Flexible fitting of atomic structures into electron microscopy maps using molecular dynamics

Visualizing the transfer-messenger RNA as the ribosome resumes translation

Bacterial ribosomes stalled by truncated mRNAs are rescued by transfer-messenger RNA (tmRNA), a dual-function molecule that contains a tRNA-like domain (TLD) and an internal open reading frame (ORF). Occupying the empty A site with its TLD, the tmRNA enters the ribosome with the help of elongation factor Tu and a protein factor called small protein B (SmpB), and switches the translation to its own ORF. In this study, using cryo-electron microscopy, we obtained the first structure of an in vivo-formed complex containing ribosome and the tmRNA at the point where the TLD is accommodated into the ribosomal P site. We show that tmRNA maintains a stable 'arc and fork' structure on the ribosome when its TLD moves to the ribosomal P site and translation resumes on its ORF. Based on the density map, we built an atomic model, which suggests that SmpB interacts with the five nucleotides immediately upstream of the resume codon, thereby determining the correct selection of the reading frame on the ORF of tmRNA.

Localization of eukaryote-specific ribosomal proteins in a 5.5-Å cryo-EM map of the 80S eukaryotic ribosome

Protein synthesis in all living organisms occurs on ribonucleoprotein particles, called ribosomes. Despite the universality of this process, eukaryotic ribosomes are significantly larger in size than their bacterial counterparts due in part to the presence of 80 r proteins rather than 54 in bacteria. Using cryoelectron microscopy reconstructions of a translating plant (Triticum aestivum) 80S ribosome at 5.5-Å resolution, together with a 6.1-Å map of a translating Saccharomyces cerevisiae 80S ribosome, we have localized and modeled 74/80 (92.5%) of the ribosomal proteins, encompassing 12 archaeal/eukaryote-specific small subunit proteins as well as the complete complement of the ribosomal proteins of the eukaryotic large subunit. Near-complete atomic models of the 80S ribosome provide insights into the structure, function, and evolution of the eukaryotic translational apparatus.

Formation of salt bridges mediates internal dimerization of myosin VI medial tail domain

The unconventional motor protein, myosin VI, is known to dimerize upon cargo binding to its C-terminal end. It has been shown that one of its tail domains, called the medial tail domain, is a dimerization region. The domain contains an unusual pattern of alternating charged residues and a few hydrophobic residues. To reveal the unknown dimerization mechanism of the medial tail domain, we employed molecular dynamics and single-molecule experimental techniques. Both techniques suggest that the formation of electrostatic-based interhelical salt bridges between oppositely charged residues is a key dimerization factor. For the dimerization to occur, the two identical helices within the dimer do not bind in a symmetric fashion, but rather with an offset of about one helical repeat. Calculations of the dimer-dissociation energy find the contribution of hydrophobic residues to the dimerization process to be minor; they also find that the asymmetric homodimer state is energetically favorable over a state of separate helices.

Structural basis for translational stalling by human cytomegalovirus and fungal arginine attenuator peptide

Specific regulatory nascent chains establish direct interactions with the ribosomal tunnel, leading to translational stalling. Despite a wealth of biochemical data, structural insight into the mechanism of translational stalling in eukaryotes is still lacking. Here we use cryo-electron microscopy to visualize eukaryotic ribosomes stalled during the translation of two diverse regulatory peptides: the fungal arginine attenuator peptide (AAP) and the human cytomegalovirus (hCMV) gp48 upstream open reading frame 2 (uORF2). The C terminus of the AAP appears to be compacted adjacent to the peptidyl transferase center (PTC). Both nascent chains interact with ribosomal proteins L4 and L17 at tunnel constriction in a distinct fashion. Significant changes at the PTC were observed: the eukaryotic-specific loop of ribosomal protein L10e establishes direct contact with the CCA end of the peptidyl-tRNA (P-tRNA), which may be critical for silencing of the PTC during translational stalling. Our findings provide direct structural insight into two distinct eukaryotic stalling processes.

Mechanical properties of a complete microtubule revealed through molecular dynamics simulation

Microtubules (MTs) are the largest type of cellular filament, essential in processes ranging from mitosis and meiosis to flagellar motility. Many of the processes depend critically on the mechanical properties of the MT, but the elastic moduli, notably the Young's modulus, are not directly revealed in experiment, which instead measures either flexural rigidity or response to radial deformation. Molecular dynamics (MD) is a method that allows the mechanical properties of single biomolecules to be investigated through computation. Typically, MD requires an atomic resolution structure of the molecule, which is unavailable for many systems, including MTs. By combining structural information from cryo-electron microscopy and electron crystallography, we have constructed an all-atom model of a complete MT and used MD to determine its mechanical properties. The simulations revealed nonlinear axial stress-strain behavior featuring a pronounced softening under extension, a possible plastic deformation transition under radial compression, and a distinct asymmetry in response to the two senses of twist. This work demonstrates the possibility of combining different levels of structural information to produce all-atom models suitable for quantitative MD simulations, which extends the range of systems amenable to the MD method and should enable exciting advances in our microscopic knowledge of biology.

Cryo-EM structure and rRNA model of a translating eukaryotic 80S ribosome at 5.5-A resolution

Protein biosynthesis, the translation of the genetic code into polypeptides, occurs on ribonucleoprotein particles called ribosomes. Although X-ray structures of bacterial ribosomes are available, high-resolution structures of eukaryotic 80S ribosomes are lacking. Using cryoelectron microscopy and single-particle reconstruction, we have determined the structure of a translating plant (Triticum aestivum) 80S ribosome at 5.5-Å resolution. This map, together with a 6.1-Å map of a Saccharomyces cerevisiae 80S ribosome, has enabled us to model ∼98% of the rRNA. Accurate assignment of the rRNA expansion segments (ES) and variable regions has revealed unique ES-ES and r-protein-ES interactions, providing insight into the structure and evolution of the eukaryotic ribosome.

tmRNA-SmpB: a journey to the centre of the bacterial ribosome

Ribosomes mediate protein synthesis by decoding the information carried by messenger RNAs (mRNAs) and catalysing peptide bond formation between amino acids. When bacterial ribosomes stall on incomplete messages, the trans-translation quality control mechanism is activated by the transfer-messenger RNA bound to small protein B (tmRNA-SmpB ribonucleoprotein complex). Trans-translation liberates the stalled ribosomes and triggers degradation of the incomplete proteins. Here, we present the cryo-electron microscopy structures of tmRNA-SmpB accommodated or translocated into stalled ribosomes. Two atomic models for each state are proposed. This study reveals how tmRNA-SmpB crosses the ribosome and how, as the problematic mRNA is ejected, the tmRNA resume codon is placed onto the ribosomal decoding site by new contacts between SmpB and the nucleotides upstream of the tag-encoding sequence. This provides a structural basis for the transit of the large tmRNA-SmpB complex through the ribosome and for the means by which the tmRNA internal frame is set for translation to resume.

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

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

Consistent refinement of submitted models at CASP using a knowledge-based potential

Protein structure refinement is an important but unsolved problem; it must be solved if we are to predict biological function that is very sensitive to structural details. Specifically, critical assessment of techniques for protein structure prediction (CASP) shows that the accuracy of predictions in the comparative modeling category is often worse than that of the template on which the homology model is based. Here we describe a refinement protocol that is able to consistently refine submitted predictions for all categories at CASP7. The protocol uses direct energy minimization of the knowledge-based potential of mean force that is based on the interaction statistics of 167 atom types (Summa and Levitt, Proc Natl Acad Sci USA 2007; 104:3177-3182). Our protocol is thus computationally very efficient; it only takes a few minutes of CPU time to run typical protein models (300 residues). We observe an average structural improvement of 1% in GDT_TS, for predictions that have low and medium homology to known PDB structures (Global Distance Test score or GDT_TS between 50 and 80%). We also observe a marked improvement in the stereochemistry of the models. The level of improvement varies amongst the various participants at CASP, but we see large improvements (>10% increase in GDT_TS) even for models predicted by the best performing groups at CASP7. In addition, our protocol consistently improved the best predicted models in the refinement category at CASP7 and CASP8. These improvements in structure and stereochemistry prove the usefulness of our computationally inexpensive, powerful and automatic refinement protocol.

The role of L1 stalk-tRNA interaction in the ribosome elongation cycle

The ribosomal L1 stalk is a mobile structure implicated in directing tRNA movement during translocation through the ribosome. This article investigates three aspects of L1 stalk-tRNA interaction. First, by combining data from cryo electron microscopy, X-ray crystallography, and molecular dynamics simulations through the molecular dynamics flexible fitting method, we obtained atomic models of different tRNAs occupying the hybrid P/E state interacting with the L1 stalk. These models confirm the assignment of fluorescence resonance energy transfer states from previous single-molecule investigations of L1 stalk dynamics. Second, the models reconcile how initiator tRNA(fMet) interacts less strongly with the L1 stalk compared to elongator tRNA(Phe), as seen in previous single-molecule experiments. Third, results from a simulation of the entire ribosome in which the L1 stalk is moved from a half-closed conformation to its open conformation are found to support the hypothesis that L1 stalk opening is involved in tRNA release from the ribosome.

Energy Transfer Dynamics in an RC-LH1-PufX Tubular Photosynthetic Membrane

Light absorption and the subsequent transfer of excitation energy are the first two steps of the photosynthetic process, carried out by protein-bound pigments, mainly bacteriochlorophylls (BChls), in photosynthetic bacteria. BChls are anchored in light-harvesting (LH) complexes, such as light-harvesting complex I (LH1), which directly associates with the reaction center (RC), forming the RC-LH1 core complex. In Rhodobacter sphaeroides, RC-LH1 core complexes contain an additional protein, PufX, and assemble into dimeric RC-LH1-PufX core complexes. In the absence of light-harvesting complexes II, the former complexes can aggregate into a helically ordered tubular photosynthetic membrane. We examined the excitation transfer dynamics in a single RC-LH1-PufX core complex dimer using the hierarchical equations of motion for dissipative quantum dynamics that accurately, yet computationally costly, treat the coupling between BChls and their protein environment. A widely employed description, generalized Förster theory, was also used to calculate the transfer rates of the same excitonic system in order to verify the accuracy of this computationally cheap method. Additionally, in light of the structural uncertainties in the Rhodobacter sphaeroides RC-LH1-PufX core complex, geometrical alterations were introduced in the BChl organization. It is shown that the energy transfer dynamics is not affected by the considered changes in the BChl organization, and that generalized Förster theory provides accurate transfer rates. An all-atom model for a tubular photosynthetic membrane is then constructed on the basis of electron microscopy data, and the overall energy transfer properties of this membrane are computed.

The actin-myosin interface

In order to understand the mechanism of muscle contraction at the atomic level, it is necessary to understand how myosin binds to actin in a reversible way. We have used a novel molecular dynamics technique constrained by an EM map of the actin-myosin complex at 13-A resolution to obtain an atomic model of the strong-binding (rigor) actin-myosin interface. The constraining force resulting from the EM map during the molecular dynamics simulation was sufficient to convert the myosin head from the initial weak-binding state to the strong-binding (rigor) state. Our actin-myosin model suggests extensive contacts between actin and the myosin head (S1). S1 binds to two actin monomers. The contact surface between actin and S1 has increased dramatically compared with previous models. A number of loops in S1 and actin are involved in establishing the interface. Our model also suggests how the loop carrying the critical Arg 405 Glu mutation in S1 found in a familial cardiomyopathy might be functionally involved.

The Ribosome Comes Alive

This essay is a reflection on the ways the X-ray structures of the ribosome are helping in the interpretation of cryogenic electron microscopy (cryo-EM) density maps showing the translating ribosome in motion. Through advances in classification methods, cryo-EM and single-particle reconstruction methods have recently evolved to the point where they can yield an array of structures from a single sample ("story in a sample"), providing snapshots of an entire subprocess of translation, such as translocation or decoding.

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.

Recognition of the regulatory nascent chain TnaC by the ribosome

Regulatory nascent chains interact with the ribosomal exit tunnel and modulate their own translation. To characterize nascent chain recognition by the ribosome at the atomic level, extensive molecular dynamics simulations of TnaC, the leader peptide of the tryptophanase operon, inside the exit tunnel were performed for an aggregate time of 2.1 mus. The simulations, complemented by quantum chemistry calculations, suggest that the critical TnaC residue W12 is recognized by the ribosome via a cation-pi interaction, whereas TnaC's D16 forms salt bridges with ribosomal proteins. The simulations also show that TnaC-mediated translational arrest does not involve a swinging of ribosomal protein L22, as previously proposed. Furthermore, bioinformatic analyses and simulations suggest nascent chain elements that may prevent translational arrest in various organisms. Altogether, the current study unveils atomic-detail interactions that explain the role of elements of TnaC and the ribosome essential for translational arrest.

Self-assembly of photosynthetic membranes

Bacterial photosynthetic membranes, also known as chromatophores, are tightly packed with integral membrane proteins that work together to carry out photosynthesis. Chromatophores display a wide range of cellular morphologies; spherical, tubular, and lamellar chromatophores have all been observed in different bacterial species, or with different protein constituents. Through recent computational modeling and simulation, it has been demonstrated that the light-harvesting complexes abundant in chromatophores induce local membrane curvature via multiple mechanisms. These protein complexes assemble to generate a global curvature and sculpt the chromatophores into various cellular-scale architectures.

From DNA to proteins via the ribosome: structural insights into the workings of the translation machinery

Understanding protein synthesis in bacteria and humans is important for understanding the origin of many human diseases and devising treatments for them. Over the past decade, the field of structural biology has made significant advances in the visualisation of the molecular machinery involved in protein synthesis. It is now possible to discern, at least in outline, the way that interlocking ribosomal components and factors adapt their conformations throughout this process. The determination of structures in various functional contexts, along with the application of kinetic and fluorescent resonance energy transfer approaches to the problem, has given researchers the frame of reference for what remains as the greatest challenge: the complete dynamic portrait of protein synthesis in the cell.

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

Many three-dimensional density maps of 70S ribosome at various functional states are available now in the Electron Microscopy DataBank at EBI. We used our new flexible-fitting approach to systematically analyze these maps to reveal the global conformational differences between the EM structures. Large-scale ratchet-like deformations were observed in an EM structure of the initiation complex and in some EM structures bound with EFG, RF3, and RRF. In most of them, the L1 stalk, which interacts with the tRNA molecule at the E site of ribosome and is considered to be involved in the release of the tRNA, was in "the blocking state" for the E-tRNA. Furthermore, we found that the EM structures bound with EFG or RRF were aligned in the conformational space, suggesting that the large-scale conformational changes of the 70S ribosome bound with these factors occur along a specific pathway in a concerted manner.

Specific arrangement of alpha-helical coiled coils in the core domain of the bacterial flagellar hook for the universal joint function

The bacterial flagellar hook is a short, highly curved tubular structure connecting the rotary motor to the filament acting as a helical propeller. The bending flexibility of the hook allows it to work as a universal joint. A partial atomic model of the hook revealed a sliding intersubunit domain interaction along the protofilament to produce bending flexibility. However, it remained unclear how the tightly packed inner core domains can still permit axial extension and compression. We report advances in cryoEM image analysis for high-resolution, high-throughput structural analysis and a density map of the hook that reveals most of the secondary structures, including the terminal alpha helices forming a coiled coil. The orientations and axial packing interactions of these two alpha helices are distinctly different from those of the filament, allowing them to have a room for axial compression and extension for bending flexibility without impairing the mechanical stability of the hook.

Regulation of the protein-conducting channel by a bound ribosome

During protein synthesis, it is often necessary for the ribosome to form a complex with a membrane-bound channel, the SecY/Sec61 complex, in order to translocate nascent proteins across a cellular membrane. Structural data on the ribosome-channel complex are currently limited to low-resolution cryo-electron microscopy maps, including one showing a bacterial ribosome bound to a monomeric SecY complex. Using that map along with available atomic-level models of the ribosome and SecY, we have determined, through molecular dynamics flexible fitting (MDFF), an atomic-resolution model of the ribosome-channel complex. We characterized computationally the sites of ribosome-SecY interaction within the complex and determined the effect of ribosome binding on the SecY channel. We also constructed a model of a ribosome in complex with a SecY dimer by adding a second copy of SecY to the MDFF-derived model. The study involved 2.7-million-atom simulations over altogether nearly 50 ns.

Elongation in translation as a dynamic interaction among the ribosome, tRNA, and elongation factors EF-G and EF-Tu

The ribosome is a complex macromolecular machine that translates the message encoded in the messenger RNA and synthesizes polypeptides by linking the individual amino acids carried by the cognate transfer RNAs (tRNAs). The protein elongation cycle, during which the tRNAs traverse the ribosome in a coordinated manner along a path of more than 100 A, is facilitated by large-scale rearrangements of the ribosome. These rearrangements go hand in hand with conformational changes of tRNA as well as elongation factors EF-Tu and EF-G - GTPases that catalyze tRNA delivery and translocation, respectively. This review focuses on the structural data related to the dynamics of the ribosomal machinery, which are the basis, in conjunction with existing biochemical, kinetic, and fluorescence resonance energy transfer data, of our knowledge of the decoding and translocation steps of protein elongation.

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.

Discovery through the computational microscope

All-atom molecular dynamics simulations have become increasingly popular as a tool to investigate protein function and dynamics. However, researchers are concerned about the short time scales covered by simulations, the apparent impossibility to model large and integral biomolecular systems, and the actual predictive power of the molecular dynamics methodology. Here we review simulations that were in the past both hotly disputed and considered key successes, namely of proteins with mainly mechanical functions (titin, fibrinogen, ankyrin, and cadherin). The simulation work covered shows how state-of-the-art modeling alleviates some of the prior concerns and how unrefuted discoveries are made through the "computational microscope."

Multiparticle cryo-EM of ribosomes

As the resolution of cryo-EM reconstructions has improved to the subnanometer range, conformational and compositional heterogeneity have become increasing problems in cryo-EM, limiting the resolution of reconstructions. Since further purification is not feasible, the presence of several conformational states of ribosomal complexes in thermodynamic equilibrium requires methods for separating these states in silico. We describe a procedure for generating subnanometer resolution cryo-EM structures from large sets of projection images of ribosomal complexes. The incremental K-means-like method of unsupervised 3D sorting discussed here allows separation of classes in the dataset by exploiting intrinsic divisions in the data. The classification procedure is described in detail and its effectiveness is illustrated using current examples from our work. Through a good separation of conformational modes, higher resolution reconstructions can be calculated. This increases information gained from single states, while exploiting the coexistence of multiple states to gather comprehensive mechanistic insight into biological processes like ribosomal translocation.

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.

Semi-isometric registration of line features for flexible fitting of protein structures

In this paper, we study a registration problem that is motivated by a practical biology problem - fitting protein structures to low-resolution density maps. We consider registration between two sets of lines features (e.g., helices in the proteins) that have undergone not a single, but multiple isometric transformations (e.g., hinge-motions). The problem is further complicated by the presence of symmetry in each set. We formulate the problem as a clique-finding problem in a product graph, and propose a heuristic solution that includes a fast clique-finding algorithm unique to the structure of this graph. When tested on a suite of real protein structures, the algorithm achieved high accuracy even for very large inputs containing hundreds of helices.

An automated procedure for detecting protein folds from sub-nanometer resolution electron density

The use of sub-nanometer resolution electron density as spatial constraints for de novo and ab initio structure prediction requires knowledge of protein boundaries to accurately segment the electron density for the prediction algorithms. Here we present a procedure where even poorly segmented density can be used to determine the fold of the protein. The method is automated, fast, capable of searching for multiple copies of a protein fold, and accessible to densities encompassing more than a thousand residues. The automation is particularly powerful as it allows the procedure to take full advantage of the expanding repository in the Protein Data Bank. We have tested the method on nine segmented sub-nanometer image reconstruction electron densities. The method successfully identifies the correct fold for the six densities for which an atomic structure is known, identifies a fold that agrees with prior structural data, a fold that agrees with predictions from the Fold & Function Assignment server, and a fold that correlates with secondary structure prediction. The identified folds in the last three examples can be used as templates for comparative modeling of the bacteriophage P22 tail-machine (a 3MDa complex composed of 39 protein subunits).

EM-fold: De novo folding of alpha-helical proteins guided by intermediate-resolution electron microscopy density maps

In medium-resolution (7-10 A) cryo-electron microscopy (cryo-EM) density maps, alpha helices can be identified as density rods whereas beta-strand or loop regions are not as easily discerned. We are proposing a computational protein structure prediction algorithm "EM-Fold" that resolves the density rod connectivity ambiguity by placing predicted alpha helices into the density rods and adding missing backbone coordinates in loop regions. In a benchmark of 11 mainly alpha-helical proteins of known structure a native-like model is identified in eight cases (rmsd 3.9-7.9 A). The three failures can be attributed to inaccuracies in the secondary structure prediction step that precedes EM-Fold. EM-Fold has been applied to the approximately 6 A resolution cryo-EM density map of protein IIIa from human adenovirus. We report the first topological model for the alpha-helical 400 residue N-terminal region of protein IIIa. EM-Fold also has the potential to interpret medium-resolution density maps in X-ray crystallography.

Molecular cryo-electron tomography of vitreous tissue sections: current challenges

Electron tomography of vitreous tissue sections (tissue TOVIS) allows the study of the three-dimensional structure of molecular complexes in a near-native cellular context. Its usage is, however, limited by an unfortunate combination of noisy and incomplete data, by a technically demanding sample preparation procedure, and by a disposition for specimen degradation during data collection. Here we outline some major challenges as experienced from the application of TOVIS to human skin. We further consider a number of practical measures as well as theoretical approaches for its future development.

Structure of monomeric yeast and mammalian Sec61 complexes interacting with the translating ribosome

The trimeric Sec61/SecY complex is a protein-conducting channel (PCC) for secretory and membrane proteins. Although Sec complexes can form oligomers, it has been suggested that a single copy may serve as an active PCC. We determined subnanometer-resolution cryo-electron microscopy structures of eukaryotic ribosome-Sec61 complexes. In combination with biochemical data, we found that in both idle and active states, the Sec complex is not oligomeric and interacts mainly via two cytoplasmic loops with the universal ribosomal adaptor site. In the active state, the ribosomal tunnel and a central pore of the monomeric PCC were occupied by the nascent chain, contacting loop 6 of the Sec complex. This provides a structural basis for the activity of a solitary Sec complex in cotranslational protein translocation.

Structural insight into nascent polypeptide chain-mediated translational stalling

Expression of the Escherichia coli tryptophanase operon depends on ribosome stalling during translation of the upstream TnaC leader peptide, a process for which interactions between the TnaC nascent chain and the ribosomal exit tunnel are critical. We determined a 5.8 angstrom-resolution cryo-electron microscopy and single-particle reconstruction of a ribosome stalled during translation of the tnaC leader gene. The nascent chain was extended within the exit tunnel, making contacts with ribosomal components at distinct sites. Upon stalling, two conserved residues within the peptidyltransferase center adopted conformations that preclude binding of release factors. We propose a model whereby interactions within the tunnel are relayed to the peptidyltransferase center to inhibit translation. Moreover, we show that nascent chains adopt distinct conformations within the ribosomal exit tunnel.

Molecular dynamics flexible fitting: a practical guide to combine cryo-electron microscopy and X-ray crystallography

Hybrid computational methods for combining structural data from different sources and resolutions are becoming an essential part of structural biology, especially as the field moves toward the study of large macromolecular assemblies. We have developed the molecular dynamics flexible fitting (MDFF) method for combining high-resolution atomic structures with cryo-electron microscopy (cryo-EM) maps, that results in atomic models representing the conformational state captured by cryo-EM. The method has been applied successfully to the ribosome, a ribonucleoprotein complex responsible for protein synthesis. MDFF involves a molecular dynamics simulation in which a guiding potential, based on the cryo-EM map, is added to the standard force field. Forces proportional to the gradient of the density map guide an atomic structure, available from X-ray crystallography, into high-density regions of a cryo-EM map. In this paper we describe the necessary steps to set up, run, and analyze MDFF simulations and the software packages that implement the corresponding functionalities.

Single-particle reconstruction of biological macromolecules in electron microscopy--30 years

This essay gives the autho's personal account on the development of concepts underlying single-particle reconstruction, a technique in electron microscopy of macromolecular assemblies with a remarkable record of achievements as of late. The ribosome proved to be an ideal testing ground for the development of specimen preparation methods, cryo-EM techniques, and algorithms, with discoveries along the way as a rich reward. Increasingly, cryo-EM and single-particle reconstruction, in combination with classification techniques, is revealing dynamic information on functional molecular machines uninhibited by molecular contacts.

Protein-induced membrane curvature investigated through molecular dynamics flexible fitting

In the photosynthetic purple bacterium Rhodobacter (Rba.) sphaeroides, light is absorbed by membrane-bound light-harvesting (LH) proteins LH1 and LH2. LH1 directly surrounds the reaction center (RC) and, together with PufX, forms a dimeric (RC-LH1-PufX)2 protein complex. In LH2-deficient Rba. sphaeroides mutants, RC-LH1-PufX dimers aggregate into tubular vesicles with a radius of approximately 250-550 A, making RC-LH1-PufX one of the few integral membrane proteins known to actively induce membrane curvature. Recently, a three-dimensional electron microscopy density map showed that the Rba. sphaeroides RC-LH1-PufX dimer exhibits a prominent bend at its dimerizing interface. To investigate the curvature properties of this highly bent protein, we employed molecular dynamics simulations to fit an all-atom structural model of the RC-LH1-PufX dimer within the electron microscopy density map. The simulations reveal how the dimer produces a membrane with high local curvature, even though the location of PufX cannot yet be determined uniquely. The resulting membrane curvature agrees well with the size of RC-LH1-PufX tubular vesicles, and demonstrates how the local curvature properties of the RC-LH1-PufX dimer propagate to form the observed long-range organization of the Rba. sphaeroides tubular vesicles.

Refinement of protein structures into low-resolution density maps using rosetta

We describe a method based on Rosetta structure refinement for generating high-resolution, all-atom protein models from electron cryomicroscopy density maps. A local measure of the fit of a model to the density is used to directly guide structure refinement and to identify regions incompatible with the density that are then targeted for extensive rebuilding. Over a range of test cases using both simulated and experimentally generated data, the method consistently increases the accuracy of starting models generated either by comparative modeling or by hand-tracing the density. The method can achieve near-atomic resolution starting from density maps at 4-6 A resolution.

Interpretation of very low resolution X-ray electron-density maps using core objects

A novel approach to obtaining structural information from macromolecular X-ray data extending to resolutions as low as 20 A is presented. Following a simple map-segmentation procedure, the approximate shapes of the domains forming the structure are identified. A pattern-recognition comparative analysis of these shapes and those derived from the structures of domains from the PDB results in candidate structural models that can be used for a fit into the density map. It is shown that the placed candidate models can be employed for subsequent phase extension to higher resolution.

Confidence intervals for fitting of atomic models into low-resolution densities

The fitting of high-resolution structures into low-resolution densities obtained from techniques such as electron microscopy or small-angle X-ray scattering can yield powerful new insights. While several algorithms for achieving optimal fits have recently been developed, relatively little effort has been devoted to developing objective measures for judging the quality of the resulting fits, in particular with regard to the danger of overfitting. Here, a general method is presented for obtaining confidence intervals for atomic coordinates resulting from fitting of atomic resolution domain structures into low-resolution densities using well established statistical tools. It is demonstrated that the resulting confidence intervals are sufficiently accurate to allow meaningful statistical tests and to provide tools for detecting potential overfitting.

Inferential optimization for simultaneous fitting of multiple components into a CryoEM map of their assembly

Models of macromolecular assemblies are essential for a mechanistic description of cellular processes. Such models are increasingly obtained by fitting atomic-resolution structures of components into a density map of the whole assembly. Yet, current density-fitting techniques are frequently insufficient for an unambiguous determination of the positions and orientations of all components. Here, we describe MultiFit, a method used for simultaneously fitting atomic structures of components into their assembly density map at resolutions as low as 25 A. The component positions and orientations are optimized with respect to a scoring function that includes the quality-of-fit of components in the map, the protrusion of components from the map envelope, and the shape complementarity between pairs of components. The scoring function is optimized by our exact inference optimizer DOMINO (Discrete Optimization of Multiple INteracting Objects) that efficiently finds the global minimum in a discrete sampling space. MultiFit was benchmarked on seven assemblies of known structure, consisting of up to seven proteins each. The input atomic structures of the components were obtained from the Protein Data Bank, as well as by comparative modeling based on a 16-99% sequence identity to a template structure. A near-native configuration was usually found as the top-scoring model. Therefore, MultiFit can provide initial configurations for further refinement of many multicomponent assembly structures described by electron microscopy.

Hybrid approaches: applying computational methods in cryo-electron microscopy

Recent advances in cryo-electron microscopy have led to an increasing number of high (3-5A) to medium (5-10A) resolution cryoEM density maps. These density maps contain valuable information about the protein structure but frequently require computational algorithms to aid their structural interpretation. It is these hybrid approaches between cryoEM and computational protein structure prediction algorithms that will shape protein structure elucidation from density maps.

Structural model and excitonic properties of the dimeric RC-LH1-PufX complex from Rhodobacter sphaeroides

The light-harvesting apparatus of the purple bacterial photosynthetic unit consists of a pool of peripheral light-harvesting complexes that transfer excitation energy to a reaction center (RC) via the surrounding pigment-protein complex LH1. Recent electron microscopy and atomic force microscopy studies have revealed that RC-LH1 units of Rhodobacter sphaeroides form membrane-bending dimeric complexes together with the polypeptide PufX. We present a structural model for these RC-LH1-PufX dimeric complexes constructed using the molecular dynamics flexible fitting method based on an EM density map. The arrangement of the LH1 BChls displays a distortion near the proposed location of the PufX polypeptide. The resulting atomic model for BChl arrays is used to compute the excitonic properties of the dimeric RC-LH1 complex. A comparison is presented between the structural and excitonic features of the S-shaped dimeric BChl array of Rhodobacter sphaeroides and the circular BChl arrangement found in other purple bacteria.

Introduction. Biomolecular simulation

'Everything that living things do can be understood in terms of the jigglings and wigglings of atoms' as Richard Feynman provocatively stated nearly 50 years ago. But how can we 'see' this wiggling and jiggling and understand how it drives biology? Increasingly, computer simulations of biological macromolecules are helping to meet this challenge.

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

In translation, elongation factor Tu (EF-Tu) molecules deliver aminoacyl-tRNAs to the mRNA-programmed ribosome. The GTPase activity of EF-Tu is triggered by ribosome-induced conformational changes of the factor that play a pivotal role in the selection of the cognate aminoacyl-tRNAs. We present a 6.7-A cryo-electron microscopy map of the aminoacyl-tRNA x EF-Tu x GDP x kirromycin-bound Escherichia coli ribosome, together with an atomic model of the complex obtained through molecular dynamics flexible fitting. The model reveals the conformational changes in the conserved GTPase switch regions of EF-Tu that trigger hydrolysis of GTP, along with key interactions, including those between the sarcin-ricin loop and the P loop of EF-Tu, and between the effector loop of EF-Tu and a conserved region of the 16S rRNA. Our data suggest that GTP hydrolysis on EF-Tu is controlled through a hydrophobic gate mechanism.

The advent of near-atomic resolution in single-particle electron microscopy

Single-particle electron microscopy (EM) can provide structural information for a large variety of biological molecules, ranging from small proteins to large macromolecular assemblies, without the need to produce crystals. The year 2008 has become a landmark year for single-particle EM as for the first time density maps have been produced at a resolution that made it possible to trace protein backbones or even to build atomic models. In this review, we highlight some of the recent successes achieved by single-particle EM and describe the individual steps involved in producing a density map by this technique. We also discuss some of the remaining challenges and areas, in which further advances would have a great impact on the results that can be achieved by single-particle EM.

Flexible fitting of high-resolution x-ray structures into cryoelectron microscopy maps using biased molecular dynamics simulations

A methodology for flexible fitting of all-atom high-resolution structures into low-resolution cryoelectron microscopy (cryo-EM) maps is presented. Flexibility of the modeled structure is simulated by classical molecular dynamics and an additional effective potential is introduced to enhance the fitting process. The additional potential is proportional to the correlation coefficient between the experimental cryo-EM map and a synthetic map generated for an all-atom structure being fitted to the map. The additional forces are calculated as a gradient of the correlation coefficient. During the molecular dynamics simulations under the additional forces, the molecule undergoes a conformational transition that maximizes the correlation coefficient, which results in a high-accuracy fit of all-atom structure into a cryo-EM map. Using five test proteins that exhibit structural rearrangement during their biological activity, we demonstrate performance of our method. We also test our method on the experimental cryo-EM of elongation factor G and show that the model obtained is comparable to previous studies. In addition, we show that overfitting can be avoided by assessing the quality of the fitted model in terms of correlation coefficient and secondary structure preservation.

Exploration of parameters in cryo-EM leading to an improved density map of the E. coli ribosome

A number of image processing parameters in the 3D reconstruction of a ribosome complex from a cryo-EM data set were varied to test their effects on the final resolution. The parameters examined were pixel size, window size, and mode of Fourier amplitude enhancement at high spatial frequencies. In addition, the strategy of switching from large to small pixel size during angular refinement was explored. The relationship between resolution (in Fourier space) and the number of particles was observed to follow a lin-log dependence, a relationship that appears to hold for other data, as well. By optimizing the above parameters, and using a lin-log extrapolation to the full data set in the estimation of resolution from half-sets, we obtained a 3D map from 131,599 ribosome particles at 6.7A resolution (FSC=0.5).

Problems in fitting high resolution structures into electron microscopic reconstructions

Great advances have been made in electron microscopy (EM) over the past decade, with the result that a number of protein complexes have been solved at near-atomic resolution using EM imaging. However, only a limited number of such complexes are expected to have the high degree of internal order needed to achieve this type of resolution. Many other complexes and polymers will be visualized and reconstructed by EM at an intermediate level of resolution, where the polypeptide chain cannot be directly traced. Crystal and nuclear magnetic resonance structures for components or subunits of these higher-order assemblies are frequently available. One of the greatest strengths of EM continues to be the ability to dock high-resolution structures of components into low or intermediate resolution reconstructions of assemblies to build pseudoatomic models for quaternary structure. This review discusses the strengths and limitations of this approach, with particular emphasis on protein polymers. I discuss how limitations in resolution can lead to ambiguities in building models, and these cannot be always be resolved with available data. The use of homology models for quaternary structure are particularly problematic, given accumulating evidence for the divergence of quaternary structures at the same time that tertiary structure can be conserved.

Biomolecular pleiomorphism probed by spatial interpolation of coarse models

In low resolution structures of biological assemblies one can often observe conformational deviations that require a flexible rearrangement of structural domains fitted at the atomic level. We are evaluating interpolation methods for the flexible alignment of atomic models based on coarse models. Spatial interpolation is well established in image-processing and visualization to describe the overall deformation or warping of an object or an image. Combined with a coarse representation of the biological system by feature vectors, such methods can provide a flexible approximation of the molecular structure. We have compared three well-known interpolation techniques and evaluated the results by comparing them with constrained molecular dynamics. One method, inverse distance weighting interpolation, consistently produced models that were nearly indistinguishable on the alpha carbon level from the molecular dynamics results. The method is simple to apply and enables flexing of structures by non-expert modelers. This is useful for the basic interpretation of volumetric data in biological applications such as electron microscopy. The method can be used as a general interpretation tool for sparsely sampled motions derived from coarse models.