Global conformational changes of ribosome observed by normal mode fitting for 3D Cryo-EM structures

The energetics of subunit rotation in the ribosome

Protein synthesis in the cell is controlled by an elaborate sequence of conformational rearrangements in the ribosome. The composition of a ribosome varies by species, though they typically contain ∼ 50-100 RNA and protein molecules. While advances in structural techniques have revolutionized our understanding of long-lived conformational states, a vast range of transiently visited configurations can not be directly observed. In these cases, computational/simulation methods can be used to understand the mechanical properties of the ribosome. Insights from these approaches can then help guide next-generation experimental measurements. In this short review, we discuss theoretical strategies that have been deployed to quantitatively describe the energetics of collective rearrangements in the ribosome. We focus on efforts to probe large-scale subunit rotation events, which involve the coordinated displacement of large numbers of atoms (tens of thousands). These investigations are revealing how the molecular structure of the ribosome encodes the mechanical properties that control large-scale dynamics.

State-dependent sequential allostery exhibited by chaperonin TRiC/CCT revealed by network analysis of Cryo-EM maps

The eukaryotic chaperonin TRiC/CCT plays a major role in assisting the folding of many proteins through an ATP-driven allosteric cycle. Recent structures elucidated by cryo-electron microscopy provide a broad view of the conformations visited at various stages of the chaperonin cycle, including a sequential activation of its subunits in response to nucleotide binding. But we lack a thorough mechanistic understanding of the structure-based dynamics and communication properties that underlie the TRiC/CCT machinery. In this study, we present a computational methodology based on elastic network models adapted to cryo-EM density maps to gain a deeper understanding of the structure-encoded allosteric dynamics of this hexadecameric machine. We have analysed several structures of the chaperonin resolved in different states toward mapping its conformational landscape. Our study indicates that the overall architecture intrinsically favours cooperative movements that comply with the structural variabilities observed in experiments. Furthermore, the individual subunits CCT1-CCT8 exhibit state-dependent sequential events at different states of the allosteric cycle. For example, in the ATP-bound state, subunits CCT5 and CCT4 selectively initiate the lid closure motions favoured by the overall architecture; whereas in the apo form of the heteromer, the subunit CCT7 exhibits the highest predisposition to structural change. The changes then propagate through parallel fluxes of allosteric signals to neighbours on both rings. The predicted state-dependent mechanisms of sequential activation provide new insights into TRiC/CCT intra- and inter-ring signal transduction events.

Normal mode analysis calculation for a full-atom model with a smaller number of degrees of freedom for huge protein molecules

The number of degrees of freedom (DOF), N, in normal mode analysis (NMA) calculations of proteins is a crucial problem in huge systems because the eigenvalue problem of an N-by-N matrix must be solved. If it were possible to perform the analysis with a smaller number of DOF for the same system with minimal deterioration in accuracy, this would make a significant impact on the computational study of protein dynamics. We examined two models in which the number of DOF was reduced. Both of them adopted a full-atom model with dihedral angles as independent variables. In one model, side-chain dihedral angles, χ's, and a main-chain dihedral angle, ω, were fixed and only the main-chain dihedral angles, ϕ and ψ, were variable. In another model, the dihedral angles around virtual bonds that connect neighboring Cα atoms were tested. The number of DOFs for the two models was two and one per residue, respectively. The residue-by-residue fluctuation profiles for atoms and dihedral angles were well reproduced in both models. The motion of atoms in the individual lowest-frequency normal modes of the two models was also very similar to those of the original model in which all rotatable dihedral angles were variable. Consequently, these models could predict large-amplitude concerted motion. These results also imply that proteins in a full-atom model can undergo only limited large-scale conformational changes around the native conformation, and consequently, NMA results do not strongly depend on the independent variables adopted.

Resolution-exchanged structural modeling and simulations jointly unravel that subunit rolling underlies the mechanism of programmed ribosomal frameshifting

MOTIVATION

Programmed ribosomal frameshifting (PRF) is widely used by viruses and bacteria to produce different proteins from a single mRNA template. How steric hindrance of a PRF-stimulatory mRNA structure transiently modifies the conformational dynamics of the ribosome, and thereby allows tRNA slippage, remains elusive.

RESULTS

Here, we leverage linear response theories and resolution-exchanged simulations to construct a structural/dynamics model that connects and rationalizes existing structural, single-molecule and mutagenesis data by resolution-exchanged structural modelling and simulations. Our combined theoretical techniques provide a temporal and spatial description of PRF with unprecedented mechanistic details. We discover that ribosomal unfolding of the PRF-stimulating pseudoknot exerts resistant forces on the mRNA entrance of the ribosome, and thereby drives 30S subunit rolling. Such motion distorts tRNAs, leads to tRNA slippage, and in turn serves as a delicate control of cis-element's unwinding forces over PRF.

AVAILABILITY AND IMPLEMENTATION

All the simulation scripts and computational implementations of our methods/analyses (including linear response theory) are included in the bioStructureM suite, provided through GitHub at https://github.com/Yuan-Yu/bioStructureM.

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.

Dynamic analysis of ribosome by a movie made from many three-dimensional electron-microscopy density maps

The atomic models of the 70S ribosome including the bound molecules were built from many 3D-EM density maps. The positions and conformations of the bound molecules were determined by fitting them to the regions in the density maps which remained after fitting the 70S ribosome. Then, using these atomic models, a movie for the elongation cycle was made. For determining the sequential order in which the models appeared in the movie, the knowledge about the bound molecules and the ratchet angles were used. The movie revealed several interesting points which were not apparent from each density map, suggesting the usefulness of a movie made from many 3D-EM density maps.

2DHybrid Analysis

We have developed an approach termed '2D hybrid analysis' for building three-dimensional (3D) structures from electron microscopy (EM) images of biological molecules. The key advantage is that it is applicable to flexible molecules, which are difficult to analyze by the approach in which 3DEM maps are reconstructed. In the proposed approach, a large number of atomic models with different conformations are first built by computer simulation. Then, simulated EM images are produced from each atomic model. Finally, these images are compared with an experimental EM image to identify the best-fitting atomic model. Two kinds of models are used to simulate the EM images: the negative-stain model and the simple projection model. Although the former is more realistic, the latter permits faster computation. We applied this approach to the averaged EM images of integrin. Although many of these were reproduced well by the best-fitting atomic models, others did not closely resemble any of the simulated EM images. However, the latter group were well reproduced by averaging multiple simulated EM images originating from atomic models with rather different conformations or orientations. This indicated that our approach is capable of detecting mixtures of conformations in the averaged EM images, which should assist in their correct interpretation.

Acceleration of cryo-EM Flexible Fitting for Large Biomolecular Systems by Efficient Space Partitioning

Flexible fitting is a powerful technique to build the 3D structures of biomolecules from cryoelectron microscopy (cryo-EM) density maps. One popular method is a cross-correlation coefficient-based approach, where the molecular dynamics (MD) simulation is carried out with the biasing potential that includes the cross-correlation coefficient between the experimental and simulated density maps. Here, we propose efficient parallelization schemes for the calculation of the cross-correlation coefficient to accelerate flexible fitting. Our schemes are tested for small, medium, and large biomolecules using CPU and hybrid CPU + GPU architectures. The scheme for the atomic decomposition MD is suitable for small proteins such as Ca²⁺-ATPase with the all-atom Go model, while that for the domain decomposition MD is better for larger systems such as ribosome with the all-atom Go or the all-atom explicit solvent models. Our methods allow flexible fitting for various biomolecules with reasonable computational cost. This approach also connects high-resolution structure refinements with investigation of protein structure-function relationship.

Quantitative Characterization of Domain Motions in Molecular Machines

The work of molecular machines such as the ribosome is accompanied by conformational changes, often characterized by relative motions of their domains. The method we have developed seeks to quantify these motions in a general way, facilitating comparisons of results obtained by different researchers. Typically there are multiple snapshots of a structure in the form of pdb coordinates resulting from flexible fitting of low-resolution density maps, from X-ray crystallography, or from molecular dynamics simulation trajectories. Our objective is to characterize the motion of each domain as a coordinate transformation using moments of inertia tensor, a method we developed earlier. What has been missing until now are ancillary tools that make this task practical, general, and biologically informative. We have provided a comprehensive solution to this task with a set of tools implemented on the VMD platform. These tools address the need for reproducible segmentation of domains, and provide a generalized description of their motions using principal axes of inertia. Although this methodology has been specifically developed for studying ribosome motion, it is applicable to any molecular machine.

2D hybrid analysis: Approach for building three-dimensional atomic model by electron microscopy image matching

In this study, we develop an approach termed “2D hybrid analysis” for building atomic models by image matching from electron microscopy (EM) images of biological molecules. The key advantage is that it is applicable to flexible molecules, which are difficult to analyze by 3DEM approach. In the proposed approach, first, a lot of atomic models with different conformations are built by computer simulation. Then, simulated EM images are built from each atomic model. Finally, they are compared with the experimental EM image. Two kinds of models are used as simulated EM images: the negative stain model and the simple projection model. Although the former is more realistic, the latter is adopted to perform faster computations. The use of the negative stain model enables decomposition of the averaged EM images into multiple projection images, each of which originated from a different conformation or orientation. We apply this approach to the EM images of integrin to obtain the distribution of the conformations, from which the pathway of the conformational change of the protein is deduced.

New generation of elastic network models

The intrinsic flexibility of proteins and nucleic acids can be grasped from remarkably simple mechanical models of particles connected by springs. In recent decades, Elastic Network Models (ENMs) combined with Normal Model Analysis widely confirmed their ability to predict biologically relevant motions of biomolecules and soon became a popular methodology to reveal large-scale dynamics in multiple structural biology scenarios. The simplicity, robustness, low computational cost, and relatively high accuracy are the reasons behind the success of ENMs. This review focuses on recent advances in the development and application of ENMs, paying particular attention to combinations with experimental data. Successful application scenarios include large macromolecular machines, structural refinement, docking, and evolutionary conservation.

Functional Importance of Mobile Ribosomal Proteins

Although the dynamic motions and peptidyl transferase activity seem to be embedded in the rRNAs, the ribosome contains more than 50 ribosomal proteins (r-proteins), whose functions remain largely elusive. Also, the precise forms of some of these r-proteins, as being part of the ribosome, are not structurally solved due to their high flexibility, which hinders the efforts in their functional elucidation. Owing to recent advances in cryo-electron microscopy, single-molecule techniques, and theoretical modeling, much has been learned about the dynamics of these r-proteins. Surprisingly, allosteric regulations have been found in between spatially separated components as distant as those in the opposite sides of the ribosome. Here, we focus on the functional roles and intricate regulations of the mobile L1 and L12 stalks and L9 and S1 proteins. Conformational flexibility also enables versatile functions for r-proteins beyond translation. The arrangement of r-proteins may be under evolutionary pressure that fine-tunes mass distributions for optimal structural dynamics and catalytic activity of the ribosome.

Giant cadherins Fat and Dachsous self-bend to organize properly spaced intercellular junctions

The cadherins Fat and Dachsous regulate cell polarity and proliferation via their heterophilic interactions at intercellular junctions. Their ectodomains are unusually large because of repetitive extracellular cadherin (EC) domains, which raises the question of how they fit in regular intercellular spaces. Cadherins typically exhibit a linear topology through the binding of Ca(2+) to the linker between the EC domains. Our electron-microscopic observations of mammalian Fat4 and Dachsous1 ectodomains, however, revealed that, although their N-terminal regions exhibit a linear configuration, the C-terminal regions are kinked with multiple hairpin-like bends. Notably, certain EC-EC linkers in Fat4 and Dachsous1 lost Ca(2+)-binding amino acids. When such non-Ca(2+)-binding linkers were substituted for a normal linker in E-cadherin, the mutant E-cadherins deformed more extensively than the wild-type molecule. To simulate cadherin structures with non-Ca(2+)-binding linkers, we used an elastic network model and confirmed that bent configurations can be generated by deformation of non-Ca(2+)-binding linkers. These findings suggest that Fat and Dachsous self-bend due to the loss of Ca(2+)-binding amino acids from specific EC-EC linkers, and can therefore adapt to confined spaces.

Free-energy landscape of reverse tRNA translocation through the ribosome analyzed by electron microscopy density maps and molecular dynamics simulations

To understand the mechanism of reverse tRNA translocation in the ribosome, all-atom molecular dynamics simulations of the ribosome-tRNAs-mRNA-EFG complex were performed. The complex at the post-translocational state was directed towards the translocational and pre-translocational states by fitting the complex into cryo-EM density maps. Between a series of the fitting simulations, umbrella sampling simulations were performed to obtain the free-energy landscape. Multistep structural changes, such as a ratchet-like motion and rotation of the head of the small subunit were observed. The free-energy landscape showed that there were two main free-energy barriers: one between the post-translocational and intermediate states, and the other between the pre-translocational and intermediate states. The former corresponded to a clockwise rotation, which was coupled to the movement of P-tRNA over the P/E-gate made of G1338, A1339 and A790 in the small subunit. The latter corresponded to an anticlockwise rotation of the head, which was coupled to the location of the two tRNAs in the hybrid state. This indicates that the coupled motion of the head rotation and tRNA translocation plays an important role in opening and closing of the P/E-gate during the ratchet-like movement in the ribosome. Conformational change of EF-G was interpreted to be the result of the combination of the external motion by L12 around an axis passing near the sarcin-ricin loop, and internal hinge-bending motion. These motions contributed to the movement of domain IV of EF-G to maintain its interaction with A/P-tRNA.

Structured mRNA induces the ribosome into a hyper-rotated state

During protein synthesis, mRNA and tRNA are moved through the ribosome by the process of translocation. The small diameter of the mRNA entrance tunnel only permits unstructured mRNA to pass through. However, there are structured elements within mRNA that present a barrier for translocation that must be unwound. The ribosome has been shown to unwind RNA in the absence of additional factors, but the mechanism remains unclear. Here, we show using single molecule Förster resonance energy transfer and small angle X-ray scattering experiments a new global conformational state of the ribosome. In the presence of the frameshift inducing dnaX hairpin, the ribosomal subunits are driven into a hyper-rotated state and the L1 stalk is predominantly in an open conformation. This previously unobserved conformational state provides structural insight into the helicase activity of the ribosome and may have important implications for understanding the mechanism of reading frame maintenance.

Fast fitting to low resolution density maps: elucidating large-scale motions of the ribosome

Determining the conformational rearrangements of large macromolecules is challenging experimentally and computationally. Case in point is the ribosome; it has been observed by high-resolution crystallography in several states, but many others are known only from low-resolution methods including cryo-electron microscopy. Combining these data into dynamical trajectories that may aid understanding of its largest-scale conformational changes has so far remained out of reach of computational methods. Most existing methods either model all atoms explicitly, resulting in often prohibitive cost, or use approximations that lose interesting structural and dynamical detail. In this work, I introduce Internal Coordinate Flexible Fitting, which uses full atomic forces and flexibility in limited regions of a model, capturing extensive conformational rearrangements at low cost. I use it to turn multiple low-resolution density maps, crystallographic structures and biochemical information into unified all-atoms trajectories of ribosomal translocation. Internal Coordinate Flexible Fitting is three orders of magnitude faster than the most comparable existing method.

Revealing -1 programmed ribosomal frameshifting mechanisms by single-molecule techniques and computational methods

Programmed ribosomal frameshifting (PRF) serves as an intrinsic translational regulation mechanism employed by some viruses to control the ratio between structural and enzymatic proteins. Most viral mRNAs which use PRF adapt an H-type pseudoknot to stimulate −1 PRF. The relationship between the thermodynamic stability and the frameshifting efficiency of pseudoknots has not been fully understood. Recently, single-molecule force spectroscopy has revealed that the frequency of −1 PRF correlates with the unwinding forces required for disrupting pseudoknots, and that some of the unwinding work dissipates irreversibly due to the torsional restraint of pseudoknots. Complementary to single-molecule techniques, computational modeling provides insights into global motions of the ribosome, whose structural transitions during frameshifting have not yet been elucidated in atomic detail. Taken together, recent advances in biophysical tools may help to develop antiviral therapies that target the ubiquitous −1 PRF mechanism among viruses.

Solution structure of RNase P RNA

The ribonucleoprotein enzyme ribonuclease P (RNase P) processes tRNAs by cleavage of precursor-tRNAs. RNase P is a ribozyme: The RNA component catalyzes tRNA maturation in vitro without proteins. Remarkable features of RNase P include multiple turnovers in vivo and ability to process diverse substrates. Structures of the bacterial RNase P, including full-length RNAs and a ternary complex with substrate, have been determined by X-ray crystallography. However, crystal structures of free RNA are significantly different from the ternary complex, and the solution structure of the RNA is unknown. Here, we report solution structures of three phylogenetically distinct bacterial RNase P RNAs from Escherichia coli, Agrobacterium tumefaciens, and Bacillus stearothermophilus, determined using small angle X-ray scattering (SAXS) and selective 2'-hydroxyl acylation analyzed by primer extension (SHAPE) analysis. A combination of homology modeling, normal mode analysis, and molecular dynamics was used to refine the structural models against the empirical data of these RNAs in solution under the high ionic strength required for catalytic activity.

Unified data resource for cryo-EM

Three-dimensional (3D) cryoelectron microscopy reconstruction methods are uniquely able to reveal structures of many important macromolecules and macromolecular complexes. EMDataBank.org, a joint effort of the Protein Databank in Europe (PDBe), the Research Collaboratory for Structural Bioinformatics (RCSB), and the National Center for Macromolecular Imaging (NCMI), is a "one-stop shop" resource for global deposition and retrieval of cryo-EM map, model, and associated metadata. The resource unifies public access to the two major EM Structural Data archives: EM Data Bank (EMDB) and Protein Data Bank (PDB), and facilitates use of EM structural data of macromolecules and macromolecular complexes by the wider scientific community.

EMDataBank.org: unified data resource for CryoEM

Cryo-electron microscopy reconstruction methods are uniquely able to reveal structures of many important macromolecules and macromolecular complexes. EMDataBank.org, a joint effort of the Protein Data Bank in Europe (PDBe), the Research Collaboratory for Structural Bioinformatics (RCSB) and the National Center for Macromolecular Imaging (NCMI), is a global ‘one-stop shop’ resource for deposition and retrieval of cryoEM maps, models and associated metadata. The resource unifies public access to the two major archives containing EM-based structural data: EM Data Bank (EMDB) and Protein Data Bank (PDB), and facilitates use of EM structural data of macromolecules and macromolecular complexes by the wider scientific community.