The Structural Coupling between ATPase Activation and Recovery Stroke in the Myosin II Motor

UNC-45 assisted myosin folding depends on a conserved FX3HY motif implicated in Freeman Sheldon Syndrome

Myosin motors are critical for diverse motility functions, ranging from cytokinesis and endocytosis to muscle contraction. The UNC-45 chaperone controls myosin function mediating the folding, assembly, and degradation of the muscle protein. Here, we analyze the molecular mechanism of UNC-45 as a hub in myosin quality control. We show that UNC-45 forms discrete complexes with folded and unfolded myosin, forwarding them to downstream chaperones and E3 ligases. Structural analysis of a minimal chaperone:substrate complex reveals that UNC-45 binds to a conserved FX3HY motif in the myosin motor domain. Disrupting the observed interface by mutagenesis prevents myosin maturation leading to protein aggregation in vivo. We also show that a mutation in the FX3HY motif linked to the Freeman Sheldon Syndrome impairs UNC-45 assisted folding, reducing the level of functional myosin. These findings demonstrate that a faulty myosin quality control is a critical yet unexplored cause of human myopathies.

A weak coupling mechanism for the early steps of the recovery stroke of myosin VI: A free energy simulation and string method analysis

Myosin motors use the energy of ATP to produce force and directed movement on actin by a swing of the lever-arm. ATP is hydrolysed during the off-actin re-priming transition termed recovery stroke. To provide an understanding of chemo-mechanical transduction by myosin, it is critical to determine how the reverse swing of the lever-arm and ATP hydrolysis are coupled. Previous studies concluded that the recovery stroke of myosin II is initiated by closure of the Switch II loop in the nucleotide-binding site. Recently, we proposed that the recovery stroke of myosin VI starts with the spontaneous re-priming of the converter domain to a putative pre-transition state (PTS) intermediate that precedes Switch II closing and ATPase activation. Here, we investigate the transition from the pre-recovery, post-rigor (PR) state to PTS in myosin VI using geometric free energy simulations and the string method. First, our calculations rediscover the PTS state agnostically and show that it is accessible from PR via a low free energy transition path. Second, separate path calculations using the string method illuminate the mechanism of the PR to PTS transition with atomic resolution. In this mechanism, the initiating event is a large movement of the converter/lever-arm region that triggers rearrangements in the Relay-SH1 region and the formation of the kink in the Relay helix with no coupling to the active site. Analysis of the free-energy barriers along the path suggests that the converter-initiated mechanism is much faster than the one initiated by Switch II closure, which supports the biological relevance of PTS as a major on-pathway intermediate of the recovery stroke in myosin VI. Our analysis suggests that lever-arm re-priming and ATP hydrolysis are only weakly coupled, so that the myosin recovery stroke is initiated by thermal fluctuations and stabilised by nucleotide consumption via a ratchet-like mechanism.

Binding pocket dynamics along the recovery stroke of human β-cardiac myosin

The druggability of small-molecule binding sites can be significantly affected by protein motions and conformational changes. Ligand binding, protein dynamics and protein function have been shown to be closely interconnected in myosins. The breakthrough discovery of omecamtiv mecarbil (OM) has led to an increased interest in small molecules that can target myosin and modulate its function for therapeutic purposes (myosin modulators). In this work, we use a combination of computational methods, including steered molecular dynamics, umbrella sampling and binding pocket tracking tools, to follow the evolution of the OM binding site during the recovery stroke transition of human β-cardiac myosin. We found that steering two internal coordinates of the motor domain can recapture the main features of the transition and in particular the rearrangements of the binding site, which shows significant changes in size, shape and composition. Possible intermediate conformations were also identified, in remarkable agreement with experimental findings. The differences in the binding site properties observed along the transition can be exploited for the future development of conformation-selective myosin modulators.

Insights into the Mechanism of the Cardiac Drug Omecamtiv Mecarbil─A Computational Study

Omecamtiv mecarbil (OM) is a positive inotrope that is thought to bind directly to an allosteric site of the β-cardiac myosin. The drug is under investigation for the treatment of systolic heart failure. The drug is classified as a cardiac myosin modulator and has been observed to affect multiple vital steps of the cross-bridge cycle including the recovery stroke and the chemical step. We explored the free-energy surface of the recovery stroke of the human cardiac β-myosin in the presence of OM to determine its influence on this process. We also investigated the effects of OM on the recovery stroke in the presence of genetic cardiomyopathic mutations R712L, F764L, and P710R using metadynamics. We also utilized the method of transition path sampling to generate an unbiased ensemble of reactive trajectories for the ATP hydrolysis step in the presence of OM that were able to provide insight into the differences observed due to OM in the dynamics and mechanism of the decomposition of ATP to ADP and HPO₄^2-, a central part of the power generation in cardiac muscle. We studied chemistry in the presence of the same three mutations to further elucidate the effect of OM, and its use in the treatment of cardiac disease.

Binding Networks Identify Targetable Protein Pockets for Mechanism-Based Drug Design

The human genome codes only a few thousand druggable proteins, mainly receptors and enzymes. While this pool of available drug targets is limited, there is an untapped potential for discovering new drug-binding mechanisms and modes. For example, enzymes with long binding cavities offer numerous prerequisite binding sites that may be visited by an inhibitor during migration from a bulk solution to the destination site. Drug design can use these prerequisite sites as new structural targets. However, identifying these ephemeral sites is challenging. Here, we introduce a new method called NetBinder for the systematic identification and classification of prerequisite binding sites at atomic resolution. NetBinder is based on atomistic simulations of the full inhibitor binding process and provides a networking framework on which to select the most important binding modes and uncover the entire binding mechanism, including previously undiscovered events. NetBinder was validated by a study of the binding mechanism of blebbistatin (a potent inhibitor) to myosin 2 (a promising target for cancer chemotherapy). Myosin 2 is a good test enzyme because, like other potential targets, it has a long internal binding cavity that provides blebbistatin with numerous potential prerequisite binding sites. The mechanism proposed by NetBinder of myosin 2 structural changes during blebbistatin binding shows excellent agreement with experimentally determined binding sites and structural changes. While NetBinder was tested on myosin 2, it may easily be adopted to other proteins with long internal cavities, such as G-protein-coupled receptors or ion channels, the most popular current drug targets. NetBinder provides a new paradigm for drug design by a network-based elucidation of binding mechanisms at an atomic resolution.

Actomyosin Complex

Formation of cross-bridges between actin and myosin occurs ubiquitously in eukaryotic cells and mediates muscle contraction, intracellular cargo transport, and cytoskeletal remodeling. Myosin motors repeatedly bind to and dissociate from actin filaments in a cycle that transduces the chemical energy from ATP hydrolysis into mechanical force generation. While the general layout of surface elements within the actin-binding interface is conserved among myosin classes, sequence divergence within these motifs alters the specific contacts involved in the actomyosin interaction as well as the kinetics of mechanochemical cycle phases. Additionally, diverse lever arm structures influence the motility and force production of myosin molecules during their actin interactions. The structural differences generated by myosin's molecular evolution have fine-tuned the kinetics of its isoforms and adapted them for their individual cellular roles. In this chapter, we will characterize the structural and biochemical basis of the actin-myosin interaction and explain its relationship with myosin's cellular roles, with emphasis on the structural variation among myosin isoforms that enables their functional specialization. We will also discuss the impact of accessory proteins, such as the troponin-tropomyosin complex and myosin-binding protein C, on the formation and regulation of actomyosin cross-bridges.

Investigation of the Recovery Stroke and ATP Hydrolysis and Changes Caused Due to the Cardiomyopathic Point Mutations in Human Cardiac β Myosin

Human cardiac β myosin undergoes the cross-bridge cycle as part of the force-generating mechanism of cardiac muscle. The recovery stroke is considered one of the key steps of the kinetic cycle as it is the conformational rearrangement required to position the active site residues for hydrolysis of ATP and interaction with actin. We explored the free-energy surface of the transition and investigated the effect of the genetic cardiomyopathy causing mutations R453C, I457T, and I467T on this step using metadynamics. This work extends previous studies on Dictyostelium myosin II with engineered mutations. Here, like previously, we generated an unbiased thermodynamic ensemble of reactive trajectories for the chemical step using transition path sampling. Our methodologies were able to predict the changes to the dynamics of the recovery stroke as well as predict the pathway of breakdown of ATP to ADP and HPO₄^2- with the stabilization of the metaphosphate intermediate. We also observed clear differences between the Dictyostelium myosin II and human cardiac β myosin for ATP hydrolysis as well as predict the effect of the mutation I467T on the chemical step.

Myosin dilated cardiomyopathy mutation S532P disrupts actomyosin interactions, leading to altered muscle kinetics, reduced locomotion, and cardiac dilation in Drosophila

Dilated cardiomyopathy (DCM), a life-threatening disease characterized by pathological heart enlargement, can be caused by myosin mutations that reduce contractile function. To better define the mechanistic basis of this disease, we employed the powerful genetic and integrative approaches available in Drosophila melanogaster. To this end, we generated and analyzed the first fly model of human myosin-induced DCM. The model reproduces the S532P human β-cardiac myosin heavy chain DCM mutation, which is located within an actin-binding region of the motor domain. In concordance with the mutation's location at the actomyosin interface, steady-state ATPase and muscle mechanics experiments revealed that the S532P mutation reduces the rates of actin-dependent ATPase activity and actin binding and increases the rate of actin detachment. The depressed function of this myosin form reduces the number of cross-bridges during active wing beating, the power output of indirect flight muscles, and flight ability. Further, S532P mutant hearts exhibit cardiac dilation that is mutant gene dose-dependent. Our study shows that Drosophila can faithfully model various aspects of human DCM phenotypes and suggests that impaired actomyosin interactions in S532P myosin induce contractile deficits that trigger the disease.

Electrostatic interactions in the SH1-SH2 helix of human cardiac myosin modulate the time of strong actomyosin binding

Two single mutations, R694N and E45Q, were introduced in the beta isoform of human cardiac myosin to remove permanent salt bridges E45:R694 and E98:R694 in the SH1-SH2 helix of the myosin head. Beta isoform-specific bridges E45:R694 and E98:R694 were discovered in the molecular dynamics simulations of the alpha and beta myosin isoforms. Alpha and beta isoforms exhibit different kinetics, ADP dissociates slower from actomyosin containing beta myosin isoform, therefore, beta myosin stays strongly bound to actin longer. We hypothesize that the electrostatic interactions in the SH1-SH2 helix modulate the affinity of ADP to actomyosin, and therefore, the time of the strong actomyosin binding. Wild type and the mutants of the myosin head construct (1-843 amino acid residues) were expressed in differentiated C2C12 cells, and the duration of the strongly bound state of actomyosin was characterized using transient kinetics spectrophotometry. All myosin constructs exhibited a fast rate of ATP binding to actomyosin and a slow rate of ADP dissociation, showing that ADP release limits the time of the strongly bound state of actomyosin. The mutant R694N showed a faster rate of ADP release from actomyosin, compared to the wild type and the E45Q mutant, thus indicating that electrostatic interactions within the SH1-SH2 helix region of human cardiac myosin modulate ADP release and thus, the duration of the strongly bound state of actomyosin.

Switch-1 instability at the active site decouples ATP hydrolysis from force generation in myosin II

Myosin active site elements (i.e., switch‐1) bind both ATP and a divalent metal to coordinate ATP hydrolysis. ATP hydrolysis at the active site is linked via allosteric communication to the actin polymer binding site and lever arm movement, thus coupling the free energy of ATP hydrolysis to force generation. How active site motifs are functionally linked to actin binding and the power stroke is still poorly understood. We hypothesize that destabilizing switch‐1 movement at the active site will negatively affect the tight coupling of the ATPase catalytic cycle to force production. Using a metal‐switch system, we tested the effect of interfering with switch‐1 coordination of the divalent metal cofactor on force generation. We found that while ATPase activity increased, motility was inhibited. Our results demonstrate that a single atom change that affects the switch‐1 interaction with the divalent metal directly affects actin binding and productive force generation. Even slight modification of the switch‐1 divalent metal coordination can decouple ATP hydrolysis from motility. Switch‐1 movement is therefore critical for both structural communication with the actin binding site, as well as coupling the energy of ATP hydrolysis to force generation.

Mechanochemical Function of Myosin II: Investigation into the Recovery Stroke and ATP Hydrolysis

Myosin regulates muscle function through a complex cycle of conformational rearrangements coupled with the hydrolysis of adenosine triphosphate (ATP). The recovery stroke reorganizes the myosin active site to hydrolyze ATP and cross bridge with the thin filament to produce muscle contraction. Engineered mutations K84M and R704E in Dictyostelium myosin have been designed to specifically inhibit the recovery stroke and have been shown to indirectly affect the ATPase activity of myosin. We investigated these mutagenic perturbations to the recovery stroke and generated thermodynamically correct and unbiased trajectories for native ATP hydrolysis with computationally enhanced sampling methods. Our methodology was able to resolve experimentally observed changes to kinetic and equilibrium dynamics for the recovery stroke with the correct prediction in the severity of these changes. For ATP hydrolysis, the sequential nature along with the stabilization of a metaphosphate intermediate was observed in agreement with previous studies. However, we observed glutamate 459 being utilized as a proton abstractor to prime the attacking water instead of a lytic water, a phenomenon not well categorized in myosin but has in other ATPases. Both rare event methodologies can be extended to human myosin to investigate isoformic differences from Dictyostelium and scan cardiomyopathic mutations to see differential perturbations to kinetics of other conformational changes in myosin such as the power stroke.

Structural and Computational Insights into a Blebbistatin-Bound Myosin•ADP Complex with Characteristics of an ADP-Release Conformation along the Two-Step Myosin Power Stoke

The motor protein myosin drives a wide range of cellular and muscular functions by generating directed movement and force, fueled through adenosine triphosphate (ATP) hydrolysis. Release of the hydrolysis product adenosine diphosphate (ADP) is a fundamental and regulatory process during force production. However, details about the molecular mechanism accompanying ADP release are scarce due to the lack of representative structures. Here we solved a novel blebbistatin-bound myosin conformation with critical structural elements in positions between the myosin pre-power stroke and rigor states. ADP in this structure is repositioned towards the surface by the phosphate-sensing P-loop, and stabilized in a partially unbound conformation via a salt-bridge between Arg131 and Glu187. A 5 Å rotation separates the mechanical converter in this conformation from the rigor position. The crystallized myosin structure thus resembles a conformation towards the end of the two-step power stroke, associated with ADP release. Computationally reconstructing ADP release from myosin by means of molecular dynamics simulations further supported the existence of an equivalent conformation along the power stroke that shows the same major characteristics in the myosin motor domain as the resolved blebbistatin-bound myosin-II·ADP crystal structure, and identified a communication hub centered on Arg232 that mediates chemomechanical energy transduction.

Plasmodium myosin A drives parasite invasion by an atypical force generating mechanism

Plasmodium parasites are obligate intracellular protozoa and causative agents of malaria, responsible for half a million deaths each year. The lifecycle progression of the parasite is reliant on cell motility, a process driven by myosin A, an unconventional single-headed class XIV molecular motor. Here we demonstrate that myosin A from Plasmodium falciparum (PfMyoA) is critical for red blood cell invasion. Further, using a combination of X-ray crystallography, kinetics, and in vitro motility assays, we elucidate the non-canonical interactions that drive this motor’s function. We show that PfMyoA motor properties are tuned by heavy chain phosphorylation (Ser19), with unphosphorylated PfMyoA exhibiting enhanced ensemble force generation at the expense of speed. Regulated phosphorylation may therefore optimize PfMyoA for enhanced force generation during parasite invasion or for fast motility during dissemination. The three PfMyoA crystallographic structures presented here provide a blueprint for discovery of specific inhibitors designed to prevent parasite infection.

Here, Robert-Paganin et al. show that myosin A from Plasmodium falciparum is critical for red blood cell invasion and that non-canonical interactions and regulated phosphorylation are important for force generation during parasite invasion.

Molecular parallelism in fast-twitch muscle proteins in echolocating mammals

Detecting associations between genomic changes and phenotypic differences is fundamental to understanding how phenotypes evolved. By systematically screening for parallel amino acid substitutions, we detected known as well as novel cases (Strc, Tecta, and Cabp2) of parallelism between echolocating bats and toothed whales in proteins that could contribute to high-frequency hearing adaptations. Our screen also showed that echolocating mammals exhibit an unusually high number of parallel substitutions in fast-twitch muscle fiber proteins. Both echolocating bats and toothed whales produce an extremely rapid call rate when homing in on their prey, which was shown in bats to be powered by specialized superfast muscles. We show that these genes with parallel substitutions (Casq1, Atp2a1, Myh2, and Myl1) are expressed in the superfast sound-producing muscle of bats. Furthermore, we found that the calcium storage protein calsequestrin 1 of the little brown bat and the bottlenose dolphin functionally converged in its ability to form calcium-sequestering polymers at lower calcium concentrations, which may contribute to rapid calcium transients required for superfast muscle physiology. The proteins that our genomic screen detected could be involved in the convergent evolution of vocalization in echolocating mammals by potentially contributing to both rapid Ca²⁺ transients and increased shortening velocities in superfast muscles.

Long-range coupling between ATP-binding and lever-arm regions in myosin via dielectric allostery

A protein molecule is a dielectric substance, so the binding of a ligand is expected to induce dielectric response in the protein molecule, considering that ligands are charged or polar in general. We previously reported that binding of adenosine triphosphate (ATP) to molecular motor myosin actually induces such a dielectric response in myosin due to the net negative charge of ATP. By this dielectric response, referred to as "dielectric allostery," spatially separated two regions in myosin, the ATP-binding region and the actin-binding region, are allosterically coupled. In this study, from the statistically stringent analyses of the extensive molecular dynamics simulation data obtained in the ATP-free and the ATP-bound states, we show that there exists the dielectric allostery that transmits the signal of ATP binding toward the distant lever-arm region. The ATP-binding-induced electrostatic potential change observed on the surface of the main domain induced a movement of the converter subdomain from which the lever arm extends. The dielectric response was found to be caused by an underlying large-scale concerted rearrangement of the electrostatic bond network, in which highly conserved charged/polar residues are involved. Our study suggests the importance of the dielectric property for molecular machines in exerting their function.

An intermediate along the recovery stroke of myosin VI revealed by X-ray crystallography and molecular dynamics

Myosins form a class of actin-based, ATPase motor proteins that mediate important cellular functions such as cargo transport and cell motility. Their functional cycle involves two large-scale swings of the lever arm: the force-generating powerstroke, which takes place on actin, and the recovery stroke during which the lever arm is reprimed into an armed configuration. Previous analyses of the prerecovery (postrigor) and postrecovery (prepowerstroke) states predicted that closure of switch II in the ATP binding site precedes the movement of the converter and the lever arm. Here, we report on a crystal structure of myosin VI, called pretransition state (PTS), which was solved at 2.2 Å resolution. Structural analysis and all-atom molecular dynamics simulations are consistent with PTS being an intermediate along the recovery stroke, where the Relay/SH1 elements adopt a postrecovery conformation, and switch II remains open. In this state, the converter appears to be largely uncoupled from the motor domain and explores an ensemble of partially reprimed configurations through extensive, reversible fluctuations. Moreover, we found that the free energy cost of hydrogen-bonding switch II to ATP is lowered by more than 10 kcal/mol compared with the prerecovery state. These results support the conclusion that closing of switch II does not initiate the recovery stroke transition in myosin VI. Rather, they suggest a mechanism in which lever arm repriming would be mostly driven by thermal fluctuations and eventually stabilized by the switch II interaction with the nucleotide in a ratchet-like fashion.

Mechanistic insights into the active site and allosteric communication pathways in human nonmuscle myosin-2C

Despite a generic, highly conserved motor domain, ATP turnover kinetics and their activation by F-actin vary greatly between myosin-2 isoforms. Here, we present a 2.25 Å pre-powerstroke state (ADP⋅VO₄) crystal structure of the human nonmuscle myosin-2C motor domain, one of the slowest myosins characterized. In combination with integrated mutagenesis, ensemble-solution kinetics, and molecular dynamics simulation approaches, the structure reveals an allosteric communication pathway that connects the distal end of the motor domain with the active site. Disruption of this pathway by mutation of hub residue R788, which forms the center of a cluster of interactions connecting the converter, the SH1-SH2 helix, the relay helix, and the lever, abolishes nonmuscle myosin-2 specific kinetic signatures. Our results provide insights into structural changes in the myosin motor domain that are triggered upon F-actin binding and contribute critically to the mechanochemical behavior of stress fibers, actin arcs, and cortical actin-based structures.

Allosteric modulation of cardiac myosin dynamics by omecamtiv mecarbil

New promising avenues for the pharmacological treatment of skeletal and heart muscle diseases rely on direct sarcomeric modulators, which are molecules that can directly bind to sarcomeric proteins and either inhibit or enhance their activity. A recent breakthrough has been the discovery of the myosin activator omecamtiv mecarbil (OM), which has been shown to increase the power output of the cardiac muscle and is currently in clinical trials for the treatment of heart failure. While the overall effect of OM on the mechano-chemical cycle of myosin is to increase the fraction of myosin molecules in the sarcomere that are strongly bound to actin, the molecular basis of its action is still not completely clear. We present here a Molecular Dynamics study of the motor domain of human cardiac myosin bound to OM, where the effects of the drug on the dynamical properties of the protein are investigated for the first time with atomistic resolution. We found that OM has a double effect on myosin dynamics, inducing a) an increased coupling of the motions of the converter and lever arm subdomains to the rest of the protein and b) a rewiring of the network of dynamic correlations, which produces preferential communication pathways between the OM binding site and distant functional regions. The location of the residues responsible for these effects suggests possible strategies for the future development of improved drugs and the targeting of specific cardiomyopathy-related mutations.

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.

Highly selective inhibition of myosin motors provides the basis of potential therapeutic application

Direct inhibition of smooth muscle myosin (SMM) is a potential means to treat hypercontractile smooth muscle diseases. The selective inhibitor CK-2018571 prevents strong binding to actin and promotes muscle relaxation in vitro and in vivo. The crystal structure of the SMM/drug complex reveals that CK-2018571 binds to a novel allosteric pocket that opens up during the "recovery stroke" transition necessary to reprime the motor. Trapped in an intermediate of this fast transition, SMM is inhibited with high selectivity compared with skeletal muscle myosin (IC₅₀ = 9 nM and 11,300 nM, respectively), although all of the binding site residues are identical in these motors. This structure provides a starting point from which to design highly specific myosin modulators to treat several human diseases. It further illustrates the potential of targeting transition intermediates of molecular machines to develop exquisitely selective pharmacological agents.

ATP-dependent interplay between local and global conformational changes in the myosin motor

The ATPase active site of myosin is located at the core of the motor head. During the Lymn-Taylor actomyosin contractile cycle, small conformational changes in the active site upon ATP binding, ATP hydrolysis and ADP/P_i release are accompanied by large conformational transitions of the motor domains, such as opening and closing of the actin binding cleft and the movement of lever arm. Here, our previous computational studies of myosin are summarized in a comprehensive model at the level of atomic detail. Molecular movies show how the successive domain motions during the ATP induced actin dissociation and the recovery stroke are coupled with the precise positioning of the key catalytic groups in the active site. This leads to a precise timing of the activation of the ATPase function: it allows ATP hydrolysis only after unbinding from actin and the priming of the lever arm, both pre-requisites for an efficient functioning of the motor during the subsequent power stroke. These coupling mechanisms constitute essential principles of every myosin motor, of which the ATP-site can be seen as the central allosteric control unit. © 2016 Wiley Periodicals, Inc.

Piezoelectric allostery of protein

Allostery is indispensable for a protein to work, where a locally applied stimulus is transmitted to a distant part of the molecule. While the allostery due to chemical stimuli such as ligand binding has long been studied, the growing interest in mechanobiology prompts the study of the mechanically stimulated allostery, the physical mechanism of which has not been established. By molecular dynamics simulation of a motor protein myosin, we found that a locally applied mechanical stimulus induces electrostatic potential change at distant regions, just like the piezoelectricity. This novel allosteric mechanism, "piezoelectric allostery", should be of particularly high value for mechanosensor/transducer proteins.

The Relay/Converter Interface Influences Hydrolysis of ATP by Skeletal Muscle Myosin II

The interface between relay and converter domain of muscle myosin is critical for optimal myosin performance. Using Drosophila melanogaster indirect flight muscle S1, we performed a kinetic analysis of the effect of mutations in the converter and relay domain. Introduction of a mutation (R759E) in the converter domain inhibits the steady-state ATPase of myosin S1, whereas an additional mutation in the relay domain (N509K) is able to restore the ATPase toward wild-type values. The R759E S1 construct showed little effect on most steps of the actomyosin ATPase cycle. The exception was a 25-30% reduction in the rate constant of the hydrolysis step, the step coupled to the cross-bridge recovery stroke that involves a change in conformation at the relay/converter domain interface. Significantly, the double mutant restored the hydrolysis step to values similar to the wild-type myosin. Modeling the relay/converter interface suggests a possible interaction between converter residue 759 and relay residue 509 in the actin-detached conformation, which is lost in R759E but is restored in N509K/R759E. This detailed kinetic analysis of Drosophila myosin carrying the R759E mutation shows that the interface between the relay loop and converter domain is important for fine-tuning myosin kinetics, in particular ATP binding and hydrolysis.

Structural basis for drug-induced allosteric changes to human β-cardiac myosin motor activity

Omecamtiv Mecarbil (OM) is a small molecule allosteric effector of cardiac myosin that is in clinical trials for treatment of systolic heart failure. A detailed kinetic analysis of cardiac myosin has shown that the drug accelerates phosphate release by shifting the equilibrium of the hydrolysis step towards products, leading to a faster transition from weak to strong actin-bound states. The structure of the human β-cardiac motor domain (cMD) with OM bound reveals a single OM-binding site nestled in a narrow cleft separating two domains of the human cMD where it interacts with the key residues that couple lever arm movement to the nucleotide state. In addition, OM induces allosteric changes in three strands of the β-sheet that provides the communication link between the actin-binding interface and the nucleotide pocket. The OM-binding interactions and allosteric changes form the structural basis for the kinetic and mechanical tuning of cardiac myosin.

Poorly understood aspects of striated muscle contraction

Muscle contraction results from cyclic interactions between the contractile proteins myosin and actin, driven by the turnover of adenosine triphosphate (ATP). Despite intense studies, several molecular events in the contraction process are poorly understood, including the relationship between force-generation and phosphate-release in the ATP-turnover. Different aspects of the force-generating transition are reflected in the changes in tension development by muscle cells, myofibrils and single molecules upon changes in temperature, altered phosphate concentration, or length perturbations. It has been notoriously difficult to explain all these events within a given theoretical framework and to unequivocally correlate observed events with the atomic structures of the myosin motor. Other incompletely understood issues include the role of the two heads of myosin II and structural changes in the actin filaments as well as the importance of the three-dimensional order. We here review these issues in relation to controversies regarding basic physiological properties of striated muscle. We also briefly consider actomyosin mutation effects in cardiac and skeletal muscle function and the possibility to treat these defects by drugs.

Omecamtiv Mecarbil modulates the kinetic and motile properties of porcine β-cardiac myosin

We determined the effect of Omecamtiv Mecarbil, a novel allosteric effector of cardiac muscle myosin, on the kinetic and "in vitro" motility properties of the porcine ventricular heavy meromyosin (PV-HMM). Omecamtiv Mecarbil increases the equilibrium constant of the hydrolysis step (M-ATP ⇄ M-ADP-Pi) from 2.4 to 6 as determined by quench flow, but the maximal rates of both the hydrolysis step and tryptophan fluorescence increase are unchanged by the drug. OM also increases the amplitude of the fast phase of phosphate dissociation (AM-ADP-Pi → AM-ADP + Pi) that is associated with force production in muscle by 4-fold. These results suggest a mechanism in which hydrolysis of M-ATP to M-ADP-Pi occurs both before and after the recovery stroke, but rapid acceleration of phosphate dissociation by actin occurs only on post-recovery stroke A-M-ADP-Pi. One of the more dramatic effects of OM on PV-HMM is a 14-fold decrease in the unloaded shortening velocity measured by the in vitro motility assay. The increase in flux through phosphate dissociation and the unchanged rate of ADP dissociation (AM-ADP → AM + ADP) by the drug produce a higher duty ratio motor in which a larger fraction of myosin heads are strongly bound to actin filaments. The increased internal load produced by a larger fraction of strongly attached crossbridges explains the reduced rate of in vitro motility velocity in the presence of OM and predicts that the drug will produce slower and stronger contraction of cardiac muscle.

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.

Morphing methods to visualize coarse-grained protein dynamics

Morphing was initially developed as a cinematic effect, where one image is seamlessly transformed into another image. The technique was widely adopted by biologists to visualize the transition between protein conformational states, generating an interpolated pathway from an initial to a final protein structure. Geometric morphing seeks to create visually suggestive movies that illustrate structural changes between conformations but do not necessarily represent a biologically relevant pathway, while minimum energy path (MEP) interpolations aim at describing the true transition state between the crystal structure minima in the energy landscape.

Stabilization of the ADP/metaphosphate intermediate during ATP hydrolysis in pre-power stroke myosin: quantitative anatomy of an enzyme

It has been proposed recently that ATP hydrolysis in ATPase enzymes proceeds via an initial intermediate in which the dissociated γ-phosphate of ATP is bound in the protein as a metaphosphate (PγO3(-)). A combined quantum/classical analysis of this dissociated nucleotide state inside myosin provides a quantitative understanding of how the enzyme stabilizes this unusual metaphosphate. Indeed, in vacuum, the energy of the ADP(3-) · PγO3(-) · Mg(2+) complex is much higher than that of the undissociated ATP(4-). The protein brings it to a surprisingly low value. Energy decomposition reveals how much each interaction in the protein stabilizes the metaphosphate state; backbone peptides of the P-loop contribute 50% of the stabilization energy, and the side chain of Lys-185(+) contributes 25%. This can be explained by the fact that these groups make strong favorable interactions with the α- and β-phosphates, thus favoring the charge distribution of the metaphosphate state over that of the ATP state. Further stabilization (16%) is achieved by a hydrogen bond between the backbone C=O of Ser-237 (on loop Switch-1) and a water molecule perfectly positioned to attack the PγO3(-) in the subsequent hydrolysis step. The planar and singly negative PγO3(-) is a much better target for the subsequent nucleophilic attack by a negatively charged OH(-) than the tetrahedral and doubly negative PγO4(2-) group of ATP. Therefore, we argue that the present mechanism of metaphosphate stabilization is common to the large family of nucleotide-hydrolyzing enzymes. Methodologically, this work presents a computational approach that allows us to obtain a truly quantitative conception of enzymatic strategy.

The myosin start-of-power stroke state and how actin binding drives the power stroke

We propose that on binding to actin at the start of the power stroke the myosin cross-bridge takes on the rigor configuration at the actin interface. Starting from the prepower stroke state, this can be achieved by a small movement (16° rotation) of the lower 50K domain without twisting the central β-sheet or opening switch-1 or switch-2. The movement of the lower 50K domain puts a strain on the W-helix. This strain tries to twist the β-sheet, which could drive the power stroke. This would provide a coupling between actin binding and the execution of the power stroke. During the power stroke the β-sheet twists, moving the P-loop away from switch-2, which opens the nucleotide binding pocket and separates ADP from Pi . The power stroke is different from the recovery stroke because the upper and lower 50K domains are tethered in the rigor configuration.

The superfast human extraocular myosin is kinetically distinct from the fast skeletal IIa, IIb, and IId isoforms

Humans express five distinct myosin isoforms in the sarcomeres of adult striated muscle (fast IIa, IId, the slow/cardiac isoform I/β, the cardiac specific isoform α, and the specialized extraocular muscle isoform). An additional isoform, IIb, is present in the genome but is not normally expressed in healthy human muscles. Muscle fibers expressing each isoform have distinct characteristics including shortening velocity. Defining the properties of the isoforms in detail has been limited by the availability of pure samples of the individual proteins. Here we study purified recombinant human myosin motor domains expressed in mouse C₂C₁₂ muscle cells. The results of kinetic analysis show that among the closely related adult skeletal isoforms, the affinity of ADP for actin·myosin (K_AD) is the characteristic that most readily distinguishes the isoforms. The three fast muscle myosins have K_AD values of 118, 80, and 55 μm for IId, IIa, and IIb, respectively, which follows the speed in motility assays from fastest to slowest. Extraocular muscle is unusually fast with a far weaker K_AD = 352 μm. Sequence comparisons and homology modeling of the structures identify a few key areas of sequence that may define the differences between the isoforms, including a region of the upper 50-kDa domain important in signaling between the nucleotide pocket and the actin-binding site.

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.

The mechanism of the converter domain rotation in the recovery stroke of myosin motor protein

Upon ATP binding, myosin motor protein is found in two alternative conformations, prerecovery state M* and postrecovery state M**. The transition from one state to the other, known as the recovery stroke, plays a key role in the myosin functional cycle. Despite much recent research, the microscopic details of this transition remain elusive. A critical step in the recovery stroke is the rotation of the converter domain from "up" position in prerecovery state to "down" position in postrecovery state that leads to the swing of the lever arm attached to it. In this work, we demonstrate that the two rotational states of the converter domain are determined by the interactions within a small structural motif in the force-generating region of the protein that can be accurately modeled on computers using atomic representation and explicit solvent. Our simulations show that the transition between the two states is controlled by a small helix (SH1) located next to the relay helix and relay loop. A small translation in the position of SH1 away from the relay helix is seen to trigger the transition from "up" state to "down" state. The transition is driven by a cluster of hydrophobic residues I687, F487, and F506 that make significant contributions to the stability of both states. The proposed mechanism agrees well with the available structural and mutational studies.

Interactions between relay helix and Src homology 1 (SH1) domain helix drive the converter domain rotation during the recovery stroke of myosin II

Myosin motor protein exists in two alternative conformations, prerecovery state M* and postrecovery state M**, on adenosine triphosphate binding. The details of the M*-to-M** transition, known as the recovery stroke to reflect its role as the functional opposite of the force-generating power stroke, remain elusive. The defining feature of the postrecovery state is a kink in the relay helix, a key part of the protein involved in force generation. In this article, we determine the interactions that are responsible for the appearance of the kink. We design a series of computational models that contain three other segments, relay loop, converter domain, and Src homology 1 (SH1) domain helix, with which relay helix interacts and determine their structure in accurate replica exchange molecular dynamics simulations in explicit solvent. By conducting an exhaustive combinatorial search among different models, we find that: (1) the converter domain must be attached to the relay helix during the transition, so it does not interfere with other parts of the protein and (2) the structure of the relay helix is controlled by SH1 helix. The kink is strongly coupled to the position of SH1 helix. It arises as a result of direct interactions between SH1 and the relay helix and leads to a rotation of the C-terminal part of the relay helix, which is subsequently transmitted to the converter domain.

How myosin motors power cellular functions: an exciting journey from structure to function: based on a lecture delivered at the 34th FEBS Congress in Prague, Czech Republic, July 2009

Molecular motors such as myosins are allosteric enzymes that power essential motility functions in the cell. Structural biology is an important tool for deciphering how these motors work. Myosins produce force upon the actin-driven conformational changes controlling the sequential release of the hydrolysis products of ATP (Pi followed by ADP). These conformational changes are amplified by a 'lever arm', which includes the region of the motor known as the converter and the adjacent elongated light chain binding region. Analysis of four structural states of the motor provides a detailed understanding of the rearrangements and pathways of communication in the motor that are necessary for detachment from the actin track and repriming of the motor. However, the important part of the cycle in which force is produced remains enigmatic and awaits new high-resolution structures. The value of a structural approach is particularly evident from clues provided by the structural states of the reverse myosin VI motor. Crystallographic structures have revealed that rearrangements within the converter subdomain occur, which explains why this myosin can produce a large stroke in the opposite direction to all other myosins, despite a very short lever arm. By providing a detailed understanding of the motor rearrangements, structural biology will continue to reveal essential information and help solve current enigma, such as how actin promotes force production, how motors are tuned for specific cellular roles or how motor/cargo interactions regulate the function of myosin in the cell.

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.

Disrupting the myosin converter-relay interface impairs Drosophila indirect flight muscle performance

Structural interactions between the myosin converter and relay domains have been proposed to be critical for the myosin power stroke and muscle power generation. We tested this hypothesis by mutating converter residue 759, which interacts with relay residues I508, N509, and D511, to glutamate (R759E) and determined the effect on Drosophila indirect flight muscle mechanical performance. Work loop analysis of mutant R759E indirect flight muscle fibers revealed a 58% and 31% reduction in maximum power generation (P(WL)) and the frequency at which maximum power (f(WL)) is generated, respectively, compared to control fibers at 15 °C. Small amplitude sinusoidal analysis revealed a 30%, 36%, and 32% reduction in mutant elastic modulus, viscous modulus, and mechanical rate constant 2πb, respectively. From these results, we infer that the mutation reduces rates of transitions through work-producing cross-bridge states and/or force generation during strongly bound states. The reductions in muscle power output, stiffness, and kinetics were physiologically relevant, as mutant wing beat frequency and flight index decreased about 10% and 45% compared to control flies at both 15 °C and 25 °C. Thus, interactions between the relay loop and converter domain are critical for lever-arm and catalytic domain coordination, high muscle power generation, and optimal Drosophila flight performance.

Modeling of relay helix functional dynamics and feasibility of experimental verification by neutron scattering

Cellular long-range transport involves motor proteins (MPs) (especially, kinesin and myosin) which contain a so-called relay helix. Its motion is of crucial importance to the conversion of chemical energy released in ATP hydrolysis into the coordinated mechanical movement of the entire motor protein. In this paper, we propose two combined nonlinear mechanisms for this particular functional activity and suggest the application of neutron scattering assays to experimentally determine the incoherent dynamic structure factor S(q,ω). We argue that this type of experiment is not only feasible but it could offer significant insights into the mechanism of MP function at a molecular level.

Early stages of the recovery stroke in myosin II studied by molecular dynamics simulations

The recovery stroke is a key step in the functional cycle of muscle motor protein myosin, during which pre-recovery conformation of the protein is changed into the active post-recovery conformation, ready to exersice force. We study the microscopic details of this transition using molecular dynamics simulations of atomistic models in implicit and explicit solvent. In more than 2 μs of aggregate simulation time, we uncover evidence that the recovery stroke is a two-step process consisting of two stages separated by a time delay. In our simulations, we directly observe the first stage at which switch II loop closes in the presence of adenosine triphosphate at the nucleotide binding site. The resulting configuration of the nucleotide binding site is identical to that detected experimentally. Distribution of inter-residue distances measured in the force generating region of myosin is in good agreement with the experimental data. The second stage of the recovery stroke structural transition, rotation of the converter domain, was not observed in our simulations. Apparently it occurs on a longer time scale. We suggest that the two parts of the recovery stroke need to be studied using separate computational models.

Free energy simulations of a GTPase: GTP and GDP binding to archaeal initiation factor 2

Archaeal initiation factor 2 (aIF2) is a protein involved in the initiation of protein biosynthesis. In its GTP-bound, "ON" conformation, aIF2 binds an initiator tRNA and carries it to the ribosome. In its GDP-bound, "OFF" conformation, it dissociates from tRNA. To understand the specific binding of GTP and GDP and its dependence on the ON or OFF conformational state of aIF2, molecular dynamics free energy simulations (MDFE) are a tool of choice. However, the validity of the computed free energies depends on the simulation model, including the force field and the boundary conditions, and on the extent of conformational sampling in the simulations. aIF2 and other GTPases present specific difficulties; in particular, the nucleotide ligand coordinates a divalent Mg(2+) ion, which can polarize the electronic distribution of its environment. Thus, a force field with an explicit treatment of electronic polarizability could be necessary, rather than a simpler, fixed charge force field. Here, we begin by comparing a fixed charge force field to quantum chemical calculations and experiment for Mg(2+):phosphate binding in solution, with the force field giving large errors. Next, we consider GTP and GDP bound to aIF2 and we compare two fixed charge force fields to the recent, polarizable, AMOEBA force field, extended here in a simple, approximate manner to include GTP. We focus on a quantity that approximates the free energy to change GTP into GDP. Despite the errors seen for Mg(2+):phosphate binding in solution, we observe a substantial cancellation of errors when we compare the free energy change in the protein to that in solution, or when we compare the protein ON and OFF states. Finally, we have used the fixed charge force field to perform MDFE simulations and alchemically transform GTP into GDP in the protein and in solution. With a total of about 200 ns of molecular dynamics, we obtain good convergence and a reasonable statistical uncertainty, comparable to the force field uncertainty, and somewhat lower than the predicted GTP/GDP binding free energy differences. The sign and magnitudes of the differences can thus be interpreted at a semiquantitative level, and are found to be consistent with the experimental binding preferences of ON- and OFF-aIF2.

Coarse-grained modeling of conformational transitions underlying the processive stepping of myosin V dimer along filamentous actin

To explore the structural basis of processive stepping of myosin V along filamentous actin, we have performed comprehensive modeling of its key conformational states and transitions with an unprecedented residue level of details. We have built structural models for a myosin V monomer complexed with filamentous actin at four biochemical states [adenosine diphosphate (ATP)-, adenosine diphosphate (ADP)-phosphate-, ADP-bound or nucleotide-free]. Then we have modeled a myosin V dimer (consisting of lead and rear head) at various two-head-bound states with nearly straight lever arms rotated by intramolecular strain. Next, we have performed transition pathway modeling to determine the most favorable sequence of transitions (namely, phosphate release at the lead head followed by ADP release at the rear head, while ADP release at the lead head is inhibited), which underlie the kinetic coordination between the two heads. Finally, we have used transition pathway modeling to reveal the order of structural changes during three key biochemical transitions (phosphate release at the lead head, ADP release and ATP binding at the rear head), which shed lights on the strain-dependence of the allosterically coupled motions at various stages of myosin V's work cycle. Our modeling results are in agreement with and offer structural insights to many results of kinetic, single-molecule and structural studies of myosin V.

Structural mechanism of the ATP-induced dissociation of rigor myosin from actin

Myosin is a true nanomachine, which produces mechanical force from ATP hydrolysis by cyclically interacting with actin filaments in a four-step cycle. The principle underlying each step is that structural changes in separate regions of the protein must be mechanically coupled. The step in which myosin dissociates from tightly bound actin (the rigor state) is triggered by the 30 Å distant binding of ATP. Large conformational differences between the crystal structures make it difficult to perceive the coupling mechanism. Energetically accessible transition pathways computed at atomic detail reveal a simple coupling mechanism for the reciprocal binding of ATP and actin.

Unsuspected pathway of the allosteric transition in hemoglobin

Large conformational transitions play an essential role in the function of many proteins, but experiments do not provide the atomic details of the path followed in going from one end structure to the other. For the hemoglobin tetramer, the transition path between the unliganded (T) and tetraoxygenated (R) structures is not known, which limits our understanding of the cooperative mechanism in this classic allosteric system, where both tertiary and quaternary changes are involved. The conjugate peak refinement algorithm is used to compute an unbiased minimum energy path at atomic detail between the two end states. Although the results confirm some of the proposals of Perutz [Perutz MF (1970) Stereochemistry of cooperative effects in haemoglobin. Nature 228:726-734], the subunit motions do not follow the textbook description of a simple rotation of one αβ-dimer relative to the other. Instead, the path consists of two sequential quaternary rotations, each involving different subdomains and axes. The quaternary transitions are preceded and followed by phases of tertiary structural changes. The results explain the recent photodissociation measurements, which suggest that the quaternary transition has a fast (2 μs) as well as a slow (20 μs) component and provide a testable model for single molecule FRET experiments.

Predicting order of conformational changes during protein conformational transitions using an interpolated elastic network model

The decryption of sequence of structural events during protein conformational transitions is essential to a detailed understanding of molecular functions of various biological nanomachines. Coarse-grained models have proven useful by allowing highly efficient simulations of protein conformational dynamics. By combining two coarse-grained elastic network models constructed based on the beginning and end conformations of a transition, we have developed an interpolated elastic network model to generate a transition pathway between the two protein conformations. For validation, we have predicted the order of local and global conformational changes during key ATP-driven transitions in three important biological nanomachines (myosin, F(1) ATPase and chaperonin GroEL). We have found that the local conformational change associated with the closing of active site precedes the global conformational change leading to mechanical motions. Our finding is in good agreement with the distribution of intermediate experimental structures, and it supports the importance of local motions at active site to drive or gate various conformational transitions underlying the workings of a diverse range of biological nanomachines.

Peripheral fatigue: high-energy phosphates and hydrogen ions

Peripheral fatigue results from an overactivity-induced decline in muscle function that originates from non-central nervous system mechanisms. A common symptom of fatigue is a feeling of tiredness or weariness because of overexertion, such as that associated with intense or prolonged physical exercise. Fatigue is worsened by low physical fitness and chronic illnesses. These conditions may intensify fatigue to levels that limit physical and social functioning and severely diminish health-related quality of life. Although etiologic aspects of peripheral fatigue are often associated with regulatory system (neurologic, endocrine, immunologic, muscular) and support system (cardiovascular, pulmonary, metabolic, renal, digestive, skeletal) limitations, final mediation occurs in muscle cells as a result of altered crossbridge functioning. Specifically, the final product and ionic metabolite accumulation that result from adenosine triphosphate hydrolysis appear to inhibit crossbridge formation and activation. Thus, clinical manifestations of peripheral fatigue often can be observed as limitations placed upon muscle or cardiorespiratory endurance, here defined as fatigue resistance. An overview of the common pathways by which peripheral fatigue can be mediated is provided. Product inhibition of contractile chemistry is brought into focus as a common pathway through which the mechanisms of peripheral fatigue often act.

Alpha-cardiac myosin heavy chain (MYH6) mutations affecting myofibril formation are associated with congenital heart defects

Congenital heart defects (CHD) are collectively the most common form of congenital malformation. Studies of human cases and animal models have revealed that mutations in several genes are responsible for both familial and sporadic forms of CHD. We have previously shown that a mutation in MYH6 can cause an autosomal dominant form of atrial septal defect (ASD), whereas others have identified mutations of the same gene in patients with hypertrophic and dilated cardiomyopathy. In the present study, we report a mutation analysis of MYH6 in patients with a wide spectrum of sporadic CHD. The mutation analysis of MYH6 was performed in DNA samples from 470 cases of isolated CHD using denaturing high-performance liquid chromatography and sequence analysis to detect point mutations and small deletions or insertions, and multiplex amplifiable probe hybridization to detect partial or complete copy number variations. One non-sense mutation, one splicing site mutation and seven non-synonymous coding mutations were identified. Transfection of plasmids encoding mutant and non-mutant green fluorescent protein-MYH6 fusion proteins in mouse myoblasts revealed that the mutations A230P and A1366D significantly disrupt myofibril formation, whereas the H252Q mutation significantly enhances myofibril assembly in comparison with the non-mutant protein. Our data indicate that functional variants of MYH6 are associated with cardiac malformations in addition to ASD and provide a novel potential mechanism. Such phenotypic heterogeneity has been observed in other genes mutated in CHD.

Targeting of the hair cell proteins cadherin 23, harmonin, myosin XVa, espin, and prestin in an epithelial cell model

We have developed an advantageous epithelial cell transfection model for examining the targeting, interactions, and mutations of hair cell proteins. When expressed in LLC-PK1-CL4 epithelial cells (CL4 cells), the outer hair cell protein prestin showed faithful domain-specific targeting to the basolateral plasma membrane. We examined the consequences of mutations affecting prestin activity and assigned a targeting role to the cytoplasmic tail. The stereociliary link protein cadherin 23 (Cdh23) was targeted to the plasma membrane of CL4 cell microvilli, the topological equivalent of stereocilia. In cells coexpressing the Cdh23 cytoplasmic binding protein harmonin, a large fraction of harmonin became colocalized with Cdh23 in microvilli. Using this assay and in vitro protein binding assays, we formulated an alternative model for Cdh23-harmonin binding, in which the primary interaction is between the harmonin N-domain and a 35-residue internal peptide in the Cdh23 cytoplasmic tail. Contrary to a previous model, we found no role for the Cdh23 C-terminal PDZ (PSD-95/Dlg/ZO-1)-binding motif and observed that Cdh23 bound similar levels of harmonin with or without the exon 68 peptide. We also examined two proteins involved in stereocilium elongation. The stereociliary actin-bundling protein espin was targeted to CL4 cell microvilli and caused microvillar elongation, whereas espin with the c.2469delGTCA or c.1988delAGAG human deafness mutation showed defects in microvillar targeting and elongation. The unconventional myosin motor myosin XVa accumulated at the tips of espin-elongated microvilli, by analogy to its location in stereocilia, whereas myosin XVa with the c.4351G>A or c.4669A>G human deafness mutation did not, revealing functional deficits in motor activity.

Mutating the converter-relay interface of Drosophila myosin perturbs ATPase activity, actin motility, myofibril stability and flight ability

We used an integrative approach to probe the significance of the interaction between the relay loop and converter domain of the myosin molecular motor from Drosophila melanogaster indirect flight muscle. During the myosin mechanochemical cycle, ATP-induced twisting of the relay loop is hypothesized to reposition the converter, resulting in cocking of the contiguous lever arm into the pre-power stroke configuration. The subsequent movement of the lever arm through its power stroke generates muscle contraction by causing myosin heads to pull on actin filaments. We generated a transgenic line expressing myosin with a mutation in the converter domain (R759E) at a site of relay loop interaction. Molecular modeling suggests that the interface between the relay loop and converter domain of R759E myosin would be significantly disrupted during the mechanochemical cycle. The mutation depressed calcium as well as basal and actin-activated MgATPase (V(max)) by approximately 60% compared to wild-type myosin, but there is no change in apparent actin affinity (K(m)). While ATP or AMP-PNP (adenylyl-imidodiphosphate) binding to wild-type myosin subfragment-1 enhanced tryptophan fluorescence by approximately 15% or approximately 8%, respectively, enhancement does not occur in the mutant. This suggests that the mutation reduces lever arm movement. The mutation decreases in vitro motility of actin filaments by approximately 35%. Mutant pupal indirect flight muscles display normal myofibril assembly, myofibril shape, and double-hexagonal arrangement of thick and thin filaments. Two-day-old fibers have occasional "cracking" of the crystal-like array of myofilaments. Fibers from 1-week-old adults show more severe cracking and frayed myofibrils with some disruption of the myofilament lattice. Flight ability is reduced in 2-day-old flies compared to wild-type controls, with no upward mobility but some horizontal flight. In 1-week-old adults, flight capability is lost. Thus, altered myosin function permits myofibril assembly, but results in a progressive disruption of the myofilament lattice and flight ability. We conclude that R759 in the myosin converter domain is essential for normal ATPase activity, in vitro motility and locomotion. Our results provide the first mutational evidence that intramolecular signaling between the relay loop and converter domain is critical for myosin function both in vitro and in muscle.

Multiscale modeling of structural dynamics underlying force generation and product release in actomyosin complex

To decrypt the mechanistic basis of myosin motor function, it is essential to probe the conformational changes in actomyosin with high spatial and temporal resolutions. In a computational effort to meet this challenge, we have performed a multiscale modeling of the allosteric couplings and transition pathway of actomyosin complex by combining coarse-grained modeling of the entire complex with all-atom molecular dynamics simulations of the active site. Our modeling of allosteric couplings at the pre-powerstroke state has pinpointed key actin-activated couplings to distant myosin parts which are critical to force generation and the sequential release of phosphate and ADP. At the post-powerstroke state, we have identified isoform-dependent couplings which underlie the reciprocal coupling between actin binding and nucleotide binding in fast Myosin II, and load-dependent ADP release in Myosin V. Our modeling of transition pathway during powerstroke has outlined a clear sequence of structural events triggered by actin binding, which lead to subsequent force generation, twisting of central beta-sheet, and the sequential release of phosphate and ADP. Finally we have performed atomistic simulations of active-site dynamics based on an on-path "transition-state" myosin conformation, which has revealed significantly weakened coordination of phosphate by Switch II, and a disrupted key salt bridge between Switch I and II. Meanwhile, the coordination of MgADP by Switch I and P loop is less perturbed. As a result, the phosphate can be released prior to MgADP. This study has shed new lights on the controversy over the structural mechanism of actin-activated phosphate release and force generation in myosin motor.