The 13Å Structure of a Chaperonin GroEL–Protein Substrate Complex by Cryo-electron Microscopy

Friends in need: How chaperonins recognize and remodel proteins that require folding assistance

Chaperonins are biological nanomachines that help newly translated proteins to fold by rescuing them from kinetically trapped misfolded states. Protein folding assistance by the chaperonin machinery is obligatory in vivo for a subset of proteins in the bacterial proteome. Chaperonins are large oligomeric complexes, with unusual seven fold symmetry (group I) or eight/nine fold symmetry (group II), that form double-ring constructs, enclosing a central cavity that serves as the folding chamber. Dramatic large-scale conformational changes, that take place during ATP-driven cycles, allow chaperonins to bind misfolded proteins, encapsulate them into the expanded cavity and release them back into the cellular environment, regardless of whether they are folded or not. The theory associated with the iterative annealing mechanism, which incorporated the conformational free energy landscape description of protein folding, quantitatively explains most, if not all, the available data. Misfolded conformations are associated with low energy minima in a rugged energy landscape. Random disruptions of these low energy conformations result in higher free energy, less folded, conformations that can stochastically partition into the native state. Two distinct mechanisms of annealing action have been described. Group I chaperonins (GroEL homologues in eubacteria and endosymbiotic organelles), recognize a large number of misfolded proteins non-specifically and operate through highly coordinated cooperative motions. By contrast, the less well understood group II chaperonins (CCT in Eukarya and thermosome/TF55 in Archaea), assist a selected set of substrate proteins. Sequential conformational changes within a CCT ring are observed, perhaps promoting domain-by-domain substrate folding. Chaperonins are implicated in bacterial infection, autoimmune disease, as well as protein aggregation and degradation diseases. Understanding the chaperonin mechanism and the specific proteins they rescue during the cell cycle is important not only for the fundamental aspect of protein folding in the cellular environment, but also for effective therapeutic strategies.

Electron scattering properties of biological macromolecules and their use for cryo-EM map sharpening

Resolution-dependent loss of contrast in cryo-EM maps may obscure features at high resolution that are critical for map interpretation. Post-processing of cryo-EM maps can improve the interpretability by adjusting the resolution-dependence of structure factor amplitudes through map sharpening. Traditionally this has been done by rescaling the relative contribution of low and high-resolution frequencies globally. More recently, the realisation that molecular motion and heterogeneity cause non-uniformity of resolution throughout the map has inspired the development of techniques that optimise sharpening locally. We previously developed LocScale, a method that utilises the radial structure factor from a refined atomic model as a restraint for local map sharpening. While this method has proved beneficial for the interpretation of cryo-EM maps, the dependence on the availability of (partial) model information limits its general applicability. Here, we review the basic assumptions of resolution-dependent contrast loss in cryo-EM maps and propose a route towards a robust alternative for local map sharpening that utilises information on expected scattering properties of biological macromolecules, but requires no detailed knowledge of the underlying molecular structure. We examine remaining challenges for implementation and discuss possible applications.

Chaperonin-assisted protein folding: a chronologue

This chronologue seeks to document the discovery and development of an understanding of oligomeric ring protein assemblies known as chaperonins that assist protein folding in the cell. It provides detail regarding genetic, physiologic, biochemical, and biophysical studies of these ATP-utilizing machines from both in vivo and in vitro observations. The chronologue is organized into various topics of physiology and mechanism, for each of which a chronologic order is generally followed. The text is liberally illustrated to provide firsthand inspection of the key pieces of experimental data that propelled this field. Because of the length and depth of this piece, the use of the outline as a guide for selected reading is encouraged, but it should also be of help in pursuing the text in direct order.

A two-domain folding intermediate of RuBisCO in complex with the GroEL chaperonin

The chaperonins (GroEL and GroES in Escherichia coli) are ubiquitous molecular chaperones that assist a subset of essential substrate proteins to undergo productive folding to the native state. Using single particle cryo EM and image processing we have examined complexes of E. coli GroEL with the stringently GroE-dependent substrate enzyme RuBisCO from Rhodospirillum rubrum. Here we present snapshots of non-native RuBisCO - GroEL complexes. We observe two distinct substrate densities in the binary complex reminiscent of the two-domain structure of the RuBisCO subunit, so that this may represent a captured form of an early folding intermediate. The occupancy of the complex is consistent with the negative cooperativity of GroEL with respect to substrate binding, in accordance with earlier mass spectroscopy studies.

Automated map sharpening by maximization of detail and connectivity

An algorithm for automatic map sharpening is presented that is based on optimization of the detail and connectivity of the sharpened map. The detail in the map is reflected in the surface area of an iso-contour surface that contains a fixed fraction of the volume of the map, where a map with high level of detail has a high surface area. The connectivity of the sharpened map is reflected in the number of connected regions defined by the same iso-contour surfaces, where a map with high connectivity has a small number of connected regions. By combining these two measures in a metric termed the `adjusted surface area', map quality can be evaluated in an automated fashion. This metric was used to choose optimal map-sharpening parameters without reference to a model or other interpretations of the map. Map sharpening by optimization of the adjusted surface area can be carried out for a map as a whole or it can be carried out locally, yielding a locally sharpened map. To evaluate the performance of various approaches, a simple metric based on map-model correlation that can reproduce visual choices of optimally sharpened maps was used. The map-model correlation is calculated using a model with B factors (atomic displacement factors; ADPs) set to zero. This model-based metric was used to evaluate map sharpening and to evaluate map-sharpening approaches, and it was found that optimization of the adjusted surface area can be an effective tool for map sharpening.

Model-based local density sharpening of cryo-EM maps

Atomic models based on high-resolution density maps are the ultimate result of the cryo-EM structure determination process. Here, we introduce a general procedure for local sharpening of cryo-EM density maps based on prior knowledge of an atomic reference structure. The procedure optimizes contrast of cryo-EM densities by amplitude scaling against the radially averaged local falloff estimated from a windowed reference model. By testing the procedure using six cryo-EM structures of TRPV1, β-galactosidase, γ-secretase, ribosome-EF-Tu complex, 20S proteasome and RNA polymerase III, we illustrate how local sharpening can increase interpretability of density maps in particular in cases of resolution variation and facilitates model building and atomic model refinement.

GroEL actively stimulates folding of the endogenous substrate protein PepQ

Many essential proteins cannot fold without help from chaperonins, like the GroELS system of Escherichia coli. How chaperonins accelerate protein folding remains controversial. Here we test key predictions of both passive and active models of GroELS-stimulated folding, using the endogenous E. coli metalloprotease PepQ. While GroELS increases the folding rate of PepQ by over 15-fold, we demonstrate that slow spontaneous folding of PepQ is not caused by aggregation. Fluorescence measurements suggest that, when folding inside the GroEL-GroES cavity, PepQ populates conformations not observed during spontaneous folding in free solution. Using cryo-electron microscopy, we show that the GroEL C-termini make physical contact with the PepQ folding intermediate and help retain it deep within the GroEL cavity, resulting in reduced compactness of the PepQ monomer. Our findings strongly support an active model of chaperonin-mediated protein folding, where partial unfolding of misfolded intermediates plays a key role.

In the prevailing model for assisted protein folding, chaperonins act passively by preventing protein aggregation. Here, the authors use single-molecule fluorescence measurements and cryo-electron microscopy and show that the E. coli GroELS chaperonin system also has an active role in folding the endogenous bacterial protein PepQ.

ATP-driven molecular chaperone machines

This review is focused on the mechanisms by which ATP binding and hydrolysis drive chaperone machines assisting protein folding and unfolding. A survey of the key, general chaperone systems Hsp70 and Hsp90, and the unfoldase Hsp100 is followed by a focus on the Hsp60 chaperonin machine which is understood in most detail. Cryo-electron microscopy analysis of the E. coli Hsp60 GroEL reveals intermediate conformations in the ATPase cycle and in substrate folding. These structures suggest a mechanism by which GroEL can forcefully unfold and then encapsulate substrates for subsequent folding in isolation from all other binding surfaces.

Allosteric Mechanisms in Chaperonin Machines

Chaperonins are nanomachines that facilitate protein folding by undergoing energy (ATP)-dependent movements that are coordinated in time and space owing to complex allosteric regulation. They consist of two back-to-back stacked oligomeric rings with a cavity at each end where protein substrate folding can take place. Here, we focus on the GroEL/GroES chaperonin system from Escherichia coli and, to a lesser extent, on the more poorly characterized eukaryotic chaperonin CCT/TRiC. We describe their various functional (allosteric) states and how they are affected by substrates and allosteric effectors that include ATP, ADP, nonfolded protein substrates, potassium ions, and GroES (in the case of GroEL). We also discuss the pathways of intra- and inter-ring allosteric communication by which they interconvert and the coupling between allosteric transitions and protein folding reactions.

The GroEL-GroES Chaperonin Machine: A Nano-Cage for Protein Folding

The bacterial chaperonin GroEL and its cofactor GroES constitute the paradigmatic molecular machine of protein folding. GroEL is a large double-ring cylinder with ATPase activity that binds non-native substrate protein (SP) via hydrophobic residues exposed towards the ring center. Binding of the lid-shaped GroES to GroEL displaces the bound protein into an enlarged chamber, allowing folding to occur unimpaired by aggregation. GroES and SP undergo cycles of binding and release, regulated allosterically by the GroEL ATPase. Recent structural and functional studies are providing insights into how the physical environment of the chaperonin cage actively promotes protein folding, in addition to preventing aggregation. Here, we review different models of chaperonin action and discuss issues of current debate.

Following Natures Lead: On the Construction of Membrane-Inserted Toxins in Lipid Bilayer Nanodiscs

Bacterial toxin or viral entry into the cell often requires cell surface binding and endocytosis. The endosomal acidification induces a limited unfolding/refolding and membrane insertion reaction of the soluble toxins or viral proteins into their translocation competent or membrane inserted states. At the molecular level, the specific orientation and immobilization of the pre-transitioned toxin on the cell surface is often an important prerequisite prior to cell entry. We propose that structures of some toxin membrane insertion complexes may be observed through procedures where one rationally immobilizes the soluble toxin so that potential unfolding ↔ refolding transitions that occur prior to membrane insertion orientate away from the immobilization surface in the presence of lipid micelle pre-nanodisc structures. As a specific example, the immobilized prepore form of the anthrax toxin pore translocon or protective antigen can be transitioned, inserted into a model lipid membrane (nanodiscs), and released from the immobilized support in its membrane solubilized form. This particular strategy, although unconventional, is a useful procedure for generating pure membrane-inserted toxins in nanodiscs for electron microscopy structural analysis. In addition, generating a similar immobilized platform on label-free biosensor surfaces allows one to observe the kinetics of these acid-induced membrane insertion transitions. These platforms can facilitate the rational design of inhibitors that specifically target the toxin membrane insertion transitions that occur during endosomal acidification. This approach may lead to a new class of direct anti-toxin inhibitors.

Probing the kinetic stabilities of Friedreich's ataxia clinical variants using a solid phase GroEL chaperonin capture platform

Numerous human diseases are caused by protein folding defects where the protein may become more susceptible to degradation or aggregation. Aberrant protein folding can affect the kinetic stability of the proteins even if these proteins appear to be soluble in vivo. Experimental discrimination between functional properly folded and misfolded nonfunctional conformers is not always straightforward at near physiological conditions. The differences in the kinetic behavior of two initially folded frataxin clinical variants were examined using a high affinity chaperonin kinetic trap approach at 25 °C. The kinetically stable wild type frataxin (FXN) shows no visible partitioning onto the chaperonin. In contrast, the clinical variants FXN-p.Asp122Tyr and FXN-p.Ile154Phe kinetically populate partial folded forms that tightly bind the GroEL chaperonin platform. The initially soluble FXN-p.Ile154Phe variant partitions onto GroEL more rapidly and is more kinetically liable. These differences in kinetic stability were confirmed using differential scanning fluorimetry. The kinetic and aggregation stability differences of these variants may lead to the distinct functional impairments described in Friedreich's ataxia, the neurodegenerative disease associated to frataxin functional deficiency. This chaperonin platform approach may be useful for identifying small molecule stabilizers since stabilizing ligands to frataxin variants should lead to a concomitant decrease in chaperonin binding.

High-resolution modeling of protein structures based on flexible fitting of low-resolution structural data

To circumvent the difficulty of directly solving high-resolution biomolecular structures, low-resolution structural data from Cryo-electron microscopy (EM) and small angle solution X-ray scattering (SAXS) are increasingly used to explore multiple conformational states of biomolecular assemblies. One promising venue to obtain high-resolution structural models from low-resolution data is via data-constrained flexible fitting. To this end, we have developed a new method based on a coarse-grained Cα-only protein representation, and a modified form of the elastic network model (ENM) that allows large-scale conformational changes while maintaining the integrity of local structures including pseudo-bonds and secondary structures. Our method minimizes a pseudo-energy which linearly combines various terms of the modified ENM energy with an EM/SAXS-fitting score and a collision energy that penalizes steric collisions. Unlike some previous flexible fitting efforts using the lowest few normal modes, our method effectively utilizes all normal modes so that both global and local structural changes can be fully modeled with accuracy. This method is also highly efficient in computing time. We have demonstrated our method using adenylate kinase as a test case which undergoes a large open-to-close conformational change. The EM-fitting method is available at a web server (http://enm.lobos.nih.gov), and the SAXS-fitting method is available as a pre-compiled executable upon request.

Bayesian analysis of individual electron microscopy images: towards structures of dynamic and heterogeneous biomolecular assemblies

We develop a method to extract structural information from electron microscopy (EM) images of dynamic and heterogeneous molecular assemblies. To overcome the challenge of disorder in the imaged structures, we analyze each image individually, avoiding information loss through clustering or averaging. The Bayesian inference of EM (BioEM) method uses a likelihood-based probabilistic measure to quantify the consistency between each EM image and given structural models. The likelihood function accounts for uncertainties in the molecular position and orientation, variations in the relative intensities and noise in the experimental images. The BioEM formalism is physically intuitive and mathematically simple. We show that for experimental GroEL images, BioEM correctly identifies structures according to the functional state. The top-ranked structure is the corresponding X-ray crystal structure, followed by an EM structure generated previously from a superset of the EM images used here. To analyze EM images of highly flexible molecules, we propose an ensemble refinement procedure, and validate it with synthetic EM maps of the ESCRT-I-II supercomplex. Both the size of the ensemble and its structural members are identified correctly. BioEM offers an alternative to 3D-reconstruction methods, extracting accurate population distributions for highly flexible structures and their assemblies. We discuss limitations of the method, and possible applications beyond ensemble refinement, including the cross-validation and unbiased post-assessment of model structures, and the structural characterization of systems where traditional approaches fail. Overall, our results suggest that the BioEM framework can be used to analyze EM images of both ordered and disordered molecular systems.

Visualizing GroEL/ES in the act of encapsulating a folding protein

The GroEL/ES chaperonin system is required for the assisted folding of many proteins. How these substrate proteins are encapsulated within the GroEL-GroES cavity is poorly understood. Using symmetry-free, single-particle cryo-electron microscopy, we have characterized a chemically modified mutant of GroEL (EL43Py) that is trapped at a normally transient stage of substrate protein encapsulation. We show that the symmetric pattern of the GroEL subunits is broken as the GroEL cis-ring apical domains reorient to accommodate the simultaneous binding of GroES and an incompletely folded substrate protein (RuBisCO). The collapsed RuBisCO folding intermediate binds to the lower segment of two apical domains, as well as to the normally unstructured GroEL C-terminal tails. A comparative structural analysis suggests that the allosteric transitions leading to substrate protein release and folding involve concerted shifts of GroES and the GroEL apical domains and C-terminal tails.

Consensus among multiple approaches as a reliability measure for flexible fitting into cryo-EM data

Cryo-electron microscopy (cryo-EM) can provide low-resolution density maps of large macromolecular assemblies. As the number of structures deposited in the Protein Data Bank by fitting a high-resolution structure into a low-resolution cryo-EM map is increasing, there is a need to revise the protocols and improve the measures for fitting. A recent study suggested using a combination of multiple automated flexible fitting approaches to improve the interpretation of cryo-EM data. The current work further explores the use of multiple approaches by validating this "consensus" fitting approach and deriving a local reliability measure. Here four different flexible fitting approaches are applied for fitting an initial structure into a simulated density map of known target structure from a dataset of proteins. It is found that the models produced from different approaches often have a consensus in conformation and are also near to the target structure, whereas cases not showing consensus are away from the target. A high correlation is also observed between the RMSF profiles calculated with respect to the average and the target structures, which indicates that the relation between consensus and accuracy can also be extended to a per-residue level. Therefore, the RMSF among the fitted models is proposed as a local reliability measure, which can be used to assess the reliability of the fit at specific regions. Hence, we encourage the community to use consensus flexible fitting with different methods to report on local reliability of the resulting models and improve the interpretation of cryo-EM data.

Go hybrid: EM, crystallography, and beyond

A mechanistic understanding of the molecular transactions that govern cellular function requires knowledge of the dynamic organization of the macromolecular machines involved in these processes. Structural biologists employ a variety of biophysical methods to study large macromolecular complexes, but no single technique is likely to provide a complete description of the structure-function relationship of all the constituent components. Since structural studies generally only provide snapshots of these dynamic machines as they accomplish their molecular functions, combining data from many methodologies is crucial to our understanding of molecular function.

Fitting multimeric protein complexes into electron microscopy maps using 3D Zernike descriptors

A novel computational method for fitting high-resolution structures of multiple proteins into a cryoelectron microscopy map is presented. The method named EMLZerD generates a pool of candidate multiple protein docking conformations of component proteins, which are later compared with a provided electron microscopy (EM) density map to select the ones that fit well into the EM map. The comparison of docking conformations and the EM map is performed using the 3D Zernike descriptor (3DZD), a mathematical series expansion of three-dimensional functions. The 3DZD provides a unified representation of the surface shape of multimeric protein complex models and EM maps, which allows a convenient, fast quantitative comparison of the three-dimensional structural data. Out of 19 multimeric complexes tested, near native complex structures with a root-mean-square deviation of less than 2.5 Å were obtained for 14 cases while medium range resolution structures with correct topology were computed for the additional 5 cases.

Hybrid approaches to molecular simulation

Molecular dynamics (MD) simulation is an established method for studying the conformational changes that are important for protein function. Recent advances in hardware and software have allowed MD simulations over the same timescales as experiment, improving the agreement between theory and experiment to a large extent. However, running such simulations are costly, in terms of resources, storage, and trajectory analysis. There is still a place for techniques that involve short MD simulations. In order to overcome the sampling paucity of short time-scales, hybrid methods that include some form of MD simulation can exploit certain features of the system of interest, often combining experimental information in surprising ways. Here, we review some recent hybrid approaches to the simulation of proteins.

Accurate flexible fitting of high-resolution protein structures to small-angle x-ray scattering data using a coarse-grained model with implicit hydration shell

Small-angle x-ray scattering (SAXS) is a powerful technique widely used to explore conformational states and transitions of biomolecular assemblies in solution. For accurate model reconstruction from SAXS data, one promising approach is to flexibly fit a known high-resolution protein structure to low-resolution SAXS data by computer simulations. This is a highly challenging task due to low information content in SAXS data. To meet this challenge, we have developed what we believe to be a novel method based on a coarse-grained (one-bead-per-residue) protein representation and a modified form of the elastic network model that allows large-scale conformational changes while maintaining pseudobonds and secondary structures. Our method optimizes a pseudoenergy that combines the modified elastic-network model energy with a SAXS-fitting score and a collision energy that penalizes steric collisions. Our method uses what we consider a new implicit hydration shell model that accounts for the contribution of hydration shell to SAXS data accurately without explicitly adding waters to the system. We have rigorously validated our method using five test cases with simulated SAXS data and three test cases with experimental SAXS data. Our method has successfully generated high-quality structural models with root mean-squared deviation of 1 ∼ 3 Å from the target structures.

Consensus among flexible fitting approaches improves the interpretation of cryo-EM data

Cryo-elecron microscopy (cryo-EM) can provide important structural information of large macromolecular assemblies in different conformational states. Recent years have seen an increase in structures deposited in the Protein Data Bank (PDB) by fitting a high-resolution structure into its low-resolution cryo-EM map. A commonly used protocol for accommodating the conformational changes between the X-ray structure and the cryo-EM map is rigid body fitting of individual domains. With the emergence of different flexible fitting approaches, there is a need to compare and revise these different protocols for the fitting. We have applied three diverse automated flexible fitting approaches on a protein dataset for which rigid domain fitting (RDF) models have been deposited in the PDB. In general, a consensus is observed in the conformations, which indicates a convergence from these theoretically different approaches to the most probable solution corresponding to the cryo-EM map. However, the result shows that the convergence might not be observed for proteins with complex conformational changes or with missing densities in cryo-EM map. In contrast, RDF structures deposited in the PDB can represent conformations that not only differ from the consensus obtained by flexible fitting but also from X-ray crystallography. Thus, this study emphasizes that a "consensus" achieved by the use of several automated flexible fitting approaches can provide a higher level of confidence in the modeled configurations. Following this protocol not only increases the confidence level of fitting, but also highlights protein regions with uncertain fitting. Hence, this protocol can lead to better interpretation of cryo-EM data.

Initial evaluation of a direct detection device detector for single particle cryo-electron microscopy

We report on initial results of using a new direct detection device (DDD) for single particle reconstruction of vitreous ice embedded specimens. Images were acquired on a Tecnai F20 at 200keV and a nominal magnification of 29,000×. This camera has a significantly improved signal to noise ratio and modulation transfer function (MTF) at 200keV compared to a standard CCD camera installed on the same microscope. Control of the DDD has been integrated into Leginon, an automated data collection system. Using GroEL as a test specimen, we obtained images of ∼30K particles with the CCD and the DDD from the same specimen sample using essentially identical imaging conditions. Comparison of the maps reconstructed from the CCD images and the DDD images demonstrates the improved performance of the DDD. We also obtained a 3D reconstruction from ∼70K GroEL particles acquired using the DDD; the quality of the density map demonstrates the potential of this new recording device for cryoEM data acquisition.

Analysis of peptides and proteins in their binding to GroEL

The GroEL-GroES is an essential molecular chaperon system that assists protein folding in cell. Binding of various substrate proteins to GroEL is one of the key aspects in GroEL-assisted protein folding. Small peptides may mimic segments of the substrate proteins in contact with GroEL and allow detailed structural analysis of the interactions. A model peptide SBP has been shown to bind to a region in GroEL that is important for binding of substrate proteins. Here, we investigated whether the observed GroEL-SBP interaction represented those of GroEL-substrate proteins, and whether SBP was able to mimic various aspects of substrate proteins in GroE-assisted protein folding cycle. We found that SBP competed with substrate proteins, including α-lactalbumin, rhodanese, and malate dehydrogenase, in binding to GroEL. SBP stimulated GroEL ATP hydrolysis rate in a manner similar to that of α-lactalbumin. SBP did not prevent GroES from binding to GroEL, and GroES association reduced the ATPase rates of GroEL/SBP and GroEL/α-lactalbumin to a comparable extent. Binding of both SBP and α-lactalbumin to apo GroEL was dominated by hydrophobic interaction. Interestingly, association of α-lactalbumin to GroEL/GroES was thermodynamically distinct from that to GroEL with reduced affinity and decreased contribution from hydrophobic interaction. However, SBP did not display such differential binding behaviors to apo GroEL and GroEL/GroES, likely due to the lack of a contiguous polypeptide chain that links all of the bound peptide fragments. Nevertheless, studies using peptides provide valuable information on the nature of GroEL-substrate protein interaction, which is central to understand the mechanism of GroEL-assisted protein folding.

Accurate flexible fitting of high-resolution protein structures into cryo-electron microscopy maps using coarse-grained pseudo-energy minimization

Cryo-electron microscopy (cryo-EM) has been widely used to explore conformational states of large biomolecular assemblies. The detailed interpretation of cryo-EM data requires the flexible fitting of a known high-resolution protein structure into a low-resolution cryo-EM map. To this end, we have developed what we believe is a new method based on a two-bead-per-residue protein representation, and a modified form of the elastic network model that allows large-scale conformational changes while maintaining pseudobonds and secondary structures. Our method minimizes a pseudo-energy which linearly combines various terms of the modified elastic network model energy with a cryo-EM-fitting score and a collision energy that penalizes steric collisions. Unlike previous flexible fitting efforts using the lowest few normal modes, our method effectively utilizes all normal modes so that both global and local structural changes can be fully modeled. We have validated our method for a diverse set of 10 pairs of protein structures using simulated cryo-EM maps with a range of resolutions and in the absence/presence of random noise. We have shown that our method is both accurate and efficient compared with alternative techniques, and its performance is robust to the addition of random noise. Our method is also shown to be useful for the flexible fitting of three experimental cryo-EM maps.

Cooperative dynamics of proteins unraveled by network models

Recent years have seen a significant increase in the number of computational studies that adopted network models for investigating biomolecular systems dynamics and interactions. In particular, elastic network models have proven useful in elucidating the dynamics and allosteric signaling mechanisms of proteins and their complexes. Here we present an overview of two most widely used elastic network models, the Gaussian Network Model (GNM) and Anisotropic Network Model (ANM). We illustrate their use in (i) explaining the anisotropic response of proteins observed in external pulling experiments, (ii) identifying residues that possess high allosteric potentials, and demonstrating in this context the propensity of catalytic sites and metal-binding sites for enabling efficient signal transduction, and (iii) assisting in structure refinement, molecular replacement and comparative modeling of ligand-bound forms via efficient sampling of energetically favored conformers.

Exploring the spatial and temporal organization of a cell's proteome

To increase our current understanding of cellular processes, such as cell signaling and division, knowledge is needed about the spatial and temporal organization of the proteome at different organizational levels. These levels cover a wide range of length and time scales: from the atomic structures of macromolecules for inferring their molecular function, to the quantitative description of their abundance, and spatial distribution in the cell. Emerging new experimental technologies are greatly increasing the availability of such spatial information on the molecular organization in living cells. This review addresses three fields that have significantly contributed to our understanding of the proteome's spatial and temporal organization: first, methods for the structure determination of individual macromolecular assemblies, specifically the fitting of atomic structures into density maps generated from electron microscopy techniques; second, research that visualizes the spatial distributions of these complexes within the cellular context using cryo electron tomography techniques combined with computational image processing; and third, methods for the spatial modeling of the dynamic organization of the proteome, specifically those methods for simulating reaction and diffusion of proteins and complexes in crowded intracellular fluids. The long-term goal is to integrate the varied data about a proteome's organization into a spatially explicit, predictive model of cellular processes.

Characterisation of a GroEL single-ring mutant that supports growth of Escherichia coli and has GroES-dependent ATPase activity

Binding and folding of substrate proteins by the molecular chaperone GroEL alternates between its two seven-membered rings in an ATP-regulated manner. The association of ATP and GroES to a polypeptide-bound ring of GroEL encapsulates the folding proteins in the central cavity of that ring (cis ring) and allows it to fold in a protected environment where the risk of aggregation is reduced. ATP hydrolysis in the cis ring changes the potentials within the system such that ATP binding to the opposite (trans) ring triggers the release of all ligands from the cis ring of GroEL through a complex network of allosteric communication between the rings. Inter-ring allosteric communication thus appears indispensable for the function of GroEL, and an engineered single-ring version (SR1) cannot substitute for GroEL in vivo. We describe here the isolation and characterisation of an active single-ring form of the GroEL protein (SR-A92T), which has an exceptionally low ATPase activity that is strongly stimulated by the addition of GroES. Dissection of the kinetic pathway of the ATP-induced structural changes in this active single ring can be explained by the fact that the mutation effectively blocks progression through the full allosteric pathway of the GroEL reaction cycle, thus trapping an early allosteric intermediate. Addition of GroES is able to overcome this block by binding this intermediate and pulling the allosteric pathway to completion via mass action, explaining how bacterial cells expressing this protein as their only chaperonin are viable.

Chaperonin-mediated protein folding: using a central cavity to kinetically assist polypeptide chain folding

The chaperonin ring assembly GroEL provides kinetic assistance to protein folding in the cell by binding non-native protein in the hydrophobic central cavity of an open ring and subsequently, upon binding ATP and the co-chaperonin GroES to the same ring, releasing polypeptide into a now hydrophilic encapsulated cavity where productive folding occurs in isolation. The fate of polypeptide during binding, encapsulation, and folding in the chamber has been the subject of recent experimental studies and is reviewed and considered here. We conclude that GroEL, in general, behaves passively with respect to its substrate proteins during these steps. While binding appears to be able to rescue non-native polypeptides from kinetic traps, such rescue is most likely exerted at the level of maximizing hydrophobic contact, effecting alteration of the topology of weakly structured states. Encapsulation does not appear to involve 'forced unfolding', and if anything, polypeptide topology is compacted during this step. Finally, chamber-mediated folding appears to resemble folding in solution, except that major kinetic complications of multimolecular association are prevented.

Cryo-EM structure of the native GroEL-GroES complex from thermus thermophilus encapsulating substrate inside the cavity

The chaperonin GroEL interacts with various proteins, leading them to adopt their correct conformations with the aid of GroES and ATP. The actual mechanism is still being debated. In this study, by use of cryo-electron microscopy, we determined the solution structure of the Thermus thermophilus GroEL-GroES complex encapsulating its substrate proteins. We observed the averaged density of substrate proteins in the center of the GroEL-GroES cavity. The position of the averaged substrate density in the cavity suggested a repulsive interaction between a majority of the substrate proteins and the interior wall of the cavity, which is suitable for substrate release. In addition, we observed a distortion of the cis-GroEL ring, especially at the position near the substrate, which indicated that the interaction between the encapsulated proteins and the GroEL ring results in an adjustment in the cavity's shape to accommodate the substrate.

Chaperonin complex with a newly folded protein encapsulated in the folding chamber

A subset of essential cellular proteins requires the assistance of chaperonins (in Escherichia coli, GroEL and GroES), double-ring complexes in which the two rings act alternately to bind, encapsulate and fold a wide range of nascent or stress-denatured proteins. This process starts by the trapping of a substrate protein on hydrophobic surfaces in the central cavity of a GroEL ring. Then, binding of ATP and co-chaperonin GroES to that ring ejects the non-native protein from its binding sites, through forced unfolding or other major conformational changes, and encloses it in a hydrophilic chamber for folding. ATP hydrolysis and subsequent ATP binding to the opposite ring trigger dissociation of the chamber and release of the substrate protein. The bacteriophage T4 requires its own version of GroES, gp31, which forms a taller folding chamber, to fold the major viral capsid protein gp23 (refs 16-20). Polypeptides are known to fold inside the chaperonin complex, but the conformation of an encapsulated protein has not previously been visualized. Here we present structures of gp23-chaperonin complexes, showing both the initial captured state and the final, close-to-native state with gp23 encapsulated in the folding chamber. Although the chamber is expanded, it is still barely large enough to contain the elongated gp23 monomer, explaining why the GroEL-GroES complex is not able to fold gp23 and showing how the chaperonin structure distorts to enclose a large, physiological substrate protein.

A method for the alignment of heterogeneous macromolecules from electron microscopy

We propose a feature-based image alignment method for single-particle electron microscopy that is able to accommodate various similarity scoring functions while efficiently sampling the two-dimensional transformational space. We use this image alignment method to evaluate the performance of a scoring function that is based on the Mutual Information (MI) of two images rather than one that is based on the cross-correlation function. We show that alignment using MI for the scoring function has far less model-dependent bias than is found with cross-correlation based alignment. We also demonstrate that MI improves the alignment of some types of heterogeneous data, provided that the signal-to-noise ratio is relatively high. These results indicate, therefore, that use of MI as the scoring function is well suited for the alignment of class-averages computed from single-particle images. Our method is tested on data from three model structures and one real dataset.

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.

Structural frameworks for considering microbial protein- and nucleic acid-dependent motor ATPases

Many fundamental cellular processes depend on enzymes that utilize chemical energy to catalyse unfavourable reactions. Certain classes of ATPases provide a particularly vivid example of the process of energy conversion, employing cycles of nucleotide turnover to move and/or rearrange biological polymers such as proteins and nucleic acids. Four well-characterized classes of ATP-dependent protein/nucleic acid translocases and remodelling factors are found in all three domains of life (bacteria, archaea and eukarya): additional strand catalytic 'E' (ASCE) P-loop NTPases, GHL proteins, actin-fold enzymes and chaperonins. These unrelated protein superfamilies have each evolved the ability to couple ATP binding and hydrolysis to the generation of motion and force along or within their substrates. The past several years have witnessed the emergence of a wealth of structural data that help explain how such molecular engines link nucleotide turnover to conformational change. In this review, we highlight several recent advances to illustrate some of the mechanisms by which each family of ATP-dependent motors facilitates the rearrangement and movement of proteins, protein complexes and nucleic acids.

Chaperonins: The hunt for the Group II mechanism

Chaperonins are multi-subunit complexes that enhance the efficiency of protein-folding reactions by capturing protein substrates in their central cavities. They occur in all prokaryotic and eukaryotic cell types and, alone amongst molecular chaperones, chaperonin knockouts are always lethal. Chaperonins come in two forms; the Group I are found in bacteria, mitochondria and plastids [W.A. Fenton, A.L. Horwich, Q. Rev. Biophys. 36 (2003) 229-256, [1]] and the Group II in the eukaryotic cytoplasm and in archaea [N.J. Cowan, S.A. Lewis, Adv. Protein Chem. 59 (2001) 73-104, [2]]. Both use energy derived from ATP binding and hydrolysis to drive a series of structural rearrangements that enable them to capture, engulf and then release polypeptide chains that have either not yet acquired the native, biologically active state or have been denatured in the cell.

Combining efficient conformational sampling with a deformable elastic network model facilitates structure refinement at low resolution

Structural studies of large proteins and protein assemblies are a difficult and pressing challenge in molecular biology. Experiments often yield only low-resolution or sparse data that are not sufficient to fully determine atomistic structures. We have developed a general geometry-based algorithm that efficiently samples conformational space under constraints imposed by low-resolution density maps obtained from electron microscopy or X-ray crystallography experiments. A deformable elastic network (DEN) is used to restrain the sampling to prior knowledge of an approximate structure. The DEN restraints dramatically reduce over-fitting, especially at low resolution. Cross-validation is used to optimally weight the structural information and experimental data. Our algorithm is robust even for noise-added density maps and has a large radius of convergence for our test case. The DEN restraints can also be used to enhance reciprocal space simulated annealing refinement.

Tandem mass spectrometry of intact GroEL-substrate complexes reveals substrate-specific conformational changes in the trans ring

It has been suggested that the bacterial GroEL chaperonin accommodates only one substrate at any given time, due to conformational changes to both the cis and trans ring that are induced upon substrate binding. Using electrospray ionization mass spectrometry, we show that indeed GroEL binds only one molecule of the model substrate Rubisco. In contrast, the capsid protein of bacteriophage T4, a natural GroEL substrate, can occupy both rings simultaneously. As these substrates are of similar size, the data indicate that each substrate induces distinct conformational changes in the GroEL chaperonin. The distinctive binding behavior of Rubisco and the capsid protein was further investigated using tandem mass spectrometry on the intact 800-914 kDa GroEL-substrate complexes. Our data suggest that even in the gas phase the substrates remain bound inside the GroEL cavity. The analysis revealed further that binding of Rubisco to the GroEL oligomer stabilizes the chaperonin complex significantly, whereas binding of one capsid protein did not have the same effect. However, addition of a second capsid protein molecule to GroEL resulted in a similar stabilizing effect to that obtained after the binding of a single Rubisco. On the basis of the stoichiometry of the GroEL chaperonin-substrate complex and the dissociation behavior of the two different substrates, we hypothesize that the binding of a single capsid polypeptide does not induce significant conformational changes in the GroEL trans ring, and hence the unoccupied GroEL ring remains accessible for a second capsid molecule.

3D electron microscopy of biological nanomachines: principles and applications

Transmission electron microscopy is a powerful technique for studying the three-dimensional (3D) structure of a wide range of biological specimens. Knowledge of this structure is crucial for fully understanding complex relationships among macromolecular complexes and organelles in living cells. In this paper, we present the principles and main application domains of 3D transmission electron microscopy in structural biology. Moreover, we survey current developments needed in this field, and discuss the close relationship of 3D transmission electron microscopy with other experimental techniques aimed at obtaining structural and dynamical information from the scale of whole living cells to atomic structure of macromolecular complexes.

Topologies of a substrate protein bound to the chaperonin GroEL

The chaperonin GroEL assists polypeptide folding through sequential steps of binding nonnative protein in the central cavity of an open ring, via hydrophobic surfaces of its apical domains, followed by encapsulation in a hydrophilic cavity. To examine the binding state, we have classified a large data set of GroEL binary complexes with nonnative malate dehydrogenase (MDH), imaged by cryo-electron microscopy, to sort them into homogeneous subsets. The resulting electron density maps show MDH associated in several characteristic binding topologies either deep inside the cavity or at its inlet, contacting three to four consecutive GroEL apical domains. Consistent with visualization of bound polypeptide distributed over many parts of the central cavity, disulfide crosslinking could be carried out between a cysteine in a bound substrate protein and cysteines substituted anywhere inside GroEL. Finally, substrate binding induced adjustments in GroEL itself, observed mainly as clustering together of apical domains around sites of substrate binding.

Inter-ring communication allows the GroEL chaperonin complex to distinguish between different substrates

The productive folding of substrate proteins by the GroEL complex of Escherichia coli requires the activity of both the chaperonin rings. These heptameric rings were shown to regulate the chaperonins' affinity for substrates and co-chaperonin via inter-ring communications; however, the molecular details of the interactions are not well understood. We have investigated the effect of substrate binding on inter-ring communications of the chaperonin complex, both the double-ring GroEL as well as the single-ring SR1 chaperonin in complex with four different substrates by using mass spectrometry. This approach shows that whereas SR1 is unable to distinguish between Rubisco, gp23, gp5, and MDH, GroEL shows clear differences upon binding these substrates. The most distinctive binding behavior is observed for Rubisco, which only occupies one GroEL ring. Both bacteriophage capsid proteins (gp23 and gp5) as well as MDH are able to bind to the two GroEL rings simultaneously. Our data suggest that inter-ring communication allows the chaperonin complex to differentiate between substrates. Using collision induced dissociation in the gas phase, differences between the chaperonin(substrate) complexes are observed only when both rings are present. The data indicate that the size of the substrate is an important factor that determines the degree of stabilization of the chaperonin complex.

Cross-correlation of common lines: a novel approach for single-particle reconstruction of a structure containing a flexible domain

We describe a novel approach to sorting class averages of a structure in multiple conformational states in order to generate 3D reconstructions that account for conformational variability present in the sample. The method assumes that the relative Euler angles between class averages are known, then uses a common lines approach to match any given class against a set of distinct conformations from a selected view of the structure. We show the effectiveness of the method both on model data and on an experimental dataset for which the conformational variability is limited to a defined region within the structure. During our studies of hepatitis C virus (HCV) internal ribosome entry site (IRES) interaction with the human translation initiation factor eIF3, we observed that the IRES RNA included a flexible region holding multiple conformations. While current classification methods were used to produce two-dimensional averages of the complex showing these different conformations, no method existed for relating these averages in three dimensions. Our approach overcame these limitations, giving us structural insight that was previously not possible.

Multiscale modeling of biomolecular systems: in serial and in parallel

Considerable progress has been recently achieved in the multiscale modeling of complex biological processes. Multiscale models have now investigated the structure and dynamics of lipid membranes, proteins, peptides and DNA over length and time scales ranging from the atomic to the macroscopic. Serial multiscale methods that parameterize low-resolution coarse-grained models with data from high-resolution models have studied long time or length scale phenomena that cannot be investigated with atomically detailed models. Parallel multiscale methods that directly couple high- and low-resolution models have efficiently explored slow structural transitions and the importance of long-wavelength fluctuations for biological molecules. The success of such models relies upon new theories and methods for constructing accurate multiscale bridges that transfer information between models with different resolutions.

An expanded conformation of single-ring GroEL-GroES complex encapsulates an 86 kDa substrate

Electron cryomicroscopy reveals an unprecedented conformation of the single-ring mutant of GroEL (SR398) bound to GroES in the presence of Mg-ATP. This conformation exhibits a considerable expansion of the folding cavity, with approximately 80% more volume than the X-ray structure of the equivalent cis cavity in the GroEL-GroES-(ADP)(7) complex. This expanded conformation can encapsulate an 86 kDa heterodimeric (alphabeta) assembly intermediate of mitochondrial branched-chain alpha-ketoacid dehydrogenase, the largest substrate ever observed to be cis encapsulated. The SR398-GroES-Mg-ATP complex is found to exist as a mixture of standard and expanded conformations, regardless of the absence or presence of the substrate. However, the presence of even a small substrate causes a pronounced bias toward the expanded conformation. Encapsulation of the large assembly intermediate is supported by a series of electron cryomicroscopy studies as well as the protection of both alpha and beta subunits of the substrate from tryptic digestion.

GroEL-mediated protein folding: making the impossible, possible

Protein folding is a spontaneous process that is essential for life, yet the concentrated and complex interior of a cell is an inherently hostile environment for the efficient folding of many proteins. Some proteins-constrained by sequence, topology, size, and function-simply cannot fold by themselves and are instead prone to misfolding and aggregation. This problem is so deeply entrenched that a specialized family of proteins, known as molecular chaperones, evolved to assist in protein folding. Here we examine one essential class of molecular chaperones, the large, oligomeric, and energy utilizing chaperonins or Hsp60s. The bacterial chaperonin GroEL, along with its co-chaperonin GroES, is probably the best-studied example of this family of protein-folding machine. In this review, we examine some of the general properties of proteins that do not fold well in the absence of GroEL and then consider how folding of these proteins is enhanced by GroEL and GroES. Recent experimental and theoretical studies suggest that chaperonins like GroEL and GroES employ a combination of protein isolation, unfolding, and conformational restriction to drive protein folding under conditions where it is otherwise not possible.

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

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