Analysis of spontaneous oscillations for a three-state power-stroke model.

Abstract

Our study considers the mechanism of the spontaneous oscillations of molecular motors that are driven by the power stroke principle by applying linear stability analysis around the stationary solution. By representing the coupling equation of microscopic molecular motor dynamics and mesoscopic sarcomeric dynamics by a rank-1 updated matrix system, we derived the analytical representations of the eigenmodes of the Jacobian matrix that cause the oscillation. Based on these analytical representations, we successfully derived the essential conditions for the oscillation in terms of the rate constants of the power stroke and the reversal stroke transitions of the molecular motor. Unlike the two-state model, in which the dependence of the detachment rates on the motor coordinates or the applied forces on the motors plays a key role for the oscillation, our three-state power stroke model demonstrates that the dependence of the rate constants of the power and reversal strokes on the strains in the elastic elements in the motor molecules plays a key role, where these rate constants are rationally determined from the free energy available for the power stroke, the stiffness of the elastic element in the molecular motor, and the working stroke size. By applying the experimentally confirmed values to the free energy, the stiffness, and the working stroke size, our numerical model reproduces well the experimentally observed oscillatory behavior. Furthermore, our analysis shows that two eigenmodes with real positive eigenvalues characterize the oscillatory behavior, where the eigenmode with the larger eigenvalue indicates the transient of the system of the quick sarcomeric lengthening induced by the collective reversal strokes, and the smaller eigenvalue correlates with the speed of sarcomeric shortening, which is much slower than lengthening. Applying the perturbation analyses with primal physical parameters, we find that these two real eigenvalues occur on two branches derived from a merge point of a pair of complex-conjugate eigenvalues generated by Hopf bifurcation.