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.     

Spectroscopic studies of the unfolding of a multimeric protein α‐crystallin

α‐Crystallin is a multimeric eye lens protein having molecular chaperone‐like function which is crucial for lens transparency. The stability and unfolding‐refolding properties of α‐crystallin plays important roles for its function. We undertook a multi probe based fluorescence spectroscopic approach to explore the changes in the various levels of organization of this protein at different urea concentration. Steady state fluorescence studies reveal that at 0.6M urea a compact structural intermediate is formed which has a native‐like secondary structure with enhanced surface exposure of hydrophobic groups. At 2.8M urea the tertiary interactions are largely collapsed with partial collapse of secondary and quaternary structure. The surface solvation probed by picosecond time resolved fluorescence of acrylodan labeled α‐crystallin revealed dry native‐like core of α‐crystallin at 0.6M urea compared to enhanced water penetration at 2.8M urea and extensive solvation at 6M urea. Activation energy for the subunit exchange decreased by 22 kJ mol^?1 on changing urea concentration from 0 to 0.6M compared with over 75 kJ mol^?1 on changing urea concentration from 0 to 2.8M. Light scattering and analytical ultracentrifugation techniques were used to determine size and oligomerization of the unfolding intermediates. The data indicated swelling but no oligomer breakdown at 0.6M urea. At 2.8M urea the oligomeric size is considerably reduced and a monomer is produced at 6M urea. The data clearly reveals that structural breakdown of α‐crystallin does not follow hierarchical sequence as tertiary structure dissolution takes place before complete oligomeric dissociation. © 2013 Wiley Periodicals, Inc. Biopolymers 101: 549–560, 2014.     

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.     

Protein molecular dynamics with electrostatic force entirely determined by a single Poisson-Boltzmann calculation

The electrostatic force including the intramolecular Coulombic interactions and the electrostatic contribution of solvation effect were entirely calculated by using the finite difference Poisson-Boltzmann method (FDPB), which was incorporated into the GROMOS96 force field to complete a new finite difference stochastic dynamics procedure (FDSD). Simulations were performed on an insulin dimer. Different relative dielectric constants were successively assigned to the protein interior; a value of 17 was selected as optimal for our system. The simulation data were analyzed and compared with those obtained from 500-ps molecular dynamics (MD) simulation with explicit water and a 500-ps conventional stochastic dynamics (SD) simulation without the mean solvent force. The results indicate that the FDSD method with GROMOS96 force field is suitable to study the dynamics and structure of proteins in solution if used with the optimal protein dielectric constant.     

Comparison of Guanidine Hydrochloride (GdnHCl) and Urea Denaturation on Inactivation and Unfolding of Human Placental Cystatin (HPC)

The activity and conformational change of human placental cystatin (HPC), a low molecular weight thiol proteinase inhibitor (12,500) has been investigated in presence of guanidine hydrochloride (GdnHCl) and urea. The denaturation of HPC was followed by activity measurements, fluorescence spectroscopy and Circular Dichroism (CD) studies. Increasing the denaturant concentration significantly enhanced the inactivation and unfolding of HPC. The enzyme was 50% inactivated at 1.5 M GdnHCl or 3 M urea. Up to 1.5 M GdnHCl concentration there was quenching of fluorescence intensity compared to native form however at 2 M concentration intensity increased and emission maxima had 5 nm red shift with complete unfolding in 4–6 M range. The mid point of transition was in the region of 1.5–2 M. In case of urea denaturation, the fluorescence intensity increased gradually with increase in the concentration of denaturant. The protein unfolded completely in 6–8 M concentration of urea with a mid-point of transition at 3 M. CD spectroscopy shows that the ellipticity of HPC has increased compared to that of native up to 1.5 M GdnHCl and then there is gradual decrease in ellipticity from 2 to 5 M concentration. At 6 M GdnHCl the protein had random coil conformation. For urea the ellipticity decreases with increase in concentration showing a sigmoidal shaped transition curve with little change up to 1 M urea. The protein greatly loses its structure at 6 M urea and at 8 M it is a random coil. The urea induced denaturation follows two-state rule in which Native→Denatured state transition occurs in a single step whereas in case of GdnHCl, intermediates or non-native states are observed at lower concentrations of denaturant. These intermediate states are possibly due to stabilizing properties of guanidine cation (Gdn⁺) at lower concentrations, whereas at higher concentrations it acts as a classical denaturant.     

Microsecond simulations of the folding/unfolding thermodynamics of the Trp‐cage miniprotein

We study the unbiased folding/unfolding thermodynamics of the Trp‐cage miniprotein using detailed molecular dynamics simulations of an all‐atom model of the protein in explicit solvent using the Amberff99SB force field. Replica‐exchange molecular dynamics simulations are used to sample the protein ensembles over a broad range of temperatures covering the folded and unfolded states at two densities. The obtained ensembles are shown to reach equilibrium in the 1 μs/replica timescale. The total simulation time used in the calculations exceeds 100 μs. Ensemble averages of the fraction folded, pressure, and energy differences between the folded and unfolded states as a function of temperature are used to model the free energy of the folding transition, ΔG(P, T), over the whole region of temperatures and pressures sampled in the simulations. The ΔG(P, T) diagram describes an ellipse over the range of temperatures and pressures sampled, predicting that the system can undergo pressure‐induced unfolding and cold denaturation at low temperatures and high pressures, and unfolding at low pressures and high temperatures. The calculated free energy function exhibits remarkably good agreement with the experimental folding transition temperature (T_f = 321 K), free energy, and specific heat changes. However, changes in enthalpy and entropy are significantly different than the experimental values. We speculate that these differences may be due to the simplicity of the semiempirical force field used in the simulations and that more elaborate force fields may be required to describe appropriately the thermodynamics of proteins. Proteins 2010. © 2010 Wiley‐Liss, Inc.     

Biophysical characterization and unfolding of LEF4 factor of RNA polymerase from AcNPV

Late expression factor 4 (LEF4) is one of the four subunits of Autographa californica nuclear polyhedrosis virus (AcNPV) RNA polymerase. LEF4 was overexpressed in Escherichia coli and recombinant protein was subjected to structural characterization. Chemical induced unfolding of LEF4 was investigated using intrinsic fluorescence, hydrophobic dye binding, fluorescence quenching, and circular dichroism (CD) techniques. The unfolding of LEF4 was found to be a non‐two state, biphasic transition. Intermediate states of LEF4 at 2M GnHCl and 4M urea shared some common structural features and hence may lie on the same pathway of protein folding. Steady‐state fluorescence and far‐UV CD showed that while there was considerable shift in the wavelength of emission maximum (λ_max), the secondary structure of LEF4 intermediates at 2M GnHCl and 4M urea remained intact. Further, temperature induced denaturation of LEF4 was monitored using far‐UV CD. This study points to the structural stability of LEF4 under the influence of denaturants like urea and temperature. Although LEF4 is an interesting model protein to study protein folding intermediates, in terms of functional significance the robust nature of this protein might reflect one of the several strategies adapted by the virus to survive under very adverse environmental and physiological conditions. © 2009 Wiley Periodicals, Inc. Biopolymers 91: 574–582, 2009. This article was originally published online as an accepted preprint. The “Published Online” date corresponds to the preprint version. You can request a copy of the preprint by emailing the Biopolymers editorial office at biopolymers@wiley.com     

Dielectric study of the urea denaturation of delipidated and relipidated bovine serum albumin

Dielectric relaxation and viscosity measurements were performed on delipidated and relipidated samples of bovine serum albumin (BSA) at urea concentrations between O and 6M. By the combined interpretation of these two hydrodynamic methods the characterization of conformational changes of the molecule during urea denaturation is possible. The denaturation of delipidated BSA results from two mechanisms. The first one is a slow, time-dependent elongation of the molecule; the second one is a rapid swelling which becomes most pronounced at urea concentrations higher than 4M. For relipidated albumin, the slow elongation mechanism occurs but the presence of fatty acids protects the protein aganist molecular swelling. In both cases these conformational changes are accompanied by an increased disymmetry of charge repartition and a concomitant increase of the dipole moment. From these results it follows that lipidated albumin (as occurs under physiological conditions) is less sensitive to denaturation than delipidated albumin.     

Thermodynamics of the Trp‐cage miniprotein unfolding in urea

The thermodynamic properties of unfolding of the Trp‐cage mini protein in the presence of various concentrations of urea have been characterized using temperature‐induced unfolding monitored by far‐UV circular dichroism spectroscopy. Analysis of the data using a two‐state model allowed the calculation of the Gibbs energy of unfolding at 25°C as a function of urea concentration. This in turn was analyzed by the linear extrapolation model that yielded the dependence of Gibbs energy on urea concentration, i.e. the m‐value for Trp‐cage unfolding. The m‐value obtained from the experimental data, as well as the experimental heat capacity change upon unfolding, were correlated with the structural parameters derived from the three dimensional structure of Trp‐cage. It is shown that the m‐value can be predicted well using a transfer model, while the heat capacity changes are in very good agreement with the empirical models based on model compounds studies. These results provide direct evidence that Trp‐cage, despite its small size, is an excellent model for studies of protein unfolding and provide thermodynamic data that can be used to compare with atomistic computer simulations. Proteins 2010. © 2009 Wiley‐Liss, Inc.     

Beyond Lumry-Eyring: an unexpected pattern of operational reversibility/irreversibility in protein denaturation

We have found that, contrary to naïve intuition, the degree of operational reversibility in the thermal denaturation of lipase from Thermomyces lanuginosa (an important industrial enzyme) in urea solutions is maximum when the protein is heated several degrees above the end of the temperature‐induced denaturation transition. Upon cooling to room temperature, the protein seems to reach a state with enzymatic activity similar to that of the initial native state, but with higher denaturation temperature and radically different behavior in terms of susceptibility to irreversible denaturation. These results show that patterns of operational reversibility/irreversibility in protein denaturation may be more complex than the often‐taken‐for‐granted, two‐situation classification (reversible vs. irreversible). Furthermore, they are consistent with the possibility of existence of different native or native‐like states separated by high kinetic barriers under native conditions and they suggest experimental procedures to reach and study such “alternative” native states. Proteins 2008. © 2007 Wiley‐Liss, Inc.     

Temperature and pressure denaturation of chignolin: Folding and unfolding simulation by multibaric-multithermal molecular dynamics method

A multibaric‐multithermal molecular dynamics (MD) simulation of a 10‐residue protein, chignolin, was performed. All‐atom model with the Amber parm99SB force field was used for the protein and the TIP3P model was used for the explicit water molecules. This MD simulation covered wide ranges of temperature between 260 and 560 K and pressure between 0.1 and 600 MPa and sampled many conformations without getting trapped in local‐minimum free‐energy states. Folding events to the native β‐hairpin structure occurred five times and unfolding events were observed four times. As the temperature and/or pressure increases, fraction of folded chignolin decreases. The partial molar enthalpy change ΔH and partial molar volume change ΔV of unfolding were calculated as ΔH = 24.1 ± 4.9 kJ/mol and ΔV = ?5.6 ± 1.5 cm³/mol, respectively. These values agree well with recent experimental results. Illustrating typical local‐minimum free‐energy conformations, folding and unfolding pathways were revealed. When chignolin unfolds from the β‐hairpin structure, only the C terminus or both C and N termini open first. It may undergo an α‐helix or 3₁₀‐helix structure and finally unfolds to the extended structure. Difference of the mechanism between temperature denaturation and pressure denaturation is also discussed. Temperature denaturation is caused by making the protein transferred to a higher entropy state and making it move around more with larger space. The reason for pressure denaturation is that water molecules approach the hydrophobic residues, which are not well hydrated at the folded state, and some hydrophobic contacts are broken. Proteins 2012;. © 2012 Wiley Periodicals, Inc.     

Tryptophan fluorescence studies of plasma fibronectin: effects of environmental factors

Steady state fluorescence measurements have been used to study tryptophan fluorescence of plasma fibronectin. The native protein has an emission maximum at 337 nm with a quantum yield of 0.03. A red shift of emission maximum was observed in 3–5M urea and a further red shift in 7–8M urea. The emission maximum shifted from 337 to 345 nm when the temperature was changed from 30 to 80°C, with a midpoint of thermal denaturation at 58°C. Similarly, the emission maximum shifted from 337 to 345 nm when the solution pH was increased from 9 to 12, with a midpoint of pH transition at 10.6. The results obtained from difference absorption spectroscopy studies suggest that the unfolding of fibronectin at alkaline pH is related at least in part to ionization of tyrosine residues. Since most of the tryptophan residues are in invariant positions in homology sequences, it is suggested here that tryptophan residues are useful intrinsic probes for elucidating fibronectin structure in solution.     

Can current force fields reproduce ring puckering in 2-O-sulfo-α-l-iduronic acid? A molecular dynamics simulation study

The monosaccharide 2-O-sulfo-α-l-iduronic acid (IdoA2S) is one of the major components of glycosaminoglycans. The ability of molecular mechanics force fields to reproduce ring-puckering conformational equilibrium is important for the successful prediction of the free energies of interaction of these carbohydrates with proteins. Here we report unconstrained molecular dynamics simulations of IdoA2S monosaccharide that were carried out to investigate the ability of commonly used force fields to reproduce its ring conformational flexibility in aqueous solution. In particular, the distribution of ring conformer populations of IdoA2S was determined. The GROMOS96 force field with the SPC/E water potential can predict successfully the dominant skew-boat to chair conformational transition of the IdoA2S monosaccharide in aqueous solution. On the other hand, the GLYCAM06 force field with the TIP3P water potential sampled transitional conformations between the boat and chair forms. Simulations using the GROMOS96 force field showed no pseudorotational equilibrium fluctuations and hence no inter-conversion between the boat and twist boat ring conformers. Calculations of theoretical proton NMR coupling constants showed that the GROMOS96 force field can predict the skew-boat to chair conformational ratio in good agreement with the experiment, whereas GLYCAM06 shows worse agreement. The omega rotamer distribution about the C5-C6 bond was predicted by both force fields to have torsions around 10°, 190°, and 360°.     

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.     

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.     

Why the OPLS-AA force field cannot produce the β-hairpin structure of H1 peptide in solution when comparing with the GROMOS 43A1 force field?

The optimal combination of force field and water model is an essential problem that is able to increase molecular dynamics simulation quality for different types of proteins and peptides. In this work, an attempt has been made to explore the problem by studying H1 peptide using four different models based on different force fields, water models and electrostatic schemes. The driving force for H1 peptide conformation transition and the reason why the OPLS-AA force field cannot produce the β-hairpin structure of H1 peptide in solution while the GROMOS 43A1 force field can do were investigated by temperature replica exchange molecular dynamics simulation (T-REMD). The simulation using the GROMOS 43A1 force field preferred to adopt a β-hairpin structure, which was in good agreement with the several other simulations and the experimental evidences. However, the simulation using the OPLS-AA force field has a significant difference from the simulations with the GROMOS 43A1 force field simulation. The results show that the driving force in H1 peptide conformation transition is solvent exposure of its hydrophobic residues. However, the subtle balances between residue-residue interactions and residue-solvent interaction are disrupted by using the OPLS-AA force field, which induced the reduction in the number of residue-residue contact. Similar solvent exposure of the hydrophobic residues is observed for all the conformations sampled using the OPLS-AA force field. For H1 peptide which exhibits large solvent exposure of the hydrophobic residues, the GROMOS 43A1 force field with the SPC water model can provide more accurate results.     

Molecular dynamics simulation of hen egg white lysozyme: a test of the GROMOS96 force field against nuclear magnetic resonance data

Biomolecular force fields for use in molecular dynamics (MD) simulations of proteins, DNA, or membranes are generally parametrized against ab initio quantum-chemical and experimental data for small molecules. The application of a force field in a simulation of a biomolecular system, such as a protein in solution, may then serve as a test of the quality and transferability of the force field. Here, we compare various properties obtained from two MD simulations of the protein hen egg white lysozyme (HEWL) in aqueous solution using the latest version, GROMOS96, of the GROMOS force field and an earlier version, GROMOS87+, with data derived from nuclear magnetic resonance (NMR) experiments: NOE atom-atom distance bounds, (3)J(HNalpha)-coupling constants, and backbone and side-chain order parameters. The convergence of these quantities over a 2-ns period is considered, and converged values are compared to experimental ones. The GROMOS96 simulation shows better agreement with the NMR data and also with the X-ray crystal structure of HEWL than the GROMOS87+ simulation, which was based on an earlier version of the GROMOS force field.     

Conformational space comparison of GnRH and lGnRH-III using molecular dynamics, cluster analysis and Monte Carlo thermodynamic integration

The conformational space available to GnRH and lGnRH-III was compared using 5.2 ns constant temperature and pressure molecular dynamics simulations with explicit TIP3P solvation and the AMBER v. 5.0 force field. Cluster analysis of both trajectories resulted in two groups of conformations. Results of free energy calculations, in agreement with previous experimental data, indicate that a conformation with a turn from residues 5 through 8 is preferred for GnRH in an aqueous environment. By contrast, a conformation with a helix from residues 2 through 7 with a bend from residues 6 through 10 is preferred for lGnRH-III in an aqueous environment. The side chains of His2 and Trp3 in lGnRH-III occupy different regions of phase space and participate in weakly polar interactions different from those in GnRH. The unique conformational properties of lGnRH-III may account for its specific anti cancer activity.     

Engineering out motion: introduction of a de novo disulfide bond and a salt bridge designed to close a dynamic cleft on the surface of cytochrome b5

A previous molecular dynamics (MD) simulation of cytochrome b5 (cyt b5) at 25 degrees C displayed localized dynamics on the surface of the protein giving rise to the periodic formation of a cleft that provides access to the heme through a protected hydrophobic channel [Storch and Daggett (1995) Biochemistry 34, 9682]. Here we describe the production and testing of mutants designed to prevent the cleft from opening using a combination of experimental and theoretical techniques. Two mutants have been designed to close the surface cleft: S18D to introduce a salt bridge and S18C:R47C to incorporate a disulfide bond. The putative cleft forms between two separate cores of the protein: one is structural in nature and can be monitored through the fluorescence of Trp 22, and the other binds the heme prosthetic group and can be tracked via heme absorbance. An increase in motion localized to the cleft region was observed for each protein, except for the disulfide-containing variant, in MD simulations at 50 degrees C compared to simulations at 25 degrees C. For the disulfide-containing variant, the cleft remained closed. Both urea and temperature denaturation curves were nearly identical for wild-type and mutant proteins when heme absorbance was monitored. In contrast, fluorescence studies revealed oxidized S18C:R47C to be considerably more stable based on the midpoints of the denaturation transitions, Tm and U1/2. Moreover, the fluorescence changes for each protein were complete at approximately 50 degrees C and a urea concentration of approximately 3.9 M, significantly below the temperature and urea concentration (62 degrees C, 5 M urea) required to observe heme release. In addition, solvent accessibility based on acrylamide quenching of Trp 22 was lower in the S18C:R47C mutant, particularly at 50 degrees C, before heme release [presented in the accompanying paper (58)]. The results suggest that a constraining disulfide bond can be designed to inhibit dynamic cleft formation on the surface of cyt b5. Located near the heme, the native dynamics of the cleft may be functionally important for protein-protein recognition and/or complex stabilization.