Packing is a key selection factor in the evolution of protein hydrophobic cores.

The energy derived from optimized van der Waals interactions in closely packed, folded proteins has been proposed to be of similar energetic magnitude to hydrophobicity in stabilizing the native state. If packing is this energetically important, it should influence the evolution of protein core sequences. To test this hypothesis, the occurrence of various amino acid side chains in the major hydrophobic core of staphylococcal nuclease and 42 homologous proteins was determined. Most such positions in this protein family are usually isoleucine, leucine, or valine. Previously we have constructed and measured the stabilities of 12 single, 44 double, 64 triple, and 32 quadruple mutants, representing all possible permutations of these three side chains at two overlapping sets of four positions in the core of staphylococcal nuclease. The stabilities and interaction energies of those mutants with various combinations of the most common, or consensus, sequence were compared to the stabilities of all other mutants. Mutants which had the consensus side-chain combinations were not necessarily the most stable, but usually were found to have the best interaction energies. In other words, these proteins were far more stable than would be predicted from simply summing the observed energetic effects of the component single mutations, apparently reflecting particularly favorable packing interactions that are possible for the most common side chains. An additional 12 mutants which tested possible alternate explanations of the results were constructed. The stabilities and interaction energies of these mutants also support the conclusion that packing is a crucial determinant guiding the sequence evolution of protein cores.     

Increasing the thermostability of staphylococcal nuclease: implications for the origin of protein thermostability

Seven hyper-stable multiple mutants have been constructed in staphylococcal nuclease by various combinations of eight different stabilizing single mutants. The stabilities of these multiple mutants determined by guanidine hydrochloride denaturation were 3.4 to 5.6 kcal/mol higher than that of the wild-type. Their thermal denaturation midpoint temperatures were 12.6 to 22.9 deg. C higher than that of the wild-type. These are among the greatest increases in protein stability and thermal denaturation midpoint temperature relative to the wild-type yet attained. There has been great interest in understanding how proteins found in thermophilic organisms are stabilized. One frequently cited theory is that the packing of hydrophobic side-chains is improved in the cores of proteins isolated from thermophiles when compared to proteins from mesophiles. The crystal structures of four single and five multiple stabilizing mutants of staphylococcal nuclease were solved to high resolution. No large overall structural change was found, with most changes localized around the sites of mutation. Rearrangements were observed in the packing of side-chains in the major hydrophobic core, although none of the mutations was in the core. It is surprising that detailed structural analysis showed that packing had improved, with the volume of the mutant protein's hydrophobic cores decreasing as protein stability increased. Further, the number of van der Waals interactions in the entire protein showed an experimentally significant increase correlated with increasing stability. These results indicate that optimization of packing follows as a natural consequence of increased protein thermostability and that good packing is not necessarily the proximate cause of high stability. Another popular theory is that thermostable proteins have more electrostatic and hydrogen bonding interactions and these are responsible for the high stabilities. The mutants here show that increased numbers of electrostatic and hydrogen bonding interactions are not obligatory for large increases in protein stability.     

Energetics of side chain packing in staphylococcal nuclease assessed by exchange of valines, isoleucines, and leucines

To examine the importance of side chain packing to protein stability, each of the 11 leucines in staphylococcal nuclease was substituted with isoleucine and valine. The nine valines were substituted with leucine and isoleucine, while the five isoleucines, previously substituted with valine, were substituted with leucine and methionine. These substitutions conserve the hydrophobic character of these side chains but alter side chain geometry and, in some cases, size. In addition, eight threonine residues, previously substituted with valine, were substituted with isoleucine to test the importance of packing at sites normally not occupied by a hydrophobic residue. The stabilities of these 58 mutant proteins were measured by guanidine hydrochloride denaturation. To the best of our knowledge, this is the largest library of single packing mutants yet characterized. As expected, repacking stability effects are tied to the degree of side chain burial. The average energetic cost of moving a single buried methyl group was 0.9 kcal/mol, albeit with a standard deviation of 0.8 kcal/mol. This average is actually slightly greater than the value of 0.7-0.8 kcal/mol estimated for the hydrophobic transfer energy of a methylene from octanol to water. These results appear to indicate that van der Waals interactions gained from optimal packing are at least as important in stabilizing the native state of proteins as hydrophobic transfer effects.     

Replacement of staphylococcal nuclease hydrophobic core residues with those from thermophilic homologues indicates packing is improved in some thermostable proteins

The importance of tight hydrophobic core packing in stabilizing proteins found in thermophilic organisms has been vigorously disputed. Here, portions of the cores found in three thermophilic homologues were transplanted into the core of staphylococcal nuclease, a protein of modest stability. Packing of the core was evaluated by comparing interaction energy of the three mutants to the comprehensive mutant library built up previously at these same sites in staphylococcal nuclease. It was found that the interaction energy of one thermophilic sequence is extraordinarily favorable and the interaction energies of other two transplanted thermophilic sequences are good, comparable to the interaction energies of mutant cores based on cores found in mesophilic homologues. As expected when transferring just a portion of the core sequence, the mutant proteins were destabilized overall relative to wild-type staphylococcal nuclease. The overall conclusion is that improvement of packing interactions is a mechanism to confer stability employed in some proteins from thermophiles, but not all.     

Structure and dynamics of staphylococcal nuclease mutants as studied by fluorescence quenching techniques

Fluorescence techniques have been used to investigate the effect of mutations on the structure and dynamics of staphylococcal nuclease. An estimate of the accessibility to acrylamide of the enzyme's single tryptophan residue (Trp140) was obtained from the Stern-Volmer constant for fluorescence quenching. This was indicative of a partially buried tryptophan in the wild-type nuclease. Five single-site mutant nucleases (H124L, V66L, G88V, G79S and F76V) and one double mutant (V66L + G88V), with widely differing stabilities to denaturants, gave Stern-Volmer constants which were very similar to that of their parent enzyme. Studies of the temperature- and viscosity-dependence of quenching suggest that access by acrylamide to Trp140 is limited by diffusion rather than by protein structural fluctuations. Lifetime-resolved fluorescence anisotropy studies using steady-state instrumentation suggest that there is very little segmental motion of the Trp140; most of the anisotropy therefore decays due to protein rotation in the solution. Rotational correlation times for several nuclease mutants have been determined and these are very similar to that of the native nuclease. Thus it appears that these substitutions in the primary amino acid sequence, which have significant effects on the stability of the folded proteins, do not cause a significant change in the protein structure or dynamics around Trp140.     

Higher-order packing interactions in triple and quadruple mutants of staphylococcal nuclease

Sixty-four triple and 32 quadruple mutants were constructed in the core of staphylococcal nuclease. This is the first time that a large number of multiple mutants with all possible variations and all possible lower-order mutants has been systematically constructed in any protein core. Stabilities were determined by solvent denaturation. The energetic effects of these multiple mutants have been analyzed in combination with the stability data from the component single and double mutants. It was found that most of the stability changes in triple and quadruple mutants cannot be correctly predicted from stability effects of component single mutants. However, if the interaction energy between pairs of side chains in the component double mutants is taken into account, correct stability prediction can be made for most triple and quadruple mutants. The data further show that while packing interactions unique to triple and quadruple mutations do occur, they are of much less energetic significance than interactions between pairs of residues. The results presented here show that the packing of a protein interior can be closely approximated in most cases as a series of short-range, nearest-neighbor interactions. This has profound implications for rational protein design and structure prediction.     

Proteins with simplified hydrophobic cores compared to other packing mutants

Efforts to design proteins with greatly reduced sequence diversity have often resulted in proteins with so-called molten globule properties. Substitutions were made at six neighboring sites in the major hydrophobic core of staphylococcal nuclease to create variants with all leucine, all isoleucine or all valine at these sites. The mutant proteins with simplified cores constructed here are quite unstable and have poorly packed cores, attested to by interaction energies. Eight related mutants with greater sequence diversity were also constructed. Comparison to these mutants and 159 other permutations of these 3 aliphatic side chains at these same 6 sites previously constructed shows that the simplified cores are not unusual in their stabilities or interaction energies. Further, crystal structures of the two mutants with the worst packing, as measured by interaction energies, showed no unusual disorder in the core. Therefore, reduction of sequence diversity is not necessarily incompatible with a single stable native structure. Other factors must also contribute to previous protein design failures.     

Structural and energetic differences between insertions and substitutions in staphylococcal nuclease.

In a previous study, the small protein staphylococcal nuclease was shown to readily accommodate single alanine and glycine insertions, with average losses in stability comparable to substitutions at the same sites (PROT. 7:299-305, 1990). To more fully explore this unexpected adaptability to changes in residue spacing, 2 double amino acid insertions (alanyl-glycine, glycyl-glycine) and 3 additional single amino acid insertions with dissimilar side chains (proline, leucine, and glutamine) were constructed at 10 of the sites previously studied. At 8 of these sites, the type of amino acid side chain on the inserted residue significantly influenced the stability of the mutant protein. However, at 9 of the 10 sites, the double insertions were found to be no more destabilizing than the single alanine or glycine insertions. In contrast, double substitution mutations of staphylococcal nuclease, which replace two adjacent residues with alanine, do not show this striking degree of non-additivity. A comparison of the effects of single glutamine and single glycine insertions with alanyl-glycine insertions indicates that insertion of alanine into the peptide backbone is, on average, less destabilizing than appending the equivalent atoms onto the side chain of a glycine insertion. To explain their very different energetic effects, we propose that, unlike most substitutions, the inserted residue(s) must induce lateral displacements of the polypeptide chain, forcing the folded conformation away from that of wild type. The resulting obligatory shifts in the positioning of residues flanking the insertion generate a large number of degrees of freedom around which the mutant structure can relax.(ABSTRACT TRUNCATED AT 250 WORDS)     

Complex interactions at the helix-helix interface stabilize the glycophorin A transmembrane dimer

To explore the residue interactions in the glycophorin A dimerization motif, an alanine scan double mutant analysis at the helix-helix interface was carried out. These data reveal a combination of additive and coupled effects. The majority of the double mutants are found to be equally or slightly more stable than would be predicted by the sum of the energetic cost of the single-point mutants. The proximity of the mutated sites is not related to the presence of coupling between those sites. Previous studies reveal that a single face of the glycophorin A monomer contains a specific glycine-containing motif (GxxxG) that is thought to be a driving force for the association of transmembrane helices. Double mutant cycles suggest that the relationship of the GxxxG motif to the remainder of the helix-helix interface is complex. Sequences containing mutations that abolish the GxxxG motif retain an ability to dimerize, while a sequence containing a GxxxG motif appears unable to form dimers. The energetic effects of weakly coupled and additive double mutants can be explained by changes in van der Waals interactions at the dimer interface. These results emphasize the fact that the sequence context of the dimer interface modulates the strength of the glycophorin A GxxxG-mediated transmembrane dimerization reaction.     

Structural mechanisms of the degenerate sequence recognition by Bse634I restriction endonuclease

Restriction endonuclease Bse634I recognizes and cleaves the degenerate DNA sequence 5'-R/CCGGY-3' (R stands for A or G; Y for T or C, '/' indicates a cleavage position). Here, we report the crystal structures of the Bse634I R226A mutant complexed with cognate oligoduplexes containing ACCGGT and GCCGGC sites, respectively. In the crystal, all potential H-bond donor and acceptor atoms on the base edges of the conserved CCGG core are engaged in the interactions with Bse634I amino acid residues located on the α6 helix. In contrast, direct contacts between the protein and outer base pairs are limited to van der Waals contact between the purine nucleobase and Pro203 residue in the major groove and a single H-bond between the O2 atom of the outer pyrimidine and the side chain of the Asn73 residue in the minor groove. Structural data coupled with biochemical experiments suggest that both van der Waals interactions and indirect readout contribute to the discrimination of the degenerate base pair by Bse634I. Structure comparison between related enzymes Bse634I (R/CCGGY), NgoMIV (G/CCGGC) and SgrAI (CR/CCGGYG) reveals how different specificities are achieved within a conserved structural core.     

Kinetic folding and unfolding of staphylococcal nuclease and its six mutants studied by stopped-flow circular dichroism

Kinetics of refolding and unfolding of staphylococcal nuclease and its six mutants, each carrying single or double amino acid substitutions, are studied by stopped-flow circular dichroism measurements. A transient kinetic intermediate formed within 10 ms after refolding starts possesses a substantial part of the N-domain core β-structure, whereas helices are formed at the later stages. The structure of the kinetic intermediate is less organized than the structure that is known to be formed by a nuclease 1-136 fragment. Only the refolding kinetics are affected by the mutations in all the mutants except two in which the mutations have changed the native structure. From this result and also from the locations of the mutation sites, the major N-terminal domain of the nuclease in the transition state of folding has a structure nearly identical to the native one. On the other hand, the minor C-terminal domain has previously been shown to be still disorganized in the transition state. The effects of the amino acid substitutions on the stability of the native and the transition states are in good agreement with the changes in the hydration free energy, expected for the corresponding amino acid replacements in the unfolded polypeptide. Since side chains of all the mutated residues are not accessible to solvent in the native structure, the result suggests that it is the unfolded state that is mainly affected by the mutations. © 1995 Wiley-Liss, Inc.     

Structural response to mutation at a protein-protein interface

We have crystallised three mutants of the barnase-barstar complex in which interactions across the interface have been deleted by simultaneous mutation of both residues involved in the interaction. Each mutant deletes a different type of interaction at the interface: the first complex bnHis102-->Ala-bsTyr29-->Phe (bn, barnase; bs, barstar), deletes a van der Waals packing interaction; the second complex, bnLys27-->Ala-bsThr42-->Ala, deletes a hydrogen bond; the third, bnLys27-->Ala-bsAsp35-->Ala, deletes a long-range charge-charge interaction. The contribution of each of these side-chains to the stability of the complex is known; the coupling energy between the deleted side-chains is also known. Despite each of the double mutants being significantly destabilised compared with the wild-type, the effects of mutation are local. Only small movements in the main-chain surrounding the sites of mutation and some larger movements of neighbouring side-chains are observed in the mutant complexes. The exact response to mutation is context-dependent and for the same mutant can vary depending upon the environment within the crystal. In some double mutant complexes, interfacial pockets, which are accessible to bulk solvent are formed, whereas interfacial cavities which are isolated from bulk solvent, are formed in others. In all double mutants, water molecules fill the created pockets and cavities. These water molecules mimic the deleted side-chains by occupying positions close to the non-carbon atoms of truncated side-chains and re-making many hydrogen bonds made by the truncated side-chains in the wild-type. It remains extremely difficult, however, to correlate energetic and structural responses to mutation because of unknown changes in entropy and entropy-enthalpy compensation.     

Residue cluster additivity of thermodynamic stability in the hydrophobic core of mesophile vs. hyperthermophile rubredoxins

The branched sidechain residues 24 and 33 in the hydrophobic core of rubredoxin differ between the Clostridium pasteurianum (Cp) and Pyrococcus furiosus (Pf) sequences. Their X-ray structures indicate that these two sidechains are in van der Waals contact with each other, while neither appears to significantly interact with the other nonconserved residues. The simultaneous interchange of residues 24 and 33 between the Cp and Pf rubredoxin sequences yield a complementary pair of hybrid proteins for which the sum of their thermodynamic stabilities equals that of the parental rubredoxins. The 1.2 kcal/mol change arising from this two residues interchange accounts for 21% of the differential thermodynamic stability between the mesophile and hyperthermophile proteins. The additional interchange of the sole nonconserved aromatic residue in the hydrophobic core yields a 0.78 kcal/mol deviation from thermodynamic additivity.     

Core side-chain packing and backbone conformation in Lpp-56 coiled-coil mutants

Native proteins exhibit precise geometric packing of atoms in their hydrophobic interiors. Nonetheless, controversy remains about the role of core side-chain packing in specifying and stabilizing the folded structures of proteins. Here we investigate the role of core packing in determining the conformation and stability of the Lpp-56 trimerization domain. The X-ray crystal structures of Lpp-56 mutants with alanine substitutions at two and four interior core positions reveal trimeric coiled coils in which the twist of individual helices and the helix-helix spacing vary significantly to achieve the most favored superhelical packing arrangement. Introduction of each alanine "layer" into the hydrophobic core destabilizes the superhelix by 1.4 kcal mol(-1). Although the methyl groups of the alanine residues pack at their optimum van der Waals contacts in the coiled-coil trimer, they provide a smaller component of hydrophobic interactions than bulky hydrophobic side-chains to the thermodynamic stability. Thus, specific side-chain packing in the hydrophobic core of coiled coils are important determinants of protein main-chain conformation and stability.     

The constraints protein-protein interactions place on sequence divergence

IL4-BP, the extracellular binding domain of the IL-4 receptor alpha chain, contains a high-affinity binding epitope for IL-4 (K(D) 150 pM). Previous results on the crystal structure of the IL-4/IL4-BP complex and on the functional epitope of IL-4 suggested that this contact comprises a mosaic of two binding clusters. The present mutational analysis of IL4-BP supports this view and demonstrates that the energetically most important group is the receptor carboxylate group of D72 forming an ion pair with IL-4 R88 in cluster II. The second main receptor determinant is the hydroxyl group of Y183 forming a hydrogen bond with IL-4 E9 in cluster I. The latter is engaged in additional hydrogen bonds with Y13 and also in van der Waals contacts with Y127. Receptor residue D72 as well as Y183 are each surrounded by a shell of hydrophobic groups from residues that upon mutation lead to smaller decreases in binding affinity. Analysis of IL4-BP double mutants showed that receptor side-chains within one cluster but not those of different clusters cooperate. Interaction analysis of IL-4 and IL4-BP single mutants also revealed additivity in binding of side-chains between clusters and cooperativity within each cluster I or II.These results show that the high-affinity IL-4/IL4-BP contact is constituted by two independent binding units, each containing a central polar or charged side-chain surrounded by hydrophobic groups (avocado cluster).     

Comparison of straight chain and cyclic unnatural amino acids embedded in the core of staphylococcal nuclease.

We have determined by X-ray crystallography the structures of several variants of staphylococcal nuclease with long flexible straight chain and equivalent length cyclic unnatural amino acid side chains embedded in the protein core. The terminal atoms in the straight side chains are not well defined by the observed electron density even though they remain buried within the protein interior. We have previously observed this behavior and have suggested that it may arise from the addition of side-chain vibrational and oscillational motions with each bond as a side chain grows away from the relatively rigid protein main chain and/or the population of multiple rotamers (Wynn R, Harkins P, Richards FM. Fox RO. 1996. Mobile unnatural amino acid side chains in the core of staphylococcal nuclease. Protein Sci 5:1026-1031). Reduction of the number of degrees of freedom by cyclization of a side chain would be expected to constrain these motions. These side chains are in fact well defined in the structures described here. Over-packing of the protein core results in a 1.0 A shift of helix 1 away from the site of mutation. Additionally, we have determined the structure of a side chain containing a single hydrogen to fluorine atom replacement on a methyl group. A fluorine atom is intermediate in size between methyl group and a hydrogen atom. The fluorine atom is observed in a single position indicating it does not rotate like methyl hydrogen atoms. This change also causes subtle differences in the packing interactions.     

Effect of mutations involving charged residues on the stability of staphylococcal nuclease: a continuum electrostatics study

A continuum electrostatics model is used to calculate the relative stabilities of 117 mutants of staphylococcal nuclease (SNase) involving the mutation of a charged residue to an uncharged residue. The calculations are based on the crystallographic structure of the wild-type protein and attempt to take implicitly into account the effect of the mutations in the denatured state by assuming a linear relationship between the free energy changes caused by the mutation in the native and denatured states. A good correlation (linear correlation coefficient of approximately 0.8) is found with published experimental relative stabilities of these mutants. The results suggest that in the case of SNase (i) charged residues contribute to the stability of the native state mainly through electrostatic interactions, and (ii) native-like electrostatic interactions may persist in the denatured state. The continuum electrostatics method is only moderately sensitive to model parameters and leads to quasi-predictive results for the relative mutant stabilities (error of 2-3 kJ mol(-1) or of the order of k(B)T), except for mutants in which a charged residue is mutated to glycine.     

Atoms-in-molecules study of the genetically encoded amino acids. III. Bond and atomic properties and their correlations with experiment including mutation-induced changes in protein stability and genetic coding

This article presents a study of the molecular charge distributions of the genetically encoded amino acids (AA), one that builds on the previous determination of their equilibrium geometries and the demonstrated transferability of their common geometrical parameters. The properties of the charge distributions are characterized and given quantitative expression in terms of the bond and atomic properties determined within the quantum theory of atoms-in-molecules (QTAIM) that defines atoms and bonds in terms of the observable charge density. The properties so defined are demonstrated to be remarkably transferable, a reflection of the underlying transferability of the charge distributions of the main chain and other groups common to the AA. The use of the atomic properties in obtaining an understanding of the biological functions of the AA, whether free or bound in a polypeptide, is demonstrated by the excellent statistical correlations they yield with experimental physicochemical properties. A property of the AA side chains of particular importance is the charge separation index (CSI), a quantity previously defined as the sum of the magnitudes of the atomic charges and which measures the degree of separation of positive and negative charges in the side chain of interest. The CSI values provide a correlation with the measured free energies of transfer of capped side chain analogues, from the vapor phase to aqueous solution, yielding a linear regression equation with r2 = 0.94. The atomic volume is defined by the van der Waals isodensity surface and it, together with the CSI, which accounts for the electrostriction of the solvent, yield a linear regression (r2 = 0.98) with the measured partial molar volumes of the AAs. The changes in free energies of transfer from octanol to water upon interchanging 153 pairs of AAs and from cyclohexane to water upon interchanging 190 pairs of AAs, were modeled using only three calculated parameters (representing electrostatic and volume contributions) yielding linear regressions with r2 values of 0.78 and 0.89, respectively. These results are a prelude to the single-site mutation-induced changes in the stabilities of two typical proteins: ubiquitin and staphylococcal nuclease. Strong quadratic correlations (r2 approximately 0.9) were obtained between DeltaCSI upon mutation and each of the two terms DeltaDeltaH and TDeltaDeltaS taken from recent and accurate differential scanning calorimetry experiments on ubiquitin. When the two terms are summed to yield DeltaDeltaG, the quadratic terms nearly cancel, and the result is a simple linear fit between DeltaDeltaG and DeltaCSI with r2 = 0.88. As another example, the change in the stability of staphylococcal nuclease upon mutation has been fitted linearly (r2 = 0.83) to the sum of a DeltaCSI term and a term representing the change in the van der Waals volume of the side chains upon mutation. The suggested correlation of the polarity of the side chain with the second letter of the AA triplet genetic codon is given concrete expression in a classification of the side chains in terms of their CSI values and their group dipole moments. For example, all amino acids with a pyrimidine base as their second letter in mRNA possess side-chain CSI < or = 2.8 (with the exception of Cys), whereas all those with CSI > 2.8 possess an purine base. The article concludes with two proposals for measuring and predicting molecular complementarity: van der Waals complementarity expressed in terms of the van der Waals isodensity surface and Lewis complementarity expressed in terms of the local charge concentrations and depletions defined by the topology of the Laplacian of the electron density. A display of the experimentally accessible Laplacian distribution for a folded protein would offer a clear picture of the operation of the "stereochemical code" proposed as the determinant in the folding process.     

Theoretical studies of the response of a protein structure to cavity-creating mutations

We have investigated the response of a protein structure to cavity-creating mutations by molecular dynamics (MD) simulations for the wild-type and the five mutants of phage T4 lysozyme. Essential dynamics (ED) analysis and the methods for calculating different components of local interaction energies are used to examine the structural and energetic characteristics associated with the mutations. In agreement with the x-ray results, it is found that the structural changes due to the replacements of a bulky side chain such as Leu or Phe with Ala within the hydrophobic core can be characterized as slight adjustments rather than substantial reorganization of the protein. The relative stability of different mutant structures can be related with the extent of structural readjustments in response to the mutation. The destabilization of the mutant Leu-->Ala proteins relative to the wild-type is closely related with the loss of van der Waals contacts due to the cavity-creating mutations.