Estimating the contribution of engineered surface electrostatic interactions to protein stability by using double-mutant cycles

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.     

Variable contributions of tyrosine residues to the structural and spectroscopic properties of the factor for inversion stimulation

The diverse roles of tyrosine residues in proteins may be attributed to their dual hydrophobic and polar nature, which can result in hydrophobic and ring stacking interactions, as well as hydrogen bonding. The small homodimeric DNA binding protein, factor for inversion stimulation (FIS), contains four tyrosine residues located at positions 38, 51, 69, and 95, each involved in specific intra- or intermolecular interactions. To investigate their contributions to the stability, flexibility, and spectroscopic properties of FIS, each one was independently mutated to phenylalanine. Equilibrium denaturation experiments show that Tyr95 and Tyr51 stabilize FIS by about 2 and 1 kcal/mol, respectively, as a result of their involvement in a hydrogen bond-salt bridge network. In contrast, Tyr38 destabilizes FIS by about 1 kcal/mol due to the placement of a hydroxyl group in a hydrophobic environment. The stability of FIS was not altered when the solvent-exposed Tyr69 was mutated. Limited proteolysis with trypsin and V8 proteases was used to monitor the flexibility of the C-terminus (residues 71-98) and the dimer core (residues 26-70), respectively. The results for Y95F and Y51F FIS revealed a different proteolytic susceptibility of the dimer core compared to the C-terminus, suggesting an increased flexibility of the latter. DNA binding affinity of the various FIS mutants was only modestly affected and correlated inversely with the C-terminal flexibility probed by trypsin proteolysis. Deconvolution of the fluorescence contribution of each mutant revealed that it varies in intensity and direction for each tyrosine in WT FIS, highlighting the role of specific interactions and the local environment in determining the fluorescence of tyrosine residues. The significant changes in stability, flexibility, and signals observed for the Y51F and Y95F mutations are attributed to their coupled participation in the hydrogen bond-salt bridge network. These results highlight the importance of tyrosine hydrogen-bonding and packing interactions for the stability of FIS and demonstrate the varying roles that tyrosine residues can play on the structural and spectroscopic properties of even small proteins.     

A helix propensity scale based on experimental studies of peptides and proteins.

The average globular protein contains 30% alpha-helix, the most common type of secondary structure. Some amino acids occur more frequently in alpha-helices than others; this tendency is known as helix propensity. Here we derive a helix propensity scale for solvent-exposed residues in the middle positions of alpha-helices. The scale is based on measurements of helix propensity in 11 systems, including both proteins and peptides. Alanine has the highest helix propensity, and, excluding proline, glycine has the lowest, approximately 1 kcal/mol less favorable than alanine. Based on our analysis, the helix propensities of the amino acids are as follows (kcal/mol): Ala = 0, Leu = 0.21, Arg = 0.21, Met = 0.24, Lys = 0.26, Gln = 0.39, Glu = 0.40, Ile = 0.41, Trp = 0.49, Ser = 0.50, Tyr = 0. 53, Phe = 0.54, Val = 0.61, His = 0.61, Asn = 0.65, Thr = 0.66, Cys = 0.68, Asp = 0.69, and Gly = 1.     

Contributions of aromatic pairs to the folding and stability of long-lived human γD-crystallin

Human γD-crystallin (HγD-Crys) is a highly stable protein that remains folded in the eye lens for the majority of an individual's lifetime. HγD-Crys exhibits two homologous crystallin domains, each containing two Greek key motifs with eight β-strands. Six aromatic pairs (four Tyr/Tyr, one Tyr/Phe and one Phe/Phe) are present in the β-hairpin sequences of the Greek keys. Ultraviolet damage to the aromatic residues in lens crystallins may contribute to the genesis of cataract. Mutant proteins with these aromatic residues substituted with alanines were constructed and expressed in E. coli. All mutant proteins except F115A and F117A had lower thermal stability than the WT protein. In equilibrium experiments in guanidine hydrochloride (GuHCl), all mutant proteins had lower thermodynamic stability than the WT protein. N-terminal domain (N-td) substitutions shifted the N-td transition to lower GuHCl concentration, but the C-terminal domain (C-td) transition remained unaffected. C-td substitutions led to a more cooperative unfolding/refolding process, with both the N-td and C-td transitions shifted to lower GuHCl concentration. The aromatic pairs conserved for each Greek key motif (Greek key pairs) had larger contributions to both thermal stability and thermodynamic stability than the other pairs. Aromatic-aromatic interaction was estimated as 1.5-2.0 kcal/mol. In kinetic experiments, N-td substitutions accelerated the early phase of unfolding, while C-td substitutions accelerated the late phase, suggesting independent domain unfolding. Only substitutions of the second Greek key pair of each crystallin domain slowed refolding. The second Greek keys may provide nucleation sites during the folding of the double-Greek-key crystallin domains.     

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.     

A survey of aspartate-phenylalanine and glutamate-phenylalanine interactions in the protein data bank: searching for anion-π pairs

Protein structures are stabilized using noncovalent interactions. In addition to the traditional noncovalent interactions, newer types of interactions are thought to be present in proteins. One such interaction, an anion-π pair, in which the positively charged edge of an aromatic ring interacts with an anion, forming a favorable anion-quadrupole interaction, has been previously proposed [Jackson, M. R., et al. (2007) J. Phys. Chem. B111, 8242-8249]. To study the role of anion-π interactions in stabilizing protein structure, we analyzed pairwise interactions between phenylalanine (Phe) and the anionic amino acids, aspartate (Asp) and glutamate (Glu). Particular emphasis was focused on identification of Phe-Asp or -Glu pairs separated by less than 7 ? in the high-resolution, nonredundant Protein Data Bank. Simplifying Phe to benzene and Asp or Glu to formate molecules facilitated in silico analysis of the pairs. Kitaura-Morokuma energy calculations were performed on roughly 19000 benzene-formate pairs and the resulting energies analyzed as a function of distance and angle. Edgewise interactions typically produced strongly stabilizing interaction energies (-2 to -7.3 kcal/mol), while interactions involving the ring face resulted in weakly stabilizing to repulsive interaction energies. The strongest, most stabilizing interactions were identified as preferentially occurring in buried residues. Anion-π pairs are found throughout protein structures, in helices as well as β strands. Numerous pairs also had nearby cation-π interactions as well as potential π-π stacking. While more than 1000 structures did not contain an anion-π pair, the 3134 remaining structures contained approximately 2.6 anion-π pairs per protein, suggesting it is a reasonably common motif that could contribute to the overall structural stability of a protein.     

Role of hydrogen bond networks and dynamics in positive and negative cooperative stabilization of a protein

Cooperativity mediated through hydrogen bond networks in yeast iso-1-cytochrome c was studied using a thermodynamic triple mutant cycle. Three known stabilizing mutations, Asn 26 to His, Asn 52 to Ile, and Tyr 67 to Phe, were used to construct the triple mutant cycle. The side chain of His 26, a wild-type residue, forms two hydrogen bonds that bridge two substructures of the wild-type protein, and Tyr 67 and Asn 52 are part of an extensive buried hydrogen bond network. The stabilities of all variants in the triple mutant cycle were determined by guanidine hydrochloride denaturation methods and used to determine the pairwise, Delta(2)G(int), and triple interaction energies. His 26 and Ile 52 interact cooperatively (Delta(2)G(int) is 1-2 kcal/mol), whereas the two other pairs of mutations interact anticooperatively (Delta(2)G(int) is -0.5 to -1.5 kcal/mol). Previously reported structural data for iso-1-cytochrome c variants containing these mutations show that changes in the strength of the His 26 to Glu 44 hydrogen bond, apparently caused by changes in main chain dynamics, provide a mechanism for the long distance (His 26 to Phe 67 and His 26 to Ile 52) propagation of pairwise interaction energies. Opposing changes in the strength of the His 26 to Glu 44 hydrogen bond caused by the N52I and Y67F mutations generate a negative triple interaction energy (-0.9 +/-0.7 kcal/mol) that combined with cancellation of cooperative and anticooperative pairwise interactions produce apparent additivity for the stabilizing effects of the single mutations in the triple mutant variant.     

Probing weakly polar interactions in cytochrome c.

Theoretical, statistical, and model studies suggest that proteins are stabilized by weakly polar attractions between sulfur atoms and properly oriented aromatic rings. The two sulfur-containing amino acids, methionine and cysteine, occur frequently among functional alleles in random mutant libraries of Saccharomyces cerevisiae iso-1-cytochrome c genes at positions that form a weakly polar aromatic-aromatic interaction, the wild-type protein. To determine if a weakly polar sulfur-aromatic interaction replaced the aromatic-aromatic interaction, the structure and stability of two variants were examined. Phenylalanine 10, which interacts with tyrosine 97, was replaced by methionine and cysteine. The cysteine was modified to form the methionine and cysteine analog, S-methyl cysteine (CysSMe). Proton NMR studies indicate that changing Phe 10 to Met or CysSMe affects only local structure and that the structures of sulfur-containing variants are nearly identical. Analysis of chemical shifts and nuclear Overhauser effect data indicates that both sulfur-containing side chains are in position to form a weakly polar interaction with Tyr 97. The F10M and F10CSMe variants are 2-3 kcal mol-1 less stable than iso-1-cytochrome c at 300 K. Comparison of the stabilities of the F10M and F10CSMe variants allows evaluation of the potential weakly polar interaction between the additional sulfur atom of F10CSMe and the aromatic moiety of Tyr 97. The F10CSMe;C102T variant is 0.7 +/- 0.3 kcal mol-1 more stable than the F10M;C102T protein. The increased stability is explained by the difference in hydrophobicity of the sulfur-containing side chains. We conclude that any weakly polar interaction between the additional sulfur and the aromatic ring is too weak to detect or is masked by destabilizing contributions to the free energy of denaturation.     

Surface salt bridges,double-mutant cycles,and protein stability: an experimental and computational analysis of the interaction of the Asp 23 side chain with the N-terminus of the N-terminal domain of the ribosomal protein l9

Experimental and theoretical double-mutant cycles have been used to investigate a salt bridge in the N-terminal domain of the protein L9. Aspartic acid 23 is the only acidic residue involved in a well-defined pairwise interaction, namely, a partially solvent-exposed salt bridge with the protonated N-terminus of the protein. Mutations were studied in which Asp 23 was substituted by alanine, asparagine, and nitrile alanine. Interactions with the N-terminus were probed by comparisons between proteins with a protonated and acetylated N-terminus. The mutants were all folded, and the structures were unchanged from wild type as judged by CD and 2-D NMR. The coupling free energy between the N-terminus and the side chain of Asp 23 measured through double-mutant cycle analysis was favorable and ranged from -0.7 to -1.7 kcal mol(-)(1), depending upon the set of mutants used. This relatively large coupling free energy for a surface salt bridge likely arises from geometric factors that reduce the entropy loss associated with salt-bridge formation and from structural relaxation in the mutants. Coupling free energies computed with continuum electrostatic calculations agreed well with the experimental values when full account was taken of all potential interactions, particularly those involving Asp 23 and the acetylated N-terminus as well as interactions with solvent. The measured and calculated coupling free energy decreased only slightly when the salt concentration was increased from 100 to 750 mM NaCl. The calculations suggest that the coupling free energy between D23 and the N-terminus measured through the experimental double-mutant cycle analysis is significantly smaller than the actual interaction free energy between the groups in the wild-type structure because of the inapplicability of assumptions frequently used to interpret double-mutant cycles.     

Pressure-dependent ionization of Tyr 9 in glutathione S-transferase A1-1: contribution of the C-terminal helix to a "soft" active site.

The glutathione S-transferase (GST) isozyme A1-1 contains at its active site a catalytic tyrosine, Tyr9, which hydrogen bonds to, and stabilizes, the thiolate form of glutathione, GS-. In the substrate-free GST A1-1, the Tyr 9 has an unusually low pKa, approximately 8.2, for which the ionization to tyrosinate is monitored conveniently by UV and fluorescence spectroscopy in the tryptophan-free mutant, W21F. In addition, a short alpha-helix, residues 208-222, provides part of the GSH and hydrophobic ligand binding sites, and the helix becomes "disordered" in the absence of ligands. Here, hydrostatic pressure has been used to probe the conformational dynamics of the C-terminal helix, which are apparently linked to Tyr 9 ionization. The extent of ionization of Tyr 9 at pH 7.6 is increased dramatically at low pressures (p1/2 = 0.52 kbar), based on fluorescence titration of Tyr 9. The mutant protein W21F:Y9F exhibits no changes in tyrosine fluorescence up to 1.2 kbar; pressure specifically ionizes Tyr 9. The volume change, delta V, for the pressure-dependent ionization of Tyr 9 at pH 7.6, 19 degrees C, was -33 +/- 3 mL/mol. In contrast, N-acetyl tyrosine exhibits a delta V for deprotonation of -11 +/- 1 mL/mol, beginning from the same extent of initial ionization, pH 9.5. The pressure-dependent ionization is completely reversible for both Tyr 9 and N-acetyl tyrosine. Addition of S-methyl GSH converted the "soft" active site to a noncompressible site that exhibited negligible pressure-dependent ionization of Tyr 9 below 0.8 kbar. In addition, Phe 220 forms part of an "aromatic cluster" with Tyr 9 and Phe 10, and interactions among these residues were hypothesized to control the order of the C-terminal helix. The amino acid substitutions F220Y, F2201, and F220L afford proteins that undergo pressure-dependent ionization of Tyr 9 with delta V values of 31 +/- 2 mL/mol, 43 +/- 3 mL/mol, and 29 +/- 2 mL/mol, respectively. The p1/2 values for Tyr 9 ionization were 0.61 kbar, 0.41 kbar, and 0.46 kbar for F220Y, F220I, and F220L, respectively. Together, the results suggest that the C-terminal helix is conformationally heterogeneous in the absence of ligands. The conformations differ little in free energy, but they are significantly different in volume, and mutations at Phe 220 control the conformational distribution.     

Tyr74 is essential for the formation, stability and function of Plasmodium falciparum triosephosphate isomerase dimer

Plasmodium falciparum triosephosphate isomerase (PfTIM) is known to be functional only as a homodimer. Although many studies have shown that the interface Cys13 plays a major role in the stability of the dimer, a few reports have demonstrated that structurally conserved Tyr74 may be essential for the stability of PfTIM dimer. To understand the role of Tyr74, we have performed molecular dynamics (MD) simulations of monomeric and dimeric PfTIM mutated to glycine and cysteine at position 74. Simulations of the monomer revealed that mutant Tyr74Gly does not produce changes in folding and stability of the monomer. Interestingly, comparison of the flexibility of Tyr74 in the monomer and dimer revealed that this residue possesses an intrinsic restricted mobility, indicating that Tyr74 is an anchor residue required for homodimerization. Tyr74 also appears to play an important role in binding by facilitating the disorder-to-order transitions of loops 1 and 3, which allows Cys13 to form favorable interactions with loop 3 and Lys12 to be locked in a favorable position for catalysis. High-temperature MD simulations of the wild-type and Tyr74Gly PfTIM dimers showed that the aromatic moiety of Tyr74 is necessary to preserve the geometry and native contacts between loops 1 and 3 at the interface of the dimer. Disulfide cross-linking between mutant Tyr74Cys and Cys13 further revealed that Tyr74 stabilizes the geometry of loop 1 (which contains the catalytic residue Lys12) and the interactions between loops 1 and 3 via aromatic-aromatic interactions with residues Phe69, Tyr101, and Phe102. Principal component analysis showed that Tyr74 is also necessary to preserve the collective motions in the dimer that contribute to the catalytic efficiency of PfTIM dimer. We conclude that Tyr74 not only plays a role in the stability of the dimer, but also participates in the dimerization process and collective motions via coupled disorder-to-order transitions of intrinsically disordered regions, necessary for efficiency in the catalytic function of PfTIM.     

Contributions of engineered surface salt bridges to the stability of T4 lysozyme determined by directed mutagenesis.

Six designed mutants of T4 lysozyme were created in an attempt to create putative salt bridges on the surface of the protein. The first three of the mutants, T115E (Thr 115 to Glu), Q123E, and N144E, were designed to introduce a new charged side chain close to one or more existing charged groups of the opposite sign on the surface of the protein. In each of these cases the putative electrostatic interactions introduced by the mutation include possible salt bridges between residues within consecutive turns of an alpha-helix. Effects of the mutations ranged from no change in stability to a 1.5 degrees C (0.5 kcal/mol) increase in melting temperature. In two cases, secondary (double) mutants were constructed as controls in which the charge partner was removed from the primary mutant structure. These controls proteins indicate that the contributions to stability from each of the engineered salt bridges is very small (about 0.1-0.25 kcal/mol in 0.15 M KCl). The structures of the three primary mutants were determined by X-ray crystallography and shown to be essentially the same as the wild-type structure except at the site of the mutation. Although the introduced charges in the T115E and Q123E structures are within 3-5 A of their intended partner, the introduced side chains and their intended partners were observed to be quite mobile. It has been shown that the salt bridge between His 31 and Asp 70 in T4 lysozyme stabilizes the protein by 3-5 kcal/mol [Anderson, D. E., Becktel, W. J., & Dahlquist, F. W. (1990) Biochemistry 29, 2403-2408]. To test the effectiveness of His...Asp interactions in general, three additional double mutants, K60H/L13D, K83H/A112D, and S90H/Q122D, were created in order to introduce histidine-aspartate charge pairs on the surface of the protein. Each of these mutants destabilizes the protein by 1-3 kcal/mol in 0.15 M KCl at pH values from 2 to 6.5. The X-ray crystallographic structure of the mutant K83H/A112D has been determined and shows that there are backbone conformational changes of 0.3-0.6 A extending over several residues. The introduction of the histidine and aspartate presumably introduces strain into the folded protein that destabilizes this variant. It is concluded that pairs of oppositely charged residues that are on the surface of a protein and have freedom to adopt different conformations do not tend to come together to form structurally localized salt bridges. Rather, such residues tend to remain mobile, interact weakly if at all, and do not contribute significantly to protein stability. It is argued that the entropic cost of localizing a pair of solvent-exposed charged groups on the surface of a protein largely offsets the interaction energy expected from the formation of a defined salt bridge. There are examples of strong salt bridges in proteins, but such interactions require that the folding of the protein provides the requisite driving energy to hold the interacting partners in the correct rigid alignment.     

Probing aromatic, hydrophobic, and steric effects on the self-assembly of an amyloid-β fragment peptide

Aromatic amino acids have been shown to promote self-assembly of amyloid peptides, although the basis for this amyloid-inducing behavior is not understood. We adopted the amyloid-β 16-22 peptide (Aβ(16-22), Ac-KLVFFAE-NH(2)) as a model to study the role of aromatic amino acids in peptide self-assembly. Aβ(16-22) contains two consecutive Phe residues (19 and 20) in which Phe 19 side chains form interstrand contacts in fibrils while Phe 20 side chains interact with the side chain of Va l18. The kinetic and thermodynamic effect of varying the hydrophobicity and aromaticity at positions 19 and 20 by mutation with Ala, Tyr, cyclohexylalanine (Cha), and pentafluorophenylalanine (F(5)-Phe) (order of hydrophobicity is Ala < Tyr < Phe < F(5)-Phe < Cha) was characterized. Ala and Tyr position 19 variants failed to undergo fibril formation at the peptide concentrations studied, but Cha and F(5)-Phe variants self-assembled at dramatically enhanced rates relative to wild-type. Cha mutation was thermodynamically stabilizing at position 20 (ΔΔG = -0.2 kcal mol(-1) relative to wild-type) and destabilizing at position 19 (ΔΔG = +0.2 kcal mol(-1)). Conversely, F(5)-Phe mutations were strongly stabilizing at both positions (ΔΔG = -1.3 kcal mol(-1) at 19, ΔΔG = -0.9 kcal mol(-1) at 20). The double Cha and F(5)-Phe mutants showed that the thermodynamic effects were additive (ΔΔG = 0 kcal mol(-1) for Cha 19,20 and -2.1 kcal mol(-1) for F(5)-Phe 19,20). These results indicate that sequence hydrophobicity alone does not dictate amyloid potential, but that aromatic, hydrophobic, and steric considerations collectively influence fibril formation.     

Superadditive and subadditive effects of "hot spot" mutations within the interfaces of placental ribonuclease inhibitor with angiogenin and ribonuclease A.

Previous single-site mutagenesis studies on the complexes of ribonuclease inhibitor (RI) with angiogenin (Ang) and RNase A suggested that in both cases a substantial fraction of the binding energy is concentrated within one small part of the crystallographically observed interface, involving RI residues 434-438. Such energetic "hot spots" are common in protein-protein complexes, but their physical meaning is generally unclear. Here we have investigated this question by examining the detailed interactions within the RI.ligand hot spots and the extent to which they function independently. The effects of Phe versus Ala substitutions show that the key residue Tyr434 interacts with both ligands primarily through its phenyl ring; for Tyr437, the OH group forms the important contacts with RNase A, whereas the phenyl group interacts with Ang. Kinetic characterization of complexes containing multiple substitutions reveals striking, but distinctive, cooperativity in the interactions of RI with the two ligands. The losses in binding energy for the RNase complex associated with replacements of Tyr434 and Asp435, and Tyr434 and Tyr437, are markedly less than additive (i.e., by 2.4 and 1.3 kcal/mol, respectively). In contrast, the energetic effects of the 434 and 435, and 434 and 437, substitution pairs on binding of Ang are fully additive and 2.5 kcal/mol beyond additive, respectively. Superadditivities (0.9-2.4 kcal/mol) are also observed for several multisite replacements involving these inhibitor residues and two Ang residues, Arg5 and Lys40, from this part of the interface. Consequently, the decreases in binding energy for some triple-variant complexes are as large as 8.5-10.1 kcal/mol (compared to a total DeltaG of -21.0 kcal/mol for the wild-type complex). Potential explanations for these functional couplings, many of which occur over distances of >13 A and are not mediated by direct or triangulated contacts, are proposed. These findings show that the basis for the generation of hot spots can be complex, and that these sites can assume significantly more (as with Ang) or less (as with RNase) importance than indicated from the effects of single-site mutations.     

Contribution of hydrophobic interactions to protein stability

Our goal was to gain a better understanding of the contribution of hydrophobic interactions to protein stability. We measured the change in conformational stability, Δ(ΔG), for hydrophobic mutants of four proteins: villin headpiece subdomain (VHP) with 36 residues, a surface protein from Borrelia burgdorferi (VlsE) with 341 residues, and two proteins previously studied in our laboratory, ribonucleases Sa and T1. We compared our results with those of previous studies and reached the following conclusions: (1) Hydrophobic interactions contribute less to the stability of a small protein, VHP (0.6 ± 0.3 kcal/mol per -CH₂- group), than to the stability of a large protein, VlsE (1.6 ± 0.3 kcal/mol per -CH₂- group). (2) Hydrophobic interactions make the major contribution to the stability of VHP (40 kcal/mol) and the major contributors are (in kilocalories per mole) Phe18 (3.9), Met13 (3.1), Phe7 (2.9), Phe11 (2.7), and Leu21 (2.7). (3) Based on the Δ(ΔG) values for 148 hydrophobic mutants in 13 proteins, burying a -CH₂- group on folding contributes, on average, 1.1 ± 0.5 kcal/mol to protein stability. (4) The experimental Δ(ΔG) values for aliphatic side chains (Ala, Val, Ile, and Leu) are in good agreement with their ΔG_tr values from water to cyclohexane. (5) For 22 proteins with 36 to 534 residues, hydrophobic interactions contribute 60 ± 4% and hydrogen bonds contribute 40 ± 4% to protein stability. (6) Conformational entropy contributes about 2.4 kcal/mol per residue to protein instability. The globular conformation of proteins is stabilized predominantly by hydrophobic interactions.     

Influence of cation-pi interactions in different folding types of membrane proteins

Cation-pi interactions play an important role to the stability of protein structures. In this work, we analyze the influence of cation-pi interactions in three-dimensional structures of membrane proteins. We found that transmembrane strand (TMS) proteins have more number of cation-pi interactions than transmembrane helical (TMH) proteins. In TMH proteins, both the positively charged residues Lys and Arg equally experience favorable cation-pi interactions whereas in TMS proteins, Arg is more likely than Lys to be in such interactions. There is no relationship between number of cation-pi interactions and number of residues in TMH proteins whereas a good correlation was observed in TMS proteins. The average cation-pi interaction energy for TMH proteins is -16 kcal/mol and that for TMS proteins is -27 kcal/mol. The pair-wise cation-pi interaction energy between aromatic and positively charged residues showed that Lys-Trp energy is stronger in TMS proteins than TMH proteins; Arg-Phe, Arg-Tyr and Lys-Phe have higher energy in TMH proteins than TMS proteins. The decomposition of energies into electrostatic and van der Waals revealed that the contribution from electrostatic energy is twice as that from van der Waals energy in both TMH and TMS proteins. The results obtained in the present study would be helpful to understand the contribution of cation-pi interactions to the stability of membrane proteins.     

Binding and release of cholesterol in the Osh4 protein of yeast

Sterols have been shown experimentally to bind to the Osh4 protein (a homolog of the oxysterol binding proteins) of Saccharomyces cerevisiae within a binding tunnel, which consists of antiparallel beta-sheets that resemble a beta-barrel and three alpha-helices of the N-terminus. This and other Osh proteins are essential for intracellular transport of sterols and ultimately cell life. Molecular dynamics (MD) simulations are used to study the binding of cholesterol to Osh4 at the atomic level. The structure of the protein is stable during the course of all MD simulations and has little deviation from the experimental crystal structure. The conformational stability of cholesterol within the binding tunnel is aided in part by direct or water-mediated interactions between the 3-hydroxyl (3-OH) group of cholesterol and Trp(46), Gln(96), Tyr(97), Asn(165), and/or Gln(181) as well as dispersive interactions with Phe(42), Leu(24), Leu(39), Ile(167), and Ile(203). These residues along with other nonpolar residues in the binding tunnel and lid contribute nearly 75% to the total binding energy. The strongest and most populated interaction is between Gln(96) and 3-OH with a cholesterol/Gln(96) interaction energy of -4.5 +/- 1.0 kcal/mol. Phe(42) has a similar level of attraction to cholesterol with -4.1 +/- 0.3 kcal/mol. A MD simulation without the N-terminus lid that covers the binding tunnel resulted in similar binding conformations and binding energies when compared with simulations with the full-length protein. Steered MD was used to determine details of the mechanism used by Osh4 to release cholesterol to the cytoplasm. Phe(42), Gln(96), Asn(165), Gln(181), Pro(211), and Ile(206) are found to direct the cholesterol as it exits the binding tunnel as well as Lys(109). The mechanism of sterol release is conceptualized as a molecular ladder with the rungs being amino acids or water-mediated amino acids that interact with 3-OH.     

An essential residue in the flexible peptide linking the two idiosynchratic domains of bacterial tyrosyl-tRNA synthetases

Tyrosyl-tRNA synthetase (TyrRS) from Bacillus stearothermophilus comprises three sequential domains: an N-terminal catalytic domain, an alpha-helical domain with unknown function, and a C-terminal tRNA binding domain (residues 320-419). The properties of the polypeptide segment that links the alpha-helical and C-terminal domains, were analyzed by measuring the effects of sequence changes on the aminoacylation of tRNA(Tyr) with tyrosine. Mutations F323A (Phe323 into Ala), S324A, and G325A showed that the side chain of Phe323 was essential but not those of Ser324 and Gly325. Insertions of Gly residues between Leu322 and Phe323 and the point mutation L322P showed that the position and precise orientation of Phe323 relative to the alpha-helical domain were important. Insertions of Gly residues between Gly325 and Asp326 and deletion of residues 330-339 showed that the length and flexibility of the sequence downstream from Gly325 were unimportant but that this sequence could not be deleted. Mutations F323A, -L, -Y, and -W showed that the essential property of Phe323 was its aromaticity. The Phe323 side chain contributed to the stability of the initial complex between TyrRS and tRNA(Tyr) for 2.0 +/- 0.2 kcal x mol(-1) and to the stability of their transition state complex for 4.2 +/- 0.1 kcal x mol(-1), even though it is located far from the catalytic site. The results indicate that the disorder of the C-terminal domain in the crystals of TyrRS is due to the flexibility of the peptide that links it to the helical domain. They identified Phe323 as an essential residue for the recognition of tRNA(Tyr).     

Electrostatic stabilization in four-helix bundle proteins.

Charge substitutions generated by site-directed mutagenesis at the termini of adjacent anti-parallel alpha-helices in a four-helix bundle protein were used to determine a precise value for the contribution of indirect charge-charge interactions to overall protein stability, and to simulate the electrostatic effects of alpha-helix macrodipoles. Thermodynamic double mutant cycles were constructed to measure the interaction energy between such charges on adjacent anti-parallel helices in the four-helix bundle cytochrome b562 from Escherichia coli. Previously, theoretical calculations of helix macrodipole interactions using modeled four-helix bundle proteins have predicted values ranging over an order of magnitude from 0.2 to 2.5 kcal/mol. Our system represents the first experimental evidence for electrostatic interactions such as those between partial charges due to helix macrodipole charges. At the positions mutated, we have measured a favorable interaction energy of 0.6 kcal/mol between opposite charges simulating an anti-parallel helix pair. Pairs of negative or positive charges simulating a parallel orientation of helices produce an unfavorable interaction of similar magnitude. The interaction energies show a strong dependence upon ionic strength, consistent with an electrostatic effect. Indirect electrostatic contacts do appear to confer a limited stabilization upon the association of anti-parallel packing of helices, favoring this orientation by as much as 1 kcal/mol at 20 mM K phosphate.     

Tyrosine hydrogen bonds make a large contribution to protein stability

The aim of this study was to gain a better understanding of the contribution of hydrogen bonds by tyrosine -OH groups to protein stability. The amino acid sequences of RNases Sa and Sa3 are 69 % identical and each contains eight Tyr residues with seven at equivalent structural positions. We have measured the stability of the 16 tyrosine to phenylalanine mutants. For two equivalent mutants, the stability increases by 0.3 kcal/mol (RNase Sa Y30F) and 0.5 kcal/mol (RNase Sa3 Y33F) (1 kcal=4.184 kJ). For all of the other mutants, the stability decreases with the greatest decrease being 3.6 kcal/mol for RNase Sa Y52F. Seven of the 16 tyrosine residues form intramolecular hydrogen bonds and the average decrease in stability for these is 2.0(+/-1.0) kcal/mol. For the nine tyrosine residues that do not form intramolecular hydrogen bonds, the average decrease in stability is 0.4(+/-0.6) kcal/mol. Thus, most tyrosine -OH groups contribute favorably to protein stability even if they do not form intramolecular hydrogen bonds. Generally, the stability changes for equivalent positions in the two proteins are remarkably similar. Crystal structures were determined for two of the tyrosine to phenylalanine mutants of RNase Sa: Y80F (1.2 A), and Y86F (1.7 A). The structures are very similar to that of wild-type RNase Sa, and the hydrogen bonding partners of the tyrosine residues always form intermolecular hydrogen bonds to water in the mutants. These results provide further evidence that the hydrogen bonding and van der Waals interactions of polar groups in the tightly packed interior of folded proteins are more favorable than similar interactions with water in the unfolded protein, and that polar group burial makes a substantial contribution to protein stability.