Molecular dynamics simulations of the protein chymotrypsin inhibitor 2 in 8 M urea at 60°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 Tm of 75°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 ␤-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.S mall organic molecules in aqueous solution can have profound effects on protein stability, structure, and function. The use of these solutions to stabilize or destabilize proteins, depending on the cosolvent, is commonplace. In fact, protein studies are conducted almost exclusively in complex solutions. Chemical denaturation, with an agent such as urea, is one of the primary ways to assess protein stability, the effects of mutations on stability, and protein unfolding (1). Despite its widespread use, the molecular basis for urea's ability to denature proteins remains unknown. Urea may exert its effect directly, by binding to the protein, or indirectly, by altering the solvent environment (2-20). Most versions of the direct interaction model posit that urea binds to, and stabilizes, the denatured state (D), thereby favoring unfolding. But this interpretation does not explain how the protein surmounts the kinetic barrier to unfolding. In this regard, urea could bind to the protein and compete with native interactions, thereby actively participating in the unfolding process. Alternatively, it has been proposed that urea acts indirectly by altering the solvent environment, thereby mitigating the hydrophobic effect and facilitating the exposure of residues in the hydrophobic core. It is also possible that the mechanism of urea-promoted unfolding depends on the urea concentration. Unfortunately, it seems unlikely that experimental approaches will provide the molecular details of how urea denatures proteins, so we are employing atomic-resolution molecular dynamics (MD) simulations to address this issue....

Trimethylamine n-oxide (TMAO) is a naturally occurring osmolyte that stabilizes proteins and offsets the destabilizing effects of urea. To investigate the molecular mechanism of these effects, we have studied the thermodynamics of interaction between TMAO and protein functional groups. The solubilities of a homologous series of cyclic dipeptides were measured by differential refractive index and the dissolution heats were determined calorimetrically as a function of TMAO concentration at 25 degrees C. The transfer free energy of the amide unit (-CONH-) from water to 1 M TMAO is large and positive, indicating an unfavorable interaction between the TMAO solution and the amide unit. This unfavorable interaction is enthalpic in origin. The interaction between TMAO and apolar groups is slightly favorable. The transfer free energy of apolar groups from water to TMAO consists of favorable enthalpic and unfavorable entropic contributions. This is in contrast to the contributions for the interaction between urea and apolar groups. Molecular dynamics simulations were performed to provide a structural framework for the interpretation of these results. The simulations show enhancement of water structure by TMAO in the form of a slight increase in the number of hydrogen bonds per water molecule, stronger water hydrogen bonds, and long-range spatial ordering of the solvent. These findings suggest that TMAO stabilizes proteins via enhancement of water structure, such that interactions with the amide unit are discouraged.

The objective of this work is to obtain a water model for simulations of biological macromolecules in solution. A pragmatic approach is taken in which we use the same type of force field for the water as used for the solute and derive the water potential as an integral part of the ENCAD macromolecular potential. , Here we describe a flexible three-centered water model (F3C), which has already been used for many large-scale biological simulations, and compare it with other water models. The model is further tested by comparing calculated energetic, structural, and dynamic properties of liquid water, at several temperatures and pressures, with experiment. The F3C model is extremely simple and fits experimental data well for different temperatures, pressures, system sizes, and integration time steps. Because the F3C model works well with short-range truncation, it is well-suited to high-speed computation of long molecular dynamics trajectories of macromolecules in solution.

Combining experimental and simulation data to describe all of the structures and the pathways involved in folding a protein is problematical. Transition states can be mapped experimentally by phi values, but the denatured state is very difficult to analyse under conditions that favour folding. Also computer simulation at atomic resolution is currently limited to about a microsecond or less. Ultrafast-folding proteins fold and unfold on timescales accessible by both approaches, so here we study the folding pathway of the three-helix bundle protein Engrailed homeodomain. Experimentally, the protein collapses in a microsecond to give an intermediate with much native alpha-helical secondary structure, which is the major component of the denatured state under conditions that favour folding. A mutant protein shows this state to be compact and contain dynamic, native-like helices with unstructured side chains. In the transition state between this and the native state, the structure of the helices is nearly fully formed and their docking is in progress, approximating to a classical diffusion-collision model. Molecular dynamics simulations give rate constants and structural details highly consistent with experiment, thereby completing the description of folding at atomic resolution.

No abstract

Experiment and simulation are now conspiring to give atomic-level descriptions of protein folding relevant to folding, misfolding, trafficking, and degradation in the cell. We are on the threshold of predicting those protein folding events using simulation that has been carefully benchmarked by experiment.

No abstract

We compare the folding of representative members of a protein superfamily by experiment and simulation to investigate common features in folding mechanisms. The homeodomain superfamily of three-helical, single-domain proteins exhibits a spectrum of folding processes that spans the complete transition from concurrent secondary and tertiary structure formation (nucleation-condensation mechanism) to sequential secondary and tertiary formation (framework mechanism). The unifying factor in their mechanisms is that the transition state for (un)folding is expanded and very native-like, with the proportion and degree of formation of secondary and tertiary interactions varying. There is a transition, or slide, from the framework to nucleation-condensation mechanism with decreasing stability of the secondary structure. Thus, framework and nucleation-condensation are different manifestations of an underlying common mechanism.two-state ͉ three-state ͉ framework ͉ nucleation ͉ homeodomain A Holy Grail of protein folding is to find a single mechanism. Given the diversity of protein structure and the evolutionary pressure on function and not on folding rates, a unique mechanism for folding would seem unlikely. If there are simplifying features, then small, single-domain proteins may be the most likely to exhibit them. But such proteins seem to fold by two distinct mechanisms. The 6-85 repressor fragment (1) and the engrailed homeodomain (En-HD; ref.2) seem to fold by a classical diffusion-collision mechanism (3-5) whereby secondary structural elements form independently and then dock to form the tertiary structure. Chymotrypsin inhibitor 2, on the other hand, folds by nucleation-condensation, which is characterized by concerted consolidation of secondary and tertiary interactions as the whole domain collapses around an extended nucleus (6). It has been argued on general grounds that nucleation-condensation and diffusion-collision are different manifestations of a common mechanism in which secondary structure and tertiary structure form in parallel (7,8). Nucleationcondensation reflects the situation when secondary structure is inherently unstable in the absence of tertiary interactions whereas diffusion-collision becomes more probable with increasing stability of secondary structure.Studies of the folding of point mutants of a prototype protein are essential for discovering atomic level details of folding mechanisms and kinetics. Single-point mutants may even cause gross changes in the kinetics of folding, such as the transition from three-state to two-state folding (9). But, to extrapolate a general understanding of folding mechanisms, studies on members of the same fold family (different homologues sharing the same overall topology but with different primary structures) can be useful in finding correlations between amino acid sequences and three-dimensional structures (10-16). Although there can be different folding routes through different transition states for some proteins (17), it seems that mechanisms of folding are oft...