Ribosome-induced changes in elongation factor Tu conformation control GTP hydrolysis

Supercomputing in the biological sciences: Toward Zettascale and Yottascale simulations

Molecular simulations of biological systems tend to be significantly more compute-intensive than those in materials science and astrophysics, due to important contributions of long-range electrostatic forces and large numbers of time steps (>1E9) required. Simulations of biomolecular complexes of microseconds to milliseconds are considered state-of-the-art today. However, these time scales are miniscule in comparison to physiological time scales relevant to molecular machine activity, drug action, and elongation cycles for protein synthesis, RNA synthesis, and DNA synthesis (seconds to days). While an exascale supercomputer has simulated an entire virus for nanoseconds, this supercomputer would need to be 10 billion times faster to simulate that virus for 3 hours of physiological time, demonstrating the insatiable need for computing power. With growing interest in computational drug design from the pharmaceutical sector, the biological sciences are positioned to be an industry driver in computing.

The central role of transfer RNAs in mistranslation

Transfer RNAs (tRNA) are essential small non-coding RNAs that enable the translation of genomic information into proteins in all life forms. The principal function of tRNAs is to bring amino acid building blocks to the ribosomes for protein synthesis. In the ribosome, tRNAs interact with messenger RNA (mRNA) to mediate the incorporation of amino acids into a growing polypeptide chain following the rules of the genetic code. Accurate interpretation of the genetic code requires tRNAs to carry amino acids matching their anticodon identity and decode the correct codon on mRNAs. Errors in these steps cause the translation of codons with the wrong amino acids (mistranslation), compromising the accurate flow of information from DNA to proteins. Accumulation of mutant proteins due to mistranslation jeopardizes proteostasis and cellular viability. However, the concept of mistranslation is evolving, with increasing evidence indicating that mistranslation can be used as a mechanism for survival and acclimatization to environmental conditions. In this review, we discuss the central role of tRNAs in modulating translational fidelity through their dynamic and complex interplay with translation factors. We summarize recent discoveries of mistranslating tRNAs and describe the underlying molecular mechanisms and the specific conditions and environments that enable and promote mistranslation.

β-Amino Acids Reduce Ternary Complex Stability and Alter the Translation Elongation Mechanism

Templated synthesis of proteins containing non-natural amino acids (nnAAs) promises to expand the chemical space available to biological therapeutics and materials, but existing technologies are still limiting. Addressing these limitations requires a deeper understanding of the mechanism of protein synthesis and how it is perturbed by nnAAs. Here we examine the impact of nnAAs on the formation and ribosome utilization of the central elongation substrate: the ternary complex of native, aminoacylated tRNA, thermally unstable elongation factor, and GTP. By performing ensemble and single-molecule fluorescence resonance energy transfer measurements, we reveal that both the (R)- and (S)-β2 isomers of phenylalanine (Phe) disrupt ternary complex formation to levels below in vitro detection limits, while (R)- and (S)-β3-Phe reduce ternary complex stability by 1 order of magnitude. Consistent with these findings, (R)- and (S)-β2-Phe-charged tRNAs were not utilized by the ribosome, while (R)- and (S)-β3-Phe stereoisomers were utilized inefficiently. (R)-β3-Phe but not (S)-β3-Phe also exhibited order of magnitude defects in the rate of translocation after mRNA decoding. We conclude from these findings that non-natural amino acids can negatively impact the translation mechanism on multiple fronts and that the bottlenecks for improvement must include the consideration of the efficiency and stability of ternary complex formation.

β-amino acids reduce ternary complex stability and alter the translation elongation mechanism

Templated synthesis of proteins containing non-natural amino acids (nnAAs) promises to vastly expand the chemical space available to biological therapeutics and materials. Existing technologies limit the identity and number of nnAAs than can be incorporated into a given protein. Addressing these bottlenecks requires deeper understanding of the mechanism of messenger RNA (mRNA) templated protein synthesis and how this mechanism is perturbed by nnAAs. Here we examine the impact of both monomer backbone and side chain on formation and ribosome-utilization of the central protein synthesis substate: the ternary complex of native, aminoacylated transfer RNA (aa-tRNA), thermally unstable elongation factor (EF-Tu), and GTP. By performing ensemble and single-molecule fluorescence resonance energy transfer (FRET) measurements, we reveal the dramatic effect of monomer backbone on ternary complex formation and protein synthesis. Both the (R) and (S)-β2 isomers of Phe disrupt ternary complex formation to levels below in vitro detection limits, while (R)- and (S)-β3-Phe reduce ternary complex stability by approximately one order of magnitude. Consistent with these findings, (R)- and (S)-β2-Phe-charged tRNAs were not utilized by the ribosome, while (R)- and (S)-β3-Phe stereoisomers were utilized inefficiently. The reduced affinities of both species for EF-Tu ostensibly bypassed the proofreading stage of mRNA decoding. (R)-β3-Phe but not (S)-β3-Phe also exhibited order of magnitude defects in the rate of substrate translocation after mRNA decoding, in line with defects in peptide bond formation that have been observed for D-α-Phe. We conclude from these findings that non-natural amino acids can negatively impact the translation mechanism on multiple fronts and that the bottlenecks for improvement must include consideration of the efficiency and stability of ternary complex formation.

Geometric alignment of aminoacyl-tRNA relative to catalytic centers of the ribosome underpins accurate mRNA decoding

Accurate protein synthesis is determined by the two-subunit ribosome's capacity to selectively incorporate cognate aminoacyl-tRNA for each mRNA codon. The molecular basis of tRNA selection accuracy, and how fidelity can be affected by antibiotics, remains incompletely understood. Using molecular simulations, we find that cognate and near-cognate tRNAs delivered to the ribosome by Elongation Factor Tu (EF-Tu) can follow divergent pathways of motion into the ribosome during both initial selection and proofreading. Consequently, cognate aa-tRNAs follow pathways aligned with the catalytic GTPase and peptidyltransferase centers of the large subunit, while near-cognate aa-tRNAs follow pathways that are misaligned. These findings suggest that differences in mRNA codon-tRNA anticodon interactions within the small subunit decoding center, where codon-anticodon interactions occur, are geometrically amplified over distance, as a result of this site's physical separation from the large ribosomal subunit catalytic centers. These insights posit that the physical size of both tRNA and ribosome are key determinants of the tRNA selection fidelity mechanism.

Visualizing translation dynamics at atomic detail inside a bacterial cell

Translation is the fundamental process of protein synthesis and is catalysed by the ribosome in all living cells1. Here we use advances in cryo-electron tomography and sub-tomogram analysis2,3 to visualize the structural dynamics of translation inside the bacterium Mycoplasma pneumoniae. To interpret the functional states in detail, we first obtain a high-resolution in-cell average map of all translating ribosomes and build an atomic model for the M. pneumoniae ribosome that reveals distinct extensions of ribosomal proteins. Classification then resolves 13 ribosome states that differ in their conformation and composition. These recapitulate major states that were previously resolved in vitro, and reflect intermediates during active translation. On the basis of these states, we animate translation elongation inside native cells and show how antibiotics reshape the cellular translation landscapes. During translation elongation, ribosomes often assemble in defined three-dimensional arrangements to form polysomes4. By mapping the intracellular organization of translating ribosomes, we show that their association into polysomes involves a local coordination mechanism that is mediated by the ribosomal protein L9. We propose that an extended conformation of L9 within polysomes mitigates collisions to facilitate translation fidelity. Our work thus demonstrates the feasibility of visualizing molecular processes at atomic detail inside cells.

Proteome Analysis of the Soybean Nodule Phosphorus Response Mechanism and Characterization of Stress-Induced Ribosome Structural and Protein Expression Changes

In agroecosystems, a plant-usable form of nitrogen is mainly generated by legume-based biological nitrogen fixation, a process that requires phosphorus (P) as an essential nutrient. To investigate the physiological mechanism whereby phosphorus influences soybean nodule nitrogen fixation, soybean root nodules were exposed to four phosphate levels: 1 mg/L (P stress), 11 mg/L (P stress), 31 mg/L (Normal P), and 61 mg/L (High P) then proteome analysis of nodules was conducted to identify phosphorus-associated proteome changes. We found that phosphorus stress-induced ribosomal protein structural changes were associated with altered key root nodule protein synthesis profiles. Importantly, up-regulated expression of peroxidase was observed as an important phosphorus stress-induced nitrogen fixation-associated adaptation that supported two nodule-associated activities: scavenging of reactive oxygen species (ROS) and cell wall growth. In addition, phosphorus transporter (PT) and purple acid phosphatase (PAPs) were up-regulated that regulated phosphorus transport and utilization to maintain phosphorus balance and nitrogen fixation function in phosphorus-stressed root nodules.

Genome-Wide Investigation of Pasteurella multocida Identifies the Stringent Response as a Negative Regulator of Hyaluronic Acid Capsule Production

The bacterial pathogen P. multocida can cause serious disease in production animals, including fowl cholera in poultry, hemorrhagic septicemia in cattle and buffalo, atrophic rhinitis in pigs, and respiratory diseases in a range of livestock. P. multocida produces a capsule that is essential for systemic disease, but the complete mechanisms underlying synthesis and regulation of capsule production are not fully elucidated. A whole-genome analysis using TraDIS was undertaken to identify genes essential for growth in rich media and to obtain a comprehensive characterization of capsule production.

Regulation of translation by site-specific ribosomal RNA methylation

Ribosomes are complex ribozymes that interpret genetic information by translating messenger RNA (mRNA) into proteins. Natural variation in ribosome composition has been documented in several organisms and can arise from several different sources. A key question is whether specific control over ribosome heterogeneity represents a mechanism by which translation can be regulated. We used RiboMeth-seq to demonstrate that differential 2'-O-methylation of ribosomal RNA (rRNA) represents a considerable source of ribosome heterogeneity in human cells, and that modification levels at distinct sites can change dynamically in response to upstream signaling pathways, such as MYC oncogene expression. Ablation of one prominent methylation resulted in altered translation of select mRNAs and corresponding changes in cellular phenotypes. Thus, differential rRNA 2'-O-methylation can give rise to ribosomes with specialized function. This suggests a broader mechanism where the specific regulation of rRNA modification patterns fine tunes translation.

Snapshots of native pre-50S ribosomes reveal a biogenesis factor network and evolutionary specialization

Ribosome biogenesis is a fundamental multi-step cellular process that culminates in the formation of ribosomal subunits, whose production and modification are regulated by numerous biogenesis factors. In this study, we analyze physiologic prokaryotic ribosome biogenesis by isolating bona fide pre-50S subunits from an Escherichia coli strain with the biogenesis factor ObgE, affinity tagged at its native gene locus. Our integrative structural approach reveals a network of interacting biogenesis factors consisting of YjgA, RluD, RsfS, and ObgE on the immature pre-50S subunit. In addition, our study provides mechanistic insight into how the GTPase ObgE, in concert with other biogenesis factors, facilitates the maturation of the 50S functional core and reveals both conserved and divergent evolutionary features of ribosome biogenesis between prokaryotes and eukaryotes.

Scalable molecular dynamics on CPU and GPU architectures with NAMD

NAMDis a molecular dynamics program designed for high-performance simulations of very large biological objects on CPU- and GPU-based architectures. NAMD offers scalable performance on petascale parallel supercomputers consisting of hundreds of thousands of cores, as well as on inexpensive commodity clusters commonly found in academic environments. It is written in C++ and leans on Charm++ parallel objects for optimal performance on low-latency architectures. NAMD is a versatile, multipurpose code that gathers state-of-the-art algorithms to carry out simulations in apt thermodynamic ensembles, using the widely popular CHARMM, AMBER, OPLS, and GROMOS biomolecular force fields. Here, we review the main features of NAMD that allow both equilibrium and enhanced-sampling molecular dynamics simulations with numerical efficiency. We describe the underlying concepts utilized by NAMD and their implementation, most notably for handling long-range electrostatics; controlling the temperature, pressure, and pH; applying external potentials on tailored grids; leveraging massively parallel resources in multiple-copy simulations; and hybrid quantum-mechanical/molecular-mechanical descriptions. We detail the variety of options offered by NAMD for enhanced-sampling simulations aimed at determining free-energy differences of either alchemical or geometrical transformations and outline their applicability to specific problems. Last, we discuss the roadmap for the development of NAMD and our current efforts toward achieving optimal performance on GPU-based architectures, for pushing back the limitations that have prevented biologically realistic billion-atom objects to be fruitfully simulated, and for making large-scale simulations less expensive and easier to set up, run, and analyze. NAMD is distributed free of charge with its source code at www.ks.uiuc.edu.

Cryo-EM of elongating ribosome with EF-Tu•GTP elucidates tRNA proofreading

Ribosomes accurately decode mRNA by proofreading each aminoacyl-tRNA delivered by elongation factor EF-Tu¹. Understanding the molecular mechanism of proofreading requires visualizing GTP-catalyzed elongation, which has remained a challenge^2–4. Here, time-resolved cryo-EM revealed 33 states following aminoacyl-tRNA delivery by EF-Tu•GTP. Instead of locking cognate tRNA upon initial recognition, the ribosomal decoding center (DC) dynamically monitors codon-anticodon interactions before and after GTP hydrolysis. GTP hydrolysis allows EF-Tu’s GTPase domain to extend away, releasing EF-Tu from tRNA. Then, the 30S subunit locks cognate tRNA in the DC, and rotates, enabling the tRNA to bypass 50S protrusions during accommodation into the peptidyl transferase center. By contrast, the DC fails to lock near-cognate tRNA, allowing dissociation of near-cognate tRNA during both initial selection (before GTP hydrolysis) and proofreading (after GTP hydrolysis). These findings reveal structural similarity between initial selection^5,6 and the previously unseen proofreading, which together govern efficient rejection of incorrect tRNA.

Elongation Factor Tu Switch I Element is a Gate for Aminoacyl-tRNA Selection

Selection of correct aminoacyl (aa)-tRNA at the ribosomal A site is fundamental to maintaining translational fidelity. Aa-tRNA selection is a multistep process facilitated by the guanosine triphosphatase elongation factor (EF)-Tu. EF-Tu delivers aa-tRNA to the ribosomal A site and participates in tRNA selection. The structural mechanism of how EF-Tu is involved in proofreading remains to be fully resolved. Here, we provide evidence that switch I of EF-Tu facilitates EF-Tu's involvement during aa-tRNA selection. Using structure-based and explicit solvent molecular dynamics simulations based on recent cryo-electron microscopy reconstructions, we studied the conformational change of EF-Tu from the guanosine triphosphate to guanine diphosphate conformation during aa-tRNA accommodation. Switch I of EF-Tu rapidly converts from an α-helix into a β-hairpin and moves to interact with the acceptor stem of the aa-tRNA. In doing so, switch I gates the movement of the aa-tRNA during accommodation through steric interactions with the acceptor stem. Pharmacological inhibition of the aa-tRNA accommodation pathway prevents the proper positioning of switch I with the aa-tRNA acceptor stem, suggesting that the observed interactions are specific for cognate aa-tRNA substrates, and thus capable of contributing to the fidelity mechanism.

Elongation factor-Tu can repetitively engage aminoacyl-tRNA within the ribosome during the proofreading stage of tRNA selection

Elongation factor Tu (EF-Tu) facilitates rapid and accurate selection of aminoacyl-tRNA (aa-tRNA) by the bacterial ribosome during protein synthesis. We show that EF-Tu dissociates from the ribosome as aa-tRNA navigates the accommodation corridor en route to peptide bond formation. We find that EF-Tu’s release from the ribosome during aa-tRNA selection can be reversible. We also demonstrate that new ternary complex formation, accompanied by futile cycles of GTP hydrolysis, can occur on aa-tRNA bound within the ribosome. These findings inform on the decoding mechanism, the contributions of EF-Tu to the fidelity of translation, and the potential consequences of reduced rates of peptide bond formation on cellular physiology.

The substrate for ribosomes actively engaged in protein synthesis is a ternary complex of elongation factor Tu (EF-Tu), aminoacyl-tRNA (aa-tRNA), and GTP. EF-Tu plays a critical role in mRNA decoding by increasing the rate and fidelity of aa-tRNA selection at each mRNA codon. Here, using three-color single-molecule fluorescence resonance energy transfer imaging and molecular dynamics simulations, we examine the timing and role of conformational events that mediate the release of aa-tRNA from EF-Tu and EF-Tu from the ribosome after GTP hydrolysis. Our investigations reveal that conformational changes in EF-Tu coordinate the rate-limiting passage of aa-tRNA through the accommodation corridor en route to the peptidyl transferase center of the large ribosomal subunit. Experiments using distinct inhibitors of the accommodation process further show that aa-tRNA must at least partially transit the accommodation corridor for EF-Tu⋅GDP to release. aa-tRNAs failing to undergo peptide bond formation at the end of accommodation corridor passage after EF-Tu release can be reengaged by EF-Tu⋅GTP from solution, coupled to GTP hydrolysis. These observations suggest that additional rounds of ternary complex formation can occur on the ribosome during proofreading, particularly when peptide bond formation is slow, which may serve to increase both the rate and fidelity of protein synthesis at the expense of GTP hydrolysis.

Insights into the base-pairing preferences of 8-oxoguanosine on the ribosome

Of the four bases, guanine is the most susceptible to oxidation, which results in the formation of 8-oxoguanine (8-oxoG). In protein-free DNA, 8-oxodG adopts the syn conformation more frequently than the anti one. In the syn conformation, 8-oxodG base pairs with dA. The equilibrium between the anti and syn conformations of the adduct are known to be altered by the enzyme recognizing 8-oxodG. We previously showed that 8-oxoG in mRNA severely disrupts tRNA selection, but the underlying mechanism for these effects was not addressed. Here, we use miscoding antibiotics and ribosome mutants to probe how 8-oxoG interacts with the tRNA anticodon in the decoding center. Addition of antibiotics and introduction of error-inducing mutations partially suppressed the effects of 8-oxoG. Under these conditions, rates and/or endpoints of peptide-bond formation for the cognate (8-oxoG•C) and near-cognate (8-oxoG•A) aminoacyl-tRNAs increased. In contrast, the antibiotics had little effect on other mismatches, suggesting that the lesion restricts the nucleotide from forming other interactions. Our findings suggest that 8-oxoG predominantly adopts the syn conformation in the A site. However, its ability to base pair with adenosine in this conformation is not sufficient to promote the necessary structural changes for tRNA selection to proceed.

Synthetic 5' UTRs Can Either Up- or Downregulate Expression upon RNA-Binding Protein Binding

The construction of complex gene-regulatory networks requires both inhibitory and upregulatory modules. However, the vast majority of RNA-based regulatory "parts" are inhibitory. Using a synthetic biology approach combined with SHAPE-seq, we explored the regulatory effect of RNA-binding protein (RBP)-RNA interactions in bacterial 5' UTRs. By positioning a library of RNA hairpins upstream of a reporter gene and co-expressing them with the matching RBP, we observed a set of regulatory responses, including translational stimulation, translational repression, and cooperative behavior. Our combined approach revealed three distinct states in vivo: in the absence of RBPs, the RNA molecules can be found in either a molten state that is amenable to translation or a structured phase that inhibits translation. In the presence of RBPs, the RNA molecules are in a semi-structured phase with partial translational capacity. Our work provides new insight into RBP-based regulation and a blueprint for designing complete gene-regulatory circuits at the post-transcriptional level.

Cryo-EM shows stages of initial codon selection on the ribosome by aa-tRNA in ternary complex with GTP and the GTPase-deficient EF-TuH84A

The GTPase EF-Tu in ternary complex with GTP and aminoacyl-tRNA (aa-tRNA) promotes rapid and accurate delivery of cognate aa-tRNAs to the ribosomal A site. Here we used cryo-EM to study the molecular origins of the accuracy of ribosome-aided recognition of a cognate ternary complex and the accuracy-amplifying role of the monitoring bases A1492, A1493 and G530 of the 16S rRNA. We used the GTPase-deficient EF-Tu variant H84A with native GTP, rather than non-cleavable GTP analogues, to trap a near-cognate ternary complex in high-resolution ribosomal complexes of varying codon-recognition accuracy. We found that ribosome complexes trapped by GTPase-deficicent ternary complex due to the presence of EF-Tu^H84A or non-cleavable GTP analogues have very similar structures. We further discuss speed and accuracy of initial aa-tRNA selection in terms of conformational changes of aa-tRNA and stepwise activation of the monitoring bases at the decoding center of the ribosome.

An Overview of Molecular Modeling for Drug Discovery with Specific Illustrative Examples of Applications

In this paper we review the current status of high-performance computing applications in the general area of drug discovery. We provide an introduction to the methodologies applied at atomic and molecular scales, followed by three specific examples of implementation of these tools. The first example describes in silico modeling of the adsorption of small molecules to organic and inorganic surfaces, which may be applied to drug delivery issues. The second example involves DNA translocation through nanopores with major significance to DNA sequencing efforts. The final example offers an overview of computer-aided drug design, with some illustrative examples of its usefulness.

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.

Biogenic regions of cyanobacterial thylakoids form contact sites with the plasma membrane

Little is known about how the photosynthetic machinery is arranged in time and space during the biogenesis of thylakoid membranes. Using in situ cryo-electron tomography to image the three-dimensional architecture of the cyanobacterium Synechocystis, we observed that the tips of multiple thylakoids merge to form a substructure called the 'convergence membrane'. This high-curvature membrane comes into close contact with the plasma membrane at discrete sites. We generated subtomogram averages of 70S ribosomes and array-forming phycobilisomes, then mapped these structures onto the native membrane architecture as markers for protein synthesis and photosynthesis, respectively. This molecular localization identified two distinct biogenic regions in the thylakoid network: thylakoids facing the cytosolic interior of the cell that were associated with both marker complexes, and convergence membranes that were decorated by ribosomes but not phycobilisomes. We propose that the convergence membranes perform a specialized biogenic function, coupling the synthesis of thylakoid proteins with the integration of cofactors from the plasma membrane and the periplasmic space.

Joachim Frank's Binding with the Ribosome

With recent technological advancements, single-particle cryogenic electron microscopy (cryo-EM) is now the technique of choice to study structure and function of biological macromolecules at near-atomic resolution. Many single-particle EM reconstruction methods necessary for these advances were pioneered by Joachim Frank, and were optimized using the ribosome as a benchmark specimen. In doing so, he made several landmark contributions to the understanding of the structure and function of ribosomes. These include the first 3D visualization of ribosome-bound transfer RNAs, the first experimentally derived structures of the primary complexes formed during the bacterial translation elongation cycle, and the critical ribosomal conformational transitions required for translation. Over the years, his laboratory studied many important functional complexes of the ribosome from both eubacterial and eukaryotic systems, including ribosomes from pathogenic organisms. This article presents a brief account of the contributions made by Joachim Frank to the ribosome field.

Hydroxylation of protein constituents of the human translation system: structural aspects and functional assignments

During the current decade, data on the post-translational hydroxylation of specific amino acid residues of some ribosomal proteins and translation factors in both eukaryotes and eubacteria have accumulated. The reaction is catalyzed by dedicated oxygenases (so-called ribosomal oxygenases), whose action is impaired under hypoxia conditions. The modification occurs at amino acid residues directly involved in the formation of the main functional sites of ribosomes and factors. This review summarizes currently available data on the specific hydroxylation of protein constituents of eukaryotic and eubacterial translation systems with a special emphasis on the human system, as well as on the links between hypoxia impacts on the operation of ribosomal oxygenases, the functioning of the translational apparatus and human health problems.

New Insights into Ribosome Structure and Function

In the past 4 years, because of the advent of new cameras, many ribosome structures have been solved by cryoelectron microscopy (cryo-EM) at high, often near-atomic resolution, bringing new mechanistic insights into the processes of translation initiation, peptide elongation, termination, and recycling. Thus, cryo-EM has joined X-ray crystallography as a powerful technique in structural studies of translation. The significance of this new development is that structures of ribosomes in complex with their functional binding partners can now be determined to high resolution in multiple states as they perform their work. The aim of this article is to provide an overview of these new studies and assess the contributions they have made toward an understanding of translation and translational control.

Visualization of biological macromolecules at near-atomic resolution: cryo-electron microscopy comes of age

Structural biology is going through a revolution as a result of transformational advances in the field of cryo-electron microscopy (cryo-EM) driven by the development of direct electron detectors and ultrastable electron microscopes. High-resolution cryo-EM images of isolated biomolecules (single particles) suspended in a thin layer of vitrified buffer are subjected to powerful image-processing algorithms, enabling near-atomic resolution structures to be determined in unprecedented numbers. Prior to these advances, electron crystallography of two-dimensional crystals and helical assemblies of proteins had established the feasibility of atomic resolution structure determination using cryo-EM. Atomic resolution single-particle analysis, without the need for crystals, now promises to resolve problems in structural biology that were intractable just a few years ago.

Exploring Dynamics and Structure of Biomolecules, Cryoprotectants, and Water Using Molecular Dynamics Simulations: Implications for Biostabilization and Biopreservation

Successful stabilization and preservation of biological materials often utilize low temperatures and dehydration to arrest molecular motion. Cryoprotectants are routinely employed to help the biological entities survive the physicochemical and mechanical stresses induced by cold or dryness. Molecular interactions between biomolecules, cryoprotectants, and water fundamentally determine the outcomes of preservation. The optimization of assays using the empirical approach is often limited in structural and temporal resolution, whereas classical molecular dynamics simulations can provide a cost-effective glimpse into the atomic-level structure and interaction of individual molecules that dictate macroscopic behavior. Computational research on biomolecules, cryoprotectants, and water has provided invaluable insights into the development of new cryoprotectants and the optimization of preservation methods. We describe the rapidly evolving state of the art of molecular simulations of these complex systems, summarize the molecular-scale protective and stabilizing mechanisms, and discuss the challenges that motivate continued innovation in this field.

Constructing atomic structural models into cryo-EM densities using molecular dynamics - Pros and cons

Accurate structure determination from electron density maps at 3-5 Å resolution necessitates a balance between extensive global and local sampling of atomistic models, yet with the stereochemical correctness of backbone and sidechain geometries. Molecular Dynamics Flexible Fitting (MDFF), particularly through a resolution-exchange scheme, ReMDFF, provides a robust way of achieving this balance for hybrid structure determination. Employing two high-resolution density maps, namely that of β-galactosidase at 3.2 Å and TRPV1 at 3.4 Å, we showcase the quality of ReMDFF-generated models, comparing them against ones submitted by independent research groups for the 2015-2016 Cryo-EM Model Challenge. This comparison offers a clear evaluation of ReMDFF's strengths and shortcomings, and those of data-guided real-space refinements in general. ReMDFF results scored highly on the various metric for judging the quality-of-fit and quality-of-model. However, some systematic discrepancies are also noted employing a Molprobity analysis, that are reproducible across multiple competition entries. A space of key refinement parameters is explored within ReMDFF to observe their impact within the final model. Choice of force field parameters and initial model seem to have the most significant impact on ReMDFF model-quality. To this end, very recently developed CHARMM36m force field parameters provide now more refined ReMDFF models than the ones originally submitted to the Cryo-EM challenge. Finally, a set of good-practices is prescribed for the community to benefit from the MDFF developments.

Peptide Mimics of the Ribosomal P Stalk Inhibit the Activity of Ricin A Chain by Preventing Ribosome Binding

Ricin A chain (RTA) depurinates the sarcin/ricin loop (SRL) by interacting with the C-termini of the ribosomal P stalk. The ribosome interaction site and the active site are located on opposite faces of RTA. The interaction with P proteins allows RTA to depurinate the SRL on the ribosome at physiological pH with an extremely high activity by orienting the active site towards the SRL. Therefore, if an inhibitor disrupts RTA⁻ribosome interaction by binding to the ribosome binding site of RTA, it should inhibit the depurination activity. To test this model, we synthesized peptides mimicking the last 3 to 11 amino acids of P proteins and examined their interaction with wild-type RTA and ribosome binding mutants by Biacore. We measured the inhibitory activity of these peptides on RTA-mediated depurination of yeast and rat liver ribosomes. We found that the peptides interacted with the ribosome binding site of RTA and inhibited depurination activity by disrupting RTA⁻ribosome interactions. The shortest peptide that could interact with RTA and inhibit its activity was four amino acids in length. RTA activity was inhibited by disrupting its interaction with the P stalk without targeting the active site, establishing the ribosome binding site as a new target for inhibitor discovery.

Single-Particle Reconstruction of Biological Molecules-Story in a Sample (Nobel Lecture)

Pictures tell a thousand words: The development of single-particle cryo-electron microscopy set the stage for high-resolution structure determination of biological molecules. In his Nobel lecture, J. Frank describes the ground-breaking discoveries that have enabled the development of cryo-EM. The method has taken biochemistry into a new era.

CryoEM-based hybrid modeling approaches for structure determination

Recent advances in cryo-electron microscopy (cryoEM) have dramatically improved the resolutions at which vitrified biological specimens can be studied, revealing new structural and mechanistic insights over a broad range of spatial scales. Bolstered by these advances, much effort has been directed toward the development of hybrid modeling methodologies for the construction and refinement of high-fidelity atomistic models from cryoEM data. In this brief review, we will survey the key elements of cryoEM-based hybrid modeling, providing an overview of available computational tools and strategies as well as several recent applications.

RNA Structural Dynamics As Captured by Molecular Simulations: A Comprehensive Overview

With both catalytic and genetic functions, ribonucleic acid (RNA) is perhaps the most pluripotent chemical species in molecular biology, and its functions are intimately linked to its structure and dynamics. Computer simulations, and in particular atomistic molecular dynamics (MD), allow structural dynamics of biomolecular systems to be investigated with unprecedented temporal and spatial resolution. We here provide a comprehensive overview of the fast-developing field of MD simulations of RNA molecules. We begin with an in-depth, evaluatory coverage of the most fundamental methodological challenges that set the basis for the future development of the field, in particular, the current developments and inherent physical limitations of the atomistic force fields and the recent advances in a broad spectrum of enhanced sampling methods. We also survey the closely related field of coarse-grained modeling of RNA systems. After dealing with the methodological aspects, we provide an exhaustive overview of the available RNA simulation literature, ranging from studies of the smallest RNA oligonucleotides to investigations of the entire ribosome. Our review encompasses tetranucleotides, tetraloops, a number of small RNA motifs, A-helix RNA, kissing-loop complexes, the TAR RNA element, the decoding center and other important regions of the ribosome, as well as assorted others systems. Extended sections are devoted to RNA-ion interactions, ribozymes, riboswitches, and protein/RNA complexes. Our overview is written for as broad of an audience as possible, aiming to provide a much-needed interdisciplinary bridge between computation and experiment, together with a perspective on the future of the field.

The translation elongation cycle-capturing multiple states by cryo-electron microscopy

During the work cycle of elongation, the ribosome, a molecular machine of vast complexity, exists in a large number of states distinguished by constellation of its subunits, its subunit domains and binding partners. Single-particle cryogenic electron microscopy (cryo-EM), developed over the past 40 years, is uniquely suited to determine the structure of molecular machines in their native states. With the emergence, 10 years ago, of unsupervised clustering techniques in the analysis of single-particle data, it has been possible to determine multiple structures from a sample containing ribosomes equilibrating in different thermally accessible states. In addition, recent advances in detector technology have made it possible to reach near-atomic resolution for some of these states. With these capabilities, single-particle cryo-EM has been at the forefront of exploring ribosome dynamics during its functional cycle, along with single-molecule fluorescence resonance energy transfer and molecular dynamics computations, offering insights into molecular architecture uniquely honed by evolution to capitalize on thermal energy in the ambient environment.This article is part of the themed issue 'Perspectives on the ribosome'.

The potential of cryo-electron microscopy for structure-based drug design

Structure-based drug design plays a central role in therapeutic development. Until recently, protein crystallography and NMR have dominated experimental approaches to obtain structural information of biological molecules. However, in recent years rapid technical developments in single particle cryo-electron microscopy (cryo-EM) have enabled the determination to near-atomic resolution of macromolecules ranging from large multi-subunit molecular machines to proteins as small as 64 kDa. These advances have revolutionized structural biology by hugely expanding both the range of macromolecules whose structures can be determined, and by providing a description of macromolecular dynamics. Cryo-EM is now poised to similarly transform the discipline of structure-based drug discovery. This article reviews the potential of cryo-EM for drug discovery with reference to protein ligand complex structures determined using this technique.

Tracking fluctuation hotspots on the yeast ribosome through the elongation cycle

Chemical modification was used to quantitatively determine the flexibility of nearly the entire rRNA component of the yeast ribosome through 8 discrete stages of translational elongation, revealing novel observations at the gross and fine-scales. These include (i) the bulk transfer of energy through the intersubunit bridges from the large to the small subunit after peptidyltransfer, (ii) differences in the interaction of the sarcin ricin loop with the two elongation factors and (iii) networked information exchange pathways that may functionally facilitate intra- and intersubunit coordination, including the 5.8S rRNA. These analyses reveal hot spots of fluctuations that set the stage for large-scale conformational changes essential for translocation and enable the first molecular dynamics simulation of an 80S complex. Comprehensive datasets of rRNA base flexibilities provide a unique resource to the structural biology community that can be computationally mined to complement ongoing research toward the goal of understanding the dynamic ribosome.

Elfamycins: inhibitors of elongation factor-Tu

Elfamycins are a relatively understudied group of antibiotics that target the essential process of translation through impairment of EF-Tu function. For the most part, the utility of these compounds has been as laboratory tools for the study of EF-Tu and the ribosome, as their poor pharmacokinetic profile and solubility has prevented implementation as therapeutic agents. However, due to the slowing of the antibiotic pipeline and the rapid emergence of resistance to approved antibiotics, this group is being reconsidered. Some researchers are using screens for novel naturally produced variants, while others are making directed, systematic chemical improvements on publically disclosed compounds. As an example of the latter approach, a GE2270 A derivative, LFF571, has completed phase 2 clinical trials, thus demonstrating the potential for elfamycins to become more prominent antibiotics in the future.

Structure of the 40S-ABCE1 post-splitting complex in ribosome recycling and translation initiation

The essential ATP-binding cassette protein ABCE1 splits 80S ribosomes into 60S and 40S subunits after canonical termination or quality-control-based mRNA surveillance processes. However, the underlying splitting mechanism remains enigmatic. Here, we present a cryo-EM structure of the yeast 40S-ABCE1 post-splitting complex at 3.9-Å resolution. Compared to the pre-splitting state, we observe repositioning of ABCE1's iron-sulfur cluster domain, which rotates 150° into a binding pocket on the 40S subunit. This repositioning explains a newly observed anti-association activity of ABCE1. Notably, the movement implies a collision with A-site factors, thus explaining the splitting mechanism. Disruption of key interactions in the post-splitting complex impairs cellular homeostasis. Additionally, the structure of a native post-splitting complex reveals ABCE1 to be part of the 43S initiation complex, suggesting a coordination of termination, recycling, and initiation.

YcaO-Dependent Posttranslational Amide Activation: Biosynthesis, Structure, and Function

With advances in sequencing technology, uncharacterized proteins and domains of unknown function (DUFs) are rapidly accumulating in sequence databases and offer an opportunity to discover new protein chemistry and reaction mechanisms. The focus of this review, the formerly enigmatic YcaO superfamily (DUF181), has been found to catalyze a unique phosphorylation of a ribosomal peptide backbone amide upon attack by different nucleophiles. Established nucleophiles are the side chains of Cys, Ser, and Thr which gives rise to azoline/azole biosynthesis in ribosomally synthesized and posttranslationally modified peptide (RiPP) natural products. However, much remains unknown about the potential for YcaO proteins to collaborate with other nucleophiles. Recent work suggests potential in forming thioamides, macroamidines, and possibly additional post-translational modifications. This review covers all knowledge through mid-2016 regarding the biosynthetic gene clusters (BGCs), natural products, functions, mechanisms, and applications of YcaO proteins and outlines likely future research directions for this protein superfamily.

The mechanism of translation

Translation of the genetic code on the ribosome into protein is a process of extraordinary complexity, and understanding its mechanism has remained one of the major challenges even though x-ray structures have been available since 2000. In the past two decades, single-particle cryo-electron microscopy has contributed a major share of information on structure, binding modes, and conformational changes of the ribosome during its work cycle, but the contributions of this technique in the translation field have recently skyrocketed after the introduction of a new recording medium capable of detecting individual electrons. As many examples in the recent literature over the past three years show, the impact of this development on the advancement of knowledge in this field has been transformative and promises to be lasting.

CryoEM Structure Refinement by Integrating NMR Chemical Shifts with Molecular Dynamics Simulations

Single particle cryoEM has emerged as a powerful method for structure determination of proteins and complexes, complementing X-ray crystallography and NMR spectroscopy. Yet, for many systems, the resolution of cryoEM density map has been limited to 4-6 Å, which only allows for resolving bulky amino acids side chains, thus hindering accurate model building from the density map. On the other hand, experimental chemical shifts (CS) from solution and solid state MAS NMR spectra provide atomic level data for each amino acid within a molecule or a complex; however, structure determination of large complexes and assemblies based on NMR data alone remains challenging. Here, we present a novel integrated strategy to combine the highly complementary experimental data from cryoEM and NMR computationally by molecular dynamics simulations to derive an atomistic model, which is not attainable by either approach alone. We use the HIV-1 capsid protein (CA) C-terminal domain as well as the large capsid assembly to demonstrate the feasibility of this approach, termed NMR CS-biased cryoEM structure refinement.

Advances in the field of single-particle cryo-electron microscopy over the last decade

In single-particle cryo-electron microscopy (cryo-EM), molecules suspended in a thin aqueous layer are rapidly frozen and imaged at cryogenic temperature in the transmission electron microscope. From the random projection views, a three-dimensional image is reconstructed, enabling the structure of the molecule to be obtained. In this article I discuss technological progress over the past decade, which has, in my own field of study, culminated in the determination of ribosome structure at 2.5-Å resolution. I also discuss likely future improvements in methodology.

Partially Assembled Nucleosome Structures at Atomic Detail

The evidence is now overwhelming that partially assembled nucleosome states (PANS) are as important as the canonical nucleosome structure for the understanding of how accessibility to genomic DNA is regulated in cells. We use a combination of molecular dynamics simulation and atomic force microscopy to deliver, in atomic detail, structural models of three key PANS: the hexasome (H2A·H2B)·(H3·H4)₂, the tetrasome (H3·H4)₂, and the disome (H3·H4). Despite fluctuations of the conformation of the free DNA in these structures, regions of protected DNA in close contact with the histone core remain stable, thus establishing the basis for the understanding of the role of PANS in DNA accessibility regulation. On average, the length of protected DNA in each structure is roughly 18 basepairs per histone protein. Atomistically detailed PANS are used to explain experimental observations; specifically, we discuss interpretation of atomic force microscopy, Förster resonance energy transfer, and small-angle x-ray scattering data obtained under conditions when PANS are expected to exist. Further, we suggest an alternative interpretation of a recent genome-wide study of DNA protection in active chromatin of fruit fly, leading to a conclusion that the three PANS are present in actively transcribing regions in a substantial amount. The presence of PANS may not only be a consequence, but also a prerequisite for fast transcription in vivo.

Unravelling biological macromolecules with cryo-electron microscopy

Knowledge of the three-dimensional structures of proteins and other biological macromolecules often aids understanding of how they perform complicated tasks in the cell. Because many such tasks involve the cleavage or formation of chemical bonds, structural characterization at the atomic level is most useful. Developments in the electron microscopy of frozen hydrated samples (cryo-electron microscopy) are providing unprecedented opportunities for the structural characterization of biological macromolecules. This is resulting in a wave of information about processes in the cell that were impossible to characterize with existing techniques in structural biology.

The C-terminal Helix of Pseudomonas aeruginosa Elongation Factor Ts Tunes EF-Tu Dynamics to Modulate Nucleotide Exchange

Little is known about the conservation of critical kinetic parameters and the mechanistic strategies of elongation factor (EF) Ts-catalyzed nucleotide exchange in EF-Tu in bacteria and particularly in clinically relevant pathogens. EF-Tu from the clinically relevant pathogen Pseudomonas aeruginosa shares over 84% sequence identity with the corresponding elongation factor from Escherichia coli Interestingly, the functionally closely linked EF-Ts only shares 55% sequence identity. To identify any differences in the nucleotide binding properties, as well as in the EF-Ts-mediated nucleotide exchange reaction, we performed a comparative rapid kinetics and mutagenesis analysis of the nucleotide exchange mechanism for both the E. coli and P. aeruginosa systems, identifying helix 13 of EF-Ts as a previously unnoticed regulatory element in the nucleotide exchange mechanism with species-specific elements. Our findings support the base side-first entry of the nucleotide into the binding pocket of the EF-Tu·EF-Ts binary complex, followed by displacement of helix 13 and rapid binding of the phosphate side of the nucleotide, ultimately leading to the release of EF-Ts.

Molecular dynamics-based refinement and validation for sub-5 Å cryo-electron microscopy maps

Two structure determination methods, based on the molecular dynamics flexible fitting (MDFF) paradigm, are presented that resolve sub-5 Å cryo-electron microscopy (EM) maps with either single structures or ensembles of such structures. The methods, denoted cascade MDFF and resolution exchange MDFF, sequentially re-refine a search model against a series of maps of progressively higher resolutions, which ends with the original experimental resolution. Application of sequential re-refinement enables MDFF to achieve a radius of convergence of ~25 Å demonstrated with the accurate modeling of β-galactosidase and TRPV1 proteins at 3.2 Å and 3.4 Å resolution, respectively. The MDFF refinements uniquely offer map-model validation and B-factor determination criteria based on the inherent dynamics of the macromolecules studied, captured by means of local root mean square fluctuations. The MDFF tools described are available to researchers through an easy-to-use and cost-effective cloud computing resource on Amazon Web Services.

DOI:http://dx.doi.org/10.7554/eLife.16105.001

To understand the roles that proteins and other large molecules play inside cells, it is important to determine their structures. One of the techniques that researchers can use to do this is called cryo-electron microscopy (cryo-EM), which rapidly freezes molecules to fix them in position before imaging them in fine detail.

The cryo-EM images are like maps that show the approximate position of atoms. These images must then be processed in order to build a three-dimensional model of the protein that shows how its atoms are arranged relative to each other. One computational approach called Molecular Dynamics Flexible Fitting (MDFF) works by flexibly fitting possible atomic structures into cryo-EM maps. Although this approach works well with relatively undetailed (or ‘low resolution’) cryo-EM images, it struggles to handle the high-resolution cryo-EM maps now being generated.

Singharoy, Teo, McGreevy et al. have now developed two MDFF methods – called cascade MDFF and resolution exchange MDFF – that help to resolve atomic models of biological molecules from cryo-EM images. Each method can refine poorly guessed models into ones that are consistent with the high-resolution experimental images. The refinement is achieved by interpreting a range of images that starts with a ‘fuzzy’ image. The contrast of the image is then progressively improved until an image is produced that has a resolution that is good enough to almost distinguish individual atoms.

The method works because each cryo-EM image shows not just one, but a collection of atomic structures that the molecule can take on, with the fuzzier parts of the image representing the more flexible parts of the molecule. By taking into account this flexibility, the large-scale features of the protein structure can be determined first from the fuzzier images, and increasing the contrast of the images allows smaller-scale refinements to be made to the structure.

The MDFF tools have been designed to be easy to use and are available to researchers at low cost through cloud computing platforms. They can now be used to unravel the structure of many different proteins and protein complexes including those involved in photosynthesis, respiration and protein synthesis.

DOI:http://dx.doi.org/10.7554/eLife.16105.002

Computational Methodologies for Real-Space Structural Refinement of Large Macromolecular Complexes

The rise of the computer as a powerful tool for model building and refinement has revolutionized the field of structure determination for large biomolecular systems. Despite the wide availability of robust experimental methods capable of resolving structural details across a range of spatiotemporal resolutions, computational hybrid methods have the unique ability to integrate the diverse data from multimodal techniques such as X-ray crystallography and electron microscopy into consistent, fully atomistic structures. Here, commonly employed strategies for computational real-space structural refinement are reviewed, and their specific applications are illustrated for several large macromolecular complexes: ribosome, virus capsids, chemosensory array, and photosynthetic chromatophore. The increasingly important role of computational methods in large-scale structural refinement, along with current and future challenges, is discussed.

Advances in the molecular dynamics flexible fitting method for cryo-EM modeling

Molecular Dynamics Flexible Fitting (MDFF) is an established technique for fitting all-atom structures of molecules into corresponding cryo-electron microscopy (cryo-EM) densities. The practical application of MDFF is simple but requires a user to be aware of and take measures against a variety of possible challenges presented by each individual case. Some of these challenges arise from the complexity of a molecular structure or the limited quality of available structural models and densities to be interpreted, while others stem from the intricacies of MDFF itself. The current article serves as an overview of the strategies that have been developed since MDFF's inception to overcome common challenges and successfully perform MDFF simulations.

Characterizing RNA structures in vitro and in vivo with selective 2'-hydroxyl acylation analyzed by primer extension sequencing (SHAPE-Seq)

RNA molecules adopt a wide variety of structures that perform many cellular functions, including, among others, catalysis, small molecule sensing, and cellular defense. Our ability to characterize, predict, and design RNA structures are key factors for understanding and controlling the biological roles of RNAs. Fortunately, there has been rapid progress in this area, especially with respect to experimental methods that can characterize RNA structures in a high throughput fashion using chemical probing and next-generation sequencing. Here, we describe one such method, selective 2'-hydroxyl acylation analyzed by primer extension sequencing (SHAPE-Seq), which measures nucleotide resolution flexibility information for RNAs in vitro and in vivo. We outline the process of designing and performing a SHAPE-Seq experiment and describe methods for using experimental SHAPE-Seq data to restrain computational folding algorithms to generate more accurate predictions of RNA secondary structure. We also provide a number of examples of SHAPE-Seq reactivity spectra obtained in vitro and in vivo and discuss important considerations for performing SHAPE-Seq experiments, both in terms of collecting and analyzing data. Finally, we discuss improvements and extensions of these experimental and computational techniques that promise to deepen our knowledge of RNA folding and function.