Methods for computer simulation of proteins in solution developed and tested during the first funding period will be applied to systems under active experimental investigation. Our calculation protocol uses a periodic box, includes solvent molecules explicitly, allows full bond length and angle flexibility, truncates long-range interactions smoothly and conserves total energy. It gives well-behaved trajectories that can be computed efficiently. Our newly gained ability to run much longer, more realistic simulations, will help us form strong interactions with experimental collaborators at Standard, NIH, MIT and Los Alamos.
Specific aims are: 1. Can we simulate urea and trifluoroethanol? We will molecular dynamics to simulate these denaturing and renaturing solutions at different concentrations and temperatures and calibrate them against experimental data to obtain realistic solutions for use in unfolding and refolding simulations. 2. How do alpha-helix and Beta-hairpin structures unfold? How long does helix unfolding take? What is the nucleation center for helix and hairpin formation? How does Beta-hairpin unfolding, with its long-rang hydrogen bonds, differ from alpha-helix unfolding? These questions will be answered in unfolding simulations lasting 1 to 10 nanoseconds (ns) to be done at a range of temperatures. By reducing temperature slowly, we will also attempt to refold partially denatured alpha-helices and Beta- hairpins. 3. What are the earliest events in protein unfolding? We have carefully chosen a set of five proteins that are (a) small, (b) built according to different architectural principles, and (c) under experimental study. Short simulations (<100ps) will be run at a range of temperatures and densities. Water interactions, packing contacts and packing volume changes will be analyzed to reveal how structure is disrupted. The effects of concentrated urea and TFE will also be investigated. 4. How do unfolding pathways depend on protein architecture? Running prolonged (1-5 ns) simulations at elevated temperatures should caused each of the small proteins studied here to unfold. The nature of the intermediates in unfolding and the way they are formed will be analyzed and compared. 5. How mobile are intermediate states? Dynamics of intermediate states generated by unfolding these five proteins will be run at room temperature. Distances between cysteine, hydrophobic residues and charge pairs will be monitored. 6. Can we restore a native fold to a part of alpha-lactalbumin? We will use our packing methods and simulations to predict stabilizing mutations for a part known to be a molten globule at low temperature.

Agency
National Institute of Health (NIH)
Institute
National Institute of General Medical Sciences (NIGMS)
Type
Research Project (R01)
Project #
5R01GM041455-07
Application #
2180879
Study Section
Molecular and Cellular Biophysics Study Section (BBCA)
Project Start
1989-07-01
Project End
1998-06-30
Budget Start
1995-07-01
Budget End
1996-06-30
Support Year
7
Fiscal Year
1995
Total Cost
Indirect Cost
Name
Stanford University
Department
Anatomy/Cell Biology
Type
Schools of Medicine
DUNS #
800771545
City
Stanford
State
CA
Country
United States
Zip Code
94305
Sim, Adelene Y L; Minary, Peter; Levitt, Michael (2012) Modeling nucleic acids. Curr Opin Struct Biol 22:273-8
Moreno-Hernández, Sergio; Levitt, Michael (2012) Comparative modeling and protein-like features of hydrophobic-polar models on a two-dimensional lattice. Proteins 80:1683-93
Sim, Adelene Y L; Levitt, Michael; Minary, Peter (2012) Modeling and design by hierarchical natural moves. Proc Natl Acad Sci U S A 109:2890-5
Sim, Adelene Y L; Levitt, Michael (2011) Clustering to identify RNA conformations constrained by secondary structure. Proc Natl Acad Sci U S A 108:3590-5
Bernauer, Julie; Huang, Xuhui; Sim, Adelene Y L et al. (2011) Fully differentiable coarse-grained and all-atom knowledge-based potentials for RNA structure evaluation. RNA 17:1066-75
Samson, Abraham O; Levitt, Michael (2011) Normal modes of prion proteins: from native to infectious particle. Biochemistry 50:2243-8
Bray, Jenelle K; Weiss, Dahlia R; Levitt, Michael (2011) Optimized torsion-angle normal modes reproduce conformational changes more accurately than cartesian modes. Biophys J 101:2966-9
Chopra, Gaurav; Levitt, Michael (2011) Remarkable patterns of surface water ordering around polarized buckminsterfullerene. Proc Natl Acad Sci U S A 108:14455-60
Huang, Xuhui; Wang, Dong; Weiss, Dahlia R et al. (2010) RNA polymerase II trigger loop residues stabilize and position the incoming nucleotide triphosphate in transcription. Proc Natl Acad Sci U S A 107:15745-50
Minary, Peter; Levitt, Michael (2010) Conformational optimization with natural degrees of freedom: a novel stochastic chain closure algorithm. J Comput Biol 17:993-1010

Showing the most recent 10 out of 46 publications