The myosin cardiac loop participates functionally in the actomyosin interaction

The R369 Myosin Residue within Loop 4 Is Critical for Actin Binding and Muscle Function in Drosophila

The myosin molecular motor interacts with actin filaments in an ATP-dependent manner to yield muscle contraction. Myosin heavy chain residue R369 is located within loop 4 at the actin-tropomyosin interface of myosin's upper 50 kDa subdomain. To probe the importance of R369, we introduced a histidine mutation of that residue into Drosophila myosin and implemented an integrative approach to determine effects at the biochemical, cellular, and whole organism levels. Substituting the similarly charged but bulkier histidine residue reduces maximal actin binding in vitro without affecting myosin ATPase activity. R369H mutants exhibit impaired flight ability that is dominant in heterozygotes and progressive with age in homozygotes. Indirect flight muscle ultrastructure is normal in mutant homozygotes, suggesting that assembly defects or structural deterioration of myofibrils are not causative of reduced flight. Jump ability is also reduced in homozygotes. In contrast to these skeletal muscle defects, R369H mutants show normal heart ultrastructure and function, suggesting that this residue is differentially sensitive to perturbation in different myosin isoforms or muscle types. Overall, our findings indicate that R369 is an actin binding residue that is critical for myosin function in skeletal muscles, and suggest that more severe perturbations at this residue may cause human myopathies through a similar mechanism.

Discovery of ultrafast myosin, its amino acid sequence, and structural features

Cytoplasmic streaming with extremely high velocity (∼70 μm s^-1) occurs in cells of the characean algae (Chara). Because cytoplasmic streaming is caused by myosin XI, it has been suggested that a myosin XI with a velocity of 70 μm s^-1, the fastest myosin measured so far, exists in Chara cells. However, the velocity of the previously cloned Chara corallina myosin XI (CcXI) was about 20 μm s^-1, one-third of the cytoplasmic streaming velocity in Chara Recently, the genome sequence of Chara braunii has been published, revealing that this alga has four myosin XI genes. We cloned these four myosin XI (CbXI-1, 2, 3, and 4) and measured their velocities. While the velocities of CbXI-3 and CbXI-4 motor domains (MDs) were similar to that of CcXI MD, the velocities of CbXI-1 and CbXI-2 MDs were 3.2 times and 2.8 times faster than that of CcXI MD, respectively. The velocity of chimeric CbXI-1, a functional, full-length CbXI-1 construct, was 60 μm s^-1 These results suggest that CbXI-1 and CbXI-2 would be the main contributors to cytoplasmic streaming in Chara cells and show that these myosins are ultrafast myosins with a velocity 10 times faster than fast skeletal muscle myosins in animals. We also report an atomic structure (2.8-Å resolution) of myosin XI using X-ray crystallography. Based on this crystal structure and the recently published cryo-electron microscopy structure of acto-myosin XI at low resolution (4.3-Å), it appears that the actin-binding region contributes to the fast movement of Chara myosin XI. Mutation experiments of actin-binding surface loops support this hypothesis.

Natural variant frequencies across domains from different sarcomere proteins cross-correlate to identify inter-protein contacts associated with cardiac muscle function and disease

Coordinated sarcomere proteins produce contraction force for muscle shortening. In human ventriculum they include the cardiac myosin motor (βmys), repetitively converting ATP free energy into work, and myosin binding protein C (MYBPC3) that in complex with βmys is regulatory. Single nucleotide variants (SNVs) causing hereditary heart diseases frequently target this protein pair. The βmys/MYBPC3 complex models a regulated motor and is used here to study how the proteins couple. SNVs in βmys or MYBPC3 survey human populations worldwide. Their protein expression modifies domain structure affecting phenotype and pathogenicity outcomes. When the SNV modified domain locates to inter-protein contacts it could affect complex coordination. Domains involved, one in βmys the other in MYBPC3, form coordinated domains (co-domains). Co-domain bilateral structure implies the possibility for a shared impact from SNV modification in either domain suggesting a correlated response to a common perturbation could identify their location. Genetic divergence over human populations is proposed to perturb SNV probability coupling that is detected by cross-correlation in 2D correlation genetics (2D-CG). SNV probability data and 2D-CG identify three critical sites, two in MYBPC3 with links to several domains across the βmys motor, and, one in βmys with links to the MYBPC3 regulatory domain. MYBPC3 sites are hinges sterically enabling regulatory interactions with βmys. The βmys site is the actin binding C-loop (residues 359-377). The C-loop is a trigger for actin-activated myosin ATPase and a contraction velocity modulator. Co-domain identification implies their spatial proximity suggesting a novel approach for in vivo protein complex structure determination.

Demographic model for inheritable cardiac disease

The cardiac muscle proteins, generating and regulating energy transduction during a heartbeat, assemble in the sarcomere into a cyclical machine repetitively translating actin relative to myosin filaments. Myosin is the motor transducing ATP free energy into actin movement against resisting force. Cardiac myosin binding protein C (mybpc3) regulates shortening velocity probably by transient N-terminus binding to actin while its C-terminus strongly binds the myosin filament. Inheritable heart disease associated mutants frequently modify these proteins involving them in disease mechanisms. Nonsynonymous single nucleotide polymorphisms (SNPs) cause single residue substitutions with independent characteristics (sequence location, residue substitution, human demographic, and allele frequency) hypothesized to decide dependent phenotype and pathogenicity characteristics in a feed-forward neural network model. Trial models train and validate on a dynamic worldwide SNP database for cardiac muscle proteins then predict phenotype and pathogenicity for any single residue substitution in myosin, mybpc3, or actin. A separate Bayesian model formulates conditional probabilities for phenotype or pathogenicity given independent SNP characteristics. Neural/Bayes forecasting tests SNP pathogenicity vs (in)dependent SNP characteristics to assess individualized disease risk and in particular to elucidate gender and human subpopulation bias in disease. Evident subpopulation bias in myosin SNP pathogenicities imply myosin normally engages multiple sarcomere proteins functionally. Consistent with this observation, mybpc3 forms a third actomyosin interaction competing with myosin essential light chain N-terminus suggesting a novel strain-dependent mechanism adapting myosin force-velocity to load dynamics. The working models, and the integral myosin/mybpc3 motor concept, portends the wider considerations involved in understanding heart disease as a systemic maladaptation.

Cardiac and skeletal actin substrates uniquely tune cardiac myosin strain-dependent mechanics

Cardiac ventricular myosin (βmys) translates actin by transducing ATP free energy into mechanical work during muscle contraction. Unitary βmys translation of actin is the step-size. In vitro and in vivo βmys regulates contractile force and velocity autonomously by remixing three different step-sizes with adaptive stepping frequencies. Cardiac and skeletal actin isoforms have a specific 1 : 4 stoichiometry in normal adult human ventriculum. Human adults with inheritable hypertrophic cardiomyopathy (HCM) upregulate skeletal actin in ventriculum probably compensating the diseased muscle's inability to meet demand by adjusting βmys force-velocity characteristics. βmys force-velocity characteristics were compared for skeletal versus cardiac actin substrates using ensemble in vitro motility and single myosin assays. Two competing myosin strain-sensitive mechanisms regulate step-size choices dividing single βmys mechanics into low- and high-force regimes. The actin isoforms alter myosin strain-sensitive regulation such that onset of the high-force regime, where a short step-size is a large or major contributor, is offset to higher loads probably by the unique cardiac essential light chain (ELC) N-terminus/cardiac actin contact at Glu6/Ser358. It modifies βmys force-velocity by stabilizing the ELC N-terminus/cardiac actin association. Uneven onset of the high-force regime for skeletal versus cardiac actin modulates force-velocity characteristics as skeletal/cardiac actin fractional content increases in diseased muscle.

Neural/Bayes network predictor for inheritable cardiac disease pathogenicity and phenotype

The cardiac muscle sarcomere contains multiple proteins contributing to contraction energy transduction and its regulation during a heartbeat. Inheritable heart disease mutants affect most of them but none more frequently than the ventricular myosin motor and cardiac myosin binding protein c (mybpc3). These co-localizing proteins have mybpc3 playing a regulatory role to the energy transducing motor. Residue substitution and functional domain assignment of each mutation in the protein sequence decides, under the direction of a sensible disease model, phenotype and pathogenicity. The unknown model mechanism is decided here using a method combing neural and Bayes networks. Missense single nucleotide polymorphisms (SNPs) are clues for the disease mechanism summarized in an extensive database collecting mutant sequence location and residue substitution as independent variables that imply the dependent disease phenotype and pathogenicity characteristics in 4 dimensional data points (4ddps). The SNP database contains entries with the majority having one or both dependent data entries unfulfilled. A neural network relating causes (mutant residue location and substitution) and effects (phenotype and pathogenicity) is trained, validated, and optimized using fulfilled 4ddps. It then predicts unfulfilled 4ddps providing the implicit disease model. A discrete Bayes network interprets fulfilled and predicted 4ddps with conditional probabilities for phenotype and pathogenicity given mutation location and residue substitution thus relating the neural network implicit model to explicit features of the motor and mybpc3 sequence and structural domains. Neural/Bayes network forecasting automates disease mechanism modeling by leveraging the world wide human missense SNP database that is in place and expanding.

Biochemical and bioinformatic analysis of the myosin-XIX motor domain

Mitochondrial dynamics are dependent on both the microtubule and actin cytoskeletal systems. Evidence for the involvement of myosin motors has been described in many systems, and until recently a candidate mitochondrial myosin transport motor had not been described in vertebrates. Myosin-XIX (MYO19) was predicted to represent a novel class of myosin and had previously been shown to bind to mitochondria and increase mitochondrial network dynamics when ectopically expressed. Our analyses comparing ∼40 MYO19 orthologs to ∼2000 other myosin motor domain sequences identified instances of homology well-conserved within class XIX myosins that were not found in other myosin classes, suggesting MYO19-specific mechanochemistry. Steady-state biochemical analyses of the MYO19 motor domain indicate that Homo sapiens MYO19 is a functional motor. Insect cell-expressed constructs bound calmodulin as a light chain at the predicted stoichiometry and displayed actin-activated ATPase activity. MYO19 constructs demonstrated high actin affinity in the presence of ATP in actin-co-sedimentation assays, and translocated actin filaments in gliding assays. Expression of GFP-MYO19 containing a mutation impairing ATPase activity did not enhance mitochondrial network dynamics, as occurs with wild-type MYO19, indicating that myosin motor activity is required for mitochondrial motility. The measured biochemical properties of MYO19 suggest it is a high-duty ratio motor that could serve to transport mitochondria or anchor mitochondria, depending upon the cellular microenvironment.

Myosin individualized: single nucleotide polymorphisms in energy transduction

BACKGROUND

Myosin performs ATP free energy transduction into mechanical work in the motor domain of the myosin heavy chain (MHC). Energy transduction is the definitive systemic feature of the myosin motor performed by coordinating in a time ordered sequence: ATP hydrolysis at the active site, actin affinity modulation at the actin binding site, and the lever-arm rotation of the power stroke. These functions are carried out by several conserved sub-domains within the motor domain. Single nucleotide polymorphisms (SNPs) affect the MHC sequence of many isoforms expressed in striated muscle, smooth muscle, and non-muscle tissue. The purpose of this work is to provide a rationale for using SNPs as a functional genomics tool to investigate structurefunction relationships in myosin. In particular, to discover SNP distribution over the conserved sub-domains and surmise what it implies about sub-domain stability and criticality in the energy transduction mechanism.

RESULTS

An automated routine identifying human nonsynonymous SNP amino acid missense substitutions for any MHC gene mined the NCBI SNP data base. The routine tested 22 MHC genes coding muscle and non-muscle isoforms and identified 89 missense mutation positions in the motor domain with 10 already implicated in heart disease and another 8 lacking sequence homology with a skeletal MHC isoform for which a crystallographic model is available. The remaining 71 SNP substitutions were found to be distributed over MHC with 22 falling outside identified functional sub-domains and 49 in or very near to myosin sub-domains assigned specific crucial functions in energy transduction. The latter includes the active site, the actin binding site, the rigid lever-arm, and regions facilitating their communication. Most MHC isoforms contained SNPs somewhere in the motor domain.

CONCLUSIONS

Several functional-crucial sub-domains are infiltrated by a large number of SNP substitution sites suggesting these domains are engineered by evolution to be too-robust to be disturbed by otherwise intrusive sequence changes. Two functional sub-domains are SNP-free or relatively SNP-deficient but contain many disease implicated mutants. These sub-domains are apparently highly sensitive to any missense substitution suggesting they have failed to evolve a robust sequence paradigm for performing their function.

The myosin C-loop is an allosteric actin contact sensor in actomyosin

Actin and myosin form the molecular motor in muscle. Myosin is the enzyme performing ATP hydrolysis under the allosteric control of actin such that actin binding initiates product release and force generation in the myosin power stroke. Non-equilibrium Monte Carlo molecular dynamics simulation of the power stroke suggested that a structured surface loop on myosin, the C-loop, is the actin contact sensor initiating actin activation of the myosin ATPase. Previous experimental work demonstrated C-loop binds actin and established the forward and reverse allosteric link between the C-loop and the myosin active site. Here, smooth muscle heavy meromyosin C-loop chimeras were constructed with skeletal (sCl) and cardiac (cCl) myosin C-loops substituted for the native sequence. In both cases, actin-activated ATPase inhibition is indicated mainly by the lower V(max). In vitro motility was also inhibited in the chimeras. Motility data were collected as a function of myosin surface density, with unregulated actin, and with skeletal and cardiac isoforms of tropomyosin-bound actin for the wild type, cCl, and sCl. Slow and fast subpopulations of myosin velocities in the wild-type species were discovered and represent geometrically unfavorable and favorable actomyosin interactions, respectively. Unfavorable interactions are detected at all surface densities tested. Favorable interactions are more probable at higher myosin surface densities. Cardiac tropomyosin-bound actin promotes the favorable actomyosin interactions by lowering the inhibiting geometrical constraint barriers with a structural effect on actin. Neither higher surface density nor cardiac tropomyosin-bound actin can accelerate motility velocity in cCl or sCl, suggesting the element initiating maximal myosin activation by actin resides in the C-loop.

Human myosin Vc is a low duty ratio nonprocessive motor

There are three distinct members of the myosin V family in vertebrates, and each isoform is involved in different membrane trafficking pathways. Both myosin Va and Vb have demonstrated that they are high duty ratio motors that are consistent with the processive nature of these motors. Here we report that the ATPase cycle mechanism of the single-headed construct of myosin Vc is quite different from those of other vertebrate myosin V isoforms. K(ATPase) of the actin-activated ATPase was 62 microm, which is much higher than that of myosin Va ( approximately 1 mum). The rate of ADP release from actomyosin Vc was 12.7 s(-1), which was 2 times greater than the entire ATPase cycle rate, 6.5 s(-1). P(i) burst size was 0.31, indicating that the equilibrium of the ATP hydrolysis step is shifted to the prehydrolysis form. Our kinetic model, based on all kinetic data we determined in this study, suggests that myosin Vc spends the majority of the ATPase cycle time in the weak actin binding state in contrast to myosin Va and Vb. Consistently, the two-headed myosin Vc construct did not show processive movement in total internal reflection fluorescence microscope analysis, demonstrating that myosin Vc is a nonprocessive motor. Our findings suggest that myosin Vc fulfills its function as a cargo transporter by different mechanisms from other myosin V isoforms.

Kinetic characterization of the function of myosin loop 4 in the actin-myosin interaction

Myosin interacts with actin during its enzymatic cycle, and actin stimulates myosin's ATPase activity. There are extensive interaction surfaces on both actin and myosin. Several surface loops of myosin play different roles in actomyosin interaction. However, the functional role of loop 4 in actin binding is still ambiguous. We explored the role of loop 4 by either mutating its conserved acidic group, Glu-365, to Gln (E365Q), or by replacing the entire loop with three glycines (DeltaAL) in a Dictyostelium discoideum myosin II motor domain (MD) containing a single tryptophan residue. This native tryptophan (Trp-501) is located in the relay loop and is sensitive to nucleotide binding and lever-arm movement. Fluorescence and fast kinetic measurements showed that the mutations in loop 4 do not alter the enzymatic steps of the ATPase cycle in the absence of actin. By contrast, actin binding was significantly weakened in the absence and presence of ADP and ATP in both mutants. Because the strength of actin-myosin interaction increases in the order of rigor, ADP, and ATP complex, we conclude that loop 4 is a functional actin-binding region that stabilizes actomyosin complex, particularly in weak actin-binding states.

An unusual transduction pathway in human tonic smooth muscle myosin

The motor protein myosin binds actin and ATP, producing work by causing relative translation of the proteins while transducing ATP free energy. Smooth muscle myosin has one of four heavy chains encoded by the MYH11 gene that differ at the C-terminus and in the active site for ATPase due to alternate splicing. A seven-amino-acid active site insert in phasic muscle myosin is absent from the tonic isoform. Fluorescence increase in the nucleotide sensitive tryptophan (NST) accompanies nucleotide binding and hydrolysis in several myosin isoforms implying it results from a common origin within the motor. A wild-type tonic myosin (smA) construct of the enzymatic head domain (subfragment 1 or S1) has seven tryptophan residues and nucleotide-induced fluorescence enhancement like other myosins. Three smA mutants probe the molecular basis for the fluorescence enhancement. W506+ contains one tryptophan at position 506 homologous to the NST in other myosins. W506F has the native tryptophans except phenylalanine replaces W506, and W506+(Y499F) is W506+ with phenylalanine replacing Y499. W506+ lacks nucleotide-induced fluorescence enhancement probably eliminating W506 as the NST. W506F has impaired ATPase activity but retains nucleotide-induced fluorescence enhancement. Y499F replacement in W506+ partially rescues nucleotide sensitivity demonstrating the role of Y499 as an NST facilitator. The exceptional response of W506 to active site conformation opens the possibility that phasic and tonic isoforms differ in how influences from active site ATPase propagate through the protein network.

A closer look at energy transduction in muscle

Muscular force is the sum of unitary force interactions generated as filaments of myosins move forcibly along parallel filaments of actins, understanding that the free energy required comes from myosin-catalyzed ATP hydrolysis. Using results from conventional biochemistry, our own mutational studies, and diffraction images from others, we attempt, in molecular detail, an account of a unitary interaction, i.e., what happens after a traveling myosin head, bearing an ADP-P(i), reaches the next station of an actin filament in its path. We first construct a reasonable model of the myosin head and actin regions that meet to form the "weakly bound state". Separately, we consider Holmes' model of the rigor state [Holmes, K. C., Angert, I., Kull, F. J., Jahn, W. & Schröder, R. R. (2003) Nature 425, 423-427], supplemented with several heretofore missing residues, thus realizing the "strongly bound state." Comparing states suggests how influences initiated at the interface travel elsewhere in myosin to discharge various functions, including striking the actins. Overall, state change seems to occur by attachment of a hydrophobic triplet (Trp-546, Phe-547, and Pro-548) of myosin to an actin conduit with a hydrophobic guiding rail (Ile-341, Ile-345, Leu-349, and Phe-352) and the subsequent linear movement of the triplet along the rail.

Myosin Surface Loop 4 Modulates Inhibition of Actomyosin 1b ATPase Activity by Tropomyosin

Structural studies of the class I myosin, MyoE, led to the predictions that loop 4, a surface loop near the actin-binding region that is longer in class I myosins than in other myosin subclasses, might limit binding of myosins I to actin when actin-binding proteins, like tropomyosin, are present, and might account for the exclusion of myosin I from stress fibers. To test these hypotheses, mutant molecules of the related mammalian class I myosin, Myo1b, in which loop 4 was truncated (from an amino acid sequence of RMNGLDES to NGLD) or replaced with the shorter and distinct loop 4 found in Dictyostelium myosin II (GAGEGA), were expressed in vitro and their interaction with actin and with actin-tropomyosin was tested. Saturating amounts of expressed fibroblast tropomyosin-2 resulted in a decrease in the maximum actin-activated Mg2+-ATPase activity of wild-type Myo1b but had little or no effect on the actin-activated Mg2+-ATPase activity of the two mutants. In motility assays, few actin filaments bound tightly to Myo1b-WT-coated cover slips when tropomyosin-2 was present, whereas actin filaments both bound and were translocated by Myo1b-NGLD or Myo1b-GAGEGA in both the presence and absence of tropomyosin-2. When expressed in mammalian cells, like the wild type, the mutant myosins were largely excluded from tropomyosin-containing actin filaments, indicating that in the cell additional factors besides loop 4 determine targeting of myosins I to specific subpopulations of actin filaments.

Fluorescence labeling and computational analysis of the strut of myosin’s 50kDa cleft

A new fluorescent labeling procedure specific for the strut sequence of myosin subfragment-1's 50kDa cleft was developed using CY3 N-hydroxy succinimidyl ester as a hydrophobic tag and hydrophobic interaction chromatography to purify the major labeled species which retained actin-activated ATPase activity. Stern-Volmer analysis suggests that the CY3 is in close proximity to basic residues, consistent with inspection of the mapped labeling site in the atomic model. Fluorescence polarization indicates that the CY3 becomes more mobile upon actin binding, supporting a location near the actomyosin interface. In contrast, nucleotide binding to myosin had little impact on the CY3. Molecular mechanics and stochastic dynamics simulations suggest that this labeling site is sensitive to forced cleft opening and closure, but the upper 50kDa cleft does not move easily. In addition, there appear to be some long-range effects of forced cleft opening and closing that could impact the lever arm position.

Protein-protein interactions in actin-myosin binding and structural effects of R405Q mutation: a molecular dynamics study

Detailed residue-wise interactions involved in the binding of myosin to actin in the rigor conformation without nucleotides have been examined using molecular dynamics simulations of the chicken skeletal myosin head complexed with two actin monomers, based on the cryo-microscopic model of Holmes et al. (Nature 2003;425:423-427). The overall interaction is largely electrostatic in nature, because of the charged residues in the four loops surrounding the central primary binding site. The 50k/20k loop, disordered in crystal structures and in simulations of free myosin in solution, was found to be in a conformation stabilized with 1 - 2 internal salt bridges. The cardiomyopathy loop forms 2 - 3 interprotein salt bridges with actin monomers upon binding, whereas its Arg405 residue, the mutation site associated with the hypertrophic cardiomyopathy, forms a strong salt bridge with Glu605 in the neighboring helix away from actin in the actin-bound myosin. The myopathy loop of the R405Q mutant maintains a high degree of two-strand beta-sheet character when bound to actin with the corresponding salt bridges broken.

Isolation of myosin XI genes from the Closterium peracerosum-strigosum-littorale complex and analysis of their expression during sexual reproduction

Myosins comprise a large superfamily of molecular motors that generate mechanical force in ATP-dependent interactions with actin filaments. On the basis of their conserved motor domain sequences, myosins can be divided into at least 17 classes, 3 of which (VIII, XI, XIII) are found in plants. Although full sequences of myosins are available from several species of green plants, little is known about the functions of these proteins. Additionally, sequence information for algal myosin is incomplete, and little attention has been given to the molecular evolution of myosin from green plants. In the present study, the Closterium peracerosum-strigosum-littorale complex was used as a model system for investigating a unicellular basal charophycean alga. This organism has been well studied with respect to sexual reproduction between its two mating types. Three types of partial sequences belonging to class XI myosins were obtained using degenerate primers designed to amplify motor domain sequences. Real-time polymerase chain reaction analysis of the respective myosin genes during various stages of the algal life cycle showed that one of the genes was more highly expressed during sexual reproduction, and that expression was cell-cycle-dependent in vegetatively grown cells.