The unfolding action of GroEL on a protein substrate

From Microstates to Macrostates in the Conformational Dynamics of GroEL: A Single-Molecule Förster Resonance Energy Transfer Study

The chaperonin GroEL is a multisubunit molecular machine that assists in protein folding in the Escherichia coli cytosol. Past studies have shown that GroEL undergoes large allosteric conformational changes during its reaction cycle. Here, we report single-molecule Förster resonance energy transfer measurements that directly probe the conformational transitions of one subunit within GroEL and its single-ring variant under equilibrium conditions. We find that four microstates span the conformational manifold of the protein and interconvert on the submillisecond time scale. A unique set of relative populations of these microstates, termed a macrostate, is obtained by varying solution conditions, e.g., adding different nucleotides or the cochaperone GroES. Strikingly, ATP titration studies demonstrate that the partition between the apo and ATP-ligated conformational macrostates traces a sigmoidal response with a Hill coefficient similar to that obtained in bulk experiments of ATP hydrolysis. These coinciding results from bulk measurements for an entire ring and single-molecule measurements for a single subunit provide new evidence for the concerted allosteric transition of all seven subunits.

How soluble misfolded proteins bypass chaperones at the molecular level

Subpopulations of soluble, misfolded proteins can bypass chaperones within cells. The extent of this phenomenon and how it happens at the molecular level are unknown. Through a meta-analysis of the experimental literature we find that in all quantitative protein refolding studies there is always a subpopulation of soluble but misfolded protein that does not fold in the presence of one or more chaperones, and can take days or longer to do so. Thus, some misfolded subpopulations commonly bypass chaperones. Using multi-scale simulation models we observe that the misfolded structures that bypass various chaperones can do so because their structures are highly native like, leading to a situation where chaperones do not distinguish between the folded and near-native-misfolded states. More broadly, these results provide a mechanism by which long-time scale changes in protein structure and function can persist in cells because some misfolded states can bypass components of the proteostasis machinery.

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.

Atomic-Level Description of Protein Folding inside the GroEL Cavity

Chaperonins (ubiquitous facilitators of protein folding) sequester misfolded proteins within an internal cavity, thus preventing protein aggregation during the process of refolding. GroEL, a tetradecameric bacterial chaperonin, is one of the most studied chaperonins, but the role of the internal cavity in the refolding process is still unclear. It has been suggested that rather than simply isolating proteins while they refold, the GroEL cavity actively promotes protein folding. A detailed characterization of the folding dynamics and thermodynamics of protein substrates encapsulated within the cavity, however, has been difficult to obtain by experimental means, due to the system's complexity and the many steps in the folding cycle. Here, we examine the influence of the GroEL cavity on protein folding based on the results of unbiased, atomistic molecular dynamics simulations. We first verified that the computational setup, which uses a recently developed state-of-the-art force field that more accurately reproduces the aggregation propensity of unfolded states, could recapitulate the essential structural dynamics of GroEL. In these simulations, the GroEL tetradecamer was highly dynamic, transitioning among states corresponding to most of the structures that have been observed experimentally. We then simulated a small, unfolded protein both in the GroEL cavity and in bulk solution and compared the protein's folding process within these two environments. Inside the GroEL cavity, the unfolded protein interacted strongly with the disordered residues in GroEL's C-terminal tails. These interactions stabilized the protein's unfolded states relative to its compact states and increased the roughness of its folding free-energy surface, resulting in slower folding compared to the rate in solution. For larger proteins, which are more typical GroEL substrates, we speculate that these interactions may allow substrates to more quickly escape kinetic traps associated with compact, misfolded states, thereby actively promoting folding.

Chaperones convert the energy from ATP into the nonequilibrium stabilization of native proteins

During and after protein translation, molecular chaperones require ATP hydrolysis to favor the native folding of their substrates and, under stress, to avoid aggregation and revert misfolding. Why do some chaperones need ATP, and what are the consequences of the energy contributed by the ATPase cycle? Here, we used biochemical assays and physical modeling to show that the bacterial chaperones GroEL (Hsp60) and DnaK (Hsp70) both use part of the energy from ATP hydrolysis to restore the native state of their substrates, even under denaturing conditions in which the native state is thermodynamically unstable. Consistently with thermodynamics, upon exhaustion of ATP, the metastable native chaperone products spontaneously revert to their equilibrium non-native states. In the presence of ATPase chaperones, some proteins may thus behave as open ATP-driven, nonequilibrium systems whose fate is only partially determined by equilibrium thermodynamics.

Electric-Field-Induced Protein Translocation via a Conformational Transition in SecDF: An MD Study

SecDF is an important component of the Sec protein translocation machinery embedded in the bacterial membrane, which is associated with many functions, such as stabilizing other Sec translocon components within the membrane, maintaining the transmembrane (TM) potential, and facilitating the ATP-independent stage of the translocation mechanism. Related studies suggest that SecDF undergoes functionally important conformational changes that involve mainly its P1-head domain and that these changes are coupled with the proton motive force (Δp). However, there still is not a clear understanding of how SecDF functions, its exact role in the translocation machinery, and how its function is related to Δp. Here, using all-atom molecular dynamics simulations combined with umbrella sampling, we study the P1-head conformational change and how it is coupled to the proton motive force. We report potentials of mean force along a root-mean-square-distance-based reaction coordinate obtained in the presence and absence of the TM electrical potential. Our results show that the interaction of the P1 domain dipole moment with the TM electrical field considerably lowers the free-energy barrier in the direction of F-form to I-form transition.

Substrate binding accelerates the conformational transitions and substrate dissociation in multidrug efflux transporter AcrB

The tripartite efflux pump assembly AcrAB-TolC is the major multidrug resistance transporter in E. coli. The inner membrane transporter AcrB is a homotrimer, energized by the proton movement down the transmembrane electrochemical gradient. The asymmetric crystal structures of AcrB with three monomers in distinct conformational states [access (A), binding (B) and extrusion (E)] support a functional rotating mechanism, in which each monomer of AcrB cycles among the three states in a concerted way. However, the relationship between the conformational changes during functional rotation and drug translocation has not been totally understood. Here, we explored the conformational changes of the AcrB homotrimer during the ABE to BEA transition in different substrate-binding states using targeted MD simulations. It was found that the dissociation of substrate from the distal binding pocket of B monomer is closely related to the concerted conformational changes in the translocation pathway, especially the side chain reorientation of Phe628 and Tyr327. A second substrate binding at the proximal binding pocket of A monomer evidently accelerates the conformational transitions as well as substrate dissociation in B monomer. The acceleration effect of the multi-substrate binding mode provides a molecular explanation for the positive cooperativity observed in the kinetic studies of substrate efflux and deepens our understanding of the functional rotating mechanism of AcrB.

Parallel Cascade Selection Molecular Dynamics (PaCS-MD) to generate conformational transition pathway

Parallel Cascade Selection Molecular Dynamics (PaCS-MD) is proposed as a molecular simulation method to generate conformational transition pathway under the condition that a set of "reactant" and "product" structures is known a priori. In PaCS-MD, the cycle of short multiple independent molecular dynamics simulations and selection of the structures close to the product structure for the next cycle are repeated until the simulated structures move sufficiently close to the product. Folding of 10-residue mini-protein chignolin from the extended to native structures and open-close conformational transition of T4 lysozyme were investigated by PaCS-MD. In both cases, tens of cycles of 100-ps MD were sufficient to reach the product structures, indicating the efficient generation of conformational transition pathway in PaCS-MD with a series of conventional MD without additional external biases. Using the snapshots along the pathway as the initial coordinates, free energy landscapes were calculated by the combination with multiple independent umbrella samplings to statistically elucidate the conformational transition pathways.

Proteomic data from human cell cultures refine mechanisms of chaperone-mediated protein homeostasis

In the crowded environment of human cells, folding of nascent polypeptides and refolding of stress-unfolded proteins is error prone. Accumulation of cytotoxic misfolded and aggregated species may cause cell death, tissue loss, degenerative conformational diseases, and aging. Nevertheless, young cells effectively express a network of molecular chaperones and folding enzymes, termed here "the chaperome," which can prevent formation of potentially harmful misfolded protein conformers and use the energy of adenosine triphosphate (ATP) to rehabilitate already formed toxic aggregates into native functional proteins. In an attempt to extend knowledge of chaperome mechanisms in cellular proteostasis, we performed a meta-analysis of human chaperome using high-throughput proteomic data from 11 immortalized human cell lines. Chaperome polypeptides were about 10% of total protein mass of human cells, half of which were Hsp90s and Hsp70s. Knowledge of cellular concentrations and ratios among chaperome polypeptides provided a novel basis to understand mechanisms by which the Hsp60, Hsp70, Hsp90, and small heat shock proteins (HSPs), in collaboration with cochaperones and folding enzymes, assist de novo protein folding, import polypeptides into organelles, unfold stress-destabilized toxic conformers, and control the conformal activity of native proteins in the crowded environment of the cell. Proteomic data also provided means to distinguish between stable components of chaperone core machineries and dynamic regulatory cochaperones.

Analysis and elimination of a bias in targeted molecular dynamics simulations of conformational transitions: application to calmodulin

The popular targeted molecular dynamics (TMD) method for generating transition paths in complex biomolecular systems is revisited. In a typical TMD transition path, the large-scale changes occur early and the small-scale changes tend to occur later. As a result, the order of events in the computed paths depends on the direction in which the simulations are performed. To identify the origin of this bias, and to propose a method in which the bias is absent, variants of TMD in the restraint formulation are introduced and applied to the complex open ↔ closed transition in the protein calmodulin. Due to the global best-fit rotation that is typically part of the TMD method, the simulated system is guided implicitly along the lowest-frequency normal modes, until the large spatial scales associated with these modes are near the target conformation. The remaining portion of the transition is described progressively by higher-frequency modes, which correspond to smaller-scale rearrangements. A straightforward modification of TMD that avoids the global best-fit rotation is the locally restrained TMD (LRTMD) method, in which the biasing potential is constructed from a number of TMD potentials, each acting on a small connected portion of the protein sequence. With a uniform distribution of these elements, transition paths that lack the length-scale bias are obtained. Trajectories generated by steered MD in dihedral angle space (DSMD), a method that avoids best-fit rotations altogether, also lack the length-scale bias. To examine the importance of the paths generated by TMD, LRTMD, and DSMD in the actual transition, we use the finite-temperature string method to compute the free energy profile associated with a transition tube around a path generated by each algorithm. The free energy barriers associated with the paths are comparable, suggesting that transitions can occur along each route with similar probabilities. This result indicates that a broad ensemble of paths needs to be calculated to obtain a full description of conformational changes in biomolecules. The breadth of the contributing ensemble suggests that energetic barriers for conformational transitions in proteins are offset by entropic contributions that arise from a large number of possible paths.

A dynamic model of long-range conformational adaptations triggered by nucleotide binding in GroEL-GroES

The molecular chaperone, GroEL, essential for correct protein folding in E. coli, is composed of 14 identical subunits organized in two interacting rings, each providing a folding chamber for non-native substrate proteins. The oligomeric assembly shows positive cooperativity within each ring and negative cooperativity between the rings. Although it is well known that ATP and long-range allosteric interactions drive the functional cycle of GroEL, an atomic resolution view of how ligand binding modulates conformational adaptations over long distances remains a major challenge. Moreover, little is known on the relation between equilibrium dynamics at physiological temperatures and the allosteric transitions in GroEL. Here we present multiple all-atom molecular dynamics simulations of the GroEL-GroES assemblies at different stages of the functional cycle. Combined with an extensive analysis of the complete set of experimentally available structures, principal component analysis and conformer plots, we provide an explicit evaluation of the accessible conformational space of unliganded GroEL. Our results suggest the presence of pre-existing conformers at the equatorial domain level, and a shift of the conformational ensemble upon ATP-binding. At the inter-ring interface the simulations capture a remarkable offset motion of helix D triggered by ATP-binding to the folding active ring. The reorientation of helix D, previously only observed upon GroES association, correlates with a change of the internal dynamics in the opposite ring. This work contributes to the understanding of the molecular mechanisms in GroEL and highlights the ability of all-atom MD simulations to model long-range structural changes and allosteric events in large systems.

Unfolding and translocation pathway of substrate protein controlled by structure in repetitive allosteric cycles of the ClpY ATPase

Clp ATPases are ring-shaped AAA+ motors in the degradation pathway that perform critical actions of unfolding and translocating substrate proteins (SPs) through narrow pores to deliver them to peptidase components. These actions are effected by conserved diaphragm-forming loops found in the central channel of the Clp ATPase hexamer. Conformational changes, that take place in the course of repetitive ATP-driven cycles, result in mechanical forces applied by the central channel loops onto the SP. We use coarse-grained simulations to elucidate allostery-driven mechanisms of unfolding and translocation of a tagged four-helix bundle protein by the ClpY ATPase. Unfolding is initiated at the tagged C-terminal region via an obligatory intermediate. The resulting nonnative conformation is competent for translocation, which proceeds on a different time scale than unfolding and involves sharp stepped transitions. Completion of the translocation process requires assistance from the ClpQ peptidase. These mechanisms contrast nonallosteric mechanical unfolding of the SP. In atomic force microscopy experiments, multiple unfolding pathways are available and large mechanical forces are required to unravel the SP relative to those exerted by the central channel loops of ClpY. SP threading through a nonallosteric ClpY nanopore involves simultaneous unfolding and translocation effected by strong pulling forces.

Predicting order of conformational changes during protein conformational transitions using an interpolated elastic network model

The decryption of sequence of structural events during protein conformational transitions is essential to a detailed understanding of molecular functions of various biological nanomachines. Coarse-grained models have proven useful by allowing highly efficient simulations of protein conformational dynamics. By combining two coarse-grained elastic network models constructed based on the beginning and end conformations of a transition, we have developed an interpolated elastic network model to generate a transition pathway between the two protein conformations. For validation, we have predicted the order of local and global conformational changes during key ATP-driven transitions in three important biological nanomachines (myosin, F(1) ATPase and chaperonin GroEL). We have found that the local conformational change associated with the closing of active site precedes the global conformational change leading to mechanical motions. Our finding is in good agreement with the distribution of intermediate experimental structures, and it supports the importance of local motions at active site to drive or gate various conformational transitions underlying the workings of a diverse range of biological nanomachines.

Equilibrium sampling for biomolecules under mechanical tension

In the studies of force-induced conformational transitions of biomolecules, the large timescale difference from experiments presents the challenge of obtaining convergent sampling for molecular dynamics simulations. To circumvent this fundamental problem, an approach combining the replica-exchange method and umbrella sampling (REM-US) was developed to simulate mechanical stretching of biomolecules under equilibrium conditions. Equilibrium properties of conformational transitions can be obtained directly from simulations without further assumptions. To test the performance, we carried out REM-US simulations of atomic force microscope (AFM) stretching and relaxing measurements on the polysaccharide pustulan, a (1-->6)-beta-D-glucan, which undergoes well-characterized rotameric transitions in the backbone bonds. With significantly enhanced sampling convergence and efficiency, the REM-US approach closely reproduced the equilibrium force-extension curves measured in AFM experiments. Consistent with the reversibility in the AFM measurements, the new approach generated identical force-extension curves in both stretching and relaxing simulations-an outcome not reported in previous studies, proving that equilibrium conditions were achieved in the simulations. REM-US may provide a robust approach to modeling of mechanical stretching on polysaccharides and even nucleic acids.

The conformational transition pathway of ATP binding cassette transporter MsbA revealed by atomistic simulations

ATP binding cassette transporters are integral membrane proteins that use the energy released from ATP hydrolysis at the two nucleotide binding domains (NBDs) to translocate a wide variety of substrates through a channel at the two transmembrane domains (TMDs) across the cell membranes. MsbA from Gram-negative bacteria is a lipid and multidrug resistance ATP binding cassette exporter that can undergo large scale conformational changes between the outward-facing and the inward-facing conformations revealed by crystal structures in different states. Here, we use targeted molecular dynamics simulation methods to explore the atomic details of the conformational transition from the outward-facing to the inward-facing states of MsbA. The molecular dynamics trajectories revealed a clear spatiotemporal order of the conformational movements. The disruption of the nucleotide binding sites at the NBD dimer interface is the very first event that initiates the following conformational changes, verifying the assumption that the conformational conversion is triggered by ATP hydrolysis. The conserved x-loops of the NBDs were identified to participate in the interaction network that stabilizes the cytoplasmic tetrahelix bundle of the TMDs and play an important role in mediating the cross-talk between the NBD and TMD. The movement of the NBD dimer is transmitted through x-loops to break the tetrahelix bundle, inducing the packing rearrangements of the transmembrane helices at the cytoplasmic side and the periplasmic side sequentially. The packing rearrangement within each periplasmic wing of TMD that results in exposure of the substrate binding sites occurred at the end stage of the trajectory, preventing the wrong timing of the binding site accessibility.

CHARMM: the biomolecular simulation program

CHARMM (Chemistry at HARvard Molecular Mechanics) is a highly versatile and widely used molecular simulation program. It has been developed over the last three decades with a primary focus on molecules of biological interest, including proteins, peptides, lipids, nucleic acids, carbohydrates, and small molecule ligands, as they occur in solution, crystals, and membrane environments. For the study of such systems, the program provides a large suite of computational tools that include numerous conformational and path sampling methods, free energy estimators, molecular minimization, dynamics, and analysis techniques, and model-building capabilities. The CHARMM program is applicable to problems involving a much broader class of many-particle systems. Calculations with CHARMM can be performed using a number of different energy functions and models, from mixed quantum mechanical-molecular mechanical force fields, to all-atom classical potential energy functions with explicit solvent and various boundary conditions, to implicit solvent and membrane models. The program has been ported to numerous platforms in both serial and parallel architectures. This article provides an overview of the program as it exists today with an emphasis on developments since the publication of the original CHARMM article in 1983.

Allosteric transitions of supramolecular systems explored by network models: application to chaperonin GroEL

Identification of pathways involved in the structural transitions of biomolecular systems is often complicated by the transient nature of the conformations visited across energy barriers and the multiplicity of paths accessible in the multidimensional energy landscape. This task becomes even more challenging in exploring molecular systems on the order of megadaltons. Coarse-grained models that lend themselves to analytical solutions appear to be the only possible means of approaching such cases. Motivated by the utility of elastic network models for describing the collective dynamics of biomolecular systems and by the growing theoretical and experimental evidence in support of the intrinsic accessibility of functional substates, we introduce a new method, adaptive anisotropic network model (aANM), for exploring functional transitions. Application to bacterial chaperonin GroEL and comparisons with experimental data, results from action minimization algorithm, and previous simulations support the utility of aANM as a computationally efficient, yet physically plausible, tool for unraveling potential transition pathways sampled by large complexes/assemblies. An important outcome is the assessment of the critical inter-residue interactions formed/broken near the transition state(s), most of which involve conserved residues.

Most proteins are biomolecular machines. They perform their function by undergoing changes between different structures. Understanding the mechanism of transition between these structures is of major importance to design methods for controlling such transitions, and thereby modulating protein function. Although there are many computational methods for exploring the transitions of small proteins, the task of exploring the transition pathways becomes prohibitively expensive in the case of supramolecular systems. The bacterial chaperonin GroEL is such a supramolecular machine. It plays an important role in assisting protein folding. During its function, GroEL undergoes structural transitions between multiple forms. Here, we are introducing a new methodology, based on elastic network models, for elucidating the transition mechanisms in such supramolecular systems. Application to GroEL provides us with biologically significant information on critical interactions and sequence of events occurring during the chaperonin machinery and key contacts that make and break at the transition. The method can be readily applied to other supramolecular machines.

Polar or apolar--the role of polarity for urea-induced protein denaturation

Urea-induced protein denaturation is widely used to study protein folding and stability; however, the molecular mechanism and driving forces of this process are not yet fully understood. In particular, it is unclear whether either hydrophobic or polar interactions between urea molecules and residues at the protein surface drive denaturation. To address this question, here, many molecular dynamics simulations totalling ca. 7 µs of the CI2 protein in aqueous solution served to perform a computational thought experiment, in which we varied the polarity of urea. For apolar driving forces, hypopolar urea should show increased denaturation power; for polar driving forces, hyperpolar urea should be the stronger denaturant. Indeed, protein unfolding was observed in all simulations with decreased urea polarity. Hyperpolar urea, in contrast, turned out to stabilize the native state. Moreover, the differential interaction preferences between urea and the 20 amino acids turned out to be enhanced for hypopolar urea and suppressed (or even inverted) for hyperpolar urea. These results strongly suggest that apolar urea–protein interactions, and not polar interactions, are the dominant driving force for denaturation. Further, the observed interactions provide a detailed picture of the underlying molecular driving forces. Our simulations finally allowed characterization of CI2 unfolding pathways. Unfolding proceeds sequentially with alternating loss of secondary or tertiary structure. After the transition state, unfolding pathways show large structural heterogeneity.

To perform their physiological function, proteins have to fold into their characteristic three-dimensional structure. While the folded state is stable under physiological conditions, changes in the solvent can destabilize the folded state and even induce denaturation. One of the most commonly used denaturants is urea. Despite its widespread use to study protein folding and stability, however, the molecular mechanism and particularly the driving forces of urea-induced protein denaturation are not yet understood. Two mechanisms have been suggested, according to which denaturation is driven either by polar interactions via hydrogen bonds or by hydrophobic interactions with apolar amino acids. By systematically varying urea polarity and quantifying the interactions of the solvent molecules with all amino acids of the protein, the present simulation study reveals that it is mainly the apolar interactions that drive denaturation. Our results suggest a coherent microscopic picture for urea-induced denaturation and bear more general implications for protein stability in other environments, e.g., in chaperone-assisted folding.

Allostery and cooperativity revisited

Although phenomenlogical models that account for cooperativity in allosteric systems date back to the early and mid-60's (e.g., the KNF and MWC models), there is resurgent interest in the topic due to the recent experimental and computational studies that attempted to reveal, at an atomistic level, how allostery actually works. In this review, using systems for which atomistic simulations have been carried out in our groups as examples, we describe the current understanding of allostery, how the mechanisms go beyond the classical MWC/Pauling-KNF descriptions, and point out that the "new view" of allostery, emphasizing "population shifts," is, in fact, an "old view." The presentation offers not only an up-to-date description of allostery from a theoretical/computational perspective, but also helps to resolve several outstanding issues concerning allostery.

Proteins: coexistence of stability and flexibility

We introduce an equation for protein native topology based on recent analysis of data from the Protein Data Bank and on a generalization of the Landau-Peierls instability criterion for fractals. The equation relates the protein fractal dimension df, the spectral dimension ds, and the number of amino acids N. Deviations from the equation may render a protein unfolded. The fractal nature of proteins is shown to bridge their seemingly conflicting properties of stability and flexibility. Over 500 proteins have been analyzed (df, ds, and N) and found to obey this equation of state.

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.

Allosteric transitions in the chaperonin GroEL are captured by a dominant normal mode that is most robust to sequence variations

The Escherichia coli chaperonin GroEL, which helps proteins to fold, consists of two heptameric rings stacked back-to-back. During the reaction cycle GroEL undergoes a series of allosteric transitions triggered by ligand (substrate protein, ATP, and the cochaperonin GroES) binding. Based on an elastic network model of the bullet-shaped double-ring chaperonin GroEL-(ADP)(7)-GroES structure (R''T state), we perform a normal mode analysis to explore the energetically favorable collective motions encoded in the R''T structure. By comparing each normal mode with the observed conformational changes in the R''T --> TR'' transition, a single dominant normal mode provides a simple description of this highly intricate allosteric transition. A detailed analysis of this relatively high-frequency mode describes the structural and dynamic changes that underlie the positive intra-ring and negative inter-ring cooperativity. The dynamics embedded in the dominant mode entails highly concerted structural motions with approximate preservation of sevenfold symmetry within each ring and negatively correlated ones between the two rings. The dominant normal mode (in comparison with the other modes) is robust to parametric perturbations caused by sequence variations, which validates its functional importance. Response of the dominant mode to local changes that mimic mutations using the structural perturbation method technique leads to a wiring diagram that identifies a network of key residues that regulate the allosteric transitions. Many of these residues are located in intersubunit interfaces, and may therefore play a critical role in transmitting allosteric signals between subunits.

Coupling between allosteric transitions in GroEL and assisted folding of a substrate protein

Escherichia coli chaperonin, GroEL, helps proteins fold under nonpermissive conditions. During the reaction cycle, GroEL undergoes allosteric transitions in response to binding of a substrate protein (SP), ATP, and the cochaperonin GroES. Using coarse-grained representations of the GroEL and GroES structures, we explore the link between allosteric transitions and the folding of a model SP, a de novo-designed four-helix bundle protein, with low spontaneous yield. The ensemble of GroEL-bound SP is less structured than the bulk misfolded structures. Upon binding, which kinetically occurs in two stages, the SP loses not only native tertiary contacts but also experiences a decrease in helical content. During multivalent binding and the subsequent ATP-driven transition of GroEL the SP undergoes force-induced stretching. Upon encapsulation, which occurs upon GroES binding, the SP finds itself in a "hydrophilic" cavity in which it can reach the folded conformation. Surprisingly, we find that the yield of the native state in the expanded GroEL cavity is relatively small even after it remains in it for twice the spontaneous folding time. Thus, in accord with the iterative annealing mechanism, multiple rounds of binding, partial unfolding, and release of the SP are required to enhance the yield of the folded SP.

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.

Spinach on the ceiling: a theoretical chemist's return to biology

I was born in Vienna and came to the United States as a refugee in October 1938. This experience played an important role in my view of the world and my approach to science: It contributed to my realization that it was safe to stop working in fields that I felt I understood and to focus on different areas of research by asking questions that would teach me and others something new. I describe my experiences that led me from chemistry and physics back to my first love, biology, and outline some of the contributions I have made as part of my ongoing learning experience.

Simulation of conformational transitions by the restricted perturbation-targeted molecular dynamics method

A method for the simulation of conformational transitions is presented. The method, based on targeted molecular dynamics (TMD), limits the conformational change at each molecular dynamics step to a fixed size, that minimizes the root mean square deviation from the target. The method is more efficient than standard TMD and yields lower energy pathways, but, like the TMD method, requires only a single molecular dynamics simulation. Test calculations and comparisons with standard TMD calculations for the alanine dipeptide with the analytic continuum electrostatics implicit solvent model are presented.

Glu257 in GroEL is a sensor involved in coupling polypeptide substrate binding to stimulation of ATP hydrolysis

The ATPase activity of many types of molecular chaperones is stimulated by polypeptide substrate binding via molecular mechanisms that are, for the most part, unknown. Here, we report that such stimulation of the ATPase activity of GroEL is abolished when its conserved apical domain residue Glu257 is replaced by alanine. This mutation is also found to convert the ATPase profile of GroEL, a group I chaperonin, into one that is characteristic of group II chaperonins. Steady-state and transient kinetic analysis indicate that both effects are due, at least in part, to a reduction of the affinity of GroEL for ADP. This finding indicates that nonfolded proteins stimulate ATP hydrolysis by accelerating the off-rate of the ADP formed, thereby allowing more rapid cycles of ATP binding and hydrolysis.

Folding behavior of chaperonin-mediated substrate protein

Chaperonin-mediated protein folding is complex. There have been diverse results on folding behavior, and the chaperonin molecules have been investigated as enhancing or retarding the folding rate. To understand the diversity of chaperonin-mediated protein folding, we report a study based on simulations using a simplified Gō-type model. By considering effects of affinity between the substrate protein and the chaperonin wall and spatial confinement of the chaperonin cavity, we study the thermodynamics and kinetics of folding of an unfrustrated substrate protein encapsulated in a chaperonin cavity. The affinity makes the hydrophobic residues of the protein bind to the chaperonin wall, and a strong (or weak) affinity results in a large (or small) effect of binding. Compared with the folding in bulk, the folding in chaperonin cavity with different strengths of affinity shows two kinds of behaviors: one with less dependence on the affinity but more reliance on the spatial confinement effect and the other relying strongly on the affinity. It is found that the enhancement or retardation of the folding rate depends on the competition between the spatial confinement and the affinity due to the chaperonin cavity, and a strong affinity produces a slow folding while a weak affinity induces a fast folding. The crossover between two kinds of folding behaviors happens in the case that the favorable effect of confinement is balanced by the unfavorable effect of the affinity, and a critical affinity strength is roughly defined. By analyzing the contacts formed between the residues of the protein and the chaperonin wall and between the residues of the protein themselves, the role of the affinity in the folding processes is studied. The binding of the residues with the chaperonin wall reduces the formation of both native contacts and nonnative contact or mis-contacts, providing a loose structure for further folding after allosteric change of the chaperonin cavity. In addition, 15 single-site-mutated mutants are simulated in order to test the validity of our model and to investigate the importance of affinity. Inspiringly, our results of the folding rates have a good correlation with those obtained from experiments. The folding rates are inversely correlated with the strength of the binding interactions, i.e., the weaker the binding, the faster the folding. We also find that the inner hydrophobic residues have larger effects on the folding kinetics than those of the exterior hydrophobic residues. We suggest that, besides the confinement effect, the affinity acts as another important factor to affect the folding of the substrate proteins in chaperonin systems, providing an understanding of the folding mechanism of the molecular chaperonin systems.

Translocation boost protein-folding efficiency of double-barreled chaperonins

Incorrect folding of proteins in living cells may lead to malfunctioning of the cell machinery. To prevent such cellular disasters from happening, all cells contain molecular chaperones that assist nonnative proteins in folding into the correct native structure. One of the most studied chaperone complexes is the GroEL-GroES complex. The GroEL part has a "double-barrel" structure, which consists of two cylindrical chambers joined at the bottom in a symmetrical fashion. The hydrophobic rim of one of the GroEL chambers captures nonnative proteins. The GroES part acts as a lid that temporarily closes the filled chamber during the folding process. Several capture-folding-release cycles are required before the nonnative protein reaches its native state. Here we report molecular simulations that suggest that translocation of the nonnative protein through the equatorial plane of the complex boosts the efficiency of the chaperonin action. If the target protein is correctly folded after translocation, it is released. However, if it is still nonnative, it is likely to remain trapped in the second chamber, which then closes to start a reverse translocation process. This shuttling back and forth continues until the protein is correctly folded. Our model provides a natural explanation for the prevalence of double-barreled chaperonins. Moreover, we argue that internal folding is both more efficient and safer than a scenario where partially refolded proteins escape from the complex before being recaptured.

Probing the "annealing" mechanism of GroEL minichaperone using molecular dynamics simulations

Although the intact chaperonin machinery is needed to rescue natural substrate proteins (SPs) under non-permissive conditions the "minichaperone" alone, containing only the isolated apical domain of GroEL, can assist folding of a certain class of proteins. To understand the annealing function of the minichaperone, we have carried out molecular dynamics simulations in the NPT ensemble totaling 300ns for four systems; namely, the isolated strongly binding peptide (SBP), the minichaperone, and the SBP and a weakly binding peptide (WBP) in complex with the minichaperone. The SBP, which is structureless in isolation, adopts a beta-hairpin conformation in complex with the minichaperone suggesting that favorable non-specific interactions of the SPs confined to helices H and I of the apical domains can induce local secondary structures. Comparison of the dynamical fluctuations of the apo and the liganded forms of the minichaperone shows that the stability (needed for SP capture) involves favorable hydrophobic interactions and hydrogen bond network formation between the SBP and WBP, and helices H and I. The release of the SP, which is required for the annealing action, involves water-mediated interactions of the charged residues at the ends of H and I helices. The simulation results are consistent with a transient binding release (TBR) model for the annealing action of the minichaperone. According to the TBR model, SP annealing occurs in two stages. In the first stage the SP is captured by the apical domain. This is followed by SP release (by thermal fluctuations) that places it in a different region of the energy landscape from which it can partition rapidly to the native state with probability Phi or be trapped in another misfolded state. The process of binding and release can result in enhancement of the native state yield. The TBR model suggests "that any cofactor that can repeatedly bind and release SPs can be effective in assisting protein folding." By comparing the structures of the non-chaperone alpha-casein (which has no sequence similarity with the apical domain) and the minichaperone and the hydrophobicity profiles we show that alpha-casein has a pair of helices that have similar sequence and structural profiles as H and I. Based on this comparison we identify residues that stabilize (destabilize) alpha-casein-protein complexes. This suggests that alpha-casein assists folding by the TBR mechanism.

Protein structural transitions and their functional role

Living cells are a collection of molecular machines which carry out many of the functions essential for the cell's existence, differentiation and reproduction. Most, though not all, of these machines are made up of proteins. Because of their complexity, an understanding of how they work requires a synergistic combination of experimental and theoretical studies. In this paper we outline our studies of two such protein machines. One is GroEL, the chaperone from Escherichia coli, which aids in protein folding; the other is F(1)-ATPase, a motor protein which synthesizes and hydrolyses ATP.