Mechanochemical model for myosin II dimer that can explain the spontaneous oscillatory contraction of muscle

Sun Wei^{1, 2}, Zhao Xiao-Yang^{3, 4}, Zhang Jun-Ping¹, Tala ¹, Guo Wei-Sheng^{1, †}

(color online) Proportion of myosin II dimer in different states. The values of rate constants are k₁ = 0.1 s⁻¹, k₂ = 0.5 s⁻¹, , k₄ = 22 s⁻¹, k₅ = 10 s⁻¹, k₆ = 800 s⁻¹, k₇ = 200 s⁻¹, k₈ = 200 s⁻¹, , a = 0.2 mM, b = 0.1, and c = 0.5 mM. The initial values of x_i are x₁(0) = 0.105, x₂(0) = 0.041, x₃(0) = 0.029, x₄(0) = 0.061, x₅(0) = 0.295, x₆(0) = 0.388, x₇(0) = 0.021, x₈(0) = 0.029, and x₉(0) = 0.031; y(0) = 0.039 mM.