Additive transfer free energies of the peptide backbone unit that are independent of the model compound and the choice of concentration scale

With knowledge of individual transfer free energies of chemical groups that become newly exposed on protein denaturation and assuming the group transfer free energy contributions are additive, it should be possible to predict the stability of a protein in the presence of denaturant. Unfortunately, several unresolved issues have seriously hampered quantitative development of this transfer model for protein folding/unfolding. These issues include the lack of adequate demonstration that group transfer free energies (DeltaG(tr)) are additive and independent of the choice of model compound, the problem arising from dependence of DeltaG(tr) on concentration scales, the lack of knowledge of activity coefficients, and the validity of the mathematical constructs used in obtaining DeltaG(tr) values. Regarding transfer from water to 1 M concentrations of the naturally occurring osmolytes, trimethylamine-N-oxide (TMAO), sarcosine, betaine, proline, glycerol, sorbitol, sucrose, trehalose, and urea, using cyclic glycylglycine, zwitterionic glycine peptides, and N-acetylglycine amide peptides as models for the peptide backbone of proteins, we set out to address these issues and obtain DeltaG(tr) values for the peptide backbone unit. We demonstrate experimental approaches that obviate the choice of concentration scale and demonstrate additivity in DeltaG(tr) of the peptide backbone unit for all solvent systems studied. Evidence is presented to show that the DeltaG(tr) values are independent of the chemical model studied, and experimental conditions are given to illustrate when the mathematical constructs are valid and when activity coefficients can be ignored. Resolution of the long-standing issues that have stymied development of the transfer model now make it possible to design transfer experiments that yield reliable and quantitative values for the interactions between osmolyte-containing solvents and native and unfolded protein.     

Effect of guanidine hydrochloride and urea on the interaction of 6-thioguanine with human serum albumin: a spectroscopic and molecular dynamics based study

6-thioguanine (6-TG) is an antineoplastic, nucleobase guanine, purine analog drug belongs to thiopurine drug-family of antimetabolites. In the present study, we report an experimental approach towards interaction mechanism of 6-TG with human serum albumin (HSA) and examine the chemical stability of HSA in the presence of denaturants such as guanidine hydrochloride (GdnHCl) and urea. Interaction of 6-TG with HSA has been studied by various spectroscopic and spectropolarimeteric methods to investigate what short of binding occurs at physiological conditions. 6-TG binds in the hydrophobic cavity of subdomain IIA of HSA by static quenching mechanism which induces conformation alteration in the protein structure. That helpful for further study of denaturation process where change in secondary structures causes unfolding of protein that also responsible for severance of domain III from rest of the protein part. We have also performed molecular simulation and molecular docking study in the presence of denaturating agents to determine the binding property of 6-TG and the effect of denaturating agents on the structural activity of HSA. We had found that GdnHCl is more effective denaturating agent when compared to urea. Hence, this study provides straight evidence of the binding mechanism of 6-TG with HSA and the formation of intermediate or unfolding transition that causes unfolding of HSA.     

Refolding kinetics of cytochrome c(551) reveals a mechanistic difference between urea and guanidine

The energetic parameters for the folding of small globular proteins can be very different if derived from guanidine hydrochloride (GdnHCl) or urea denaturation experiments. A study of the equilibrium and kinetics of the refolding of wild-type (wt) cytochrome c(551) (cyt c(551)) from Pseudomonas aeruginosa and of two site-directed mutants (E70Q and E70V) shows that the nonionic nature of urea reveals the role of a salt bridge between residues E70 and K10 on the transition state, which is otherwise completely masked in GdnHCl experiments. Mixed denaturant refolding experiments allow us to conclude that the masking effect of GdnHCl is complete at fairly low GdnHCl concentrations ( congruent with 0.1 M). The fact that potassium chloride is unable to reproduce this quenching effect, together with the results obtained on the mutants, suggests a specific binding of the Gdn(+) cation, which involves the E70-K10 ion pair in wt cyt c(551).We propose, therefore, a simple kinetic test to obtain a mechanistic interpretation of nonlinear dependences of DeltaG(w) on GdnHCl concentration on the basis of kinetic refolding experiments in the presence of both denaturants.     

New evidence for the denaturant binding model.

Denaturation with guanidine hydrochloride (GdnHCl) or urea is one of the primary ways of measuring the conformational stability of proteins and comparing the stability of mutant proteins. Despite the widespread use of these two denaturants to provide quantitative data for the free energies of unfolding, the mode of action of these agents is not well understood. We are not even certain whether the action of these agents on proteins is direct and can be regarded as ligand binding, or indirect and involves a change in the properties of solvent (water) in the presence of GdnHCl and urea. In this paper, an extensive kinetic study of the inhibition of ribonuclease A and papain by urea has been performed. The results suggest that the effect of urea on activities of these enzymes can be well described by the denaturant binding model. The binding constants of urea determined by the present method are nearly identical to that determined from a variety of different studies on model compounds and proteins.     

Preferential interactions of urea with lysozyme and their linkage to protein denaturation

The interactions involved in the denaturation of lysozyme in the presence of urea were examined by thermal transition studies and measurements of preferential interactions of urea with the protein at pH 7.0, where it remains native up to 9.3 M urea, and at pH 2.0, where it undergoes a transition between 2.5 and 5.0 M urea. The destabilization of lysozyme by urea was found to follow the linear dependence on urea molar concentration, M(u), DeltaG(u)(o)=DeltaG(w)(o)-2.1 M(u), over the combined data, where DeltaG(u)(o) and DeltaG(w)(o) are the standard free energy changes of the N right harpoon over left harpoon D reaction in urea and water, respectively. Combination with the measured preferential binding gave the result that the increment of preferential binding, deltaGamma(23)=Gamma(23)(D)-Gamma(23)(N), is also linear in M(u). A temperature dependence study of preferential interactions permitted the evaluation of the transfer enthalpy, DeltaHmacr;(2,tr)(o), and entropy, DeltaSmacr;(2,tr)(o) of lysozyme from water into urea in both the native and denatured states. These values were found to be consistent with the enthalpy and entropy of formation of inter urea hydrogen bonds (Schellman, 1955; Kauzmann, 1959), with estimated values of DeltaHmacr;(2,tr)(o)=ca. -2.5 kcal mol(-1) and DeltaSmacr;(2,tr)(o)=ca. -7.0 e.u. per site. Analysis of the results led to the conclusion that the stabilization of the denatured form was predominantly by preferential binding to newly exposed peptide groups. Combination with the knowledge that stabilizing osmolytes act by preferential exclusion from peptide groups (Liu and Bolen, 1995) has led to the general conclusion that both the stabilization and destabilization of proteins by co-solvents are controlled predominantly by preferential interactions with peptide groups newly exposed on denaturation.     

Equilibrium unfolding of RNase Rs from Rhizopus stolonifer: pH dependence of chemical and thermal denaturation

The conformational stability of RNase Rs was determined with chemical and thermal denaturants over the pH range of 1-10. Equilibrium unfolding with urea showed that values of D(1/2) (5.7 M) and DeltaG(H(2)O) (12.8 kcal/mol) were highest at pH 5.0, its pI and the maximum conformational stability of RNase Rs was observed near pH 5.0. Denaturation with guanidine hydrochloride (GdnHCl), at pH 5.0, gave similar values of DeltaG(H(2)O) although GdnHCl was 2-fold more potent denaturant with D(1/2) value of 3.1 M. The curves of fraction unfolded (f(U)) obtained with fluorescence and CD measurements overlapped at pH 5.0. Denaturation of RNase Rs with urea in the pH range studied was reversible but the enzyme denatured irreversibly >pH 11.0. Thermal denaturation of RNase Rs was reversible in the pH range of 2.0-3.0 and 6.0-9.0. Thermal denaturation in the pH range 4.0-5.5 resulted in aggregation and precipitation of the protein above 55 degrees C. The aggregate was amorphous or disordered precipitate as observed in TE micrographs. Blue shift in emission lambda(max) and enhancement of fluorescence intensity of ANS at 70 degrees C indicated the presence of solvent exposed hydrophobic surfaces as a result of heat treatment. Aggregation could be prevented partially with alpha-cyclodextrin (0.15 M) and completely with urea at concentrations >3 M. Aggregation was probably due to intermolecular hydrophobic interaction favored by minimum charge-charge repulsion at the pI of the enzyme. Both urea and temperature-induced denaturation studies showed that RNase Rs unfolds through a two-state F right arrow over left arrow U mechanism. The pH dependence of stability described by DeltaG(H(2)O) (urea) and DeltaG (25 degrees C) suggested that electrostatic interactions among the charged groups make a significant contribution to the conformational stability of RNase Rs. Since RNase Rs is a disulfide-containing protein, the major element for structural stability are the covalent disulfide bonds.     

Measuring the stability of partly folded proteins using TMAO

Standard methods for measuring free energy of protein unfolding by chemical denaturation require complete folding at low concentrations of denaturant so that a native baseline can be observed. Alternatively, proteins that are completely unfolded in the absence of denaturant can be folded by addition of the osmolyte trimethylamine N-oxide (TMAO), and the unfolding free energy can then be calculated through analysis of the refolding transition. However, neither chemical denaturation nor osmolyte-induced refolding alone is sufficient to yield accurate thermodynamic unfolding parameters for partly folded proteins, because neither method produces both native and denatured baselines in a single transition. Here we combine urea denaturation and TMAO stabilization as a means to bring about baseline-resolved structural transitions in partly folded proteins. For Barnase and the Notch ankyrin domain, which both show two-state equilibrium unfolding, we found that DeltaG degrees for unfolding depends linearly on TMAO concentration, and that the sensitivity of DeltaG degrees to urea (the m-value) is TMAO independent. This second observation confirms that urea and TMAO exert independent effects on stability over the range of cosolvent concentrations required to bring about baseline-resolved structural transitions. Thermodynamic parameters calculated using a global fit that assumes additive, linear dependence of DeltaG degrees on each cosolvent are similar to those obtained by standard urea-induced unfolding in the absence of TMAO. Finally, we demonstrate the applicability of this method to measurement of the free energy of unfolding of a partly folded protein, a fragment of the full-length Notch ankyrin domain.     

The peculiar nature of the guanidine hydrochloride-induced two-state denaturation of staphylococcal nuclease: a calorimetric study.

This work determines the ratio of DeltaH(vH) /DeltaH(cal) for staphylococcal nuclease (SN) denaturation in guanidine hydrochloride (GdnHCl) to test whether GdnHCl-induced denaturation is two-state. Heats of mixing of SN as a function of [GdnHCl] were determined at pH 7.0 and 25 degrees C. The resulting plot of DeltaH(mix) vs [GdnHCl] exhibits a sigmoid shaped curve with linear pre- and post-denaturational base lines. Extending the pre- and post-denaturational lines to zero [GdnHCl] gives a calorimetric DeltaH (DeltaH(cal)) of 24.1 +/- 1.0 kcal/mol, for SN denaturation in the limit of zero GdnHCl concentration. Guanidine hydrochloride-induced denaturation Gibbs energy changes in the limit of zero denaturant concentration (DeltaG degrees (N)(-)(D)) at pH 7. 0 were determined for SN from fluorescence measurements at fixed temperatures over the range from 15 to 35 degrees C. Analysis of the resulting temperature-dependent DeltaG degrees (N)(-)(D) data defines a van't Hoff denaturation enthalpy change (DeltaH(vH)) of 26. 4 +/- 2.8 kcal/mol. The model-dependent van't Hoff DeltaH(vH) divided by the model-independent DeltaH(cal) gives a ratio of 1.1 +/- 0.1 for DeltaH(vH)/DeltaH(cal), a result that rules out the presence of thermodynamically important intermediate states in the GdnHCl-induced denaturation of SN. The likelihood that GdnHCl-induced SN denaturation involves a special type of two-state denaturation, known as a variable two-state process, is discussed in terms of the thermodynamic implications of the process.     

Interactions of myoglobin with urea and some alkylureas. II. Calorimetric and circular dichroic studies

The interactions of myoglobin with urea, methyl-, N,N'-dimethyl- and ethylurea were studied by means of calorimetry and circular dichroism (CD). The enthalpies of transfer from water to aqueous denaturant solutions are positive for the alkylureas and negative for urea. The difference is due to the presence of hydrophobic groups in the alkylureas. Gibbs free energies of transfer for urea solutions were obtained from preferential binding data determined previously. An attempt is made to interpret the values of the thermodynamic quantities in terms of various interactions between protein and denaturant. Analysis of the far-ultraviolet CD spectra reveals some differences in the denaturing activity of urea and the alkylureas, the latter being stronger denaturants than urea. Myoglobin displays relatively high stability towards these denaturants since concentrations above 5 M are needed for achieving major conformational changes.     

The effect of base sequence on the stability of RNA and DNA single base bulges

Forty-eight RNA duplexes were constructed that contained all common single base bulges at six different locations. The stabilities of the RNAs were determined by temperature gradient gel electrophoresis (TGGE). The relative stability of a single base bulge was dependent on both base identity and the nearest neighbor context. The single base bulges were placed into two categories. A bulged base with no identical neighboring base was defined as a Group I base bulge. Group II-bulged bases had at least one neighboring base identical to it. Group II bulges were generally more stable than Group I bulges in the same nearest neighbor environments. This indicates that position degeneracy of an unpaired base enhances stability. Differences in the mobility transition temperatures between the RNA fragments with bulges and the completely base-paired reference RNAs were related to free energy differences. Simple models for estimating the free energy contribution of single base bulges were evaluated from the free energy difference data. The contribution of a Group I bulge 5'-(XNZ)-3'.5'-(Z'-X')-3' where N is the unpaired base and X.X' and Z.Z' the neighboring base pairs, could be well-represented (+/-0.34 kcal/mol) by the equation, DeltaG((X)(N)()(Z))(.)((Z)(')(-)(X)(')()) = 3.11 + 0. 40DeltaG(s)()((XZ))(.)((Z)(')(X)(')()). DeltaG(s)()((XZ))(. )((Z)(')(X)(')()) is the stacking energy of the closing base pair doublet. By adding a constant term, delta = -0.3 kcal/mol, to the right side of the above equation, free energies of Group II bulges could also be predicted with the same accuracy. The term delta represents the stabilizing effect due to position degeneracy. A similar equation/model was applied to previous data from 32 DNA fragments with single base bulges. It predicted the free energy differences with a similar standard deviation.     

A stability transition at mildly acidic pH in the alpha-hemolysin (alpha-toxin) from Staphylococcus aureus.

The effects of mildly acidic conditions on the free energy of unfolding (DeltaG(u)(buff)) of the pore-forming alpha-hemolysin (alphaHL) from Staphylococcus aureus were assessed between pH 5.0 and 7.5 by measuring intrinsic tryptophan fluorescence, circular dichroism and elution time in size exclusion chromatography during urea denaturation. Decreasing the pH from 7.0 to 5.0 reduced the calculated DeltaG(u)(buff) from 8.9 to 4.2 kcal mol(-1), which correlates with an increased rate of pore formation previously observed over the same pH range. It is proposed that the lowered surface pH of biological membranes reduces the stability of alphaHL thereby modulating the rate of pore formation.     

A computer simulation of functional group contributions to free energy in water and a DPPC lipid bilayer

A series of all-atom molecular dynamics simulations has been performed to evaluate the contributions of various functional groups to the free energy of solvation in water and a dipalmitoylphospatidylcholine lipid bilayer membrane and to the free energies of solute transfer (Delta(DeltaG(o))X) from water into the ordered-chain interior of the bilayer. Free energies for mutations of the alpha-H atom in p-toluic acid to six different substituents (-CH3, -Cl, -OCH3, -CN, -OH, -COOH) were calculated by a combined thermodynamic integration and perturbation method and compared to literature results from vapor pressure measurements, partition coefficients, and membrane transport experiments. Convergence of the calculated free energies was indicated by substantial declines in standard deviations for the calculated free energies with increased simulation length, by the independence of the ensemble-averaged Boltzmann factors to simulation length, and the weak dependence of hysteresis effects on simulation length over two different simulation lengths and starting from different initial configurations. Calculated values of Delta(DeltaG(o))X correlate linearly with corresponding values obtained from lipid bilayer transport experiments with a slope of 1.1 and from measurements of partition coefficients between water and hexadecane or decadiene, with slopes of 1.1 and 0.9, respectively. Van der Waals interactions between the functional group of interest and the acyl chains in the ordered chain region account for more than 95% of the overall potential energy of interaction. These results support the view that the ordered chain region within the bilayer interior is the barrier domain for transport and that solvation interactions within this region resemble those occurring in a nonpolar hydrocarbon.     

Kinetic inactivation study of mushroom tyrosinase: intermediate detection by denaturants

The unfolding and inhibition study of mushroom tyrosinase have been studied in the presence of different denaturants such as sodium dodecyl sulfate (SDS), guanidine hydrochloride (GdnHCl), and urea. The kinetic two-phase rate constants were commonly measured from semilogarithmic plots of the activity versus time, which resolved into two straight lines, indicating that the inactivation process consisted of fast and slow phases as a first-order reaction. This result also implied that transient partially folded intermediate existed during tyrosinase unfolding pathway. Mushroom tyrosinase had different behaviors to denaturants in regard with: noncooperative binding manner by SDS while cooperative interactions by GdnHCl and urea; in equilibrium state, SDS-micelle never completely inactivated enzyme activity while GdnHCl has single step denaturation and urea induced a typical transition-like process. Various kinetic parameters for each denaturant were calculated and the possible unfolding pathway scheme was discussed.     

Measurement of thermodynamic parameters for hydrophobic mismatch 2: intermembrane transfer of a transmembrane helix

Hydrophobic matching between proteins and lipids is essential for the thermodynamic stability of integral membrane proteins. However, there is no direct thermodynamic information available about the intermembrane transfer of proteins between membranes with different hydrophobic thicknesses, which is crucial for understanding hydrophobic mismatch. This article reports the complete set of thermodynamic parameters (DeltaG, DeltaH, DeltaS, and DeltaC(p)) for the intermembrane transfer of the inert hydrophobic model transmembrane helix NBD-(AALALAA)(3)-NH(2) (NBD: 7-nitro-2-1,3-benzoxadiazol-4-yl), which is exchangeable between vesicles, from 1-palmitoyl-2-oleoyl-sn-glycero-3-phosphocholine (POPC) to dimonounsaturated-phosphocholine lipid bilayers with different hydrophobic thicknesses (C14-C22) at 37-58 degrees C. The transfer free energies were calculated from equilibrium values of the extent of helix transfer from donor to acceptor lipid vesicles, as monitored by a decrease in fluorescence resonance energy transfer from the NBD group to a lipid-labeled Rhodamine in the donor upon transfer to the quencher-free acceptor. Under hydrophobic mismatch conditions up to approximately 7 A, the helix partitioning became unfavorable up to +7 kJ mol(-)(1), hampered by an increase in entropic (up to +20 kJ mol(-)(1)) and enthalpic (up to +66 kJ mol(-)(1)) terms in thinner and thicker membranes, respectively. Together with the results that H/D exchange at the membrane interface was accelerated in thinner membranes the obtained thermodynamic parameters were reasonably explained assuming that hydrophobic mismatch induces aqueous exposure or membrane burial of the helix termini, resulting in excess energies originating from the hydration of terminal hydrophobic residues or the unfavorable Born energy of terminal partial charges of the helix macrodipole.     

Characterization of intra-molecular distances and site-specific dynamics in chemically unfolded barstar: evidence for denaturant-dependent non-random structure

The structure and dynamics of the unfolded form of a protein are expected to play critical roles in determining folding pathways. In this study, the urea and guanidine hydrochloride (GdnHCl)-unfolded forms of the small protein barstar were explored by time-resolved fluorescence techniques. Barstar was labeled specifically with thionitrobenzoate (TNB), by coupling it to the thiol side-chain of a cysteine residue at one of the following positions on the sequence: 14, 25, 40, 42, 62, 82 and 89, in single cysteine-containing mutant proteins. Seven intra-molecular distances (R(DA)) under unfolding conditions were estimated from measurements of time-resolved fluorescence resonance energy transfer between the donor Trp53 and the non-fluorescent acceptor TNB coupled to one of the seven cysteine side-chains. The unfolded protein chain expands with an increase in the concentration of the denaturants. The extent of expansion was found to be non-uniform, with different intra-molecular distances expanding to different extents. In general, shorter distances were found to expand less when compared to longer spans. The extent of expansion was higher in the case of GdnHCl when compared to urea. A comparison of the measured values of R(DA) with those derived from a model based on excluded volume, revealed that while shorter spans showed good agreement, the experimental values of R(DA) of longer spans were smaller when compared to the theoretical values. Sequence-specific flexibility of the polypeptide was determined by time-resolved fluorescence anisotropy decay measurements on acrylodan or 1,5-IAEDANS labeled single cysteine-containing proteins under unfolding conditions. Rotational dynamics derived from these measurements indicated that the level of flexibility increased with increase in the concentration of denaturants and showed a graded increase towards the C-terminal end. Taken together, these results appear to indicate the presence of specific non-random coil structures and show that the deviation from random coil structure is different for the two denaturants.     

Bromophenol blue binding as a probe to study urea and guanidine hydrochloride denaturation of bovine serum albumin

Urea and guanidine hydrochloride (GdnHCl) denaturation of bovine serum albumin (BSA) were investigated using bromophenol blue (BPB) binding as a probe. Addition of BPB to BSA produced an absorption difference spectrum in the wavelength range, 525-675 nm with a minimum at 587 nm and a maximum at 619 nm. The magnitude of absorption difference (DeltaAbs.) at 619 nm decreased on increasing urea/GdnHCl concentration and followed the denaturation curve. The denaturation was found to be a two-state, single-step transition. The transitions started at 1.75 and 0.875 M and completed at 6.5 and 3.25 M with the mid point occurring around 4.0 and 1.5 M urea and GdnHCl concentrations, respectively. The value of free energy of stabilization, DeltaGDH2O as determined from urea and GdnHCl denaturation curves was found to be 4041 and 4602 cal/mol, respectively. Taken together, these results suggest that BPB binding can be used as a probe to study urea and GdnHCl denaturation of BSA.     

Characterization of the unfolding process of lipocalin-type prostaglandin D synthase

We found that low concentrations of guanidine hydrochloride (GdnHCl, <0.75 M) or urea (<1.5 M) enhanced the enzyme activity of lipocalin-type prostaglandin (PG) D synthase (L-PGDS) maximally 2.5- and 1.6-fold at 0.5 M GdnHCl and 1 M urea, respectively. The catalytic constants in the absence of denaturant and in the presence of 0.5 M GdnHCl or 1 m urea were 22, 57, and 30 min(-1), respectively, and the K(m) values for the substrate, PGH(2), were 2.8, 8.3, and 2.3 microm, respectively, suggesting that the increase in the catalytic constant was mainly responsible for the activation of L-PGDS. The intensity of the circular dichroism (CD) spectrum at 218 nm, reflecting the beta-sheet content, was also increased by either denaturant in a concentration-dependent manner, with the maximum at 0.5 M GdnHCl or 1 M urea. By plotting the enzyme activities against the ellipticities at 218 nm of the CD spectra of L-PGDS in the presence or absence of GdnHCl or urea, we found two states in the reversible folding process of L-PGDS: one is an activity-enhanced state and the other, an inactive state. The NMR analysis of L-PGDS revealed that the hydrogen-bond network was reorganized to be increased in the activity-enhanced state formed in the presence of 0.5 M GdnHCl or 1 m urea and to be decreased but still remain in the inactive intermediate observed in the presence of 2 M GdnHCl or 4 M urea. Furthermore, binding of the nonsubstrate ligands, bilirubin or 13-cis-retinal, to L-PGDS changed from a multistate mode in the native form of L-PGDS to a simple two-state mode in the activity-enhanced form, as monitored by CD spectra of the bound ligands. Therefore, L-PGDS is a unique protein whose enzyme activity and ligand-binding property are biphasically altered during the unfolding process by denaturants.