Thermodynamics of barnase unfolding.

The thermodynamics of barnase denaturation has been studied calorimetrically over a broad range of temperature and pH. It is shown that in acidic solutions the heat denaturation of barnase is well approximated by a 2-state transition. The heat denaturation of barnase proceeds with a significant increase of heat capacity, which determines the temperature dependencies of the enthalpy and entropy of its denaturation. The partial specific heat capacity of denatured barnase is very close to that expected for the completely unfolded protein. The specific denaturation enthalpy value extrapolated to 130 degrees C is also close to the value expected for the full unfolding. Therefore, the calorimetrically determined thermodynamic characteristics of barnase denaturation can be considered as characteristics of its complete unfolding and can be correlated with structural features--the number of hydrogen bonds, extent of van der Waals contacts, and the surface areas of polar and nonpolar groups. Using this information and thermodynamic information on transfer of protein groups into water, the contribution of various factors to the stabilization of the native structure of barnase has been estimated. The main contributors to the stabilization of the native state of barnase appear to be intramolecular hydrogen bonds. The contributions of van der Waals interactions between nonpolar groups and those of hydration effects of these groups are not as large if considered separately, but the combination of these 2 factors, known as hydrophobic interactions, is of the same order of magnitude as the contribution of hydrogen bonding.     

Characterization of the denaturation of human alpha-lactalbumin in urea by molecular dynamics simulations

Molecular dynamics (MD) simulations were used to characterize the non-cooperative denaturation of the molten globule A-state of human alpha-lactalbumin by urea. A solvent of explicit urea and water molecules was used, corresponding to a urea concentration of approximately 6M. Three simulations were performed at temperatures of 293K, 360K and 400K, with lengths of 2 ns, 8 ns and 8 ns respectively. The results of the simulations were compared with experimental data from NMR studies of human alpha-lactalbumin and related peptides. During the simulations, hydrogen bonds were formed from the protein to both urea and water molecules as intra-protein hydrogen bonds were lost. Urea was shown to compete efficiently with water as both a hydrogen bond donor and acceptor. Radial distribution functions of water and urea around hydrophobic side chain atoms showed a significant increase in urea molecules in the solvation shell as the side chains became exposed during denaturation. A considerable portion of the native-like secondary structure persisted throughout the simulations. However, in the simulations at 360K and 400K, there were substantial changes in the packing of aromatic and other hydrophobic side chains in the protein, and many native contacts were lost. The results suggest that during the non-cooperative denaturation of the molten globule, secondary structure elements are stabilized by non-specific, non-native interactions.     

Molecular dynamics simulations of urea and thermal-induced denaturation of S-peptide analogue

Molecular dynamics simulations of the S-peptide analogue AETAAAKFLREHMDS in water at 278 and 358 K, and in 8 M urea at 278 K were performed. The results show agreement with experiments. The helix is stable at low temperature (278 K), while at 358 K, unfolding is observed. The effects of urea on protein stability have been studied. The data support a model in which urea denatures proteins by: (1) diminishing the hydrophobic effect by displacing water molecules from the solvent shell around nonpolar groups; and (2) binding directly to amide units (NH and CO groups) via hydrogen bonds. The results of cluster analysis and essential dynamics analysis suggest that the mechanism of urea and thermal-induced denaturation may not be the same.     

Effect of urea on peptide conformation in water: molecular dynamics and experimental characterization

Molecular dynamics simulations of a ribonuclease A C-peptide analog and a sequence variant were performed in water at 277 and 300 K and in 8 M urea to clarify the molecular denaturation mechanism induced by urea and the early events in protein unfolding. Spectroscopic characterization of the peptides showed that the C-peptide analog had a high alpha-helical content, which was not the case for the variant. In the simulations, interdependent side-chain interactions were responsible for the high stability of the alpha-helical C-peptide analog in the different solvents. The other peptide displayed alpha-helical unwinding that propagated cooperatively toward the N-terminal. The conformations sampled by the peptides depended on their sequence and on the solvent. The ability of water molecules to form hydrogen bonds to the peptide as well as the hydrogen bond lifetimes increased in the presence of urea, whereas water mobility was reduced near the peptide. Urea accumulated in excess around the peptide, to which it formed long-lived hydrogen bonds. The unfolding mechanisms induced by thermal denaturation and by urea are of a different nature, with urea-aqueous solutions providing a better peptide solvation than pure water. Our results suggest that the effect of urea on the chemical denaturation process involves both the direct and indirect mechanisms.     

Polar or Apolar—The Role of Polarity for Urea-Induced Protein Denaturation

Urea-induced protein denaturation is widely used to study protein folding and stability; however, the molecular mechanism and driving forces of this process are not yet fully understood. In particular, it is unclear whether either hydrophobic or polar interactions between urea molecules and residues at the protein surface drive denaturation. To address this question, here, many molecular dynamics simulations totalling ca. 7 µs of the CI2 protein in aqueous solution served to perform a computational thought experiment, in which we varied the polarity of urea. For apolar driving forces, hypopolar urea should show increased denaturation power; for polar driving forces, hyperpolar urea should be the stronger denaturant. Indeed, protein unfolding was observed in all simulations with decreased urea polarity. Hyperpolar urea, in contrast, turned out to stabilize the native state. Moreover, the differential interaction preferences between urea and the 20 amino acids turned out to be enhanced for hypopolar urea and suppressed (or even inverted) for hyperpolar urea. These results strongly suggest that apolar urea–protein interactions, and not polar interactions, are the dominant driving force for denaturation. Further, the observed interactions provide a detailed picture of the underlying molecular driving forces. Our simulations finally allowed characterization of CI2 unfolding pathways. Unfolding proceeds sequentially with alternating loss of secondary or tertiary structure. After the transition state, unfolding pathways show large structural heterogeneity.     

分子动力学模拟尿素对水溶液中蛋白质构象转变的影响

采用分子动力学方法和全原子模型研究尿素和水分子对模型蛋白S-肽链结构转化的影响。模拟结果显示S-肽链的变性速率常数k值随着尿素浓度的增加而先降低后升高,在尿素浓度为2.9 mol/L时达到最低值。模拟了不同尿素浓度下尿素-肽链、水-肽链以及肽链分子氢键的形成状况。结果表明:尿素浓度较低时,尿素分子与S-肽链的极性氨基酸侧链形成氢键,但不破坏其分子内的骨架氢键,尿素在S-肽链水化层外形成限制性空间,增强了S-肽链的稳定性。随着尿素的升高,尿素分子进入S-肽链内部并与其内部氨基酸残基形成氢键,导致S-肽链的骨架氢键丧失,S-肽链发生去折叠。上述模拟结果与文献报道的实验结果一致,从分子水平上揭示了尿素对蛋白质分子结构变化的影响机制,对于研究和发展蛋白质折叠及稳定化技术具有指导意义。     

The structural properties of magainin in water, TFE/water, and aqueous urea solutions: molecular dynamics simulations

Here, the MD simulations and comparative structural analysis of Magainin in water, TFE/water, and 2M, 4M, and BM urea solutions is reported. For MAG-TFE/water and MAG-2M urea the largely alpha helical conformation of the peptide is maintained throughout the 9-ns simulation. While in water, 4M urea, and 8M urea, the helix length decreases and at the same time helix radius increases. This suggests a more destabilized magainin secondary structure. Our simulation data reveals that the stabilizing effect of TFE is induced by preferential accumulation of TFE molecules around the alpha helical peptide. These results indicate that an aqueous urea solution solvates the surface of polypeptide chain more favorably than pure water. Urea molecules interact more favorably with nonpolar groups of the peptide in comparison with water, and the presence of urea improves the interactions of water molecules with the hydrophilic groups of the peptide. At 8M urea, there are more direct interactions between the urea and solute, and the helix is destabilized. At 2M urea, the interaction of urea molecules and nonpolar residues are weak, therefore, the presence of urea molecules decreases the interactions of water molecules with hydrophilic groups. Urea could not deteriorate the peptide secondary structure with time from an initial helix structure.     

Apolar and polar solvation thermodynamics related to the protein unfolding process.

Thermodynamics related to hydrated water upon protein unfolding is studied over a broad temperature range (5-125 degrees C). The hydration effect arising from the apolar interior is modeled as an increased number of hydrogen bonds between water molecules compared with bulk water. The corresponding contribution from the polar interior is modeled as a two-step process. First, the polar interior breaks hydrogen bonds in bulk water upon unfolding. Second, due to strong bonds between the polar surface and the nearest water molecules, we assume quantization using a simplified two-state picture. The heat capacity change upon hydration is compared with model compound data evaluated previously for 20 different proteins. We obtain good correspondence with the data for both the apolar and the polar interior. We note that the effective coupling constants for both models have small variations among the proteins we have investigated.     

A Gaussian statistical mechanical model for the equilibrium thermodynamics of barnase folding

Given an all non-hydrogen-atom potential function that implicitly includes solvation effects, it is possible to adjust its parameters to favor the correct native structure for several proteins over decoys produced by ungapped threading. It is also possible to further train it to reproduce the experimental free energy of unfolding in aqueous solution at 298 K for wild-type barnase and 66 mutants. For this, the native state is represented by the crystal structure at a single energy level with a calculated low degeneracy; the denatured state is represented by the extended conformation and a high calculated degeneracy. The same two-state model can be extended to account for the stability of all 67 sequences toward urea denaturation at 298 K by building in a solvation term that depends on urea concentration. With the addition of one more parameter set to give the correct heat capacity of unfolded barnase in solution, it is possible to approximate the experimental thermodynamics of barnase thermal denaturation: melting temperature, width of thermal transition, deltaG, deltaH, deltaS, and deltaCp. This requires a novel sort of statistical mechanical model where the two states each have a Gaussian density of microscopic state distribution as a function of energy.     

Replica exchange simulation of reversible folding/unfolding of the Trp-cage miniprotein in explicit solvent: on the structure and possible role of internal water

We simulate the folding/unfolding equilibrium of the 20-residue miniprotein Trp-cage. We use replica exchange molecular dynamics simulations of the AMBER94 atomic detail model of the protein explicitly solvated by water, starting from a completely unfolded configuration. We employ a total of 40 replicas, covering the temperature range between 280 and 538 K. Individual simulation lengths of 100 ns sum up to a total simulation time of about 4 micros. Without any bias, we observe the folding of the protein into the native state with an unfolding-transition temperature of about 440 K. The native state is characterized by a distribution of root mean square distances (RMSD) from the NMR data that peaks at 1.8A, and is as low as 0.4A. We show that equilibration times of about 40 ns are required to yield convergence. A folded configuration in the entire extended ensemble is found to have a lifetime of about 31 ns. In a clamp-like motion, the Trp-cage opens up during thermal denaturation. In line with fluorescence quenching experiments, the Trp-residue sidechain gets hydrated when the protein opens up, roughly doubling the number of water molecules in the first solvation shell. We find the helical propensity of the helical domain of Trp-cage rather well preserved even at very high temperatures. In the folded state, we can identify states with one and two buried internal water molecules interconnecting parts of the Trp-cage molecule by hydrogen bonds. The loss of hydrogen bonds of these buried water molecules in the folded state with increasing temperature is likely to destabilize the folded state at elevated temperatures.     

Hydrophobic tendencies of polar groups as a major force in molecular recognition

Proteins and nucleic acids are able to adopt their native conformation and perform their biological role only in the presence of water with which they actively interact in a mutually modifying way. Traditionally, hydrophobic effect has been considered to be the major factor stabilizing biopolymeric structures. However, solvent reorganization around polar groups is an event thermodynamically more unfavorable than solvent reorganization around nonpolar groups. Consequently, burial of polar groups with formation of complementary solute-solute hydrogen bonds out of contact with water is an energetically favorable process that also provides a major force driving macromolecular association and folding. In contrast to nonpolar groups, polar groups may form their complementary intra- or intersolute hydrogen bonds out of contact with water only provided that an appropriate solute structure has been formed with properly positioned hydrogen bond donors and acceptors. Formation of such structures is disfavored entropically and may not be possible due to steric reasons. However, the interior of a folded protein, alpha-helices and beta-sheets, double helical nucleic acid structures, and protein-ligand interfaces all provide rigid matrices where polar groups may form their complementary hydrogen bonds. For these structures, the inward drive of polar groups represents a considerable stabilizing factor.     

Role of hydration water in protein unfolding

In this paper, following our work on the two-state outer neighbor mixed bonding model of water, it is proposed that polar groups promote the formation of the low density ice Ih-type bonding in their neighborhood, whereas nonpolar groups tend to promote the higher density ice II-type structure. In a protein, because of the large numbers of exposed polar and nonpolar groups, large changes in the neighboring water structure can occur. These changes, of course, depend on whether the protein is in its native or its unfolded state and will be shown here to have a direct impact on the thermodynamics of protein unfolding at both high and low temperatures. For example, it is known that the polar hydration entropies become rapidly more negative with increasing temperature. This very unusual behavior can be directly related to the promotion in the outer bulk liquid of the more stable Ih-type bonding at the expense of II-type bonding by polar groups of the protein. In contrast, nonpolar groups have an opposite effect on the thermodynamics. It is the delicate balance created by these outer hydration contributions, mixed with ordinary thermodynamic contributions from the inner hydration shell and those from hydrogen-bond and van der Waals forces within the protein molecule itself that is responsible for both heat and cold denaturation of proteins.     

Exploring the differences and similarities between urea and thermally driven denaturation of bovine serum albumin: intermolecular forces and solvation preferences

The interactions of bovine serum albumin (BSA) with urea/water were investigated by computer simulation. It was revealed that the BSA-hydrophobic residues in urea solutions favored contact with urea more than with water. Energy decomposition analysis showed that van der Waals energy was the dominant driving force behind urea affinity for hydrophobic residues, whereas coulombic attraction was largely responsible for water affinity for these residues. Meanwhile, urea–BSA hydrogen bond energies were found to be weaker than water–BSA hydrogen bond energies. The greater strength of water–BSA hydrogen bonds than urea–BSA hydrogen bonds, and the opposing preferential interaction between the BSA and urea suggest that disruption of hydrophobic interaction predominates urea–protein denaturation. In pure water, hydrophobic residues showed aggregation tendencies at 323 K, suggesting an increase in hydrophobicity, while at 353 K the residues were partly denatured due to loss of hydrogen bonds; thus, disruption of hydrophobic interactions appeared to contribute less to thermal denaturation.     

Contribution of polar groups in the interior of a protein to the conformational stability

It has been generally believed that polar residues are usually located on the surface of protein structures. However, there are many polar groups in the interior of the structures in reality. To evaluate the contribution of such buried polar groups to the conformational stability of a protein, nonpolar to polar mutations (L8T, A9S, A32S, I56T, I59T, I59S, A92S, V93T, A96S, V99T, and V100T) in the interior of a human lysozyme were examined. The thermodynamic parameters for denaturation were determined using a differential scanning calorimeter, and the crystal structures were analyzed by X-ray crystallography. If a polar group had a heavy energy cost to be buried, a mutant protein would be remarkably destabilized. However, the stability (Delta G) of the Ala to Ser and Val to Thr mutant human lysozymes was comparable to that of the wild-type protein, suggesting a low-energy penalty of buried polar groups. The structural analysis showed that all polar side chains introduced in the mutant proteins were able to find their hydrogen bond partners, which are ubiquitous in protein structures. The empirical structure-based calculation of stability change (Delta Delta G) [Takano et al. (1999) Biochemistry 38, 12698--12708] revealed that the mutant proteins decreased the hydrophobic effect contributing to the stability (Delta G(HP)), but this destabilization was recovered by the hydrogen bonds newly introduced. The present study shows the favorable contribution of polar groups with hydrogen bonds in the interior of protein molecules to the conformational stability.     

Thermodynamic properties of globular proteins and the principle of stabilization of their native structure

A semi-empirical method has been used to estimate the thermodynamic parameters of hydration of buried surface areas of ribonuclease S, lysozyme and myoglobin from the model of complete unfolding according to Ooi et al. ((1987) Proc. Natl. Acad. Sci. USA 84, 3086-3090). The buried surface area of proteins is considered as the difference between the accessible surface area of native protein and the completely extended polypeptide chain according to Lee and Richards ((1971) J. Mol. Biol. 55, 379-400). The contributions of nonpolar and polar protein groups to the general value of Gibbs energy, enthalpy, entropy and heat capacity of hydration have been determined. The obtained results on the thermodynamic behavior of proteins in the process of complete unfolding are in good agreement with the results of microcalorimetric studies of thermal denaturation.     

Collagen denaturation in spread monolayers at the air-water interface: Experiments and a possible model of the process

Spread monolayers of the fibril protein collagen were studied at the air-water interface in the presence of denaturants, urea and thiourea. The most prominent feature of spread collagen monolayers at the air-water interface is the ability to form supramolecular structures (fibrils), which themselves can form monolayers with collapse points of their own. The surface pressure isotherms of collagen monolayers have two “quasi-linear” centers, which are separated by a plateau and correspond to liquid-expanded and liquid-condensed states; this unique capability makes collagen different from other proteins. When in monolayer, collagen acquires the same level of structural organization as in the bulk. In the presence of denaturants, subphase characteristics of collagen monolayers change rapidly and irreversibly. Thiourea exerts more pronounced denaturing action on collagen monolayers than urea; this effect increases with exposure time and denaturant concentration. A hypothetical mechanism of thiourea-induced denaturation of fibril proteins is proposed according to which interactions between hydrophobic C=S groups of thiourea and nonpolar surface groups of the protein lead to reorientation of carbonyl groups to formation of intrinsic hydrogen bonds with NH₂-groups of thiourea eventually resulting in the rupture of intrinsic hydrogen bonds and denaturation of the protein.     

On the edge of the denaturation process: application of X-ray diffraction to barnase and lysozyme cross-linked crystals with denaturants in molar concentrations

Structural data about the early step of protein denaturation were obtained from cross-linked crystals for two small proteins: barnase and lysozyme. Several denaturant agents like urea, bromoethanol or thiourea were used at increasing concentrations up to a limit leading to crystal disruption (>or=2 to 6 M). Before the complete destruction of the crystal order started, specific binding sites were observed at the protein surfaces, an indication that the preliminary step of denaturation is the disproportion of intermolecular polar bonds to the benefit of the agent "parasiting" the surface. The analysis of the thermal factors first agree with a stabilization effect at low or moderate concentration of denaturants rapidly followed by a destabilization at specific weak points when the number of sites increase (overflooding effect).     

Geofold: topology-based protein unfolding pathways capture the effects of engineered disulfides on kinetic stability

Protein unfolding is modeled as an ensemble of pathways, where each step in each pathway is the addition of one topologically possible conformational degree of freedom. Starting with a known protein structure, GeoFold hierarchically partitions (cuts) the native structure into substructures using revolute joints and translations. The energy of each cut and its activation barrier are calculated using buried solvent accessible surface area, side chain entropy, hydrogen bonding, buried cavities, and backbone degrees of freedom. A directed acyclic graph is constructed from the cuts, representing a network of simultaneous equilibria. Finite difference simulations on this graph simulate native unfolding pathways. Experimentally observed changes in the unfolding rates for disulfide mutants of barnase, T4 lysozyme, dihydrofolate reductase, and factor for inversion stimulation were qualitatively reproduced in these simulations. Detailed unfolding pathways for each case explain the effects of changes in the chain topology on the folding energy landscape. GeoFold is a useful tool for the inference of the effects of disulfide engineering on the energy landscape of protein unfolding.     

Understanding the role of hydrogen bonds in water dynamics and protein stability

The mechanisms of cold and pressure denaturation of proteins are a matter of debate, but it is commonly accepted that water plays a fundamental role in the process. It has been proposed that the denaturation process is related to an increase of hydrogen bonds among hydration water molecules. Other theories suggest that the causes of denaturation are the density fluctuations of surface water, or the destabilization of hydrophobic contacts as a consequence of water molecule inclusions inside the protein, especially at high pressures. We review some theories that have been proposed to give insight into this problem, and we describe a coarse-grained model of water that compares well with experiments for proteins’ hydration water. We introduce its extension for a homopolymer in contact with the water monolayer and study it by Monte Carlo simulations in an attempt to understand how the interplay of water cooperativity and interfacial hydrogen bonds affects protein stability.     

Native hydrogen bonds in a molten globule: the apoflavodoxin thermal intermediate

The structure and energetics of protein-folding intermediates are poorly understood. We have identified, in the thermal unfolding of the apoflavodoxin from Anabaena PCC 7119, an equilibrium intermediate with spectroscopic properties of a molten globule and substantial enthalpy and heat capacity of unfolding. The structure of the intermediate is probed by mutagenesis (and phi analysis) of polar residues involved in surface-exposed hydrogen bonds connecting secondary-structure elements in the native protein. All hydrogen bonds analysed are formed in the molten globule intermediate, either with native strength or debilitated. This suggests the overall intermediate's topology and surface tertiary interactions are close to native, and indicates that hydrogen bonding may contribute significantly to shape the conformation and energetics of folding intermediates.