Coulombic interactions between charges on the surface of proteins contribute to stability. It is difficult, however, to estimate their importance by protein engineering methods because mutation of one residue in an ion pair alters the energetics of many interactions in addition to the coulombic energy between the two components. We have estimated the interaction energy between two charged residues, Asp-12 and Arg-16, in an alpha-helix on the surface of a barnase mutant by invoking a double-mutant cycle involving wild-type enzyme (Asp-12, Thr-16), the single mutants Thr----Arg-16 and Asp----Ala-12, and the double mutant Asp----Ala-12, Thr----Arg-16. The changes in free energy of unfolding of the single mutants are not additive because of the coulombic interaction energy. Additivity is restored at high concentrations of salt that shield electrostatic interactions. The geometry of the ion pair in the mutant was assumed to be the same as that in the highly homologous ribonuclease from Bacillus intermedius, binase, which has Asp-12 and Arg-16 in the native enzyme. The ion pair does not form a hydrogen-bonded salt bridge, but the charges are separated by 5-6 A. The mutant barnase containing the ion pair Asp-12/Arg-16 is more stable than wild type by 0.5 kcal/mol, but only a part of the increased stability is attributable to the electrostatic interaction. We present a formal analysis of how double-mutant cycles can be used to measure the energetics of pairwise interactions.

No abstract

No abstract

A double-mutant cycle involves wild-type protein, two single mutants and the corresponding double mutant protein. If the change in free energy associated with a structural or functional property of the protein upon a double mutation differs from the sum of changes in free energy due to the single mutations, then the residues at the two positions are coupled. Such coupling reflects either direct or indirect interactions between these residues. Double-mutant cycle analysis can be used to measure the strength of intramolecular and intermolecular pairwise interactions in proteins or protein-ligand complexes with known structure. Double-mutant cycles can also be employed to characterize structures that are inaccessible to NMR and X-ray crystallography, such as those of transition states for protein folding, ligand binding and enzyme catalysis, or of membrane proteins. Multidimensional mutant cycle analysis can be used to measure higher-order cooperativity between intramolecular or intermolecular interactions. In the absence of coupling between residues, prediction of mutational effects is possible by assuming their additivity.

Initial rates of ATP hydrolysis by wild-type GroEL were measured as a function of ATP concentration from 0 to 0.8 mM. Two allosteric transitions are observed: one at relatively low ATP concentrations (< or = 100 microM) and the second at higher concentrations of ATP with respective midpoints of about 16 and 160 microM. Two allosteric transitions were previously observed also in the case of the Arg-196-->Ala GroEL mutant [Yifrach, O., & Horovitz, A. (1994) J. Mol. Biol. 243, 397-401]. On the basis of these observations a mathematical model for nested cooperativity in ATP hydrolysis by GroEL is developed in which there are two levels of allostery: one within each ring and the second between rings. In the first level, each hepatameric ring is in equilibrium between the T and R states, in accordance with the Monod-Wyman-Changeux (MWC) model of cooperativity [Monod et al. (1965) J. Mol. Biol. 12, 88-118]. A second level of allostery is between the rings of the GroEL particle which undergoes sequential Koshland-Némethy-Filmer (KNF)-type transitions from the TT state via the TR state to the RR state [Koshland et al. (1966) Biochemistry 5, 365-385]. Using our model, we estimate the values of the Hill coefficient for the negative cooperativity between rings in wild-type GroEL and the Arg-196-->Ala mutant to be 0.003 (+/- 0.001) and 0.07 (+/- 0.02), respectively. The inter-ring coupling free energies in wild-type GroEL and the Arg-196-->Ala mutant are -7.5 (+/- 0.4) and -3.9 (+/- 0.3) kcal mol-1, respectively.

Allosteric regulation plays an important role in many biological processes, such as signal transduction, transcriptional regulation, and metabolism. Allostery is rooted in the fundamental physical properties of macromolecular systems, but its underlying mechanisms are still poorly understood. A collection of contributions to a recent interdisciplinary CECAM (Center Européen de Calcul Atomique et Moléculaire) workshop is used hereto provide an overview of the progress and remaining limitations in the understanding of the mechanistic foundations of allostery gained from computational and experimental analyses of real protein systems and model systems. The main conceptual frameworks instrumental in driving the field are discussed. We illustrate the role of these frameworks in illuminating molecular mechanisms and explaining cellular processes, and describe some of their promising practical applications in engineering molecular sensors and informing drug design efforts.

An interesting example of an allosteric protein is the chaperonin GroEL. It undergoes adenosine 5'-triphosphate-induced conformational changes that are reflected in binding of adenosine 5'-triphosphate with positive cooperativity within rings and negative cooperativity between rings. Herein, correlated mutations in chaperonins are analyzed to unravel routes of allosteric communication in GroEL and in its complex with its co-chaperonin GroES. It is shown that analysis of correlated mutations in the chaperonin family can provide information about pathways of allosteric communication within GroEL and between GroEL and GroES. The results are discussed in the context of available structural, genetic, and biochemical data concerning short- and long-range interactions in the GroE system.

The activity of many proteins, including metabolic enzymes, molecular machines, and ion channels, is often regulated by conformational changes that are induced or stabilized by ligand binding. In cases of multimeric proteins, such allosteric regulation has often been described by the concerted Monod-Wyman-Changeux and sequential Koshland-Némethy-Filmer classic models of cooperativity. Despite the important functional implications of the mechanism of cooperativity, it has been impossible in many cases to distinguish between these various allosteric models using ensemble measurements of ligand binding in bulk protein solutions. Here, we demonstrate that structural MS offers a way to break this impasse by providing the full distribution of ligand-bound states of a protein complex. Given this distribution, it is possible to determine all the binding constants of a ligand to a highly multimeric cooperative system, and thereby infer its allosteric mechanism. Our approach to the dissection of allosteric mechanisms relies on advances in MSwhich provide the required resolution of ligand-bound states-and in data analysis. We validated our approach using the well-characterized Escherichia coli chaperone GroEL, a double-heptameric ring containing 14 ATP binding sites, which has become a paradigm for molecular machines. The values of the 14 binding constants of ATP to GroEL were determined, and the ATP-loading pathway of the chaperone was characterized. The methodology and analyses presented here are directly applicable to numerous other cooperative systems and are therefore expected to promote further research on allosteric systems.chaperonins | Hill coefficient M ultimeric proteins are often subject to allosteric regulation that is achieved by conformational changes induced or stabilized by ligand binding (1). Such allosteric regulation has been described by two classic models: (i) the Monod-WymanChangeux (MWC) model (2), in which conformational changes occur in a concerted manner and symmetry is conserved, and (ii) the Koshland-Némethy-Filmer (KNF) model (3), in which conformational changes take place in a sequential manner and symmetry is broken. In addition, it has been proposed more recently that conformational changes can take place in a probabilistic manner (4). The allosteric control of protein activity is frequently manifested in sigmoidal plots of initial reaction velocity or fractional saturation as a function of the ligand (substrate) concentration that indicates positive cooperativity in ligand binding. It has been impossible, however, to extract any mechanistic insights from these plots (5) because they only show how an average property of the ensemble (e.g., fractional saturation) changes with ligand concentration and do not reveal how the distribution of ligand-bound states changes with ligand concentration. Thus, for example, it is not possible to determine from such sigmoidal plots whether an allosteric transition takes place in a concerted MWC-like fashion (2) or via a sequential KNF-like mechanism (3). T...