Bottom-up mechanokinetic models predict ensemble function of actin and myosin based on parameter values derived from studies using isolated proteins. To be generally useful, e.g., to analyze disease effects, such models must also be able to predict ensemble function when actomyosin interaction kinetics are modified differently from normal. Here, we test this capability for a model recently shown to predict several physiological phenomena along with the effects of the small molecular compound blebbistatin. We demonstrate that this model also qualitatively predicts effects of other well-characterized drugs as well as varied concentrations of MgATP. However, the effects of one compound, amrinone, are not well accounted for quantitatively. We therefore systematically varied key model parameters to address this issue, leading to the increased amplitude of the second sub-stroke of the power stroke from 1 nm to 2.2 nm, an unchanged first sub-stroke (5.3−5.5 nm), and an effective cross-bridge attachment rate that more than doubled. In addition to better accounting for the effects of amrinone, the modified model also accounts well for normal physiological ensemble function. Moreover, a Monte Carlo simulation-based version of the model was used to evaluate force−velocity data from small myosin ensembles. We discuss our findings in relation to key aspects of actin−myosin operation mechanisms causing a non-hyperbolic shape of the force−velocity relationship at high loads. We also discuss remaining limitations of the model, including uncertainty of whether the cross-bridge elasticity is linear or not, the capability to account for contractile properties of very small actomyosin ensembles (<20 myosin heads), and the mechanism for requirements of a higher cross-bridge attachment rate during shortening compared to during isometric contraction.
Keywords: actin; mechanokinetic model; muscle contraction; myosin; myosin-active compounds; statistical model.