Search

Although bovine beta-lactoglobulin assumes a monomeric native structure at pH 3 in the absence of salt, the addition of salts stabilizes the dimer. Thermodynamics of the monomer-dimer equilibrium dependent on the salt concentration were studied by sedimentation equilibrium. The addition of NaCl, KCl, or guanidine hydrochloride below 1 M stabilized the dimer in a similar manner. On the other hand, NaClO(4) was more effective than other salts by about 20-fold, suggesting that anion binding is responsible for the salt-induced dimer formation, as observed for acid-unfolded proteins. The addition of guanidine hydrochloride at 5 M dissociated the dimer into monomers because of the denaturation of protein structure. In the presence of either NaCl or NaClO(4), the dimerization constant decreased with an increase in temperature, indicating that the enthalpy change (DeltaH(D)) of dimer formation is negative. The heat effect of the dimer formation was directly measured with an isothermal titration calorimeter by titrating the monomeric beta-lactoglobulin at pH 3.0 with NaClO(4). The net heat effects after subtraction of the heat of salt dilution, corresponding to DeltaH(D), were negative, and were consistent with those obtained by the sedimentation equilibrium. From the dependence of dimerization constant on temperature measured by sedimentation equilibrium, we estimated the DeltaH(D) value at 20 degrees C and the heat capacity change (DeltaC(p)) of dimer formation. In both NaCl and NaClO(4), the obtained DeltaC(p) value was negative, indicating the dominant role of burial of the hydrophobic surfaces upon dimer formation. The observed DeltaC(p) values were consistent with the calculated value from the X-ray dimeric structure using a method of accessible surface area. These results indicated that monomer-dimer equilibrium of beta-lactoglobulin at pH 3 is determined by a subtle balance of hydrophobic and electrostatic effects, which are modulated by the addition of salts or by changes in temperature.

Conformational analyses of a recombinant mouse tooth enamel amelogenin (rM179) were performed using circular dichroism (CD), fluorescence, differential scanning calorimetry, and sedimentation equilibrium studies. The results show that the far-UV CD spectra of rM179 at acidic pH and 10 degrees C are different from the spectra of random coil in 6 M GdnHCl. A near-UV CD spectrum of rM179 at 10 degrees C is similar to that of rM179 in 6 M GdnHCl, which indicates that aromatic residues of native structure are exposed to solvent and rotate freely. Far-UV CD values of rM179 at 80 degrees C are different from that of random-coil structure in 6 M GdnHCl, which suggests that rM179 at 80 degrees C has specific secondary structures. A gradual thermal transition was observed by far-UV CD, which is interpreted as a weak cooperative transition from specific secondary structures to other specific secondary structures. The fluorescence emission maximum for the spectrum due to Trp residues in rM179 at 10 degrees C shows the same fluorescence emission maximum as rM179 in 6 M GdnHCl and amino acid Trp, which indicates that the three Trp in rM179 are exposed to solvent. Deconvolution of differential scanning calorimetry curve gives the population of three states (A, I, and C states). These results indicate that three states (A, I, and C) have specific secondary structures, in which hydrophobic and Trp residues are exposed to the solvent. The thermodynamic characteristics of rM179 are unique and different from a typical globular protein, proline-rich peptides, and a molten globule state.

Recognition of specific DNA sequences by proteins is essential for regulation of gene expression. To fully understand the recognition mechanism, it is necessary to understand not only the structure of the specific protein-DNA interactions but also the energetics. We therefore performed a computer analysis in which a phage DNA-binding protein, lambda repressor, was used to examine the changes in binding free energy (DeltaDeltaG) and its energy components caused by single base mutations. We then determined which of the calculated energy components best correlated with the experimental data. The experimental DeltaDeltaG values were well reproduced by the calculations. Component analysis revealed that the electrostatic and hydrogen bond energies were most strongly correlated with the experimental data. Among the 51 single base-substitution mutants examined, positive DeltaDeltaG values, corresponding to weakened binding, were caused by the loss of favorable electrostatic interactions and hydrogen bonds, the introduction of steric collisions and electrostatic repulsion, the loss of favorable interactions with a thymine methyl group, and the increase of unfavorable hydration energy from isolated DNA. This analysis also showed distinct patterns of recognition at A-T and G-C positions, as different combinations of energy components were involved in DeltaDeltaG caused by the two substitution types. We have thus been able to identify the energy components that most strongly correlate with sequence-dependent DeltaDeltaG and determine their contribution to the specificity of DNA sequence recognition by the lambda repressor. Application of this method to other systems should provide additional insight into the molecular mechanism of protein-DNA recognition.

Acid-induced denaturation of the ribonuclease HI protein from Escherichia coli was analyzed by CD and NMR spectroscopies. The CD measurement revealed that the acid denaturation at 10°C proceeds from the native state (N-state) to a molten globule-like state (A-state), through an apparently more unfolded state (UA-state). In 1H–15N heteronuclear single-quantum coherence (HSQC) spectra, cross peaks from the N-state and those from the other two states are distinctively observed, while the UA-state and A-state are not distinguished from each other. Cross peaks from the UA/A-states showed a small pH dependence, which suggests a similarity in the backbone structure between the two states. The direct hydrogen–deuterium (H–D) exchange measurement at pH with the largest population of UA-state revealed that at least α-helix I is highly protected in the structure of the UA-state. A pH-jump H–D exchange analysis showed that the protection of α-helix I is highest also in the A-state. The profile of hydrogen-bond protection indicated that the structure of the A-state is closely related to that of the kinetic folding intermediate. © 2003 Wiley Periodicals, Inc. Biopolymers 69: 176–188, 2003

The activities of the eight mutant proteins of Escherichia coli RNase HI, in which the four carboxylic amino acids (Asp(10), Glu(48), Asp(70), and Asp(134)) involved in catalysis are changed to Asn (Gln) or Ala, were examined in the presence of Mn(2+). Of these proteins, the E48A, E48Q, D134A, and D134N proteins exhibited the activity, indicating that Glu(48) and Asp(134) are dispensable for Mn(2+)-dependent activity. The maximal activities of the E48A and D134A proteins were comparable to that of the wild-type protein. However, unlike the wild-type protein, these mutant proteins exhibited the maximal activities in the presence of100 microM MnCl(2), and their activities were not inhibited at higher Mn(2+) concentrations (up to 10 mM). The wild-type protein contains two Mn(2+) binding sites and is activated upon binding of one Mn(2+) ion at site 1 at low ( approximately 1 microM) Mn(2+) concentrations. This activity is attenuated upon binding of a second Mn(2+) ion at site 2 at high (10 microM) Mn(2+) concentrations. The cleavage specificities of the mutant proteins, which were examined using oligomeric substrates at high Mn(2+) concentrations, were identical to that of the wild-type protein at low Mn(2+) concentrations but were different from that of the wild-type protein at high Mn(2+) concentrations. These results suggest that one Mn(2+) ion binds to the E48A, E48Q, D134A, and D134N proteins at site 1 or a nearby site with weaker affinities. The binding analyses of the Mn(2+) ion to these proteins in the absence of the substrate support this hypothesis. When Mn(2+) ion is used as a metal cofactor, the Mn(2+) ion itself, instead of Glu(48) and Asp(134), probably holds water molecules required for activity.

The high resolution dielectric spectra of semidilute solutions of apomyoglobin in native (N, pH = 5), acid-induced molten globule (A, pH = 4), and unfolded (U(A), pH = 3) states have been measured in the range from 0.2 to 20 GHz. Based on a two-component mixture theory, we obtained the following hydration numbers per protein molecule: 590 +/- 65 for N, 630 +/- 73 for A, and 1110 +/- 67 for U(A). There was no clear difference between N and A states in contrast to the 25% reduction of helix content and the 50% reduction of heat capacity change upon unfolding. This suggests that the association of hydrophobic moieties might follow the disruption of secondary structures from N to A states. The measured hydration number of U(A) was close to that of the accessible water number (1340) of a protein molecule calculated for a fully extended structure, indicating that the structure of U(A) is extended but somewhat more compact than that of a fully extended state.

Thermodynamic data regarding proteins and their interactions are important for understanding the mechanisms of protein folding, protein stability, and molecular recognition. Although there are several structural databases available for proteins and their complexes with other molecules, databases for experimental thermodynamic data on protein stability and interactions are rather scarce. Thus, we have developed two electronically accessible thermodynamic databases. ProTherm, Thermodynamic Database for Proteins and Mutants, contains numerical data of several thermodynamic parameters of protein stability, experimental methods and conditions, along with structural, functional, and literature information. ProNIT, Thermodynamic Database for Protein-Nucleic Acid Interactions, contains thermodynamic data for protein-nucleic acid binding, experimental conditions, structural information of proteins, nucleic acids and the complex, and literature information. These data have been incorporated into 3DinSight, an integrated database for structure, function, and properties of biomolecules. A WWW interface allows users to search for data based on various conditions, with different display and sorting options, and to visualize molecular structures and their interactions. These thermodynamic databases, together with structural databases, help researchers gain insight into the relationship among structure, function, and thermodynamics of proteins and their interactions, and will become useful resources for studying proteins in the postgenomic era.

For understanding the factors influencing protein stability, we have analyzed the relationship between changes in protein stability caused by partially buried mutations and changes in 48 physico-chemical, energetic and conformational properties of amino acid residues. Multiple regression equations were derived to predict the stability of protein mutants and the efficiency of the method has been verified with both back-check and jack-knife tests. We observed a good agreement between experimental and computed stabilities. Further, we have analyzed the effect of sequence window length from 1 to 12 residues on each side of the mutated residue to include the sequence information for predicting protein stability and we found that the preferred window length for obtaining the highest correlation is different for each secondary structure; the preferred window length for helical, strand and coil mutations are, respectively, 0, 9 and 4 residues on both sides of the mutant residues. However, all the secondary structures have significant correlation for a window length of one residue on each side of the mutant position, implying the role of short-range interactions. Extraction of surrounding residue information for various distances (3 to 20A) around the mutant position showed the highest correlation at 8A, 6A and 7A, respectively, for mutations in helical, strand and coil segments. Overall, the information about the surrounding residues within the sphere of 7 to 8A, may explain better the stability in all subsets of partially buried mutations implying that this distance is sufficient to accommodate the residues influenced by major intramolecular interactions for the stability of protein structures.

Pressure denaturation of Escherichia coli ribonuclease HI (RNase HI) was studied by Fourier transform infrared (FTIR) and two-dimensional NMR spectroscopy at pD* 3.0 and 25 degrees C. A reversible transition in the pressure range of 0.1-1090 MPa was observed with second-derivative FTIR experiments. A cooperative and gradual denaturation, involving both the secondary and tertiary structures, was observed between 240 and 450 MPa. The two peaks at 1629 and 1652 cm(-1), due to beta-strands and alpha-helices, respectively, did not fully disappear after the denaturation, and are different from the spectra of the random coil peptides. The hydrogen-deuterium exchange rates of the individual backbone amide protons were determined by heteronuclear NMR combined with the pressure-jump technique at 500, 650, and 850 MPa. Although most of the amides protected in the native structure are also highly protected in the pressure-denatured state, the rate constants (0.048 +/- 0.007 min(-1)) for the amide protons at 500 MPa are similar regardless of their locations, which is an indication of the EX1 mechanism of hydrogen-deuterium exchange. The pressure-denatured state of RNase HI at 500 MPa represents a novel denatured state, which is different from a typical molten globule state at atmospheric pressure (0.1 MPa), from the viewpoint of the homogeneous rate constants. The observations at 650 MPa are essentially the same as those at 500 MPa. However, at 850 MPa, the amide exchange rates for the highly hydrophobic C-terminal half of alpha-helix I are significantly slower than those for the other part of the protein, which can be interpreted as a hydrophobic collapse centered at the C-terminal half of alpha-helix I.

Escherichia coli ribonuclease HI, which requires divalent cations (Mg2+ or Mn2+) for activity, was thermostabilized by 2.6-3.0 kcal/mol in the presence of the Mg2+, Mn2+, or Ca2+ ion, probably because the negative charge repulsion around the active site was canceled upon the binding of these metal ions. The dissociation constants were determined to be 0.71 mM for Mg2+, 0.035 mM for Mn2+, and 0.16 mM for Ca2+. Likewise, various active site mutants at Asp10, Glu48, Asp70, or Asp134 were thermostabilized by 0.4-3.0 kcal/mol in the presence of the Mg2+ ion, suggesting that this ion binds to these mutant proteins as well. The dissociation constants of Mg2+ were determined to be 9.8 mM for D10N, 1.1 mM for E48Q, 18.8 mM for D70N, and 1.8 mM for D134N. Thus, the mutation of Asp10 or Asp70 to Asn considerably impairs the Mg2+ binding, whereas the mutation of Glu48 to Gln or Asp134 to Asn does not. Comparison of the thermal stability of the mutant proteins with that of the wild-type protein in the absence of the Mg2+ ion suggests that the negative charge repulsion between Asp10 and Asp70 is responsible for the binding of the metal cofactor. Glu48 may be required to anchor a water molecule, which functions as a general acid.

Circular dichroism has been used to investigate the thermodynamics of the thermal unfolding of Escherichia coli ribonuclease HI as a function of pH over the pH range 0–4. This protein undergoes a reversible thermal conformational change from the native state to the denatured state with an isodichroic point. The calculated thermodynamic values at pH 3.0 are as follows: tm = 50.2°C, ΔHm (tm) = 93.8 kcal mol−1, ΔG (25°) = 6.08 kcal mol−1, and ΔG (10°C) = 8.14 kcal mol−1. At pH 4, ΔG (25°C) = 10 kcal mol−1 and ΔG (10°C) = 12 kcal mol−1. At a pH below 2, this protein denatures at 25°C with ΔG (25°C) = −1 kcal mol−1, but it is stable at 10°C with ΔG (10°C) = 2 kcal mol−1. The ΔCp value determined from the ΔHm(Tm) versus Tm plot is 1.4 kcal mol−1 K−1. The thermal unfolding curves at pH values above 2.18 showed a highly cooperative thermal transition. Cold denaturation was observed at temperatures below 10°C between pH 2.03 and 1.55. A cooperative heat denaturation was also observed over this pH range. At pH values lower than 1.45, cold denaturation was not observed. A broad thermal transition was observed between pH 0.83 and 0.53.

In order to explain the large permittivity of liquid cyanogen consisting of centrosymmetric molecules, a molecular dynamics simulation which includes all degrees of freedom of vibration, rotation, and translation has been performed. The intramolecular potential energy is expressed in terms of the internal coordinates with quadratic, cubic, and quartic force constants. The intermolecular potential consists of an atom–atom pair potential with Coulombic and Buckingham terms. The simulation has given an adequate value of self‐diffusion coefficient and a typical form of radial distribution function for the liquid state. The simulation has reproduced well the large value of the permittivity observed in experiments, suggesting that the value results from the extraordinarily large atomic polarization brought by the intramolecular bending motion in an antisymmetric mode. The simulation has also suggested that the nitrogen atoms collide with each other more frequently than with carbon atoms and the collision shifts the normal frequencies related to the N≡C bond to higher frequency. The infrared and Raman spectra together with the real and imaginary parts of the refractive index have also been calculated.

A strategy to genetically select Escherichia coli ribonuclease HI mutants with enhanced thermostability is described. E. coli strain MIC3001, which shows an RNase H-dependent, temperature-sensitive growth phenotype, was used for this purpose. Introduction of the rnhA gene permits the growth of this temperature-sensitive strain, whereas the gene for the truncated protein, 142-RNase HI, which lacks the carboxyl-terminal 13 residues, cannot. Analyses of the production levels and the stability of a series of mutant proteins with COOH-terminal truncations suggested that 142-RNase HI is nonfunctional in vivo because of a dramatic decrease in the protein stability. Polymerase chain reaction-mediated random mutagenesis of the rnhA142 gene, encoding 142-RNase HI, followed by selection of revertants, allowed us to isolate 11 single amino acid substitutions that render 142-RNase HI functional in vivo. Of them, eight substitutions were shown to enhance the thermal stability of the wild-type RNase HI protein, and of these, six were novel. The genetic selection strategy employed in this experiment was thus shown to be effective for identifying amino acid substitutions that enhance the thermal stability of E. coli RNase HI. Such a strategy would be versatile if a protein of interest could be destabilized by a deletion or a truncation and a conditional-lethal strain were available.

Protein de novo design has been performed, as an exercise of the inverse folding problem. A beta/alpha-barrel protein was designed and synthesized using the Escherichia coli expression system for the structural characterization. A tertiary model with a two-fold symmetry was built, based upon the geometrical parameters extracted from X-ray crystal structures of several beta/alpha-barrel proteins. Amino acid frequencies at each position on the alpha- and beta-structures were investigated, and an amino acid sequence with 201 residues was designed. The associated gene was chemically synthesized and the fusion protein with human growth hormone was expressed in Escherichia coli. The purified protein after being cleaved and refolded was found to be stable and globular with the large amount of secondary structures. However, it has similar characteristics to the molten globules of natural proteins, with loose packing of side-chains. The approach for the tight packing is discussed.

To investigate the structural stability of proteins, we analyzed the thermodynamics of an artificially designed 30-residue peptide. The designed peptide, NH2-EELLPLAEALAPLLEALLPLAEALAPLLKK-COOH (PERI COIL-1), with prolines at i + 7 positions, forms a pentameric alpha-helical structure in aqueous solution. The thermal denaturation curves of the CD at 222 nm (pH 7.5) show an unusual cold denaturation occurring well above 0 degrees C and no thermal denaturation is observable under 90 degrees C. This conformational change is reversible and depends on peptide concentration. A 2-state model between the monomeric denatured state (5D) and the pentameric helical state (H5) was sufficient to analyze 5 thermal denaturation curves of PERI COIL-1 with concentrations between 23 and 286 microM. The analysis was carried out by a nonlinear least-squares method using 3 fitting parameters: the midpoint temperature, Tm, the enthalpy change, delta H(Tm), and the heat capacity change, delta Cp. The association number (n = 5) was determined by sedimentation equilibrium and was not used as a fitting parameter. The heat capacity change suggests that the hydrophobic residues are buried in the helical state and exposed in the denatured one, as it occurs normally for natural globular proteins. On the other hand, the enthalpy and the entropy changes have values close to those found for coiled-coils and are quite distinct from typical values reported for natural globular proteins. In particular, the enthalpy change extrapolated at 110 degrees C is about 3 kJ/mol per amino acid residue, i.e., half of the value found for globular proteins.(ABSTRACT TRUNCATED AT 250 WORDS)

Molecular dynamics simulations were performed for the dissociation and association (D/A) reactions N2O4■2 NO2 in the gas phase and in liquid N2O4. The trajectory was initialized from an equilibrium distribution of all variables in liquid N2O4, except the reactive mode, the NN distance of a reactant NO2 pair, was excited above the dissociation limit of the Morse‐like potential between NO2 fragments, and the dynamics were calculated for 500 fs both forward and backward in time. Characteristics of the translational and vibrational energy relaxations of the reactant were studied in detail. Energy ERT, which is defined to be the sum of the potential and kinetic energies of interfragment motion, is found to play a key role in the D/A dynamics; a reactant pair is associated when ERT

In order to study the cation dependence of the ionic dynamics in molten nitrates, molecular dynamics simulations including vibrational degrees of freedom were carried out for molten LiNO3, NaNO3, and RbNO3. Coulomb pair potential with Born‐type repulsion was adopted for the interionic interaction. The simulated diffusion coefficient was smaller for a larger cation, and that of nitrate ions did not change with changing cation species. The mean squared charge displacements showed that the static conductivity decreased considerably as the cation size increased from Li+ to Rb+. The simulated orientational correlation function of nitrate ions decayed more quickly as the cation size increased. Far infrared absorption spectrum simulated from the time evolution of the dipole moment (or the current) of the system showed that the peak shifted to the low energy side and the intensity decreased as the cation size increased. Results of the simulation were compared with the experimental diffusion constants, static and ...

In summary, the thermal denaturation of proteins has been elucidated in terms of the chain free energy and the hydration free energy as follows. (1) Method to calculate the unfolding free energy. The free energy of unfolding consists of two contributions: the hydration around the molecule, and the intramolecular interactions. A method to calculate the free energy of hydration from the accessible surface area (ASA) of the constituent atomic groups in a protein has been developed. This assumes a proportionality between the free energy and the ASA, where the proportional constants were determined by least-squares fitting to the experimentally derived thermodynamic data on small molecules. Similarly, the free energy of unfolding for the chain in vacuo can be also calculated from the ASA, using the unfolding thermodynamics derived from the experimental data of the ten proteins. (2) Thermodynamics of protein unfolding predicted from the three-dimensional structures and from the amino acid content in proteins. First, our method is applied to predict the thermodynamics of protein unfolding from the X-ray structure. The predicted values of four test proteins agree well with the experimentally derived values. It also accounts for the temperature dependence of the free energy and of the enthalpy upon unfolding for 14 proteins. Second, this method is applied to the helix-coil transition of short peptides of poly(L-Ala)20 and Ac-(AAAAK)3A-NH2. The calculated enthalpy change is close to the experimental values for poly-L-Lys and poly-L-Glu. Since delta Hcu at 25 degrees C significantly contributes to delta Gu, the helix formation is enthalpy-driven through interactions in the chain. Third, the method is applied to predict the unfolding thermodynamics of a globular protein from its amino acid content. It also accounts for the temperature dependence of the free energy of unfolding for the 14 proteins. The agreement between the experimental and the calculated values by this method for the 14 proteins is not so different from those obtained with the three-dimensional structures. Fourth, the values of delta Cpu for 14 proteins may be closely approximated to the predicted values of delta Cp,hu. The delta Cp,hu value in a protein consists of the major contribution from the hydrophobic and the aromatic residues, and the minor one from the hydrophilic residues. (3) Dominant free energies in protein folding.(ABSTRACT TRUNCATED AT 400 WORDS)

To study the microdynamics of translational, orientational and vibrational motions in liquid N2O4, molecular dynamics simulations, including all degrees of freedom of vibrating N2O4 molecules, were carried out. Atom-atom Lennard-Jones potentials with and without Coulomb potentials (LJ + C and LJ, respectively), were adopted in the simulation. When Coulomb interaction between atoms was included, mean-squared force and torque increased, the velocity and angular momentum correlation functions decayed faster to more evident negative parts, and both the translational and orientational diffusion slowed down. Simulated rotational correlation functions for LJ + C agreed reasonably with those derived from Raman spectroscopy, while those for LJ were found to decay too fast. Vibrational correlation functions for the NN stretching v 3 mode were calculated and the vibrational spectra were obtained from their Fourier transforms. Two simulations for the harmonic and the anharmonic intramolecular potentials, where anharm...

A thirty-residue peptide (PERI COIL-1) has been designed with a new type of α-helical structure, which is capable of folding into an amphiphilic helix bending at 4 periodic prolines in the sequence. Two such helices should form a dimer by supercoiling about one another in an antiparallel direction in the design. With this arrangement, close packing between them is maintained through the hydrophobic interaction pattern called ‘leucine zipper’. PERI COIL-1 has been obtained by solid-phase peptide synthesis, and characterized by circular dichroic spectroscopy, sedimentation equilibrium experiments and NMR. The result of the analyses shows that it preferentially forms a helical tetramer in aqueous solution.

Molecular dynamics simulation has been used to investigate the infrared spectra of ordered and disordered benzoic acid crystals consisting of various proportions of the A and B configurations, whose C–O and C=O bonds are approximately parallel to the a axis, respectively. The simulation reproduced well the splitting and the temperature dependence in intensity for the C=O stretching bands observed. The absorption coefficient of the lower frequency band was simulated as several times larger than that of the higher frequency band in accord with the observed. To interpret this peculiar result of the molecular dynamics simulation, a simple calculation of the band splitting has been made in terms of the transition dipole–transition dipole interaction. It has been shown that the interaction is an essential factor to produce the large difference in intensity and the splitting for the C=O stretching bands. The higher and lower bands are interpreted as due to the A configuration reduced in intensity by the B form, and due to the B form enhanced in intensity by the A form, respectively. Compression effect in the O⋅⋅⋅O distance of the hydrogen bond is also studied in comparison with the observed value.

In order to study the effect of the potential well on the single ion dynamics and the charge transport in molten LiNO3, a molecular dynamics simulation, which included all degrees of freedom of vibrating nitrate ions, was carried out. Simple Coulomb pair potential with Born‐type repulsion was adopted as a standard, and a potential well of varying depth between Li+ and O of NO−3 was added to the standard potential. As the well depth increased, both the translational and orientational diffusion were found to slow down, leading to solid like dynamics. The potential dependence of the ionic dynamics was studied from the translational, orientational, and current correlation functions. Results of the simulation were compared with the experimental diffusion constants, static conductivity, and rotational behavior revealed by Raman spectroscopy. Using the relation between the infrared absorption spectrum and the autocorrelation function of the dipole moment (or the current) of the system, we carried out the first s...

Understanding the mechanisms by which mutations affect protein stability is one of the most important problems in molecular biology. In this work, we analyzed the relationship between changes in protein stability caused by surface mutations and changes in 49 physicochemical, energetic, and conformational properties of amino acid residues. We found that the hydration entropy was the major contributor to the stability of surface mutations in helical segments; other properties responsible for size and volume of molecule also correlated significantly with stability. Classification of coil mutations based on their locations in the (ϕ–ψ) map improved the correlation significantly, demonstrating the existence of a relationship between stability and strain energy, which indicates that the role of strain energy is very important for the stability of surface mutations. We observed that the inclusion of sequence and structural information raised the correlation, indicating the influence of surrounding residues on the stability of surface mutations. Further, we examined the previously reported “inverse relationship” between stability and hydrophobicity, and observed that the inverse hydrophobic effect was generally applicable only to coil mutations. The present study leads to a simple method for predicting protein stability changes caused by amino acid substitutions, which will be useful for protein engineering in designing novel proteins with increased stability and altered function. © 2002 Wiley Periodicals, Inc. Biopolymers 64: 210–220, 2002

ProTherm 2.0 is the second release of the Thermo-dynamic Database for Proteins and Mutants that includes numerical data for several thermodynamic parameters, structural information, experimental methods and conditions, functional and literature information. The present release contains >5500 entries, an approximately 67% increase over the previous version. In addition, we have included information about reversibility of data, details about buffer and ion concentrations and the surrounding residues in space for all mutants. A WWW interface enables users to search data based on various conditions with different sorting options for outputs. Further, ProTherm has links with other structural and literature databases, and the mutation sites and surrounding residues are automatically mapped on the structures and can be directly viewed through 3DinSight developed in our laboratory. The ProTherm database is freely available through the WWW at http://www.rtc.riken.go.jp/protherm.html

Understanding the role of various interactions in enhancing the thermostability of proteins is important not only for clarifying the mechanism of protein stability but also for designing stable proteins. In this work, we have analyzed the thermostability of 16 different families by comparing mesophilic and thermophilic proteins with 48 various physicochemical, energetic and conformational properties. We found that the increase in shape, s (location of branch point in side chain) increases the thermostability, whereas, an opposite trend is observed for Gibbs free energy change of hydration for native proteins, G hN , in 14 families. A good correlation is observed between these two properties and the simultaneous increases of − G hN and s is necessary to enhance the thermostability from mesophile to thermophile. The increase in shape, which tends to increase with increasing number of carbon atoms both for polar and non-polar residues, may generate more packing and compactness, and the position of β and higher order branches may be important for better packing. On the other hand, the increase in − G hN in thermophilic proteins increases the solubility of the proteins. This tendency counterbalances the increases in insolubility and unfolding heat capacity change due to the increase in the number of carbon atoms. Thus, the present results suggest that the stability of thermophilic proteins may be achieved by a balance between better packing and solubility.

The first release of the Thermodynamic Database for Proteins and Mutants (ProTherm) contains more than 3300 data of several thermodynamic parameters for wild type and mutant proteins. Each entry includes numerical data for unfolding Gibbs free energy change, enthalpy change, heat capacity change, transition temperature, activity etc., which are important for understanding the mechanism of protein stability. ProTherm also includes structural information such as secondary structure and solvent accessibility of wild type residues, and experimental methods and other conditions. A WWW interface enables users to search data based on various conditions with different sorting options for outputs. Further, ProTherm is cross-linked with NCBI PUBMED literature database, Protein Mutant Database, Enzyme Code and Protein Data Bank structural database. Moreover, all the mutation sites associated with each PDB structure are automatically mapped and can be directly viewed through 3DinSight developed in our laboratory. The database is available at the URL, http://www.rtc.riken.go.jp/protherm.htm l

Escherichia coli ribonuclease HI has a cavity within the hydrophobic core. Two core residues, Ala52 and Val74, resided at both ends of this cavity. We have constructed a series of single mutant proteins at Ala52, and double mutant proteins, in which Ala52 was replaced by Gly, Val, Ile, Leu, or Phe, and Val74 was replaced by Ala or Leu. All of these mutant proteins, except for A52W, A52R, and A52G/V74A, were overproduced and purified. Measurement of the thermal denaturations of the proteins at pH 3.2 by CD suggests that the cavity is large enough to accommodate three methyl or methylene groups without creating serious strains. A correlation was observed between the protein stability and the hydrophobicity of the substituted residue. As a result, a number of the mutant proteins were more stable than the wild-type protein. The stabilities of the mutant proteins with charged or extremely bulky residues at the cavity were lower than those expected from the hydrophobicities of the substituted residues, suggesting that considerable strains are created at the mutation sites in these mutant proteins. However, examination of the far- and near-UV CD spectra and the enzymatic activities suggest that all of the mutant proteins have structures similar to that of the wild-type protein. These results suggest that the cavity in the hydrophobic core of E. coli RNase HI is conformationally fairly stable. This may be the reason why the cavity-filling mutations effectively increase the thermal stability of this protein.

The unfolding and refolding processes of Escherichia coli ribonuclease HI at 25 degrees C, induced by concentration jumps of either guanidine hydrochloride (GuHCl) or urea, were investigated using stopped-flow circular dichroism (CD), stopped-flow fluorescence, and NMR spectroscopies. Only a single exponential process was detected for the fast time scale unfolding (rate constants from 0.014 to 0.54 s-1, depending on the final denaturant concentration). For refolding, the far-UV CD value largely recovered within 50 ms of the stopped-flow mixing dead time (burst phase). This phase was followed by either one or two phases, with rate constants from 0.035 to 2.45 s-1 as detected by CD and fluorescence, respectively. Although this protein has a single cis-Pro residue, a very slow phase due to proline isomerization was not observed, for either unfolding or refolding. The difference in the amplitudes of the burst phases for refolding in the far- and near-UV CD spectra revealed that an intermediate state exists, with the characteristics of a molten globule. Because the one-phased fast exponential process detected by CD corresponds to the slower of the two phases detected by fluorescence, the intermediate detected by CD might be the most stable. GuHCl denaturation experiments revealed that this intermediate cooperatively unfolds, with a transition midpoint of 1.33 +/- 0.03 M. The Gibbs free energy difference (delta G) between the intermediate and the unfolded states, under physiological conditions (25 degrees C, pH 5.5, and 0 M GuHCl), was estimated to be 20.0 +/- 2.3 kJ mol-1. Therefore, it is reasonable to assume that the refolding intermediate, rather than the unfolded state, is the latent denatured state under physiological conditions. Approximately linear relationships between the GuHCl concentration and the logarithm of the microscopic rate constants determined by CD and fluorescence were also observed. By extrapolation to a GuHCl concentration of 0 M, activation Gibbs free energies of 98.5 +/- 1.1 kJ mol-1 for unfolding and 69.5 +/- 0.2 kJ mol-1 for refolding under physiological conditions were obtained. The hydrogen-exchange-refolding competition combined with two-dimensional NMR revealed that the amide protons of alpha-helix I are the most highly protected, suggesting that alpha-helix I is the initial site of protein folding. The CD and NMR data showed that the intermediate state has a structure similar to that of the acid-denatured molten globule.

In order to understand the contribution of hydrophobic residues to the conformational stability of human lysozyme, five Ile mutants (Ile --Val) in the interior of the protein were constructed. The thermodynamic parameters characterizing the denaturation of these mutant proteins were determined by scanning calorimetry, and the three-dimensional structure of each mutant protein was solved at high resolution by X-ray crystallography. The thermodynamic analyses at 64.9 degrees C and at pH 2.7 revealed the following. (1) The stabilities of all the mutant proteins were decreased as compared with that of the wild-type protein. (2) The changes in the calorimetric enthalpies were larger than those in the Gibbs energies, and were compensated by entropy changes. (3) The destabilization mechanism of the mutant proteins differs, depending on the location of the mutation sites. X-ray analyses showed that the overall structures of all the mutant human lysozymes examined were identical to that of the wild-type protein, and only small structural rearrangements were observed locally around some of the mutation sites. The most striking change among the mutant proteins was found in the mutant protein, 159V, which contains a new water molecule in the cavity created by the mutation. The thermodynamic stabilities of the mutant proteins are discussed in light of the high-resolution X-ray structures of the wild-type and five mutant human lysozymes examined.

To test whether the combination of multiple thermostabilizing mutations is a useful strategy to generate a hyperstable mutant protein, five mutations, Gly23-->Ala, His62-->Pro, Val74-->Leu, Lys95-->Gly, and Asp134-->His or Asn, were simultaneously introduced into Escherichia coli ribonuclease HI. The enzymatic activities of the resultant quintuple mutant proteins, 5H- and 5N-RNases HI, which have His and Asn at position 134, respectively, were 35 and 55% of that of the wild-type protein. The far-UV and near-UV CD spectra of these mutant proteins were similar to those of the wild-type protein, suggesting that the mutations did not seriously affect the tertiary structure of the protein. The differences in the free energy change of unfolding between the wild-type and mutant proteins, delta delta G, were estimated by analyzing the thermal denaturation of the proteins by CD. The 5H-RNase HI protein, which was slightly more stable than the 5N-RNase HI, was more stable than the wild-type protein by 20.2 degrees C in Tm and 5.6 kcal/mol in delta G at pH 5.5. In addition, the 5H-RNase HI was highly resistant to proteolysis and acid denaturation. The effects of each mutation on the thermal stability and the susceptibility to chymotryptic digestion were nearly cumulative, and the 5H-RNase HI undergoes chymotryptic digestion at a rate that is 41 times slower than that of the wild-type protein. Good correlation was observed between the thermal stability and the resistance to chymotryptic digestion for all proteins examined.(ABSTRACT TRUNCATED AT 250 WORDS)

The role of the conserved Asp134 residue in Escherichia coli ribonuclease HI, which is located at the center of the alpha V helix and lies close to the active site, was analyzed by means of site-directed random mutagenesis. Mutant rnhA genes encoding proteins with ribonuclease H activities were screened by their ability to suppress the ribonuclease-H-dependent, temperature-sensitive growth phenotype of E. coli strain MIC3001. Based on the DNA sequences, nine mutant proteins were predicted to have ribonuclease H activity in vivo. All of these mutant proteins were purified to homogeneity and examined for enzymic activity and protein stability. Among them, only the mutant proteins [D134H]RNase H and [D134N]RNase H were shown to have considerable ribonuclease H activities. Determination of the kinetic parameters revealed that replacement of Asp134 by amino acid residues other than asparagine and histidine dramatically decreased the enzymic activity without seriously affecting the substrate binding. Determination of the CD spectra indicated that none of the mutations seriously affected secondary and tertiary structure. The protein stability was determined from the thermal denaturation curves. All mutant proteins were more stable than the wild-type protein. Such stabilization effects would be a result of a reduction in the negative charge repulsion between Asp134 and the active-site residues, and/or an enhancement of the stability of the alpha V helix. These results strongly suggest that Asp134 does not contribute to the maintenance of the molecular architecture but the carboxyl oxygen at its delta 1 position impacts catalysis.

Thermal stabilities of mutant ribonuclease HI proteins from Escherichia coli , in which each of five histidine residues was replaced with alanine, were examined at various pHs. Increases in the T m values were observed at pH 3.0 for four of the mutant proteins, in which each of the four histidine residues exposed to the solvent was mutated, as compared to the T m of the wild-type protein. The thermostabilization of three of the mutant proteins was dependent on pH, and only observed at low pH. The thermostabilizing effects of the His → Ala substitutions were cumulative. The temperature of the midpoint of the transition in the thermal unfolding curves, T m , of the most stable mutant enzyme, in which His 62, His 83, His 124, and His 127 were replaced by Ala, was 5.5°C higher than that of the wild-type enzyme at pH 3.0. The stability of the wild-type protein decreased as the pH was lowered below pH 4, a condition favoring the protonation of carboxyl groups, probably due to unfavorable electrostatic interactions introduced by the increase in positive charges on the protein. Since imidazole groups are positively charged at pH 3.0, it seems likely that thermal stabilization at pH 3.0 by a His → Ala substitution would be the result of a reduction in such unfavorable electrostatic interactions. These results suggest that amino acid substitutions that cause a decrease in the number of positive charges on the surface of a protein can be used as a general strategy to enhance protein stability at pH values below pH 4.

An α/β-barrel protein was designed de novo and synthesized. Based upon the characteristic structural parameters extracted from the native proteins, a model structure composed of the typical parallel β-barrel with eight β-strands and eight α-helices was built, accompanying with the unique short turns between α and β structures. A unique amino acid sequence of the total 201 amino acid was designed. The associated gene was chemically synthesized and expressed in Eschericia coli as a fusion protein with human growth hormone. The protein obtained after cleaved and refolded was found to be stable and globular with a large amount of secondary structures.

In summary, the thermal denaturation of proteins has been elucidated in terms of the chain free energy and the hydration free energy as follows. (1) Method to calculate the unfolding free energy. The free energy of unfolding consists of two contributions: the hydration around the molecule, and the intramolecular interactions. A method to calculate the free energy of hydration from the accessible surface area (ASA) of the constituent atomic groups in a protein has been developed. This assumes a proportionality between the free energy and the ASA, where the proportional constants were determined by least-squares fitting to the experimentally derived thermodynamic data on small molecules. Similarly, the free energy of unfolding for the chain in vacuo can be also calculated from the ASA, using the unfolding thermodynamics derived from the experimental data of the ten proteins. (2) Thermodynamics of protein unfolding predicted from the three-dimensional structures and from the amino acid content in proteins. First, our method is applied to predict the thermodynamics of protein unfolding from the X-ray structure. The predicted values of four test proteins agree well with the experimentally derived values. It also accounts for the temperature dependence of the free energy and of the enthalpy upon unfolding for 14 proteins. Second, this method is applied to the helix-coil transition of short peptides of poly(L-Ala)20 and Ac-(AAAAK)3A-NH2. The calculated enthalpy change is close to the experimental values for poly-L-Lys and poly-L-Glu. Since delta Hcu at 25 degrees C significantly contributes to delta Gu, the helix formation is enthalpy-driven through interactions in the chain. Third, the method is applied to predict the unfolding thermodynamics of a globular protein from its amino acid content. It also accounts for the temperature dependence of the free energy of unfolding for the 14 proteins. The agreement between the experimental and the calculated values by this method for the 14 proteins is not so different from those obtained with the three-dimensional structures. Fourth, the values of delta Cpu for 14 proteins may be closely approximated to the predicted values of delta Cp,hu. The delta Cp,hu value in a protein consists of the major contribution from the hydrophobic and the aromatic residues, and the minor one from the hydrophilic residues. (3) Dominant free energies in protein folding.(ABSTRACT TRUNCATED AT 400 WORDS)

The method given earlier for predicting the thermodynamics of protein unfolding from the x-ray structure of a protein is applied here to the poly(L-alanine) helix. First, the fitting parameters derived earlier from a data base of 10 proteins were used to predict the unfolding thermodynamics of 4 other proteins. The agreement between the observed and predicted values is comparable to that found for the 10 proteins studied initially. Next, the temperature dependences of the Gibbs energy and enthalpy changes for unfolding of bacteriophage T4 lysozyme were predicted and compared with data in the literature. The predicted and observed temperature dependences are similar and the predicted results indicate that cold denaturation should be observed at low temperatures, as observed recently for a T4 lysozyme mutant. The fitting parameters derived from thermodynamic data for protein unfolding and for hydration of model compounds were used to predict the unfolding thermodynamics of the poly(L-alanine) helix. The results predict that helix formation is enthalpy-driven, and the predicted enthalpy change for unfolding (0.86 kcal per mol per residue) is close to the value found in a recent calorimetric study of a 50-residue alanine-rich helix.

We have developed the database, "ProTherm:Thermodynamic Database for Proteins and Mutants," which contains more than 11,000 numerical data for several important thermodynamic properties, structural information of proteins and mutants, experimental methods and conditions, and, functional and literature information. ProTherm can be searched through WWW on various conditions with different options for output. In addition, it is cross-linked with structure, function and literature databases such as PDB, PMD, ENZYME, SWISS-PROT, PIR and PubMed. The structure of mutation sites and surrounding residues can be directly viewed through 3DinSight, which is developed in our laboratory. ProTherm is freely available at above URL.

Non-bonded energy of 16 proteins was calculated using the atomic co-ordinates obtained by X-ray crystallography. The curve of total energy against the number of atoms in proteins is approximately linear with a slight concaved shape. According to a linear equation to fit the curve, the extrapolated length of a polypeptide chain of a globular shape is expected to be 18 residues, which corresponds conceivably to an approximate size of nucleus for a folding of the polypeptide chain. Contributions from short-range and medium-range energies are always much greater than those from long-range energy for all the proteins and there seems to exist a change of each contribution in a range from 1200 to 1700 atoms. The energies with a lag less than four residues are a major part of the total energy and the contribution of energy from main-chain atoms is considerably higher than that from side-chain atoms. Side-chain atoms of a residue have a tendency to interact more strongly with main-chain atoms of N-terminal, than with those of C-terminal side of the residue, indicating asymmetry of the interaction in a protein. Amino acid residues in proteins may be divided into three groups by the order of strength of average energy. The first group exhibiting strong interaction consists mainly of hydrophobic amino acids and the third group consists of hydrophilic ones corresponding to the location in a protein molecule. Cys, val, leu and met are important for medium-range and long-range energies; gly and ala for medium-range energy; ile, trp, phe, tyr and arg for long-range energy. One simple application of the average energy of amino acid residues is illustrated to estimate local energy of a segment of nine residues given by a protein sequence. There is a good correlation between the curve computed by the average energy and the experimental curve for myoglobin.

To reveal the molecular orientation of plastocyanin (PC) on spinach thylakoid membranes, the position of Lys residues modified by acetic anhydride was compared between thylakoid-bound PC and the isolated one. Digestion of the isolated PC by a trypsin yielded a peptide map with seven spots prior to the acetylation of the protein; none of the spots appeared after the isolated PC was acetylated. On the other hand, there were two spots on the peptide map of the PC acetylated when it was bound to the thylakoids. Those spots were revealed by their amino acid compositions to correspond to the peptide fragments Phe 82-Lys 95 and Val 96-Asn 99. Thus, the Lys residues 81 and 95 of the thylakoid-bound PC were not acetylated. These results suggest that the PC molecule binds to the thylakoids with a specific region including the Lys's 81 and 95 in contact with the membranes. The Lys's 81 and 95 are located near Tyr 83, which has been thought to be the delivery site of electrons from the Cu2+ center.

The conformational stability of a protein in aqueous solution is described in terms of the thermodynamic properties such as unfolding Gibbs free energy, which is the difference in the free energy (Gibbs function) between the native and random conformations in solution. The properties are composed of two contributions, one from enthalpy due to intramolecular interactions among constituent atoms and chain entropy of the backbone and side chains, and the other from the hydrated water around a protein molecule. The hydration free energy and enthalpy at a given temperature for a protein of known three-dimensional structure can be calculated from the accessible surface areas of constituent atoms according to a method developed recently. Since the hydration free energy and enthalpy for random conformations are computed from those for an extended conformation, the thermodynamic properties of unfolding are evaluated quantitatively. The evaluated hydration properties for proteins of known transition temperature (Tm) and unfolding enthalpy (delta Hm) show an approximately linear dependence on the number of constituent heavy atoms. Since the unfolding free energy is zero at Tm, the enthalpy originating from interatomic interactions of a polypeptide chain and the chain entropy are evaluated from an experimental value of delta Hm and computed properties due to the hydrated water around the molecule at Tm. The chain enthalpy and entropy thus estimated are largely compensated by the hydration enthalpy and entropy, respectively, making the unfolding free energy and enthalpy relatively small. The computed temperature dependences of the unfolding free energy and enthalpy for RNase A, T4 lysozyme, and myoglobin showed a good agreement with the experimental ones.(ABSTRACT TRUNCATED AT 250 WORDS)

The native conformation of a protein may be expressed in terms of the dihedral angles, phi's and psi's for the backbone, and kappa's for the side chains, for a given geometry (bond lengths and bond angles). We have developed a method to obtain the dihedral angles for a low-energy structure of a protein, starting with the X-ray structure; it is applied here to examine the degree of flexibility of bovine pancreatic trypsin inhibitor. Minimization of the total energy of the inhibitor (including nonbonded, electrostatic, torsional, hydrogen bonding, and disulfide loop energies) yields a conformation having a total energy of -221 kcal/mol and a root mean square deviation between all atoms of the computed and experimental structures of 0.63 A. The optimal conformation is not unique, however, there being at least two other conformations of low-energy (-222 and -220 kcal/mol), which resemble the experimental one (root mean square deviations of 0.66 and 0.64 A, respectively). These three conformations are located in different positions in phi, psi space, i.e., with a total deviation of 81 degrees, 100 degrees and 55 degrees from each other (with a root mean square deviation of several degrees per dihedral angle from each other). The nonbonded energies of the backbones, calculated along lines in phi, psi space connecting these three conformations, are all negative, without any intervening energy barriers (on an energy contour map in the phi, psi plane). Side chains were attached at several representative positions in this plane, and the total energy was minimized by varying the kappa's. The energies were of approximately the same magnitude as the previous ones, indicating that the conformation of low energy is flexible to some extent in a restricted region of phi, psi space. Interestingly, the difference delta phi i+1 in phi i+1 for the (i + 1)th residue from one conformation to another is approximately the same as -delta psi i for the ith residue; i.e., the plane of the peptide group between the ith and (i + 1)th residues re-orient without significant changes in the positions of the other atoms. The flexibility of the orientations of the planes of the peptide groups is probably coupled in a cooperative manner to the flexibility of the positions of the backbone and side-chain atoms.

A method is described for the inclusion of the effects of hydration in empirical conformational energy computations on polypeptides. The free energy of hydration is composed of additive contributions of various functional groups. The hydration of each group is assumed to be proportional to the accessible surface area of the group. The constants of proportionality, representing the free energy of hydration per unit area of accessible surface, have been evaluated for seven classes of groups (occurring in peptides) by least-squares fitting to experimental free energies of solution of small monofunctional aliphatic and aromatic molecules. The same method has also been applied to the modeling of the enthalpy and heat capacity of hydration, each of which is computed from the accessible surface area.