Valerie Daggett

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°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 proteins 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 ureas 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.

Protein Folding and Unfolding at Atomic Resolution

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.

The complete folding pathway of a protein from nanoseconds to microseconds

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 φ 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 α-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.

Is there a unifying mechanism for protein folding

Proteins appear to fold by diverse pathways, but variations of a simple mechanism - nucleation-condensation - describe the overall features of folding of most domains. In general, secondary structure is inherently unstable and its stability is enhanced by tertiary interactions. Consequently, an extensive interplay of secondary and tertiary interactions determines the transition-state for folding, which is structurally similar to the native state, being formed in a general collapse (condensation) around a diffuse nucleus. As the propensity for stable secondary structure increases, folding becomes more hierarchical and eventually follows a framework mechanism where the transition state is assembled from pre-formed secondary structural elements.

Potential energy function and parameters for simulations of the molecular dynamics of proteins and nucleic acids in solution

We present the complete set of energy parameters used in the ENCAD (Energy Calculation and Dynamics) simulation program [J. Mol. Biol. 168 (1983) 595]. Full details are given of the form of the potential, which has been designed for efficient simulation of trajectories of macromolecules in solution. Emphasis is placed on energy conservation and the nonbonded truncation schemes needed to achieve it. Simulations of macromolecules in solution with ENCAD are both very stable in that the native structure is preserved at room temperature and efficient in that nanosecond simulations take a few weeks on an ordinary workstation.

The present view of the mechanism of protein folding

We can track the positions and movements of all the atoms in small proteins as they fold and unfold by combining experimental studies with atomic-resolution molecular dynamics simulations. General principles as to how such complex architectures form so rapidly are now emerging from in-depth studies of a few proteins.

Increasing temperature accelerates protein unfolding without changing the pathway of unfolding.

We have traditionally relied on extremely elevated temperatures (498K, 225 degrees C) to investigate the unfolding process of proteins within the timescale available to molecular dynamics simulations with explicit solvent. However, recent advances in computer hardware have allowed us to extend our thermal denaturation studies to much lower temperatures. Here we describe the results of simulations of chymotrypsin inhibitor 2 at seven temperatures, ranging from 298K to 498K. The simulation lengths vary from 94ns to 20ns, for a total simulation time of 344ns, or 0.34 micros. At 298K, the protein is very stable over the full 50ns simulation. At 348K, corresponding to the experimentally observed melting temperature of CI2, the protein unfolds over the first 25ns, explores partially unfolded conformations for 20ns, and then refolds over the last 35ns. Above its melting temperature, complete thermal denaturation occurs in an activated process. Early unfolding is characterized by sliding or breathing motions in the protein core, leading to an unfolding transition state with a weakened core and some loss of secondary structure. After the unfolding transition, the core contacts are rapidly lost as the protein passes on to the fully denatured ensemble. While the overall character and order of events in the unfolding process are well conserved across temperatures, there are substantial differences in the timescales over which these events take place. We conclude that 498K simulations are suitable for elucidating the details of protein unfolding at a minimum of computational expense.

Can Non-Mechanical Proteins Withstand Force? Stretching Barnase by Atomic Force Microscopy and Molecular Dynamics Simulation

Atomic force microscopy (AFM) experiments have provided intriguing insights into the mechanical unfolding of proteins such as titin I27 from muscle, but will the same be possible for proteins that are not physiologically required to resist force? We report the results of AFM experiments on the forced unfolding of barnase in a chimeric construct with I27. Both modules are independently folded and stable in this construct and have the same thermodynamic and kinetic properties as the isolated proteins. I27 can be identified in the AFM traces based on its previous characterization, and distinct, irregular low-force peaks are observed for barnase. Molecular dynamics simulations of barnase unfolding also show that it unfolds at lower forces than proteins with mechanical function. The unfolding pathway involves the unraveling of the protein from the termini, with much more native-like secondary and tertiary structure being retained in the transition state than is observed in simulations of thermal unfolding or experimentally, using chemical denaturant. Our results suggest that proteins that are not selected for tensile strength may not resist force in the same way as those that are, and that proteins with similar unfolding rates in solution need not have comparable unfolding properties under force.

Unifying features in protein-folding mechanisms

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.

Selective Neuronal Targeting in Prion Disease

The pattern of scrapie prion protein (PrP(Sc)) accumulation in the brain is different for each prion strain. We tested whether the PrP(Sc) deposition pattern is influenced by the Asn-linked oligosaccharides of PrP(C) in transgenic mice. Deletion of the first oligosaccharide altered PrP(C) trafficking and prevented infection with two prion strains. Deletion of the second did not alter PrP(C) trafficking, permitted infection with one prion strain, and had a profound effect on the PrP(Sc) deposition pattern. Our data raise the possibility that glycosylation can modify the conformation of PrP(C). Glycosylation could affect the affinity of PrP(C) for a particular conformer of PrP(Sc), thereby determining the rate of nascent PrP(Sc) formation and the specific patterns of PrP(Sc) deposition.