The requirement for mechanical coupling between head and S2 domains in smooth muscle myosin ATPase regulation and its implications for dimeric motor function

The myosin II coiled-coil domain atomic structure in its native environment

The atomic structure of the complete myosin tail within thick filaments isolated from Lethocerus indicus flight muscle is described and compared to crystal structures of recombinant, human cardiac myosin tail segments. Overall, the agreement is good with three exceptions: the proximal S2, in which the filament has heads attached but the crystal structure doesn't, and skip regions 2 and 4. At the head-tail junction, the tail α-helices are asymmetrically structured encompassing well-defined unfolding of 12 residues for one myosin tail, ∼4 residues of the other, and different degrees of α-helix unwinding for both tail α-helices, thereby providing an atomic resolution description of coiled-coil "uncoiling" at the head-tail junction. Asymmetry is observed in the nonhelical C termini; one C-terminal segment is intercalated between ribbons of myosin tails, the other apparently terminating at Skip 4 of another myosin tail. Between skip residues, crystal and filament structures agree well. Skips 1 and 3 also agree well and show the expected α-helix unwinding and coiled-coil untwisting in response to skip residue insertion. Skips 2 and 4 are different. Skip 2 is accommodated in an unusual manner through an increase in α-helix radius and corresponding reduction in rise/residue. Skip 4 remains helical in one chain, with the other chain unfolded, apparently influenced by the acidic myosin C terminus. The atomic model may shed some light on thick filament mechanosensing and is a step in understanding the complex roles that thick filaments of all species undergo during muscle contraction.

Cryo-EM structure of the inhibited (10S) form of myosin II

Myosin II is the motor protein responsible for contractility in muscle and nonmuscle cells¹. The molecule has two identical heads attached to an elongated tail, and can exist in two conformations: 10S and 6S, named for their sedimentation coefficients^2,3. The 6S conformation has an extended tail and assembles into polymeric filaments, which pull on actin filaments to generate force and motion. In 10S myosin, the tail is folded into three segments and the heads bend back and interact with each other and the tail^3-7, creating a compact conformation, in which ATPase activity, actin activation and filament assembly are all highly inhibited^7,8. This switched-off structure appears to function as a key energy-conserving storage molecule in muscle and nonmuscle cells^9-12, which can be activated to form functional filaments as needed¹³ – but the mechanism of its inhibition is not understood. Here we have solved the structure of smooth muscle 10S myosin to a global resolution of 4.3 Å by cryo-EM, revealing near-atomic level details of its structure for the first time. The reconstruction provides a new understanding of the head and tail regions of the molecule and of the key intramolecular contacts that cause inhibition. Our results suggest an atomic model for the off-state of myosin II, for its activation and unfolding by phosphorylation, and for understanding the clustering of disease-causing mutations near sites of intramolecular interaction.

An interphase contractile ring reshapes primordial germ cells to allow bulk cytoplasmic remodeling

Some cells discard undesired inherited components in bulk by forming large compartments that are subsequently eliminated. Caenorhabditis elegans primordial germ cells (PGCs) jettison mitochondria and cytoplasm by forming a large lobe that is cannibalized by intestinal cells. Although PGCs are nonmitotic, we find that lobe formation is driven by constriction of a contractile ring and requires the RhoGEF ECT-2, a RhoA activator also essential for cytokinesis. Whereas centralspindlin activates ECT-2 to promote cytokinetic contractile ring formation, we show that the ECT-2 regulator NOP-1, but not centralspindlin, is essential for PGC lobe formation. We propose that lobe contractile ring formation is locally inhibited by the PGC nucleus, which migrates to one side of the cell before the cytokinetic ring assembles on the opposite cortex. Our findings reveal how components of the cytokinetic contractile ring are reemployed during interphase to create compartments used for cellular remodeling, and they reveal differences in the spatial cues that dictate where the contractile ring will form.

Cell-intrinsic and -extrinsic mechanisms promote cell-type-specific cytokinetic diversity

Cytokinesis, the physical division of one cell into two, is powered by constriction of an actomyosin contractile ring. It has long been assumed that all animal cells divide by a similar molecular mechanism, but growing evidence suggests that cytokinetic regulation in individual cell types has more variation than previously realized. In the four-cell Caenorhabditis elegans embryo, each blastomere has a distinct cell fate, specified by conserved pathways. Using fast-acting temperature-sensitive mutants and acute drug treatment, we identified cell-type-specific variation in the cytokinetic requirement for a robust formin^CYK-1-dependent filamentous-actin (F-actin) cytoskeleton. In one cell (P2), this cytokinetic variation is cell-intrinsically regulated, whereas in another cell (EMS) this variation is cell-extrinsically regulated, dependent on both Src^SRC-1 signaling and direct contact with its neighbor cell, P2. Thus, both cell-intrinsic and -extrinsic mechanisms control cytokinetic variation in individual cell types and can protect against division failure when the contractile ring is weakened.

The successful division of one cell into two is essential for all organisms to live, grow and reproduce. For an animal cell, the nucleus – the compartment containing the genetic material – must divide before the surrounding material. The rest of the cell, called the cytoplasm, physically separates later in a process known as cytokinesis.

Cytokinesis in animal cells is driven by the formation of a ring in the middle of the dividing cell. The ring is composed of myosin motor proteins and filaments made of a protein called actin. The movements of the motor proteins along the filaments cause the ring to contract and tighten. This pulls the cell membrane inward and physically pinches the cell into two. For a long time, the mechanism of cytokinesis was assumed to be same across different types of animal cell, but later evidence suggested otherwise. For example, in liver, heat and bone cells, cytokinesis naturally fails during development to create cells with two or more nuclei. If a similar ‘failure’ happened in other cell types, it could lead to diseases such as cancers or blood disorders. This raised the question: what are the molecular mechanisms that allow cytokinesis to happen differently in different cell types?

Davies et al. investigated this question using embryos of the worm Caenorhabditis elegans at a stage in their development when they consist of just four cells. The proteins forming the contractile ring in this worm are the same as those in humans. However, in the worm, the contractile ring can easily be damaged using chemical inhibitors or by mutating the genes that encode its proteins.

Davies et al. show that when the contractile ring was damaged, two of the four cells in the worm embryo still divided successfully. This result indicates the existence of new mechanisms to divide the cytoplasm that allow division even with a weak contractile ring. In a further experiment, the embryos were dissected to isolate each of the four cells. Davies et al. saw that one of the two dividing cells could still divide on its own, while the other cell could not. This shows that this new method of cytokinesis is regulated both by factors inherent to the dividing cell and by external signals from other cells. Moreover, one of these extrinsic signals was found to be a signaling protein that had previously been implicated in human cancers.

Future work will determine if these variations in cytokinesis between the different cell types found in the worm apply to humans too; and, more importantly from a therapeutic standpoint, if these new mechanisms exist in human cancers.

Coupling between myosin head conformation and the thick filament backbone structure

The recent high-resolution structure of the thick filament from Lethocerus asynchronous flight muscle shows aspects of thick filament structure never before revealed that may shed some light on how striated muscles function. The phenomenon of stretch activation underlies the function of asynchronous flight muscle. It is most highly developed in flight muscle, but is also observed in other striated muscles such as cardiac muscle. Although stretch activation is likely to be complex, involving more than a single structural aspect of striated muscle, the thick filament itself, would be a prime site for regulatory function because it must bear all of the tension produced by both its associated myosin motors and any externally applied force. Here we show the first structural evidence that the arrangement of myosin heads within the interacting heads motif is coupled to the structure of the thick filament backbone. We find that a change in helical angle of 0.16° disorders the blocked head preferentially within the Lethocerus interacting heads motif. This observation suggests a mechanism for how tension affects the dynamics of the myosin heads leading to a detailed hypothesis for stretch activation and shortening deactivation, in which the blocked head preferentially binds the thin filament followed by the free head when force production occurs.

Non-muscle myosin II is required for correct fate specification in the Caenorhabditis elegans seam cell divisions

During development, cell division often generates two daughters with different developmental fates. Distinct daughter identities can result from the physical polarity and size asymmetry itself, as well as the subsequent activation of distinct fate programmes in each daughter. Asymmetric divisions are a feature of the C. elegans seam lineage, in which a series of post-embryonic, stem-like asymmetric divisions give rise to an anterior daughter that differentiates and a posterior daughter that continues to divide. Here we have investigated the role of non-muscle myosin II (nmy-2) in these asymmetric divisions. We show that nmy-2 does not appear to be involved in generating physical division asymmetry, but nonetheless is important for specifying differential cell fate. While cell polarity appears normal, and chromosome and furrow positioning remains unchanged when nmy-2 is inactivated, seam cell loss occurs through inappropriate terminal differentiation of posterior daughters. This reveals a role for nmy-2 in cell fate determination not obviously linked to the primary polarity determination mechanisms it has been previously associated with.

Structure of myosin filaments from relaxed Lethocerus flight muscle by cryo-EM at 6 Å resolution

We describe a cryo-electron microscopy three-dimensional image reconstruction of relaxed myosin II-containing thick filaments from the flight muscle of the giant water bug Lethocerus indicus. The relaxed thick filament structure is a key element of muscle physiology because it facilitates the reextension process following contraction. Conversely, the myosin heads must disrupt their relaxed arrangement to drive contraction. Previous models predicted that Lethocerus myosin was unique in having an intermolecular head-head interaction, as opposed to the intramolecular head-head interaction observed in all other species. In contrast to the predicted model, we find an intramolecular head-head interaction, which is similar to that of other thick filaments but oriented in a distinctly different way. The arrangement of myosin's long α-helical coiled-coil rod domain has been hypothesized as either curved layers or helical subfilaments. Our reconstruction is the first report having sufficient resolution to track the rod α helices in their native environment at resolutions ~5.5 Å, and it shows that the layer arrangement is correct for Lethocerus. Threading separate paths through the forest of myosin coiled coils are four nonmyosin peptides. We suggest that the unusual position of the heads and the rod arrangement separated by nonmyosin peptides are adaptations for mechanical signal transduction whereby applied tension disrupts the myosin heads as a component of stretch activation.

Self-assembly Controls Self-cleavage of HHR from ASBVd (-): a Combined SANS and Modeling Study

In the Avocado Sunblotch Viroid (ASBVd: 249-nt) from the Avsunviroidae family, a symmetric rolling-circle replication operates through an autocatalytic mechanism mediated by hammerhead ribozymes (HHR) embedded in both polarity strands. The concatenated multimeric ASBVd (+) and ASBVd (-) RNAs thus generated are processed by cleavage to unit-length where ASBVd (-) self-cleaves with more efficiency. Absolute scale small angle neutron scattering (SANS) revealed a temperature-dependent dimer association in both ASBVd (-) and its derived 79-nt HHR (-). A joint thermodynamic analysis of SANS and catalytic data indicates the rate-determining step corresponds to the dimer/monomer transition. 2D and 3D models of monomeric and dimeric HHR (-) suggest that the inter-molecular contacts stabilizing the dimer (between HI and HII domains) compete with the intra-molecular ones stabilizing the active conformation of the full-length HHR required for an efficient self-cleavage. Similar competing intra- and inter-molecular contacts are proposed in ASBVd (-) though with a remoter region from an extension of the HI domain.

Principles and Overview of Sampling Methods for Modeling Macromolecular Structure and Dynamics

Investigation of macromolecular structure and dynamics is fundamental to understanding how macromolecules carry out their functions in the cell. Significant advances have been made toward this end in silico, with a growing number of computational methods proposed yearly to study and simulate various aspects of macromolecular structure and dynamics. This review aims to provide an overview of recent advances, focusing primarily on methods proposed for exploring the structure space of macromolecules in isolation and in assemblies for the purpose of characterizing equilibrium structure and dynamics. In addition to surveying recent applications that showcase current capabilities of computational methods, this review highlights state-of-the-art algorithmic techniques proposed to overcome challenges posed in silico by the disparate spatial and time scales accessed by dynamic macromolecules. This review is not meant to be exhaustive, as such an endeavor is impossible, but rather aims to balance breadth and depth of strategies for modeling macromolecular structure and dynamics for a broad audience of novices and experts.

High-resolution temporal analysis reveals a functional timeline for the molecular regulation of cytokinesis

To take full advantage of fast-acting temperature-sensitive mutations, thermal control must be extremely rapid. We developed the Therminator, a device capable of shifting sample temperature in ~17 s while simultaneously imaging cell division in vivo. Applying this technology to six key regulators of cytokinesis, we found that each has a distinct temporal requirement in the Caenorhabditis elegans zygote. Specifically, myosin-II is required throughout cytokinesis until contractile ring closure. In contrast, formin-mediated actin nucleation is only required during assembly and early contractile ring constriction. Centralspindlin is required to maintain division after ring closure, although its GAP activity is only required until just prior to closure. Finally, the chromosomal passenger complex is required for cytokinesis only early in mitosis, but not during metaphase or cytokinesis. Together, our results provide a precise functional timeline for molecular regulators of cytokinesis using the Therminator, a powerful tool for ultra-rapid protein inactivation.

Role of the essential light chain in the activation of smooth muscle myosin by regulatory light chain phosphorylation

The activity of smooth and non-muscle myosin II is regulated by phosphorylation of the regulatory light chain (RLC) at serine 19. The dephosphorylated state of full-length monomeric myosin is characterized by an asymmetric intramolecular head-head interaction that completely inhibits the ATPase activity, accompanied by a hairpin fold of the tail, which prevents filament assembly. Phosphorylation of serine 19 disrupts these head-head interactions by an unknown mechanism. Computational modeling (Tama et al., 2005. J. Mol. Biol. 345, 837-854) suggested that formation of the inhibited state is characterized by both torsional and bending motions about the myosin heavy chain (HC) at a location between the RLC and the essential light chain (ELC). Therefore, altering relative motions between the ELC and the RLC at this locus might disrupt the inhibited state. Based on this hypothesis we have derived an atomic model for the phosphorylated state of the smooth muscle myosin light chain domain (LCD). This model predicts a set of specific interactions between the N-terminal residues of the RLC with both the myosin HC and the ELC. Site directed mutagenesis was used to show that interactions between the phosphorylated N-terminus of the RLC and helix-A of the ELC are required for phosphorylation to activate smooth muscle myosin.

Structural basis of the relaxed state of a Ca2+-regulated myosin filament and its evolutionary implications

Myosin filaments of muscle are regulated either by phosphorylation of their regulatory light chains or Ca(2+) binding to the essential light chains, contributing to on-off switching or modulation of contraction. Phosphorylation-regulated filaments in the relaxed state are characterized by an asymmetric interaction between the two myosin heads, inhibiting their actin binding or ATPase activity. Here, we have tested whether a similar interaction switches off activity in myosin filaments regulated by Ca(2+) binding. Cryo-electron microscopy and single-particle image reconstruction of Ca(2+)-regulated (scallop) filaments reveals a helical array of myosin head-pair motifs above the filament surface. Docking of atomic models of scallop myosin head domains into the motifs reveals that the heads interact in a similar way to those in phosphorylation-regulated filaments. The results imply that the two major evolutionary branches of myosin regulation--involving phosphorylation or Ca(2+) binding--share a common structural mechanism for switching off thick-filament activity in relaxed muscle. We suggest that the Ca(2+)-binding mechanism evolved from the more ancient phosphorylation-based system to enable rapid response of myosin-regulated muscles to activation. Although the motifs are similar in both systems, the scallop structure is more tilted and higher above the filament backbone, leading to different intermolecular interactions. The reconstruction reveals how the myosin tail emerges from the motif, connecting the heads to the filament backbone, and shows that the backbone is built from supramolecular assemblies of myosin tails. The reconstruction provides a native structural context for understanding past biochemical and biophysical studies of this model Ca(2+)-regulated myosin.

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.

Phosphorylated smooth muscle heavy meromyosin shows an open conformation linked to activation

Smooth muscle myosin and smooth muscle heavy meromyosin (smHMM) are activated by regulatory light chain phosphorylation, but the mechanism remains unclear. Dephosphorylated, inactive smHMM assumes a closed conformation with asymmetric intramolecular head-head interactions between motor domains. The "free head" can bind to actin, but the actin binding interface of the "blocked head" is involved in interactions with the free head. We report here a three-dimensional structure for phosphorylated, active smHMM obtained using electron crystallography of two-dimensional arrays. Head-head interactions of phosphorylated smHMM resemble those found in the dephosphorylated state but occur between different molecules, not within the same molecule. The light chain binding domain structure of phosphorylated smHMM differs markedly from that of the "blocked" head of dephosphorylated smHMM. We hypothesize that regulatory light chain phosphorylation opens the inhibited conformation primarily by its effect on the blocked head. Singly phosphorylated smHMM is not compatible with the closed conformation if the blocked head is phosphorylated. This concept has implications for the extent of myosin activation at low levels of phosphorylation in smooth muscle.

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.

Role of the tail in the regulated state of myosin 2

Myosin 2 from vertebrate smooth muscle or non-muscle sources is in equilibrium between compact, inactive monomers and thick filaments under physiological conditions. In the inactive monomer, the two heads pack compactly together, and the long tail is folded into three closely packed segments that are associated chiefly with one of the heads. The molecular basis of the folding of the tail remains unexplained. By using electron microscopy, we show that compact monomers of smooth muscle myosin 2 have the same structure in both the native state and following specific, intramolecular photo-cross-linking between Cys109 of the regulatory light chain (RLC) and segment 3 of the tail. Nonspecific cross-linking between lysine residues of the folded monomer by glutaraldehyde also does not perturb the compact conformation and stabilizes it against unfolding at high ionic strength. Sequence comparisons across phyla and myosin 2 isoforms suggest that the folding of the tail is stabilized by ionic interactions between the positively charged N-terminal sequence of the RLC and a negatively charged region near the start of tail segment 3 and that phosphorylation of the RLC could perturb these interactions. Our results support the view that interactions between the heads and the distal tail perform a critical role in regulating activity of myosin 2 molecules through stabilizing the compact monomer conformation.

Modeling smooth muscle myosin's two heads: long-lived enzymatic roles and phosphorylation-dependent equilibria

Smooth muscle myosin has two heads, each capable of interacting with actin to generate force and/or motion as it hydrolyzes ATP. These heads are inhibited when their associated regulatory light chain is unphosphorylated (0P), becoming active and hydrolyzing ATP maximally when phosphorylated (2P). Interestingly, with only one of the two regulatory light chains phosphorylated (1P), smooth muscle myosin is active but its ATPase rate is <2P. To explain published 1P single ATP turnover and steady-state ATPase activities, we propose a kinetic model in which 1P myosin exists in an equilibrium between being fully active (2P) and inhibited (0P). Based on the single ATP turnover data, we also propose that each 2P head adopts a hydrolytic role distinct from its partner at any point in time, i.e., one head strongly binds actin and hydrolyzes ATP at its actin-activated rate while the other weakly binds actin. Surprisingly, the heads switch roles slowly (<0.1 s(-1)), suggesting that their activities are not independent. The phosphorylation-dependent equilibrium between active and inhibited states and the hydrolytic role that each head adopts during its interaction with actin may have implications for understanding regulation and mechanical performance of other members of the myosin family of molecular motors.

Structure modeling from small angle X-ray scattering data with elastic network normal mode analysis

Computational algorithms to construct structural models from SAXS experimental data are reviewed. SAXS data provides a wealth of information to study the structure and dynamics of biological molecules, however it does not provide atomic details of structures. Thus combining the low-resolution data with already known X-ray structure is a common approach to study conformational transitions of biological molecules. This review provides a survey of SAXS modeling approaches. In addition, we will discuss theoretical backgrounds and performance of our approach, in which elastic network normal mode analysis is used to predict reasonable conformational transitions from known X-ray structures, and find alternative conformations that are consistent with SAXS data.

A coarse-grained molecular model for actin-myosin simulation

We describe a very coarse-grained molecular model for the simulation of myosin V on an actin filament. The molecular representation is hierarchical with the finest level representing secondary structure elements (end-points) which are grouped into domains which are then grouped into molecules. Each level moves with a Brownian-like motion both in translation and rotation. Molecular integrity is maintained by steric exclusion and inter-domain restraints. A molecular description is developed for a myosin dimer on a actin filament with binding interactions also specified between domains to simulate both loose and tight binding. The stability of the model was tested in the pre- and post-power-stroke conformations with simulations in both states being used to test the preferred binding site of the myosin on the filament. The effects of the myosin twofold symmetry and the restriction of an attached cargo were also tested. These results provide the basis for the development of a dynamic model of processive motion.

Common structural motifs for the regulation of divergent class II myosins

This minireview focuses on structural studies that have provided insights into our current understanding of thick filament regulation in muscle. We describe how different domains in the myosin molecule interact to produce an inactive "off" state; included are head-head and head-rod interactions, the role of the regulatory light chain, and the significance of the alpha-helical coiled-coil rod in regulation. Several of these interactions have now been visualized in a wide variety of native myosin filaments, testifying to the generality of these structural motifs across the phylogenetic tree.

RNA polymerase II flexibility during translocation from normal mode analysis

The structural dynamics in eukaryotic RNA polymerase II (RNAPII) is described from computational normal mode analysis based on a series of crystal structures of pre- and post-translocated states with open and closed trigger loops. Conserved modes are identified that involve translocation of the nucleic acid complex coupled to motions of the enzyme, in particular in the clamp and jaw domains of RNAPII. A combination of these modes is hypothesized to be involved during active transcription. The NMA modes indicate furthermore that downstream DNA translocation may occur separately from DNA:RNA hybrid translocation. A comparison of the modes between different states of RNAPII suggests that productive translocation requires an open trigger loop and is inhibited by the presence of an NTP in the active site. This conclusion is also supported by a comparison of the overall flexibility in terms of root mean square fluctuations.

NMY-2 maintains cellular asymmetry and cell boundaries, and promotes a SRC-dependent asymmetric cell division

The nonmuscle myosin II NMY-2 is required for cytokinesis as well as for the establishment of zygote asymmetry during embryogenesis in Caenorhabditis elegans. Here we describe two conditional nmy-2 alleles that rapidly and reversibly inactivate the protein. We show that NMY-2 has late-cell-cycle roles in maintaining embryonic asymmetries and is also required for a surprisingly late step in the maintenance of the cytokinesis furrow. Finally, during a signaling-induced asymmetric cell division, NMY-2 is required for SRC-dependent phosphotyrosine signaling and acts in parallel with WNT-signaling to specify endoderm.

Coupling between normal modes drives protein conformational dynamics: illustrations using allosteric transitions in myosin II

Structure-based elastic network models (ENMs) have been remarkably successful in describing conformational transitions in a variety of biological systems. Low-frequency normal modes are usually calculated from the ENM that characterizes elastic interactions between residues in contact in a given protein structure with a uniform force constant. To explore the dynamical effects of nonuniform elastic interactions, we calculate the robustness and coupling of the low-frequency modes in the presence of nonuniform variations in the ENM force constant. The variations in the elastic interactions, approximated here by Gaussian noise, approximately account for perturbation effects of heterogeneous residue-residue interactions or evolutionary sequence changes within a protein family. First-order perturbation theory provides an efficient and qualitatively correct estimate of the mode robustness and mode coupling for finite perturbations to the ENM force constant. The mode coupling analysis and the mode robustness analysis identify groups of strongly coupled modes that encode for protein functional motions. We illustrate the new concepts using myosin II motor protein as an example. The biological implications of mode coupling in tuning the allosteric couplings among the actin-binding site, the nucleotide-binding site, and the force-generating converter and lever arm in myosin isoforms are discussed. We evaluate the robustness of the correlation functions that quantify the allosteric couplings among these three key structural motifs.

Head-head and head-tail interaction: a general mechanism for switching off myosin II activity in cells

Intramolecular interaction between myosin heads, blocking key sites involved in actin-binding and ATPase activity, appears to be a critical mechanism for switching off vertebrate smooth-muscle myosin molecules, leading to relaxation. We have tested the hypothesis that this interaction is a general mechanism for switching off myosin II-based motile activity in both muscle and nonmuscle cells. Electron microscopic images of negatively stained myosin II molecules were analyzed by single particle image processing. Molecules from invertebrate striated muscles with phosphorylation-dependent regulation showed head-head interactions in the off-state similar to those in vertebrate smooth muscle. A similar structure was observed in nonmuscle myosin II (also phosphorylation-regulated). Surprisingly, myosins from vertebrate skeletal and cardiac muscle, which are not intrinsically regulated, undergo similar head-head interactions in relaxing conditions. In all of these myosins, we also observe conserved interactions between the 'blocked' myosin head and the myosin tail, which may contribute to the switched-off state. These results suggest that intramolecular head-head and head-tail interactions are a general mechanism both for inducing muscle relaxation and for switching off myosin II-based motile activity in nonmuscle cells. These interactions are broken when myosin is activated.

Influence of oligomerization on the dynamics of G-protein coupled receptors as assessed by normal mode analysis

The recently discovered impact of oligomerization on G-protein coupled receptor (GPCR) function further complicates the already challenging goal of unraveling the molecular and dynamic mechanisms of these receptors. To help understand the effect of oligomerization on the dynamics of GPCRs, we have compared the motion of monomeric, dimeric, and tetrameric arrangements of the prototypic GPCR rhodopsin, using an approximate-yet powerful-normal mode analysis (NMA) technique termed elastic network model (ENM). Moreover, we have used ENM to discriminate between putative dynamic mechanisms likely to account for the recently observed conformational rearrangement of the TM4,5-TM4,5 dimerization interface of GPCRs that occurs upon activation. Our results indicate: (1) significant perturbation of the normal modes (NMs) of the rhodopsin monomer upon oligomerization, which is mainly manifested at interfacial regions; (2) increased positive correlation among the transmembrane domains (TMs) and between the extracellular loop (EL) and TM regions of the rhodopsin protomer; (3) highest interresidue positive correlation at the interfaces between protomers; and (4) experimentally testable hypotheses of differential motional changes within different putative oligomeric arrangements.

Langevin network model of myosin

Langevin mode theory and the coarse-grained elastic network model (ENM) for proteins are combined to yield the Langevin network model (LNM). Hydrodynamic radii of 6 A were assigned to each alpha-carbon on the basis of matching experimental translational and rotational diffusion constants of lysozyme, myoglobin, and hemoglobin with those calculated using a rigid body bead model with hydrodynamic interactions described by the Rotne-Prager tensor. LNM analysis of myosin II indicates that all ENM-like modes are overdamped at water viscosities. The low-frequency LNM modes in the pre-power stroke structure (PDB code: 1VOM) are substantially less mixed than the corresponding modes of the post-power stroke structure (1Q5G). Results from a four-bead model of the myosin "lever arm" indicate that coupling between modes increases as the array departs from linearity and are consistent with the results for 1VOM and 1Q5G. The decay times for all overdamped Langevin modes are shorter than the calculated rotational tumbling times found for lysozyme and myosin.

An unstable head-rod junction may promote folding into the compact off-state conformation of regulated myosins

The N-terminal region of myosin's rod-like subfragment 2 (S2) joins the two heads of this dimeric molecule and is key to its function. Previously, a crystal structure of this predominantly coiled-coil region was determined for a short fragment (51 residues plus a leucine zipper) of the scallop striated muscle myosin isoform. In that study, the N-terminal 10-14 residues were found to be disordered. We have now determined the structure of the same scallop peptide in three additional crystal environments. In each of two of these structures, improved order has allowed visualization of the entire N-terminus in one chain of the dimeric peptide. We have also compared the melting temperatures of this scallop S2 peptide with those of analogous peptides from three other isoforms. Taken together, these experiments, along with examination of sequences, point to a diminished stability of the N-terminal region of S2 in regulated myosins, compared with those myosins whose regulation is thin filament linked. It seems plain that this isoform-specific instability promotes the off-state conformation of the heads in regulated myosins. We also discuss how myosin isoforms with varied thermal stabilities share the basic capacity to transmit force efficiently in order to produce contraction in their on states.

The intrinsic dynamics and function of nickel-binding regulatory protein: insights from elastic network analysis

Nickel-responsive protein NikR regulates the nickel uptake in nickel-dependent bacteria by interacting with the operator of nikABCDE and subsequently repressing the transcription of NikABCDE, an ABC-type nickel transporter system. The function of NikR and its affinity for the operator DNA is highly conformation-dependent, which has been confirmed by three independent crystallographic studies on NikR proteins from different bacteria. Depending on the intracellular nickel concentration, NikR is able to adopt either the open form or one of the two closed forms (cis and trans) that differ in the domain-domain arrangement. Only the closed cis form is optimal for DNA binding. We examined the low-resolution vibrational spectrum of NikR in each conformational form using the elastic network model and observed large-scale domain-domain vibrations that are closely related to the conformational transitions required for function, particularly the symmetric bending mode and the asymmetric twisting mode. This analysis on the intrinsic dynamics coded in the three-dimensional molecular construct allows us to examine the proposed mechanisms of NikR regulation from the standpoint of protein collective motions. Our findings further support the three-state equilibrium hypothesis proposed by others, and imply that an isolated closed cis form may be dynamically unstable but can be stabilized by DNA binding. However, we also found that the simple C(alpha)-model used in the current analysis is insufficient to capture the impact of nickel binding on the protein dynamics, for which an all-atom model with detailed atom typing is more appropriate.

The B2 alternatively spliced isoform of nonmuscle myosin II-B lacks actin-activated MgATPase activity and in vitro motility

We report the initial biochemical characterization of an alternatively spliced isoform of nonmuscle heavy meromyosin (HMM) II-B2 and compare it with HMM II-B0, the nonspliced isoform. HMM II-B2 is the HMM derivative of an alternatively spliced isoform of endogenous nonmuscle myosin (NM) II-B, which has 21-amino acids inserted into loop 2, near the actin-binding region. NM II-B2 is expressed in the Purkinje cells of the cerebellum as well as in other neuronal cells [X. Ma, S. Kawamoto, J. Uribe, R.S. Adelstein, Function of the neuron-specific alternatively spliced isoforms of nonmuscle myosin II-B during mouse brain development, Mol. Biol. Cell 15 (2006) 2138-2149]. In contrast to any of the previously described isoforms of NM II (II-A, II-B0, II-B1, II-C0 and II-C1) or to smooth muscle myosin, the actin-activated MgATPase activity of HMM II-B2 is not significantly increased from a low, basal level by phosphorylation of the 20kDa myosin light chain (MLC-20). Moreover, although HMM II-B2 can bind to actin in the absence of ATP and is released in its presence, it cannot propel actin in the sliding actin filament assay following MLC-20 phosphorylation. Unlike HMM II-B2, the actin-activated MgATPase activity of a chimeric HMM with the 21-amino acid II-B2 sequence inserted into the homologous location in the heavy chain of HMM II-C is increased following MLC-20 phosphorylation. This indicates that the effect of the II-B2 insert is myosin heavy chain specific.

Normal-mode flexible fitting of high-resolution structure of biological molecules toward one-dimensional low-resolution data

We present a method for reconstructing a 3D structure from a pair distribution function by flexibly fitting known x-ray structures toward a conformation that agrees with the low-resolution data. This method uses a linear combination of low-frequency normal modes from elastic-network description of the molecule in an iterative manner to deform the structure optimally to conform to the target pair distribution function. A simple function, pair distance distribution function between atoms, is chosen as a test model to establish computational algorithms-optimization algorithm and scoring function-that can utilize low-resolution 1D data. To select a correct structural model based on less information, we developed a scoring function that takes into account a characteristic of pair distribution functions. In addition, we employ a new optimization algorithm, the trusted region method, that relies on both first and second derivatives of the scoring function. Illustrative results of our studies on simulated 1D data from five different proteins, for which large conformational changes are known to occur, are presented.

Structures of smooth muscle myosin and heavy meromyosin in the folded, shutdown state

Remodelling the contractile apparatus within smooth muscle cells allows effective contractile activity over a wide range of cell lengths. Thick filaments may be redistributed via depolymerisation into inactive myosin monomers that have been detected in vitro, in which the long tail has a folded conformation. Using negative stain electron microscopy of individual folded myosin molecules from turkey gizzard smooth muscle, we show that they are more compact than previously described, with heads and the three segments of the folded tail closely packed. Heavy meromyosin (HMM), which lacks two-thirds of the tail, closely resembles the equivalent parts of whole myosin. Image processing reveals a characteristic head region morphology for both HMM and myosin, with features identifiable by comparison with less compact molecules. The two heads associate asymmetrically: the tip of one motor domain touches the base of the other, resembling the blocked and free heads of this HMM when it forms 2D crystals on lipid monolayers. The tail of HMM lies between the heads, contacting the blocked motor domain, unlike in the 2D crystal. The tail of whole myosin is bent sharply and consistently close to residues 1175 and 1535. The first bend position correlates with a skip in the coiled coil sequence, the second does not. Tail segments 2 and 3 associate only with the blocked head, such that the second bend is near the C-lobe of the blocked head regulatory light chain. Quantitative analysis of tail flexibility shows that the single coiled coil of HMM has an apparent Young's modulus of about 0.5 GPa. The folded tail of the whole myosin is less flexible, indicating interactions between the segments. The folded tail does not modify the compact head arrangement but stabilises it, indicating a structural mechanism for the very low ATPase activity of the folded molecule.

Electron tomography of swollen rigor fibers of insect flight muscle reveals a short and variably angled S2 domain

Subfragment 2 (S2), the segment that links the two myosin heads to the thick filament backbone, may serve as a swing-out adapter allowing crossbridge access to actin, as the elastic component of crossbridges and as part of a phosphorylation-regulated on-off switch for crossbridges in smooth muscle. Low-salt expansion increases interfilament spacing (from 52 nm to 67 nm) of rigor insect flight muscle fibers and exposes a tethering segment of S2 in many crossbridges. Docking an actoS1 atomic model into EM tomograms of swollen rigor fibers identifies in situ for the first time the location, length and angle assignable to a segment of S2. Correspondence analysis of 1831 38.7 nm crossbridge repeats grouped self-similar forms from which class averages could be computed. The full range of the variability in angles and lengths of exposed S2 was displayed by using class averages for atomic fittings of acto-S1, while S2 was modeled by fitting a length of coiled-coil to unaveraged individual repeats. This hybrid modeling shows that the average length of S2 tethers along the thick filament (except near the tapered ends) is approximately 10 nm, or 16% of S2's total length, with an angular range encompassing 90 degrees axially and 120 degrees azimuthally. The large range of S2 angles indicates that some rigor bridges produce positive force that must be balanced by others producing drag force. The short tethering segment clarifies constraints on the function of S2 in accommodating variable myosin head access to actin. We suggest that the short length of S2 may also favor intermolecular head-head interactions in IFM relaxed thick filaments.

Modeling activated states of GPCRs: the rhodopsin template

Activation of G Protein-Coupled Receptors (GPCRs) is an allosteric mechanism triggered by ligand binding and resulting in conformational changes transduced by the transmembrane domain. Models of the activated forms of GPCRs have become increasingly necessary for the development of a clear understanding of signal propagation into the cell. Experimental evidence points to a multiplicity of conformations related to the activation of the receptor, rendered important physiologically by the suggestion that different conformations may be responsible for coupling to different signaling pathways. In contrast to the inactive state of rhodopsin (RHO) for which several high quality X-ray structures are available, the structure-related information for the active states of rhodopsin and all other GPCRs is indirect. We have collected and stored such information in a repository we maintain for activation-specific structural data available for rhodopsin-like GPCRs, http://www.physiology.med.cornell.edu/GPCRactivation/gpcrindex.html . Using these data as structural constraints, we have applied Simulated Annealing Molecular Dynamics to construct a number of different active state models of RHO starting from the known inactive structure. The common features of the models indicate that TM3 and TM5 play an important role in activation, in addition to the well-established rearrangement of TM6. Some of the structural changes observed in these models occur in regions that were not involved in the constraints, and have not been previously tested experimentally; they emerge as interesting candidates for further experimental exploration of the conformational space of activated GPCRs. We show that none of the normal modes calculated from the inactive structure has a dominant contribution along the path of conformational rearrangement from inactive to the active forms of RHO in the models. This result may differentiate rhodopsin from other GPCRs, and the reasons for this difference are discussed in the context of the structural properties and the physiological function of the protein.

SYMMETRY, FORM, AND SHAPE: Guiding Principles for Robustness in Macromolecular Machines

Computational studies of large macromolecular assemblages have come a long way during the past 10 years. With the explosion of computer power and parallel computing, timescales of molecular dynamics simulations have been extended far beyond the hundreds of picoseconds timescale. However, limitations remain for studies of large-scale conformational changes occurring on timescales beyond nanoseconds, especially for large macromolecules. In this review, we describe recent methods based on normal mode analysis that have enabled us to study dynamics on the microsecond timescale for large macromolecules using different levels of coarse graining, from atomically detailed models to those employing only low-resolution structural information. Emerging from such studies is a control principle for robustness in Nature's machines. We discuss this idea in the context of large-scale functional reorganization of the ribosome, virus particles, and the muscle protein myosin.

Structure and function of myosin filaments

Myosin filaments interact with actin to generate muscle contraction and many forms of cell motility. X-ray and electron microscopy (EM) studies have revealed the general organization of myosin molecules in relaxed filaments, but technical difficulties have prevented a detailed description. Recent studies using improved ultrastructural and image analysis techniques are overcoming these problems. Three-dimensional reconstructions using single-particle methods have provided many new insights into the organization of the myosin heads and tails. Docking of atomic structures into cryo-EM density maps suggests how regulated myosin filaments are 'switched off', bringing about muscle relaxation. Additionally, sequence analysis suggests probable interactions between myosin tails in the backbone, whereas crystallographic and EM studies are starting to reveal tail interactions directly in three dimensions.

Coarse-grained normal mode analysis in structural biology

The realization that experimentally observed functional motions of proteins can be predicted by coarse-grained normal mode analysis has renewed interest in applications to structural biology. Notable applications include the prediction of biologically relevant motions of proteins and supramolecular structures driven by their structure-encoded collective dynamics; the refinement of low-resolution structures, including those determined by cryo-electron microscopy; and the identification of conserved dynamic patterns and mechanically key regions within protein families. Additionally, hybrid methods that couple atomic simulations with deformations derived from coarse-grained normal mode analysis are able to sample collective motions beyond the range of conventional molecular dynamics simulations. Such applications have provided great insight into the underlying principles linking protein structures to their dynamics and their dynamics to their functions.

Structural basis for myosin V discrimination between distinct cargoes

Myosin V molecular motors move cargoes on actin filaments. A myosin V may move multiple cargoes to distinct places at different times. The cargoes attach to the globular tail of myosin V via cargo-specific receptors. Here we report the crystal structure at 2.2 A of the myosin V globular tail. The overall tertiary structure has not been previously observed. There are several patches of highly conserved regions distributed on the surface of the tail. These are candidate attachment sites for cargo-specific receptors. Indeed, we identified a region of five conserved surface residues that are solely required for vacuole inheritance. Likewise, we identified a region of five conserved surface residues that are required for secretory vesicle movement, but not vacuole movement. These two regions are at opposite ends of the oblong-shaped cargo-binding domain, and moreover are offset by 180 degrees. The fact that the cargo-binding areas are distant from each other and simultaneously exposed on the surface of the globular tail suggests that major targets for the regulation of cargo attachment are organelle-specific myosin V receptors.

Large amplitude conformational change in proteins explored with a plastic network model: adenylate kinase

The plastic network model (PNM) is used to generate a conformational change pathway for Escherichia coli adenylate kinase based on two crystal structures, namely that of an open and a closed conformer. In this model, the energy basins corresponding to known conformers are connected at their lowest common energies. The results are used to evaluate and analyze the minimal energy pathways between these basins. The open to closed transition analysis provides an identification of hinges that is in agreement with the existing definitions based on the available X-ray structures. The elastic energy distribution and the C(alpha) pseudo-dihedral variation provide similar information on these hinges. The ensemble of the 45 published structures for this protein and closely related proteins is shown to always be within 3.0 A of the pathway, which corresponds to a conformational change between two end structures that differ by a C(alpha)-atom root-mean-squared deviation of 7.1A.