Search

Determining the network of residues that transmit allosteric signals is crucial to understanding the function of biological nanomachines. During the course of a reaction cycle, biological machines in general, and Escherichia coli chaperonin GroEL in particular, undergo large-scale conformational changes in response to ligand binding. Normal mode analyses, based on structure-based coarse-grained models where each residue is represented by an alpha carbon atom, have been widely used to describe the motions encoded in the structures of proteins. Here, we propose a new Calpha-side chain elastic network model of proteins that includes information about the physical identity of each residue and accurately accounts for the side-chain topology and packing within the structure. Using the Calpha-side chain elastic network model and the structural perturbation method, which probes the response of a local perturbation at a given site at all other sites in the structure, we determine the network of key residues (allostery wiring diagram) responsible for the T-->R and R''-->T transitions in GroEL. A number of residues, both within a subunit and at the interface of two adjacent subunits, are found to be at the origin of the positive cooperativity in the ATP-driven T-->R transition. Of particular note are residues G244, R58, D83, E209, and K327. Of these, R38, D83, and K327 are highly conserved. G244 is located in the apical domain at the interface between two subunits; E209 and K327 are located in the apical domain, toward the center of a subunit; R58 and D83 are equatorial domain residues. The allostery wiring diagram shows that the network of residues are interspersed throughout the structure. Residues D83, V174, E191, and D359 play a critical role in the R''-->T transition, which implies that mutations of these residues would compromise the ATPase activity. D83 and E191 are also highly conserved; D359 is moderately conserved. The negative cooperativity between the rings in the R''-->T transition is orchestrated through several interface residues within a single ring, including N10, E434, D435, and E451. Signal from the trans ring that is transmitted across the interface between the equatorial domains is responsible for the R''-->T transition. The cochaperonin GroES plays a passive role in the R''-->T transition. Remarkably, the binding affinity of GroES for GroEL is allosterically linked to GroEL residues 350-365 that span helices K and L. The movements of helices K and L alter the polarity of the cavity throughout the GroEL functional cycle and undergo large-scale motions that are anticorrelated with the other apical domain residues. The allostery wiring diagrams for the T-->R and R''-->T transitions of GroEL provide a microscopic foundation for the cooperativity (anticooperativity) within (between) the ring (rings). Using statistical coupling analysis, we extract evolutionarily linked clusters of residues in GroEL and GroES. We find that several substrate protein binding residues as well as sites related to ATPase activity belong to a single functional network in GroEL. For GroES, the mobile loop residues and GroES/GroES interface residues are linked.

The bacterial chaperonin GroEL and the co-chaperonin GroES assist in the folding of a number of structurally unrelated substrate proteins (SPs). In the absence of chaperonins, SP folds by the kinetic partitioning mechanism (KPM), according to which a fraction of unfolded molecules reaches the native state directly, while the remaining fraction gets trapped in a potentially aggregation-prone misfolded state. During the catalytic reaction cycle, GroEL undergoes a series of allosteric transitions (T<-->R-->R"-->T) triggered by SP capture, ATP binding and hydrolysis, and GroES binding. We developed a general kinetic model that takes into account the coupling between the rates of the allosteric transitions and the folding and aggregation of the SP. Our model, in which the GroEL allosteric rates and SP-dependent folding and aggregation rates are independently varied without prior assumption, quantitatively fits the GroEL concentration-dependent data on the yield of native ribulose bisphosphate carboxylase/oxygenase (Rubisco) as a function of time. The extracted kinetic parameters for the GroEL reaction cycle are consistent with the available values from independent experiments. In addition, we also obtained physically reasonable parameters for the kinetic steps in the reaction cycle that are difficult to measure. If experimental values for GroEL allosteric rates are used, the time-dependent changes in native-state yield at eight GroEL concentrations can be quantitatively fit using only three SP-dependent parameters. The model predicts that the differences in the efficiencies (as measured by yields of the native state) of GroEL, single-ring mutant (SR1), and variants of SR1, in the rescue of mitochondrial malate dehydrogenase, citrate synthase, and Rubisco, are related to the large variations in the allosteric transition rates. We also show that GroEL/S mutants that efficiently fold one SP at the expense of all others are due to a decrease in the rate of a key step in the reaction cycle, which implies that wild-type GroEL has evolved as a compromise between generality and specificity. We predict that, under maximum loading conditions and saturating ATP concentration, the efficiency of GroEL (using parameters for Rubisco) depends predominantly on the rate of R-->R" transition, while the equilibrium constant of the T<-->R has a small effect only. Both under sub- and superstoichiometric GroEL concentrations, enhanced efficiency is achieved by rapid turnover of the reaction cycle, which is in accord with the predictions of the iterative annealing mechanism. The effects are most dramatic at substoichiometric conditions (most relevant for in vivo situations) when SP aggregation can outcompete capture of SP by chaperonins.