LAUSR

Modeling the complete actin.myosin ATPase cycle has always been limited by the lack of experimental data concerning key steps of the cycle, because these steps can only be defined at very low ionic strength. Here, using human β-cardiac myosin-S1, we combine published data from transient and steady-state kinetics to model a minimal eight-state ATPase cycle. The model illustrates the occupancy of each intermediate around the cycle and how the occupancy is altered by changes in actin concentration for [actin] = 1-20Km. The cycle can be used to predict the maximal velocity of contraction (by motility assay or sarcomeric shortening) at different actin concentrations (which is consistent with experimental velocity data) and predict the effect of a 5 pN load on a single motor. The same exercise was repeated for human α-cardiac myosin S1 and rabbit fast skeletal muscle S1. The data illustrates how the motor domain properties can alter the ATPase cycle and hence the occupancy of the key states in the cycle. These in turn alter the predicted mechanical response of the myosin independent of other factors present in a sarcomere, such as filament stiffness and regulatory proteins. We also explore the potential of this modeling approach for the study of mutations in human β-cardiac myosin using the hypertrophic myopathy mutation R453C. Our modeling, using the transient kinetic data, predicts mechanical properties of the motor that are compatible with the single-molecule study. The modeling approach may therefore be of wide use for predicting the properties of myosin mutations.