Impacts of Ca2+ cation and temperature on bovine α-lactalbumin secondary structures and foamability - Insights from computational molecular dynamics

We used computational molecular dynamics (MD) to assess molecular conformations of apo- and holo-forms (respectively without and with Ca²⁺) of bovine α-lactalbumin (α-La) at different temperatures, and to correlate them with the protein's foaming properties. At 4 °C and 25 °C no major protein conformation changes occurred. At 75 °C, lots of changes were evidenced: the Ca²⁺ depletion triggered the complete loss of h2b, h3c helices and S1, S2 and S3 β-sheets, and partial losses of H1, H2 and H3 α-helices. The absence of Ca²⁺ in apo-α-La and its leaving from holo-α-La triggered electrostatic repulsion among Asp82, Asp84 and Asp87, leading to the formation of a hydrophobic cluster involving Phe9, Phe31, Ile1, Va42, Ile55, Phe80 and Leu81. These conformational changes were related to an interfacial tension decrease and to a foaming capacity increase, for both apo-α-La and holo-α-La. This study exemplifies how powerful MD is as a tool to provide a better understanding of the molecular origins of food proteins' techno-functionalities.

Mechanistic insights into the urea-induced denaturation of a non-seleno thiol specific antioxidant human peroxiredoxin 6

Peroxiredoxin 6 (Prdx6) is a unique enzyme among mammalian peroxiredoxins as it lacks resolving cysteine. It is found to be involved in number of different diseases including tumours and its expression level is highest in lungs as compared to other organs. It has been found that Prdx6 plays a significant role different metabolic diseases, ocular damage, neurodegeneration and male infertility. It is a bifunctional protein having phospholipase A₂ and peroxidase (also has the ability to reduce phospholipid hydroperoxides) activities. In order to complete the peroxidise reaction cycle it requires glutathione catalyzed by glutathione S-transferase. Equilibrium unfolding and conformational stability of Prdx6 was studied by using urea as a chemical denaturant to understand the changes it goes under cellular stress conditions. Three different spectroscopic methods were employed to monitor urea-induced denaturation. From the results obtained, it was found that the urea denaturation of Prdx6 follows a variable two state process due to non-coincidence of the normalized transition curves obtained from different optical probes. The different denaturation curves were normalized and thermodynamic parameters, ΔG_D^o, Gibbs free energy change related to the urea-induced denaturation, midpoint of denaturation (C_m), and m = (δΔG_D / [urea]) were obtained. The structural information of Prdx6 were further analysed by several parameters obtained by 100 ns MD simulation. The results of MD simulation clearly favour the outcome of spectroscopic studies.

Spectroscopic and MD simulation studies on unfolding processes of mitochondrial carbonic anhydrase VA induced by urea

Carbonic anhydrase VA (CAVA) is primarily expressed in the mitochondria and involved in numerous physiological processes including lipogenesis, insulin secretion from pancreatic cells, ureagenesis, gluconeogenesis and neuronal transmission. To understand the biophysical properties of CAVA, we carried out a reversible urea-induced isothermal denaturation at pH 7.0 and 25°C. Spectroscopic probes, [θ]222 (mean residue ellipticity at 222 nm), F344 (Trp-fluorescence emission intensity at 344 nm) and Δε280 (difference absorption at 280 nm) were used to monitor the effect of urea on the structure and stability of CAVA. The urea-induced reversible denaturation curves were used to estimate [Formula: see text], Gibbs free energy in the absence of urea; Cm, the mid-point of the denaturation curve, i.e. molar urea concentration ([urea]) at which ΔGD = 0; and m, the slope (=∂ΔGD/∂[urea]). Coincidence of normalized transition curves of all optical properties suggests that unfolding/refolding of CAVA is a two-state process. We further performed 40 ns molecular dynamics simulation of CAVA to see the dynamics at different urea concentrations. An excellent agreement was observed between in silico and in vitro studies.

Quantifying the co-solvent effects on trypsin from the digestive system of carp Catla catla by biophysical techniques and molecular dynamics simulations

Urea molecules locate within 0.5 nm of the surface of trypsin.

Mg²⁺ coordinating dynamics in Mg:ATP fueled motor proteins

The coordination of Mg(2+) with the triphosphate group of adenosine triphosphate (ATP) in motor proteins is investigated using data mining and molecular dynamics. The possible coordination structures available from crystal data for actin, myosin, RNA polymerase, DNA polymerase, DNA helicase, and F1-ATPase are verified and investigated further by molecular dynamics. Coordination states are evaluated using structural analysis and quantified by radial distribution functions, coordination numbers, and pair interaction energy calculations. The results reveal a diverse range of both transitory and stable coordination arrangements between Mg(2+) and ATP. The two most stable coordinating states occur when Mg(2+) coordinates two or three oxygens from the triphosphate group of ATP. Evidence for five-site coordination is also reported involving water in addition to the triphosphate group. The stable states correspond to a pair interaction energy of either ∼-2750 kJ/mol or -3500 kJ/mol. The role of water molecules in the hydration shell surrounding Mg(2+) is also reported.

Cosolvent effects on protein stability

Proteins are marginally stable, and the folding/unfolding equilibrium of proteins in aqueous solution can easily be altered by the addition of small organic molecules known as cosolvents. Cosolvents that shift the equilibrium toward the unfolded ensemble are termed denaturants, whereas those that favor the folded ensemble are known as protecting osmolytes. Urea is a widely used denaturant in protein folding studies, and the molecular mechanism of its action has been vigorously debated in the literature. Here we review recent experimental as well as computational studies that show an emerging consensus in this problem. Urea has been shown to denature proteins through a direct mechanism, by interacting favorably with the peptide backbone as well as the amino acid side chains. In contrast, the molecular mechanism by which the naturally occurring protecting osmolyte trimethylamine N-oxide (TMAO) stabilizes proteins is not clear. Recent studies have established the strong interaction of TMAO with water. Detailed molecular simulations, when used with force fields that incorporate these interactions, can provide insight into this problem. We present the development of a model for TMAO that is consistent with experimental observations and that provides physical insight into the role of cosolvent-cosolvent interaction in determining its preferential interaction with proteins.

Effects of urea induced protein conformational changes on ion exchange chromatographic behavior

Urea is widely employed to facilitate protein separations in ion exchange chromatography at various scales. In this work, five model proteins were used to examine the chromatographic effects of protein conformational changes induced by urea in ion exchange chromatography. Linear gradient experiments were carried out at various urea concentrations and the protein secondary and tertiary structures were evaluated by far UV CD and fluorescence measurements, respectively. The results indicated that chromatographic retention times were well correlated with structural changes and that they were more sensitive to tertiary structural change. Steric Mass Action (SMA) isotherm parameters were also examined and the results indicated that urea induced protein conformational changes could affect both the characteristic charge and equilibrium constants in these systems. Dynamic light scattering analysis of changes in protein size due to urea-induced unfolding indicated that the size of the protein was not correlated with SMA parameter changes. These results indicate that while urea-induced structural changes can have a marked effect on protein chromatographic behavior in IEX, this behavior can be quite complicated and protein specific. These differences in protein behavior may provide insight into how these partially unfolded proteins are interacting with the resin material.

Equilibrium study of protein denaturation by urea

Though urea is commonly used to denature proteins, the molecular mechanism of its denaturing ability is still a subject of considerable debate. Previous molecular dynamics simulation studies have sought to elucidate the mechanism of urea denaturation by focusing on the pathway of denaturation rather than examining the effect of urea on the folding/unfolding equilibrium, which is commonly measured in experiment. Here we report the reversible folding/unfolding equilibrium of Trp-cage miniprotein in the presence of urea, over a broad range of urea concentrations, using all-atom Replica exchange MD simulations. The simulations capture the experimentally observed linear dependence of unfolding free energy on urea concentration. We find that the denaturation is driven by favorable direct interaction of urea with the protein through both electrostatic and van der Waals forces and quantify their contribution. Though the magnitude of direct electrostatic interaction of urea is larger than van der Waals, the difference between unfolded and folded ensembles is dominated by the van der Waals interaction. We also find that hydrogen bonding of urea to the peptide backbone does not play a dominant role in denaturation. The unfolded ensemble sampled depends on urea concentration, with greater urea concentration favoring conformations with greater solvent exposure. The m-value is predicted to increase with temperature and more strongly so with pressure.

Urea-mediated protein denaturation: a consensus view

We have performed all-atom molecular dynamics simulations of three structurally similar small globular proteins in 8 M urea and compared the results with pure aqueous simulations. Protein denaturation is preceded by an initial loss of water from the first solvation shell and consequent in-flow of urea toward the protein. Urea reaches the first solvation shell of the protein mainly due to electrostatic interaction with a considerable contribution coming from the dispersion interaction. Urea shifts the equilibrium from the native to denatured ensemble by making the protein-protein contact less stable than protein-urea contact, which is just the reverse of the condition in pure water, where protein-protein contact is more stable than protein-water contact. We have also seen that water follows urea and reaches the protein interior at later stages of denaturation, while urea preferentially and efficiently solvates different parts of the protein. Solvation of the protein backbone via hydrogen bonding, favorable electrostatic interaction with hydrophilic residues, and dispersion interaction with hydrophobic residues are the key steps through which urea intrudes the core of the protein and denatures it. Why urea is preferred over water for binding to the protein backbone and how urea orients itself toward the protein backbone have been identified comprehensively. All the key components of intermolecular forces are found to play a significant part in urea-induced protein denaturation and also toward the stability of the denatured state ensemble. Changes in water network/structure and dynamical properties and higher degree of solvation of the hydrophobic residues validate the presence of "indirect mechanism" along with the "direct mechanism" and reinforce the effect of urea on protein.

Probing the urea dependence of residual structure in denatured human alpha-lactalbumin

Backbone (15)N relaxation parameters and (15)N-(1)H(N) residual dipolar couplings (RDCs) have been measured for a variant of human alpha-lactalbumin (alpha-LA) in 4, 6, 8 and 10 M urea. In the alpha-LA variant, the eight cysteine residues in the protein have been replaced by alanines (all-Ala alpha-LA). This protein is a partially folded molten globule at pH 2 and has been shown previously to unfold in a stepwise non-cooperative manner on the addition of urea. (15)N R(2) values in some regions of all-Ala alpha-LA show significant exchange broadening which is reduced as the urea concentration is increased. Experimental RDC data are compared with RDCs predicted from a statistical coil model and with bulkiness, average area buried upon folding and hydrophobicity profiles in order to identify regions of non-random structure. Residues in the regions corresponding to the B, D and C-terminal 3(10) helices in native alpha-LA show R(2) values and RDC data consistent with some non-random structural propensities even at high urea concentrations. Indeed, for residues 101-106 the residual structure persists in 10 M urea and the RDC data suggest that this might include the formation of a turn-like structure. The data presented here allow a detailed characterization of the non-cooperative unfolding of all-Ala alpha-LA at higher concentrations of denaturant and complement previous studies which focused on structural features of the molten globule which is populated at lower concentrations of denaturant.

Stability of the core domain of p53: insights from computer simulations

Background

The tumour suppressor protein p53 protein has a core domain that binds DNA and is the site for most oncogenic mutations. This domain is quite unstable compared to its homologs p63 and p73. Two key residues in the core domain of p53 (Tyr236, Thr253), have been mutated in-silico, to their equivalent residues in p63 (Phe238 and Ile255) and p73 (Phe238 and Ile255), with subsequent increase in stability of p53. Computational studies have been performed to examine the basis of instability in p53.

Results

Molecular dynamics simulations suggest that mutations in p53 lead to increased conformational sampling of the phase space which stabilizes the system entropically. In contrast, reverse mutations, where p63 and p73 were mutated by replacing the Phe238 and Ile255 by Tyr and Thr respectively (as in p53), showed reduced conformational sampling although the change for p63 was much smaller than that for p73. Barriers to the rotation of sidechains containing aromatic rings at the core of the proteins were reduced several-fold when p53 was mutated; in contrast they increased when p73 was mutated and decreased by a small amount in p63. The rate of ring flipping of a Tyrosine residue at the boundary of two domains can be correlated with the change in stability, with implications for possible pathways of entry of agents that induce unfolding.

Conclusion

A double mutation at the core of the DNA binding domain of p53 leads to enhanced stability by increasing the softness of the protein. A change from a highly directional polar interaction of the core residues Tyr236 and Thr253 to a non-directional apolar interaction between Phe and Ile respectively may enable the system to adapt more easily and thus increase its robustness to structural perturbations, giving it increased stability. This leads to enhanced conformational sampling which in turn is associated with an increased "softness" of the protein core. However the system seems to become more rigid at the periphery. The success of this methodology in reproducing the experimental trends in the stability of p53 suggests that it has the potential to complement structural studies for rapidly estimating changes in stability upon mutations and could be an additional tool in the design of specific classes of proteins.

Urea and guanidinium chloride denature protein L in different ways in molecular dynamics simulations

In performing protein-denaturation experiments, it is common to employ different kinds of denaturants interchangeably. We make use of molecular dynamics simulations of Protein L in water, in urea, and in guanidinium chloride (GdmCl) to ascertain if there are any structural differences in the associated unfolding processes. The simulation of proteins in solutions of GdmCl is complicated by the large number of charges involved, making it difficult to set up a realistic force field. Furthermore, at high concentrations of this denaturant, the motion of the solvent slows considerably. The simulations show that the unfolding mechanism depends on the denaturing agent: in urea the beta-sheet is destabilized first, whereas in GdmCl, it is the alpha-helix. Moreover, whereas urea interacts with the protein accumulating in the first solvation shell, GdmCl displays a longer-range electrostatic effect that does not perturb the structure of the solvent close to the protein.

Characterization of the protein unfolding processes induced by urea and temperature

Correct folding is critical for the biological activities of proteins. As a contribution to a better understanding of the protein (un)folding problem, we studied the effect of temperature and of urea on peptostreptococcal Protein L destructuration. We performed standard molecular dynamics simulations at 300 K, 350 K, 400 K, and 480 K, both in 10 M urea and in water. Protein L followed at least two alternative unfolding pathways. Urea caused the loss of secondary structure acting preferentially on the beta-sheets, while leaving the alpha-helices almost intact; on the contrary, high temperature preserved the beta-sheets and led to a complete loss of the alpha-helices. These data suggest that urea and high temperature act through different unfolding mechanisms, and protein secondary motives reveal a differential sensitivity to various denaturant treatments. As further validation of our results, replica-exchange molecular dynamics simulations of the temperature-induced unfolding process in the presence of urea were performed. This set of simulations allowed us to compute the thermodynamical parameters of the process and confirmed that, in the configurational space of Protein L unfolding, both of the above pathways are accessible, although to a different relative extent.

Urea-amide preferential interactions in water: quantitative comparison of model compound data with biopolymer results using water accessible surface areas

Two fundamentally different thermodynamic approaches are in use to interpret or predict the effects of urea on biopolymer processes: one is a synthesis of transfer free energies obtained from measurements of the effects of urea on the solubilities of small, model compounds; the other is an analysis of preferential interactions of urea with a range of folded and unfolded biopolymer surfaces. Here, we compare the predictions of these two approaches for the contribution of urea-amide (peptide) interactions to destabilization of folded proteins by urea. For these comparisons, we develop independent thermodynamic analyses of osmometric and solubility data characterizing interactions of a model compound with urea (or any other solute) and apply them to all five model compounds (glycine, alanine, diglycine, glycylalanine, and triglycine) where both isopiestic distillation (ID) and solubility data in aqueous urea solutions are available. We use model-independent expressions to calculate mu ex 23, the derivative of the "excess" chemical potential of solute "2" (either a model compound or a biopolymer) with respect to the molality of solute "3" (urea). Analyses of ID data for these systems reveal significant dependences of mu ex 23 on both m2 and m3, which must be taken into account in making comparisons with values of mu ex 23 obtained from solubility studies or from analyses of urea-biopolymer preferential interactions. Values of mu ex 23 calculated from model compound ID data at low m2 and m3 are directly proportional to the amount of polar amide (N, O) surface area, and not to any other type of surface. The proportionality constant in this limit, mu ex 23 /(RT x ASA) = (1.0 +/- 0.1) x 10(-3) A(-2), is very similar to that previously obtained by analysis of urea-biopolymer preferential interactions ((1.4 +/- 0.3) x 10(-3) A(-2)). This level of agreement for amide surface in the low concentration limit, as well as the absence of any significant preferential interaction of urea with Gly and Ala, reinforces the conclusion that the primary preferential interaction of urea with protein surface is a favorable interaction (resulting in local accumulation of urea) at polar amide surface, located mostly on the peptide backbone. However, mu ex 23 for interactions of urea with these model amides is found from both ID and solubility data to be urea concentration-dependent, in contrast to the urea concentration independence of the analogous quantity for protein unfolding.

Biomolecular modeling: Goals, problems, perspectives

Computation based on molecular models is playing an increasingly important role in biology, biological chemistry, and biophysics. Since only a very limited number of properties of biomolecular systems is actually accessible to measurement by experimental means, computer simulation can complement experiment by providing not only averages, but also distributions and time series of any definable quantity, for example, conformational distributions or interactions between parts of systems. Present day biomolecular modeling is limited in its application by four main problems: 1) the force-field problem, 2) the search (sampling) problem, 3) the ensemble (sampling) problem, and 4) the experimental problem. These four problems are discussed and illustrated by practical examples. Perspectives are also outlined for pushing forward the limitations of biomolecular modeling.

Aqueous urea solutions: structure, energetics, and urea aggregation

Urea is ubiquitously used as a protein denaturant. To study the structure and energetics of aqueous urea solutions, we have carried out molecular dynamics simulations for a wide range of urea concentrations and temperatures. The hydrogen bonds between urea and water were found to be significantly weaker than those between water molecules, which drives urea self-aggregation due to the hydrophobic effect. From the reduction of the water exposed urea surface area, urea was found to exhibit an aggregation degree of ca. 20% at concentrations commonly used for protein denaturation. Structurally, three distinct urea pair conformations were identified and their populations were analyzed by translational and orientational pair distribution functions. Furthermore, urea was found to strengthen water structure in terms of hydrogen bond energies and population of solvation shells. Our findings are consistent with a direct interaction between urea and the protein as the main driving force for protein denaturation. As an additional, more indirect effect, urea was found to enhance water structure, which would suggest a weakening of the hydrophobic effect.

Electrostatic contributions to protein stability and folding energy

The ability to predict the thermal stability of proteins based on their corresponding sequence is a problem of great fundamental and practical importance. Here we report an approach for calculating the electrostatic contribution to protein stability based on the use of the semimacroscopic protein dipole Langevin dipole (PDLD/S) in its linear response approximation version for self-energy with a dielectric constant, (epsilon(p)) and an effective dielectric for charge-charge interactions (epsilon(eff)). The method is applied to the test cases of ubiquitin, lipase, dihydrofolate reductase and cold shock proteins with series of epsilon(p) and epsilon(eff). It is found that the optimal values of these dielectric constants lead to very promising results, both for the relative stability and the absolute folding energy. Consideration of the specific values of the optimal dielectric constants leads to an exciting conceptual description of the reorganization effect during the folding process. Although this description should be examined by further microscopic studies, the practical use of the current approach seems to offer a powerful tool for protein design and for studies of the energetics of protein folding.

Simulations of macromolecules in protective and denaturing osmolytes: properties of mixed solvent systems and their effects on water and protein structure and dynamics

Rarely is any solution simply solute and water. In vivo, solutes, such as proteins and nucleic acids, swim in a sea of water, salts, ions, small molecules, and lipids, not to mention other macromolecules. In vitro, virtually all solutions contain a mixture of aqueous solvents, or "cosolvents" [i.e., solvent(s) in addition to water], that can alter the dynamics, behavior, solubility, and stability of proteins and nucleic acids. We have developed models for a number of cosolvents, including the denaturant urea and the small chemical chaperone trimethylamine N-oxide (TMAO). This chapter examines the models for these two cosolvents in the context of experimental data. The direct and indirect effects of these molecules on water and protein are studied with molecular dynamics simulations. These observations and conclusions are drawn from simulations of these molecules in pure water and as a cosolvent for the protein chymotrypsin inhibitor 2. Urea-induced denaturation occurs initially through attack of the protein by water and hydration of hydrophobic protein moieties as a result of disruption of the hydrogen bonding network of water by urea. This indirect denaturing effect of urea is followed by more direct action as urea replaces some waters involved in the initial hydration of the hydrophobic core and subsequently binds to polar residues and the protein main chain to compete with the intraprotein hydrogen bonds. In the case of TMAO, we find that it encourages water-water interactions, thereby stabilizing the protein as a result of the increased penalty for the hydration of hydrophobic residues.

Modeling electrostatic effects in proteins

Electrostatic energies provide what is perhaps the most effective tool for structure-function correlation of biological molecules. This review considers the current state of simulations of electrostatic energies in macromolecules as well as the early developments of this field. We focus on the relationship between microscopic and macroscopic models, considering the convergence problems of the microscopic models and the fact that the dielectric 'constants' in semimacroscopic models depend on the definition and the specific treatment. The advances and the challenges in the field are illustrated considering a wide range of functional properties including pK(a)'s, redox potentials, ion and proton channels, enzyme catalysis, ligand binding and protein stability. We conclude by pointing out that, despite the current problems and the significant misunderstandings in the field, there is an overall progress that should lead eventually to quantitative descriptions of electrostatic effects in proteins and thus to quantitative descriptions of the function of proteins.

Effective interactions between chaotropic agents and proteins

Chaotropic agents are cosolutes that can disrupt the hydrogen bonding network between water molecules and reduce the stability of the native state of proteins by weakening the hydrophobic effect. In this work, we represent the chaotropic agent as a factor that reduces the amount of order in the structures formed by water molecules, both in the bulk and the hydration shells around hydrophobic amino acids. In this framework we show that low chaotrope concentrations lead to a destabilization of the native state of proteins, and that high concentrations induce complete denaturation. We also find that the reduction of the number of bulk ordered states of water molecules can give origin to an effective interaction between chaotropic molecules and proteins.

The GROMOS software for biomolecular simulation: GROMOS05

We present the latest version of the Groningen Molecular Simulation program package, GROMOS05. It has been developed for the dynamical modelling of (bio)molecules using the methods of molecular dynamics, stochastic dynamics, and energy minimization. An overview of GROMOS05 is given, highlighting features not present in the last major release, GROMOS96. The organization of the program package is outlined and the included analysis package GROMOS++ is described. Finally, some applications illustrating the various available functionalities are presented.