Out-of-equilibrium conformational cycling of GroEL under saturating ATP concentrations

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 fast are the motions of tertiary-structure elements in proteins?

Protein motions occur on multiple time and distance scales. Large-scale motions of protein tertiary-structure elements, i.e., domains, are particularly intriguing as they are essential for the catalytic activity of many enzymes and for the functional cycles of protein machines and motors. Theoretical estimates suggest that domain motions should be very fast, occurring on the nanosecond or microsecond time scales. Indeed, free-energy barriers for domain motions are likely to involve salt bridges, which can break in microseconds. Experimental methods that can directly probe domain motions on fast time scales have appeared only in recent years. This Perspective discusses briefly some of these techniques, including nuclear magnetic resonance and single-molecule fluorescence spectroscopies. We introduce a few recent studies that demonstrate ultrafast domain motions and discuss their potential roles. Particularly surprising is the observation of tertiary-structure element dynamics that are much faster than the functional cycles in some protein machines. These swift motions can be rationalized on a case-by-case basis. For example, fast domain closure in multi-substrate enzymes may be utilized to optimize relative substrate orientation. Whether a large mismatch in time scales of conformational dynamics vs functional cycles is a general design principle in proteins remains to be determined.

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.

Molecular chaperones maximize the native state yield on biological times by driving substrates out of equilibrium

Molecular chaperones facilitate the folding of proteins and RNA in vivo. Under physiological conditions, the in vitro folding of Tetrahymena ribozyme by the RNA chaperone CYT-19 behaves paradoxically; increasing the chaperone concentration reduces the yield of native ribozymes. In contrast, the protein chaperone GroEL works as expected; the yield of the native substrate increases with chaperone concentration. The discrepant chaperone-assisted ribozyme folding thus contradicts the expectation that it operates as an efficient annealing machine. To resolve this paradox, we propose a minimal stochastic model based on the Iterative Annealing Mechanism (IAM) that offers a unified description of chaperone-mediated folding of both proteins and RNA. Our theory provides a general relation that quantitatively predicts how the yield of native states depends on chaperone concentration. Although the absolute yield of native states decreases in the Tetrahymena ribozyme, the product of the folding rate and the steady-state native yield increases in both cases. By using energy from ATP hydrolysis, both CYT-19 and GroEL drive their substrate concentrations far out of equilibrium, thus maximizing the native yield in a short time. This also holds when the substrate concentration exceeds that of GroEL. Our findings satisfy the expectation that proteins and RNA be folded by chaperones on biologically relevant time scales, even if the final yield is lower than what equilibrium thermodynamics would dictate. The theory predicts that the quantity of chaperones in vivo has evolved to optimize native state production of the folded states of RNA and proteins in a given time.

Tracking Equilibrium and Nonequilibrium Shifts in Data with TREND

Principal component analysis (PCA) discovers patterns in multivariate data that include spectra, microscopy, and other biophysical measurements. Direct application of PCA to crowded spectra, images, and movies (without selecting peaks or features) was shown recently to identify their equilibrium or temporal changes. To enable the community to utilize these capabilities with a wide range of measurements, we have developed multiplatform software named TREND to Track Equilibrium and Nonequilibrium population shifts among two-dimensional Data frames. TREND can also carry this out by independent component analysis. We highlight a few examples of finding concurrent processes. TREND extracts dual phases of binding to two sites directly from the NMR spectra of the titrations. In a cardiac movie from magnetic resonance imaging, TREND resolves principal components (PCs) representing breathing and the cardiac cycle. TREND can also reconstruct the series of measurements from selected PCs, as illustrated for a biphasic, NMR-detected titration and the cardiac MRI movie. Fidelity of reconstruction of series of NMR spectra or images requires more PCs than needed to plot the largest population shifts. TREND reads spectra from many spectroscopies in the most common formats (JCAMP-DX and NMR) and multiple movie formats. The TREND package thus provides convenient tools to resolve the processes recorded by diverse biophysical methods.

Quantifying kinetics from time series of single-molecule Förster resonance energy transfer efficiency histograms

Single-molecule fluorescence spectroscopy is a powerful approach for probing biomolecular structure and dynamics, including protein folding. For the investigation of nonequilibrium kinetics, Förster resonance energy transfer combined with confocal multiparameter detection has proven particularly versatile, owing to the large number of observables and the broad range of accessible timescales, especially in combination with rapid microfluidic mixing. However, a comprehensive kinetic analysis of the resulting time series of transfer efficiency histograms and complementary observables can be challenging owing to the complexity of the data. Here we present and compare three different methods for the analysis of such kinetic data: singular value decomposition, multivariate curve resolution with alternating least square fitting, and model-based peak fitting, where an explicit model of both the transfer efficiency histogram of each species and the kinetic mechanism of the process is employed. While each of these methods has its merits for specific applications, we conclude that model-based peak fitting is most suitable for a quantitative analysis and comparison of kinetic mechanisms.

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.

Non-equilibrium conformational dynamics in the function of molecular chaperones

Why do chaperones need ATP hydrolysis to help proteins reach their native, functional states? In this review, we highlight the most recent experimental and theoretical evidences suggesting that ATP hydrolysis allows molecular chaperones to escape the bounds imposed by equilibrium thermodynamics. We argue here that energy consumption must be fully taken into account to understand the mechanism of these intrinsically non-equilibrium machines and we propose a novel perspective in the way the relation between function and ATP hydrolysis is viewed.

Substrate protein switches GroE chaperonins from asymmetric to symmetric cycling by catalyzing nucleotide exchange

The complex kinetics of Pi and ADP release by the chaperonin GroEL/GroES is influenced by the presence of unfolded substrate protein (SP). Without SP, the kinetics of Pi release are described by four phases: a "lag," a "burst" of ATP hydrolysis by the nascent cis ring, a "delay" caused by ADP release from the nascent trans ring, and steady-state ATP hydrolysis. The release of Pi precedes the release of ADP. The rate-determining step of the asymmetric cycle is the release of ADP from the trans ring of the GroEL-GroES1 "bullet" complex that is, consequently, the predominant species. In the asymmetric cycle, the two rings of GroEL function alternately, 180° out of phase. In the presence of SP, a change in the kinetic mechanism occurs. With SP present, the kinetics of ADP release are also described by four phases: a lag, a "surge" of ADP release attributable to SP-induced ADP/ATP exchange, and a "pause" during which symmetrical "football" particles are formed, followed by steady-state ATP hydrolysis. SP catalyzes ADP/ATP exchange on the trans ring. Now ADP release precedes the release of Pi, and the rate-determining step of the symmetric cycle becomes the hydrolysis of ATP by the symmetric GroEL-GroES2 football complex that is, consequently, the predominant species. A FRET-based analysis confirms that asymmetric GroEL-GroES1 bullets predominate in the absence of SP, whereas symmetric GroEL-GroES2 footballs predominate in the presence of SP. This evidence suggests that symmetrical football particles are the folding functional form of the chaperonin machine in vivo.

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.

Crystal structure of a GroEL-ADP complex in the relaxed allosteric state at 2.7 Å resolution

The chaperonin proteins GroEL and GroES are cellular nanomachines driven by the hydrolysis of ATP that facilitate the folding of structurally diverse substrate proteins. In response to ligand binding, the subunits of a ring cycle in a concerted manner through a series of allosteric states (T, R, and R″), enabling work to be performed on the substrate protein. Removing two salt bridges that ordinarily break during the allosteric transitions of the WT permitted the structure of GroEL-ADP in the R state to be solved to 2.7 Å resolution. Whereas the equatorial domain displays almost perfect sevenfold symmetry, the apical domains, to which substrate proteins bind, and to a lesser extent, the intermediate domains display a remarkable asymmetry. Freed of intersubunit contacts, the apical domain of each subunit adopts a different conformation, suggesting a flexibility that permits interaction with diverse substrate proteins. This result contrasts with a previous cryo-EM study of a related allosteric ATP-bound state at lower resolution. After artificially imposing sevenfold symmetry it was concluded that a GroEL ring in the R-ATP state existed in six homogeneous but slightly different states. By imposing sevenfold symmetry on each of the subunits of the crystal structure of GroEL-ADP, we showed that the synthetic rings of (X-ray) GroEL-ADP and (cryo-EM) GroEL-ATP are structurally closely related. A deterministic model, the click stop mechanism, that implied temporal transitions between these states was proposed. Here, however, these conformational states are shown to exist as a structurally heterogeneous ensemble within a single ring.

Biomolecular robotics for chemomechanically driven guest delivery fuelled by intracellular ATP

The development of nanocarriers that selectively release guest molecules on sensing a particular biological signal is being actively pursued in nanomedicine for diagnostic and therapeutic purposes. Here we report a protein-based nanocarrier that opens in the presence of intracellular adenosine-5'-triphosphate (ATP). The nanocarrier consists of multiple barrel-shaped chaperonin units assembled through coordination with Mg(2+) into a tubular structure that protects guest molecules against biological degradation. When its surface is functionalized with a boronic acid derivative, the nanocarrier is able to enter cells. The hydrolysis of intracellular ATP into adenosine-5'-diphosphate (ADP) induces conformational changes of the chaperonin units, which in turns generate a mechanical force that leads to the disassembly of the tube and release of the guests. This scission occurs with a sigmoidal dependence on ATP concentration, which means that the nanocarrier can differentiate biological environments in terms of the concentration of ATP for selective guest release. Furthermore, biodistribution tests reveal preferential accumulation of the nanocarriers in a tumour tissue.

The multi-targeted effects of Chrysanthemum herb extract against Escherichia coli O157:H7

The Chrysanthemum lavandulifolium extract, which includes chrysoeriol, sudachitin, and acacetin, has excellent antibiotic effects on Escherichia coli O157:H7 (E. coli O157). A notable point is that the antibiotic targets of the herb extract are similar to the targets of commonly used antibiotic drugs, including bacterial cell wall biosynthesis, bacterial protein synthesis, and bacterial DNA replication and repair. In addition, the herbal antibiotic inhibits the etiological factors that contribute to the pathogenic property. The herbal sample was extracted and fractionated and then inoculated through a disk diffusion method to confirm its antibiotic effect against E. coli O157. Total RNA was isolated from the affected bacterial cells, and its expression level was analyzed through a microarray analysis. To confirm the accuracy of the microarray data, a real-time PCR was performed. Three active compounds, chrysoeriol, sudachitin, and acacetin, were identified with a high-performance liquid chromatography-electrospray ionization/mass spectrometry chromatogram, and the disk diffusion study confirmed that chrysoeriol and sudachitin contribute to the antibiotic properties of the herb extract. The results demonstrate that the multi-target efficacy of the herbal sample may indicate the potential for the development of more effective and safer drugs that will act as substitutes for existing antibiotics.

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.

Dynamics, flexibility and ligand-induced conformational changes in biological macromolecules: a computational approach

Biomolecules possess important dynamical properties that enable them to adapt and alternate their conformation as a response to environmental stimuli. Recent advancements in computational resources and methodology allow a higher capability to mimic in vitro conditions and open up the possibility of studying large systems over longer timescales. Here, we describe commonly used computational approaches for studying the dynamic properties of proteins. We review a selected set of simulation studies on ligand-induced changes in the chaperonin GroEL-GroES, a molecular folding machine, maltose-binding protein, a prototypical member of the periplasmic binding proteins, and the bacterial ribosomal A-site, focusing on aminoglycoside antibiotic recognition. We also discuss a recent quantitative reconstruction of the binding process of benzamidine and trypsin. These studies contribute to the understanding and further development of the medicinal regulation of large biomolecular systems.

Fluorescence correlation spectroscopy and allostery: the case of GroEL

Fluorescence correlation spectroscopy (FCS) is an experimental technique in which the equilibrium fluctuations of the fluorescent signal of molecules diffusing through a focused laser beam are measured. An autocorrelation analysis of these fluctuations provides information on dynamic processes, such as allosteric transitions, that the molecules undergo provided that they are fast relative to the diffusion time through the beam. In cases when the dynamics are slow relative to the diffusion time through the beam, FCS curves can yield information about the number of conformational states and their relative populations. Hence, FCS can be used to investigate allosteric systems with either slow or fast dynamics but the type of information gained in these two situations is different.Here, the utility of the FCS technique is exemplified in the case of the single-ring version of the Escherichia coli molecular chaperone GroEL that interconverts with relatively slow dynamics between two allosteric states: a T state with low affinity for ATP and an R state with high affinity for ATP. Thermodynamic analysis suggests that the T-state population should become negligible with increasing ATP concentrations, in conflict with the requirement for conformation cycling, which is essential for the operation of molecular machines. Surprisingly, FCS experiments showed that, even at ATP saturation, ~50% of the molecules still populate the T state at any instance of time, indicating constant out-of-equilibrium cycling between T and R.

Conformational sampling and nucleotide-dependent transitions of the GroEL subunit probed by unbiased molecular dynamics simulations

GroEL is an ATP dependent molecular chaperone that promotes the folding of a large number of substrate proteins in E. coli. Large-scale conformational transitions occurring during the reaction cycle have been characterized from extensive crystallographic studies. However, the link between the observed conformations and the mechanisms involved in the allosteric response to ATP and the nucleotide-driven reaction cycle are not completely established. Here we describe extensive (in total long) unbiased molecular dynamics (MD) simulations that probe the response of GroEL subunits to ATP binding. We observe nucleotide dependent conformational transitions, and show with multiple 100 ns long simulations that the ligand-induced shift in the conformational populations are intrinsically coded in the structure-dynamics relationship of the protein subunit. Thus, these simulations reveal a stabilization of the equatorial domain upon nucleotide binding and a concomitant “opening” of the subunit, which reaches a conformation close to that observed in the crystal structure of the subunits within the ADP-bound oligomer. Moreover, we identify changes in a set of unique intrasubunit interactions potentially important for the conformational transition.

Molecular machines convert chemical energy to mechanical work in the process of carrying out their specific tasks. Often these proteins are fueled by ATP binding and hydrolysis, enabling switching between different conformations. The ATP-dependent chaperone GroEL is a molecular machine that opens and closes its barrel-like structure in order to provide a folding cage for unfolded proteins. The quest to fully understand and control GroEL and other molecular machines is enhanced by complementing experimental work with computational approaches. Here, we provide a description of the molecular basis for the conformational changes in the GroEL subunit by performing extensive molecular dynamics simulations. The simulations sample the conformational population for the different nucleotide-free and bound states in the isolated subunit. The results reveal that the conformations of the subunit when isolated resemble those of the subunit integrated in the GroEL complex. Moreover, the molecular dynamics simulations allow following detailed changes in individual interatomic interactions brought about by ATP-binding.

Single-molecule spectroscopy of protein folding in a chaperonin cage

Molecular chaperones are known to be essential for avoiding protein aggregation in vivo, but it is still unclear how they affect protein folding mechanisms. We use single-molecule Förster resonance energy transfer to follow the folding of a protein inside the GroEL/GroES chaperonin cavity over a time range from milliseconds to hours. Our results show that confinement in the chaperonin decelerates the folding of the C-terminal domain in the substrate protein rhodanese, but leaves the folding rate of the N-terminal domain unaffected. Microfluidic mixing experiments indicate that strong interactions of the substrate with the cavity walls impede the folding process, but the folding hierarchy is preserved. Our results imply that no universal chaperonin mechanism exists. Rather, a competition between intra- and intermolecular interactions determines the folding rates and mechanisms of a substrate inside the GroEL/GroES cage.