Structural details of urea binding to barnase: a molecular dynamics analysis.

BACKGROUND: The molecular mechanism of urea-induced protein unfolding has not been established. It is generally thought that denaturation results from the stabilizing interactions of urea with portions of the protein that are buried in the native state and become exposed upon unfolding of the protein. RESULTS: We have performed molecular dynamics simulations of barnase (a 110 amino acid RNase from Bacillus amyloliquefaciens) with explicit water and urea molecules at 300 K and 360 K. The native conformation was unaffected in the 300 K simulations at neutral and low pH. Two of the three runs at 360 K and low pH showed some denaturation, with partial unfolding of the hydrophobic core 2. The first solvation shell has a much higher density of urea molecules (water/urea ratio ranging from 2.07 to 2.73) than the bulk (water/urea ratio of 4.56). About one half of the first-shell urea molecules are involved in hydrogen bonds with polar or charged groups on the barnase surface, and between 15% and 18% of the first-shell urea molecules participate in multiple hydrogen bonds with barnase. The more stably bound urea molecules tend to be in crevices or pockets on the barnase surface. CONCLUSIONS: The simulation results indicate that an aqueous urea solution solvates the surface of a polypeptide chain more favorably than pure water. Urea molecules interact more favorably with nonpolar groups of the protein than water does, and the presence of urea improves the interactions of water molecules with the hydrophilic groups of the protein. The results suggest that urea denaturation involves effects on both nonpolar and polar groups of proteins.     

Cold denaturation of alpha-lactalbumin

We have investigated the thermal unfolding of bovine alpha-lactalbumin by means of circular dichroism spectroscopy in the far- and near-ultraviolet regions, and shown that the native alpha-lactalbumin undergoes heat and cold denaturation. The guanidine hydrochloride-induced unfolding of alpha-lactalbumin was also investigated by circular dichroism spectroscopy at various temperatures from 261 to 318 K. It is shown that the population of the molten globule state is strongly dependent on temperature and that the molten globule state does not accumulate during the guanidine hydrochloride-induced unfolding transition at 261 K. Our results indicate that the molten globule state of alpha-lactalbumin undergoes cold denaturation as the native alpha-lactalbumin does, and that the heat capacity change of unfolding from the molten globule to the unfolded state is positive and significant. The present results further support the idea that the molten globule and the unfolded states do not belong to the same thermodynamic state, and that the native, molten globule and unfolded states are sufficient for interpreting the guanidine hydrochloride-induced unfolding behavior of alpha-lactalbumin.     

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.     

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.     

Temperature and urea induced denaturation of the TRP‐cage mini protein TC5b: A simulation study consistent with experimental observations

The effects of temperature and urea denaturation (6M urea) on the dominant structures of the 20‐residue Trp‐cage mini‐protein TC5b are investigated by molecular dynamics simulations of the protein at different temperatures in aqueous and in 6M urea solution using explicit solvent degrees of freedom and the GROMOS force‐field parameter set 45A3. In aqueous solution at 278 K, TC5b is stable throughout the 20 ns of MD simulation and the trajectory structures largely agree with the NMR‐NOE atom–atom distance data available. Raising the temperature to 360 K and to 400 K, the protein denatures within 22 ns and 3 ns, showing that the denaturation temperature is well below 360 K using the GROMOS force field. This is 40–90 K lower than the denaturation temperatures observed in simulations using other much used protein force fields. As the experimental denaturation temperature is about 315 K, the GROMOS force field appears not to overstabilize TC5b, as other force fields and the use of continuum solvation models seem to do. This feature may directly stem from the GROMOS force‐field parameter calibration protocol, which primarily involves reproduction of condensed phase thermodynamic quantities such as energies, densities, and solvation free energies of small compounds representative for protein fragments. By adding 6M urea to the solution, the onset of denaturation is observed in the simulation, but is too slow to observe a particular side‐chain side‐chain contact (Trp6‐Ile4) that was experimentally observed to be characteristic for the denatured state. Interestingly, using temperature denaturation, the process is accelerated and the experimental data are reproduced.     

Structural characterization of the molten globule of alpha-lactalbumin by solution X-ray scattering.

A compact denatured state is often observed under a mild denaturation condition for various proteins. A typical example is the alpha-lactalbumin molten globule. Although the molecular compactness and shape are the essential properties for defining the molten globule, there have been ambiguities of these properties for the molten globule of alpha-lactalbumin. Using solution X-ray scattering, we have examined the structural properties of two types of molten globule of alpha-lactalbumin, the apo-protein at neutral pH and the acid molten globule. The radius of gyration for the native holo-protein was 15.7 A, but the two different molten globules both had a radius of gyration of 17.2 A. The maximum dimension of the molecule was also increased from 50 A for the native state to 60 A for the molten globule. These values clearly indicate that the molten globule is not as compact as the native state. The increment in the radius of gyration was less than 10% for the alpha-lactalbumin molten globule, compared with up to 30% for the molten globules of other globular proteins. Intramolecular disulfide bonds restrict the molecular expansion of the molten globule. The distance distribution function of the alpha-lactalbumin molten globule is composed of a single peak suggesting a globular shape, which is simply swollen from the native state. The scattering profile in the high Q region of the molten globule indicates the presence of a significant amount of tertiary fold. Based on the structural properties obtained by solution X-ray scattering, general and conceptual structural images for the molten globules of various proteins are described and compared with the individual, detailed structural model obtained by nuclear magnetic resonance.     

Spin diffusion in globular proteins: alpha-lactalbumin

The structural properties of globular proteins analyzed by two different methods: high-resolution NMR and circular dichroism were compared. We established that the spin diffusion method shows changes in the secondary structure during the unfolding of the alpha-lactalbumin molten globule by urea. It was shown that the spin diffusion method is extremely effective in studies of interactions of water and denaturant molecules with the protein both in the native and the molten globule states.     

1,1,1,3,3,3-hexafluoroisopropanol induced thermal unfolding and molten globule state of bovine alpha-lactalbumin: calorimetric and spectroscopic studies

The thermal denaturation of alpha-lactalbumin was studied at pH 7.0 and 9.0 in aqueous 1,1,1,3,3,3-hexafluoroisopropanol (HFIP) by high-sensitivity differential scanning calorimetry. The conformation of the protein was analyzed by a combination of fluorescence and circular dichroism measurements. The most obvious effect of HFIP was lowering of the transition temperature with an increase in the concentration of the alcohol up to 0.30M, beyond which no calorimetric transition was observed. Up to 0.30M HFIP the calorimetric and van't Hoff enthalpy remained the same, indicating the validity of the two-state approximation for the thermal unfolding of alpha-lactalbumin. The quantitative thermodynamic parameters accompanying the thermal transitions have been evaluated. Spectroscopic observations confirm that alpha-lactalbumin is in the molten globule state in the presence of 0.50M HFIP at pH 7.0 and 0.75M HFIP at pH 9.0. The results also demonstrate that alpha-lactalbumin in the molten globule state undergoes a noncooperative thermal transition to the denatured state. It is observed that two of four tryptophans are exposed to the solvent in the HFIP induced molten globule state of alpha-lactalbumin compared to four in the 8.5M urea induced denatured state of the protein. It is also observed that the HFIP induced molten globule states at the two pH values are different from the acid induced molten globule state (A state) of alpha-lactalbumin.     

The structure of denatured alpha-lactalbumin elucidated by the technique of disulfide scrambling: fractionation of conformational isomers of alpha-lactalbumin

The structure of denatured alpha-lactalbumin (alpha-LA) has been characterized using the method of disulfide scrambling. Under denaturing conditions (urea, guanidine hydrochloride, guanidine thiocyanate, organic solvent or elevated temperature) and in the presence of thiol initiator, alpha-LA denatures by shuffling its four native disulfide bonds and converts to a mixture of fully oxidized scrambled structures. Analysis by reversed-phase HPLC reveals that the denatured alpha-LA comprises a minimum of 45 fractions of scrambled isomers. Among them, six well populated isomers have been isolated and structurally characterized. Their relative concentrations, which represent the fingerprinting of the denatured alpha-LA, vary substantially under different denaturing conditions. These results permit independent plotting of the denaturation and unfolding curves of alpha-LA. Most importantly, unique isomers of partially unfolded alpha-LA were shown to populate at mild and selected denaturing conditions. Organic solvent disrupts preferentially the hydrophobic alpha-helical domain, generating a predominant isomer containing two native disulfide bonds at the beta-sheet domain and two scrambled disulfide bonds at the alpha-helical region. Thermal denaturation selectively unfolds the beta-sheet domain of alpha-LA, producing a prevalent isomer that exhibits structural characteristics of the molten globule state of alpha-LA.     

Entropy calculations on the molten globule state of a protein: side-chain entropies of alpha-lactalbumin.

We present entropy estimates based on molecular dynamics simulations of models of the molten globule state of the protein alpha-lactalbumin at low pH. The entropy calculations use the covariance matrix of atom-positional fluctuations and yield the complete configurational entropy. The configurational entropy of the entire protein and of each of its side chains is calculated. Exposed side chains show a larger entropy compared to buried side chains. A comparison to data from rotamer counting is made and significant differences are found.     

Increased exposure of hydrophobic surface in molten globule state of alpha-lactalbumin. Fluorescence and hydrophobic photolabeling studies.

The involvement of molten globule state as a distinct intermediate in the denaturation process in proteins is well documented. However, the structural characterization of such an intermediate is far from complete. We have, using fluorescence and fluorescence quenching, studied the molten globule state of bovine alpha-lactalbumin. Unlike the native state, where all the 4 tryptophans are buried in the protein, 2 tryptophans are exposed in the molten globule state. Using the hydrophobic photoactivable reagent [3H]diazofluorene, we observe an increased hydrophobic exposure in the molten globule state. These structural characteristics conform to the current views on the molten globule state, i.e. it has similar secondary structure but a poorly defined tertiary structure. Our fluorescence studies indicate the involvement of a premolten globule state in the native to molten globule state transition. This premolten globule state exists at pH 5.0 and has a very compact structure involving increased hydrophobic interactions in the protein interior. These results are also supported by circular dichroism studies.     

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.     

A stabilized molten globule protein

A predominant conformational isomer of non-native alpha-lactalbumin (alpha-LA) has been purified by thermal denaturation of the native alpha-LA using the technique of disulfide scrambling. This unique isomer retains a substantial content of alpha-helical structure. It is stabilized by two native disulfide bonds within the alpha-helical domain and two scrambled non-native disulfide bonds at the beta-sheet domain. This denatured isomer of alpha-LA exhibits structural characteristics that are consistent with the well-documented molten globule state. The ability to prepare a stabilized and structurally defined molten globule provides a useful model for studying the folding and unfolding pathways of proteins.     

Effects of the stabilization of the molten globule state on the folding mechanism of alpha-lactalbumin: a study of a chimera of bovine and human alpha-lactalbumin

The N-terminal half of the alpha-domain (residues 1 to 34) is more important for the stability of the acid-induced molten globule state of alpha-lactalbumin than the C-terminal half (residues 86 to 123). The refolding and unfolding kinetics of a chimera, in which the amino acid sequence of residues 1 to 34 was from human alpha-lactalbumin and the remainder of the sequence from bovine alpha-lactalbumin, were studied by stopped-flow tryptophan fluorescence spectroscopy. The chimeric protein refolded and unfolded substantially faster than bovine alpha-lactalbumin. The stability of the molten globule state formed by the chimera was greater than that of bovine alpha-lactalbumin, and the hydrophobic surface area buried inside of the molecule in the molten globule state was increased by the substitution of residues 1 to 34. Peptide fragments corresponding to the A- and B-helix of the chimera showed higher helix propensity than those of the bovine protein, indicating the contribution of local interactions to the high stability of the molten globule state of the chimera. Moreover, the substitution of residues 1-34 decreased the free energy level of the transition state and increased hydrophobic surface area buried inside of the molecule in the transition state. Our results indicate that local interactions as well as hydrophobic interactions formed in the molten globule state are important in guiding the subsequent structural formation of alpha-lactalbumin.     

Defining the core structure of the alpha-lactalbumin molten globule state.

Molten globules are partially folded states of proteins which are generally believed to mimic structures formed during the folding process. In order to determine the minimal requirements for the formation of a molten globule state, we have prepared a set of peptide models of the molten globule state of human alpha-lactalbumin (alphaLA). A peptide consisting of residues 1-38 crosslinked, via the native 28-111 disulfide bond, to a peptide corresponding to residues 95-120 forms a partially folded state at pH 2.8 which has all of the characteristics of the molten globule state of alphaLA as judged by near and far UV CD, fluorescence, ANS binding and urea denaturation experiments. The structure of the peptide construct is the same at pH 7.0. Deletion of residues 95-100 from the construct has little effect. Thus, less than half the sequence is required to form a molten globule. Further truncation corresponding to the selective deletion of the A (residues 1-19) or D (residues 101-110) helices or the C-terminal 310 helix (residues 112-120) leads to a significant loss of structure. The loss of structure which results from the deletion of any of these three regions is much greater than that which would be expected based upon the non-cooperative loss of local helical structure. Deletion of residues corresponding to the region of the D helix or C-terminal 310 helix region results in a peptide construct which is largely unfolded and contains no more helical structure than is expected from the sum of the helicity of the two reduced peptides. These experiments have defined the minimum core structure of the alphaLA molten globule state.     

Energetic basis of structural stability in the molten globule state: alpha-lactalbumin

The denatured states of alpha-lactalbumin, which have features of a molten globule state, have been studied to elucidate the energetics of the molten globule state and its contribution to the stability of the native conformation. Analysis of calorimetric and CD data shows that the heat capacity increment of alpha-lactalbumin denaturation highly correlates with the degree of disorder of the residual structure of the state. As a result, the denaturational transition of alpha-lactalbumin from the native to a highly ordered compact denatured state, and from the native to the disordered unfolded state are described by different thermodynamic functions. The enthalpy and entropy of the denaturation of alpha-lactalbumin to compact denatured state are always greater than the enthalpy and entropy of its unfolding. This difference represents the unfolding of the molten globule state. Calorimetric measurements of the heat effect associated with the unfolding of the molten globule state reveal that it is negative in sign over the temperature range of molten globule stability. This observation demonstrates the energetic specificity of the molten globule state, which, in contrast to a protein with unique tertiary structure, is stabilized by the dominance of negative entropy and enthalpy of hydration over the positive conformational entropy and enthalpy of internal interactions. It is concluded that at physiological temperatures the entropy of dehydration is the dominant factor providing stability for the compact intermediate state on the folding pathway, while for the stability of the native state, the conformational enthalpy is the dominant factor.     

The molten globule state of a chimera of human alpha-lactalbumin and equine lysozyme.

The molten globule state of equine lysozyme is more stable than that of alpha-lactalbumin and is stabilized by non-specific hydrophobic interactions and native-like hydrophobic interactions. We constructed a chimeric protein which is produced by replacing the flexible loop (residues 105-110) in human alpha-lactalbumin with the helix D (residues 109-114) in equine lysozyme to investigate the possible role of the helix D for the high stability and native-like packing interaction in the molten globule state of equine lysozyme. The stability of the molten globule state formed by the chimeric protein to guanidine hydrochloride-induced unfolding is the same as that of equine lysozyme and is substantially greater than that of human alpha-lactalbumin, although only six residues come from equine lysozyme. Our results also suggest that the non-native interaction in the molten globule state of alpha-lactalbumin changes to the native-like packing interaction due to helix substitution. The solvent-accessibility of the Trp residues in the molten globule state of the chimeric protein is similar to that in the molten globule state of equine lysozyme in which packing interaction around the Trp residues in the native state is partially preserved. Therefore, the helix D in equine lysozyme is one of the contributing factors to the high stability and native-like packing interaction in the molten globule state of equine lysozyme. Our results indicate that the native-like packing interaction can stabilize the rudimentary intermediate which is stabilized by the non-specific hydrophobic interactions.     

Differences in the pathways of proteins unfolding induced by urea and guanidine hydrochloride: molten globule state and aggregates

It was shown that at low concentrations guanidine hydrochloride (GdnHCl) can cause aggregation of proteins in partially folded state and that fluorescent dye 1-anilinonaphthalene-8-sulfonic acid (ANS) binds with these aggregates rather than with hydrophobic clusters on the surface of protein in molten globule state. That is why the increase in ANS fluorescence intensity is often recorded in the pathway of protein denaturation by GdnHCl, but not by urea. So what was previously believed to be the molten globule state in the pathway of protein denaturation by GdnHCl, in reality, for some proteins represents the aggregates of partially folded molecules.     

Hydrogen exchange study of canine milk lysozyme: stabilization mechanism of the molten globule

The native state (1)H, (15)N resonance assignment of 123 of the 128 nonproline residues of canine milk lysozyme has enabled measurements of the amide hydrogen exchange of over 70 amide hydrogens in the molten globule state. To elucidate the mechanism of protein folding, the molten globule state has been studied as a model of the folding intermediate state. Lysozyme and alpha-lactalbumin are homologous to each other, but their equilibrium unfolding mechanisms differ. Generally, the folding mechanism of lysozyme obeys a two-state model, whereas that of alpha-lactalbumin follows a three-state model. Exceptions to this rule are equine and canine milk lysozymes, which exhibit a partially unfolded state during the equilibrium unfolding; this state resembles the molten globule state of alpha-lactalbumin but with extreme stability. Study of the molten globules of alpha-lactalbumin and equine milk lysozyme showed that the stabilities of their alpha-helices are similar, despite the differences in the thermodynamic stability of their molten globule states. On the other hand, our hydrogen exchange study of the molten globule of canine milk lysozyme showed that the alpha-helices are more stabilized than in alpha-lactalbumin or equine milk lysozyme and that this enhanced stability is caused by the strengthened cooperative interaction between secondary structure elements. Thus, our results underscore the importance of the cooperative interaction in the stability of the molten globule state.     

Calorimetric determination of the energetics of the molten globule intermediate in protein folding: apo-alpha-lactalbumin

High-sensitivity differential scanning calorimetry has been used to characterize the energetics of the molten globule state of apo-alpha-lactalbumin. This characterization has been possible by performing temperature scans at different guanidine hydrochloride (GuHCl) concentrations in order to experimentally define the temperature-GuHCl stability surface of the protein. Multidimensional analysis of the heat capacity surface has allowed simultaneous resolution of the energetics of the unfolded and molten globule states. These experiments indicate that the intrinsic enthalpy difference (i.e., excluding additional contributions such as those arising from differential GuHCl binding) between the unfolded and native states is 31.8 kcal/mol at 25 degrees C whereas that of the molten globule and native states is only 7.7 kcal/mol. At the same temperature, the entropy changes are 99.2 and 23.7 cal/K.mol and the heat capacity changes are 1821 and 326 cal/K.mol, respectively. Analysis of the thermodynamic data indicates that in passing from the native to the molten globule state only approximately 19% of the hydrogen bonds are broken. In addition, the magnitude of delta Cp for the molten globule suggests that water does not largely penetrate into the interior of the molten globule, implying that significant hydrophobic interactions are still present in this state. These parameters provide precise energetic constraints to the allowed structural conformations of the molten globule.