Average assignment method for predicting the stability of protein mutants

Prediction of protein stability upon amino acid substitutions is an important problem in molecular biology and it will be helpful for designing stable mutants. In this work, we have analyzed the stability of protein mutants using three different data sets of 1791, 1396, and 2204 mutants, respectively, for thermal stability (DeltaTm), free energy change due to thermal (DeltaDeltaG), and denaturant denaturations (DeltaDeltaGH2O), obtained from the ProTherm database. We have classified the mutants into 380 possible substitutions and assigned the stability of each mutant using the information obtained with similar type of mutations. We observed that this assignment could distinguish the stabilizing and destabilizing mutants to an accuracy of 70-80% at different measures of stability. Further, we have classified the mutants based on secondary structure and solvent accessibility (ASA) and observed that the classification significantly improved the accuracy of prediction. The classification of mutants based on helix, strand, and coil distinguished the stabilizing/destabilizing mutants at an average accuracy of 82% and the correlation is 0.56; information about the location of residues at the interior, partially buried, and surface regions of a protein correctly identified the stabilizing/destabilizing residues at an average accuracy of 81% and the correlation is 0.59. The nine subclassifications based on three secondary structures and solvent accessibilities improved the accuracy of assigning stabilizing/destabilizing mutants to an accuracy of 84-89% for the three data sets. Further, the present method is able to predict the free energy change (DeltaDeltaG) upon mutations within a deviation of 0.64 kcal/mol. We suggest that this method could be used for predicting the stability of protein mutants.     

Determination of protein secondary structure and solvent accessibility using site-directed fluorescence labeling. Studies of T4 lysozyme using the fluorescent probe monobromobimane

We report an investigation of how much protein structural information could be obtained using a site-directed fluorescence labeling (SDFL) strategy. In our experiments, we used 21 consecutive single-cysteine substitution mutants in T4 lysozyme (residues T115-K135), located in a helix-turn-helix motif. The mutants were labeled with the fluorescent probe monobromobimane and subjected to an array of fluorescence measurements. Thermal stability measurements show that introduction of the label is substantially perturbing only when it is located at buried residue sites. At buried sites (solvent surface accessibility of <40 A(2)), the destabilizations are between 3 and 5.5 kcal/mol, whereas at more exposed sites, DeltaDeltaG values of < or = 1.5 kcal/mol are obtained. Of all the fluorescence parameters that were explored (excitation lambda(max), emission lambda(max), fluorescence lifetime, quantum yield, and steady-state anisotropy), the emission lambda(max) and the steady-state anisotropy values most accurately reflect the solvent surface accessibility at each site as calculated from the crystal structure of cysteine-less T4 lysozyme. The parameters we identify allow the classification of each site as buried, partially buried, or exposed. We find that the variations in these parameters as a function of residue number reflect the sequence-specific secondary structure, the determination of which is a key step for modeling a protein of unknown structure.     

Revisiting “reverse hydrophobic effect”: Applicable only to coil mutations at the surface

In a seminal paper, Pakula and Sauer (Nature, 1990, 344, 363–364) demonstrated that the increase in side‐chain hydrophobicity has a reverse relationship with protein stability. We have addressed this problem with several examples of mutants that span at different locations in protein structure based on secondary structure and solvent accessibility. We confirmed that the stability change upon single coil mutation at exposed region is reversely correlated with hydrophobicity with a single exception. In addition, we found the existence of such relationship in partially buried coil mutants. The stability of exposed helical mutants is governed by conformational properties. In buried and partially buried helical and strand mutants properties reflecting hydrophobicity have direct relationship with stability, whereas an opposite relationship was obtained with entropy and flexibility. The structural analysis of partially buried/exposed mutants showed that the surrounding residues are important for the stability change upon mutation. These results provide insights to understand the general behavior for the stability of proteins upon amino acid substitutions. © 2009 Wiley Periodicals, Inc. Biopolymers 91: 591–599, 2009. This article was originally published online as an accepted preprint. The “Published Online” date corresponds to the preprint version. You can request a copy of the preprint by emailing the Biopolymers editorial office at biopolymers@wiley.com     

Contribution of the hydrophobic effect to globular protein stability.

The decrease in conformational stability, delta(delta G), has been measured for 72 aliphatic side-chain mutants from four proteins in which a larger side-chain is replaced by a smaller side-chain so that steric effects are minimal. When these delta(delta G) values are corrected to the same accessibility, namely 100% buried, then the following -delta(delta G) values per -CH2- group (in kcal/mol) are obtained: Ile----Val (1.26), Ala (1.26), Gly (1.26); Leu----Ala (1.16), Gly (1.21); Val----Ala (1.23), Gly (1.53). The average of these values is 1.27(+/- 0.07) kcal/mol. The 72 individual values range from 0 to 2.4 kcal/mol with an average value of 1.27(+/- 0.51) (standard deviation) kcal/mol. When the delta Gtr values from n-octanol to water are corrected for the difference in volume between the solutes and the solvents, the average value for the same substitutions is 1.25(+/- 0.05) kcal/mol. This suggests that proteins gain 1.3(+/- 0.5) kcal/mol in stability for each -CH2- group buried in folding, and, furthermore, that the volume corrected delta Gtr values for n-octanol for the amino acid side-chains provide good estimates of the contribution of the hydrophobic effect to globular protein stability.     

Fluorous effect in proteins: de novo design and characterization of a four-alpha-helix bundle protein containing hexafluoroleucine

Several studies have demonstrated that proteins incorporating fluorinated analogues of hydrophobic amino acids such as leucine and valine into their hydrophobic cores exhibit increased stability toward thermal denaturation and unfolding by guanidinium chloride. However, estimates for the increase in the thermodynamic stability of a protein (DeltaDeltaG(unfold)) afforded by the substitution of a hydrophobic amino acid with its fluorinated analogue vary quite significantly. To address this, we have designed a peptide that adopts an antiparallel four-helix bundle structure in which the hydrophobic core is packed with leucine, and investigated the effects of substituting the central two layers of the core with L-5,5,5,5',5',5'-hexafluoroleucine (hFLeu). We find that DeltaDeltaG(unfold) is increased by 0.3 kcal/mol per hFLeu residue. This is in good agreement with the predicted increase in DeltaDeltaG(unfold) of 0.4 kcal/mol per residue arising from the increased hydrophobicity of the hFLeu side chain, which we determined experimentally from partitioning measurements on hFLeu and leucine. The increased stability of this fluorinated protein may therefore be ascribed to simple hydrophobic effects, rather than specific "fluorous" interactions between the hFLeu residues.     

Toward the active conformation of insulin: stereospecific modulation of a structural switch in the B chain

How insulin binds to the insulin receptor has long been a subject of speculation. Although the structure of the free hormone has been extensively characterized, a variety of evidence suggests that a conformational change occurs upon receptor binding. Here, we employ chiral mutagenesis, comparison of corresponding d and l amino acid substitutions, to investigate a possible switch in the B-chain. To investigate the interrelation of structure, function, and stability, isomeric analogs have been synthesized in which an invariant glycine in a beta-turn (Gly(B8)) is replaced by d- or l-Ser. The d substitution enhances stability (DeltaDeltaG(u) 0.9 kcal/mol) but impairs receptor binding by 100-fold; by contrast, the l substitution markedly impairs stability (DeltaDeltaG(u) -3.0 kcal/mol) with only 2-fold reduction in receptor binding. Although the isomeric structures each retain a native-like overall fold, the l-Ser(B8) analog exhibits fewer helix-related and long range nuclear Overhauser effects than does the d-Ser(B8) analog or native monomer. Evidence for enhanced conformational fluctuations in the unstable analog is provided by its attenuated CD spectrum. The inverse relationship between stereospecific stabilization and receptor binding strongly suggests that the B7-B10 beta-turn changes conformation on receptor binding.     

Contribution of cation-pi interactions to protein stability

Calculations predict that cation- interactions make an important contribution to protein stability. While there have been some attempts to experimentally measure strengths of cation-pi interactions using peptide model systems, much less experimental data are available for globular proteins. We have attempted to determine the magnitude of cation-pi interactions of Lys with aromatic amino acids in four different proteins (LIVBP, MBP, RBP, and Trx). In each case, Lys was replaced with Gln and Met. In a separate series of experiments, the aromatic amino acid in each cation-pi pair was replaced by Leu. Stabilities of wild-type (WT) and mutant proteins were characterized by both thermal and chemical denaturation. Gln and aromatic --> Leu mutants were consistently less stable than corresponding Met mutants, reflecting the nonisosteric nature of these substitutions. The strength of the cation-pi interaction was assessed by the value of the change in the free energy of unfolding [DeltaDeltaG(degrees) = DeltaG(degrees)(Met) - DeltaG(degrees)(WT)]. This ranged from +1.1 to -1.9 kcal/mol (average value -0.4 kcal/mol) at 298 K and +0.7 to -2.6 kcal/mol (average value -1.1 kcal/mol) at the Tm of each WT. It therefore appears that the strength of cation-pi interactions increases with temperature. In addition, the experimentally measured values are appreciably smaller in magnitude than calculated values with an average difference /DeltaG(degrees)expt - DeltaG(degrees)calc/av of 2.9 kcal/mol. At room temperature, the data indicate that cation-pi interactions are at best weakly stabilizing and in some cases are clearly destabilizing. However, at elevated temperatures, close to typical Tm's, cation-pi interactions are generally stabilizing.     

Are acidic and basic groups in buried proteins predicted to be ionized?

Ionizable residues play essential roles in proteins, modulating protein stability, fold and function. Asp, Glu, Arg, and Lys make up about a quarter of the residues in an average protein. Multi-conformation continuum electrostatic (MCCE) calculations were used to predict the ionization states of all acidic and basic residues in 490 proteins. Of all 36,192 ionizable residues, 93.5% were predicted to be ionized. Thirty-five percent have lost 4.08 kcal/mol solvation energy (DeltaDeltaG(rxn)) sufficient to shift a pK(a) by three pH units in the absence of other interactions and 17% have DeltaDeltaG(rxn) sufficient to shift pK(a) by five pH units. Overall 85% of these buried residues (DeltaDeltaG(rxn)>5DeltapK units) are ionized, including 92% of the Arg, 86% of the Asp, 77% of the Glu, and 75% of the Lys. Ion-pair interactions stabilize the ionization of both acids and bases. The backbone dipoles stabilize anions more than cations. The interactions with polar side-chains are also different for acids and bases. Asn and Gln stabilize all charges, Ser and Thr stabilize only acids while Tyr rarely stabilize Lys. Thus, hydroxyls are better hydrogen bond donors than acceptors. Buried ionized residues are more likely to be conserved than those on the surface. There are 3.95 residues buried per 100 residues in an average protein.     

Thermodynamic stability of the bacteriorhodopsin lattice as measured by lipid dilution

To determine the strength of noncovalent interactions that stabilize a membrane protein complex, we have developed an in vitro method for quantifying the dissociation of the bacteriorhodopsin (BR) lattice, a naturally occurring two-dimensional crystal. A lattice suspension was titrated with a short- and long-chain phosphatidylcholine mixture to dilute BR within the lipid bilayer. The fraction of BR in the lattice form as a function of added lipid was determined by visible circular dichroism spectroscopy and fit with a cooperative self-assembly model to obtain a critical concentration for lattice assembly. Critical concentration values of wild-type and mutant proteins were used to calculate the change in lattice stability upon mutation (DeltaDeltaG). By using this method, a series of mutant proteins was examined in which residues at the BR-BR interface were replaced with smaller amino acids, either Ala or Gly. Most of the mutant lattices were destabilized, with DeltaDeltaG values of 0.2-1.1 kcal/mol at 30 degrees C, consistent with favorable packing of apolar residues in the membrane. One mutant, I45A, was stabilized by approximately 1.0 kcal/mol, possibly due to increased lipid entropy. The DeltaDeltaG values agreed well with previous in vivo measurements, except in the case of I45A. The ability to measure the change in stability of mutant protein complexes in a lipid bilayer may provide a means of determining the contributions of specific protein-protein and protein-lipid interactions to membrane protein structure.     

Role of proline residues in human lysozyme stability: a scanning calorimetric study combined with X-ray structure analysis of proline mutants.

It has been shown that protein stability can be modulated from site-directed mutations that affect the entropy of protein unfolding [Matthews, B. W., Nicholson, H., & Becktel, W. J. (1987) Proc. Natl. Acad. Sci. U.S.A. 84, 6663-6667]. However, the effect of a specific amino acid replacement on stability highly depends on the location of the mutation site and its environment in the protein structure [Yutani, K., Hayashi, S., Sugisaki, Y., & Ogasahara, K. (1991) Proteins Struct., Funct., Genet. 9, 90-98). To clarify the role of specific proline residues in the thermostability of human lysozyme (h-lysozyme), a series of proline mutants were investigated by means of scanning calorimetry and high-resolution X-ray crystallography. The thermodynamic properties of the mutant and wild-type h-lysozymes are compared and discussed on the basis of their three-dimensional structure. h-Lysozyme contains two proline residues at positions 71 and 103. The Pro71----Gly substitution was found to destabilize h-lysozyme by decreasing the entropic contribution of unfolding by about 2 kcal/mol at 68.8 degrees C. This is consistent with the theoretical expectations for such a substitution. However, the same substitution at position 103 (Pro103----Gly) does not affect h-lysozyme stability, and the thermodynamic properties of the P71G/P103G and P71G mutants are essentially the same. Pro71 which is conserved among lysozymes from other species, appears to be important for stability, whereas Pro103, which is not conserved, does not. These differences are explained in terms of residue accessibility to the solvent and crystallographic B-factor, which reflects the amino acid mobility.(ABSTRACT TRUNCATED AT 250 WORDS)     

Role of the minor energetic determinants of chicken egg white lysozyme (HEWL) to the stability of the HEWL.antibody scFv-10 complex

Seven of the 13 non-glycine contact amino acids in the hen (chicken) egg white lysozyme (HEWL) epitope for antibody Fab-10 each contribute < or =0.3 kcal/mol to the change in free energy (DeltaDeltaG(D)) from wild type (WT) when replaced by alanine (nullspots), and three others each give (0.7 < DeltaDeltaG(D) < or = 1. 0) kcal/mol (warm spots) (Rajpal et al. Protein Sci 1998;7:1868-1874). The low DeltaDeltaG(D) values introduced by alanine mutations present an opportunity to explore accurately their cumulative effects, as the sum of the combined DeltaDeltaG(D) values is not so large as to destabilize the complex beyond the range of accurate measurement. Substitution of six of the seven null spot residues by alanine leads to a cumulative DeltaDeltaG(D) = 2.25 +/- 0.04 kcal/mol, whereas the sum of the six individual changes is only -0.36 +/- 0.32 kcal/mol. The triple warm spot mutation generates a DeltaDeltaG(D) = 5.11 +/- 0.06 kcal/mol versus DeltaDeltaG(D) = 2.52 +/- 0.22 kcal/mol for the sum of the three individuals. The non-additivity in the individual DeltaDeltaG(D) values for the alanine mutations may indicate that these residues provide a conformationally stabilizing effect on the hot spot residues, each of which exhibits DeltaDeltaG(D) > 4.0 kcal/mol on alanine substitution.     

A new scale for side-chain contribution to protein stability based on the empirical stability analysis of mutant proteins

The hydrophobicity scales for amino acid side chains based on the transfer Gibbs energy (DeltaG(trans)) of amino acids from non-aqueous phases to water have been widely used to estimate the contribution of buried side chains to the conformational stability of proteins. In this paper, we propose a new scale for the side-chain contribution to protein stability, which is derived from data on protein denaturation experiments using systematic and comprehensive mutant proteins. In the experiments, the contribution of some physical properties were quantitatively determined as parameters in a unique equation representing the stability change (DeltaDeltaG) of mutant proteins as a function of the structural changes due to the mutations. These parameters are able conveniently to provide a scale for the side-chain contribution to protein stability. This new scale also has the advantage over the previously reported hydrophobicity scales of residues with the contributions of hydrogen bonds or secondary structural propensity. It may find practical application in algorithms for the prediction of protein structures.     

Support Vector Machine-based classification of protein folds using the structural properties of amino acid residues and amino acid residue pairs

MOTIVATION: Fold recognition is a key step in the protein structure discovery process, especially when traditional sequence comparison methods fail to yield convincing structural homologies. Although many methods have been developed for protein fold recognition, their accuracies remain low. This can be attributed to insufficient exploitation of fold discriminatory features. RESULTS: We have developed a new method for protein fold recognition using structural information of amino acid residues and amino acid residue pairs. Since protein fold recognition can be treated as a protein fold classification problem, we have developed a Support Vector Machine (SVM) based classifier approach that uses secondary structural state and solvent accessibility state frequencies of amino acids and amino acid pairs as feature vectors. Among the individual properties examined secondary structural state frequencies of amino acids gave an overall accuracy of 65.2% for fold discrimination, which is better than the accuracy by any method reported so far in the literature. Combination of secondary structural state frequencies with solvent accessibility state frequencies of amino acids and amino acid pairs further improved the fold discrimination accuracy to more than 70%, which is approximately 8% higher than the best available method. In this study we have also tested, for the first time, an all-together multi-class method known as Crammer and Singer method for protein fold classification. Our studies reveal that the three multi-class classification methods, namely one versus all, one versus one and Crammer and Singer method, yield similar predictions. AVAILABILITY: Dataset and stand-alone program are available upon request.     

Electrostatic contributions to T4 lysozyme stability: solvent-exposed charges versus semi-buried salt bridges

We carried our Poisson-Boltzmann (PB) calculations for the effects of charge reversal at five exposed sites (K16E, R119E, K135E, K147E, and R154E) and charge neutralization and proton titration of the H31-D70 semi-buried salt bridge on the stability of T4 lysozyme. Instead of the widely used solvent-exclusion (SE) surface, we used the van der Waals (vdW) surface as the boundary between the protein and solvent dielectrics (a protocol established in our earlier study on charge mutations in barnase). By including residual charge-charge interactions in the unfolded state, the five charge reversal mutations were found to have DeltaDeltaG(unfold) from -1.6 to 1.3 kcal/mol. This indicates that the variable effects of charge reversal observed by Matthews and co-workers are not unexpected. The H31N, D70N, and H31N/D70N mutations were found to destabilize the protein by 2.9, 1.3, and 1.6 kcal/mol, and the pK(a) values of H31 and D70 were shifted to 9.4 and 0.6, respectively. These results are in good accord with experimental data of Dahlquist and co-workers. In contrast, if the SE surface were used, the H31N/D70N mutant would be more stable than the wild-type protein by 1.3 kcal/mol. From these and additional results for 27 charge mutations on five other proteins, we conclude that 1) the popular view that electrostatic interactions are generally destabilizing may have been based on overestimated desolvation cost as a result of using the SE surface as the dielectric boundary; and 2) while solvent-exposed charges may not reliably contribute to protein stability, semi-buried salt bridges can provide significant stabilization.     

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.     

Calorimetric examination of high-affinity Src SH2 domain-tyrosyl phosphopeptide binding: dissection of the phosphopeptide sequence specificity and coupling energetics

SH2 domains are protein modules which interact with specific tyrosine phosphorylated sequences in target proteins. The SH2 domain of the Src kinase binds with high affinity to a tyrosine phosphorylated peptide containing the amino acids Glu, Glu, and Ile (EEI) at the positions +1, +2, and +3 C-terminal to the phosphotyrosine, respectively. To investigate the degree of selectivity of the Src SH2 domain for each amino acid of the EEI motif, the binding thermodynamics of a panel of substitutions at the +1 (Gln, Asp, Ala, Gly), +2 (Gln, Asp, Ala, Gly), and +3 (Leu, Val, Ala, Gly) positions were examined using titration microcalorimetry. It was revealed that the Src SH2 domain is insensitive (DeltaDeltaG degrees 相似文献   

SUPREX (Stability of Unpurified Proteins from Rates of H/D Exchange) analysis of the thermodynamics of synergistic anion binding by ferric-binding protein (FbpA), a bacterial transferrin

SUPREX (stability of unpurified proteins from rates of H/D exchange) is a H/D exchange- and matrix-assisted laser desorption/ionization (MALDI)-based technique for characterizing the equilibrium unfolding/refolding properties of proteins and protein-ligand complexes. Here, we describe the application of SUPREX to the thermodynamic analysis of synergistic anion binding to iron-loaded ferric-binding protein (Fe(3+)FbpA-X, X = synergistic anion). The in vivo function of FbpA is to transport unchelated Fe(3+) across the periplasmic space of certain Gram-negative bacteria, a process that requires simultaneous binding of a synergistic anion. Our results indicate that Fe(3+)FbpA-X is not a so-called "ideal" protein system for SUPREX analyses because it does not exhibit two-state folding properties and it does not exhibit EX2 H/D exchange behavior. However, despite these nonideal properties of the Fe(3+)FbpA-X protein-folding/unfolding reaction, we demonstrate that the SUPREX technique is still amenable to the quantitative thermodynamic analysis of synergistic anion binding to Fe(3+)FbpA. As part of this work, the SUPREX technique was used to evaluate the DeltaDeltaG(f) values of four synergistic anion-containing complexes of Fe(3+)FbpA (i.e., Fe(3+)FbpA-PO(4), Fe(3+)FbpA-citrate, Fe(3+)FbpA-AsO(4), and Fe(3+)FbpA-SO(4)). The DeltaDeltaG(f) value obtained for Fe(3+)FbpA-citrate relative to Fe(3+)FbpA-PO(4) (1.45 +/- 0.44 kcal/mol), is in good agreement with that reported previously (1.98 kcal/mol). The value obtained for Fe(3+)FbpA-AsO(4) (0.58 +/- 0.45 kcal/mol) was also consistent with that reported previously (0.68 kcal/mol), but the measurement error is very close to the magnitude of the value. This work (i) demonstrates the utility of the SUPREX method for studying anion binding by FbpA, (ii) provides the first evaluation of a DeltaDeltaG(f) value for Fe(3+)FbpA-SO(4), -1.43 +/- 0.17 kcal/mol, and (iii) helps substantiate our hypothesis that the synergistic anion plays a role in controlling the lability of iron bound to FbpA in the transport process.     

Effects of tryptophan microenvironment, soluble domain, and vesicle size on the thermodynamics of membrane protein folding: lessons from the transmembrane protein OmpA

Refolding curves of the integral membrane protein outer membrane protein A (OmpA) were measured to determine the conformational stabilities of this model system for membrane protein folding. Wild-type OmpA exhibits a free energy of unfolding (DeltaG degrees H2O) of 10.5 kcal/mol. Mutants, containing a single tryptophan residue at the native positions 7, 15, 57, 102, or 143, are less stable than wild-type OmpA, with DeltaG degrees H2O values of 6.7, 4.8, 2.4, 4.7, and 2.8 kcal/mol, respectively. The trend observed here is discussed in terms of noncovalent interactions, including aromatic interactions and hydrogen bonding. The effect of the soluble tail on the conformational stability of the transmembrane domain of OmpA was also investigated via truncated single-Trp mutants; DeltaG degrees H2O values for four of the five truncated mutants are greater by >2.7 kcal/mol relative to the full-length versions, suggesting that the absence of the soluble domain may destabilize the unfolded transmembrane domain. Finally, dynamic light scattering experiments were performed to measure the effects of urea and protein on vesicle size and stability. Urea concentrations greater than 1 M cause an increase in vesicle size, and these diameters are unaltered in the presence of protein. These dynamic light scattering results complement the fluorescence studies and illustrate the important effects of vesicle size on protein conformational stability.     

Context-dependent protein stabilization by methionine-to-leucine substitution shown in T4 lysozyme.

The substitution of methionines with leucines within the interior of a protein is expected to increase stability both because of a more favorable solvent transfer term as well as the reduced entropic cost of holding a leucine side chain in a defined position. Together, these two terms are expected to contribute about 1.4 kcal/mol to protein stability for each Met --> Leu substitution when fully buried. At the same time, this expected beneficial effect may be offset by steric factors due to differences in the shape of leucine and methionine. To investigate the interplay between these factors, all methionines in T4 lysozyme except at the amino-terminus were individually replaced with leucine. Of these mutants, M106L and M120L have stabilities 0.5 kcal/mol higher than wild-type T4 lysozyme, while M6L is significantly destabilized (-2.8 kcal/mol). M102L, described previously, is also destabilized (-0.9 kcal/mol). Based on this limited sample it appears that methionine-to-leucine substitutions can increase protein stability but only in a situation where the methionine side chain is fully or partially buried, yet allows the introduction of the leucine without concomitant steric interference. The variants, together with methionine-to-lysine substitutions at the same sites, follow the general pattern that substitutions at rigid, internal sites tend to be most destabilizing, whereas replacements at more solvent-exposed sites are better tolerated.     

Nine hydrophobic side chains are key determinants of the thermodynamic stability and oligomerization status of tumour suppressor p53 tetramerization domain.

The contribution of almost each amino acid side chain to the thermodynamic stability of the tetramerization domain (residues 326-353) of human p53 has been quantitated using 25 mutants with single-residue truncations to alanine (or glycine). Truncation of either Leu344 or Leu348 buried at the tetramer interface, but not of any other residue, led to the formation of dimers of moderate stability (8-9 kcal/mol of dimer) instead of tetramers. One-third of the substitutions were moderately destabilizing (<3.9 kcal/mol of tetramer). Truncations of Arg333, Asn345 or Glu349 involved in intermonomer hydrogen bonds, Ala347 at the tetramer interface or Thr329 were more destabilizing (4.1-5.7 kcal/mol). Strongly destabilizing (8.8- 11.7 kcal/mol) substitutions included those of Met340 at the tetramer interface and Phe328, Arg337 and Phe338 involved peripherally in the hydrophobic core. Truncation of any of the three residues involved centrally in the hydrophobic core of each primary dimer either prevented folding (Ile332) or allowed folding only at high protein concentration or low temperature (Leu330 and Phe341). Nine hydrophobic residues per monomer constitute critical determinants for the stability and oligomerization status of this p53 domain.