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

Modeling and Simulation of the NMDA Receptor at Coarse-Grained and Atomistic Levels

N-Methyl-D-aspartate (NMDA) receptors are glutamate-gated excitatory channels that play essential roles in brain functions. While high-resolution structures were solved for an allosterically inhibited form of functional NMDA receptor, other key functional states (particularly the active open-channel state) have not yet been resolved at atomic resolutions. To decrypt the molecular mechanism of the NMDA receptor activation, structural modeling and simulation are instrumental in providing detailed information about the dynamics and energetics of the receptor in various functional states. In this chapter, we describe coarse-grained modeling of the NMDA receptor using an elastic network model and related modeling/analysis tools (e.g., normal mode analysis, flexibility and hotspot analysis, cryo-EM flexible fitting, and transition pathway modeling) based on available structures. Additionally, we show how to build an atomistic model of the active-state receptor with targeted molecular dynamics (MD) simulation and explore its energetics and dynamics with conventional MD simulation. Taken together, these modeling and simulation can offer rich structural and dynamic information which will guide experimental studies of the activation of this key receptor.

Investigating the structural dynamics of the PIEZO1 channel activation and inactivation by coarse-grained modeling

The PIEZO channels, a family of mechanosensitive channels in vertebrates, feature a fast activation by mechanical stimuli (eg, membrane tension) followed by a slower inactivation. Although a medium-resolution structure of the trimeric form of PIEZO1 was solved by cryo-electron microscopy (cryo-EM), key structural changes responsible for the channel activation and inactivation are still unknown. Toward decrypting the structural mechanism of the PIEZO1 activation and inactivation, we performed systematic coarse-grained modeling using an elastic network model and related modeling/analysis tools (ie, normal mode analysis, flexibility and hotspot analysis, correlation analysis, and cryo-EM-based hybrid modeling and flexible fitting). We identified four key motional modes that may drive the tension-induced activation and inactivation, with fast and slow relaxation time, respectively. These modes allosterically couple the lateral and vertical motions of the peripheral domains to the opening and closing of the intra-cellular vestibule, enabling external mechanical forces to trigger, and regulate the activation/inactivation transitions. We also calculated domain-specific flexibility profiles, and predicted hotspot residues at key domain-domain interfaces and hinges. Our results offer unprecedented structural and dynamic information, which is consistent with the literature on mutational and functional studies of the PIEZO channels, and will guide future studies of this important family of mechanosensitive channels.

Probing the Structural Dynamics of the NMDA Receptor Activation by Coarse-Grained Modeling

N-Methyl-D-aspartate (NMDA) receptors are glutamate-gated excitatory channels that play essential roles in brain functions. High-resolution structures have been solved for an allosterically inhibited and agonist-bound form of a functional NMDA receptor; however, other key functional states (particularly the active open-channel state) were only resolved at moderate resolutions by cryo-electron microscopy (cryo-EM). To decrypt the mechanism of the NMDA receptor activation, structural modeling is essential to provide presently missing information about structural dynamics. We performed systematic coarse-grained modeling using an elastic network model and related modeling/analysis tools (e.g., normal mode analysis, flexibility and hotspot analysis, cryo-EM flexible fitting, and transition pathway modeling) based on an active-state cryo-EM map. We observed extensive conformational changes that allosterically couple the extracellular regulatory and agonist-binding domains to the pore-forming trans-membrane domain (TMD), and validated these, to our knowledge, new observations against known mutational and functional studies. Our results predict two key modes of collective motions featuring shearing/twisting of the extracellular domains relative to the TMD, reveal subunit-specific flexibility profiles, and identify functional hotspot residues at key domain-domain interfaces. Finally, by examining the conformational transition pathway between the allosterically inhibited form and the active form, we predict a discrete sequence of domain motions, which propagate from the extracellular domains to the TMD. In summary, our results offer rich structural and dynamic information, which is consistent with the literature on structure-function relationships in NMDA receptors, and will guide in-depth studies on the activation dynamics of this important neurotransmitter receptor.

Probing the structural dynamics of the CRISPR-Cas9 RNA-guided DNA-cleavage system by coarse-grained modeling

In the adaptive immune systems of many bacteria and archaea, the Cas9 endonuclease forms a complex with specific guide/scaffold RNA to identify and cleave complementary target sequences in foreign DNA. This DNA targeting machinery has been exploited in numerous applications of genome editing and transcription control. However, the molecular mechanism of the Cas9 system is still obscure. Recently, high-resolution structures have been solved for Cas9 in different structural forms (e.g., unbound forms, RNA-bound binary complexes, and RNA-DNA-bound tertiary complexes, corresponding to an inactive state, a pre-target-bound state, and a cleavage-competent or product state), which offered key structural insights to the Cas9 mechanism. To further probe the structural dynamics of Cas9 interacting with RNA and DNA at the amino-acid level of details, we have performed systematic coarse-grained modeling using an elastic network model and related analyses. Our normal mode analysis predicted a few key modes of collective motions that capture the observed conformational changes featuring large domain motions triggered by binding of RNA and DNA. Our flexibility analysis identified specific regions with high or low flexibility that coincide with key functional sites (such as DNA/RNA-binding sites, nuclease cleavage sites, and key hinges). We also identified a small set of hotspot residues that control the energetics of functional motions, which overlap with known functional sites and offer promising targets for future mutagenesis efforts to improve the specificity of Cas9. Finally, we modeled the conformational transitions of Cas9 from the unbound form to the binary complex and then the tertiary complex, and predicted a distinct sequence of domain motions. In sum, our findings have offered rich structural and dynamic details relevant to the Cas9 machinery, and will guide future investigation and engineering of the Cas9 systems. Proteins 2017; 85:342-353. © 2016 Wiley Periodicals, Inc.

Rigid Residue Scan Simulations Systematically Reveal Residue Entropic Roles in Protein Allostery

Intra-protein information is transmitted over distances via allosteric processes. This ubiquitous protein process allows for protein function changes due to ligand binding events. Understanding protein allostery is essential to understanding protein functions. In this study, allostery in the second PDZ domain (PDZ2) in the human PTP1E protein is examined as model system to advance a recently developed rigid residue scan method combining with configurational entropy calculation and principal component analysis. The contributions from individual residues to whole-protein dynamics and allostery were systematically assessed via rigid body simulations of both unbound and ligand-bound states of the protein. The entropic contributions of individual residues to whole-protein dynamics were evaluated based on covariance-based correlation analysis of all simulations. The changes of overall protein entropy when individual residues being held rigid support that the rigidity/flexibility equilibrium in protein structure is governed by the La Châtelier's principle of chemical equilibrium. Key residues of PDZ2 allostery were identified with good agreement with NMR studies of the same protein bound to the same peptide. On the other hand, the change of entropic contribution from each residue upon perturbation revealed intrinsic differences among all the residues. The quasi-harmonic and principal component analyses of simulations without rigid residue perturbation showed a coherent allosteric mode from unbound and bound states, respectively. The projection of simulations with rigid residue perturbation onto coherent allosteric modes demonstrated the intrinsic shifting of ensemble distributions supporting the population-shift theory of protein allostery. Overall, the study presented here provides a robust and systematic approach to estimate the contribution of individual residue internal motion to overall protein dynamics and allostery.

A combined coarse-grained and all-atom simulation of TRPV1 channel gating and heat activation

Coarse-grained modeling and all-atom molecular dynamics simulation provide insight into the mechanism for heat activation of TRPV1 gating.

The transient receptor potential (TRP) channels act as key sensors of various chemical and physical stimuli in eukaryotic cells. Despite years of study, the molecular mechanisms of TRP channel activation remain unclear. To elucidate the structural, dynamic, and energetic basis of gating in TRPV1 (a founding member of the TRPV subfamily), we performed coarse-grained modeling and all-atom molecular dynamics (MD) simulation based on the recently solved high resolution structures of the open and closed form of TRPV1. Our coarse-grained normal mode analysis captures two key modes of collective motions involved in the TRPV1 gating transition, featuring a quaternary twist motion of the transmembrane domains (TMDs) relative to the intracellular domains (ICDs). Our transition pathway modeling predicts a sequence of structural movements that propagate from the ICDs to the TMDs via key interface domains (including the membrane proximal domain and the C-terminal domain), leading to sequential opening of the selectivity filter followed by the lower gate in the channel pore (confirmed by modeling conformational changes induced by the activation of ICDs). The above findings of coarse-grained modeling are robust to perturbation by lipids. Finally, our MD simulation of the ICD identifies key residues that contribute differently to the nonpolar energy of the open and closed state, and these residues are predicted to control the temperature sensitivity of TRPV1 gating. These computational predictions offer new insights to the mechanism for heat activation of TRPV1 gating, and will guide our future electrophysiology and mutagenesis studies.

Effects of ATP and actin-filament binding on the dynamics of the myosin II S1 domain

Actin and myosin interact with one another to perform a variety of cellular functions. Central to understanding the processive motion of myosin on actin is the characterization of the individual states along the mechanochemical cycle. We present an all-atom molecular dynamics simulation of the myosin II S1 domain in the rigor state interacting with an actin filament. We also study actin-free myosin in both rigor and post-rigor conformations. Using all-atom level and coarse-grained analysis methods, we investigate the effects of myosin binding on actin, and of actin binding on myosin. In particular, we determine the domains of actin and myosin that interact strongly with one another at the actomyosin interface using a highly coarse-grained level of resolution, and we identify a number of salt bridges and hydrogen bonds at the interface of myosin and actin. Applying coarse-grained analysis, we identify differences in myosin states dependent on actin-binding, or ATP binding. Our simulations also indicate that the actin propeller twist-angle and nucleotide cleft-angles are influenced by myosin at the actomyosin interface. The torsional rigidity of the myosin-bound filament is also calculated, and is found to be increased compared to previous simulations of the free filament.

Metal switch-controlled myosin II from Dictyostelium discoideum supports closure of nucleotide pocket during ATP binding coupled to detachment from actin filaments

G-proteins, kinesins, and myosins are hydrolases that utilize a common protein fold and divalent metal cofactor (typically Mg(2+)) to coordinate purine nucleotide hydrolysis. The nucleoside triphosphorylase activities of these enzymes are activated through allosteric communication between the nucleotide-binding site and the activator/effector/polymer interface to convert the free energy of nucleotide hydrolysis into molecular switching (G-proteins) or force generation (kinesins and myosin). We have investigated the ATPase mechanisms of wild-type and the S237C mutant of non-muscle myosin II motor from Dictyostelium discoideum. The S237C substitution occurs in the conserved metal-interacting switch-1, and we show that this substitution modulates the actomyosin interaction based on the divalent metal present in solution. Surprisingly, S237C shows rapid basal steady-state Mg(2+)- or Mn(2+)-ATPase kinetics, but upon binding actin, its MgATPase is inhibited. This actin inhibition is relieved by Mn(2+), providing a direct and experimentally reversible linkage of switch-1 and the actin-binding cleft through the swapping of divalent metals in the reaction. Using pyrenyl-labeled F-actin, we demonstrate that acto·S237C undergoes slow and weak MgATP binding, which limits the rate of steady-state catalysis. Mn(2+) rescues this effect to near wild-type activity. 2'(3')-O-(N-Methylanthraniloyl)-ADP release experiments show the need for switch-1 interaction with the metal cofactor for tight ADP binding. Our results are consistent with strong reciprocal coupling of nucleoside triphosphate and F-actin binding and provide additional evidence for the allosteric communication pathway between the nucleotide-binding site and the filament-binding region.

Structure-based simulations of the translocation mechanism of the hepatitis C virus NS3 helicase along single-stranded nucleic acid

The NS3 helicase of Hepatitis C virus is an ATP-fueled molecular motor that can translocate along single-stranded (ss) nucleic acid, and unwind double-stranded nucleic acids. It makes a promising antiviral target and an important prototype system for helicase research. Despite recent progress, the detailed mechanism of NS3 helicase remains unknown. In this study, we have combined coarse-grained (CG) and atomistic simulations to probe the translocation mechanism of NS3 helicase along ssDNA. At the residue level of detail, our CG simulations have captured functionally important interdomain motions of NS3 helicase and reproduced single-base translocation of NS3 helicase along ssDNA in the 3'-5' direction, which is in good agreement with experimental data and the inchworm model. By combining the CG simulations with residue-specific perturbations to protein-DNA interactions, we have identified a number of key residues important to the translocation machinery that agree with previous structural and mutational studies. Additionally, our atomistic simulations with targeted molecular dynamics have corroborated the findings of CG simulations and further revealed key protein-DNA hydrogen bonds that break/form during the transitions. This study offers, to our knowledge, the most detailed and realistic simulations of translocation mechanism of NS3 helicase. The simulation protocol established in this study will be useful for designing inhibitors that target the translocation machinery of NS3 helicase, and for simulations of a variety of nucleic-acid-based molecular motors.

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.

Coarse-grained modeling of the structural states and transition underlying the powerstroke of dynein motor domain

This study aims to model a minimal dynein motor domain capable of motor function, which consists of the linker domain, six AAA+ modules (AAA1-AAA6), coiled coil stalk, and C-terminus domain. To this end, we have used the newly solved X-ray structures of dynein motor domain to perform a coarse-grained modeling of dynein's post- and pre-powerstroke conformation and the conformational transition between them. First, we have used normal mode analysis to identify a single normal mode that captures the coupled motions of AAA1-AAA2 closing and linker domain rotation, which enables the ATP-driven recovery stroke of dynein. Second, based on the post-powerstroke conformation solved crystallographically, we have modeled dynein's pre-powerstroke conformation by computationally inducing AAA1-AAA2 closing and sliding of coiled coil stalk, and the resulting model features a linker domain near the pre-powerstroke position and a slightly tilted stalk. Third, we have modeled the conformational transition from pre- to post-powerstroke conformation, which predicts a clear sequence of structural events that couple microtubule binding, powerstroke and product release, and supports a signaling path from stalk to AAA1 via AAA3 and AAA4. Finally, we have found that a closed AAA3-AAA4 interface (compatible with nucleotide binding) is essential to the mechano-chemical coupling in dynein. Our modeling not only offers unprecedented structural insights to the motor function of dynein as described by past single-molecule, fluorescence resonance energy transfer, and electron microscopy studies, but also provides new predictions for future experiments to test.

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.

Myosin-V as a mechanical sensor: an elastic network study

According to recent experiments, the molecular-motor myosin behaves like a strain sensor, exhibiting different functional responses when loads in opposite directions are applied to its tail. Within an elastic-network model, we explore the sensitivity of the protein to the forces acting on the tail and find, in agreement with experiments, that such forces invoke conformational changes that should affect filament binding and ADP release. Furthermore, conformational responses of myosin to the application of forces to individual residues in its principal functional regions are systematically investigated and a detailed sensitivity map of myosin-V is thus obtained. The results suggest that the strain-sensor behavior is involved in the intrinsic operation of this molecular motor.

Extension of a three-helix bundle domain of myosin VI and key role of calmodulins

The molecular motor protein myosin VI moves toward the minus-end of actin filaments with a step size of 30-36 nm. Such large step size either drastically limits the degree of complex formation between dimer subunits to leave enough length for the lever arms, or requires an extension of the lever arms' crystallographically observed structure. Recent experimental work proposed that myosin VI dimerization triggers the unfolding of the protein's proximal tail domain which could drive the needed lever-arm extension. Here, we demonstrate through steered molecular dynamics simulation the feasibility of sufficient extension arising from turning a three-helix bundle into a long α-helix. A key role is played by the known calmodulin binding that facilitates the extension by altering the strain path in myosin VI. Sequence analysis of the proximal tail domain suggests that further calmodulin binding sites open up when the domain's three-helix bundle is unfolded and that subsequent calmodulin binding stabilizes the extended lever arms.

Force and power generating mechanism(s) in active muscle as revealed from temperature perturbation studies

The basic characteristics of the process of force and power generation in active muscle that have emerged from temperature studies are examined. This is done by reviewing complementary findings from temperature-dependence studies and rapid temperature-jump (T-jump) experiments and from intact and skinned fast mammalian muscle fibres. In isometric muscle, a small T-jump leads to a characteristic rise in force showing that crossbridge force generation is endothermic (heat absorbed) and associated with increased entropy (disorder). The sensitivity of the T-jump force generation to added inorganic phosphate (Pi) indicates that a T-jump enhances an early step in the actomyosin (crossbridge) ATPase cycle before Pi-release. During muscle lengthening when steady force is increased, the T-jump force generation is inhibited. Conversely, during shortening when steady force is decreased, the T-jump force generation is enhanced in a velocity-dependent manner, showing that T-jump force generation is strain sensitive. Within the temperature range of ∼5–35◦C, the temperature dependence of steady active force is sigmoidal both in isometric and in shortening muscle. However, in shortening muscle, the endothermic character of force generation becomes more pronounced with increased velocity and this can, at least partly, account for the marked increase with warming of the mechanical power output of active muscle.