RIBOSOME DYNAMICS: Insights from Atomic Structure Modeling into Cryo-Electron Microscopy Maps

Biomotors, viral assembly, and RNA nanobiotechnology: Current achievements and future directions

The International Society of RNA Nanotechnology and Nanomedicine (ISRNN) serves to further the development of a wide variety of functional nucleic acids and other related nanotechnology platforms. To aid in the dissemination of the most recent advancements, a biennial discussion focused on biomotors, viral assembly, and RNA nanobiotechnology has been established where international experts in interdisciplinary fields such as structural biology, biophysical chemistry, nanotechnology, cell and cancer biology, and pharmacology share their latest accomplishments and future perspectives. The results summarized here highlight advancements in our understanding of viral biology and the structure-function relationship of frame-shifting elements in genomic viral RNA, improvements in the predictions of SHAPE analysis of 3D RNA structures, and the understanding of dynamic RNA structures through a variety of experimental and computational means. Additionally, recent advances in the drug delivery, vaccine design, nanopore technologies, biomotor and biomachine development, DNA packaging, RNA nanotechnology, and drug delivery are included in this critical review. We emphasize some of the novel accomplishments, major discussion topics, and present current challenges and perspectives of these emerging fields.

The dynamic, motile and deformative properties of RNA nanoparticles facilitate the third milestone of drug development

Besides mRNA, rRNA, and tRNA, cells contain many other noncoding RNA that display critical roles in the regulation of cellular functions. Human genome sequencing revealed that the majority of non-protein-coding DNA actually codes for non-coding RNAs. The dynamic nature of RNA results in its motile and deformative behavior. These conformational transitions such as the change of base-pairing, breathing within complemented strands, and pseudoknot formation at the 2D level as well as the induced-fit and conformational capture at the 3D level are important for their biological functions including regulation, translation, and catalysis. The dynamic, motile and catalytic activity has led to a belief that RNA is the origin of life. We have recently reported that the deformative property of RNA nanoparticles enhances their penetration through the leaky blood vessel of cancers which leads to highly efficient tumor accumulation. This special deformative property also enables RNA nanoparticles to pass the glomerulus, overcoming the filtration size limit, resulting in fast renal excretion and rapid body clearance, thus low or no toxicity. The biodistribution of RNA nanoparticles can be further improved by the incorporation of ligands for cancer targeting. In addition to the favorable biodistribution profiles, RNA nanoparticles possess other properties including self-assembly, negative charge, programmability, and multivalency; making it a great material for pharmaceutical applications. The intrinsic negative charge of RNA nanoparticles decreases the toxicity of drugs by preventing nonspecific binding to the negative charged cell membrane and enhancing the solubility of hydrophobic drugs. The polyvalent property of RNA nanoparticles allows the multi-functionalization which can apply to overcome drug resistance. This review focuses on the summary of these unique properties of RNA nanoparticles, which describes the mechanism of RNA dynamic, motile and deformative properties, and elucidates and prepares to welcome the RNA therapeutics as the third milestone in pharmaceutical drug development.

Probing Structural Perturbation of Biomolecules by Extracting Cryo-EM Data Heterogeneity

Single-particle cryogenic electron microscopy (cryo-EM) has become an indispensable tool to probe high-resolution structural detail of biomolecules. It enables direct visualization of the biomolecules and opens a possibility for averaging molecular images to reconstruct a three-dimensional Coulomb potential density map. Newly developed algorithms for data analysis allow for the extraction of structural heterogeneity from a massive and low signal-to-noise-ratio (SNR) cryo-EM dataset, expanding our understanding of multiple conformational states, or further implications in dynamics, of the target biomolecule. This review provides an overview that briefly describes the workflow of single-particle cryo-EM, including imaging and data processing, and new methods developed for analyzing the data heterogeneity to understand the structural variability of biomolecules.

Can DyeCycling break the photobleaching limit in single-molecule FRET?

Biomolecular systems, such as proteins, crucially rely on dynamic processes at the nanoscale. Detecting biomolecular nanodynamics is therefore key to obtaining a mechanistic understanding of the energies and molecular driving forces that control biomolecular systems. Single-molecule fluorescence resonance energy transfer (smFRET) is a powerful technique to observe in real-time how a single biomolecule proceeds through its functional cycle involving a sequence of distinct structural states. Currently, this technique is fundamentally limited by irreversible photobleaching, causing the untimely end of the experiment and thus, a narrow temporal bandwidth of ≤ 3 orders of magnitude. Here, we introduce "DyeCycling", a measurement scheme with which we aim to break the photobleaching limit in smFRET. We introduce the concept of spontaneous dye replacement by simulations, and as an experimental proof-of-concept, we demonstrate the intermittent observation of a single biomolecule for one hour with a time resolution of milliseconds. Theoretically, DyeCycling can provide > 100-fold more information per single molecule than conventional smFRET. We discuss the experimental implementation of DyeCycling, its current and fundamental limitations, and specific biological use cases. Given its general simplicity and versatility, DyeCycling has the potential to revolutionize the field of time-resolved smFRET, where it may serve to unravel a wealth of biomolecular dynamics by bridging from milliseconds to the hour range.

ELECTRONIC SUPPLEMENTARY MATERIAL

Supplementary material is available for this article at 10.1007/s12274-022-4420-5 and is accessible for authorized users.

Ribosome-membrane crosstalk: Co-translational targeting pathways of proteins across membranes in prokaryotes and eukaryotes

Ribosomes are the molecular machine of living cells designed for decoding mRNA-encoded genetic information into protein. Being sophisticated machinery, both in design and function, the ribosome not only carries out protein synthesis, but also coordinates several other ribosome-associated cellular processes. One such process is the translocation of proteins across or into the membrane depending on their secretory or membrane-associated nature. These proteins comprise a large portion of a cell's proteome and act as key factors for cellular survival as well as several crucial functional pathways. Protein transport to extra- and intra-cytosolic compartments (across the eukaryotic endoplasmic reticulum (ER) or across the prokaryotic plasma membrane) or insertion into membranes majorly occurs through an evolutionarily conserved protein-conducting channel called translocon (eukaryotic Sec61 or prokaryotic SecYEG channels). Targeting proteins to the membrane-bound translocon may occur via post-translational or co-translational modes and it is often mediated by recognition of an N-terminal signal sequence in the newly synthesizes polypeptide chain. Co-translational translocation is coupled to protein synthesis where the ribosome-nascent chain complex (RNC) itself is targeted to the translocon. Here, in the light of recent advances in structural and functional studies, we discuss our current understanding of the mechanistic models of co-translational translocation, coordinated by the actively translating ribosomes, in prokaryotes and eukaryotes.

DR-SIP: protocols for higher order structure modeling with distance restraints- and cyclic symmetry-imposed packing

MOTIVATION

Quaternary structure determination for transmembrane/soluble proteins requires a reliable computational protocol that leverages observed distance restraints and/or cyclic symmetry (Cn symmetry) found in most homo-oligomeric transmembrane proteins.

RESULTS

We survey 118 X-ray crystallographically solved structures of homo-oligomeric transmembrane proteins (HoTPs) and find that ∼97% are Cn symmetric. Given the prevalence of Cn symmetric HoTPs and the benefits of incorporating geometry restraints in aiding quaternary structure determination, we introduce two new filters, the distance-restraints (DR) and the Symmetry-Imposed Packing (SIP) filters. SIP relies on a new method that can rebuild the closest ideal Cn symmetric complex from docking poses containing a homo-dimer without prior knowledge of the number (n) of monomers. Using only the geometrical filter, SIP, near-native poses of 7 HoTPs in their monomeric states can be correctly identified in the top-10 for 71% of all cases, or 29% among 31 HoTP structures obtained through homology modeling, while ZDOCK alone returns 14 and 3%, respectively. When the n is given, the optional n-mer filter is applied with SIP and returns the near-native poses for 76% of the test set within the top-10, outperforming M-ZDOCK's 55% and Sam's 47%. While applying only SIP to three HoTPs that comes with distance restraints, we found the near-native poses were ranked 1st, 1st and 10th among 54 000 possible decoys. The results are further improved to 1st, 1st and 3rd when both DR and SIP filters are used. By applying only DR, a soluble system with distance restraints is recovered at the 1st-ranked pose.

AVAILABILITY AND IMPLEMENTATION

https://github.com/capslockwizard/drsip.

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.

Modeling Biological Complexes Using Integrative Modeling Platform

Integrative structure modeling provides 3D models of macromolecular systems that are based on information from multiple types of experiments, physical principles, statistical inferences, and prior structural models. Here, we provide a hands-on realistic example of integrative structure modeling of the quaternary structure of the actin, tropomyosin, and gelsolin protein assembly based on electron microscopy, solution X-ray scattering, and chemical crosslinking data for the complex as well as excluded volume, sequence connectivity, and rigid atomic X-ray structures of the individual subunits. We follow the general four-stage process for integrative modeling, including gathering the input information, converting the input information into a representation of the system and a scoring function, sampling alternative model configurations guided by the scoring function, and analyzing the results. The computational aspects of this approach are implemented in our open-source Integrative Modeling Platform (IMP), a comprehensive and extensible software package for integrative modeling ( https://integrativemodeling.org ). In particular, we rely on the Python Modeling Interface (PMI) module of IMP that provides facile mixing and matching of macromolecular representations, restraints based on different types of information, sampling algorithms, and analysis including validations of the input data and output models. Finally, we also outline how to deposit an integrative structure and corresponding experimental data into PDB-Dev, the nascent worldwide Protein Data Bank (wwPDB) resource for archiving and disseminating integrative structures ( https://pdb-dev.wwpdb.org ). The example application provides a starting point for a user interested in using IMP for integrative modeling of other biomolecular systems.

Analytical Approaches to Improve Accuracy in Solving the Protein Topology Problem

To take advantage of recent advances in genomics and proteomics it is critical that the three-dimensional physical structure of biological macromolecules be determined. Cryo-Electron Microscopy (cryo-EM) is a promising and improving method for obtaining this data, however resolution is often not sufficient to directly determine the atomic scale structure. Despite this, information for secondary structure locations is detectable. De novo modeling is a computational approach to modeling these macromolecular structures based on cryo-EM derived data. During de novo modeling a mapping between detected secondary structures and the underlying amino acid sequence must be identified. DP-TOSS (Dynamic Programming for determining the Topology Of Secondary Structures) is one tool that attempts to automate the creation of this mapping. By treating the correspondence between the detected structures and the structures predicted from sequence data as a constraint graph problem DP-TOSS achieved good accuracy in its original iteration. In this paper, we propose modifications to the scoring methodology of DP-TOSS to improve its accuracy. Three scoring schemes were applied to DP-TOSS and tested: (i) a skeleton-based scoring function; (ii) a geometry-based analytical function; and (iii) a multi-well potential energy-based function. A test of 25 proteins shows that a combination of these schemes can improve the performance of DP-TOSS to solve the topology determination problem for macromolecule proteins.

PEM-fitter: A Coarse-Grained Method to Validate Protein Candidate Models

The volumetric images produced by Cryo-Electron Microscopy (cryo-EM) technique are used to model macromolecular assemblies and machines. De novo protein modeling uses these images to computationally model the structure of the molecules. Many candidate conformations are usually generated during the intermediate step. Conventionally, each of these candidates is evaluated by time-consuming approaches such as potential energy. We introduce an initial version of a geometrical screening method that uses the skeleton of the cryo-EM images to evaluate candidate structures. The aim of this method is to reduce the number of native-like candidate conformations and, therefore, reduce the time required for structural evaluation by energy calculations. A test of two datasets was performed. The first dataset contains 10 proteins and shows that our method can successfully detect the correct native structure for the given skeleton among a set of different protein structures. The second dataset contains 12 proteins and shows that our method can filter slightly modified decoy conformations of the same protein. The efficiency of the method is also reported.

Integrative structure modeling with the Integrative Modeling Platform

Building models of a biological system that are consistent with the myriad data available is one of the key challenges in biology. Modeling the structure and dynamics of macromolecular assemblies, for example, can give insights into how biological systems work, evolved, might be controlled, and even designed. Integrative structure modeling casts the building of structural models as a computational optimization problem, for which information about the assembly is encoded into a scoring function that evaluates candidate models. Here, we describe our open source software suite for integrative structure modeling, Integrative Modeling Platform (https://integrativemodeling.org), and demonstrate its use.

The impact of recent improvements in cryo-electron microscopy technology on the understanding of bacterial ribosome assembly

Cryo-electron microscopy (cryo-EM) had played a central role in the study of ribosome structure and the process of translation in bacteria since the development of this technique in the mid 1980s. Until recently cryo-EM structures were limited to ∼10 Å in the best cases. However, the recent advent of direct electron detectors has greatly improved the resolution of cryo-EM structures to the point where atomic resolution is now achievable. This improved resolution will allow cryo-EM to make groundbreaking contributions in essential aspects of ribosome biology, including the assembly process. In this review, we summarize important insights that cryo-EM, in combination with chemical and genetic approaches, has already brought to our current understanding of the ribosomal assembly process in bacteria using previous detector technology. More importantly, we discuss how the higher resolution structures now attainable with direct electron detectors can be leveraged to propose precise testable models regarding this process. These structures will provide an effective platform to develop new antibiotics that target this fundamental cellular process.

Use of evolutionary information in the fitting of atomic level protein models in low resolution cryo-EM map of a protein assembly improves the accuracy of the fitting

Protein-protein interface residues, especially those at the core of the interface, exhibit higher conservation than residues in solvent exposed regions. Here, we explore the ability of this differential conservation to evaluate fittings of atomic models in low-resolution cryo-EM maps and select models from the ensemble of solutions that are often proposed by different model fitting techniques. As a prelude, using a non-redundant and high-resolution structural dataset involving 125 permanent and 95 transient complexes, we confirm that core interface residues are conserved significantly better than nearby non-interface residues and this result is used in the cryo-EM map analysis. From the analysis of inter-component interfaces in a set of fitted models associated with low-resolution cryo-EM maps of ribosomes, chaperones and proteasomes we note that a few poorly conserved residues occur at interfaces. Interestingly a few conserved residues are not in the interface, though they are close to the interface. These observations raise the potential requirement of refitting the models in the cryo-EM maps. We show that sampling an ensemble of models and selection of models with high residue conservation at the interface and in good agreement with the density helps in improving the accuracy of the fit. This study indicates that evolutionary information can serve as an additional input to improve and validate fitting of atomic models in cryo-EM density maps.

Ribosome Mechanics Informs about Mechanism

The essential aspects of the ribosome's mechanism can be extracted from coarse-grained simulations, including the ratchet motion, the movement together of critical bases at the decoding center, and movements of the peptide tunnel lining that assist in the expulsion of the synthesized peptide. Because of its large size, coarse graining helps to simplify and to aid in the understanding of its mechanism. Results presented here utilize coarse-grained elastic network modeling to extract the dynamics, and both RNAs and proteins are coarse grained. We review our previous results, showing the well-known ratchet motions and the motions in the peptide tunnel and in the mRNA tunnel. The motions of the lining of the peptide tunnel appear to assist in the expulsion of the growing peptide chain, and clamps at the ends of the mRNA tunnel with three proteins ensure that the mRNA is held tightly during decoding and essential for the helicase activity at the entrance. The entry clamp may also assist in base recognition to ensure proper selection of the incoming tRNA. The overall precision of the ribosome machine-like motions is remarkable.

Ribosomal protein uS19 mutants reveal its role in coordinating ribosome structure and function

Prior studies identified allosteric information pathways connecting functional centers in the large ribosomal subunit to the decoding center in the small subunit through the B1a and B1b/c intersubunit bridges in yeast. In prokaryotes a single SSU protein, uS13, partners with H38 (the A-site finger) and uL5 to form the B1a and B1b/c bridges respectively. In eukaryotes, the SSU component was split into 2 separate proteins during the course of evolution. One, also known as uS13, participates in B1b/c bridge with uL5 in eukaryotes. The other, called uS19 is the SSU partner in the B1a bridge with H38. Here, polyalanine mutants of uS19 involved in the uS19/uS13 and the uS19/H38 interfaces were used to elucidate the important amino acid residues involved in these intersubunit communication pathways. Two key clusters of amino acids were identified: one located at the junction between uS19 and uS13, and a second that appears to interact with the distal tip of H38. Biochemical analyses reveal that these mutations shift the ribosomal rotational equilibrium toward the unrotated state, increasing ribosomal affinity for tRNAs in the P-site and for ternary complex in the A-site, and inhibit binding of the translocase, eEF2. These defects in turn affect specific aspects of translational fidelity. These findings suggest that uS19 plays a critical role as a conduit of information exchange between the large and small ribosomal subunits directly through the B1a, and indirectly through the B1b/c bridges.

Ab initio RNA folding

RNA molecules are essential cellular machines performing a wide variety of functions for which a specific three-dimensional structure is required. Over the last several years, the experimental determination of RNA structures through x-ray crystallography and NMR seems to have reached a plateau in the number of structures resolved each year, but as more and more RNA sequences are being discovered, the need for structure prediction tools to complement experimental data is strong. Theoretical approaches to RNA folding have been developed since the late nineties, when the first algorithms for secondary structure prediction appeared. Over the last 10 years a number of prediction methods for 3D structures have been developed, first based on bioinformatics and data-mining, and more recently based on a coarse-grained physical representation of the systems. In this review we are going to present the challenges of RNA structure prediction and the main ideas behind bioinformatic approaches and physics-based approaches. We will focus on the description of the more recent physics-based phenomenological models and on how they are built to include the specificity of the interactions of RNA bases, whose role is critical in folding. Through examples from different models, we will point out the strengths of physics-based approaches, which are able not only to predict equilibrium structures, but also to investigate dynamical and thermodynamical behavior, and the open challenges to include more key interactions ruling RNA folding.

From static to dynamic: the need for structural ensembles and a predictive model of RNA folding and function

To understand RNA, it is necessary to move beyond a descriptive categorization towards quantitative predictions of its molecular conformations and functional behavior. An incisive approach to understanding the function and folding of biological RNA systems involves characterizing small, simple components that are largely responsible for the behavior of complex systems including helix-junction-helix elements and tertiary motifs. State-of-the-art methods have permitted unprecedented insight into the conformational ensembles of these elements revealing, for example, that conformations of helix-junction-helix elements are confined to a small region of the ensemble, that this region is highly dependent on the junction's topology, and that the correct alignment of tertiary motifs may be a rare conformation on the overall folding landscape. Further characterization of RNA components and continued development of experimental and computational methods with the goal of quantitatively predicting RNA folding and functional behavior will be critical to understanding biological RNA systems.

Molecular mechanics of 30S subunit head rotation

During ribosomal translocation, a process central to the elongation phase of protein synthesis, movement of mRNA and tRNAs requires large-scale rotation of the head domain of the small (30S) subunit of the ribosome. It has generally been accepted that the head rotates by pivoting around the neck helix (h28) of 16S rRNA, its sole covalent connection to the body domain. Surprisingly, we observe that the calculated axis of rotation does not coincide with the neck. Instead, comparative structure analysis across 55 ribosome structures shows that 30S head movement results from flexing at two hinge points lying within conserved elements of 16S rRNA. Hinge 1, although located within the neck, moves by straightening of the kinked helix h28 at the point of contact with the mRNA. Hinge 2 lies within a three-way helix junction that extends to the body through a second, noncovalent connection; its movement results from flexing between helices h34 and h35 in a plane orthogonal to the movement of hinge 1. Concerted movement at these two hinges accounts for the observed magnitudes of head rotation. Our findings also explain the mode of action of spectinomycin, an antibiotic that blocks translocation by binding to hinge 2.

Modeling of proteins and their assemblies with the Integrative Modeling Platform

To understand the workings of the living cell, we need to characterize protein assemblies that constitute the cell (for example, the ribosome, 26S proteasome, and the nuclear pore complex). A reliable high-resolution structural characterization of these assemblies is frequently beyond the reach of current experimental methods, such as X-ray crystallography, NMR spectroscopy, electron microscopy, footprinting, chemical cross-linking, FRET spectroscopy, small angle X-ray scattering, and proteomics. However, the information garnered from different methods can be combined and used to build models of the assembly structures that are consistent with all of the available datasets, and therefore more accurate, precise, and complete. Here, we describe a protocol for this integration, whereby the information is converted to a set of spatial restraints and a variety of optimization procedures can be used to generate models that satisfy the restraints as well as possible. These generated models can then potentially inform about the precision and accuracy of structure determination, the accuracy of the input datasets, and further data generation. We also demonstrate the Integrative Modeling Platform (IMP) software, which provides the necessary computational framework to implement this protocol, and several applications for specific use cases.

Enhanced sampling and overfitting analyses in structural refinement of nucleic acids into electron microscopy maps

Flexible fitting computational algorithms are often useful to interpret low-resolution maps of many macromolecular complexes generated by electron microscopy (EM) imaging. One such atomistic simulation technique is molecular dynamics flexible fitting (MDFF), which has been widely applied to generate structural models of large ribonucleoprotein assemblies such as the ribosome. We have previously shown that MDFF simulations of globular proteins are sensitive to the resolution of the target EM map and the strength of restraints used to preserve the secondary structure elements during fitting (Vashisth, H.; et al. Structure 2012, 20, 1453-1462). In this work, we aim to systematically examine the quality of structural models of various nucleic acids obtained via MDFF by varying the map resolution and the strength of structural restraints. We also demonstrate how an enhanced conformational sampling technique for proteins, temperature-accelerated molecular dynamics (TAMD), can be combined with MDFF for the structural refinement of nucleic acids in EM maps. Finally, we also demonstrate application of TAMD-assisted MDFF (TAMDFF) on a RNA/protein complex and suggest that TAMDFF is a viable strategy for enhanced conformational fitting in target maps of ribonucleoprotein complexes.

Structural insights into the mechanism of translational inhibition by the fungicide sordarin

The translational machinery has been found to be the target for a number of antibiotics. One such antibiotic sordarin selectively inhibits fungal translation by impairing the function of elongation factor 2 (eEF2) while being ineffective to higher eukaryotes. Surprisingly, sordarin is not even equally effective in impairing translation for all fungal species. The binding cavity of sordarin on eEF2 has been localized by X-ray crystallographic study and its unique specificity towards sordarin has been attributed to the species specific substitutions within a stretch of amino acids (sordarin specificity region, SSR) at the entrance of the cavity. In this study, we have analyzed the sordarin-binding cavity of eEF2 from different species both in isolated and ribosome-bound forms in order to decipher the mechanism of sordarin binding selectivity. Our results reveal that the molecular architecture as well as the microenvironment of the sordarin-binding cavity changes significantly from one species to another depending on the species specific substitutions within the cavity. Moreover, eEF2 binding to ribosome aggravates the effects of these substitutions. Thus, this study, while shedding light on the molecular mechanism underpinning the selective inhibitory effects of sordarin, will also be a helpful guide for future studies aiming at developing novel antifungal drugs with broader spectrum of activity.

An Introduction to the Structure and Function of the Ribosome

E. coli continues to serve as a key model for the structure and function of the ribosome, structures of ribosome from other organisms and domains of life have also greatly contributed to our knowledge of protein synthesis. Many structural models of the ribosome in a number of steps of the protein synthesis cycle have been solved by cryo-electron microscopy (cryo-EM) and x-ray crystallography. This chapter introduces the structure and dynamics of the ribosome based on these structures and ends with a brief discussion of the many questions that the structures leave unanswered. Protein synthesis is a multistep process, and the structural features of the ribosome along with the large number of cofactors reflect the complexity of translation. Numerous protein factors in addition to the ribosome contribute to translation in bacteria during the steps of initiation, elongation, termination, and recycling. These protein factors make intimate contacts to key regions of the ribosome, and this aspect is discussed in the chapter in light of our present understanding of the structure and function of the ribosome. The intact ribosome contains three binding sites for substrate tRNAs that are termed as the aminoacyl-tRNA site (A site), peptidyl-tRNA site (P site), and exit-tRNA site (E site). These three binding sites span the interface between the 30S and 50S subunits. The central activity of the ribosome is catalysis of peptide bond formation. The region of the ribosome responsible for catalyzing the reaction is called the peptidyl transferase center (PTC).

Using enhanced sampling and structural restraints to refine atomic structures into low-resolution electron microscopy maps

For a variety of problems in structural biology, low-resolution maps generated by electron microscopy imaging are often interpreted with the help of various flexible-fitting computational algorithms. In this work, we systematically analyze the quality of final models of various proteins obtained via molecular dynamics flexible fitting (MDFF) by varying the map-resolution, strength of structural restraints, and the steering forces. We find that MDFF can be extended to understand conformational changes in lower-resolution maps if larger structural restraints and lower steering forces are used to prevent overfitting. We further show that the capabilities of MDFF can be extended by combining it with an enhanced conformational sampling method, temperature-accelerated molecular dynamics (TAMD). Specifically, either TAMD can be used to generate better starting configurations for MDFF fitting or TAMD-assisted MDFF (TAMDFF) can be performed to accelerate conformational search in atomistic simulations.

Identical RNA-protein interactions in vivo and in vitro and a scheme of folding the newly synthesized proteins by ribosomes

A distinct three-dimensional shape of rRNA inside the ribosome is required for the peptidyl transfer activity of its peptidyltransferase center (PTC). In contrast, even the in vitro transcribed PTC RNA interacts with unfolded protein(s) at about five sites to let them attain their native states. We found that the same set of conserved nucleotides in the PTC interact identically with nascent and chemically unfolded proteins in vivo and in vitro, respectively. The time course of this interaction, difficult to follow in vivo, was observed in vitro. It suggested nucleation of folding of cytosolic globular proteins vectorially from hydrophilic N to hydrophobic C termini, consistent with our discovery of a regular arrangement of cumulative hydrophobic indices of the peptide segments of cytosolic proteins from N to C termini. Based on this observation, we propose a model here for the nucleation of folding of the nascent protein chain by the PTC.

Bridging the gap between ribosome structure and biochemistry by mechanistic computations

The wealth of structural and biochemical data now available for protein synthesis on the ribosome presents major new challenges for computational biochemistry. Apart from technical difficulties in modeling ribosome systems, the complexity of the overall translation cycle with a multitude of different kinetic steps presents a formidable problem for computational efforts where we have only seen the beginning. However, a range of methodologies including molecular dynamics simulations, free energy calculations, molecular docking and quantum chemical approaches have already been put to work with promising results. In particular, the combined efforts of structural biology, biochemistry, kinetics and computational modeling can lead towards a quantitative structure-based description of translation.

Site-specific platinum(II) cross-linking in a ribozyme active site

The function of RNA depends on its ability to adopt complex and dynamic structures, and the incorporation of site-specific cross-linking probes is a powerful method for providing distance constraints that are valuable in RNA structural biology. Here we describe a new RNA-RNA cross-linking strategy based on Pt(II) targeting of specific phosphorothioate substitutions. In this strategy cis-diammine Pt(II) complexes are kinetically recruited and anchored to a phosphorothioate substitution embedded within a structured RNA. Substitution of the remaining exchangeable Pt(II) ligand with a nucleophile supplied by a nearby RNA nucleobase results in metal-mediated cross-links that are stable during isolation. This type of cross-linking strategy was explored within the catalytic core of the Hammerhead ribozyme (HHRz). When a phosphorothioate substitution is installed at the scissile bond normally cleaved by the HHRz, Pt(II) cross-linking takes place to nucleotides G8 and G10 in the ribozyme active site. Both of these positions are predicted to be within ~8 Å of a phosphorothioate-bound Pt(II) metal center. Cross-linking depends on Mg(2+) ion concentration, reaching yields as high as 30%, with rates that indicate cation competition within the RNA three-helix junction. Cross-linking efficiency depends on accurate formation of the HHRz tertiary structure, and cross-links are not observed for RNA helices. Combined, these results show promise for using kinetically inert Pt(II) complexes as new site-specific cross-linking tools for exploring RNA structure and dynamics.

A-site residues move independently from P-site residues in all-atom molecular dynamics simulations of the 70S bacterial ribosome

The ribosome is a large macromolecular machine, and correlated motion between residues is necessary for coordinating function across multiple protein and RNA chains. We ran two all-atom, explicit solvent molecular dynamics simulations of the bacterial ribosome and calculated correlated motion between residue pairs by using mutual information. Because of the short timescales of our simulation (ns), we expect that dynamics are largely local fluctuations around the crystal structure. We hypothesize that residues that show coupled dynamics are functionally related, even on longer timescales. We validate our model by showing that crystallographic B-factors correlate well with the entropy calculated as part of our mutual information calculations. We reveal that A-site residues move relatively independently from P-site residues, effectively insulating A-site functions from P-site functions during translation.

Key intermolecular interactions in the E. coli 70S ribosome revealed by coarse-grained analysis

The ribosome is a very large complex that consists of many RNA and protein molecules and plays a central role in protein biosynthesis in all organisms. Extensive interactions between different molecules are critical to ribosomal functional dynamics. In this work, intermolecular interactions in the Escherichia coli 70S ribosome are investigated by coarse-grained (CG) analysis. CG models are defined to preserve dynamic domains in RNAs and proteins and to capture functional motions in the ribosome, and then the CG sites are connected by harmonic springs, and spring constants are obtained by matching the computed fluctuations to those of an all-atom molecular dynamics (MD) simulation. Those spring constants indicate how strong the interactions are between the ribosomal components, and they are in good agreement with various experimental data. Nearly all the bridges between the small and large ribosomal subunits are indicated by CG interactions with large spring constants. The head of the small subunit is very mobile because it has minimal CG interactions with the rest of the subunit; however, a large number of small subunit proteins bind to maintain the internal structure of the head. The results show a clear connection between the intermolecular interactions and the structural and functional properties of the ribosome because of the reduced complexity in domain-based CG models. The present approach also provides a useful strategy to map interactions between molecules within large biomolecular complexes since it is not straightforward to investigate these by either atomistic MD simulations or residue-based elastic network models.

Allosteric vs. spontaneous exit-site (E-site) tRNA dissociation early in protein synthesis

During protein synthesis, deacylated transfer RNAs leave the ribosome via an exit (E) site after mRNA translocation. How the ribosome regulates tRNA dissociation and whether functional linkages between the aminoacyl (A) and E sites modulate the dynamics of protein synthesis have long been debated. Using single molecule fluorescence resonance energy transfer experiments, we find that, during early cycles of protein elongation, tRNAs are often held in the E site until being allosterically released when the next aminoacyl tRNA binds to the A site. This process is regulated by the length and sequence of the nascent peptide and by the conformational state, detected by tRNA proximity, prior to translocation. In later cycles, E-site tRNA dissociates spontaneously. Our results suggest that the distribution of pretranslocation tRNA states and posttranslocation pathways are correlated within each elongation cycle via communication between distant subdomains in the ribosome, but that this correlation between elongation cycle intermediates does not persist into succeeding cycles.

Understanding ribosome assembly: the structure of in vivo assembled immature 30S subunits revealed by cryo-electron microscopy

Four decades after early in vitro assembly studies demonstrated that ribosome assembly is a controlled process, our understanding of ribosome assembly is still incomplete. Just as structure determination has been so important to understanding ribosome function, so too will it be critical to sorting out the assembly process. Here, we used a viable deletion in the yjeQ gene, a recognized ribosome assembly factor, to isolate and structurally characterize immature 30S subunits assembled in vivo. These small ribosome subunits contained unprocessed 17S rRNA and lacked some late ribosomal proteins. Cryo-electron microscopy reconstructions revealed that the presence of precursor sequences in the rRNA induces a severe distortion in the 3′ minor domain of the subunit involved in the decoding of mRNA and interaction with the large ribosome subunit. These findings suggest that rRNA processing events induce key local conformational changes directing the structure toward the mature assembly. We concluded that rRNA processing, folding, and the entry of tertiary r-proteins are interdependent events in the late stages of 30S subunit assembly. In addition, we demonstrate how studies of emerging assembly factors in ribosome biogenesis can help to elucidate the path of subunit assembly in vivo.

Predicting and modeling RNA architecture

A general approach for modeling the architecture of large and structured RNA molecules is described. The method exploits the modularity and the hierarchical folding of RNA architecture that is viewed as the assembly of preformed double-stranded helices defined by Watson-Crick base pairs and RNA modules maintained by non-Watson-Crick base pairs. Despite the extensive molecular neutrality observed in RNA structures, specificity in RNA folding is achieved through global constraints like lengths of helices, coaxiality of helical stacks, and structures adopted at the junctions of helices. The Assemble integrated suite of computer tools allows for sequence and structure analysis as well as interactive modeling by homology or ab initio assembly with possibilities for fitting within electronic density maps. The local key role of non-Watson-Crick pairs guides RNA architecture formation and offers metrics for assessing the accuracy of three-dimensional models in a more useful way than usual root mean square deviation (RMSD) values.

Distribution of dwell times of a ribosome: effects of infidelity, kinetic proofreading and ribosome crowding

Ribosome is a molecular machine that polymerizes a protein where the sequence of the amino acid residues, the monomers of the protein, is dictated by the sequence of codons (triplets of nucleotides) on a messenger RNA (mRNA) that serves as the template. The ribosome is a molecular motor that utilizes the template mRNA strand also as the track. Thus, in each step the ribosome moves forward by one codon and, simultaneously, elongates the protein by one amino acid. We present a theoretical model that captures most of the main steps in the mechanochemical cycle of a ribosome. The stochastic movement of the ribosome consists of an alternating sequence of pause and translocation; the sum of the durations of a pause and the following translocation is the time of dwell of the ribosome at the corresponding codon. We derive the analytical expression for the distribution of the dwell times of a ribosome in our model. Wherever experimental data are available, our theoretical predictions are consistent with those results. We suggest appropriate experiments to test the new predictions of our model, particularly the effects of the quality control mechanism of the ribosome and that of their crowding on the mRNA track.

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.

Modeling of proteins and their assemblies with the integrative modeling platform

To understand the workings of the living cell, we need to characterize protein assemblies that constitute the cell (for example, the ribosome, 26S proteasome, and the nuclear pore complex). A reliable high-resolution structural characterization of these assemblies is frequently beyond the reach of current experimental methods, such as X-ray crystallography, NMR spectroscopy, electron microscopy, footprinting, chemical cross-linking, FRET spectroscopy, small-angle X-ray scattering, and proteomics. However, the information garnered from different methods can be combined and used to build computational models of the assembly structures that are consistent with all of the available datasets. Here, we describe a protocol for this integration, whereby the information is converted to a set of spatial restraints and a variety of optimization procedures can be used to generate models that satisfy the restraints as much as possible. These generated models can then potentially inform about the precision and accuracy of structure determination, the accuracy of the input datasets, and further data generation. We also demonstrate the Integrative Modeling Platform (IMP) software, which provides the necessary computational framework to implement this protocol, and several applications for specific-use cases.

Quantum chemical studies of nucleic acids: can we construct a bridge to the RNA structural biology and bioinformatics communities?

In this feature article, we provide a side-by-side introduction for two research fields: quantum chemical calculations of molecular interaction in nucleic acids and RNA structural bioinformatics. Our main aim is to demonstrate that these research areas, while largely separated in contemporary literature, have substantial potential to complement each other that could significantly contribute to our understanding of the exciting world of nucleic acids. We identify research questions amenable to the combined application of modern ab initio methods and bioinformatics analysis of experimental structures while also assessing the limitations of these approaches. The ultimate aim is to attain valuable physicochemical insights regarding the nature of the fundamental molecular interactions and how they shape RNA structures, dynamics, function, and evolution.

An approach for de novo structure determination of dynamic molecular assemblies by electron cryomicroscopy

Single-particle electron cryomicroscopy is a powerful method for three-dimensional (3D) structure determination of macromolecular assemblies. Here we address the challenge of determining a 3D structure in the absence of reference models. The 3D structures are determined by alignment and weighted averaging of densities obtained by native cryo random conical tilt (RCT) reconstructions including consideration of missing data. Our weighted averaging scheme (wRCT) offers advantages for potentially heterogeneous 3D densities of low signal-to-noise ratios. Sets of aligned RCT structures can also be analyzed by multivariate statistical analysis (MSA) to provide insights into snapshots of the assemblies. The approach is used to compute 3D structures of the Escherichia coli 70S ribosome and the human U4/U6.U5 tri-snRNP under vitrified unstained cryo conditions, and to visualize by 3D MSA the L7/L12 stalk of the 70S ribosome and states of tri-snRNP. The approach thus combines de novo 3D structure determination with an analysis of compositional and conformational heterogeneity.

Structure and dynamics of a processive Brownian motor: the translating ribosome

There is mounting evidence indicating that protein synthesis is driven and regulated by mechanisms that direct stochastic, large-scale conformational fluctuations of the translational apparatus. This mechanistic paradigm implies that a free-energy landscape governs the conformational states that are accessible to and sampled by the translating ribosome. This scenario presents interdependent opportunities and challenges for structural and dynamic studies of protein synthesis. Indeed, the synergism between cryogenic electron microscopic and X-ray crystallographic structural studies, on the one hand, and single-molecule fluorescence resonance energy transfer (smFRET) dynamic studies, on the other, is emerging as a powerful means for investigating the complex free-energy landscape of the translating ribosome and uncovering the mechanisms that direct the stochastic conformational fluctuations of the translational machinery. In this review, we highlight the principal insights obtained from cryogenic electron microscopic, X-ray crystallographic, and smFRET studies of the elongation stage of protein synthesis and outline the emerging themes, questions, and challenges that lie ahead in mechanistic studies of translation.

Integrative structure modeling of macromolecular assemblies from proteomics data

Proteomics techniques have been used to generate comprehensive lists of protein interactions in a number of species. However, relatively little is known about how these interactions result in functional multiprotein complexes. This gap can be bridged by combining data from proteomics experiments with data from established structure determination techniques. Correspondingly, integrative computational methods are being developed to provide descriptions of protein complexes at varying levels of accuracy and resolution, ranging from complex compositions to detailed atomic structures.

Molecular dynamics simulations suggest that RNA three-way junctions can act as flexible RNA structural elements in the ribosome

We present extensive explicit solvent molecular dynamics analysis of three RNA three-way junctions (3WJs) from the large ribosomal subunit: the 3WJ formed by Helices 90–92 (H90–H92) of 23S rRNA; the 3WJ formed by H42–H44 organizing the GTPase associated center (GAC) of 23S rRNA; and the 3WJ of 5S rRNA. H92 near the peptidyl transferase center binds the 3′-CCA end of amino-acylated tRNA. The GAC binds protein factors and stimulates GTP hydrolysis driving protein synthesis. The 5S rRNA binds the central protuberance and A-site finger (ASF) involved in bridges with the 30S subunit. The simulations reveal that all three 3WJs possess significant anisotropic hinge-like flexibility between their stacked stems and dynamics within the compact regions of their adjacent stems. The A-site 3WJ dynamics may facilitate accommodation of tRNA, while the 5S 3WJ flexibility appears to be essential for coordinated movements of ASF and 5S rRNA. The GAC 3WJ may support large-scale dynamics of the L7/L12-stalk region. The simulations reveal that H42–H44 rRNA segments are not fully relaxed and in the X-ray structures they are bent towards the large subunit. The bending may be related to L10 binding and is distributed between the 3WJ and the H42–H97 contact.

Evolutionary tabu search strategies for the simultaneous registration of multiple atomic structures in cryo-EM reconstructions

A structural characterization of multi-component cellular assemblies is essential to explain the mechanisms governing biological function. Macromolecular architectures may be revealed by integrating information collected from various biophysical sources - for instance, by interpreting low-resolution electron cryomicroscopy reconstructions in relation to the crystal structures of the constituent fragments. A simultaneous registration of multiple components is beneficial when building atomic models as it introduces additional spatial constraints to facilitate the native placement inside the map. The high-dimensional nature of such a search problem prevents the exhaustive exploration of all possible solutions. Here we introduce a novel method based on genetic algorithms, for the efficient exploration of the multi-body registration search space. The classic scheme of a genetic algorithm was enhanced with new genetic operations, tabu search and parallel computing strategies and validated on a benchmark of synthetic and experimental cryo-EM datasets. Even at a low level of detail, for example 35-40 A, the technique successfully registered multiple component biomolecules, measuring accuracies within one order of magnitude of the nominal resolutions of the maps. The algorithm was implemented using the Sculptor molecular modeling framework, which also provides a user-friendly graphical interface and enables an instantaneous, visual exploration of intermediate solutions.

Building and refining protein models within cryo-electron microscopy density maps based on homology modeling and multiscale structure refinement

Automatic modeling methods using cryoelectron microscopy (cryoEM) density maps as constraints are promising approaches to building atomic models of individual proteins or protein domains. However, their application to large macromolecular assemblies has not been possible largely due to computational limitations inherent to such unsupervised methods. Here we describe a new method, EM-IMO (electron microscopy-iterative modular optimization), for building, modifying and refining local structures of protein models using cryoEM maps as a constraint. As a supervised refinement method, EM-IMO allows users to specify parameters derived from inspections so as to guide, and as a consequence, significantly speed up the refinement. An EM-IMO-based refinement protocol is first benchmarked on a data set of 50 homology models using simulated density maps. A multiscale refinement strategy that combines EM-IMO-based and molecular dynamics-based refinement is then applied to build backbone models for the seven conformers of the five capsid proteins in our near-atomic-resolution cryoEM map of the grass carp reovirus virion, a member of the Aquareovirus genus of the Reoviridae family. The refined models allow us to reconstruct a backbone model of the entire grass carp reovirus capsid and provide valuable functional insights that are described in the accompanying publication [Cheng, L., Zhu, J., Hui, W. H., Zhang, X., Honig, B., Fang, Q. & Zhou, Z. H. (2010). Backbone model of an aquareovirus virion by cryo-electron microscopy and bioinformatics. J. Mol. Biol. (this issue). doi:10.1016/j.jmb.2009.12.027.]. Our study demonstrates that the integrated use of homology modeling and a multiscale refinement protocol that combines supervised and automated structure refinement offers a practical strategy for building atomic models based on medium- to high-resolution cryoEM density maps.

Probing ribosome-nascent chain complexes produced in vivo by NMR spectroscopy

The means by which a polypeptide chain acquires its unique 3-D structure is a fundamental question in biology. During its synthesis on the ribosome, a nascent chain (NC) emerges vectorially and will begin to fold in a cotranslational fashion. The complex environment of the cell, coupled with the gradual emergence of the ribosome-tethered NC during its synthesis, imposes conformational restraints on its folding landscape that differ from those placed on an isolated protein when stimulated to fold following denaturation in solution. To begin to examine cotranslational folding as it would occur within a cell, we produce highly selective, isotopically labeled NCs bound to isotopically silent ribosomes in vivo. We then apply NMR spectroscopy to study, at a residue specific level, the conformation of NCs consisting of different fractional lengths of the polypeptide chain corresponding to a given protein. This combined approach provides a powerful means of generating a series of snapshots of the folding of the NC as it emerges from the ribosome. Application of this strategy to the NMR analysis of the progressive synthesis of an Ig-like domain reveals the existence of a partially folded ribosome-bound species that is likely to represent an intermediate species populated during the cotranslational folding process.

Single particle electron microscopy

Electron microscopy (EM) in combination with image analysis is a powerful technique to study protein structures at low, medium, and high resolution. Since electron micrographs of biological objects are very noisy, improvement of the signal-to-noise ratio by image processing is an integral part of EM, and this is performed by averaging large numbers of individual projections. Averaging procedures can be divided into crystallographic and non-crystallographic methods. The crystallographic averaging method, based on two-dimensional (2D) crystals of (membrane) proteins, yielded in solving atomic protein structures in the last century. More recently, single particle analysis could be extended to solve atomic structures as well. It is a suitable method for large proteins, viruses, and proteins that are difficult to crystallize. Because it is also a fast method to reveal the low-to-medium resolution structures, the impact of its application is growing rapidly. Technical aspects, results, and possibilities are presented.