In silico modeling of cardiac contractility modulation