Amnon Horovitz

Strength and co-operativity of contributions of surface salt bridges to protein stability.

Many of the interactions that stabilize proteins are co-operative and cannot be reduced to a sum of pairwise interactions. Such interactions may be analysed by protein engineering methods using multiple thermodynamic cycles comprising wild-type protein and all combinations of mutants in the interacting residues. There is a triad of charged residues on the surface of barnase, comprising residues Asp8, Asp12 and Arg110, that interact by forming two exposed salt bridges. The three residues have been mutated to alanine to give all the single, double and triple mutants. The free energies of unfolding of wild-type and the seven mutant proteins have been determined and the results analysed to give the contributions of the residues in the two salt bridges to protein stability. It is possible to isolate the energies of forming the salt bridges relative to the solvation of the separated ions by water. In the intact triad, the apparent contribution to the stabilization energy of the protein of the salt bridge between Asp12 and Arg110 is -1.25 kcal mol-1, whereas that of the salt bridge between Asp8 with Arg110 is -0.98 kcal mol-1. The strengths of the two salt bridges are coupled: the energy of each is reduced by 0.77 kcal mol-1 when the other is absent. The salt-linked triad, relative to alanine residues at the same positions, does not contribute to the stability of the protein since the favourable interactions of the salt bridges are more than offset by other electrostatic and non-electrostatic energy terms. Salt-linked triads occur in other proteins, for example, haemoglobin, where the energy of only the salt-bridge term is important and so the coupling of salt bridges could be of general importance to the stability and function of proteins.

Strategy for analysing the co-operativity of intramolecular interactions in peptides and proteins

Double mutant cycles enable the measurement of pairwise interactions in proteins. This method is extended for mutations at any number of positions in the protein. This provides a way for determining the context dependence of pairwise interactions on other neighbouring residues.

Double-mutant cycles: a powerful tool for analyzing protein structure and function

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.

α-Helix stability in proteins: II. Factors that influence stability at an internal position

The solvent-exposed residue Ala32 in the second alpha-helix of barnase was replaced by all other naturally occurring amino acids and the concomitant effects on the protein stability were determined. The results are assumed to reflect both the distinct conformational preferences of the different amino acids and also possible intrahelical interactions. The conformational preferences may be fully rationalized by invoking only a few physical principles. The results agree well with recently experimentally determined rank-order of helix-forming tendencies determined on a model peptide. There is very weak correlation between the results and the experimental host-guest values. There is a weak correlation between our results and the statistical helix propensities and a slightly better correlation with the positional-dependent statistical parameters of J. S. Richardson, and D. C. Richardson.

Mapping pathways of allosteric communication in GroEL by analysis of correlated mutations.

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. Proteins 2002;48:611–617.

Co-operative interactions during protein folding.

The theory for measuring co-operativity between interactions in proteins by protein engineering experiments is developed by introducing a procedure for analysing increasing orders of synergy in a protein with increasing numbers of residues. The (pairwise) interaction energy (delta 2Gint) between two side-chains may be measured experimentally by a double-mutant cycle consisting of the wild-type protein, the two single mutants and the double mutant. This procedure may be extended to three residues to give a value for delta 3Gint for a triple-mutant cube, and to higher orders using multi-dimensional mutant space. We now show that delta 3Gint is the excess energy of adding all three chains compared with the sum of all the pairwise values of delta 2Gint for each of the constituent double-mutant cycles and the sum of all the single addition energies. This physical interpretation extends to higher orders of mutation. delta nGint (i.e. the interaction energy for n residues), thus, reveals the layers of synergy in interactions as a protein is built up. This procedure is applied to measuring changes in synergy during the refolding of barnase for the triad of salt-linked residues Asp8, Asp12 and Arg110, which are mutated to alanine residues. The value of delta 3Gint in the folded structure is 0.77(+/- 0.06) kcal mol-1 (i.e. the triad is 0.77 kcal mol-1 more stable than expected from the sum of the individual pairwise interactions and single contributions). The value of delta 3Gint is still significant in the transition state for unfolding (0.60(+/- 0.07) kcal mol-1) and in the folding intermediate (0.60(+/- 0.13 kcal mol-1)). These results show that synergistic interactions exist in barnase, in its transition state for unfolding and in a refolding intermediate. A direct measurement of the change of co-operativity between the folded state and the transition state for unfolding shows a decrease of 0.17(+/- 0.04) kcal mol-1, suggesting that the initial stages of protein unfolding may be accompanied by some loosening of structure in parts that still interact. The similar extent of co-operativity in the transition state for unfolding and the intermediate in refolding suggests that the intermediate is homogeneous, at least in the region of the salt-linked triad, as heterogeneity would lower the co-operativity.

Sequential ATP-induced allosteric transitions of the cytoplasmic chaperonin containing TCP-1 revealed by EM analysis.

The eukaryotic cytoplasmic chaperonin containing TCP-1 (CCT) is a hetero-oligomeric complex that assists the folding of actins, tubulins and other proteins in an ATP-dependent manner. To understand the allosteric transitions that occur during the functional cycle of CCT, we imaged the chaperonin complex in the presence of different ATP concentrations. Labeling by monoclonal antibodies that bind specifically to the CCTα and CCTδ subunits enabled alignment of all the CCT subunits of a given type in different particles. The analysis shows that the apo state of CCT has considerable apparent conformational heterogeneity that decreases with increasing ATP concentration. In contrast with the concerted allosteric switch of GroEL, ATP-induced conformational changes in CCT are found to spread around the ring in a sequential fashion that may facilitate domain-by-domain substrate folding. The approach described here can be used to unravel the allosteric mechanisms of other ring-shaped molecular machines.

Allosteric mechanisms can be distinguished using structural mass spectrometry

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 MS—which 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.

Non-additivity in protein-protein interactions

The energy of binding between proteins may be seen as the sum of the contributions of the individual amino acid residues. These contributions are additive when the binding energy, due to different amino acid residues, is independent of the interactions between amino acids in the same polypeptide chain. A measure of non-additivity is the coupling free energy. In this communication it is shown that: (1) the coupling free energy is the sum of intramolecular and intermolecular contributions; and (2), when additivity exists, experimentally determined values for the free energy of transfer of amino acids from water to the hydrophobic protein-protein interface are a very good approximation of their contribution to the energy of binding. Additivity cycles can be useful in determining the precise conditions where this approximation holds.

Nested allosteric interactions in the cytoplasmic chaperonin containing TCP-1

Initial rates of ATP hydrolysis by the chaperonin containing TCP‐1 (CCT) from bovine testis were measured as a function of ATP concentration. Two allosteric transitions are observed: one at relatively low concentrations of ATP (<100 μM) and the second at higher concentrations of ATP. The data suggest that CCT has positive intra‐ring cooperativity and negative inter‐ring cooperativity in ATP hydrolysis, with respect to ATP, as previously observed in the case of GroEL. It is shown that the relatively weak positive intra‐ring cooperativity found in the case of CCT may be due to heterogeneity in its subunit composition. Our results suggest that nested allosteric behavior may be common to chaperone double‐ring systems.