Equilibrium study of protein denaturation by urea

Free-energy analysis of protein solvation with all-atom molecular dynamics simulation combined with a theory of solutions

The structure of a protein is strongly influenced by solvation. In the present review, the solvation effect is treated within the framework of statistical thermodynamics. The key quantity is the solvation free energy of protein and a fast method for its computation with explicit solvent is introduced. The applications of the method to the structural fluctuation of protein in water, structure determination of protein-protein complex, and urea effect on protein structure are then described, and the roles of solvent water and cosolvent are discussed from the standpoint of energetics.

Lectin from seeds of a Brazilian lima bean variety (Phaseolus lunatus L. var. cascavel) presents antioxidant, antitumour and gastroprotective activities

Lectins are proteins able to interact specifically and reversibly with carbohydrates. They are present in all living beings, particularly in legume seeds, which have many biological functions. The aim of this study was to isolate, characterize and verify antioxidant, anti-hemolytic, antitumor and gastroprotective activities in a lectin present in seeds of Phaseolus lunatus L. var. cascavel (PLUN). The isolation of lectin was performed by size exclusion chromatography on Sephadex G-100, which was isolated from a protein capable of agglutinating only human erythrocytes type A, being this the only inhibited haemagglutination n-acetyl-d-galactosamine. Its weight was estimated by PAGE is 128kDa. The lectin is thermostable up to 80°C and is active between pH 2-11. As 8M urea was able to denature the lectin. PLUN is a glycoprotein consisting of 2% carbohydrate and has antioxidant action with ascorbic acid equivalent antioxidant capacity (μMAA/g) of 418.20, 326 and 82.9 for total antioxidant activity, ABTS radical capture and capture of DPPH radical, respectively. The lectin has antitumor activity against melanoma derived cells at doses of 100 and 50mg/ml, reducing up to 83% tumor cells, and gastroprotective action, reducing up to 63% damaged area of ​​the stomach induced by ethanol.

Molecular basis of the osmolyte effect on protein stability: a lesson from the mechanical unfolding of lysozyme

Osmolytes are a class of small organic molecules that shift the protein folding equilibrium. For this reason, they are accumulated by organisms under environmental stress and find applications in biotechnology where proteins need to be stabilized or dissolved. However, despite years of research, debate continues over the exact mechanisms underpinning the stabilizing and denaturing effect of osmolytes. Here, we simulated the mechanical denaturation of lysozyme in different solvent conditions to study the molecular mechanism by which two biologically relevant osmolytes, denaturing (urea) and stabilizing (betaine), affect the folding equilibrium. We found that urea interacts favorably with all types of residues via both hydrogen bonds and dispersion forces, and therefore accumulates in a diffuse solvation shell around the protein. This not only provides an enthalpic stabilization of the unfolded state, but also weakens the hydrophobic effect, as hydrophobic forces promote the association of urea with nonpolar residues, facilitating the unfolding. In contrast, we observed that betaine is excluded from the protein backbone and nonpolar side chains, but is accumulated near the basic residues, yielding a nonuniform distribution of betaine molecules at the protein surface. Spatially resolved solvent–protein interaction energies further suggested that betaine behaves in a ligand- rather than solvent-like manner and its exclusion from the protein surface arises mostly from the scarcity of favorable binding sites. Finally, we found that, in the presence of betaine, the reduced ability of water molecules to solvate the protein results in an additional enthalpic contribution to the betaine-induced stabilization.

Probing the Action of Chemical Denaturant on an Intrinsically Disordered Protein by Simulation and Experiment

Chemical denaturants are the most commonly used agents for unfolding proteins and are thought to act by better solvating the unfolded state. Improved solvation is expected to lead to an expansion of unfolded chains with increasing denaturant concentration, providing a sensitive probe of the denaturant action. However, experiments have so far yielded qualitatively different results concerning the effects of chemical denaturation. Studies using Förster resonance energy transfer (FRET) and other methods found an increase in radius of gyration with denaturant concentration, but most small-angle X-ray scattering (SAXS) studies found no change. This discrepancy therefore challenges our understanding of denaturation mechanism and more generally the accuracy of these experiments as applied to unfolded or disordered proteins. Here, we use all-atom molecular simulations to investigate the effect of urea and guanidinium chloride on the structure of the intrinsically disordered protein ACTR, which can be studied by experiment over a wide range of denaturant concentration. Using unbiased molecular simulations with a carefully calibrated denaturant model, we find that the protein chain indeed swells with increasing denaturant concentration. This is due to the favorable association of urea or guanidinium chloride with the backbone of all residues and with the side-chains of almost all residues, with denaturant-water transfer free energies inferred from this association in reasonable accord with experimental estimates. Interactions of the denaturants with the backbone are dominated by hydrogen bonding, while interactions with side-chains include other contributions. By computing FRET efficiencies and SAXS intensities at each denaturant concentration, we show that the simulation trajectories are in accord with both experiments on this protein, demonstrating that there is no fundamental inconsistency between the two types of experiment. Agreement with experiment also supports the picture of chemical denaturation described in our simulations, driven by weak association of denaturant with the protein. Our simulations support some assumptions needed for each experiment to accurately reflect changes in protein size, namely, that the commonly used FRET chromophores do not qualitatively alter the results and that possible effects such as preferential solvent partitioning into the interior of the chain do not interfere with the determination of radius of gyration from the SAXS experiments.

Hydrophobic Association in Mixed Urea-TMAO Solutions

The formation of a hydrophobic core is key to the folding and resulting function of most proteins in the cell. In several organisms, as well as in many in vitro experiments, protein folding is modulated by the presence of osmolytes, but the mechanism by which hydrophobic association occurs is not well understood. We present a study of the solvation thermodynamics of hydrophobic self-association in mixed-osmolyte urea-TMAO solutions, with neopentane as a model hydrophobic molecule. Using molecular dynamics simulations and the Kirkwood-Buff theory of solutions, we show that a sensitive balance between the TMAO-water and the TMAO-urea interactions governs the osmolyte-induced changes in hydrophobic association in mixed urea-TMAO solutions. This balance must be correctly incorporated in force-field parametrization because hydrophobic association can be either enhanced or prevented all together by slightly increasing or decreasing the osmolyte-water affinity and osmolyte-osmolyte self-affinity of TMAO molecules.

The aggregation of cytochrome C may be linked to its flexibility during refolding

Large-scale expression of biopharmaceutical proteins in cellular hosts results in production of large insoluble mass aggregates. In order to generate functional product, these aggregates require further processing through refolding with denaturant, a process in itself that can result in aggregation. Using a model folding protein, cytochrome C, we show how an increase in final denaturant concentration decreases the propensity of the protein to aggregate during refolding. Using polarised fluorescence anisotropy, we show how reduced levels of aggregation can be achieved by increasing the period of time the protein remains flexible during refolding, mediated through dilution ratios. This highlights the relationship between the flexibility of a protein and its propensity to aggregate. We attribute this behaviour to the preferential urea-residue interaction, over self-association between molecules.

Using a second-order differential model to fit data without baselines in protein isothermal chemical denaturation

In vitro protein stability studies are commonly conducted via thermal or chemical denaturation/renaturation of protein. Conventional data analyses on the protein unfolding/(re)folding require well-defined pre- and post-transition baselines to evaluate Gibbs free-energy change associated with the protein unfolding/(re)folding. This evaluation becomes problematic when there is insufficient data for determining the pre- or post-transition baselines. In this study, fitting on such partial data obtained in protein chemical denaturation is established by introducing second-order differential (SOD) analysis to overcome the limitations that the conventional fitting method has. By reducing numbers of the baseline-related fitting parameters, the SOD analysis can successfully fit incomplete chemical denaturation data sets with high agreement to the conventional evaluation on the equivalent completed data, where the conventional fitting fails in analyzing them. This SOD fitting for the abbreviated isothermal chemical denaturation further fulfills data analysis methods on the insufficient data sets conducted in the two prevalent protein stability studies.

Microscopic significance of hydrophobic residues in the protein-stabilizing effect of trimethylamine N-oxide (TMAO)

Proteins with a higher hydrophobic content are better protected by TMAO against the deleterious effect of urea.

Empirical Optimization of Interactions between Proteins and Chemical Denaturants in Molecular Simulations

Chemical denaturants are the most commonly used perturbation applied to study protein stability and folding kinetics as well as the properties of unfolded polypeptides. We build on recent work balancing the interactions of proteins and water, and accurate models for the solution properties of urea and guanidinium chloride, to develop a combined force field that is able to capture the strength of interactions between proteins and denaturants. We use solubility data for a model tetraglycine peptide in each denaturant to tune the protein-denaturant interaction by a novel simulation methodology. We validate the results against data for more complex sequences: single-molecule Förster resonance energy transfer data for a 34-residue fragment of the globular protein CspTm and photoinduced electron transfer quenching data for the disordered peptides C(AGQ)nW in denaturant solution as well as the chemical denaturation of the mini-protein Trp cage. The combined force field model should aid our understanding of denaturation mechanisms and the interpretation of experiment.

Correlation analysis for heat denaturation of Trp-cage miniprotein with explicit solvent

Energetics was analyzed for Trp-cage miniprotein in water to elucidate the solvation effect in heat denaturation. The solvation free energy was computed for a set of protein structures at room and high temperatures with all-atom molecular dynamics simulation combined with the solution theory in the energy representation, and its correlations were investigated against the intramolecular (structural) energy of the protein and the average interaction energy of the protein with the solvent water. It was observed both at room and high temperatures that the solvation free energy is anticorrelated to the structural energy and varies in parallel to the electrostatic component of the protein-water interaction energy without correlations to the van der Waals and excluded-volume components. When the set of folded structures sampled at room temperature was compared with the set of unfolded ones at high temperature, it was found that the preference order of the two sets is in correspondence to the van der Waals and excluded-volume components in the sum of the protein intramolecular and protein-water intermolecular interactions and is not distinguished by the electrostatic component.

Heat, Acid and Chemically Induced Unfolding Pathways, Conformational Stability and Structure-Function Relationship in Wheat α-Amylase

Wheat α-amylase, a multi-domain protein with immense industrial applications, belongs to α+β class of proteins with native molecular mass of 32 kDa. In the present study, the pathways leading to denaturation and the relevant unfolded states of this multi-domain, robust enzyme from wheat were discerned under the influence of temperature, pH and chemical denaturants. The structural and functional aspects along with thermodynamic parameters for α-amylase unfolding were probed and analyzed using fluorescence, circular dichroism and enzyme assay methods. The enzyme exhibited remarkable stability up to 70°C with tendency to aggregate at higher temperature. Acid induced unfolding was also incomplete with respect to the structural content of the enzyme. Strong ANS binding at pH 2.0 suggested the existence of a partially unfolded intermediate state. The enzyme was structurally and functionally stable in the pH range 4.0–9.0 with 88% recovery of hydrolytic activity. Careful examination of biophysical properties of intermediate states populated in urea and GdHCl induced denaturation suggests that α-amylase unfolding undergoes irreversible and non-coincidental cooperative transitions, as opposed to previous reports of two-state unfolding. Our investigation highlights several structural features of the enzyme in relation to its catalytic activity. Since, α-amylase has been comprehensively exploited for use in a range of starch-based industries, in addition to its physiological significance in plants and animals, knowledge regarding its stability and folding aspects will promote its biotechnological applications.

A hypothesis to reconcile the physical and chemical unfolding of proteins

High pressure (HP) or urea is commonly used to disturb folding species. Pressure favors the reversible unfolding of proteins by causing changes in the volumetric properties of the protein-solvent system. However, no mechanistic model has fully elucidated the effects of urea on structure unfolding, even though protein-urea interactions are considered to be crucial. Here, we provide NMR spectroscopy and 3D reconstructions from X-ray scattering to develop the "push-and-pull" hypothesis, which helps to explain the initial mechanism of chemical unfolding in light of the physical events triggered by HP. In studying MpNep2 from Moniliophthora perniciosa, we tracked two cooperative units using HP-NMR as MpNep2 moved uphill in the energy landscape; this process contrasts with the overall structural unfolding that occurs upon reaching a threshold concentration of urea. At subdenaturing concentrations of urea, we were able to trap a state in which urea is preferentially bound to the protein (as determined by NMR intensities and chemical shifts); this state is still folded and not additionally exposed to solvent [fluorescence and small-angle X-ray scattering (SAXS)]. This state has a higher susceptibility to pressure denaturation (lower p1/2 and larger ΔVu); thus, urea and HP share concomitant effects of urea binding and pulling and water-inducing pushing, respectively. These observations explain the differences between the molecular mechanisms that control the physical and chemical unfolding of proteins, thus opening up new possibilities for the study of protein folding and providing an interpretation of the nature of cooperativity in the folding and unfolding processes.

Engineering tandem single-chain Fv as cell surface reporters with enhanced properties of fluorescence detection

A recently described fluorescence biosensor platform utilizes single-chain Fv (scFvs) that selectively bind and activate fluorogen molecules. In this report we investigated the display of tandem scFv biosensors at the surface of mammalian cells with the aim of advancing current fluorescence detection strategies. We initially screened different peptide linkers to separate each scFv unit, and discovered that tandem proteins joined by either flexible or α-helical linkers properly fold and display at the surface of mammalian cells. Accordingly, we performed a combinatorial scFv-dimer study and identified that fluorescence activation correlated with the cellular location (membrane distal versus proximal) and selections of the different scFvs. Furthermore, in vitro measurements showed that the stability of each scFv monomer unit influenced the folding and cell surface activities of tandem scFvs. Additionally, we investigated the absence or poor signals from some scFv-dimer combinations and discovered that intramolecular and intermolecular scFv chain mispairings led to protein misfolding and/or secretory-pathway-mediated degradation. Furthermore, when tandem scFvs were utilized as fluorescence reporter tags with surface receptors, the biosensor unit and target protein showed independent activities. Thus, the live cell application of tandem scFvs permitted advanced detection of target proteins via fluorescence signal amplification, Förster resonance energy transfer resulting in the increase of Stokes shift and multi-color vesicular traffic of surface receptors.

Dispersion interactions between urea and nucleobases contribute to the destabilization of RNA by urea in aqueous solution

Urea has long been used to investigate protein folding and, more recently, RNA folding. Studies have proposed that urea denatures RNA by participating in stacking interactions and hydrogen bonds with nucleic acid bases. In this study, the ability of urea to form unconventional stacking interactions with RNA bases is investigated using ab initio calculations (RI-MP2 and CCSD(T) methods with the aug-cc-pVDZ basis set). A total of 29 stable nucleobase-urea stacked complexes are identified in which the intermolecular interaction energies (up to -14 kcal/mol) are dominated by dispersion effects. Natural bond orbital (NBO) and atoms in molecules (AIM) calculations further confirm strong interactions between urea and nucleobases. Calculations on model systems with multiple urea and water molecules interacting with a guanine base lead to a hypothesis that urea molecules along with water are able to form cage-like structures capable of trapping nucleic acid bases in extrahelical states by forming both hydrogen-bonded and dispersion interactions, thereby contributing to the unfolding of RNA in the presence of urea in aqueous solution.

Mutual Exclusion of Urea and Trimethylamine N-Oxide from Amino Acids in Mixed Solvent Environment

We study the solvation of amino acids in pure-osmolyte and mixed-osmolyte urea and trimethylamine N-oxide (TMAO) solutions using molecular dynamics simulations. Analysis of Kirkwood-Buff integrals between the solution components provides evidence that in the mixed osmolytic solution, both urea and TMAO are mutually excluded from the amino acid surface, accompanied by an increase in osmolyte-osmolyte aggregation. Similar observations are made in simulations of a model protein backbone, represented by triglycine, and suggest that TMAO stabilizes proteins under urea denaturation conditions by effectively removing urea from the protein surface. The effects of the mixed osmolytes on the solvation of the amino acids and the backbone are found to be highly nonlinear in terms of the effects of the individual osmolytes and independent of differences in the strength of the TMAO-water interactions, as observed with different TMAO force fields.

Protein denaturants at aqueous-hydrophobic interfaces: self-consistent correlation between induced interfacial fluctuations and denaturant stability at the interface

The notion of direct interaction between denaturing cosolvent and protein residues has been proposed in dialogue relevant to molecular mechanisms of protein denaturation. Here we consider the correlation between free energetic stability and induced fluctuations of an aqueous-hydrophobic interface between a model hydrophobically associating protein, HFBII, and two common protein denaturants, guanidinium cation (Gdm(+)) and urea. We compute potentials of mean force along an order parameter that brings the solute molecule close to the known hydrophobic region of the protein. We assess potentials of mean force for different relative orientations between the protein and denaturant molecule. We find that in both cases of guanidinium cation and urea relative orientations of the denaturant molecule that are parallel to the local protein-water interface exhibit greater stability compared to edge-on or perpendicular orientations. This behavior has been observed for guanidinium/methylguanidinium cations at the liquid-vapor interface of water, and thus the present results further corroborate earlier findings. Further analysis of the induced fluctuations of the aqueous-hydrophobic interface upon approach of the denaturant molecule indicates that the parallel orientation, displaying a greater stability at the interface, also induces larger fluctuations of the interface compared to the perpendicular orientations. The correlation of interfacial stability and induced interface fluctuation is a recurring theme for interface-stable solutes at hydrophobic interfaces. Moreover, observed correlations between interface stability and induced fluctuations recapitulate connections to local hydration structure and patterns around solutes as evidenced by experiment (Cooper et al., J. Phys. Chem. A 2014, 118, 5657.) and high-level ab initio/DFT calculations (Baer et al., Faraday Discuss 2013, 160, 89).

On the urea induced hydrophobic collapse of a water soluble polymer

Polymer collapse despite cosolvent binding: solvation of extended coil conformations is entropically penalized, therefore stabilizing compact globular conformations in the coil-globule equilibrium of poly(N-isopropylacrylamide) in aqueous urea solution.

Infinitely dilute partial molar properties of proteins from computer simulation

A detailed understanding of temperature and pressure effects on an infinitely dilute protein's conformational equilibrium requires knowledge of the corresponding infinitely dilute partial molar properties. Established molecular dynamics methodologies generally have not provided a way to calculate these properties without either a loss of thermodynamic rigor, the introduction of nonunique parameters, or a loss of information about which solute conformations specifically contributed to the output values. Here we implement a simple method that is thermodynamically rigorous and possesses none of the above disadvantages, and we report on the method's feasibility and computational demands. We calculate infinitely dilute partial molar properties for two proteins and attempt to distinguish the thermodynamic differences between a native and a denatured conformation of a designed miniprotein. We conclude that simple ensemble average properties can be calculated with very reasonable amounts of computational power. In contrast, properties corresponding to fluctuating quantities are computationally demanding to calculate precisely, although they can be obtained more easily by following the temperature and/or pressure dependence of the corresponding ensemble averages.

The role of cross-chain ionic interactions for the stability of collagen model peptides

The contribution of ionic interactions to the stability of the collagen triple helix was studied using molecular dynamics (MD) simulations and biophysical methods. To this end, we examined the stability of a host-guest collagen model peptide, Ac-GPOGPOGPYGXOGPOGPO-NH2, substituting KGE, KGD, EGK, and DGK for the YGX sequence. All-atom, implicit solvent MD simulations show that the fraction of cross-chain ionic interactions formed is different, with the most pronounced in the KGE and KGD sequences, and the least in the DGK sequence. To test whether the fraction of cross-chain ionic interactions correlates with the stability, experimental measurements of thermostability were done using differential scanning calorimetry and circular dichroism spectroscopy. It was found that the melting temperature is very similar for KGE and KGD peptides, whereas the EGK peptide has lower thermostability and the DGK peptide is the least thermostable. A novel, to our knowledge, computational protocol termed temperature-scan MD was applied to estimate the relative stabilities of the peptides from MD simulations. We found an excellent correlation between transition temperatures obtained from temperature-scan MD and those measured experimentally. These results suggest the importance of cross-chain ionic interactions for the stability of collagen triple helix and confirm the utility of MD simulations in predicting interactions and stability in this system.

Understanding the role of H-bonding in aqueous self-assembly of two naphthalene diimide (NDI)-conjugated amphiphiles

Supramolecular architectures with the synchronized combination of various directional noncovalent forces are ubiquitous in biological systems. However, reports of such abiotic synthetic systems involving H-bonding in aqueous medium are rare due to the challenge faced in the formation of such structures by overcoming the competition from the water molecules. In this paper we have studied self-assembly of two structurally related naphthalene-diimide (NDI) conjugated bola-amphiphiles (NDI-1 and NDI-2) in water with an aim to realize the specific role of H-bonding among the hydrazide units present in one of the two building blocks (NDI-2) on the self-assembly. Both chromophores showed vesicular assembly in aqueous solution driven primarily by π-stacking among the NDI chromophores, which could be probed by UV-vis absorption spectra. Contrary to common belief, the lack of an H-bonding group in NDI-1 was found to be a boon in disguise in terms of the stability of the aggregates. Whereas NDI-2 aggregates showed LCST around 65-70 °C owing to the breaking of the H-bonds with increased temperature, the NDI-1 aggregates were found to be structurally intact until 90 °C, which may be attributed to the increased hydrophobicity introduced by the absence of the polar hydrazide group. Further concentration- and solvent-dependent UV-vis studies showed that NDI-1 formed assembled structure at greatly dilute solution and also in a solvent such as THF, confirming greater propensity for its self-assembly. As both bola-amphiphiles contain an electron-deficient NDI chromophore, interaction of their vesicles was studied with an externally added electron-rich pyrene derivative. Surprisingly, NDI-1 did not show any charge-transfer interaction with the donor, whereas NDI-2 could effectively intercalate, leading to a functional membrane with tunable surface functionalities. This was attributed to the additional stability of the intercalated state by H-bonding among the hydrazide units.

Urea induced unfolding dynamics of flavin adenine dinucleotide (FAD): spectroscopic and molecular dynamics simulation studies from femto-second to nanosecond regime

Here, we investigate the effect of urea in the unfolding dynamics of flavin adenine dinucleotide (FAD), an important enzymatic cofactor, through steady state, time-resolved fluorescence spectroscopic and molecular dynamics (MD) simulation studies. Steady state results indicate the possibility of urea induced unfolding of FAD, inferred from increasing emission intensity of FAD with urea. The TCSPC and up-conversion results suggest that the stack-unstack dynamics of FAD severely gets affected in the presence of urea and leads to an increase in the unstack conformation population from 15% in pure water to 40% in 12 M urea. Molecular dynamics simulation was employed to understand the nature of the interaction between FAD and urea at the molecular level. Results depict that urea molecules replace many of the water molecules around adenine and isoalloxazine rings of FAD. However, the major driving force for the stability of this unstack conformations arises from the favorable stacking interaction of a significant fraction of the urea molecules with adenine and isoalloxazine rings of FAD, which overcomes the intramolecular stacking interaction between themselves observed in pure water.

Folding and unfolding thermodynamics of the TC10b Trp-cage miniprotein

We examine the folding-unfolding of a variant of the Trp-cage, known as TC10b, and compare structural stability, dynamics, and thermodynamics with that of the TC5b variant, using replica exchange molecular dynamics (REMD). The TC10b variant was designed to have larger helical stability by the substitution of amino acids with greater alpha helical propensities in the N-terminal region. Experiments have shown TC10b to possess larger overall stability than TC5b. Simulations starting from unbiased, unfolded initial conditions are run for 1 μs per replica. The calculations show a higher melting temperature for TC10b than TC5b, and suggest a more ordered folded structure through the elimination of a substate found in the folded ensemble of TC5b. We model the difference in Gibbs free energy, ΔG(P,T), of folding using the bootstrap statistical method, which is used to calculate uncertainties associated with the thermodynamic parameters for both variants of the Trp-cage. We find that while the shape of the area for which the protein is stability folded is elliptical for TC5b, there is a degree of uncertainty associated with that of TC10b, with one model suggesting elliptical and another suggesting hyperbolic. This model suggests that at high pressures, TC5b can experience pressure denaturation, but TC10b may not.

Microsecond molecular dynamics simulation of guanidinium chloride induced unfolding of ubiquitin

All atom molecular dynamics simulations have been used to explore the atomic detail mechanism of guanidinium induced unfolding of the protein ubiquitin.