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

Preferential interactions between small solutes and the protein backbone: a computational analysis

To improve our understanding of the effects of small solutes on protein stability, we conducted atomistic simulations to quantitatively characterize the interactions between two broadly used small solutes, urea and glycine betaine (GB), and a triglycine peptide, which is a good model for a protein backbone. Multiple solute concentrations were analyzed, and each solute-peptide-water ternary system was studied with approximately 200-300 ns of molecular dynamics simulations with the CHARMM force field. The comparison between calculated preferential interaction coefficients (Gamma(23)) and experimentally measured values suggests that semiquantitative agreement with experiments can be obtained if care is exercised to balance interactions among the solute, protein, and water. On the other hand, qualitatively incorrect (i.e., wrong sign in Gamma(23)) results can be obtained if a solute model is constructed by directly taking parameters for chemically similar groups from an existing force field. Such sensitivity suggests that small solute thermodynamic data can be valuable in the development of accurate force field models of biomolecules. Further decomposition of Gamma(23) into group contributions leads to additional insights regarding the effects of small solutes on protein stability. For example, use of the CHARMM force field predicts that urea preferentially interacts with not only amide groups in the peptide backbone but also aliphatic groups, suggesting a role for these interactions in urea-induced protein denaturation; quantitatively, however, it is likely that the CHARMM force field overestimates the interaction between urea and aliphatic groups. The results with GB support a simple thermodynamic model that assumes additivity of preferential interaction between GB and various biomolecular surfaces.

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.

Urea, but not guanidinium, destabilizes proteins by forming hydrogen bonds to the peptide group

The mechanism by which urea and guanidinium destabilize protein structure is controversial. We tested the possibility that these denaturants form hydrogen bonds with peptide groups by measuring their ability to block acid- and base-catalyzed peptide hydrogen exchange. The peptide hydrogen bonding found appears sufficient to explain the thermodynamic denaturing effect of urea. Results for guanidinium, however, are contrary to the expectation that it might H-bond. Evidently, urea and guanidinium, although structurally similar, denature proteins by different mechanisms.

Dewetting and hydrophobic interaction in physical and biological systems

Hydrophobicity manifests itself differently on large and small length scales. This review focuses on large-length-scale hydrophobicity, particularly on dewetting at single hydrophobic surfaces and drying in regions bounded on two or more sides by hydrophobic surfaces. We review applicable theories, simulations, and experiments pertaining to large-scale hydrophobicity in physical and biomolecular systems and clarify some of the critical issues pertaining to this subject. Given space constraints, we cannot review all the significant and interesting work in this active field.

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.

To perform their physiological function, proteins have to fold into their characteristic three-dimensional structure. While the folded state is stable under physiological conditions, changes in the solvent can destabilize the folded state and even induce denaturation. One of the most commonly used denaturants is urea. Despite its widespread use to study protein folding and stability, however, the molecular mechanism and particularly the driving forces of urea-induced protein denaturation are not yet understood. Two mechanisms have been suggested, according to which denaturation is driven either by polar interactions via hydrogen bonds or by hydrophobic interactions with apolar amino acids. By systematically varying urea polarity and quantifying the interactions of the solvent molecules with all amino acids of the protein, the present simulation study reveals that it is mainly the apolar interactions that drive denaturation. Our results suggest a coherent microscopic picture for urea-induced denaturation and bear more general implications for protein stability in other environments, e.g., in chaperone-assisted folding.

Urea denaturation by stronger dispersion interactions with proteins than water implies a 2-stage unfolding

The mechanism of denaturation of proteins by urea is explored by using all-atom microseconds molecular dynamics simulations of hen lysozyme generated on BlueGene/L. Accumulation of urea around lysozyme shows that water molecules are expelled from the first hydration shell of the protein. We observe a 2-stage penetration of the protein, with urea penetrating the hydrophobic core before water, forming a "dry globule." The direct dispersion interaction between urea and the protein backbone and side chains is stronger than for water, which gives rise to the intrusion of urea into the protein interior and to urea's preferential binding to all regions of the protein. This is augmented by preferential hydrogen bond formation between the urea carbonyl and the backbone amides that contributes to the breaking of intrabackbone hydrogen bonds. Our study supports the "direct interaction mechanism" whereby urea has a stronger dispersion interaction with protein than water.

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.

Recent applications of Kirkwood-Buff theory to biological systems

The effect of cosolvents on biomolecular equilibria has traditionally been rationalized using simple binding models. More recently, a renewed interest in the use of Kirkwood-Buff (KB) theory to analyze solution mixtures has provided new information on the effects of osmolytes and denaturants and their interactions with biomolecules. Here we review the status of KB theory as applied to biological systems. In particular, the existing models of denaturation are analyzed in terms of KB theory, and the use of KB theory to interpret computer simulation data for these systems is discussed.

Single-mutation-induced stability loss in protein lysozyme

Recent NMR experiments have revealed that a single residue mutation W62G on protein hen's-egg white lysozyme can cause a dramatic loss of long-range interactions and protein stability; however, the molecular mechanism for this surprising phenomenon is not completely clear. In this mini-review, we have summarized some of our recent work on the molecular mechanism with large-scale molecular modelling, and also utilized a new wavelet method to analyse the local structural clusters present in both the wild-type and mutant folding trajectories. These extensive MD (Molecular Dynamics) simulations (10+ μs) were performed in 8 M urea, mimicking the experimental condition. Detailed analyses revealed that the Trp62 residue is the key to a co-operative long-range interaction within the wild-type protein: it acts as a bridge between neighbouring basic residues, mainly arginine residues, through π-type hydrogen bonds or π-cation interactions to form an Arg-Trp-Arg ‘sandwich-like’ local structure. The local cluster near Trp62 further extends its interaction to other clusters, such as the one near Trp111, through Arg112, which is involved in such an Arg-Trp-Arg bridging structure, thus achieving the long-range interactions for the wild-type. On the other hand, the mutant does not have this bridging effect and forms much less local clusters or contacts, and therefore results in a much less stable structure. Overall, these findings not only support the general conclusions of the experiment, but also provide a detailed but somewhat different molecular picture of the disruption of the long-range interactions.

Anatomy of energetic changes accompanying urea-induced protein denaturation

Because of its protein-denaturing ability, urea has played a pivotal role in the experimental and conceptual understanding of protein folding and unfolding. The measure of urea's ability to force a protein to unfold is given by the m value, an experimental quantity giving the free energy change for unfolding per molar urea. With the aid of Tanford's transfer model [Tanford C (1964) J Am Chem Soc 86:2050-2059], we use newly obtained group transfer free energies (GTFEs) of protein side-chain and backbone units from water to 1 M urea to account for the m value of urea, and the method reveals the anatomy of protein denaturation in terms of residue-level free energy contributions of groups newly exposed on denaturation. The GTFEs were obtained by accounting for solubility and activity coefficient ratios accompanying the transfer of glycine from water to 1 M urea. Contrary to the opinions of some researchers, the GTFEs show that urea does not denature proteins through favorable interactions with nonpolar side chains; what drives urea-induced protein unfolding is the large favorable interaction of urea with the peptide backbone. Although the m value is said to be proportional to surface area newly exposed on denaturation, only approximately 25% of the area favorably contributes to unfolding (because of newly exposed backbone units), with approximately 75% modestly opposing urea-induced denaturation (originating from side-chain exposure). Use of the transfer model and newly determined GTFEs achieves the long-sought goal of predicting urea-dependent cooperative protein unfolding energetics at the level of individual amino acid residues.

Interactions between Hydrophobic and Ionic Solutes in Aqueous Guanidinium Chloride and Urea Solutions: Lessons for Protein Denaturation Mechanism

In order to clarify the mechanism of denaturant-induced unfolding of proteins we have calculated the interactions between hydrophobic and ionic species in aqueous guanidinium chloride and urea solutions using molecular dynamics simulations. Hydrophobic association is not significantly changed in urea or guanidinium chloride solutions. The strength of interaction between ion pairs is greatly diminished by the guanidinium ion. Although the changes in electrostatic interactions in urea are small, examination of structures, using appropriate pair functions, of urea and water around the solutes show strong hydrogen bonding between urea's carbonyl oxygen and the positively charged solute. Our results strongly suggest protein denaturation occurs by the direct interaction model according to which the most commonly used denaturants unfold proteins by altering electrostatic interactions either by solvating the charged residues or by engaging in hydrogen bonds with the protein backbone. To further validate the direct interaction model we show that, in urea and guanidinium chloride solutions, unfolding of an unusually stable helix (H1) from mouse PrPC (residues 144-153) occurs by hydrogen bonding of denaturants to charged side chains and backbone carbonyl groups.

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.

The interaction of guanidinium ions with a model peptide

In addition to promoting unfolded protein states, the denaturants urea and guanidinium (Gdm(+)) accumulate at the surface of folded proteins at subdenaturing concentrations, a phenomenon that correlates with their denaturant activities. The enhanced accumulation of Gdm(+) relative to urea indicates different binding modes, or additional binding sites, for Gdm(+), and we recently proposed potential binding modes to protein functional groups for Gdm(+) based on the determination of the weak hydration properties of this complex cation. Here we describe molecular dynamics simulations of a model helical peptide, melittin, in a 3 M solution of GdmCl, to identify potential interactions with amino-acid side chains in a nondenatured polypeptide surface. The simulations indicate that Gdm(+) can interact with a number of planar amino-acid side chains (Arg, Trp, Gln) in a stacking manner, as well as more weakly with hydrophobic surfaces composed of aliphatic side chains, and that these interactions result in enhanced number densities of Gdm(+) at certain locations on the peptide surface. These observations provide molecular scale insight into the accumulation of Gdm(+) at protein surfaces that has previously been observed experimentally.

Destruction of long-range interactions by a single mutation in lysozyme

We propose a mechanism, based on a > or =10-micros molecular dynamics simulation, for the surprising misfolding of hen egg-white lysozyme caused by a single mutation (W62G). Our simulations of the wild-type and mutant lysozymes in 8 M urea solution at biological temperature (with both pH 2 and 7) reveal that the mutant structure is much less stable than that of the wild type, with the mutant showing larger fluctuations and less native-like contacts. Analysis of local contacts reveals that the Trp-62 residue is the key to a cooperative long-range interaction within the wild type, where it acts like a bridge between two neighboring basic residues. Thus, a native-like cluster or nucleation site can form near these residues in the wild type but not in the mutant. The time evolution of the secondary structure also exhibits a quicker loss of the beta-sheets in the mutant than in the wild type, whereas some of the alpha-helices persist during the entire simulation in both the wild type and the mutant in 8 M urea (even though the tertiary structures are basically all gone). These findings, while supporting the general conclusions of a recent experimental study by Dobson and coworkers [Klein-Seetharam J, Oikama M, Grimshaw SB, Wirmer J, Duchardt E, Ueda T, Imoto T, Smith LJ, Dobson CM, Schwalbe H (2002) Science 295:1719-1722], provide a detailed but different molecular picture of the misfolding mechanism.

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.

Protein under pressure: Molecular dynamics simulation of the arc repressor

Experimental nuclear magnetic resonance results for the Arc Repressor have shown that this dimeric protein dissociates into a molten globule at high pressure. This structural change is accompanied by a modification of the hydrogen-bonding pattern of the intermolecular beta-sheet: it changes its character from intermolecular to intramolecular with respect to the two monomers. Molecular dynamics simulations of the Arc Repressor, as a monomer and a dimer, at elevated pressure have been performed with the aim to study this hypothesis and to identify the major structural and dynamical changes of the protein under such conditions. The monomer appears less stable than the dimer. However, the complete dissociation has not been seen because of the long timescale needed to observe this phenomenon. In fact, the protein structure altered very little when increasing the pressure. It became slightly compressed and the dynamics of the side-chains and the unfolding process slowed down. Increasing both, temperature and pressure, a tendency of conversion of intermolecular into intramolecular hydrogen bonds in the beta-sheet region has been detected, supporting the mentioned hypothesis. Also, the onset of denaturation of the separated chains was observed.

WATER MEDIATION IN PROTEIN FOLDING AND MOLECULAR RECOGNITION

Water is essential for life in many ways, and without it biomolecules might no longer truly be biomolecules. In particular, water is important to the structure, stability, dynamics, and function of biological macromolecules. In protein folding, water mediates the collapse of the chain and the search for the native topology through a funneled energy landscape. Water actively participates in molecular recognition by mediating the interactions between binding partners and contributes to either enthalpic or entropic stabilization. Accordingly, water must be included in recognition and structure prediction codes to capture specificity. Thus water should not be treated as an inert environment, but rather as an integral and active component of biomolecular systems, where it has both dynamic and structural roles. Focusing on water sheds light on the physics and function of biological machinery and self-assembly and may advance our understanding of the natural design of proteins and nucleic acids.

On the edge of the denaturation process: Application of X-ray diffraction to barnase and lysozyme cross-linked crystals with denaturants in molar concentrations

Structural data about the early step of protein denaturation were obtained from cross-linked crystals for two small proteins: barnase and lysozyme. Several denaturant agents like urea, bromoethanol or thiourea were used at increasing concentrations up to a limit leading to crystal disruption (>or=2 to 6 M). Before the complete destruction of the crystal order started, specific binding sites were observed at the protein surfaces, an indication that the preliminary step of denaturation is the disproportion of intermolecular polar bonds to the benefit of the agent "parasiting" the surface. The analysis of the thermal factors first agree with a stabilization effect at low or moderate concentration of denaturants rapidly followed by a destabilization at specific weak points when the number of sites increase (overflooding effect).

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.

Nanomechanical force measurements of gliadin protein interactions

The strength and nature of interactions between monomeric gliadin proteins involving alpha-alpha, omega-omega, and alpha-omega interactions in 0.01M acetic acid, and the effect of urea has been investigated. It was shown by means of nanomechanical force measurements that the stretching events in the separation curve after adhesive phenomena originated from proteins. These stretching events displayed different responses of the alpha- and omega-gliadins to urea. While 2M urea caused the more globular alpha-gliadins to unfold, the beta-turn-rich omega-gliadins remained fairly stable even in 8M urea. This suggests different roles for gliadins in the formation of dough; while the omega-gliadins are still in a compact structure being responsible for the viscous flow, the alpha-gliadins have already started to participate in forming the network in dough.

Effect of urea on peptide conformation in water: molecular dynamics and experimental characterization

Molecular dynamics simulations of a ribonuclease A C-peptide analog and a sequence variant were performed in water at 277 and 300 K and in 8 M urea to clarify the molecular denaturation mechanism induced by urea and the early events in protein unfolding. Spectroscopic characterization of the peptides showed that the C-peptide analog had a high alpha-helical content, which was not the case for the variant. In the simulations, interdependent side-chain interactions were responsible for the high stability of the alpha-helical C-peptide analog in the different solvents. The other peptide displayed alpha-helical unwinding that propagated cooperatively toward the N-terminal. The conformations sampled by the peptides depended on their sequence and on the solvent. The ability of water molecules to form hydrogen bonds to the peptide as well as the hydrogen bond lifetimes increased in the presence of urea, whereas water mobility was reduced near the peptide. Urea accumulated in excess around the peptide, to which it formed long-lived hydrogen bonds. The unfolding mechanisms induced by thermal denaturation and by urea are of a different nature, with urea-aqueous solutions providing a better peptide solvation than pure water. Our results suggest that the effect of urea on the chemical denaturation process involves both the direct and indirect mechanisms.

The heme iron coordination of unfolded ferric and ferrous cytochrome c in neutral and acidic urea solutions. Spectroscopic and electrochemical studies

The heme iron coordination of unfolded ferric and ferrous cytochrome c in the presence of 7-9 M urea at different pH values has been probed by several spectroscopic techniques including magnetic and natural circular dichroism (CD), electrochemistry, UV-visible (UV-vis) absorption and resonance Raman (RR). In 7-9 M urea at neutral pH, ferric cytochrome c is found to be predominantly a low spin bis-His-ligated heme center. In acidic 9 M urea solutions the UV-vis and near-infrared (NIR) magnetic circular dichroism (MCD) measurements have for the first time revealed the formation of a high spin His/H(2)O complex. The pK(a) for the neutral to acidic conversion is 5.2. In 9 M urea, ferrous cytochrome c is shown to retain its native ligation structure at pH 7. Formation of a five-coordinate high spin complex in equilibrium with the native form of ferrous cytochrome c takes place below the pK(a) 4.8. The formal redox potential of the His/H(2)O complex of cytochrome c in 9 M urea at pH 3 was estimated to be -0.13 V, ca. 100 mV more positive than E degrees ' estimated for the bis-His complex of cytochrome c in urea solution at pH 7.

Water and urea interactions with the native and unfolded forms of a beta-barrel protein

A fundamental understanding of protein stability and the mechanism of denaturant action must ultimately rest on detailed knowledge about the structure, solvation, and energetics of the denatured state. Here, we use (17)O and (2)H magnetic relaxation dispersion (MRD) to study urea-induced denaturation of intestinal fatty acid-binding protein (I-FABP). MRD is among the few methods that can provide molecular-level information about protein solvation in native as well as denatured states, and it is used here to simultaneously monitor the interactions of urea and water with the unfolding protein. Whereas CD shows an apparently two-state transition, MRD reveals a more complex process involving at least two intermediates. At least one water molecule binds persistently (with residence time >10 nsec) to the protein even in 7.5 M urea, where the large internal binding cavity is disrupted and CD indicates a fully denatured protein. This may be the water molecule buried near the small hydrophobic folding core at the D-E turn in the native protein. The MRD data also provide insights about transient (residence time <1 nsec) interactions of urea and water with the native and denatured protein. In the denatured state, both water and urea rotation is much more retarded than for a fully solvated polypeptide. The MRD results support a picture of the denatured state where solvent penetrates relatively compact clusters of polypeptide segments.

The molecular basis for the chemical denaturation of proteins by urea

Molecular dynamics simulations of the protein chymotrypsin inhibitor 2 in 8 M urea at 60 degrees C were undertaken to investigate the molecular basis of chemical denaturation. The protein unfolded rapidly under these conditions, but it retained its native structure in a control simulation in water at the same temperature. The overall process of unfolding in urea was similar to that observed in thermal denaturation simulations above the protein's T(m) of 75 degrees C. The first step in unfolding was expansion of the hydrophobic core. Then, the core was solvated by water and later by urea. The denatured structures in both urea and at high temperature contained residual native helical structure, whereas the beta-structure was completely disrupted. The average residence time for urea around hydrophilic groups was six times greater than around hydrophobic residues and in all cases greater than the corresponding water residence times. Water self-diffusion was reduced 40% in 8 M urea. Urea altered water structure and dynamics, thereby diminishing the hydrophobic effect and encouraging solvation of hydrophobic groups. In addition, through urea's weakening of water structure, water became free to compete with intraprotein interactions. Urea also interacted directly with polar residues and the peptide backbone, thereby stabilizing nonnative conformations. These simulations suggest that urea denatures proteins via both direct and indirect mechanisms.

The core of the tetrameric mycobacterial porin MspA is an extremely stable beta-sheet domain

MspA is the major porin of Mycobacterium smegmatis mediating the exchange of hydrophilic solutes across the cell wall and is the prototype of a new family of tetrameric porins with a single central pore of 10 nm in length. Infrared and circular dichroism spectroscopy revealed that MspA consists mainly of antiparallel beta-strands organized in a coherent domain. Heating to 92 and 112 degrees C was required to dissociate the MspA tetramer and to unfold the beta-sheet domain in the monomer, respectively. The stability of the MspA tetramer exceeded the remarkable stability of the porins of Gram-negative bacteria for every condition tested and was not reduced in the presence of 2% SDS and at any pH from 2 to 14. These results indicated that the interactions between the MspA subunits are different from those in the porins of Gram-negative bacteria and are discussed in the light of a channel-forming beta-barrel as a core structure of MspA. Surprisingly, the channel activity of MspA in 2% SDS and 7.6 m urea at 50 degrees C was reduced 13- and 30-fold, respectively, although the MspA tetramer and the beta-sheet domain were stable under those conditions. Channel closure by conformational changes of extracellular loops under those conditions is discussed to explain these observations. This study presents the first experimental evidence that outer membrane proteins not only from Gram-negative bacteria but also from mycobacteria are beta-sheet proteins and demonstrates that MspA constitutes the most stable transmembrane channel protein known so far. Thus, MspA may be of special interest for biotechnological applications.

The dominant interaction between peptide and urea is electrostatic in nature: a molecular dynamics simulation study

The conformational equilibrium of a blocked valine peptide in water and aqueous urea solution is studied using molecular dynamics simulations. Pair correlation functions indicate enhanced concentration of urea near the peptide. Stronger hydrogen bonding of urea-peptide compared to water-peptide is observed with preference for helical conformation. The potential of mean force, computed using umbrella sampling, shows only small differences between urea and water solvation that are difficult to quantify. The changes in solvent structure around the peptide are explained by favorable electrostatic interactions (hydrogen bonds) of urea with the peptide backbone. There is no evidence for significant changes in hydrophobic interactions in the two conformations of the peptide in urea solution. Our simulations suggest that urea denatures proteins by preferentially forming hydrogen bonds to the peptide backbone, reducing the barrier for exposing protein residues to the solvent, and reaching the unfolded state.

Molecular dynamics simulation reveals a surface salt bridge forming a kinetic trap in unfolding of truncated Staphylococcal nuclease

Surface salt bridges are ubiquitous in globular proteins. Their contribution to protein stability has been extensively debated in the past decade. Here, molecular dynamics simulations are performed starting from a non-equilibrium state of Staphylococcal nuclease (SNase) with C-terminal truncation (SNaseDelta). The results indicate a key role in the unfolding of the surface salt bridge between arginine 105 and glutamate 135. Experimentally, SNaseDelta is known to be partially unfolded. However, in simulations over 1 ns at 300 K and over 500 ps at 400 K, SNaseDelta remains stable in the native-like folded conformation, the salt bridge hindering unfolding. When the potential function is altered so as to selectively weaken the salt bridge, which then breaks rapidly at 430 K, the protein starts to unfold. The results suggest that breaking of this salt bridge presents a significant barrier to the unfolding transition of SNaseDelta from a native-like state to the unfolded state. Potential of mean force calculations indicate that the barrier height for this transition is approximately 7 kcal/mol.

Denaturing action of urea and guanidine hydrochloride towards two thermophilic esterases

The stability of two thermophilic esterases, AFEST from Archaeoglobus fulgidus and EST2 from Alicyclobacillus acidocaldarius, against the denaturing action of urea and guanidine hydrochloride has been investigated by means of steady-state fluorescence and circular dichroism measurements. Experimental results indicate that the two enzymes, even though very resistant to temperature and urea, show a resistance to guanidine hydrochloride weaker than expected on the basis of data collected so far for a large set of globular proteins. Structural information available for AFEST and EST2 and ideas that emerged from studies on the molecular origin of the greater thermal stability of thermophiles allow the suggestion of a reliable rationale. The present results may be an indication that the optimization of charge-charge interactions on the protein surface is a key factor for the stability of the two esterases.

Flexibility of the murine prion protein and its Asp178Asn mutant investigated by molecular dynamics simulations

Inherited forms of transmissible spongiform encephalopathy, e.g. familial Creutzfeldt-Jakob disease, Gerstmann-Sträussler-Scheinker syndrome and fatal familial insomnia, segregate with specific point mutations of the prion protein. It has been proposed that the pathologically relevant Asp178Asn (D178N) mutation might destabilize the structure of the prion protein because of the loss of the Arg164-Asp178 salt bridge. Molecular dynamics simulations of the structured C-terminal domain of the murine prion protein and the D178N mutant were performed to investigate this hypothesis. The D178N mutant did not deviate from the NMR conformation more than the wild type on the nanosecond time scale of the simulations. In agreement with CD spectroscopy experiments, no major structural rearrangement could be observed for the D178N mutant, apart from the N-terminal elongation of helix 2. The region of structure around the disulfide bridge deviated the least from the NMR conformation and showed the smallest fluctuations in all simulations in agreement with hydrogen exchange data of the wild type prion protein. Large deviations and flexibility were observed in the segments which are ill-defined in the NMR conformation. Moreover, helix 1 showed an increased degree of mobility, especially at its N-terminal region. The dynamic behavior of the D178N mutant and its minor deviation from the folded conformation suggest that the salt bridge between Arg164 and Asp178 might not be crucial for the stability of the prion protein.

Calculation of protein ionization equilibria with conformational sampling: pK(a) of a model leucine zipper, GCN4 and barnase

The use of conformational ensembles provided by nuclear magnetic resonance (NMR) experiments or generated by molecular dynamics (MD) simulations has been regarded as a useful approach to account for protein motions in the context of pK(a) calculations, yet the idea has been tested occasionally. This is the first report of systematic comparison of pK(a) estimates computed from long multiple MD simulations and NMR ensembles. As model systems, a synthetic leucine zipper, the naturally occurring coiled coil GCN4, and barnase were used. A variety of conformational averaging and titration curve-averaging techniques, or combination thereof, was adopted and/or modified to investigate the effect of extensive global conformational sampling on the accuracy of pK(a) calculations. Clustering of coordinates is proposed as an approach to reduce the vast diversity of MD ensembles to a few structures representative of the average electrostatic properties of the system in solution. Remarkable improvement of the accuracy of pK(a) predictions was achieved by the use of multiple MD simulations. By using multiple trajectories the absolute error in pK(a) predictions for the model leucine zipper was reduced to as low as approximately 0.25 pK(a) units. The validity, advantages, and limitations of explicit conformational sampling by MD, compared with the use of an average structure and a high internal protein dielectric value as means to improve the accuracy of pK(a) calculations, are discussed.

Selective inactivation of guanine-nucleotide-binding regulatory protein (G-protein) alpha and betagamma subunits by urea

G-protein-coupled receptors activate signal-transducing G-proteins, which consist of an alpha subunit and a betagamma dimer. Membrane extraction with 5-7 M urea has been used to uncouple receptors from endogenous G-proteins to permit reconstitution with purified G-proteins. We show that alpha(i) subunits are inactivated with 5 M urea whereas the betagamma dimer requires at least 7 M urea for its inactivation. There is no significant loss of receptors. Surprisingly, Western-blot analysis indicates that the urea-denatured alpha(i) subunit remains mostly membrane-bound and that beta is only partially removed. After 7 M urea treatment, both alpha(i1) and betagamma subunits are required to restore high-affinity agonist binding and receptor-catalysed guanosine 5'-[gamma-thio]triphosphate binding. We demonstrate the generality of this approach for four G(i)-coupled receptors (alpha(2A)-adrenergic, adenosine A1, 5-hydroxytryptamine(1A) and mu-opioid) expressed in insect cells and two mammalian cell lines. Thus a selectivity of urea for G-protein alpha versus betagamma subunits is established in both concentration and mechanism.