Thermodynamic effects of proline introduction on protein stability

The amino acid Pro is more rigid than other naturally occurring amino acids and, in proteins, lacks an amide hydrogen. To understand the structural and thermodynamic effects of Pro substitutions, it was introduced at 13 different positions in four different proteins, leucine-isoleucine-valine binding protein, maltose binding protein, ribose binding protein, and thioredoxin. Three of the maltose binding protein mutants were characterized by X-ray crystallography to confirm that no structural changes had occurred upon mutation. In the remaining cases, fluorescence and CD spectroscopy were used to show the absence of structural change. Stabilities of wild type and mutant proteins were characterized by chemical denaturation at neutral pH and by differential scanning calorimetry as a function of pH. The mutants did not show enhanced stability with respect to chemical denaturation at room temperature. However, 6 of the 13 single mutants showed a small but significant increase in the free energy of thermal unfolding in the range of 0.3-2.4 kcal/mol, 2 mutants showed no change, and 5 were destabilized. In five of the six cases, the stabilization was because of reduced entropy of unfolding. However, the magnitude of the reduction in entropy of unfolding was typically several fold larger than the theoretical estimate of -4 cal K(-1) mol(-1) derived from the relative areas in the Ramachandran map accessible to Pro and Ala residues, respectively. Two double mutants were constructed. In both cases, the effects of the single mutations on the free energy of thermal unfolding were nonadditive.     

Correlation between mutational destabilization of phage T4 lysozyme and increased unfolding rates

The thermodynamics and kinetics of unfolding of 28 bacteriophage T4 lysozyme variants were compared by using urea gradient gel electrophoresis. The mutations studied cause a variety of sequence changes at different residues throughout the polypeptide chain and result in a wide range of thermodynamic stabilities. A striking relationship was observed between the thermodynamic and kinetic effects of the amino acid replacements: All the substitutions that destabilized the native protein by 2 kcal/mol or more also increased the rate of unfolding. The observed increases in unfolding rate corresponded to a decrease in the activation energy of unfolding (delta Gu) at least 35% as large as the decrease in thermodynamic stability (delta Gu). Thus, the destabilizing lesions bring the free energy of the native state closer to that of both the unfolded state and the transition state for folding and unfolding. Since a large fraction of the mutational destabilization is expressed between the transition state and the native conformation, the changes in folding energetics cannot be accounted for by effects on the unfolded state alone. The results also suggest that interactions throughout much of the folded structure are altered in the formation of the transition state during unfolding.     

Thermodynamic and structural studies of cavity formation in proteins suggest that loss of packing interactions rather than the hydrophobic effect dominates the observed energetics

The hydrophobic effect is widely believed to be an important determinant of protein stability. However, it is difficult to obtain unambiguous experimental estimates of the contribution of the hydrophobic driving force to the overall free energy of folding. Thermodynamic and structural studies of large to small substitutions in proteins are the most direct method of measuring this contribution. We have substituted the buried residue Phe8 in RNase S with alanine, methionine, and norleucine. Binding thermodynamics and structures were characterized by titration calorimetry and crystallography, respectively. The crystal structures of the RNase S F8A, F8M, and F8Nle mutants indicate that the protein tolerates the changes without any main chain adjustments. The correlation of structural and thermodynamic parameters associated with large to small substitutions was analyzed for nine mutants of RNase S as well as 32 additional cavity-containing mutants of T4 lysozyme, human lysozyme, and barnase. Such substitutions were typically found to result in negligible changes in DeltaC(p)() and positive values of both DeltaDeltaH degrees and DeltaDeltaS of folding. Enthalpic effects were dominant, and the sign of DeltaDeltaS is the opposite of that expected from the hydrophobic effect. Values of DeltaDeltaG degrees and DeltaDeltaH degrees correlated better with changes in packing parameters such as residue depth or occluded surface than with the change in accessible surface area upon folding. These results suggest that the loss of packing interactions rather than the hydrophobic effect is a dominant contributor to the observed energetics for large to small substitutions. Hence, estimates of the magnitude of the hydrophobic driving force derived from earlier mutational studies are likely to be significantly in excess of the actual value.     

Thermodynamics of 2'-ribose substitutions in UUCG tetraloops

The ribose 2'-hydroxyl group confers upon RNA many unique molecular properties. To better appreciate its contribution to structure and stability and to monitor how substitutions of the 2' hydroxyl can alter an RNA molecule, each loop pyrimidine ribonucleotide in the UUCG tetraloop was substituted with a nucleotide containing either a fluorine (2'-F), hydrogen (2'-H), amino (2'-NH2), or methoxy (2'-OCH3) group, in the context of both the C:G and G:C loop-closing base pair. The thermodynamic parameters of these tetraloop variants have been determined and NMR experiments used to monitor the structural changes resulting from the substitutions. The modified riboses are better tolerated in the G[UUCG]C tetraloop, which may be due to its increased loop flexibility relative to the C[UUCG]G loop. Even for these simple substitutions, the free-energy change reflects a complex interplay of hydrogen bonding, solvation effects, and intrinsic pucker preferences of the nucleotides.     

The importance of a distal hydrogen bonding group in stabilizing the transition state in subtilisin BPN'

Stabilization of an oxyanion transition state is important to catalysis of peptide bond hydrolysis in all proteases. For subtilisin BPN', a bacterial serine protease, structural data suggest that two hydrogen bonds stabilize the tetrahedral-like oxyanion intermediate: one from the main chain NH of Ser221 and another from the side chain NH2 of Asn155. Molecular dynamic studies (Rao, S., N., Singh, U., C. Bush, P. A., and Kollman, P. A. (1987) Nature 328, 551-554) have indicated the gamma-hydroxyl of Thr220 may be a third hydrogen bond donor even though it is 4A away in the static x-ray structure. We have probed the role of Thr220 by replacing it with serine, cysteine, valine, or alanine by site-directed mutagenesis. These substitutions were intended to alter the size and hydrogen bonding ability of residue 220. Removal of the gamma-hydroxyl group reduced the transition state stabilization energy (delta delta GT) by 1.8-2.1 kcal/mol depending upon the substitution. By comparison, removal of the gamma-methyl group in the Thr220 to serine mutation only decreased delta GT by 0.5 kcal/mol. The gamma-hydroxyl of Thr220 is most important for catalysis, not substrate binding, because virtually all of the effects were on kcat and not KM. The role of the Thr220 hydroxyl is functionally independent from the amide NH2 of Asn155 because the free energy effects of double alanine mutants at these two positions are additive. These data indicate that a distal hydrogen bond donor, namely the hydroxyl of Thr220, plays a functionally important role in stabilizing the oxyanion transition state in subtilisin which is independent of Asn155.     

Electrostatic interactions in the denatured state and in the transition state for protein folding: effects of denatured state interactions on the analysis of transition state structure

The development of electrostatic interactions during the folding of the N-terminal domain of the ribosomal protein L9 (NTL9) is investigated by pH-dependent rate equilibrium free energy relationships. We show that Asp8, among six acidic residues, is involved in non-native, electrostatic interactions with K12 in the transition state for folding as well as in the denatured state. The perturbed native state pK(a) of D8 (pK(a) = 3.0) appears to be maintained through non-native interactions in both the transition state and the denatured state. Mutational effects on the stability of the transition state for protein (un)folding are often analyzed in respect to change in ground states. Thus, the interpretation of transition state analysis critically depends on an understanding of mutational effects on both the native and denatured state. Increasing evidence for structurally biased denatured states under physiological conditions raises concerns about possible denatured state effects on folding studies. We show that the structural interpretation of transition state analysis can be altered dramatically by denatured state effects.     

Transfer RNA aminoacylation: identification of a critical ribose 2'-hydroxyl-base interaction.

To understand the relationship between tRNA architecture and specific aminoacylation by aminoacyl-tRNA synthetases, we performed kinetic assays of Escherichia coli tRNA(Pro) molecules containing single deoxynucleotide substitutions. We identified an important 2'-hydroxyl group at position U8 (of 22 positions probed). Chemical modification studies showed that this 2'-hydroxyl interacts with either the N1 or the exocyclic amine of G46 in a hydrogen bonding interaction that contributes 1.8 kcal/mol to the free energy of activation for aminoacylation. Molecular modeling of tRNA(Pro) supports the existence of this interaction. This is the first study to identify a specific ribose 2'-hydroxyl-base interaction in the core region of a tRNA molecule that makes a thermodynamically significant contribution to aminoacylation.     

Effects of base substitutions in an RNA hairpin from molecular dynamics and free energy simulations

Contributions of individual interactions in the GGCGCAAGCC hairpin containing a GCAA tetraloop were studied by computer simulations using base substitutions. The G in the first tetraloop position was replaced by inosine (I) or adenosine (A), and the G in the C-G basepair closing the tetraloop was replaced by I. These substitutions eliminate particular hydrogen bonds proposed in the nuclear magnetic resonance model of the GCAA tetraloop. Molecular dynamics simulations of the GCAA tetraloop in aqueous solvent displayed a well-defined hydrogen pattern between the first and last loop nucleotides (G and A) stabilized by a bridging water molecule. Substitution of G-->I in the basepair closing the tetraloop did not significantly influence the loop structure and dynamics. The ICAA loop maintained the overall structure, but displayed variation in the hydrogen-bond network within the tetraloop itself. Molecular dynamics simulations of the ACAA loop led to conformational heterogeneity of the resulting structures. Changes of hairpin formation free energy associated with substitutions of individual bases were calculated by the free energy perturbation method. The calculated decrease of the hairpin stability upon G-->I substitution in the C-G basepair closing the tetraloop was in good agreement with experimental thermodynamic data. Our theoretical estimates for G-->I and G-->A mutations located in the tetraloop suggest larger loop destabilization than corresponding experimental results. The extent of conformational sampling of the structures resulting from base substitutions and its impact on the calculated free energy was discussed.     

All-atom structural investigation of kinesin-microtubule complex constrained by high-quality cryo-electron-microscopy maps

In this study, we have performed a comprehensive structural investigation of three major biochemical states of a kinesin complexed with microtubule under the constraint of high-quality cryo-electron-microscopy (EM) maps. In addition to the ADP and ATP state which were captured by X-ray crystallography, we have also modeled the nucleotide-free or APO state for which no crystal structure is available. We have combined flexible fitting of EM maps with regular molecular dynamics simulations, hydrogen-bond analysis, and free energy calculation. Our APO-state models feature a subdomain rotation involving loop L2 and α6 helix of kinesin, and local structural changes in active site similar to a related motor protein, myosin. We have identified a list of hydrogen bonds involving key residues in the active site and the binding interface between kinesin and microtubule. Some of these hydrogen bonds may play an important role in coupling microtubule binding to ATPase activities in kinesin. We have validated our models by calculating the binding free energy between kinesin and microtubule, which quantitatively accounts for the observation of strong binding in the APO and ATP state and weak binding in the ADP state. This study will offer promising targets for future mutational and functional studies to investigate the mechanism of kinesin motors.     

Computational Modeling of Structurally Conserved Cancer Mutations in the RET and MET Kinases: The Impact on Protein Structure, Dynamics, and Stability

Structural and biochemical characterization of protein kinases that confer oncogene addiction and harbor a large number of disease-associated mutations, including RET and MET kinases, have provided insights into molecular mechanisms associated with the protein kinase activation in human cancer. In this article, structural modeling, molecular dynamics, and free energy simulations of a structurally conserved mutational hotspot, shared by M918T in RET and M1250T in MET kinases, are undertaken to quantify the molecular mechanism of activation and the functional role of cancer mutations in altering protein kinase structure, dynamics, and stability. The mechanistic basis of the activating RET and MET cancer mutations may be driven by an appreciable free energy destabilization of the inactive kinase state in the mutational forms. According to our results, the locally enhanced mobility of the cancer mutants and a higher conformational entropy are counterbalanced by a larger enthalpy loss and result in the decreased thermodynamic stability. The computed protein stability differences between the wild-type and cancer kinase mutants are consistent with circular dichroism spectroscopy and differential scanning calorimetry experiments. These results support the molecular mechanism of activation, which causes a detrimental imbalance in the dynamic equilibrium shifted toward the active form of the enzyme. Furthermore, computer simulations of the inhibitor binding with the oncogenic and drug-resistant RET mutations have also provided a plausible molecular rationale for the observed differences in the inhibition profiles, which is consistent with the experimental data. Finally, structural mapping of RET and MET cancer mutations and the computed protein stability changes suggest a similar mechanism of activation, whereby the cancer mutations which display the higher oncogenic activity tend to have the greatest destabilization effect on the inactive kinase structure.     

Stereoselective interaction with chiral phosphorothioates at the central DNA kink of the EcoRI endonuclease-GAATTC complex.

We have probed the contacts between EcoRI endonuclease and the central phosphate of its recognition site GAApTTC, using synthetic oligonucleotides containing single stereospecific Rp- or Sp-phosphorothioates (Ps). These substitutions produce subtle stereospecific effects on EcoRI endonuclease binding and cleavage. An Sp-Ps substitution in one strand of the DNA duplex improves binding free energy by -1.5 kcal/mol, whereas the Rp-Ps substitution has an unfavorable effect (+0.3 kcal/mol) on binding free energy. These effects derive principally from changes in the first order rate constants for dissociation of the enzyme-DNA complexes. The first order rate constants for strand scission are also affected, in that a strand containing Sp-Ps substitution is cleaved 2 to 3 times more rapidly than a strand containing a normal prochiral phosphate, whereas a strand containing Rp-Ps substitution is cleaved about 3 times slower than normal. As a result, single-strand substitutions produce pronounced asymmetry in the rates of cleavage of the two DNA strands, and this effect is exaggerated in an Rp,Sp-heteroduplex. Ethylation-interference footprinting indicates that none of the Ps substitutions cause any major change in contacts between endonuclease and DNA phosphates. When an Sp-Ps localizes P = O in the DNA major groove, a hydrogen-bonding interaction with the backbone amide-NH of Gly116 of the endonuclease is improved relative to that with a prochiral phosphate having intermediate P-O bond order and delocalized charge.     

Kinetic isotope effects reveal the presence of significant secondary structure in the transition state for the folding of the N-terminal domain of L9

Our present understanding of the nature of the transition state for protein folding depends predominantly on studies where individual side-chain contributions are mapped out by mutational analysis (phi value analysis). This approach, although extremely powerful, does not in general provide direct information about the formation of backbone hydrogen bonds. Here, we report the results of amide H/D isotope effect studies that probe the development of hydrogen bonded interactions in the transition state for the folding of a small alpha-beta protein, the N-terminal domain of L9. Replacement of amide protons by deuterons in a solvent of constant isotopic composition destabilized the domain, decreasing both its T(m) and Delta G(0) of unfolding. The folding rate also decreased. The parameter Phi(H/D), defined as the ratio of the effect of isotopic substitution upon the activation free energy to the equilibrium free energy was determined to be 0.6 in a D(2)O background and 0.75 in a H(2)O background, indicating that significant intraprotein hydrogen bond interactions are developed in the transition state for the folding of NTL9. The value is in remarkably good agreement with more traditional measures of the position of the transition state, which report on the relative burial of surface area. The results provide a picture of a compact folding transition state containing significant secondary structure. Indirect analysis argues that the bulk of the kinetic isotope effect arises from the beta-sheet-rich region of the protein, and suggests that the development of intraprotein hydrogen bonds in this region plays a critical role in the folding of NTL9.     

Mutational analysis of the complex of human RNase inhibitor and human eosinophil-derived neurotoxin (RNase 2)

RNase inhibitor (RI) binds diverse proteins in the pancreatic RNase superfamily with extremely high avidity. Previous studies showed that tight binding of RNase A and angiogenin (Ang) is achieved primarily through interactions of hot spot residues in the 434-460 C-terminal segment of RI with the enzymatic active site; Asp435 of RI forms key hydrogen bonds with the catalytic lysine in both complexes, whereas the other contacts are largely distinctive. Here we have investigated the structural basis for recognition of a third ligand, eosinophil-derived neurotoxin (EDN), by single-site and multisite mutagenesis. Surprisingly, Ala replacement of Asp435 decreases affinity for EDN only by 14-fold, as compared to the several hundred-fold decreases with RNase A and Ang, and individual mutations of three other hot spot residues-Tyr434, Tyr437, and Ser460-have essentially no effect. Ala substitutions of nine additional residues, selected by examining a computational model of the RI.EDN complex, also have no marked impact. Overall, the losses in affinity for the single-residue variants examined account for only approximately 25% of the free energy of binding for the complex. However, multisite mutagenesis of RI reveals strong superadditivity of mutational effects, indicating that part of this shortfall reflects negative cooperativity. Replacement of Tyr434 together with Asp435 or Tyr437 increases K(i) by 540- and 290-fold, respectively. Thus, the C-terminal region of RI again plays an important role in ligand recognition, although probably smaller than for binding RNase A and Ang. Simultaneous substitutions of three neighboring tryptophans (261, 263, and 318) on RI attenuate affinity even more dramatically (by 4900-fold), indicating that the interactions of this RI region also contribute a considerable amount of the binding energy for the EDN complex. These findings highlight the potential importance of cooperativity in protein-protein interactions and the consequent limitations of single-site mutagenesis for assessing interface energetics.     

Context-dependent mutations predominate in an engineered high-affinity single chain antibody fragment

A mutational analysis of the femtomolar-affinity anti-fluorescein antibody 4M5.3, compared to its wild-type progenitor, 4-4-20, indicates both context-dependent and -independent mutations are responsible for the 1800-fold affinity improvement. 4M5.3 was engineered from 4-4-20 by directed evolution and contains 14 mutations. The seven mutations identified as present in each of 10 final round affinity maturation clones were studied here. Affinities of the 4-4-20 single mutant addition and 4M5.3 single site reversion mutants were compared. These experiments identified four mutations, of these seven, that were context-dependent in their contribution to higher affinity. A simplified mutant containing only these seven mutations was created to analyze complete double mutant cycles of selected sets of mutations. Specific mutational sets studied included the ligand contact mutations, the heavy chain CDR3 mutations, the heavy chain CDR3 mutations plus the neighboring residue at site H108, and the early and late acquired mutations on the directed evolution pathway. The heavy chain CDR3 mutational set and the ligand-contacting mutations were shown to provide -1.4 and -2.0 kcal/mol, respectively, of the total -3.5 kcal/mol change in free energy of binding of the seven-site consensus mutant. The mutations acquired late in the directed evolution rounds provided much of the change in free energy without the earlier acquired mutations (-3.1 kcal/mol of the total -3.5 kcal/mol). Prior structural data and electrostatic calculations presented several hypotheses for the higher affinity contributions, some of which are supported by these mutational data.     

Dynamics and stability of GCAA tetraloops with 2-aminopurine and purine substitutions

The contributions of various interactions in the GGCGCAAGCC hairpin containing a GCAA tetraloop were studied by computer simulations using the substitutions of functional groups. The guanosine (G) in the first tetraloop position or in the C-G closing base pair was replaced by 2-aminopurine (AP), and the individual tetraloop's adenosines (A) were replaced by purine (PUR). These substitutions eliminated particular hydrogen bonds thought to stabilize the GCAA tetraloop. For each substitution, molecular dynamics (MD) simulations were carried out in an aqueous solution with sodium counterions, using the CHARMM27 force field. The MD simulations showed that the substitutions in the first (G-->AP) and the third (A-->PUR) position of the GCAA tetraloop did not significantly influence the conformation of the hairpin. A long-lived bridging water molecule observed in the GCAA loop was present in both modified loops. The substitutions made in the last loop position (A-->PUR) or in the C-G base pair closing the tetraloop (G-->AP) to some extent influenced the loop structure and dynamics. These loops did not display the long-lived bridging water molecules. When the second A in the GCAA loop was replaced by PUR, the first A in the loop was observed in the anti or in the syn orientation about the glycosyl bond. The G to AP substitution in C-G base pair led to a change of their arrangement from the Watson-Crick to wobble. The MD simulations of the hairpin with C-AP wobble closing base pair showed increased conformational dynamics of the hairpin. The changes of hairpin formation free energy associated with the substitutions of individual bases were calculated by the free energy perturbation method. Our theoretical estimates suggest a larger destabilization for the G to AP substitutions in GCAA loop than for the substitutions of individual A's by PUR, which is in accordance with experimental tendency. The calculations predicted a similar free energy change for G to AP substitutions in the GCAA tetraloop and in the C-G closing base pair.     

The effect of Asp-His-Ser/Thr-Trp tetrad on the thermostability of WD40-repeat proteins

We recently found that Asp-His-Ser/Thr-Trp hydrogen-bonded tetrads are widely and uniquely present in the WD40-repeat proteins. WDR5 protein is a seven WD40-repeat propeller with five such tetrads. To explore the effect of the tetrad on the structure and stability of WD40-repeat proteins, the wild-type WDR5 and its seven mutants involving the substitutions of tetrad residues have been isolated. The crystal structures of the wild-type WDR5 and its three WDR5 mutants have been determined by X-ray diffraction method. The mutations of the tetrad residues are found not to change the basic structural features. The denaturing profiles of the wild type and the seven mutants with the use of denaturant guanidine hydrochloride have been studied by circular dichroism spectroscopy to determine the folding free energies of these proteins. The folding free energies of the wild type and the S62A, S146A, S188A, D192E, W330F, W330Y, and D324E mutants are measured to be about -11.6, -2.7, -3.1, -2.9, -3.6, -7.1, -7.0, and -7.5 kcal/mol, respectively. These suggest that (1) the hydrogen bonds in these hydrogen bond networks are unusually strong; (2) each hydrogen-bonded tetrad provides over 12 kcal/mol stability to the protein; thus, the removal of any single tetrad would cause unfolding of the protein; (3) since there are five tetrads, the protein must be in a highly unstable state without the tetrads, which might be related to its biological functions.     

Substitutions in a flexible loop of horse liver alcohol dehydrogenase hinder the conformational change and unmask hydrogen transfer.

When horse liver alcohol dehydrogenase binds coenzyme, a rotation of about 10 degrees brings the catalytic domain closer to the coenzyme binding domain and closes the active site cleft. The conformational change requires that a flexible loop containing residues 293-298 in the coenzyme binding domain rearranges so that the coenzyme and some amino acid residues from the catalytic domain can be accommodated. The change appears to control the rate of dissociation of the coenzyme and to be necessary for installation of the proton relay system. In this study, directed mutagenesis produced the activated Gly293Ala/Pro295Thr enzyme. X-ray crystallography shows that the conformations of both free and complexed forms of the mutated enzyme and wild-type apoenzyme are very similar. Binding of NAD(+) and 2,2, 2-trifluoroethanol do not cause the conformational change, but the nicotinamide ribose moiety and alcohol are not in a fixed position. Although the Gly293Ala and Pro295Thr substitutions do not disturb the apoenzyme structure, molecular modeling shows that the new side chains cannot be accommodated in the closed native holoenzyme complex without steric alterations. The mutated enzyme may be active in the "open" conformation. The turnover numbers with ethanol and acetaldehyde increase 1.5- and 5.5-fold, respectively, and dissociation constants for coenzymes and other kinetic constants increase 40-2,000-fold compared to those of the native enzyme. Substrate deuterium isotope effects on the steady state V or V/K(m) parameters of 4-6 with ethanol or benzyl alcohol indicate that hydrogen transfer is a major rate-limiting step in catalysis. Steady state oxidation of benzyl alcohol is most rapid above a pK of about 9 for V and V/K(m) and is 2-fold faster in D(2)O than in H(2)O. The results are consistent with hydride transfer from a ground state zinc alkoxide that forms a low-barrier hydrogen bond with the hydroxyl group of Ser48.     

Incorporating the amino acid properties to predict the significance of missense mutations

Determining if missense mutations are deleterious is critical for the analysis of genes implicated in disease. However, the mutational effects of many missense mutations in databases like the Breast Cancer Information Core are unclassified. Several approaches have emerged recently to determine such mutational effects but none have utilized amino acid property indices. We modified a previously described phylogenetic approach by first classifying benign substitutions based on the assumption that missense mutations that are maintained in orthologs are unlikely to affect function. A consensus conservation score based on 16 amino acid properties was used to characterize the remaining substitutions. This approach was evaluated with experimentally verified T4 lysozyme missnese mutations and is shown to be able to sieve out putative biochemical and structurally important residues. The use of amino acid properties can enhance the prediction of biochemical and structurally important residues and thus also predict the significance of missense mutations.     

Stability analysis for the cavity-filling mutations of the Myb DNA-binding domain utilizing free-energy calculations

We have analyzed the effect of cavity-filling mutations on protein stability by means of free-energy calculations based on molecular dynamics simulations to identify the factors contributing to stability changes caused by the mutations. We have studied the DNA-binding domain of Myb, which has a cavity in one of three homologous repeat units, and analyzed a series of mutations with nonnatural and natural amino acids at a single site, which change the size of the cavity. We found that the calculated free-energy changes caused by the mutations are in excellent agreement with experimental data (correlation coefficient 0.98). The free-energy changes in the native and denatured states were independently compared with the unfolding free-energy change (deltadeltaG) and cavity-volume changes (deltaV), and it was found that deltadeltaG and deltaV correlate with the native-state free-energy changes but not with the denatured-state free-energy changes. Further analyses in terms of enthalpy and entropy show that compensation between entropy and enthalpy occurs in the denatured state but not in the native state. The main contribution to the native-state free energy was found to be van der Waals interactions associated with the cavity. We estimate that the decrease in free energy per methylene group, which results from filling the cavity, is about 2 to 3 kcal/mol. These results suggest that the stabilization of a protein by cavity-filling mutations be determined primarily by the free energy associated with the cavity volume in the native state.     

Probing the roles of conserved arginine-44 of Escherichia coli dihydrofolate reductase in its function and stability by systematic sequence perturbation analysis

Sequence perturbation analysis is a powerful method to reveal roles of an amino acid residue in function and stability of a protein. By using and improving this method, we studied roles of highly conserved Arg44 of Escherichia coli dihydrofolate reductase (DHFR) in its function and stability. Here, we introduced systematic amino acid substitutions at this position and found that all 19 kinds of amino acid substitutions were tolerated, but the mutations significantly reduced the enzymatic activity and the binding affinity toward the cofactor NADPH. Moreover, the mutational effects on the cofactor binding affinity were well correlated with those on the catalytic activity, indicating that the R44X mutations affect the catalytic activity mainly by modulating the cofactor binding affinity. On the other hand, thermal denaturation measurements showed that most mutations stabilized the protein. Comparison between the mutational effects and various amino acid indices taken from the AAindex database indicated that hydrophobicity and polarity are key determinants of amino acids favorable at this position. These results suggest that through electrostatic interactions Arg44 plays a functional role in retaining the cofactor binding affinity at the cost of the protein stability.