The elastic properties of the structurally characterized myosin II S2 subdomain: a molecular dynamics and normal mode analysis.

The elastic properties (stretching and bending moduli) of myosin are expected to play an important role in its function. Of particular interest is the extended alpha-helical coiled-coil portion of the molecule. Since there is no high resolution structure for the entire coiled-coil, a study is made of the scallop myosin II S2 subdomain for which an x-ray structure is available (Protein Data Bank 1nkn). We estimate the stretching and bending moduli of the S2 subdomain with an atomic level model by use of molecular simulations. Results were obtained from nonequilibrium molecular dynamics simulations in the presence of an external force, from the fluctuations in equilibrium molecular dynamics simulations and from normal modes. In addition, a poly-Ala (78 amino acid residues) alpha-helix model was examined to test the methodology and because of its interest as part of the lever arm. As expected, both the alpha-helix and coiled-coil S2 subdomain are very stiff for stretching along the main axis, with the stretching stiffness constant in the range 60-80 pN/nm (scaled to the 60 nm long S2). Both molecules are much more flexible for bending with a lateral stiffness of approximately 0.010 pN/nm for the S2 and 0.0055 pN/nm for the alpha-helix (scaled to 60 nm). These results are expected to be useful in estimating cross-bridge elasticity, which is required for understanding the strain-dependent transitions in the actomyosin cycle and for the development of three-dimensional models of muscle contraction.