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

Energy landscape remodeling mechanism of Hsp70-chaperone-accelerated protein folding

Hsp70 chaperone is one of the key protein machines responsible for the quality control of protein production in cells. Facilitating in vivo protein folding by counteracting misfolding and aggregation is the essence of its biological function. Although the allosteric cycle during its functional actions has been well characterized both experimentally and computationally, the mechanism by which Hsp70 assists protein folding is still not fully understood. In this work, we studied the Hsp70-mediated folding of model proteins with rugged energy landscape by using molecular simulations. Different from the canonical scenario of Hsp70 functioning, which assumes that folding of substrate proteins occurs spontaneously after releasing from chaperones, our results showed that the substrate protein remains in contacts with the chaperone during its folding process. The direct chaperone-substrate interactions in the open conformation of Hsp70 tend to shield the substrate sites prone to form non-native contacts, which therefore avoids the frustrated folding pathway, leading to a higher folding rate and less probability of misfolding. Our results suggest that in addition to the unfoldase and holdase functions widely addressed in previous studies, Hsp70 can facilitate the folding of its substrate proteins by remodeling the folding energy landscape and directing the folding processes, demonstrating the foldase scenario. These findings add new, to our knowledge, insights into the general molecular mechanisms of chaperone-mediated protein folding.

Signalling networks and dynamics of allosteric transitions in bacterial chaperonin GroEL: implications for iterative annealing of misfolded proteins

Signal transmission at the molecular level in many biological complexes occurs through allosteric transitions. Allostery describes the responses of a complex to binding of ligands at sites that are spatially well separated from the binding region. We describe the structural perturbation method, based on phonon propagation in solids, which can be used to determine the signal-transmitting allostery wiring diagram (AWD) in large but finite-sized biological complexes. Application to the bacterial chaperonin GroEL-GroES complex shows that the AWD determined from structures also drives the allosteric transitions dynamically. From both a structural and dynamical perspective these transitions are largely determined by formation and rupture of salt-bridges. The molecular description of allostery in GroEL provides insights into its function, which is quantitatively described by the iterative annealing mechanism. Remarkably, in this complex molecular machine, a deep connection is established between the structures, reaction cycle during which GroEL undergoes a sequence of allosteric transitions, and function, in a self-consistent manner.This article is part of a discussion meeting issue 'Allostery and molecular machines'.

Minichaperone (GroEL191-345) mediated folding of MalZ proceeds by binding and release of native and functional intermediates

The isolated apical domain of GroEL consisting of residues 191-345 (known as "minichaperone") binds and assists the folding of a wide variety of client proteins without GroES and ATP, but the mechanism of its action is still unknown. In order to probe into the matter, we have examined minichaperone-mediated folding of a large aggregation prone protein Maltodextrin-glucosidase (MalZ). The key objective was to identify whether MalZ exists free in solution, or remains bound to, or cycling on and off the minichaperone during the refolding process. When GroES was introduced during refolding process, production of the native MalZ was inhibited. We also observed the same findings with a trap mutant of GroEL, which stably captures a predominantly non-native MalZ released from minichaperone during refolding process, but does not release it. Tryptophan and ANS fluorescence measurements indicated that refolded MalZ has the same structure as the native MalZ, but that its structure when bound to minichaperone is different. Surface plasmon resonance measurements provide an estimate for the equilibrium dissociation constant KD for the MalZ-minichaperone complex of 0.21 ± 0.04 μM, which are significantly higher than for most GroEL clients. This showed that minichaperone interacts loosely with MalZ to allow the protein to change its conformation and fold while bound during the refolding process. These observations suggest that the minichaperone works by carrying out repeated cycles of binding aggregation-prone protein MalZ in a relatively compact conformation and in a partially folded but active state, and releasing them to attempt to fold in solution.

Mechanism of transient binding and release of substrate protein during the allosteric cycle of the p97 nanomachine

ATPases associated with various cellular activities (AAA+) form a superfamily of ring-shaped motor proteins that utilize cyclical allosteric motions to remodel or translocate substrate proteins (SP) through a narrow central pore. The p97 ATPase is a homohexameric, double-ring member of this superfamily that encloses a central channel with nonuniform width. A narrow compartment is present within the D1 ring and a larger cavity within the D2 ring, separated by a constriction formed by six His amino acids. We use molecular dynamics simulations to probe the interaction between p97 and an extended peptide substrate. Mechanical pulling of the substrate through the p97 pore reveals that smaller work is required for translocation from the D1 toward the D2 compartment than in the opposite direction. These distinct energetic requirements originate in structural aspects and chemical properties of the pore lining. Whereas van der Waals interactions are dominant within the D1 pore, interaction within the D2 pore are strongly electrostatic. Two charged amino acids in the D2 pore, Arg599 and Glu554, provide the largest contribution to the interaction and hinder translocation from the D2 pore. SP threading requires smaller forces when the SP is pulled from the D1 side due to lower barrier to rotation of the His side chains in the direction of the D2 pore. Based on additional simulations of SP binding to two allosteric conformations of p97, we propose that transient binding and release of SP from the pore involves a lever mechanism. Binding to the open pore conformation of p97 occurs primarily at the Arg599 side chain, where the SP backbone is engaged through electrostatic interactions and hydrogen bonds. ATP-driven conformational transitions within the D2 ring alter the chemical environment inside the p97 cavity in the closed pore state. In this state, Glu554 side chains project further into the pore and interacts strongly through van der Waals contacts with the SP backbone. Based on mutations at the two sites in each of the states we identify a specific requirement of these side chains for interaction with the substrate.

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.

Versatile substrate protein recognition mechanism of the eukaryotic chaperonin CCT

Group II chaperonins, found in eukaryotic and archaeal organisms, recognize substrate proteins through diverse mechanisms that involve either hydrophobic- or electrostatic-dominated interactions. This action is distinct from the universal substrate recognition mechanism of group I chaperonins, which bind a wide spectrum of non-native proteins primarily through hydrophobic interactions. We use computational approaches to pinpoint the substrate protein binding sites of the gamma-subunit of the eukaryotic chaperonin CCT and to identify its interactions with the stringent substrate beta-tubulin. Protein-protein docking methods reveal intrinsic binding sites of CCT comprising a helical (HL) region, homologous to the GroEL-binding site, and the helical protrusion (HP) region. We performed molecular dynamics simulations of the solvated CCTgamma apical domain, beta-tubulin peptide-CCTgamma complexes, and isolated beta-tubulin peptides. We find that tubulin binds to CCTgamma through an extensive interface that spans both the HL region and the HP region. HL interactions involve both hydrophobic and electrostatic contacts, while binding to the HP region is stabilized almost exclusively by a salt bridge network. On the basis of additional simulations of a beta-tubulin-CCTgamma complex that involves a reduced interface, centered onto the HP region, we conclude that this salt bridge network is the minimal stabilizing interaction required. Strong conservation of the charged amino acids that participate in the salt bridge network, Arg306 and Glu271, indicates a general mechanism across the nonidentical CCT subunits and group II chaperonins.

Reconciling theories of chaperonin accelerated folding with experimental evidence

For the last 20 years, a large volume of experimental and theoretical work has been undertaken to understand how chaperones like GroEL can assist protein folding in the cell. The most accepted explanation appears to be the simplest: GroEL, like most other chaperones, helps proteins fold by preventing aggregation. However, evidence suggests that, under some conditions, GroEL can play a more active role by accelerating protein folding. A large number of models have been proposed to explain how this could occur. Focused experiments have been designed and carried out using different protein substrates with conclusions that support many different mechanisms. In the current article, we attempt to see the forest through the trees. We review all suggested mechanisms for chaperonin-mediated folding and weigh the plausibility of each in light of what we now know about the most stringent, essential, GroEL-dependent protein substrates.

An unusual intrinsically disordered protein from the model legume Lotus japonicus stabilizes proteins in vitro

Intrinsic structural disorder is a prevalent feature of proteins with chaperone activity. Using a complementary set of techniques, we have structurally characterized LjIDP1 (intrinsically disordered protein 1) from the model legume Lotus japonicus, and our results provide the first structural characterization of a member of the Lea5 protein family (PF03242). Contrary to in silico predictions, we show that LjIDP1 is intrinsically disordered and probably exists as an ensemble of conformations with limited residual beta-sheet, turn/loop, and polyproline II secondary structure. Furthermore, we show that LjIDP1 has an inherent propensity to undergo a large conformational shift, adopting a largely alpha-helical structure when it is dehydrated and in the presence of different detergents and alcohols. This is consistent with an overrepresentation of order-promoting residues in LjIDP1 compared with the average of intrinsically disordered proteins. In line with functioning as a chaperone, we show that LjIDP1 effectively prevents inactivation of two model enzymes under conditions that promote protein misfolding and aggregation. The LjIdp1 gene is expressed in all L. japonicus tissues tested. A higher expression level was found in the root tip proximal zone, in roots inoculated with compatible endosymbiotic M. loti, and in functional nitrogen-fixing root nodules. We suggest that the ability of LjIDP1 to prevent protein misfolding and aggregation may play a significant role in tissues, such as symbiotic root nodules, which are characterized by high metabolic activity.

Rattling the cage: computational models of chaperonin-mediated protein folding

Chaperonins are known to maintain the stability of the proteome by facilitating the productive folding of numerous misfolded or aggregation-prone proteins and are thus essential for cell viability. Despite their established importance, the mechanism by which chaperonins facilitate protein folding remains unknown. Computer simulation techniques are now being employed to complement experimental ones in order to shed light on this mystery. Here we review previous computational models of chaperonin-mediated protein folding in the context of the two main hypotheses for chaperonin function: iterative annealing and landscape modulation. We then discuss new results pointing to the importance of solvent (a previously neglected factor) in chaperonin activity. We conclude with our views on the future role of simulation in studying chaperonin activity as well as protein folding in other biologically relevant confined contexts.

Do chaperonins boost protein yields by accelerating folding or preventing aggregation?

The GroEL chaperonin has the ability to behave as an unfoldase, repeatedly denaturing proteins upon binding, which in turn can free them from kinetic traps and increase their folding rates. The complex formed by GroEL+GroES+ATP can also act as an infinite dilution cage, enclosing proteins within a protective container where they can fold without danger of aggregation. Controversy remains over which of these two properties is more critical to the GroEL/ES chaperonin's function. We probe the importance of the unfoldase nature of GroEL under conditions where aggregation is the predominant protein degradation pathway. We consider the effect of a hypothetical mutation to GroEL which increases the cycle frequency of GroEL/ES by increasing the rate of hydrolysis of GroEL-bound ATP. Using a simple kinetic model, we show that this modified chaperonin would be self-defeating: any potential reduction in folding time would be negated by an increase in time spent in the bulk, causing an increase in aggregation and a net decrease in protein folding yields.

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.

Monomer adds to preformed structured oligomers of Abeta-peptides by a two-stage dock-lock mechanism

Nonfibrillar soluble oligomers, which are intermediates in the transition from monomers to amyloid fibrils, may be the toxic species in Alzheimer's disease. To monitor the early events that direct assembly of amyloidogenic peptides we probe the dynamics of formation of (Abeta(16-22))(n) by adding a monomer to a preformed (Abeta(16-22))(n-1) (n = 4-6) oligomer in which the peptides are arranged in an antiparallel beta-sheet conformation. All atom molecular dynamics simulations in water and multiple long trajectories, for a cumulative time of 6.9 mus, show that the oligomer grows by a two-stage dock-lock mechanism. The largest conformational change in the added disordered monomer occurs during the rapid ( approximately 50 ns) first dock stage in which the beta-strand content of the monomer increases substantially from a low initial value. In the second slow-lock phase, the monomer rearranges to form in register antiparallel structures. Surprisingly, the mobile structured oligomers undergo large conformational changes in order to accommodate the added monomer. The time needed to incorporate the monomer into the fluid-like oligomer grows even when n = 6, which suggests that the critical nucleus size must exceed six. Stable antiparallel structure formation exceeds hundreds of nanoseconds even though frequent interpeptide collisions occur at elevated monomer concentrations used in the simulations. The dock-lock mechanism should be a generic mechanism for growth of oligomers of amyloidogenic peptides.

Dynamics of allosteric transitions in GroEL

The chaperonin GroEL-GroES, a machine that helps proteins to fold, cycles through a number of allosteric states, the T state, with high affinity for substrate proteins, the ATP-bound R state, and the R" (GroEL-ADP-GroES) complex. Here, we use a self-organized polymer model for the GroEL allosteric states and a general structure-based technique to simulate the dynamics of allosteric transitions in two subunits of GroEL and the heptamer. The T --> R transition, in which the apical domains undergo counterclockwise motion, is mediated by a multiple salt-bridge switch mechanism, in which a series of salt-bridges break and form. The initial event in the R -->R" transition, during which GroEL rotates clockwise, involves a spectacular outside-in movement of helices K and L that results in K80-D359 salt-bridge formation. In both the transitions there is considerable heterogeneity in the transition pathways. The transition state ensembles (TSEs) connecting the T, R, and R" states are broad with the TSE for the T --> R transition being more plastic than the R --> R" TSE.

Folding on the chaperone: yield enhancement through loose binding

A variety of small cageless chaperones have been discovered that can assist protein folding without the consumption of ATP. These include mini-chaperones (catalytically active fragments of larger chaperones), as well as small proteins such as alpha-casein and detergents acting as "artificial chaperones." These chaperones all possess exposed hydrophobic patches on their surface that act as recognition sites for misfolded proteins. They lack the complexity of chaperonins (that encapsulate proteins in their inner rings) and their study can offer insight into the minimal requirements for chaperone function. We use molecular dynamics simulations to investigate how a cageless chaperone, modeled as a sphere of tunable hydrophobicity, can assist folding of a substrate protein. We find that under steady-state (non-stress) conditions, cageless chaperones that bind to a single substrate protein increase folding yields by reducing the time the substrate spends in an aggregation-prone state in a dual manner: (a) by competing for aggregation-prone hydrophobic sites on the surface of a protein, hence reducing the time the protein spends unprotected in the bulk and (b) by accelerating folding rates of the protein. In both cases, the chaperone must bind to and hold the protein loosely enough to allow the protein to change its conformation and fold while bound. Loose binding may enable small cageless chaperones to help proteins fold and avoid aggregation under steady-state conditions, even at low concentrations, without the consumption of ATP.

Residues in substrate proteins that interact with GroEL in the capture process are buried in the native state

We have used a bioinformatic approach to predict the natural substrate proteins for the Escherichia coli chaperonin GroEL based on two simple criteria. Natural substrate proteins should contain binding motifs similar in sequence to the mobile loop peptide of GroES that displaces the binding motif during the chaperonin cycle. Secondly, each substrate protein should contain multiple copies of the binding motif so that the chaperonin can perform "work" on the substrate protein. To validate these criteria, we have used a database of 252 proteins that have been experimentally shown to interact with the chaperonin machinery in vivo. More than 80% are identified by these criteria. The binding motifs of all 79 proteins in the database with a known three-dimensional structure are buried (<50% solvent-accessible surface area) in the native state. Our results show that the binding motifs are inaccessible in the native state but become solvent-exposed in unfolded state, thus enabling GroEL to distinguish between unfolded and native states. The structures of the binding motif in the native states of the substrate proteins include alpha-helices, beta-strands, and random coils. The diversity of secondary structures implies that there are large and varied conformational transitions in the recognition motifs after their displacement by the mobile loops of GroES.

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.