Early stages of energy transduction by myosin: roles of Arg in switch I, of Glu in switch II, and of the salt-bridge between them

The Local Environment of Loop Switch 1 Modulates the Rate of ATP-Induced Dissociation of Human Cardiac Actomyosin

Two isoforms of human cardiac myosin, alpha and beta, share significant sequence similarities but show different kinetics. The alpha isoform is a faster motor; it spends less time being strongly bound to actin during the actomyosin cycle. With alpha isoform, actomyosin dissociates faster upon ATP binding, and the affinity of ADP to actomyosin is weaker. One can suggest that the isoform-specific actomyosin kinetics is regulated at the nucleotide binding site of human cardiac myosin. Myosin is a P-loop ATPase; the nucleotide-binding site consists of P-loop and loops switch 1 and 2. All three loops position MgATP for successful hydrolysis. Loops sequence is conserved in both myosin isoforms, and we hypothesize that the isoform-specific structural element near the active site regulates the rate of nucleotide binding and release. Previously we ran molecular dynamics simulations and found that loop S291-E317 near loop switch 1 is more compact and exhibits larger fluctuations of the position of amino acid residues in beta isoform than in alpha. In alpha isoform, the loop forms a salt bridge with loop switch 1, the bridge is not present in beta isoform. Two isoleucines I303 and I313 of loop S291-E317 are replaced with valines in alpha isoform. We introduced a double mutation I303V:I313V in beta isoform background and studied how the mutation affects the rate of ATP binding and ADP dissociation from actomyosin. We found that ATP-induced actomyosin dissociation occurs faster in the mutant, but the rate of ADP release remains the same as in the wild-type beta isoform. Due to the proximity of loop S291-E317 and loop switch 1, a faster rate of ATP-induced actomyosin dissociation indicates that loop S291-E317 affects structural dynamics of loop switch 1, and that loop switch 1 controls ATP binding to the active site. A similar rate of ADP dissociation from actomyosin in the mutant and wild-type myosin constructs indicates that loop switch 1 does not control ADP release from actomyosin.

Myosin XVI

Myosin XVI (Myo16), a vertebrate-specific motor protein, is a recently discovered member of the myosin superfamily. The detailed functionality regarding myosin XVI requires elucidating or clarification; however, it appears to portray an important role in neural development and in the proper functioning of the nervous system. It is expressed in the largest amount in neural tissues in the late embryonic-early postnatal period, specifically the time in which neuronal cell migration and dendritic elaboration coincide. The impaired expression of myosin XVI has been found lurking in the background of several neuropsychiatric disorders including autism, schizophrenia and/or bipolar disorders.Two principal isoforms of class XVI myosins have been thus far described: Myo16a, the tailless cytoplasmic isoform and Myo16b, the full-length molecule featuring both cytoplasmic and nuclear localization. Both isoforms contain a class-specific N-terminal ankyrin repeat domain that binds to the protein phosphatase catalytic subunit. Myo16b, the predominant isoform, exhibits a diverse function. In the cytoplasm, it participates in the reorganization of the actin cytoskeleton through activation of the PI3K pathway and the WAVE-complex, while in the nucleus it may possess a role in cell cycle regulation. Based on the sequence, myosin XVI may have a compromised ATPase activity, implying a potential stationary role.

Analysis of Phosphoryl-Transfer Enzymes with QM/MM Free Energy Simulations

We discuss the application of quantum mechanics/molecular mechanics (QM/MM) free energy simulations to the analysis of phosphoryl transfers catalyzed by two enzymes: alkaline phosphatase and myosin. We focus on the nature of the transition state and the issue of mechanochemical coupling, respectively, in the two enzymes. The results illustrate unique insights that emerged from the QM/MM simulations, especially concerning the interpretation of experimental data regarding the nature of enzymatic transition states and coupling between global structural transition and catalysis in the active site. We also highlight a number of technical issues worthy of attention when applying QM/MM free energy simulations, and comment on a number of technical and mechanistic issues that require further studies.

Influence of Glu/Arg, Asp/Arg, and Glu/Lys Salt Bridges on α-Helical Stability and Folding Kinetics

Using a combination of ultraviolet circular dichroism, temperature-jump transient-infrared spectroscopy, and molecular dynamics simulations, we investigate the effect of salt bridges between different types of charged amino-acid residue pairs on α-helix folding. We determine the stability and the folding and unfolding rates of 12 alanine-based α-helical peptides, each of which has a nearly identical composition containing three pairs of positively and negatively charged residues (either Glu(-)/Arg(+), Asp(-)/Arg(+), or Glu(-)/Lys(+)). Within each set of peptides, the distance and order of the oppositely charged residues in the peptide sequence differ, such that they have different capabilities of forming salt bridges. Our results indicate that stabilizing salt bridges (in which the interacting residues are spaced and ordered such that they favor helix formation) speed up α-helix formation by up to 50% and slow down the unfolding of the α-helix, whereas salt bridges with an unfavorable geometry have the opposite effect. Comparing the peptides with different types of charge pairs, we observe that salt bridges between side chains of Glu(-) and Arg(+) are most favorable for the speed of folding, probably because of the larger conformational space of the salt-bridging Glu(-)/Arg(+) rotamer pairs compared to Asp(-)/Arg(+) and Glu(-)/Lys(+). We speculate that the observed impact of salt bridges on the folding kinetics might explain why some proteins contain salt bridges that do not stabilize the final, folded conformation.

Regulation and Plasticity of Catalysis in Enzymes: Insights from Analysis of Mechanochemical Coupling in Myosin

The mechanism of ATP hydrolysis in the myosin motor domain is analyzed using a combination of DFTB3/CHARMM simulations and enhanced sampling techniques. The motor domain is modeled in the pre-powerstroke state, in the post-rigor state, and as a hybrid based on the post-rigor state with a closed nucleotide-binding pocket. The ATP hydrolysis activity is found to depend on the positioning of nearby water molecules, and a network of polar residues facilitates proton transfer and charge redistribution during hydrolysis. Comparison of the observed hydrolysis pathways and the corresponding free energy profiles leads to detailed models for the mechanism of ATP hydrolysis in the pre-powerstroke state and proposes factors that regulate the hydrolysis activity in different conformational states. In the pre-powerstroke state, the scissile P_γ-O_3β bond breaks early in the reaction. Proton transfer from the lytic water to the γ-phosphate through active site residues is an important part of the kinetic bottleneck; several hydrolysis pathways that feature distinct proton transfer routes are found to have similar free energy barriers, suggesting a significant degree of plasticity in the hydrolysis mechanism. Comparison of hydrolysis in the pre-powerstroke state and the closed post-rigor model suggests that optimization of residues beyond the active site for electrostatic stabilization and preorganization is likely important to enzyme design.

Review: The ATPase mechanism of myosin and actomyosin

Myosins are a large family of molecular motors that use the common P-loop, Switch 1 and Switch 2 nucleotide binding motifs to recognize ATP, to create a catalytic site than can efficiently hydrolyze ATP and to communicate the state of the nucleotide pocket to other allosteric binding sites on myosin. The energy of ATP hydrolysis is used to do work against an external load. In this short review I will outline current thinking on the mechanism of ATP hydrolysis and how the energy of ATP hydrolysis is coupled to a series of protein conformational changes that allow a myosin, with the cytoskeleton track actin, to operate as a molecular motor of distinct types; fast movers, processive motors or strain sensors. © 2016 Wiley Periodicals, Inc. Biopolymers 105: 483-491, 2016.

Kinetic Adaptations of Myosins for Their Diverse Cellular Functions

Members of the myosin superfamily are involved in all aspects of eukaryotic life. Their function ranges from the transport of organelles and cargos to the generation of membrane tension, and the contraction of muscle. The diversity of physiological functions is remarkable, given that all enzymatically active myosins follow a conserved mechanoenzymatic cycle in which the hydrolysis of ATP to ADP and inorganic phosphate is coupled to either actin-based transport or tethering of actin to defined cellular compartments. Kinetic capacities and limitations of a myosin are determined by the extent to which actin can accelerate the hydrolysis of ATP and the release of the hydrolysis products and are indispensably linked to its physiological tasks. This review focuses on kinetic competencies that - together with structural adaptations - result in myosins with unique mechanoenzymatic properties targeted to their diverse cellular functions.

Biological Nanomotors with a Revolution, Linear, or Rotation Motion Mechanism

The ubiquitous biological nanomotors were classified into two categories in the past: linear and rotation motors. In 2013, a third type of biomotor, revolution without rotation (http://rnanano.osu.edu/movie.html), was discovered and found to be widespread among bacteria, eukaryotic viruses, and double-stranded DNA (dsDNA) bacteriophages. This review focuses on recent findings about various aspects of motors, including chirality, stoichiometry, channel size, entropy, conformational change, and energy usage rate, in a variety of well-studied motors, including FoF1 ATPase, helicases, viral dsDNA-packaging motors, bacterial chromosome translocases, myosin, kinesin, and dynein. In particular, dsDNA translocases are used to illustrate how these features relate to the motion mechanism and how nature elegantly evolved a revolution mechanism to avoid coiling and tangling during lengthy dsDNA genome transportation in cell division. Motor chirality and channel size are two factors that distinguish rotation motors from revolution motors. Rotation motors use right-handed channels to drive the right-handed dsDNA, similar to the way a nut drives the bolt with threads in same orientation; revolution motors use left-handed motor channels to revolve the right-handed dsDNA. Rotation motors use small channels (<2 nm in diameter) for the close contact of the channel wall with single-stranded DNA (ssDNA) or the 2-nm dsDNA bolt; revolution motors use larger channels (>3 nm) with room for the bolt to revolve. Binding and hydrolysis of ATP are linked to different conformational entropy changes in the motor that lead to altered affinity for the substrate and allow work to be done, for example, helicase unwinding of DNA or translocase directional movement of DNA.

The effect of novel mutations on the structure and enzymatic activity of unconventional myosins associated with autosomal dominant non-syndromic hearing loss

Mutations in five unconventional myosin genes have been associated with genetic hearing loss (HL). These genes encode the motor proteins myosin IA, IIIA, VI, VIIA and XVA. To date, most mutations in myosin genes have been found in the Caucasian population. In addition, only a few functional studies have been performed on the previously reported myosin mutations. We performed screening and functional studies for mutations in the MYO1A and MYO6 genes in Korean cases of autosomal dominant non-syndromic HL. We identified four novel heterozygous mutations in MYO6. Three mutations (p.R825X, p.R991X and Q918fsX941) produce a premature truncation of the myosin VI protein. Another mutation, p.R205Q, was associated with diminished actin-activated ATPase activity and actin gliding velocity of myosin VI in an in vitro analysis. This finding is consistent with the results of protein modelling studies and corroborates the pathogenicity of this mutation in the MYO6 gene. One missense variant, p.R544W, was found in the MYO1A gene, and in silico analysis suggested that this variant has deleterious effects on protein function. This finding is consistent with the results of protein modelling studies and corroborates the pathogenic effect of this mutation in the MYO6 gene.

Catalytic strategy used by the myosin motor to hydrolyze ATP

Myosin is a molecular motor responsible for biological motions such as muscle contraction and intracellular cargo transport, for which it hydrolyzes adenosine 5'-triphosphate (ATP). Early steps of the mechanism by which myosin catalyzes ATP hydrolysis have been investigated, but still missing are the structure of the final ADP·inorganic phosphate (Pi) product and the complete pathway leading to it. Here, a comprehensive description of the catalytic strategy of myosin is formulated, based on combined quantum-classical molecular mechanics calculations. A full exploration of catalytic pathways was performed and a final product structure was found that is consistent with all experiments. Molecular movies of the relevant pathways show the different reorganizations of the H-bond network that lead to the final product, whose γ-phosphate is not in the previously reported HPγO4(2-) state, but in the H2PγO4(-) state. The simulations reveal that the catalytic strategy of myosin employs a three-pronged tactic: (i) Stabilization of the γ-phosphate of ATP in a dissociated metaphosphate (PγO3(-)) state. (ii) Polarization of the attacking water molecule, to abstract a proton from that water. (iii) Formation of multiple proton wires in the active site, for efficient transfer of the abstracted proton to various product precursors. The specific role played in this strategy by each of the three loops enclosing ATP is identified unambiguously. It explains how the precise timing of the ATPase activation during the force generating cycle is achieved in myosin. The catalytic strategy described here for myosin is likely to be very similar in most nucleotide hydrolyzing enzymes.

Switch II mutants reveal coupling between the nucleotide- and actin-binding regions in myosin V

Conserved active-site elements in myosins and other P-loop NTPases play critical roles in nucleotide binding and hydrolysis; however, the mechanisms of allosteric communication among these mechanoenzymes remain unresolved. In this work we introduced the E442A mutation, which abrogates a salt-bridge between switch I and switch II, and the G440A mutation, which abolishes a main-chain hydrogen bond associated with the interaction of switch II with the γ phosphate of ATP, into myosin V. We used fluorescence resonance energy transfer between mant-labeled nucleotides or IAEDANS-labeled actin and FlAsH-labeled myosin V to examine the conformation of the nucleotide- and actin-binding regions, respectively. We demonstrate that in the absence of actin, both the G440A and E442A mutants bind ATP with similar affinity and result in only minor alterations in the conformation of the nucleotide-binding pocket (NBP). In the presence of ADP and actin, both switch II mutants disrupt the formation of a closed NBP actomyosin.ADP state. The G440A mutant also prevents ATP-induced opening of the actin-binding cleft. Our results indicate that the switch II region is critical for stabilizing the closed NBP conformation in the presence of actin, and is essential for communication between the active site and actin-binding region.

Minimum energy reaction profiles for ATP hydrolysis in myosin

The minimum energy reaction profiles corresponding to two possible reaction mechanisms of adenosine triphosphate (ATP) hydrolysis in myosin are computed in this work within the framework of the quantum mechanics-molecular mechanics (QM/MM) method by using the same partitioning of the model system to the QM and MM parts and the same computational protocol. On the first reaction route, one water molecule performs nucleophilic attack at the phosphorus center P(γ) from ATP while the second water molecule in the closed protein cleft serves as a catalytic base assisted by the Glu residue from the myosin salt bridge. According to the present QM/MM calculations consistent with the results of kinetic studies this reaction pathway is characterized by a low activation energy barrier about 10 kcal/mol. The computed activation energy barrier for the second mechanism, which assumes the penta-coordinated oxyphosphorane transition state upon involvement of single water molecule in the reaction, is considerably higher than that for the two-water mechanism.

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.

ATP hydrolysis in Eg5 kinesin involves a catalytic two-water mechanism

Motor proteins couple steps in ATP binding and hydrolysis to conformational switching both in and remote from the active site. In our kinesin.AMPPPNP crystal structure, closure of the active site results in structural transformations appropriate for microtubule binding and organizes an orthosteric two-water cluster. We conclude that a proton is shared between the lytic water, positioned for gamma-phosphate attack, and a second water that serves as a general base. To our knowledge, this is the first experimental detection of the catalytic base for any ATPase. Deprotonation of the second water by switch residues likely triggers subsequent large scale structural rearrangements. Therefore, the catalytic base is responsible for initiating nucleophilic attack of ATP and for relaying the positive charge over long distances to initiate mechanotransduction. Coordination of switch movements via sequential proton transfer along paired water clusters may be universal for nucleotide triphosphatases with conserved active sites, such as myosins and G-proteins.

The hydrolysis activity of adenosine triphosphate in myosin: a theoretical analysis of anomeric effects and the nature of the transition state

Combined quantum mechanical/molecular mechanical (QM/MM) calculations with density functional theory are employed to analyze two issues related to the hydrolysis activity of adenosine triphosphate (ATP) in myosin. First, we compare the geometrical properties and electronic structure of ATP in the open (post-rigor) and closed (pre-powerstroke) active sites of the myosin motor domain. Compared to both solution and the open active site cases, the scissile P(gamma)-O(3beta) bond of ATP in the closed active site is shown to be substantially elongated. Natural bond orbital (NBO) analysis clearly shows that this structural feature is correlated with the stronger anomeric effects in the closed active site, which involve charge transfers from the lone pairs in the nonbridging oxygen in the gamma-phosphate to the antibonding orbital of the scissile bond. However, an energetic analysis finds that the ATP molecule is not significantly destabilized by the P(gamma)-O(3beta) bond elongation. Therefore, despite the notable perturbations in the geometry and electronic structure of ATP as its environment changes from solution to the hydrolysis-competent active site, ground-state destabilization is unlikely to play a major role in enhancing the hydrolysis activity in myosin. Second, two-dimensional potential energy maps are used to better characterize the energetic landscape near the hydrolysis transition state. The results indicate that the transition-state region is energetically flat and a range of structures representative of different mechanisms according to the classical nomenclature (e.g., "associative", "dissociative", and "concerted") are very close in energy. Therefore, at least in the case of ATP hydrolysis in myosin, the energetic distinction between different reaction mechanisms following the conventional nomenclature is likely small. This study highlights the importance of (i) explicitly evaluating the relevant energetic properties for determining whether a factor is essential to catalysis and (ii) broader explorations of the energy landscape beyond saddle points (even on free-energy surface) for characterizing the molecular mechanism of catalysis.

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.

Electrostatics-defying interaction between arginine termini as a thermodynamic driving force in protein-protein interaction

Apparent electrostatics-defying clustering of arginines attributed as screening effect of solvent is in this study examined as a possible thermodynamic driving force in protein-protein interaction. A dataset of 266 protein dimers is found to have approximately 22% arginines mutually paired and approximately 17% pairs in interaction across interfaces and thus putative "hotspots" of protein-protein interaction. The pairing, uncorrelated with inter or intramolecular context, could be contributing in protein folding as well, and, uncorrelated with solvent access, could be driven by effects that are generic to solvent and protein structures. Mutually stacked at shorter distances but in diverse geometrical modes otherwise, the cations tend to be in gross deficit of hydrogen-bond partners, and contributing electrostatics across protein-protein interface that, on average, is repulsive for protein-protein interaction. Embedded in local environment enriched in polarizable residues, aromatic, aliphatic, and anionic, the arginines may contribute to protein-protein interaction via environmental polarization response to electrostatics of cation clustering, a possible new principle in molecular recognition.

Extensive conformational transitions are required to turn on ATP hydrolysis in myosin

Conventional myosin is representative of biomolecular motors in which the hydrolysis of adenosine triphosphate (ATP) is coupled to large-scale structural transitions both in and remote from the active site. The mechanism that underlies such "mechanochemical coupling," especially the causal relationship between hydrolysis and allosteric structural changes, has remained elusive despite extensive experimental and computational analyses. In this study, using combined quantum mechanical and molecular mechanical simulations and different conformations of the myosin motor domain, we provide evidence to support that regulation of ATP hydrolysis activity is not limited to residues in the immediate environment of the phosphate. Specifically, we illustrate that efficient hydrolysis of ATP depends not only on the proper orientation of the lytic water but also on the structural stability of several nearby residues, especially the Arg238-Glu459 salt bridge (the numbering of residues follows myosin II in Dictyostelium discoideum) and the water molecule that spans this salt bridge and the lytic water. More importantly, by comparing the hydrolysis activities in two motor conformations with very similar active-site (i.e., Switches I and II) configurations, which distinguished this work from our previous study, the results clearly indicate that the ability of these residues to perform crucial electrostatic stabilization relies on the configuration of residues in the nearby N-terminus of the relay helix and the "wedge loop." Without the structural support from those motifs, residues in a closed active site in the post-rigor motor domain undergo subtle structural variations that lead to consistently higher calculated ATP hydrolysis barriers than in the pre-powerstroke state. In other words, starting from the post-rigor state, turning on the ATPase activity requires not only displacement of Switch II to close the active site but also structural transitions in the N-terminus of the relay helix and the "wedge loop," which have been proposed previously to be ultimately coupled to the rotation of the converter subdomain 40 A away.

Impacts of Usher syndrome type IB mutations on human myosin VIIa motor function

Usher syndrome (USH) is a human hereditary disorder characterized by profound congenital deafness, retinitis pigmentosa, and vestibular dysfunction. Myosin VIIa has been identified as the responsible gene for USH type 1B, and a number of missense mutations have been identified in the affected families. However, the molecular basis of the dysfunction of USH gene, myosin VIIa, in the affected families is unknown to date. Here we clarified the effects of USH1B mutations on human myosin VIIa motor function for the first time. The missense mutations of USH1B significantly inhibited the actin activation of ATPase activity of myosin VIIa. G25R, R212C, A397D, and E450Q mutations abolished the actin-activated ATPase activity completely. P503L mutation increased the basal ATPase activity for 2-3-fold but reduced the actin-activated ATPase activity to 50% of the wild type. While all of the mutations examined, except for R302H, reduced the affinity for actin and the ATP hydrolysis cycling rate, they did not largely decrease the rate of ADP release from actomyosin, suggesting that the mutations reduce the duty ratio of myosin VIIa. Taken together, the results suggest that the mutations responsible for USH1B cause the complete loss of the actin-activated ATPase activity or the reduction of duty ratio of myosin VIIa.

Mechanism of the myosin catalyzed hydrolysis of ATP as rationalized by molecular modeling

The intrinsic chemical reaction of adenosine triphosphate (ATP) hydrolysis catalyzed by myosin is modeled by using a combined quantum mechanics and molecular mechanics (QM/MM) methodology that achieves a near ab initio representation of the entire model. Starting with coordinates derived from the heavy atoms of the crystal structure (Protein Data Bank ID code 1VOM) in which myosin is bound to the ATP analog ADP.VO(4)(-), a minimum-energy path is found for the transformation ATP + H(2)O --> ADP + P(i) that is characterized by two distinct events: (i) a low activation-energy cleavage of the P(gamma) O(betagamma) bond and separation of the gamma-phosphate from ADP and (ii) the formation of the inorganic phosphate as a consequence of proton transfers mediated by two water molecules and assisted by the Glu-459-Arg-238 salt bridge of the protein. The minimum-energy model of the enzyme-substrate complex features a stable hydrogen-bonding network in which the lytic water is positioned favorably for a nucleophilic attack of the ATP gamma-phosphate and for the transfer of a proton to stably bound second water. In addition, the P(gamma) O(betagamma) bond has become significantly longer than in the unbound state of the ATP and thus is predisposed to cleavage. The modeled transformation is viewed as the part of the overall hydrolysis reaction occurring in the closed enzyme pocket after ATP is bound tightly to myosin and before conformational changes preceding release of inorganic phosphate.

Mechanochemical coupling in the myosin motor domain. I. Insights from equilibrium active-site simulations

Although the major structural transitions in molecular motors are often argued to couple to the binding of Adenosine triphosphate (ATP), the recovery stroke in the conventional myosin has been shown to be dependent on the hydrolysis of ATP. To obtain a clearer mechanistic picture for such "mechanochemical coupling" in myosin, equilibrium active-site simulations with explicit solvent have been carried out to probe the behavior of the motor domain as functions of the nucleotide chemical state and conformation of the converter/relay helix. In conjunction with previous studies of ATP hydrolysis with different active-site conformations and normal mode analysis of structural flexibility, the results help establish an energetics-based framework for understanding the mechanochemical coupling. It is proposed that the activation of hydrolysis does not require the rotation of the lever arm per se, but the two processes are tightly coordinated because both strongly couple to the open/close transition of the active site. The underlying picture involves shifts in the dominant population of different structural motifs as a consequence of changes elsewhere in the motor domain. The contribution of this work and the accompanying paper [] is to propose the actual mechanism behind these "population shifts" and residues that play important roles in the process. It is suggested that structural flexibilities at both the small and large scales inherent to the motor domain make it possible to implement tight couplings between different structural motifs while maintaining small free-energy drops for processes that occur in the detached states, which is likely a feature shared among many molecular motors. The significantly different flexibility of the active site in different X-ray structures with variable level arm orientations supports the notation that external force sensed by the lever arm may transmit into the active site and influence the chemical steps (nucleotide hydrolysis and/or binding).

Myosin16b: The COOH-tail region directs localization to the nucleus and overexpression delays S-phase progression

Rat Myo16a and Myo16b comprise the founding members of class XVI myosin and are characterized by an N-terminal ankyrin repeat domain thought to mediate an association with protein phosphatase 1 catalytic subunits 1alpha and 1gamma. Myo16b is the principal isoform and reveals predominant expression in developing neural tissue. Here, we use COS-7 cells as a model system to develop an understanding of Myo16b function. We find that Myo16b displays predominant localization in the nucleus of cells transitioning through interphase, but is not associated with processes of mitosis. Using a panel of EGFP-Myo16b-expression plasmids in transient transfection studies, we identified the COOH-terminal residues 1616-1912 as necessary and solely sufficient to target Myo16b to the nucleus. We show that the Myo16b-tail region directs localization to a nuclear compartment containing profilin and polymerized actin, which appears to form a three-dimensional meshwork through the depth of the nucleus. Further, we demonstrate that this compartment localizes within euchromatic regions of the genome and contains proliferating cell nuclear antigen (PCNA) and cyclin A, both markers of S-phase of the cell cycle. Cells transiently expressing Myo16b or Myo16b-tail region show limited incorporation of BrdU, delayed progression through S-phase of the cell cycle, and curtailed cellular proliferation.

Quantum Chemical Studies of the Myosin ATPase Mechanism

Muscle contraction is one of the molecular motions in the living body. The elemental process of muscle contraction is a mutual sliding of two kinds of filaments of proteins, myosin and actin, using energy liberated by ATP hydrolysis. Myosin enzymatically hydrolyzes ATP and interacts with actin. This article reviews previous quantum chemical studies to help determine how the chemical energy of ATP hydrolysis is converted to mechanical energy during sliding.

Analysis of functional motions in Brownian molecular machines with an efficient block normal mode approach: myosin-II and Ca2+ -ATPase

The structural flexibilities of two molecular machines, myosin and Ca(2+)-ATPase, have been analyzed with normal mode analysis and discussed in the context of their energy conversion functions. The normal mode analysis with physical intermolecular interactions was made possible by an improved implementation of the block normal mode (BNM) approach. The BNM results clearly illustrated that the large-scale conformational transitions implicated in the functional cycles of the two motor systems can be largely captured with a small number of low-frequency normal modes. Therefore, the results support the idea that structural flexibility is an essential part of the construction principle of molecular motors through evolution. Such a feature is expected to be more prevalent in motor proteins than in simpler systems (e.g., signal transduction proteins) because in the former, large-scale conformational transitions often have to occur before the chemical events (e.g., ATP hydrolysis in myosin and ATP binding/phosphorylation in Ca(2+)-ATPase). This highlights the importance of Brownian motions associated with the protein domains that are involved in the functional transitions; in this sense, Brownian molecular machines is an appropriate description of molecular motors, although the normal mode results do not address the origin of the ratchet effect. The results also suggest that it might be more appropriate to describe functional transitions in some molecular motors as intrinsic elastic motions modulating local structural changes in the active site, which in turn gets stabilized by the subsequent chemical events, in contrast with the conventional idea of local changes somehow getting amplified into larger-scale motions. In the case of myosin, for example, we favor the idea that Brownian motions associated with the flexible converter propagates to the Switch I/II region, where the salt-bridge formation gets stabilized by ATP hydrolysis, in contrast with the textbook notion that ATP hydrolysis drives the converter motion. Another useful aspect of the BNM results is that selected low-frequency normal modes have been identified to form a set of collective coordinates that can be used to characterize the progress of a significant fraction of large-scale conformational transitions. Therefore, the present normal mode analysis has provided a stepping-stone toward more elaborate microscopic simulations for addressing critical issues in free energy conversions in molecular machines, such as the coupling and the causal relationship between collective motions and essential local changes at the catalytic active site where ATP hydrolysis occurs.

A hypothesis about myosin catalysis

When ATP binds to the active site of myosin heads, Switch II undergoes a large conformational change and the cleft surrounding the bound gamma-phosphate closes. In the closed state, Glu470 in Switch II comes together with Arg247 in Switch I to form a salt-bridge. Here, the functional significance of the two bridging residues was tested by using site-directed mutagenesis. We conclude from such tests that (a) the attractive force between Arg247 and the gamma-phosphate of ATP moves the cleft to close, and (b) during hydrolysis, Glu470 is intimately involved in positioning the lytic water for the attack on the gamma-phosphorus. We also speculate on how the salt-bridge between Arg247 and Glu470 is related to hydrolysis.

A kinesin switch I arginine to lysine mutation rescues microtubule function

Switch I and II are key active site structural elements of kinesins, myosins, and G-proteins. Our analysis of a switch I mutant (R210A) in Drosophila melanogaster kinesin showed a reduction in microtubule affinity, a loss in cooperativity between the motor domains, and an ATP hydrolysis defect leading to aberrant detachment from the microtubule. To investigate the conserved arginine in switch I further, a lysine substitution mutant was generated. The R210K dimeric motor has lost the ability to hydrolyze ATP; however, it has rescued microtubule function. Our results show that R210K has restored microtubule association kinetics, microtubule affinity, ADP release kinetics, and motor domain cooperativity. Moreover, the active site at head 1 is able to distinguish ATP, ADP, and AMP-PNP to signal head 2 to bind the microtubule and release mantADP with kinetics comparable with wild-type. Therefore, the structural pathway of communication from head 1 to head 2 is restored, and head 2 can respond to this signal by binding the microtubule and releasing mantADP. Structural modeling revealed that lysine could retain some of the hydrogen bonds made by arginine but not all, suggesting a structural hypothesis for the ability of lysine to rescue microtubule function in the Arg210 mutant.

The hydration structure of guanidinium and thiocyanate ions: implications for protein stability in aqueous solution

Neutron diffraction experiments were carried out on aqueous solutions containing either guanidinium or thiocyanate ions. The first-order difference method of neutron diffraction and isotopic substitution was applied, and the hydration structures of two of nature's strongest denaturant ions were determined. Each ion is shown to interact weakly with water: Guanidinium has no recognizable hydration shell and is one of the most weakly hydrated cations yet characterized. Hydration of thiocyanate is characterized by a low coordination number involving around one hydrogen-bonded water molecule and approximately two water molecules weakly interacting through "hydration bonds." The weak hydration of these denaturant ions strongly supports suggestions that a major contribution to the denaturant effect is the preferential interaction of the denaturant with the protein surface. By contrast, solute species such as many sugars and related polyols that stabilize proteins are strongly hydrated and are thus preferentially retained in the bulk solvent and excluded from the protein surface.