Challenges in generation of conformational ensembles for peptides and small proteins
|
|
- Myron Allison
- 5 years ago
- Views:
Transcription
1 Challenges in generation of conformational ensembles for peptides and small proteins Carlos Simmerling Stony Brook University What could (and does) go wrong? 1. Sampling: difficult to obtain converged ensembles, long-timescale dynamics 2. Force field errors 4 Functional form 4 Inaccurate parameters 4 Treatment of solvent 3. Validation against experimental data 4 Usually need to overcome (1) to improve (2) 1
2 Force field problems >1 million decoys AMBER ff94, ff99 Native not lowest energy = refit prehistoric using decoys? Okur, Strockbine, Hornak, Simmerling, J. Comp. Chem Trpzip2 indole stacking Our model 1HRX 1LE1 4us: the stacking of the outer Trp side chains differs from the published structures 41LE1: updated coordinates differ from the previous ones in the side chain orientations of Trp4 and Trp11 2
3 Model system: trpzip2 β-hairpin 1 Reference data: 150 ns standard REMD (GB, ffga) Thermal profiles independent of initial structure Good model: expt. T m = 345K 1 1 Cochran, Skelton & Starovasnik, PNAS 2001 Trp-cage tc5b mini-protein 4Small, stable and folds quickly 420 residues 498% folded at 280K 4Folding time of 4µsec 4Blind prediction using MD simulations and energy hydrophobic core N α helix i:i 5 C 3 10 helix PP ΙΙ i:i+10 NMR-based structure: Neidigh, Fesinmeyer, Andersen, NSB 2002 Folding rate: Qiu, Pabit, Roitberg, Hagen JACS 2002r Structure prediction: Simmerling, Strockbine, Roitberg, JACS 2002 MD low energy 1L2Y model 1 3
4 Comparison of Amber parameter sets 4 Gly bb is treated differently than others in Amber 4 Study capped Ala 3 and Gly 3 4 TIP3P explicit water 4 MD, REMD converged ensembles at 300K 4 Multiple existing Amber parameter sets 4 Compare to our new ff99sb based on backbone refitting with QM data (LMP2/cc-pVTZ(-f)) for multiple conformations of Ala 3 and Gly 3 Hornak, Abel, Okur, Strockbine, Roitberg and Simmerling, Proteins 2006 Gly 3 Ensembles in explicit water Ala 3, Gly 3 in explicit water Ala 3 Gly often very unreasonable Ala basins quite sensitive Hornak, Abel, Okur, Strockbine, Roitberg and Simmerling, Proteins
5 Gly: chemical shifts for simulation ensemble (SHIFTS, D. Case) Experiment: Random coil ~4ppm Folded state stereospecific HA2 ~ 0.96 HA3 ~ 3.43 Dynamics: agreement with NMR improved Hornak, Abel, Okur, Strockbine, Roitberg and Simmerling, Proteins qualitatively better than previous force fields 4 yields NMR spin relaxation parameters in near quantitative agreement with experimental values 4RDC accuracy that is comparable to or better than the best static structural models and the NMR ensemble Showalter and Brüschweiler, JCTC 2007, JACS
6 HP36: a model system for protein folding. Folded/unfolded state? 4 NMR and X-ray structure differ: 4 core packing, hydrogen bonds, length of the helices 4 Sequence (HP35 N68H vs HP36) 4 ph, temperature X-ray = Yellow NMR = Blue McKnight, Matsudaira,Kim Nat. Struct. Biol, 1996 Chiu, Kubelka, Herbst-Irmer, Eaton, Hofrichter, Davies, PNAS 2005 MD results (explicit water) 4Simulations compared to both experimental coordinate sets Xray ref NMR ref 4MD always moves closer to X-ray 4Validated using expt. double mutant cycles Blue = NMR Yellow = X-ray Green = Simulation WIckstrom, Bi, Hornak, Raleigh and Simmerling, Biochemistry
7 Is there local structure in HP36 unfolded state? HP-3 HP-2 HP-1 REMD in TIP3P Independent initial coordinates Modest native HP-1 population No helix in polyala of same length Trends match experiment Tang, Rigotti, Fairman,. & Raleigh, Biochemistry (2004) Wickstrom, Okur, Song, Hornak, Raleigh & Simmerling, JMB 2006 HP21 : see Lauren Wickstrom s poster 7
8 Decoy-based (GB) parameters ffga, 2002 Ala tetrapeptide QM parameters ff99mod2, 2002 Ala and Gly tetrapeptide QM + tetrapeptide REMD in TIP3P + decoy screening: ff99sb, 2005 All retain ff94 RESP charge model (unlike ff03) Does GB give good results? 8
9 Influence of GB on secondary structure Capped Ala 10 (no salt bridges, etc) REMD: TIP3P, GB HCT, GB OBC, GB neck G solv : PB, GB HCT, GB OBC, GB neck, TIP3P (TI) Roe and Simmerling, J. Phys. Chem. B, 2006 Ala 10 ensembles: GB vs. TIP3P Error bars reflect the difference between 2 REMD Simulations of Ala10 from different initial conformations. Different solvent models give different secondary structure populations. Roe and Simmerling, J. Phys. Chem. B,
10 Comparing Differences in G Pol, G Pol PP2 Collapsed GPol RMSD from TIP3P PE GBHCT Collapsed GBOBC GBNeck Overall PP2 Collapsed Collapsed Analogous to folding GB is qualitative (at least ones in Amber) Ion pairing, secondary structure peptides are particularly sensitive Explicit water is costly Viscosity Dependence of # REMD replicas on system size Partial REMD 1 Hybrid solvent REMD 2 1 Cheng, Cui, Hornak & Simmerling, J. Phys. Chem. B, Okur, Wickstrom, Layten, Geney, Song, Hornak & Simmerling, J. Chem. Theory Comp.,
11 Temperature Replica Exchange 4MD runs over range of T 4Periodically swap structures 4Faster convergence 4Populations as a function of T 4Drawback: E must be small small T for large systems REMD Hansmann, U., CPL 1997 Sugita & Okamoto, CPL 1999 W 300K 325K 350K 375K 4 Impose reversibility/detailed balance I ( X ) w( X X ) = W( X ') w( X ' X ) 4Impose limiting distribution on exchange calculation W = P exch ( X ) = e( β ) E X ( β β )( E E ) n m ( ( ) ) = min 1,exp i j REMD is more efficient than MD REMD MD 4Levy: systems where temperature helps folding rate! 11
12 MD vs. REMD for β-sheet, GB solvent MD N Replica exchange Roe, Hornak and Simmerling, J. Mol. Biol Improving REMD in explicit water Replica requirement Convergence rate 12
13 Better scaling using a hybrid solvent model Perform simulations in full periodic box of explicit solvent Retain only first shell of water + reaction field in exchange calculation Smaller perceived system size permits fewer replicas Avoids problems associated with hybrid solvent MD Approximate since Hamiltonian during exchange not same as during MD Okur, Wickstrom, Layten, Geney, Song, Hornak & Simmerling, J. Chem. Theory Comp., 2006 Hybrid solvent REMD: polyala Ala 10 Explicit Solvent GB OBC Hybrid GB OBC + 1 st shell α 24.9 ± ± ± 1.6 β 19.5 ± ± ± 1.6 P II 39.5 ± ± ± 0.5 α L 8.4 ± ± ± 0.3 SASA ± ± ± 2.5 Alanine dipeptide insensitive to solvent model Alanine tetrapeptide and Ala 10 reveal GB weaknesses Hybrid explicit+gb is in agreement with TIP3P: PPII Okur, Wickstrom, Layten, Geney, Song, Hornak & Simmerling, J. Chem. Theory Comp.,
14 Salt bridges: Ace-Arg-Ala-Ala-Glu-NH 2 TIP3P standard REMD : 46 replicas, 296K to 584K (15%) GB standard REMD: 6 replicas, 300K to 636K Hybrid solvent REMD: 8 replicas, 280K to 570K, 75 water molecules in hybrid shell Okur, Wickstrom & Simmerling, JCTC in press GB Salt bridge PMFs TIP3P, hybrid GB Poor correlation between backbone conformation in GB/TIP3P Improved correlation with hybrid GB salt bridge 2-3 kcal too strong Improved PMF with hybrid Restrained backbone Hybrid reproduces solvent separated mimimum Geney, Layten, Gomperts, Hornak and Simmerling, JCTC 2006 Okur, Wickstrom & Simmerling, JCTC in press 14
15 HP-1 with hybrid solvent REMD MLSDEDFKAVFGM Helical content Hybrid solvent model again significantly improves agreement with standard REMD in explicit solvent Some residual GB-like behavior in hybrid model at low T Use PB, include dispersion term Okur, Wickstrom & Simmerling, JCTC in press Hybrid solvent allows explicit water at nearly the same # replicas as GB. What can we do about REMD convergence rates? 15
16 300K REMD with a structure reservoir 325K 350K 375K 400K J-walk into reservoir Frantz, Freeman & Doll, JCP steps: (1) generate reservoir, (2) run REMD 4uncouple slow, high T sampling from many-replica reweighting Zuckerman PRL 2006, Yang JCP 2006, Simmerling, JCTC 2007 R-REMD results: trpzip2 in GB R-REMD: 10,000 structures and velocities from 400K MD Resulting data are in good agreement with standard REMD 350K ensembles Okur, Wickstrom, Simmerling, JCTC 2007 R 2 > 0.99 R 2 = 0.96 without largest cluster 16
17 Accuracy good, what about efficiency? dpdp 3-stranded β-sheet 1 6 replicas up to 400K, 50ns (300ns) 260ns MD for reservoir Melting profile matches std REMD 2 R-REMD converges much faster even when we include reservoir generation 1 Schenck & Gellman, JACS Roe, Hornak & Simmerling, JMB 2005 (ps) 17
18 R-REMD in explicit water Ala 10 HP-1 R-REMD vs. standard REMD Wickstrom and Simmerling, in prep Reservoir must be Boltzmann weighted (detailed balance) Very difficult to generate, especially in explicit water Reservoir is at high T: faster equilibration, best aspect of REMD Is this ensemble really a requirement? 18
19 W REMD with non-boltzmann reservoir W 300K 325K 350K 375K reservoir 4Impose reversibility/detailed balance I ( X ) w( X X ) = W( X ') w( X ' X ) ( X ) = e( β ) E X = β ( E ) i E j W ( X ) 1 ' = N 4Exchange with reservoir employs new equation 4All other pairs use standard exchange criterion ( β) Roitberg, Okur and Simmerling, JPC B, 2007 Trpzip2 REMD with non-boltzmann reservoir 669 structures in flat cluster reservoir (only 1 native) Correct thermal profile is obtained only with new exchange equation Roitberg, Okur and Simmerling, JPC B,
20 P Using a non-boltzmann reservoir with user-specified probabilities exchange j β ( E ) ( ) i E q, T; q, R q, T; q, R e j i Desirable to expand reservoir (or even build it during the REMD run) j j i N N i ( ) Requires assignment of reservoir weight to replica MD structures Use dihedral angle grid-based weights D. Zuckerman, black box reweighting probability density in dihedral space cluster volume is better defined than RMSD-based approach Roe and Simmerling, in prep Trpzip2/GB with weighted non-boltzmann reservoir Standard REMD Grid-weighted reservoir Standard REMD Grid-weighted reservoir Use of reservoir speeds convergence Populations in reservoir can be arbitrary and dynamic 20
21 Convergence comparison at 300K trpzip2 grid N-B R-REMD in TIP3P 60ns 360ns Standard REMD nearly impossible to converge 4 Needs >30 folding events ~10ns grid-weighted non- Boltzmann R-REMD J-walk into REMD history, in addition to high T sampling 21
22 Goal: thermodynamic ensembles in explicit solvent at relevant temperatures 1. Standard REMD: converge large # of T 2. Hybrid solvent REMD: converge smaller # of T 3. REMD with converge at 1 high T, Boltzmann reservoir get data for lower T 4. REMD with no need to sample non-boltzmann reservoir converged ensemble in MD? 5. REMD using variables other than temperature Applying what we have learned to larger systems - Which protein parameters are best? - How well do various water models perform? - And so on 22
23 HIV-1 protease and peptide substrate Model for dynamic behavior of HIV PR closed bound enzyme semi-open unbound enzyme Hornak, Okur, Rizzo and Simmerling, PNAS 2006 Layten, Hornak and Simmerling, JACS 2006 Hornak, Okur, Rizzo and Simmerling JACS 2006 Hornak and Simmerling, Drug Discovery Today,
24 Summary Force fields and solvent models continue to evolve. Both can introduce biases. Some should be retired. Sampling and accuracy are coupled, especially solvent Standard REMD converges slowly, requires many replicas Replica requirement can be reduced Sampling problem can be decoupled from REMD Structure diversity from other methods can speed convergence of MD-based sampling 4Asim Okur, Dan Roe, Lauren Wickstrom, Ding Fangyu 4Melinda Layten, Salma Rafi, Kun Song, Catherine Kelso, AJ Campbell, Christina Bergonzo Acknowledgements 4NIH, NCSA 4Adrian Roitberg (U Florida) Postdoc positions available! 24
25 Table 4. GB Effective Radii Average RMSD from Perfect (PE) Radii (Ǻ) A) All GBHCT GBOBC GBNeck F) C GBHCT GBOBC GBNeck alpha 0.25 ± ± ± 0.02 alpha 0.16 ± ± ± 0.03 hairpin 0.18 ± ± ± 0.01 hairpin 0.08 ± ± ± 0.02 left 0.20 ± ± ± 0.03 left 0.12 ± ± ± 0.04 pp ± ± ± 0.00 pp ± ± ± 0.01 B) BB GBHCT GBOBC GBNeck G) CA GBHCT GBOBC GBNeck alpha 0.35 ± ± ± 0.03 alpha 0.05 ± ± ± 0.02 hairpin 0.20 ± ± ± 0.01 hairpin 0.09 ± ± ± 0.01 left 0.27 ± ± ± 0.04 left 0.07 ± ± ± 0.03 pp ± ± ± 0.00 pp ± ± ± 0.01 C) H GBHCT GBOBC GBNeck H) CB GBHCT GBOBC GBNeck alpha 0.71 ± ± ± 0.06 alpha 0.03 ± ± ± 0.01 hairpin 0.39 ± ± ± 0.03 hairpin 0.04 ± ± ± 0.00 left 0.50 ± ± ± 0.06 left 0.01 ± ± ± 0.00 pp ± ± ± 0.00 pp ± ± ± 0.00 D) O GBHCT GBOBC GBNeck I) HA GBHCT GBOBC GBNeck alpha 0.16 ± ± ± 0.02 alpha 0.07 ± ± ± 0.01 hairpin 0.16 ± ± ± 0.02 hairpin 0.34 ± ± ± 0.03 left 0.18 ± ± ± 0.02 left 0.10 ± ± ± 0.02 pp ± ± ± 0.00 pp ± ± ± 0.00 E) N GBHCT GBOBC GBNeck alpha 0.26 ± ± ± 0.04 hairpin 0.16 ± ± ± 0.02 left 0.27 ± ± ± 0.05 pp ± ± ± 0.01 J) Overall Averages GBHCT GBOBC GBNeck GBHCT GBOBC GBNeck All C BB CA H CB O HA N for Peter 25
Limitations of temperature replica exchange (T-REMD) for protein folding simulations
Limitations of temperature replica exchange (T-REMD) for protein folding simulations Jed W. Pitera, William C. Swope IBM Research pitera@us.ibm.com Anomalies in protein folding kinetic thermodynamic 322K
More informationComputer simulations of protein folding with a small number of distance restraints
Vol. 49 No. 3/2002 683 692 QUARTERLY Computer simulations of protein folding with a small number of distance restraints Andrzej Sikorski 1, Andrzej Kolinski 1,2 and Jeffrey Skolnick 2 1 Department of Chemistry,
More informationSequential resonance assignments in (small) proteins: homonuclear method 2º structure determination
Lecture 9 M230 Feigon Sequential resonance assignments in (small) proteins: homonuclear method 2º structure determination Reading resources v Roberts NMR of Macromolecules, Chap 4 by Christina Redfield
More informationWhat makes a good graphene-binding peptide? Adsorption of amino acids and peptides at aqueous graphene interfaces: Electronic Supplementary
Electronic Supplementary Material (ESI) for Journal of Materials Chemistry B. This journal is The Royal Society of Chemistry 21 What makes a good graphene-binding peptide? Adsorption of amino acids and
More informationCan a continuum solvent model reproduce the free energy landscape of a β-hairpin folding in water?
Can a continuum solvent model reproduce the free energy landscape of a β-hairpin folding in water? Ruhong Zhou 1 and Bruce J. Berne 2 1 IBM Thomas J. Watson Research Center; and 2 Department of Chemistry,
More informationReplica Exchange with Solute Scaling: A More Efficient Version of Replica Exchange with Solute Tempering (REST2)
pubs.acs.org/jpcb Replica Exchange with Solute Scaling: A More Efficient Version of Replica Exchange with Solute Tempering (REST2) Lingle Wang, Richard A. Friesner, and B. J. Berne* Department of Chemistry,
More informationConvergence of replica exchange molecular dynamics
THE JOURNAL OF CHEMICAL PHYSICS 123, 154105 2005 Convergence of replica exchange molecular dynamics Wei Zhang and Chun Wu Department of Chemistry and Biochemistry, University of Delaware, Newark, Delaware
More informationFree Energy Landscape of Protein Folding in Water: Explicit vs. Implicit Solvent
PROTEINS: Structure, Function, and Genetics 53:148 161 (2003) Free Energy Landscape of Protein Folding in Water: Explicit vs. Implicit Solvent Ruhong Zhou* IBM T.J. Watson Research Center, Yorktown Heights,
More informationGeneralized Born Model with a Simple, Robust Molecular Volume Correction
156 J. Chem. Theory Comput. 007, 3, 156-169 Generalized Born Model with a Simple, Robust Molecular Volume Correction John Mongan,,, Carlos Simmerling, # J. Andrew McCammon,,,,,,3 David A. Case,,f and Alexey
More informationLecture 12: Solvation Models: Molecular Mechanics Modeling of Hydration Effects
Statistical Thermodynamics Lecture 12: Solvation Models: Molecular Mechanics Modeling of Hydration Effects Dr. Ronald M. Levy ronlevy@temple.edu Bare Molecular Mechanics Atomistic Force Fields: torsion
More informationPolypeptide Folding Using Monte Carlo Sampling, Concerted Rotation, and Continuum Solvation
Polypeptide Folding Using Monte Carlo Sampling, Concerted Rotation, and Continuum Solvation Jakob P. Ulmschneider and William L. Jorgensen J.A.C.S. 2004, 126, 1849-1857 Presented by Laura L. Thomas and
More informationExploring the Changes in the Structure of α-helical Peptides Adsorbed onto Carbon and Boron Nitride based Nanomaterials
Exploring the Changes in the Structure of α-helical Peptides Adsorbed onto Carbon and Boron Nitride based Nanomaterials Dr. V. Subramanian Chemical Laboratory, IPC Division CSIR-Central Leather Research
More informationProtein Structures: Experiments and Modeling. Patrice Koehl
Protein Structures: Experiments and Modeling Patrice Koehl Structural Bioinformatics: Proteins Proteins: Sources of Structure Information Proteins: Homology Modeling Proteins: Ab initio prediction Proteins:
More informationTimescales of Protein Dynamics
Timescales of Protein Dynamics From Henzler-Wildman and Kern, Nature 2007 Summary of 1D Experiment time domain data Fourier Transform (FT) frequency domain data or Transverse Relaxation Ensemble of Nuclear
More informationSupplementary Information
Supplementary Information Resveratrol Serves as a Protein-Substrate Interaction Stabilizer in Human SIRT1 Activation Xuben Hou,, David Rooklin, Hao Fang *,,, Yingkai Zhang Department of Medicinal Chemistry
More informationTimescales of Protein Dynamics
Timescales of Protein Dynamics From Henzler-Wildman and Kern, Nature 2007 Dynamics from NMR Show spies Amide Nitrogen Spies Report On Conformational Dynamics Amide Hydrogen Transverse Relaxation Ensemble
More informationRoutine access to millisecond timescale events with accelerated molecular dynamics
Routine access to millisecond timescale events with accelerated molecular dynamics Levi C.T. Pierce, Romelia Salomon-Ferrer, Cesar Augusto F. de Oliveira #, J. Andrew McCammon #, Ross C. Walker * SUPPORTING
More informationProtein Folding Prof. Eugene Shakhnovich
Protein Folding Eugene Shakhnovich Department of Chemistry and Chemical Biology Harvard University 1 Proteins are folded on various scales As of now we know hundreds of thousands of sequences (Swissprot)
More informationMolecular Mechanics. I. Quantum mechanical treatment of molecular systems
Molecular Mechanics I. Quantum mechanical treatment of molecular systems The first principle approach for describing the properties of molecules, including proteins, involves quantum mechanics. For example,
More informationIntroduction to" Protein Structure
Introduction to" Protein Structure Function, evolution & experimental methods Thomas Blicher, Center for Biological Sequence Analysis Learning Objectives Outline the basic levels of protein structure.
More informationPhysiochemical Properties of Residues
Physiochemical Properties of Residues Various Sources C N Cα R Slide 1 Conformational Propensities Conformational Propensity is the frequency in which a residue adopts a given conformation (in a polypeptide)
More informationPacking of Secondary Structures
7.88 Lecture Notes - 4 7.24/7.88J/5.48J The Protein Folding and Human Disease Professor Gossard Retrieving, Viewing Protein Structures from the Protein Data Base Helix helix packing Packing of Secondary
More informationReproducible Protein Folding with the Stochastic Tunneling Method
Reproducible Protein Folding with the Stochastic Tunneling Method A. Schug, 1 T. Herges, 1 and W. Wenzel 1 1 Forschungszentrum Karlsruhe, Institut für Nanotechnologie, 76021 Karlsruhe, Germany (Dated:
More informationProtein folding. Today s Outline
Protein folding Today s Outline Review of previous sessions Thermodynamics of folding and unfolding Determinants of folding Techniques for measuring folding The folding process The folding problem: Prediction
More informationProtein Dynamics. The space-filling structures of myoglobin and hemoglobin show that there are no pathways for O 2 to reach the heme iron.
Protein Dynamics The space-filling structures of myoglobin and hemoglobin show that there are no pathways for O 2 to reach the heme iron. Below is myoglobin hydrated with 350 water molecules. Only a small
More informationExchange frequency in replica exchange molecular dynamics
THE JOURNAL OF CHEMICAL PHYSICS 128, 024103 2008 Exchange frequency in replica exchange molecular dynamics Daniel Sindhikara Department of Physics, University of Florida, Gainesville, Florida 32611-8435,
More informationAre the native states of proteins kinetic traps? Abstract
Are the native states of proteins kinetic traps? Leonor Cruzeiro CCMAR and FCT, Universidade do Algarve, 85-139 Faro, Portugal Paulo A. Lopes CITI, Departamento de Informática, Faculdade de Ciências e
More informationStructural and mechanistic insight into the substrate. binding from the conformational dynamics in apo. and substrate-bound DapE enzyme
Electronic Supplementary Material (ESI) for Physical Chemistry Chemical Physics. This journal is the Owner Societies 215 Structural and mechanistic insight into the substrate binding from the conformational
More informationMolecular dynamics simulations of anti-aggregation effect of ibuprofen. Wenling E. Chang, Takako Takeda, E. Prabhu Raman, and Dmitri Klimov
Biophysical Journal, Volume 98 Supporting Material Molecular dynamics simulations of anti-aggregation effect of ibuprofen Wenling E. Chang, Takako Takeda, E. Prabhu Raman, and Dmitri Klimov Supplemental
More informationMulti-scale approaches in description and design of enzymes
Multi-scale approaches in description and design of enzymes Anastassia Alexandrova and Manuel Sparta UCLA & CNSI Catalysis: it is all about the barrier The inside-out protocol: Big Aim: development of
More informationProblem Set 1
2006 7.012 Problem Set 1 Due before 5 PM on FRIDAY, September 15, 2006. Turn answers in to the box outside of 68-120. PLEASE WRITE YOUR ANSWERS ON THIS PRINTOUT. 1. For each of the following parts, pick
More informationThe Dominant Interaction Between Peptide and Urea is Electrostatic in Nature: A Molecular Dynamics Simulation Study
Dror Tobi 1 Ron Elber 1,2 Devarajan Thirumalai 3 1 Department of Biological Chemistry, The Hebrew University, Jerusalem 91904, Israel 2 Department of Computer Science, Cornell University, Ithaca, NY 14853
More informationDominant Paths in Protein Folding
Dominant Paths in Protein Folding Henri Orland SPhT, CEA-Saclay France work in collaboration with P. Faccioli, F. Pederiva, M. Sega University of Trento Henri Orland Annecy meeting 2006 Outline Basic notions
More informationMonte Carlo simulations of polyalanine using a reduced model and statistics-based interaction potentials
THE JOURNAL OF CHEMICAL PHYSICS 122, 024904 2005 Monte Carlo simulations of polyalanine using a reduced model and statistics-based interaction potentials Alan E. van Giessen and John E. Straub Department
More informationSupplementary Figures:
Supplementary Figures: Supplementary Figure 1: The two strings converge to two qualitatively different pathways. A) Models of active (red) and inactive (blue) states used as end points for the string calculations
More informationEnergy landscapes of model polyalanines
JOURNAL OF CHEMICAL PHYSICS VOLUME 117, NUMBER 3 15 JULY 2002 Energy landscapes of model polyalanines Paul N. Mortenson, David A. Evans, and David J. Wales University Chemical Laboratories, Lensfield Road,
More informationCharacterizing Structural Transitions of Membrane Transport Proteins at Atomic Detail Mahmoud Moradi
Characterizing Structural Transitions of Membrane Transport Proteins at Atomic Detail Mahmoud Moradi NCSA Blue Waters Symposium for Petascale Science and Beyond Sunriver, Oregon May 11, 2015 Outline Introduction
More informationTrpzip-based beta hairpin temperature jump IR studies enhanced by sitespecific
Trpzip-based beta hairpin temperature jump IR studies enhanced by sitespecific isotope labeling Carsten Kretjschi 1, Karin auser 1, Rong uang 2, Tim Keiderling 2 1 Institute of Biophysics, University of
More informationMagnetic Resonance Lectures for Chem 341 James Aramini, PhD. CABM 014A
Magnetic Resonance Lectures for Chem 341 James Aramini, PhD. CABM 014A jma@cabm.rutgers.edu " J.A. 12/11/13 Dec. 4 Dec. 9 Dec. 11" " Outline" " 1. Introduction / Spectroscopy Overview 2. NMR Spectroscopy
More informationAdvanced sampling. fluids of strongly orientation-dependent interactions (e.g., dipoles, hydrogen bonds)
Advanced sampling ChE210D Today's lecture: methods for facilitating equilibration and sampling in complex, frustrated, or slow-evolving systems Difficult-to-simulate systems Practically speaking, one is
More informationCommon Structural Transitions in Explicit-Solvent Simulations of Villin Headpiece Folding
2338 Biophysical Journal Volume 97 October 2009 2338 2347 Common Structural Transitions in Explicit-Solvent Simulations of Villin Headpiece Folding Peter L. Freddolino and Klaus Schulten* Center for Biophysics
More informationSignificantly Improved Protein Folding Thermodynamics Using a Dispersion-Corrected Water Model and a New Residue-Specific Force Field
Significantly Improved Protein Folding Thermodynamics Using a Dispersion-Corrected Water Model and a New Residue-Specific Force Field Hao-Nan Wu, Fan Jiang, *, and Yun-Dong Wu, Laboratory of Computational
More informationTHE TANGO ALGORITHM: SECONDARY STRUCTURE PROPENSITIES, STATISTICAL MECHANICS APPROXIMATION
THE TANGO ALGORITHM: SECONDARY STRUCTURE PROPENSITIES, STATISTICAL MECHANICS APPROXIMATION AND CALIBRATION Calculation of turn and beta intrinsic propensities. A statistical analysis of a protein structure
More informationMonte Carlo simulation of proteins through a random walk in energy space
JOURNAL OF CHEMICAL PHYSICS VOLUME 116, NUMBER 16 22 APRIL 2002 Monte Carlo simulation of proteins through a random walk in energy space Nitin Rathore and Juan J. de Pablo a) Department of Chemical Engineering,
More informationβ m] E ww (X) (4) Replica Exchange with Solute Tempering: Efficiency in Large Scale Systems
J. Phys. Chem. B 2007, 111, 5405-5410 5405 Replica Exchange with Solute Tempering: Efficiency in Large Scale Systems Xuhui Huang,, Morten Hagen, Byungchan Kim, Richard A. Friesner, Ruhong Zhou,, and B.
More informationBMB/Bi/Ch 173 Winter 2018
BMB/Bi/Ch 173 Winter 2018 Homework Set 8.1 (100 Points) Assigned 2-27-18, due 3-6-18 by 10:30 a.m. TA: Rachael Kuintzle. Office hours: SFL 220, Friday 3/2 4:00-5:00pm and SFL 229, Monday 3/5 4:00-5:30pm.
More informationI690/B680 Structural Bioinformatics Spring Protein Structure Determination by NMR Spectroscopy
I690/B680 Structural Bioinformatics Spring 2006 Protein Structure Determination by NMR Spectroscopy Suggested Reading (1) Van Holde, Johnson, Ho. Principles of Physical Biochemistry, 2 nd Ed., Prentice
More informationEnergy Minimization of Protein Tertiary Structure by Parallel Simulated Annealing using Genetic Crossover
Minimization of Protein Tertiary Structure by Parallel Simulated Annealing using Genetic Crossover Tomoyuki Hiroyasu, Mitsunori Miki, Shinya Ogura, Keiko Aoi, Takeshi Yoshida, Yuko Okamoto Jack Dongarra
More informationEquilibrium Structure and Folding of a Helix-Forming Peptide: Circular Dichroism Measurements and Replica-Exchange Molecular Dynamics Simulations
3786 Biophysical Journal Volume 87 December 2004 3786 3798 Equilibrium Structure and Folding of a Helix-Forming Peptide: Circular Dichroism Measurements and Replica-Exchange Molecular Dynamics Simulations
More informationBiology Chemistry & Physics of Biomolecules. Examination #1. Proteins Module. September 29, Answer Key
Biology 5357 Chemistry & Physics of Biomolecules Examination #1 Proteins Module September 29, 2017 Answer Key Question 1 (A) (5 points) Structure (b) is more common, as it contains the shorter connection
More informationCONSTANT ph REPLICA EXCHANGE MOLECULAR DYNAMICS STUDY OF PROTEIN STRUCTURE AND DYNAMICS
CONSTANT ph REPLICA EXCHANGE MOLECULAR DYNAMICS STUDY OF PROTEIN STRUCTURE AND DYNAMICS By YILIN MENG A DISSERTATION PRESENTED TO THE GRADUATE SCHOOL OF THE UNIVERSITY OF FLORIDA IN PARTIAL FULFILLMENT
More informationQM/MM study on inhibitor of HIV-1 protease. Graduate School of Science, Kyoto University Masahiko Taguchi, Masahiko Kaneso, Shigehiko Hayashi
QM/MM study on inhibitor of HIV-1 protease Graduate School of Science, Kyoto University Masahiko Taguchi, Masahiko Kaneso, Shigehiko Hayashi HIV Human Immunodeficiency Virus Origin of AIDS HIV appearing
More informationExploring the Free Energy Surface of Short Peptides by Using Metadynamics
John von Neumann Institute for Computing Exploring the Free Energy Surface of Short Peptides by Using Metadynamics C. Camilloni, A. De Simone published in From Computational Biophysics to Systems Biology
More informationFree energy calculations and the potential of mean force
Free energy calculations and the potential of mean force IMA Workshop on Classical and Quantum Approaches in Molecular Modeling Mark Tuckerman Dept. of Chemistry and Courant Institute of Mathematical Science
More informationDetails of Protein Structure
Details of Protein Structure Function, evolution & experimental methods Thomas Blicher, Center for Biological Sequence Analysis Anne Mølgaard, Kemisk Institut, Københavns Universitet Learning Objectives
More informationBiomolecules: lecture 10
Biomolecules: lecture 10 - understanding in detail how protein 3D structures form - realize that protein molecules are not static wire models but instead dynamic, where in principle every atom moves (yet
More informationNMR study of complexes between low molecular mass inhibitors and the West Nile virus NS2B-NS3 protease
University of Wollongong Research Online Faculty of Science - Papers (Archive) Faculty of Science, Medicine and Health 2009 NMR study of complexes between low molecular mass inhibitors and the West Nile
More informationDe novo all atom folding of helical proteins
John von Neumann Institute for Computing De novo all atom folding of helical proteins A. Verma, S. Murthy, K. H. Lee, E. Starikov, Wolfgang Wenzel published in NIC Workshop 2006, From Computational Biophysics
More informationMany proteins spontaneously refold into native form in vitro with high fidelity and high speed.
Macromolecular Processes 20. Protein Folding Composed of 50 500 amino acids linked in 1D sequence by the polypeptide backbone The amino acid physical and chemical properties of the 20 amino acids dictate
More informationVibrational Spectral Simulation for Peptides of Mixed Secondary Structure: Method Comparisons with the Trpzip Model Hairpin
J. Phys. Chem. B 2005, 109, 23687-23697 23687 Vibrational Spectral Simulation for Peptides of Mixed Secondary Structure: Method Comparisons with the Trpzip Model Hairpin Petr Bouř*, and Timothy A. Keiderling*,
More informationSampling the conformation space of complex systems, such as
Replica exchange with solute tempering: A method for sampling biological systems in explicit water Pu Liu*, Byungchan Kim*, Richard A. Friesner, and B. J. Berne Department of Chemistry and Center for Biomolecular
More informationUnderstanding the Folding and Stability of a Zinc Finger-Based Full Sequence Design Protein With Replica Exchange Molecular Dynamics Simulations
67:338 349 (2007) Understanding the Folding and Stability of a Zinc Finger-Based Full Sequence Design Protein With Replica Exchange Molecular Dynamics Simulations Wenfei Li, 1 Jian Zhang, 1 and Wei Wang
More informationProtein Structure Refinement Using 13 C α Chemical. Shift Tensors. Benjamin J. Wylie, Charles D. Schwieters, Eric Oldfield and Chad M.
Protein Structure Refinement Using 13 C α Chemical Shift Tensors Benjamin J. Wylie, Charles D. Schwieters, Eric Oldfield and Chad M. Rienstra * Department of Chemistry, University of Illinois at Urbana-Champaign,
More informationIntroduction to Comparative Protein Modeling. Chapter 4 Part I
Introduction to Comparative Protein Modeling Chapter 4 Part I 1 Information on Proteins Each modeling study depends on the quality of the known experimental data. Basis of the model Search in the literature
More informationSimulating Folding of Helical Proteins with Coarse Grained Models
366 Progress of Theoretical Physics Supplement No. 138, 2000 Simulating Folding of Helical Proteins with Coarse Grained Models Shoji Takada Department of Chemistry, Kobe University, Kobe 657-8501, Japan
More informationPrediction and refinement of NMR structures from sparse experimental data
Prediction and refinement of NMR structures from sparse experimental data Jeff Skolnick Director Center for the Study of Systems Biology School of Biology Georgia Institute of Technology Overview of talk
More informationMolecular dynamics simulation of Aquaporin-1. 4 nm
Molecular dynamics simulation of Aquaporin-1 4 nm Molecular Dynamics Simulations Schrödinger equation i~@ t (r, R) =H (r, R) Born-Oppenheimer approximation H e e(r; R) =E e (R) e(r; R) Nucleic motion described
More informationMolecular Modeling lecture 2
Molecular Modeling 2018 -- lecture 2 Topics 1. Secondary structure 3. Sequence similarity and homology 2. Secondary structure prediction 4. Where do protein structures come from? X-ray crystallography
More informationUsing NMR to study Macromolecular Interactions. John Gross, BP204A UCSF. Nov 27, 2017
Using NMR to study Macromolecular Interactions John Gross, BP204A UCSF Nov 27, 2017 Outline Review of basic NMR experiment Multidimensional NMR Monitoring ligand binding Structure Determination Review:
More informationPeptides And Proteins
Kevin Burgess, May 3, 2017 1 Peptides And Proteins from chapter(s) in the recommended text A. Introduction B. omenclature And Conventions by amide bonds. on the left, right. 2 -terminal C-terminal triglycine
More informationModel Mélange. Physical Models of Peptides and Proteins
Model Mélange Physical Models of Peptides and Proteins In the Model Mélange activity, you will visit four different stations each featuring a variety of different physical models of peptides or proteins.
More informationPROTEIN STRUCTURE PREDICTION USING GAS PHASE MOLECULAR DYNAMICS SIMULATION: EOTAXIN-3 CYTOKINE AS A CASE STUDY
International Conference Mathematical and Computational Biology 2011 International Journal of Modern Physics: Conference Series Vol. 9 (2012) 193 198 World Scientific Publishing Company DOI: 10.1142/S2010194512005259
More informationMechanical Proteins. Stretching imunoglobulin and fibronectin. domains of the muscle protein titin. Adhesion Proteins of the Immune System
Mechanical Proteins F C D B A domains of the muscle protein titin E Stretching imunoglobulin and fibronectin G NIH Resource for Macromolecular Modeling and Bioinformatics Theoretical Biophysics Group,
More informationCHEM 3653 Exam # 1 (03/07/13)
1. Using phylogeny all living organisms can be divided into the following domains: A. Bacteria, Eukarya, and Vertebrate B. Archaea and Eukarya C. Bacteria, Eukarya, and Archaea D. Eukarya and Bacteria
More informationImplicit modeling of nonpolar solvation for simulating protein folding and conformational transitions
PERSPECTIVE www.rsc.org/pccp Physical Chemistry Chemical Physics Implicit modeling of nonpolar solvation for simulating protein folding and conformational transitions Jianhan Chenw and Charles L. Brooks
More informationDihedral Angles. Homayoun Valafar. Department of Computer Science and Engineering, USC 02/03/10 CSCE 769
Dihedral Angles Homayoun Valafar Department of Computer Science and Engineering, USC The precise definition of a dihedral or torsion angle can be found in spatial geometry Angle between to planes Dihedral
More informationSupplemental Material for Global Langevin model of multidimensional biomolecular dynamics
Supplemental Material for Global Langevin model of multidimensional biomolecular dynamics Norbert Schaudinnus, Benjamin Lickert, Mithun Biswas and Gerhard Stock Biomolecular Dynamics, Institute of Physics,
More informationONETEP PB/SA: Application to G-Quadruplex DNA Stability. Danny Cole
ONETEP PB/SA: Application to G-Quadruplex DNA Stability Danny Cole Introduction Historical overview of structure and free energy calculation of complex molecules using molecular mechanics and continuum
More informationarxiv: v1 [cond-mat.soft] 22 Oct 2007
Conformational Transitions of Heteropolymers arxiv:0710.4095v1 [cond-mat.soft] 22 Oct 2007 Michael Bachmann and Wolfhard Janke Institut für Theoretische Physik, Universität Leipzig, Augustusplatz 10/11,
More informationModeling Background; Donald J. Jacobs, University of North Carolina at Charlotte Page 1 of 8
Modeling Background; Donald J. Jacobs, University of North Carolina at Charlotte Page 1 of 8 Depending on thermodynamic and solvent conditions, the interrelationships between thermodynamic stability of
More information2 Computational methods. 1 Introduction. 2.1 Simulation models
Molecular Dynamics Simulations of Folding and Insertion of the Ebola Virus Fusion Peptide into a Membrane Bilayer Mark A. Olson 1, In-Chul Yeh 2, and Michael S. Lee 1,3 1 US Army Medical Research Institute
More informationOutline. The ensemble folding kinetics of protein G from an all-atom Monte Carlo simulation. Unfolded Folded. What is protein folding?
The ensemble folding kinetics of protein G from an all-atom Monte Carlo simulation By Jun Shimada and Eugine Shaknovich Bill Hawse Dr. Bahar Elisa Sandvik and Mehrdad Safavian Outline Background on protein
More informationApplications of Molecular Dynamics
June 4, 0 Molecular Modeling and Simulation Applications of Molecular Dynamics Agricultural Bioinformatics Research Unit, Graduate School of Agricultural and Life Sciences, The University of Tokyo Tohru
More informationExam I Answer Key: Summer 2006, Semester C
1. Which of the following tripeptides would migrate most rapidly towards the negative electrode if electrophoresis is carried out at ph 3.0? a. gly-gly-gly b. glu-glu-asp c. lys-glu-lys d. val-asn-lys
More informationProgramme Last week s quiz results + Summary Fold recognition Break Exercise: Modelling remote homologues
Programme 8.00-8.20 Last week s quiz results + Summary 8.20-9.00 Fold recognition 9.00-9.15 Break 9.15-11.20 Exercise: Modelling remote homologues 11.20-11.40 Summary & discussion 11.40-12.00 Quiz 1 Feedback
More informationIntroduction solution NMR
2 NMR journey Introduction solution NMR Alexandre Bonvin Bijvoet Center for Biomolecular Research with thanks to Dr. Klaartje Houben EMBO Global Exchange course, IHEP, Beijing April 28 - May 5, 20 3 Topics
More informationGenerating, sampling, and analyzing the dynamical pathways of folding proteins with rare event techniques
Generating, sampling, and analyzing the dynamical pathways of folding proteins with rare event techniques Peter Bolhuis van t Hoff institute for Molecular Sciences University of Amsterdam, The Netherlands
More informationSupporting Information
Supporting Information Micelle-Triggered b-hairpin to a-helix Transition in a 14-Residue Peptide from a Choline-Binding Repeat of the Pneumococcal Autolysin LytA HØctor Zamora-Carreras, [a] Beatriz Maestro,
More informationProtein Structure Bioinformatics Introduction
1 Swiss Institute of Bioinformatics Protein Structure Bioinformatics Introduction Basel, 27. September 2004 Torsten Schwede Biozentrum - Universität Basel Swiss Institute of Bioinformatics Klingelbergstr
More informationEnergetics and Thermodynamics
DNA/Protein structure function analysis and prediction Protein Folding and energetics: Introduction to folding Folding and flexibility (Ch. 6) Energetics and Thermodynamics 1 Active protein conformation
More informationSecondary and sidechain structures
Lecture 2 Secondary and sidechain structures James Chou BCMP201 Spring 2008 Images from Petsko & Ringe, Protein Structure and Function. Branden & Tooze, Introduction to Protein Structure. Richardson, J.
More informationMolecular Dynamics Flexible Fitting
Molecular Dynamics Flexible Fitting Ryan McGreevy Research Programmer University of Illinois at Urbana-Champaign NIH Resource for Macromolecular Modeling and Bioinformatics Molecular Dynamics Flexible
More informationProperties of amino acids in proteins
Properties of amino acids in proteins one of the primary roles of DNA (but not the only one!) is to code for proteins A typical bacterium builds thousands types of proteins, all from ~20 amino acids repeated
More informationSecondary structure stability, beta-sheet formation & stability
Protein Physics 2016 Lecture 6, February 5 Secondary structure stability, beta-sheet formation & stability Magnus Andersson magnus.andersson@scilifelab.se Theoretical & Computational Biophysics Recap of
More informationCHEM-UA 652: Thermodynamics and Kinetics
1 CHEM-UA 652: Thermodynamics and Kinetics Notes for Lecture 11 I. PHYSICAL AND CHEMICAL RELEVANCE OF FREE ENERGY In this section, we will consider some examples showing the significance of free energies.
More informationPresenter: She Zhang
Presenter: She Zhang Introduction Dr. David Baker Introduction Why design proteins de novo? It is not clear how non-covalent interactions favor one specific native structure over many other non-native
More informationAggregation of the Amyloid-β Protein: Monte Carlo Optimization Study
John von Neumann Institute for Computing Aggregation of the Amyloid-β Protein: Monte Carlo Optimization Study S. M. Gopal, K. V. Klenin, W. Wenzel published in From Computational Biophysics to Systems
More informationTHE UNIVERSITY OF MANITOBA. PAPER NO: 409 LOCATION: Fr. Kennedy Gold Gym PAGE NO: 1 of 6 DEPARTMENT & COURSE NO: CHEM 4630 TIME: 3 HOURS
PAPER NO: 409 LOCATION: Fr. Kennedy Gold Gym PAGE NO: 1 of 6 DEPARTMENT & COURSE NO: CHEM 4630 TIME: 3 HOURS EXAMINATION: Biochemistry of Proteins EXAMINER: J. O'Neil Section 1: You must answer all of
More informationCentral Dogma. modifications genome transcriptome proteome
entral Dogma DA ma protein post-translational modifications genome transcriptome proteome 83 ierarchy of Protein Structure 20 Amino Acids There are 20 n possible sequences for a protein of n residues!
More informationStructural Perspectives on Drug Resistance
Structural Perspectives on Drug Resistance Irene Weber Departments of Biology and Chemistry Molecular Basis of Disease Program Georgia State University Atlanta, GA, USA What have we learned from 20 years
More information