Requirement of domain-domain interaction for conformational change and functional ATP hydrolysis in myosin

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.

Electrostatic interactions in the force-generating region of the human cardiac myosin modulate ADP dissociation from actomyosin

Human cardiac myosin has two isoforms, alpha and beta, sharing significant sequence similarity, but different in kinetics: ADP release from actomyosin is an order of magnitude faster in the alpha myosin isoform. Apparently, small differences in the sequence are responsible for distinct local inter-residue interactions within alpha and beta isoforms, leading to such a dramatic difference in the rate of ADP release. Our analysis of structural kinetics of alpha and beta isoforms using molecular dynamics simulations revealed distinct dynamics of SH1:SH2 helix within the force-generation region of myosin head. The simulations showed that the residue R694 of the helix forms two permanent salt bridges in the beta isoform, which are not present in the alpha isoform. We hypothesized that the isoform-specific electrostatic interactions play a role in the difference of kinetic properties of myosin isoforms. We prepared R694N mutant in the beta isoform background to destabilize electrostatic interactions in the force-generating region of the myosin head. Our experimental data confirm faster ADP release from R694N actomyosin mutant, but is not as dramatic as the difference of kinetics of ADP release in the alpha and beta isoforms.

Functional Diversity of Class XI Myosins in Arabidopsis thaliana

Plant myosin XI acts as a motive force for cytoplasmic streaming through interacting with actin filaments within the cell. Arabidopsis thaliana (At) has 13 genes belonging to the myosin XI family. Previous reverse genetic approaches suggest that At myosin XIs are partially redundant, but are functionally diverse for their specific tasks within the plant. However, the tissue-specific expression and enzymatic properties of myosin XIs have to date been poorly understood, primarily because of the difficulty in cloning and expressing large myosin XI genes and proteins. In this study, we cloned full-length cDNAs and promoter regions for all 13 At myosin XIs and identified tissue-specific expression (using promoter-reporter assays) and motile and enzymatic activities (using in vitro assays). In general, myosins belonging to the same class have similar velocities and ATPase activities. However, the velocities and ATPase activities of the 13 At myosin XIs are significantly different and are classified broadly into three groups based on velocity (high group, medium group and low group). Interestingly, the velocity groups appear roughly correlated with the tissue-specific expression patterns. Generally, ubiquitously expressed At myosin XIs belong to the medium-velocity group, pollen-specific At myosin XIs belong to the high-velocity group and only one At myosin XI (XI-I) is classified as belonging to the low-velocity group. In this study, we demonstrated the diversity of the 13 myosin XIs in Arabidopsis at the molecular and tissue levels. Our results indicate that myosin XIs in higher plants have distinct motile and enzymatic activities adapted for their specific roles.

Macromolecular Crowding Modulates Actomyosin Kinetics

Actomyosin kinetics is usually studied in dilute solutions, which do not reflect conditions in the cytoplasm. In cells, myosin and actin work in a dense macromolecular environment. High concentrations of macromolecules dramatically reduce the amount of free space available for all solutes, which results in an effective increase of the solutes' chemical potential and protein stabilization. Moreover, in a crowded solution, the chemical potential depends on the size of the solute, with larger molecules experiencing a larger excluded volume than smaller ones. Therefore, since myosin interacts with two ligands of different sizes (actin and ATP), macromolecular crowding can modulate the kinetics of individual steps of the actomyosin ATPase cycle. To emulate the effect of crowding in cells, we studied actomyosin cycle reactions in the presence of a high-molecular-weight polymer, Ficoll70. We observed an increase in the maximum velocity of the actomyosin ATPase cycle, and our transient-kinetics experiments showed that virtually all individual steps of the actomyosin cycle were affected by the addition of Ficoll70. The observed effects of macromolecular crowding on the myosin-ligand interaction cannot be explained by the increase of a solute's chemical potential. A time-resolved Förster resonance energy transfer experiment confirmed that the myosin head assumes a more compact conformation in the presence of Ficoll70 than in a dilute solution. We conclude that the crowding-induced myosin conformational change plays a major role in the changed kinetics of actomyosin ATPase.

Mechanochemical Energy Transduction during the Main Rotary Step in the Synthesis Cycle of F1-ATPase

F₁-ATPase is a highly efficient molecular motor that can synthesize ATP driven by a mechanical torque. Its ability to function reversibly in either direction requires tight mechanochemical coupling between the catalytic domain and the rotating central shaft, as well as temporal control of substrate binding and product release. Despite great efforts and significant progress, the molecular details of this synchronized and fine-tuned energy conversion mechanism are not fully understood. Here, we use extensive molecular dynamics simulations to reconcile recent single-molecule experiments with structural data and provide a consistent thermodynamic, kinetic and mechanistic description of the main rotary substep in the synthetic cycle of mammalian ATP synthase. The calculated free energy profiles capture a discrete pattern in the rotation of the central γ-shaft, with a metastable intermediate located-consistently with recent experimental findings-at 70° relative to the X-ray position. We identify this rotary step as the ATP-dependent substep, and find that the associated free energy input supports the mechanism involving concurrent nucleotide binding and release. During the main substep, our simulations show no significant opening of the ATP-bound β subunit; instead, we observe that mechanical energy is transmitted to its nucleotide binding site, thus lowering the affinity for ATP. Simultaneously, the empty subunit assumes a conformation that enables the enzyme to harness the free energy of ADP binding to drive ATP release. Finally, we show that ligand exchange is regulated by a checkpoint mechanism, an apparent prerequisite for high efficiency in protein nanomotors.

Myosin XI-I is Mechanically and Enzymatically Unique Among Class-XI Myosins in Arabidopsis

Arabidopsis possesses 13 genes encoding class-XI myosins. Among these, myosin XI-I is phylogenetically distant. To examine the molecular properties of Arabidopsis thaliana myosin XI-I (At myosin XI-I), we performed in vitro mechanical and enzymatic analyses using recombinant constructs of At myosin XI-I. Unlike other biochemically studied class-XI myosins, At myosin XI-I showed extremely low actin-activated ATPase activity (Vmax = 3.7 Pi s(-1) head(-1)). The actin-sliding velocity of At myosin XI-I was 0.25 µm s(-1), >10 times lower than those of other class-XI myosins. The ADP dissociation rate from acto-At myosin XI-I was 17 s(-1), accounting for the low actin-sliding velocity. In contrast, the apparent affinity for actin in the presence of ATP, estimated from Kapp (0.61 µM) of actin-activated ATPase, was extremely high. The equilibrium dissociation constant for actin was very low in both the presence and absence of ATP, indicating a high affinity for actin. To examine At myosin XI-I motility in vivo, green fluorescent protein-fused full-length At myosin XI-I was expressed in cultured Arabidopsis cells. At myosin XI-I localized not only on the nuclear envelope but also on small dots moving slowly (0.23 µm s(-1)) along actin filaments. Our results show that the properties of At myosin XI-I differ from those of other Arabidopsis class-XI myosins. The data suggest that At myosin XI-I does not function as a driving force for cytoplasmic streaming but regulates the organelle velocity, supports processive organelle movement or acts as a tension generator.

Molecular characterization and subcellular localization of Arabidopsis class VIII myosin, ATM1

Land plants possess myosin classes VIII and XI. Although some information is available on the molecular properties of class XI myosins, class VIII myosins are not characterized. Here, we report the first analysis of the enzymatic properties of class VIII myosin. The motor domain of Arabidopsis class VIII myosin, ATM1 (ATM1-MD), and the motor domain plus one IQ motif (ATM1-1IQ) were expressed in a baculovirus system and characterized. ATM1-MD and ATM1-1IQ had low actin-activated Mg(2+)-ATPase activity (Vmax = 4 s(-1)), although their affinities for actin were high (Kactin = 4 μM). The actin-sliding velocities of ATM1-MD and ATM1-1IQ were 0.02 and 0.089 μm/s, respectively, from which the value for full-length ATM1 is calculated to be ∼0.2 μm/s. The results of actin co-sedimentation assay showed that the duty ratio of ATM1 was ∼90%. ADP dissociation from the actin·ATM1 complex (acto-ATM1) was extremely slow, which accounts for the low actin-sliding velocity, low actin-activated ATPase activity, and high duty ratio. The rate of ADP dissociation from acto-ATM1 was markedly biphasic with fast and slow phase rates (5.1 and 0.41 s(-1), respectively). Physiological concentrations of free Mg(2+) modulated actin-sliding velocity and actin-activated ATPase activity by changing the rate of ADP dissociation from acto-ATM1. GFP-fused full-length ATM1 expressed in Arabidopsis was localized to plasmodesmata, plastids, newly formed cell walls, and actin filaments at the cell cortex. Our results suggest that ATM1 functions as a tension sensor/generator at the cell cortex and other structures in Arabidopsis.

Mode coupling points to functionally important residues in myosin II

Relevance of mode coupling to energy/information transfer during protein function, particularly in the context of allosteric interactions is widely accepted. However, existing evidence in favor of this hypothesis comes essentially from model systems. We here report a novel formal analysis of the near-native dynamics of myosin II, which allows us to explore the impact of the interaction between possibly non-Gaussian vibrational modes on fluctutational dynamics. We show that an information-theoretic measure based on mode coupling alone yields a ranking of residues with a statistically significant bias favoring the functionally critical locations identified by experiments on myosin II.

Metal cation controls myosin and actomyosin kinetics

We have perturbed myosin nucleotide binding site with magnesium-, manganese-, or calcium-nucleotide complexes, using metal cation as a probe to examine the pathways of myosin ATPase in the presence of actin. We have used transient time-resolved FRET, myosin intrinsic fluorescence, fluorescence of pyrene labeled actin, combined with the steady state myosin ATPase activity measurements of previously characterized D.discoideum myosin construct A639C:K498C. We found that actin activation of myosin ATPase does not depend on metal cation, regardless of the cation-specific kinetics of nucleotide binding and dissociation. The rate limiting step of myosin ATPase depends on the metal cation. The rate of the recovery stroke and the reverse recovery stroke is directly proportional to the ionic radius of the cation. The rate of nucleotide release from myosin and actomyosin, and ATP binding to actomyosin depends on the cation coordination number.

Allosteric activation transitions in enzymes and biomolecular motors: insights from atomistic and coarse-grained simulations

The chemical step in enzymes is usually preceded by a kinetically distinct activation step that involves large-scale conformational transitions. In "simple" enzymes this step corresponds to the closure of the active site; in more complex enzymes, such as biomolecular motors, the activation step is more complex and may involve interactions with other biomolecules. These activation transitions are essential to the function of enzymes and perturbations in the scale and/or rate of these transitions are implicated in various serious human diseases; incorporating key flexibilities into engineered enzymes is also considered a major remaining challenge in rational enzyme design. Therefore it is important to understand the underlying mechanism of these transitions. This is a significant challenge to both experimental and computational studies because of the allosteric and multi-scale nature of such transitions. Using our recent studies of two enzyme systems, myosin and adenylate kinase (AK), we discuss how atomistic and coarse-grained simulations can be used to provide insights into the mechanism of activation transitions in realistic systems. Collectively, the results suggest that although many allosteric transitions can be viewed as domain displacements mediated by flexible hinges, there are additional complexities and various deviations. For example, although our studies do not find any evidence for "cracking" in AK, our results do underline the contribution of intra-domain properties (e.g., dihedral flexibility) to the rate of the transition. The study of mechanochemical coupling in myosin highlights that local changes important to chemistry require stabilization from more extensive structural changes; in this sense, more global structural transitions are needed to activate the chemistry in the active site. These discussions further emphasize the importance of better understanding factors that control the degree of co-operativity for allosteric transitions, again hinting at the intimate connection between protein stability and functional flexibility. Finally, a number of topics of considerable future interest are briefly discussed.

Alternative relay and converter domains tune native muscle myosin isoform function in Drosophila

Myosin isoforms help define muscle-specific contractile and structural properties. Alternative splicing of myosin heavy chain gene transcripts in Drosophila melanogaster yields muscle-specific isoforms and highlights alternative domains that fine-tune myosin function. To gain insight into how native myosin is tuned, we expressed three embryonic myosin isoforms in indirect flight muscles lacking endogenous myosin. These isoforms differ in their relay and/or converter domains. We analyzed isoform-specific ATPase activities, in vitro actin motility and myofibril structure/stability. We find that dorsal acute body wall muscle myosin (EMB-9c11d) shows a significant increase in MgATPase V(max) and actin sliding velocity, as well as abnormal myofibril assembly compared to cardioblast myosin (EMB-11d). These properties differ as a result of alternative exon-9-encoded relay domains that are hypothesized to communicate signals among the ATP-binding pocket, actin-binding site and the converter domain. Further, EMB-11d shows significantly reduced levels of basal Ca- and MgATPase as well as MgATPase V(max) compared to embryonic body wall muscle isoform (EMB) (expressed in a multitude of body wall muscles). EMB-11d also induces increased actin sliding velocity and stabilizes myofibril structure compared to EMB. These differences arise from exon-11-encoded alternative converter domains that are proposed to reposition the lever arm during the power and recovery strokes. We conclude that relay and converter domains of native myosin isoforms fine-tune ATPase activity, actin motility and muscle ultrastructure. This verifies and extends previous studies with chimeric molecules and indicates that interactions of the relay and converter during the contractile cycle are key to myosin-isoform-specific kinetic and mechanical functions.

Unique charge distribution in surface loops confers high velocity on the fast motor protein Chara myosin

Most myosins have a positively charged loop 2 with a cluster of lysine residues that bind to the negatively charged N-terminal segment of actin. However, the net charge of loop 2 of very fast Chara myosin is zero and there is no lysine cluster in it. In contrast, Chara myosin has a highly positively charged loop 3. To elucidate the role of these unique surface loops of Chara myosin in its high velocity and high actin-activated ATPase activity, we have undertaken mutational analysis using recombinant Chara myosin motor domain. It was found that net positive charge in loop 3 affected V(max) and K(app) of actin activated ATPase activity, while it affected the velocity only slightly. The net positive charge in loop 2 affected K(app) and the velocity, although it did not affect V(max). Our results suggested that Chara myosin has evolved to have highly positively charged loop 3 for its high ATPase activity and have less positively charged loop 2 for its high velocity. Since high positive charge in loop 3 and low positive charge in loop 2 seem to be one of the reasons for Chara myosin's high velocity, we manipulated charge contents in loops 2 and 3 of Dictyostelium myosin (class II). Removing positive charge from loop 2 and adding positive charge to loop 3 of Dictyostelium myosin made its velocity higher than that of the wild type, suggesting that the charge strategy in loops 2 and 3 is widely applicable.

Myosin-catalyzed ATP hydrolysis elucidated by 31P NMR kinetic studies and 1H PFG-diffusion measurements

We conducted (31)P NMR kinetic studies and (1)H-diffusion measurements on myosin-catalyzed hydrolysis of adenosine triphosphate (ATP) under varied conditions. The data elucidate well the overall hydrolysis rate and various factors that significantly impact the reaction. We found that the enzymatic hydrolysis of ATP to adenosine diphosphate (ADP) was followed by ADP hydrolysis, and different nucleotides such as ADP and guanosine triphosphate acted as competitors of ATP. Increasing ATP or Mg(2+) concentration resulted in decreased hydrolysis rate, and such effect can be related to the decrease of ATP diffusion constants. Below 50 degrees C, the hydrolysis was accelerated by increasing temperature following the Arrhenius' law, but the hydrolysis rate was significantly lowered at higher temperature (approximately 60 degrees C), due to the thermal-denaturation of myosin. The optimal pH range was around pH 6-8. These results are important for characterization of myosin-catalyzed ATP hydrolysis, and the method is also applicable to other enzymatic nucleotide reactions.

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.

Alternative relay domains of Drosophila melanogaster myosin differentially affect ATPase activity, in vitro motility, myofibril structure and muscle function

The relay domain of myosin is hypothesized to function as a communication pathway between the nucleotide-binding site, actin-binding site and the converter domain. In Drosophila melanogaster, a single myosin heavy chain gene encodes three alternative relay domains. Exon 9a encodes the indirect flight muscle isoform (IFI) relay domain, whereas exon 9b encodes one of the embryonic body wall isoform (EMB) relay domains. To gain a better understanding of the function of the relay domain and the differences imparted by the IFI and the EMB versions, we constructed two transgenic Drosophila lines expressing chimeric myosin heavy chains in indirect flight muscles lacking endogenous myosin. One expresses the IFI relay domain in the EMB backbone (EMB-9a), while the second expresses the EMB relay domain in the IFI backbone (IFI-9b). Our studies reveal that the EMB relay domain is functionally equivalent to the IFI relay domain when it is substituted into IFI. Essentially no differences in ATPase activity, actin-sliding velocity, flight ability at room temperature or muscle structure are observed in IFI-9b compared to native IFI. However, when the EMB relay domain is replaced with the IFI relay domain, we find a 50% reduction in actin-activated ATPase activity, a significant increase in actin affinity, abolition of actin sliding, defects in myofibril assembly and rapid degeneration of muscle structure compared to EMB. We hypothesize that altered relay domain conformational changes in EMB-9a impair intramolecular communication with the EMB-specific converter domain. This decreases transition rates involving strongly bound actomyosin states, leading to a reduced ATPase rate and loss of actin motility.

Predicting allosteric communication in myosin via a pathway of conserved residues

We present a computational method that predicts a pathway of residues that mediate protein allosteric communication. The pathway is predicted using only a combination of distance constraints between contiguous residues and evolutionary data. We applied this analysis to find pathways of conserved residues connecting the myosin ATP binding site to the lever arm. These pathway residues may mediate the allosteric communication that couples ATP hydrolysis to the lever arm recovery stroke. Having examined pre-stroke conformations of Dictyostelium, scallop, and chicken myosin II as well as Dictyostelium myosin I, we observed a conserved pathway traversing switch II and the relay helix, which is consistent with the understood need for allosteric communication in this conformation. We also examined post-rigor and rigor conformations across several myosin species. Although initial residues of these paths are more heterogeneous, all but one of these paths traverse a consistent set of relay helix residues to reach the beginning of the lever arm. We discuss our results in the context of structural elements and reported mutational experiments, which substantiate the significance of the pre-stroke pathways. Our method provides a simple, computationally efficient means of predicting a set of residues that mediate allosteric communication. We provide a refined, downloadable application and source code (on https://simtk.org) to share this tool with the wider community (https://simtk.org/home/allopathfinder).

Kinetic mechanism of the fastest motor protein, Chara myosin

Chara corallina class XI myosin is by far the fastest molecular motor. To investigate the molecular mechanism of this fast movement, we performed a kinetic analysis of a recombinant motor domain of Chara myosin. We estimated the time spent in the strongly bound state with actin by measuring rate constants of ADP dissociation from actin.motor domain complex and ATP-induced dissociation of the motor domain from actin. The rate constant of ADP dissociation from acto-motor domain was >2800 s(-1), and the rate constant of ATP-induced dissociation of the motor domain from actin at physiological ATP concentration was 2200 s(-1). From these data, the time spent in the strongly bound state with actin was estimated to be <0.82 ms. This value is the shortest among known values for various myosins and yields the duty ratio of <0.3 with a V(max) value of the actin-activated ATPase activity of 390 s(-1). The addition of the long neck domain of myosin Va to the Chara motor domain largely increased the velocity of the motility without increasing the ATP hydrolysis cycle rate, consistent with the swinging lever model. In addition, this study reveals some striking kinetic features of Chara myosin that are suited for the fast movement: a dramatic acceleration of ADP release by actin (1000-fold) and extremely fast ATP binding rate.

Mechanochemical coupling in the myosin motor domain. II. Analysis of critical residues

An important challenge in the analysis of mechanochemical coupling in molecular motors is to identify residues that dictate the tight coupling between the chemical site and distant structural rearrangements. In this work, a systematic attempt is made to tackle this issue for the conventional myosin. By judiciously combining a range of computational techniques with different approximations and strength, which include targeted molecular dynamics, normal mode analysis, and statistical coupling analysis, we are able to identify a set of important residues and propose their relevant function during the recovery stroke of myosin. These analyses also allowed us to make connections with previous experimental and computational studies in a critical manner. The behavior of the widely used reporter residue, Trp501, in the simulations confirms the concern that its fluorescence does not simply reflect the relay loop conformation or active-site open/close but depends subtly on its microenvironment. The findings in the targeted molecular dynamics and a previous minimum energy path analysis of the recovery stroke have been compared and analyzed, which emphasized the difference and complementarity of the two approaches. In conjunction with our previous studies, the current set of investigations suggest that the modulation of structural flexibility at both the local (e.g., active-site) and domain scales with strategically placed "hotspot" residues and phosphate chemistry is likely the general feature for mechanochemical coupling in many molecular motors. The fundamental strategies of examining both collective and local changes and combining physically motivated methods and informatics-driven techniques are expected to be valuable to the study of other molecular motors and allosteric systems in general.

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).

A point mutation in the SH1 helix alters elasticity and thermal stability of myosin II

Movement generated by the myosin motor is generally thought to be driven by distortion of an elastic element within the myosin molecule and subsequent release of the resulting strain. However, the location of this elastic element in myosin remains unclear. The myosin motor domain consists of four major subdomains connected by flexible joints. The SH1 helix is the joint that connects the converter subdomain to the other domains, and is thought to play an important role in arrangements of the converter relative to the motor. To investigate the involvement of the SH1 helix in elastic distortion in myosin, we have introduced a point mutation into the SH1 helix of Dictyostelium myosin II (R689H), which in human nonmuscle myosin IIA causes nonsyndromic hereditary deafness, DFNA17. The mutation resulted in a significant impairment in motile activities, whereas actin-activated ATPase activity was only slightly affected. Single molecule mechanical measurements using optical trap showed that the step size was not shortened by the mutation, suggesting that the slower motility is caused by altered kinetics. The single molecule measurements demonstrated that the mutation significantly reduced cross-bridge stiffness. Motile activities produced by mixtures of wild-type and mutant myosins also suggested that the mutation affected the elasticity of myosin. These results suggest that the SH1 helix is involved in modulation of myosin elasticity, presumably by modulating the converter flexibility. Consistent with this, the mutation was also shown to reduce thermal stability and induce thermal aggregation of the protein, which might be implicated in the disease process.

Reconciling the "old" and "new" views of protein allostery: a molecular simulation study of chemotaxis Y protein (CheY)

A combination of thirty-two 10-ns-scale molecular dynamics simulations were used to explore the coupling between conformational transition and phosphorylation in the bacteria chemotaxis Y protein (CheY), as a simple but representative example of protein allostery. Results from these simulations support an activation mechanism in which the beta4-alpha4 loop, at least partially, gates the isomerization of Tyr106. The roles of phosphorylation and the conserved Thr87 are deemed indirect in that they stabilize the active configuration of the beta4-alpha4 loop. The indirect role of the activation event (phosphorylation) and/or conserved residues in stabilizing, rather than causing, specific conformational transition is likely a feature in many signaling systems. The current analysis of CheY also helps to make clear that neither the "old" (induced fit) nor the "new" (population shift) views for protein allostery are complete, because they emphasize the kinetic (mechanistic) and thermodynamic aspects of allosteric transitions, respectively. In this regard, an issue that warrants further analysis concerns the interplay of concerted collective motion and sequential local structural changes in modulating cooperativity between distant sites in biomolecules.

Dynamic protein domains: identification, interdependence, and stability

Existing methods of domain identification in proteins usually provide no information about the degree of domain independence and stability. However, this information is vital for many areas of protein research. The recently developed hierarchical clustering of correlation patterns (HCCP) technique provides machine-based domain identification in a computationally simple and physically consistent way. Here we present the modification of this technique, which not only allows determination of the most plausible number of dynamic domains but also makes it possible to estimate the degree of their independence (the extent of correlated motion) and stability (the range of environmental conditions, where domains remain intact). With this technique we provided domain assignments and calculated intra- and interdomain correlations and interdomain energies for >2500 test proteins. It is shown that mean intradomain correlation of motions can serve as a quantitative criterion of domain independence, and the HCCP stability gap is a measure of their stability. Our data show that the motions of domains with high stability are usually independent. In contrast, the domains with moderate stability usually exhibit a substantial degree of correlated motions. It is shown that in multidomain proteins the domains are most stable if they are of similar size, and this correlates with the observed abundance of such proteins.

Comparison of mode analyses at different resolutions applied to nucleic acid systems

More than two decades of different types of mode analyses has shown that these techniques can be useful in describing large-scale motions in protein systems. A number of mode analyses are available and include quasiharmonics, classical normal mode, block normal mode, and the elastic network model. Each of these methods has been validated for protein systems and this variety allows researchers to choose the technique that gives the best compromise between computational cost and the level of detail in the calculation. These same techniques have not been systematically tested for nucleic acid systems, however. Given the differences in interactions and structural features between nucleic acid and protein systems, the validity of these techniques in the protein regime cannot be directly translated into validity in the nucleic acid realm. In this work, we investigate the usefulness of the above mode analyses as applied to two RNA systems, i.e., the hammerhead ribozyme and a guanine riboswitch. We show that classical normal-mode analysis can match the magnitude and direction of residue fluctuations from the more detailed, anharmonic technique, quasiharmonic analysis of a molecular dynamics trajectory. The block normal-mode approximation is shown to hold in the nucleic acid systems studied. Only the mode analysis at the lowest level of detail, the elastic network model, produced mixed results in our calculations. We present data that suggest that the elastic network model, with the popular parameterization, is not best suited for systems that do not have a close packed structure; this observation also hints at why the elastic network model has been found to be valid for many globular protein systems. The different behaviors of block normal-mode analysis and the elastic network model, which invoke similar degrees of coarse-graining to the dynamics but use different potentials, suggest the importance of applying a heterogeneous potential function in a robust analysis of the dynamics of biomolecules, especially those that are not closely packed. In addition to these comparisons, we briefly discuss insights into the conformational space available to the hammerhead ribozyme.

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.