Free energy simulation to investigate the effect of amino acid sequence environment on the severity of osteogenesis imperfecta by glycine mutations in collagen

Molecular dynamics simulations were carried out to calculate the free energy change difference of two collagen-like peptide models for Gly --> Ser mutations causing two different osteogenesis imperfecta phenotypes. These simulations were performed to investigate the impact of local amino acid sequence environment adjacent to a mutation site on the stability of the collagen. The average free energy differences for a Gly --> Ser mutant relative to a wild type are 3.4 kcal/mol and 8.2 kcal/mol for a nonlethal site and a lethal site, respectively. The free energy change differences of mutant containing two Ser residues relative to the wild type at the nonlethal and lethal mutation sites are 4.6 and 9.8 kcal/mol, respectively. Although electrostatic interactions stabilize mutants containing one or two Ser residues at both mutation sites, van der Waals interactions are of sufficient magnitude to cause a net destabilization. The presence of Gln and Arg near the mutation site, which contain large and polar side chains, provide more destabilization than amino acids containing small and nonpolar side chains.     

The severity of osteogenesis imperfecta: a comparison to the relative free energy differences of collagen model peptides

Molecular dynamics simulations were carried out to calculate free energy differences between the folded and unfolded states of wild type and mutant collagen model peptides. The calculated stability of the collagen models was compared with the severity of osteogenesis imperfecta. Free energy differences of Gly → Xaa (Xaa: Ser, Cys, Glu, and Asp) mutations between the wild type and the mutants at position 15 of the model peptide were 3.8, 4.2, 5.6, and 8.8 kcal/mol, respectively. The corresponding free energy differences of a second Gly mutation at the same position in different chains were, on average, 1.3, 1.5, 2.9, and 5.4 kcal/mol, respectively. Free energy simulations were also performed to estimate the relative stability between an oxidized form and a reduced form of the mutants containing two Cys residues, which indicated that the mutant of the collagen-like peptide containing an intramolecular disulfide bond was more stable than the mutant containing one Cys residue but less stable than the wild type. The calculated free energy differences between an oxidized and a reduced form of the mutants containing two Cys residues are 0.8 and 2.6 kcal/mol for the disulfide bonds between Chains A and B and between Chains A and C, respectively.     

Molecular Dynamics/Free Energy Perturbation Studies of the Thermostable V74I Mutant of Ribonuclease HI

Abstract

Recent site-directed mutagenesis and thermodynamic studies have shown that the V74I mutant of Escherichia coli ribonuclease HI (RNase HI) is more stable than the wild type protein [Ishikawa et al., Biochemistry 32, 6171 (1993)]. In order to clarify the stabilization mechanism of this mutant, we calculated the free energy change due to the mutation Val 74→Ile in both the native and denatured states by free energy perturbations based on molecular dynamics (MD) simulations. We carried out inclusive MD simulations for the protein in water; i.e., fully solvated, no artificial constraints applied, and all long-range Coulomb interactions included. We found that the free energy of the mutant increased slightly relative to the wild type, in the native state by 1.60 kcal/mol, and in the denatured state by 2.25 kcal/mol. The unfolding free energy increment of the mutant (0.66 ± 0.19 kcal/mol) was in good agreement with the experimental value (0.6 kcal/mol). The hysteresis error in the free energy calculations, i.e., forward and reverse perturbations, was only ±0.19 kcal/mol. These results show that the V74I mutant is stabilized relative to the wild type by the increased free energy of the denatured state and not by a decrease in the free energy of the native state as had been proposed earlier based on the mutant X-ray structure. It was found that the stabilization was caused by a loss of solvation energy in the mutant denatured state and not by improved packing interactions inside the native protein.     

Thermal denaturation of iso-1-cytochrome c variants: comparison with solvent denaturation.

Thermal denaturation studies as a function of pH were carried out on wild-type iso-1-cytochrome c and three variants of this protein at the solvent-exposed position 73 of the sequence. By examining the enthalpy and Tm at various pH values, the heat capacity increment (delta Cp), which is dominated by the degree of change in nonpolar hydration upon protein unfolding, was found for the wild type where lysine 73 is normally present and for three variants. For the Trp 73 variant, the delta Cp value (1.15 +/- 0.17 kcal/mol K) decreased slightly relative to wild-type iso-1-cytochrome c (1.40 +/- 0.06 kcal/mol K), while for the Ile 73 (1.65 +/- 0.07 kcal/mol K) and the Val 73 (1.50 +/- 0.06 kcal/mol K) variants, delta Cp increased slightly. In previous studies, the Trp 73, Ile 73, and Val 73 variants have been shown to have decreased m-values in guanidine hydrochloride denaturations relative to the wild-type protein (Hermann L, Bowler BE, Dong A, Caughey WS. 1995. The effects of hydrophilic to hydrophobic surface mutations on the denatured state of iso-1-cytochrome c: Investigation of aliphatic residues. Biochemistry 34:3040-3047). Both the m-value and delta Cp are related to the change in solvent exposure upon unfolding and other investigators have shown a correlation exists between these two parameters. However, for this subset of variants of iso-1-cytochrome c, a lack of correlation exists which implies that there may be basic differences between the guanidine hydrochloride and thermal denaturations of this protein. Spectroscopic data are consistent with different denatured states for thermal and guanidine hydrochloride unfolding. The different response of m-values and delta Cp for these variants will be discussed in this context.     

Free energy perturbation study on a Trp-binding mutant (Ser88-->Cys) of the trp-repressor.

The Ser88-->Cys mutant of the trp-repressor showed a lower affinity for the corepressor than the wild-type repressor [delta delta G = 1.7 +/- 0.3 kcal/mol, Chou and Matthews (1989) J. Biol. Chem., 264, 18314-18319]. A molecular dynamics/free energy cycle perturbation study was performed to understand the origin of the decreased affinity. A value (delta delta G = 1.58 +/- 0.28 kcal/mol) comparable with the experimental value was obtained by the simulation. Free energy component analysis revealed that destabilization of the van der Waals interaction between Ser88 and Trp109 (corepressor) mainly contributed to the decreased affinity of the mutant. The rotational transition of the hydroxyl (sulfhydryl) group of Ser88 (Cys88) during the simulations affected the contributions of Arg84 and water to the free energy change in the aporepressor and those of Arg84 and Trp109 to that in the holorepressor. However, the contributions from different residues compensated each other, and the total free energy changes were almost invariable in the various simulations.     

Binding mechanism of spinosine and venenatine molecules with p300 HAT enzyme: Molecular screening,molecular dynamics and free-energy analysis

The chromatin modification is regulated by the histone acetyltransferase (HAT) and histone deacetyltransferase (HDAC) enzymes; abnormal function of these enzymes leads to several malignant diseases. The inhibition of these enzymes using natural ligand molecules is an emerging technique to cure these diseases. The in vitro analysis of natural molecules, venenatine, spinosine, palmatine and taxodione are giving the best inhibition rate against p300 HAT enzyme. However, the detailed understanding of binding and the stability of these molecules with p300 HAT is not yet known. The aim of the present study is focused to determine the binding strength of the molecules from molecular dynamics simulation analysis. The docking analysis confirms that, the venenatine (−6.97 kcal/mol - conformer 8), spinosine (−6.52 kcal/mol conformer −10), palmatine (−5.72 kcal/mol conformer-3) and taxodione (−4.99 kcal/mol conformer-4) molecules form strong hydrogen bonding interactions with the key amino acid residues (Arg1410, Thr1411 and Trp1466) present in the active site of p300. In the molecular dynamics (MD) simulation, the spinosine retain these key interactions with the active site amino acid residues (Arg1410, Thr1411, and Trp1466) than venenatine and are stable throughout the simulation. The RMSD value of spinosine (0.5 to 1.3 Å) and venenatine (0.3 to 1.3 Å) are almost equal during the MD simulation. However, during the MD simulation, the intermolecular interaction between venenatine and the active site amino acid residues (Arg1410, Thr1411, and Trp1466) decreased on comparing with the spinosine-p300 interaction. The binding free energy of the spinosine (−15.30 kcal/mol) is relatively higher than the venenatine (−11.8 kcal/mol); this increment is attributed to the strong hydrogen bonding interactions of spinosine molecule with the active site amino acid residues of p300.     

Functional role of putative critical residues in Mycobacterium tuberculosis RNase P protein

RNase P is involved in processing the 5⿲ end of pre-tRNA molecules. Bacterial RNase P contains a catalytic RNA subunit and a protein subunit. In this study, we have analyzed the residues in RNase P protein of M. tuberculosis that differ from the residues generally conserved in other bacterial RNase Ps. The residues investigated in the current study include the unique residues, Val27, Ala70, Arg72, Ala77, and Asp124, and also Phe23 and Arg93 which have been found to be important in the function of RNase P protein components of other bacteria. The selected residues were individually mutated either to those present in other bacterial RNase P protein components at respective positions or in some cases to alanine. The wild type and mutant M. tuberculosis RNase P proteins were expressed in E. coli, purified, used to reconstitute holoenzymes with wild type RNA component in vitro, and functionally characterized. The Phe23Ala and Arg93Ala mutants showed very poor catalytic activity when reconstituted with the RNA component. The catalytic activity of holoenzyme with Val27Phe, Ala70Lys, Arg72Leu and Arg72Ala was also significantly reduced, whereas with Ala77Phe and Asp124Ser the activity of holoenzyme was similar to that with the wild type protein. Although the mutants did not suffer from any binding defects, Val27Phe, Ala70Lys, Arg72Ala and Asp124Ser were less tolerant towards higher temperatures as compared to the wild type protein. The K_m of Val27Phe, Ala70Lys, Arg72Ala and Ala77Phe were >2-fold higher than that of the wild type, indicating the substituted residues to be involved in substrate interaction. The study demonstrates that residues Phe23, Val27 and Ala70 are involved in substrate interaction, while Arg72 and Arg93 interact with other residues within the protein to provide it a functional conformation.     

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.     

Mutation of the hydrophobic residue on helix alpha5 of the Bacillus thuringiensis Cry4B affects structural stability

Cry4B toxin is a mosquito-larvicidal protein from the Bacillus thuringiensis subsp. israelensis. We have investigated the role of two conserved hydrophobic residues of Cry4B in structural stabilization. Substitutions of the leucine-175 and isoleucine-189 on helix alpha5 with valine and leucine did not affect the expression level, solubility and proteolytic processing. Steady state analysis of an unfolding experiment as monitored by circular dichroism and fluorescence spectroscopy demonstrated a typical two-state transition. The determined unfolding free energy for the L175V mutant revealed a structural destabilization of 10.49 kcal/mol relative to the wild type. However unfolding kinetic analysis gave identical activation energy for wild type and both mutants. Our findings suggested that a perturbation on the close packing of the hydrophobic side chains in protein interior could lead to a significant destabilization of the native conformation.     

Stabilizing nonpolar/polar side-chain interactions in the alpha-helix.

A simplistic, yet often used, view of protein stability is that amino acids attract other amino acids with similar polarity, whereas nonpolar and polar side chains repel. Here we show that nonpolar/polar interactions, namely Val or Ile bonding to Lys or Arg in alpha-helices, can in fact be stabilizing. Residues spaced i, i + 4 in alpha-helices are on the same face of the helix, with potential to favorably interact and stabilize the structure. We observe that the nonpolar/polar pairs Ile-Lys, Ile-Arg, and Val-Lys occur in protein helices more often than expected when spaced i, i + 4. Partially helical peptides containing pairs of nonpolar/polar residues were synthesized. Controls with i, i + 5 spacing have the residues on opposite faces of the helix and are less helical than the test peptides with the i, i + 4 interactions. Experimental circular dichroism results were analyzed with helix-coil theory to calculate the free energy for the interactions. All three stabilize the helix with DeltaG between -0.14 and -0.32 kcal x mol(-1). The interactions are hydrophobic with contacts between Val or Ile and the alkyl groups in Arg or Lys. Side chains such as Lys and Arg can thus interact favorably with both polar and nonpolar residues.     

The role of tryptophan residues in an integral membrane protein: diacylglycerol kinase

Tryptophan residues are thought to play special roles in integral membrane proteins, anchoring transmembrane alpha-helices into the lipid bilayer. We have studied the effect of mutating the five Trp residues in the diacylglycerol kinase (DGK) of Escherichia coli to Leu residues. The fluorescence emission maxima for DGK and a variety of Trp mutants in bilayers of dioleoylphosphatidylcholine [di(C18:1)PC] are all centered at ca. 327 nm, suggesting that all five Trp residues are located close to the glycerol backbone region of the bilayer. This is also consistent with fluorescence quenching experiments, measuring the separation between the Trp residues and the bromine atoms in a bilayer of dibromostearoylphosphatidylcholine. Mutation of Trp residues in DGK was found to have significant effects on activity for DGK reconstituted into bilayers of di(C18:1)PC containing 30 mol % 1,2-dihexanoylglycerol (DHG). Of the mutants containing a single Trp residue, only that containing Trp-112 was found to give active protein. The presence of both Trp-25 and Trp-112 gave higher activity than Trp-112 alone. Trp-25 and Trp-112 are the most important Trp residues in DGK as far as activity is concerned. Effects of mutations on K(m) for DHG were generally greater than effects on v(max). The activity of wild-type and mutant DHGs reconstituted into bilayers of phosphatidylcholines was sensitive to the chain length of the phospholipid, with highest activities for chain lengths of C18 or C20 and lower activities in phosphatidylcholines with shorter or longer chains. Compared to wild-type DGK, the Trp mutants were less affected by long-chain phosphatidylcholines but more affected by short-chain phospholipids. In mutants lacking Trp-25, low activities in short-chain phospholipids followed from a decrease in v(max) compared to wild type, combined with an increase in K(m) value for DHG, as observed in the wild type. It is suggested that Trp-25 plays a role in maintaining the alignment of ATP and DHG at the active site. Fluorescence emission spectra for the Trp mutants do not change significantly with changing fatty acyl chain length from C14 to C24, showing efficient hydrophobic matching between DGK and the surrounding lipid bilayer. It is suggested that hydrophobic matching is achieved by tilting of the transmembrane alpha-helix or rotation of residues at the ends of the helices about the Calpha-Cbeta bond linking the residue to the helix backbone. As well as any structural effects, the presence of Trp residues in DGK has a clear effect on thermal stability.     

Structural insight into the kinetics and DeltaCp of interactions between TEM-1 beta-lactamase and beta-lactamase inhibitory protein (BLIP)

In a previous study, we examined thermodynamic parameters for 20 alanine mutants in beta-lactamase inhibitory protein (BLIP) for binding to TEM-1 beta-lactamase. Here we have determined the structures of two thermodynamically distinctive complexes of BLIP mutants with TEM-1 beta-lactamase. The complex BLIP Y51A-TEM-1 is a tight binding complex with the most negative binding heat capacity change (DeltaG = approximately -13 kcal mol(-1) and DeltaCp = approximately -0.8 kcal mol(-1) K(-1)) among all of the mutants, whereas BLIP W150A-TEM-1 is a weak complex with one of the least negative binding heat capacity changes (DeltaG = approximately -8.5 kcal mol(-1) and DeltaCp = approximately -0.27 kcal mol(-1) K(-1)). We previously determined that BLIP Tyr51 is a canonical and Trp150 an anti-canonical TEM-1-contact residue, where canonical refers to the alanine substitution resulting in a matched change in the hydrophobicity of binding free energy. Structure determination indicates a rearrangement of the interactions between Asp49 of the W150A BLIP mutant and the catalytic pocket of TEM-1. The Asp49 of W150A moves more than 4 angstroms to form two new hydrogen bonds while losing four original hydrogen bonds. This explains the anti-canonical nature of the Trp150 to alanine substitution, and also reveals a strong long distance coupling between Trp150 and Asp49 of BLIP, because these two residues are more than 25 angstroms apart. Kinetic measurements indicate that the mutations influence the dissociation rate but not the association rate. Further analysis of the structures indicates that an increased number of interface-trapped water molecules correlate with poor interface packing in a mutant. It appears that the increase of interface-trapped water molecules is inversely correlated with negative binding heat capacity changes.     

Analysis of the interactions of human ribonuclease inhibitor with angiogenin and ribonuclease A by mutagenesis: importance of inhibitor residues inside versus outside the C-terminal "hot spot"

Ribonuclease inhibitor (RI) binds diverse mammalian RNases with extraordinary avidity. Here, we have investigated the structural basis for this tight binding and broad specificity by mutational analysis of the complexes of RI with angiogenin (Ang) and RNase A (K(D)=0.5 fM and 43 fM, respectively). Both crystal structures are known; the interfaces are large, and the ligands dock similarly, although few of the specific interactions formed are analogous. Our previous mutagenesis studies focused primarily on one contact region, containing RI 434-438 and the enzymatic active site. Many single-residue replacements produced extensive losses of binding energy (2.3-5.9 kcal/mol), suggesting that this region constitutes a "hot spot" in both cases. We have now explored the roles of most of the remaining RI residues that interact with Ang and/or RNase A. One major cluster in each complex lies in a Trp-rich area of RI, containing Trp261, Trp263, Trp318, and Trp375. Although the energy losses from individual replacements in this portion of the Ang complex were small-to-moderate (0-1.5 kcal/mol), the changes from multiple substitutions were much greater than additive, and the binding energy provided by this region is estimated to be approximately 6 kcal/mol (30 % of total). Effects of replacing combinations of hot spot components had also been found to be superadditive, and this negative cooperativity is now shown to extend to the neighboring contact residue RI Ser460. The overall contribution of the hot spot, taking superadditivity into account, is then approximately 14-15 kcal/mol. The hot spot and Trp-rich regions, although spatially well separated, are themselves functionally linked. No other parts of the RI-Ang interface appear to be energetically important. Binding of RNase A is more sensitive to substitutions throughout the interface, with free energy losses>/=1 kcal/mol produced by nearly all replacements examined, so that the sum of losses greatly exceeds the binding energy of the complex. This discrepancy can be explained, in part, by positive cooperativity, as evident from the subadditive effects observed when combinations of residues in either the hot spot or Trp-rich region are replaced. These findings suggest that the binding energy may be more widely distributed in the RNase A complex than in the Ang complex.     

Importance of stalk segment S5 for intramolecular communication in the sarcoplasmic reticulum Ca2+-ATPase

Sixteen residues in stalk segment S5 of the Ca(2+)-ATPase of sarcoplasmic reticulum were studied by site-directed mutagenesis. The rate of the Ca(2+) binding transition, determined at 0 degrees C, was enhanced relative to wild type in mutants Ile(743) --> Ala, Val(747) --> Ala, Glu(748) --> Ala, Glu(749) --> Ala, Met(757) --> Gly, and Gln(759) --> Ala and reduced in mutants Asp(737) --> Ala, Asp(738) --> Ala, Ala(752) --> Leu, and Tyr(754) --> Ala. In mutant Arg(762) --> Ile, the rate of the Ca(2+) binding transition was wild type like at 0 degrees C, whereas it was 3.5-fold reduced relative to wild type at 25 degrees C. The rate of dephosphorylation of the ADP-insensitive phosphoenzyme was increased conspicuously in mutants Ile(743) --> Ala and Tyr(754) --> Ala (close to 20-fold in the absence of K(+)) and increased to a lesser extent in Asn(739) --> Ala, Glu(749) --> Ala, Gly(750) --> Ala, Ala(752) --> Gly, Met(757) --> Gly, and Arg(762) --> Ile, whereas it was reduced in mutants Asp(737) --> Ala, Val(744) --> Gly, Val(744) --> Ala, Val(747) --> Ala, and Ala(752) --> Leu. In mutants Ile(743) --> Ala, Tyr(754) --> Ala, and Arg(762) --> Ile, the apparent affinities for vanadate were enhanced 23-, 30-, and 18-fold, respectively, relative to wild type. The rate of Ca(2+) dissociation was 11-fold increased in Gly(750) --> Ala and 2-fold reduced in Val(747) --> Ala. Mutants with alterations to Arg(751) either were not expressed at a significant level or were completely nonfunctional. The findings show that S5 plays a crucial role in mediating communication between the Ca(2+) binding pocket and the catalytic domain and that Arg(751) is important for both structural and functional integrity of the enzyme.     

Crevice-forming mutants of bovine pancreatic trypsin inhibitor: stability changes and new hydrophobic surface.

Four mutants of bovine pancreatic trypsin inhibitor (BPTI) with replacements in the rigid core result in the creation of deep crevices on the surface of the protein. Other than crevices at the site of the mutation, few other differences are observed in the crystal structures of wild-type BPTI and the mutants F22A, Y23A, N43G, and F45A. These mutants are highly destabilized relative to wild type (WT). The differences between WT and mutants in the free energy change associated with cooperative folding/unfolding, delta delta G0 (WT-->mut), have been measured by calorimetry, and they are in good agreement with delta delta G0(WT-->mut) values from hydrogen exchange rates. For F22A the change in free energy difference is about 1.7 kcal/mol at 25 degrees C; for the other three mutants it is in the range of 5-7 kcal/mol at 25 degrees C. The experimental delta delta G0(WT-->mut) values of F22A, Y23A, and F45A are reasonably well accounted for as the sum of two terms: the difference in transfer free energy change, and a contribution from exposure to solvent of new surface (Eriksson, A.E., et al., 1992, Science 255, 178-183), if the recently corrected transfer free energies and surface hydrophobicities (De Young, L. & Dill, K., 1990, J. Phys. Chem. 94, 801-809; Sharp, K.A., et al., 1991a, Science 252, 106-109) are used and only nonpolar surface is taken into account. In N43G, three protein-protein hydrogen bonds are replaced by protein-water hydrogen bonds.(ABSTRACT TRUNCATED AT 250 WORDS)     

Prediction of protein mutant stability using classification and regression tool

Prediction of protein stability upon amino acid substitutions is an important problem in molecular biology and the solving of which would help for designing stable mutants. In this work, we have analyzed the stability of protein mutants using two different datasets of 1396 and 2204 mutants obtained from ProTherm database, respectively for free energy change due to thermal (DeltaDeltaG) and denaturant denaturations (DeltaDeltaG(H(2)O)). We have used a set of 48 physical, chemical energetic and conformational properties of amino acid residues and computed the difference of amino acid properties for each mutant in both sets of data. These differences in amino acid properties have been related to protein stability (DeltaDeltaG and DeltaDeltaG(H(2)O)) and are used to train with classification and regression tool for predicting the stability of protein mutants. Further, we have tested the method with 4 fold, 5 fold and 10 fold cross validation procedures. We found that the physical properties, shape and flexibility are important determinants of protein stability. The classification of mutants based on secondary structure (helix, strand, turn and coil) and solvent accessibility (buried, partially buried, partially exposed and exposed) distinguished the stabilizing/destabilizing mutants at an average accuracy of 81% and 80%, respectively for DeltaDeltaG and DeltaDeltaG(H(2)O). The correlation between the experimental and predicted stability change is 0.61 for DeltaDeltaG and 0.44 for DeltaDeltaG(H(2)O). Further, the free energy change due to the replacement of amino acid residue has been predicted within an average error of 1.08 kcal/mol and 1.37 kcal/mol for thermal and chemical denaturation, respectively. The relative importance of secondary structure and solvent accessibility, and the influence of the dataset on prediction of protein mutant stability have been discussed.     

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.     

Effects of five-tryptophan mutations on structure, stability and function of Escherichia coli dihydrofolate reductase

To elucidate the roles of tryptophan residues in the structure, stability, and function of Escherichia coli dihydrofolate reductase (DHFR), its five tryptophan residues were replaced by site-directed mutagenesis with leucine, phenylalanine or valine (W22F, W22L, W30L, W47L, W74F, W74L, W133F, and W133V). Far-ultraviolet circular dichroism (CD) spectra of these mutants reveal that exciton coupling between Trp47 and Trp74 strongly affects the peptide CD of wild-type DHFR, and that Trp133 also contributes appreciably. No additivity was observed in the contributions of individual tryptophan residues to the fluorescence spectrum of wild-type DHFR, Trp74 having a dominant effect. These single-tryptophan mutations induce large changes in the free energy of urea unfolding, which showed values of 1.79-7.14 kcal/mol, compared with the value for wild-type DHFR of 6.08 kcal/mol. Analysis of CD and fluorescence spectra suggests that thermal unfolding involves an intermediate with the native-like secondary structure, the disrupted Trp47-Trp74 exciton coupling, and the solvent-exposed Trp30 and Trp47 side chains. All the mutants except W22L (13%) retain more than 50% of the enzyme activity of wild-type DHFR. These results demonstrate that the five tryptophan residues of DHFR play important roles in its structure and stability but do not crucially affect its enzymatic function.     

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.     

Insights into protein-protein binding by binding free energy calculation and free energy decomposition for the Ras-Raf and Ras-RalGDS complexes

Absolute binding free energy calculations and free energy decompositions are presented for the protein-protein complexes H-Ras/C-Raf1 and H-Ras/RalGDS. Ras is a central switch in the regulation of cell proliferation and differentiation. In our study, we investigate the capability of the molecular mechanics (MM)-generalized Born surface area (GBSA) approach to estimate absolute binding free energies for the protein-protein complexes. Averaging gas-phase energies, solvation free energies, and entropic contributions over snapshots extracted from trajectories of the unbound proteins and the complexes, calculated binding free energies (Ras-Raf: -15.0(+/-6.3)kcal mol(-1); Ras-RalGDS: -19.5(+/-5.9)kcal mol(-1)) are in fair agreement with experimentally determined values (-9.6 kcal mol(-1); -8.4 kcal mol(-1)), if appropriate ionic strength is taken into account. Structural determinants of the binding affinity of Ras-Raf and Ras-RalGDS are identified by means of free energy decomposition. For the first time, computationally inexpensive generalized Born (GB) calculations are applied in this context to partition solvation free energies along with gas-phase energies between residues of both binding partners. For selected residues, in addition, entropic contributions are estimated by classical statistical mechanics. Comparison of the decomposition results with experimentally determined binding free energy differences for alanine mutants of interface residues yielded correlations with r(2)=0.55 and 0.46 for Ras-Raf and Ras-RalGDS, respectively. Extension of the decomposition reveals residues as far apart as 25A from the binding epitope that can contribute significantly to binding free energy. These "hotspots" are found to show large atomic fluctuations in the unbound proteins, indicating that they reside in structurally less stable regions. Furthermore, hotspot residues experience a significantly larger-than-average decrease in local fluctuations upon complex formation. Finally, by calculating a pair-wise decomposition of interactions, interaction pathways originating in the binding epitope of Raf are found that protrude through the protein structure towards the loop L1. This explains the finding of a conformational change in this region upon complex formation with Ras, and it may trigger a larger structural change in Raf, which is considered to be necessary for activation of the effector by Ras.