Protein Folding by Robotics

Size: px
Start display at page:

Download "Protein Folding by Robotics"

Transcription

1 Protein Folding by Robotics 1 TBI Winterseminar 2006 February 21, 2006

2 Protein Folding by Robotics 1 TBI Winterseminar 2006 February 21, 2006

3 Protein Folding by Robotics Probabilistic Roadmap Planning (PRM) [ Thomas, Song, Amato. Protein folding by motion planning. Phys. Biol., 2005 ] 1 TBI Winterseminar 2006 February 21, 2006

4 Aims Find good quality folding paths (into given native structure) no structure prediction! Predict formation orders (of secondary structure) 2 TBI Winterseminar 2006 February 21, 2006

5 Motion planning Motion planning 3 TBI Winterseminar 2006 February 21, 2006

6 Motion planning Motion planning 3 TBI Winterseminar 2006 February 21, 2006

7 Motion planning Motion planning Probabilistic roadmap planing 3 TBI Winterseminar 2006 February 21, 2006

8 Motion planning Motion planning Probabilistic roadmap planing Sampling of configuration space Q 3 TBI Winterseminar 2006 February 21, 2006

9 Motion planning Motion planning Probabilistic roadmap planing Sampling of configuration space Q Connecting nearest configurations by a (simple) local planner 3 TBI Winterseminar 2006 February 21, 2006

10 Motion planning Motion planning Probabilistic roadmap planing Sampling of configuration space Q Connecting nearest configurations by a (simple) local planner Apply graph algorithms to roadmap : Find shortest path 3 TBI Winterseminar 2006 February 21, 2006

11 More on PRM for motion planning tree-like robots 4 TBI Winterseminar 2006 February 21, 2006

12 More on PRM for motion planning tree-like robots (articulated robots) Articulated Joint 4 TBI Winterseminar 2006 February 21, 2006

13 More on PRM for motion planning tree-like robots (articulated robots) Articulated Joint 4 TBI Winterseminar 2006 February 21, 2006

14 More on PRM for motion planning tree-like robots (articulated robots) Articulated Joint configuration = vector of angles configuration space Q = {q q S n } S set of angles n number of angles = degrees of freedom (dof) 4 TBI Winterseminar 2006 February 21, 2006

15 Proteins are Robots (aren t they?) Obvious similarity == 5 TBI Winterseminar 2006 February 21, 2006

16 Proteins are Robots (aren t they?) Obvious similarity ;-) == Our model O O N C C N C C N C C O 5 TBI Winterseminar 2006 February 21, 2006

17 Proteins are Robots (aren t they?) Obvious similarity ;-) == Our model O O N C C N C C N C C O 5 TBI Winterseminar 2006 February 21, 2006

18 Proteins are Robots (aren t they?) Obvious similarity ;-) == Our model O O N C phi C psi N C C N C C O 5 TBI Winterseminar 2006 February 21, 2006

19 Proteins are Robots (aren t they?) Obvious similarity ;-) /== Our model O O N C phi C psi N C C N C C O 5 TBI Winterseminar 2006 February 21, 2006

20 Proteins are Robots (aren t they?) Obvious similarity ;-) /== Our model O O N C phi C psi N C C N C C O Protein == vector of phi and psi angles (treelike robot with 2n dof) possible models range from only backbone up to full atom 5 TBI Winterseminar 2006 February 21, 2006

21 Differences to usual PRM no external obstacles, but self-avoidingness torsion angles quality of paths low energy intermediate states kinetically prefered paths highly probable paths 6 TBI Winterseminar 2006 February 21, 2006

22 Energy Function method can use any potential 7 TBI Winterseminar 2006 February 21, 2006

23 Energy Function method can use any potential Our coarse potential [Levitt. J.Mol.Biol., ] each sidechain by only one atom (zero dof) U tot = 7 TBI Winterseminar 2006 February 21, 2006

24 Energy Function method can use any potential Our coarse potential [Levitt. J.Mol.Biol., ] each sidechain by only one atom (zero dof) U tot = K d {[(d i d 0 ) 2 + dc 2 ] 2 1 d c } restraints first term favors known secondary structure through main chain hydrogen bonds and disulphide bonds 7 TBI Winterseminar 2006 February 21, 2006

25 Energy Function method can use any potential Our coarse potential [Levitt. J.Mol.Biol., ] each sidechain by only one atom (zero dof) U tot = K d {[(d i d 0 ) 2 + dc 2 ] 2 1 d c } + E hp restraints first term favors known secondary structure through main chain hydrogen bonds and disulphide bonds second term hydrophobic effect 7 TBI Winterseminar 2006 February 21, 2006

26 Energy Function method can use any potential Our coarse potential [Levitt. J.Mol.Biol., ] each sidechain by only one atom (zero dof) U tot = K d {[(d i d 0 ) 2 + dc 2 ] 2 1 d c } + E hp restraints first term favors known secondary structure through main chain hydrogen bonds and disulphide bonds second term hydrophobic effect Van der Waals interaction modeled by step function 7 TBI Winterseminar 2006 February 21, 2006

27 Energy Function method can use any potential Our coarse potential [Levitt. J.Mol.Biol., ] each sidechain by only one atom (zero dof) U tot = K d {[(d i d 0 ) 2 + dc 2 ] 2 1 d c } + E hp restraints first term favors known secondary structure through main chain hydrogen bonds and disulphide bonds second term hydrophobic effect Van der Waals interaction modeled by step function All-atom potential: EEF1 [Lazaridis, Karplus. Proteins, ] 7 TBI Winterseminar 2006 February 21, 2006

28 PRM method for Proteins 8 TBI Winterseminar 2006 February 21, 2006

29 PRM method for Proteins Sampling 8 TBI Winterseminar 2006 February 21, 2006

30 PRM method for Proteins Sampling Connecting 8 TBI Winterseminar 2006 February 21, 2006

31 PRM method for Proteins Sampling Connecting Extracting 8 TBI Winterseminar 2006 February 21, 2006

32 Sampling Node Generation Sampling Connecting Extracting 9 TBI Winterseminar 2006 February 21, 2006

33 Node Generation No uniform sampling configuration space too large need biased sampling strategy 10 TBI Winterseminar 2006 February 21, 2006

34 Node Generation No uniform sampling configuration space too large need biased sampling strategy Gaussian sampling centered around native conformation with different STDs 5,10,...,160 ensure representants for different numbers of native contacts 10 TBI Winterseminar 2006 February 21, 2006

35 Node Generation No uniform sampling configuration space too large need biased sampling strategy Gaussian sampling centered around native conformation with different STDs 5,10,...,160 ensure representants for different numbers of native contacts Selection by energy P(accept q) = 1 if E(q) < E min E max E(q) E max E if E min min E(q) E max 0 if E(q) > E max 10 TBI Winterseminar 2006 February 21, 2006

36 More on Node Generation Visualization of Sampling Strategy 11 TBI Winterseminar 2006 February 21, 2006

37 More on Node Generation Visualization of Sampling Strategy Distribution Psi and Phi angles RMSD vs. Energy 11 TBI Winterseminar 2006 February 21, 2006

38 Node Connection Sampling Connecting Extracting 12 TBI Winterseminar 2006 February 21, 2006

39 Connecting Nodes by Local Planner connect configurations in close distance generate N intermediary nodes by local planner 13 TBI Winterseminar 2006 February 21, 2006

40 Connecting Nodes by Local Planner connect configurations in close distance generate N intermediary nodes by local planner 13 TBI Winterseminar 2006 February 21, 2006

41 Connecting Nodes by Local Planner connect configurations in close distance generate N intermediary nodes by local planner 13 TBI Winterseminar 2006 February 21, 2006

42 Connecting Nodes by Local Planner connect configurations in close distance generate N intermediary nodes by local planner 13 TBI Winterseminar 2006 February 21, 2006

43 Connecting Nodes by Local Planner connect configurations in close distance generate N intermediary nodes by local planner 13 TBI Winterseminar 2006 February 21, 2006

44 Connecting Nodes by Local Planner connect configurations in close distance generate N intermediary nodes by local planner 13 TBI Winterseminar 2006 February 21, 2006

45 Connecting Nodes by Local Planner connect configurations in close distance generate N intermediary nodes by local planner P5 P4 P1 P2 P3 assign weights to edges { e E kt if E > 0 P i = 1 if E 0 13 TBI Winterseminar 2006 February 21, 2006

46 Connecting Nodes by Local Planner connect configurations in close distance generate N intermediary nodes by local planner P5 Weight P4 P1 P2 P3 assign weights to edges { e E kt if E > 0 P i = 1 if E 0 Weight = N log(p i ) i=0 13 TBI Winterseminar 2006 February 21, 2006

47 Connecting Nodes by Local Planner connect configurations in close distance generate N intermediary nodes by local planner P5 Weight P4 P1 P2 P3 assign weights to edges { e E kt if E > 0 P i = 1 if E 0 Weight = N log(p i ) i=0 13 TBI Winterseminar 2006 February 21, 2006

48 Extracting Paths Sampling Connecting Extracting 14 TBI Winterseminar 2006 February 21, 2006

49 Extracting Paths Shortest Path extract one shortest path from some starting conformation, one path at a time 15 TBI Winterseminar 2006 February 21, 2006

50 Extracting Paths Shortest Path extract one shortest path from some starting conformation, one path at a time Single Source Shortest Paths (SSSP) extract shortest paths from all starting conformation compute paths simultaneously generate tree of shortest paths (SSSP tree) 15 TBI Winterseminar 2006 February 21, 2006

51 Big Picture Sampling Connecting Extracting 16 TBI Winterseminar 2006 February 21, 2006

52 Studied Proteins Overview of studied proteins, roadmap size, and construction times 17 TBI Winterseminar 2006 February 21, 2006

53 Formation orders formation order of secondary structure for verifying method formation orders can be determined experimentally [ Li, Woodward. Protein Science, ] Pulse labeling Out-exchange prediction of formation orders single paths averaging over multiple paths (SSSP-tree) 18 TBI Winterseminar 2006 February 21, 2006

54 Timed Contact Maps 19 TBI Winterseminar 2006 February 21, 2006

55 Formation Order 20 TBI Winterseminar 2006 February 21, 2006

56 Formation Order no (reported) contradictions between prediction and validation different kind of information from experiment and prediction 20 TBI Winterseminar 2006 February 21, 2006

57 The Proteins G and L Studied in more detail good test case structurally similar: 1α + 4β fold differently Protein G: β-turn 2 forms first Protein L: β-turn 1 forms first 21 TBI Winterseminar 2006 February 21, 2006

58 Comparison of Analysis Techniques β-turn Formation 22 TBI Winterseminar 2006 February 21, 2006

59 Conclusion 23 TBI Winterseminar 2006 February 21, 2006

60 Conclusion PRM can be applied to realistic protein models 23 TBI Winterseminar 2006 February 21, 2006

61 Conclusion PRM can be applied to realistic protein models Introduced method makes verifiable prediction 23 TBI Winterseminar 2006 February 21, 2006

62 Conclusion PRM can be applied to realistic protein models Introduced method makes verifiable prediction Coarse potential is sufficient 23 TBI Winterseminar 2006 February 21, 2006

63 Conclusion PRM can be applied to realistic protein models Introduced method makes verifiable prediction Coarse potential is sufficient Predictions are in good accordance to experimental data 23 TBI Winterseminar 2006 February 21, 2006

64 Conclusion PRM can be applied to realistic protein models Introduced method makes verifiable prediction Coarse potential is sufficient Predictions are in good accordance to experimental data Interesting relations to e.g. computation of barrier trees 23 TBI Winterseminar 2006 February 21, 2006

A Path Planning-Based Study of Protein Folding with a Case Study of Hairpin Formation in Protein G and L

A Path Planning-Based Study of Protein Folding with a Case Study of Hairpin Formation in Protein G and L A Path Planning-Based Study of Protein Folding with a Case Study of Hairpin Formation in Protein G and L G. Song, S. Thomas, K.A. Dill, J.M. Scholtz, N.M. Amato Pacific Symposium on Biocomputing 8:240-251(2003)

More information

Protein folding by motion planning

Protein folding by motion planning INSTITUTE OFPHYSICS PUBLISHING PHYSICAL BIOLOGY Phys. Biol. 2 () S148 S5 doi:1.188/1478-3975/2/4/s9 Protein folding by motion planning Shawna Thomas 1, Guang Song 2 and Nancy M Amato 1,3 1 Department of

More information

A Motion Planning Approach to Folding:

A Motion Planning Approach to Folding: A Motion Planning Approach to Folding: From Paper Craft to Protein Folding Λ Guang Song Nancy M. Amato Department of Computer Science Texas A&M University College Station, TX 77843-32 fgsong,amatog@cs.tamu.edu

More information

Using Motion Planning to Study Protein Folding Pathways Λ

Using Motion Planning to Study Protein Folding Pathways Λ Using Motion Planning to Study Protein Folding Pathways Λ Guang Song Department of Computer Science Texas A&M University College Station, TX 77843-3 gsong@cs.tamu.edu Nancy M. Amato Department of Computer

More information

DISCRETE TUTORIAL. Agustí Emperador. Institute for Research in Biomedicine, Barcelona APPLICATION OF DISCRETE TO FLEXIBLE PROTEIN-PROTEIN DOCKING:

DISCRETE TUTORIAL. Agustí Emperador. Institute for Research in Biomedicine, Barcelona APPLICATION OF DISCRETE TO FLEXIBLE PROTEIN-PROTEIN DOCKING: DISCRETE TUTORIAL Agustí Emperador Institute for Research in Biomedicine, Barcelona APPLICATION OF DISCRETE TO FLEXIBLE PROTEIN-PROTEIN DOCKING: STRUCTURAL REFINEMENT OF DOCKING CONFORMATIONS Emperador

More information

Molecular Modelling. part of Bioinformatik von RNA- und Proteinstrukturen. Sonja Prohaska. Leipzig, SS Computational EvoDevo University Leipzig

Molecular Modelling. part of Bioinformatik von RNA- und Proteinstrukturen. Sonja Prohaska. Leipzig, SS Computational EvoDevo University Leipzig part of Bioinformatik von RNA- und Proteinstrukturen Computational EvoDevo University Leipzig Leipzig, SS 2011 Protein Structure levels or organization Primary structure: sequence of amino acids (from

More information

Outline. The ensemble folding kinetics of protein G from an all-atom Monte Carlo simulation. Unfolded Folded. What is protein folding?

Outline. 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 information

Conformational Geometry of Peptides and Proteins:

Conformational Geometry of Peptides and Proteins: Conformational Geometry of Peptides and Proteins: Before discussing secondary structure, it is important to appreciate the conformational plasticity of proteins. Each residue in a polypeptide has three

More information

The protein folding problem consists of two parts:

The protein folding problem consists of two parts: Energetics and kinetics of protein folding The protein folding problem consists of two parts: 1)Creating a stable, well-defined structure that is significantly more stable than all other possible structures.

More information

Supporting Online Material for

Supporting Online Material for www.sciencemag.org/cgi/content/full/309/5742/1868/dc1 Supporting Online Material for Toward High-Resolution de Novo Structure Prediction for Small Proteins Philip Bradley, Kira M. S. Misura, David Baker*

More information

A Motion Planning Approach to Flexible Ligand Binding

A Motion Planning Approach to Flexible Ligand Binding From: ISMB-99 Proceedings. Copyright 1999, AAAI (www.aaai.org). All rights reserved. A Motion Planning Approach to Flexible Ligand Binding Amit P. Singh 1, Jean-Claude Latombe 2, Douglas L. Brutlag 3 1

More information

Many proteins spontaneously refold into native form in vitro with high fidelity and high speed.

Many 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 information

Lecture 2 and 3: Review of forces (ctd.) and elementary statistical mechanics. Contributions to protein stability

Lecture 2 and 3: Review of forces (ctd.) and elementary statistical mechanics. Contributions to protein stability Lecture 2 and 3: Review of forces (ctd.) and elementary statistical mechanics. Contributions to protein stability Part I. Review of forces Covalent bonds Non-covalent Interactions: Van der Waals Interactions

More information

Protein Folding experiments and theory

Protein Folding experiments and theory Protein Folding experiments and theory 1, 2,and 3 Protein Structure Fig. 3-16 from Lehninger Biochemistry, 4 th ed. The 3D structure is not encoded at the single aa level Hydrogen Bonding Shared H atom

More information

Useful background reading

Useful background reading Overview of lecture * General comment on peptide bond * Discussion of backbone dihedral angles * Discussion of Ramachandran plots * Description of helix types. * Description of structures * NMR patterns

More information

Can 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? 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 information

Lecture 2-3: Review of forces (ctd.) and elementary statistical mechanics. Contributions to protein stability

Lecture 2-3: Review of forces (ctd.) and elementary statistical mechanics. Contributions to protein stability Lecture 2-3: Review of forces (ctd.) and elementary statistical mechanics. Contributions to protein stability Part I. Review of forces Covalent bonds Non-covalent Interactions Van der Waals Interactions

More information

= (-22) = +2kJ /mol

= (-22) = +2kJ /mol Lecture 8: Thermodynamics & Protein Stability Assigned reading in Campbell: Chapter 4.4-4.6 Key Terms: DG = -RT lnk eq = DH - TDS Transition Curve, Melting Curve, Tm DH calculation DS calculation van der

More information

HOMOLOGY MODELING. The sequence alignment and template structure are then used to produce a structural model of the target.

HOMOLOGY MODELING. The sequence alignment and template structure are then used to produce a structural model of the target. HOMOLOGY MODELING Homology modeling, also known as comparative modeling of protein refers to constructing an atomic-resolution model of the "target" protein from its amino acid sequence and an experimental

More information

Protein Folding & Stability. Lecture 11: Margaret A. Daugherty. Fall How do we go from an unfolded polypeptide chain to a

Protein Folding & Stability. Lecture 11: Margaret A. Daugherty. Fall How do we go from an unfolded polypeptide chain to a Lecture 11: Protein Folding & Stability Margaret A. Daugherty Fall 2004 How do we go from an unfolded polypeptide chain to a compact folded protein? (Folding of thioredoxin, F. Richards) Structure - Function

More information

Lecture 18 Generalized Belief Propagation and Free Energy Approximations

Lecture 18 Generalized Belief Propagation and Free Energy Approximations Lecture 18, Generalized Belief Propagation and Free Energy Approximations 1 Lecture 18 Generalized Belief Propagation and Free Energy Approximations In this lecture we talked about graphical models and

More information

Lecture 11: Protein Folding & Stability

Lecture 11: Protein Folding & Stability Structure - Function Protein Folding: What we know Lecture 11: Protein Folding & Stability 1). Amino acid sequence dictates structure. 2). The native structure represents the lowest energy state for a

More information

Protein Folding & Stability. Lecture 11: Margaret A. Daugherty. Fall Protein Folding: What we know. Protein Folding

Protein Folding & Stability. Lecture 11: Margaret A. Daugherty. Fall Protein Folding: What we know. Protein Folding Lecture 11: Protein Folding & Stability Margaret A. Daugherty Fall 2003 Structure - Function Protein Folding: What we know 1). Amino acid sequence dictates structure. 2). The native structure represents

More information

Prediction and refinement of NMR structures from sparse experimental data

Prediction 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 information

A Multi-Directional Rapidly Exploring Random Graph (mrrg) for Protein Folding

A Multi-Directional Rapidly Exploring Random Graph (mrrg) for Protein Folding A Multi-Directional Rapidly Exploring Random Graph (mrrg) for Protein Folding Shuvra Kanti Nath, Shawna Thomas, Chinwe Ekenna, and Nancy M. Amato Parasol Lab, Department of Computer Science and Engineering

More information

Protein Structure Prediction

Protein Structure Prediction Page 1 Protein Structure Prediction Russ B. Altman BMI 214 CS 274 Protein Folding is different from structure prediction --Folding is concerned with the process of taking the 3D shape, usually based on

More information

Packing of Secondary Structures

Packing 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 information

Presenter: She Zhang

Presenter: 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 information

Polypeptide Folding Using Monte Carlo Sampling, Concerted Rotation, and Continuum Solvation

Polypeptide 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 information

Biomolecules: lecture 9

Biomolecules: lecture 9 Biomolecules: lecture 9 - understanding further why amino acids are the building block for proteins - understanding the chemical properties amino acids bring to proteins - realizing that many proteins

More information

PRM Roadmap Query. PRM Roadmap after Node Generation. goal. C obst. C obst. C obst. C obst. C obst. C obst. C obst. C obst. C Space. start.

PRM Roadmap Query. PRM Roadmap after Node Generation. goal. C obst. C obst. C obst. C obst. C obst. C obst. C obst. C obst. C Space. start. Using Motion Planning to Map Protein Folding Landscapes and Analyze Folding Kinetics of Known Native Structures 1 Nancy M. Amato amato@cs.tamu.edu Guang Song gsong@cs.tamu.edu Technical Report TR1-1 PARASOL

More information

From Amino Acids to Proteins - in 4 Easy Steps

From Amino Acids to Proteins - in 4 Easy Steps From Amino Acids to Proteins - in 4 Easy Steps Although protein structure appears to be overwhelmingly complex, you can provide your students with a basic understanding of how proteins fold by focusing

More information

Protein Folding Prof. Eugene Shakhnovich

Protein 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 information

NMR, X-ray Diffraction, Protein Structure, and RasMol

NMR, X-ray Diffraction, Protein Structure, and RasMol NMR, X-ray Diffraction, Protein Structure, and RasMol Introduction So far we have been mostly concerned with the proteins themselves. The techniques (NMR or X-ray diffraction) used to determine a structure

More information

arxiv: v1 [cond-mat.soft] 22 Oct 2007

arxiv: 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 information

STRUCTURAL BIOINFORMATICS. Barry Grant University of Michigan

STRUCTURAL BIOINFORMATICS. Barry Grant University of Michigan STRUCTURAL BIOINFORMATICS Barry Grant University of Michigan www.thegrantlab.org bjgrant@umich.edu Bergen, Norway 28-Sep-2015 Objective: Provide an introduction to the practice of structural bioinformatics,

More information

Protein function prediction via graph kernels

Protein function prediction via graph kernels Protein function prediction via graph kernels Karsten Borgwardt Joint work with Cheng oon Ong and.v.n. Vishwanathan, tefan chönauer, Hans-Peter Kriegel and Alex mola appeared in IMB 2005 1 Content Introduction

More information

Introduction to Computational Structural Biology

Introduction to Computational Structural Biology Introduction to Computational Structural Biology Part I 1. Introduction The disciplinary character of Computational Structural Biology The mathematical background required and the topics covered Bibliography

More information

Molecular Mechanics, Dynamics & Docking

Molecular Mechanics, Dynamics & Docking Molecular Mechanics, Dynamics & Docking Lawrence Hunter, Ph.D. Director, Computational Bioscience Program University of Colorado School of Medicine Larry.Hunter@uchsc.edu http://compbio.uchsc.edu/hunter

More information

Dihedral 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 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 information

Protein Structure Basics

Protein Structure Basics Protein Structure Basics Presented by Alison Fraser, Christine Lee, Pradhuman Jhala, Corban Rivera Importance of Proteins Muscle structure depends on protein-protein interactions Transport across membranes

More information

CMPS 6630: Introduction to Computational Biology and Bioinformatics. Tertiary Structure Prediction

CMPS 6630: Introduction to Computational Biology and Bioinformatics. Tertiary Structure Prediction CMPS 6630: Introduction to Computational Biology and Bioinformatics Tertiary Structure Prediction Tertiary Structure Prediction Why Should Tertiary Structure Prediction Be Possible? Molecules obey the

More information

LS1a Midterm Exam 1 Review Session Problems

LS1a Midterm Exam 1 Review Session Problems LS1a Midterm Exam 1 Review Session Problems 1. n aqueous mixture of a weak acid and its conjugate base is often used in the laboratory to prepare solutions referred to as buffers. ne commonly used acid

More information

Short Announcements. 1 st Quiz today: 15 minutes. Homework 3: Due next Wednesday.

Short Announcements. 1 st Quiz today: 15 minutes. Homework 3: Due next Wednesday. Short Announcements 1 st Quiz today: 15 minutes Homework 3: Due next Wednesday. Next Lecture, on Visualizing Molecular Dynamics (VMD) by Klaus Schulten Today s Lecture: Protein Folding, Misfolding, Aggregation

More information

CMPS 3110: Bioinformatics. Tertiary Structure Prediction

CMPS 3110: Bioinformatics. Tertiary Structure Prediction CMPS 3110: Bioinformatics Tertiary Structure Prediction Tertiary Structure Prediction Why Should Tertiary Structure Prediction Be Possible? Molecules obey the laws of physics! Conformation space is finite

More information

Template Free Protein Structure Modeling Jianlin Cheng, PhD

Template Free Protein Structure Modeling Jianlin Cheng, PhD Template Free Protein Structure Modeling Jianlin Cheng, PhD Associate Professor Computer Science Department Informatics Institute University of Missouri, Columbia 2013 Protein Energy Landscape & Free Sampling

More information

arxiv:cond-mat/ v1 2 Feb 94

arxiv:cond-mat/ v1 2 Feb 94 cond-mat/9402010 Properties and Origins of Protein Secondary Structure Nicholas D. Socci (1), William S. Bialek (2), and José Nelson Onuchic (1) (1) Department of Physics, University of California at San

More information

CS273: Algorithms for Structure Handout # 2 and Motion in Biology Stanford University Thursday, 1 April 2004

CS273: Algorithms for Structure Handout # 2 and Motion in Biology Stanford University Thursday, 1 April 2004 CS273: Algorithms for Structure Handout # 2 and Motion in Biology Stanford University Thursday, 1 April 2004 Lecture #2: 1 April 2004 Topics: Kinematics : Concepts and Results Kinematics of Ligands and

More information

Multi-Scale Hierarchical Structure Prediction of Helical Transmembrane Proteins

Multi-Scale Hierarchical Structure Prediction of Helical Transmembrane Proteins Multi-Scale Hierarchical Structure Prediction of Helical Transmembrane Proteins Zhong Chen Dept. of Biochemistry and Molecular Biology University of Georgia, Athens, GA 30602 Email: zc@csbl.bmb.uga.edu

More information

BIOCHEMISTRY Unit 2 Part 4 ACTIVITY #6 (Chapter 5) PROTEINS

BIOCHEMISTRY Unit 2 Part 4 ACTIVITY #6 (Chapter 5) PROTEINS BIOLOGY BIOCHEMISTRY Unit 2 Part 4 ACTIVITY #6 (Chapter 5) NAME NAME PERIOD PROTEINS GENERAL CHARACTERISTICS AND IMPORTANCES: Polymers of amino acids Each has unique 3-D shape Vary in sequence of amino

More information

Stochastic Conformational Roadmaps for Computing Ensemble Properties of Molecular Motion

Stochastic Conformational Roadmaps for Computing Ensemble Properties of Molecular Motion Stochastic Conformational Roadmaps for Computing Ensemble Properties of Molecular Motion Mehmet Serkan Apaydın, Douglas L. Brutlag, Carlos Guestrin, David Hsu, and Jean-Claude Latombe Stanford University,

More information

Biology Chemistry & Physics of Biomolecules. Examination #1. Proteins Module. September 29, Answer Key

Biology 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 information

MCB100A/Chem130 MidTerm Exam 2 April 4, 2013

MCB100A/Chem130 MidTerm Exam 2 April 4, 2013 MCBA/Chem Miderm Exam 2 April 4, 2 Name Student ID rue/false (2 points each).. he Boltzmann constant, k b sets the energy scale for observing energy microstates 2. Atoms with favorable electronic configurations

More information

Protein Structure Determination Using NMR Restraints BCMB/CHEM 8190

Protein Structure Determination Using NMR Restraints BCMB/CHEM 8190 Protein Structure Determination Using NMR Restraints BCMB/CHEM 8190 Programs for NMR Based Structure Determination CNS - Brünger, A. T.; Adams, P. D.; Clore, G. M.; DeLano, W. L.; Gros, P.; Grosse-Kunstleve,

More information

Biochemistry,530:,, Introduc5on,to,Structural,Biology, Autumn,Quarter,2015,

Biochemistry,530:,, Introduc5on,to,Structural,Biology, Autumn,Quarter,2015, Biochemistry,530:,, Introduc5on,to,Structural,Biology, Autumn,Quarter,2015, Course,Informa5on, BIOC%530% GraduateAlevel,discussion,of,the,structure,,func5on,,and,chemistry,of,proteins,and, nucleic,acids,,control,of,enzyma5c,reac5ons.,please,see,the,course,syllabus,and,

More information

Template Free Protein Structure Modeling Jianlin Cheng, PhD

Template Free Protein Structure Modeling Jianlin Cheng, PhD Template Free Protein Structure Modeling Jianlin Cheng, PhD Professor Department of EECS Informatics Institute University of Missouri, Columbia 2018 Protein Energy Landscape & Free Sampling http://pubs.acs.org/subscribe/archive/mdd/v03/i09/html/willis.html

More information

Proteins are not rigid structures: Protein dynamics, conformational variability, and thermodynamic stability

Proteins are not rigid structures: Protein dynamics, conformational variability, and thermodynamic stability Proteins are not rigid structures: Protein dynamics, conformational variability, and thermodynamic stability Dr. Andrew Lee UNC School of Pharmacy (Div. Chemical Biology and Medicinal Chemistry) UNC Med

More information

Design of a Novel Globular Protein Fold with Atomic-Level Accuracy

Design of a Novel Globular Protein Fold with Atomic-Level Accuracy Design of a Novel Globular Protein Fold with Atomic-Level Accuracy Brian Kuhlman, Gautam Dantas, Gregory C. Ireton, Gabriele Varani, Barry L. Stoddard, David Baker Presented by Kate Stafford 4 May 05 Protein

More information

Alpha-helical Topology and Tertiary Structure Prediction of Globular Proteins Scott R. McAllister Christodoulos A. Floudas Princeton University

Alpha-helical Topology and Tertiary Structure Prediction of Globular Proteins Scott R. McAllister Christodoulos A. Floudas Princeton University Alpha-helical Topology and Tertiary Structure Prediction of Globular Proteins Scott R. McAllister Christodoulos A. Floudas Princeton University Department of Chemical Engineering Program of Applied and

More information

Unfolding CspB by means of biased molecular dynamics

Unfolding CspB by means of biased molecular dynamics Chapter 4 Unfolding CspB by means of biased molecular dynamics 4.1 Introduction Understanding the mechanism of protein folding has been a major challenge for the last twenty years, as pointed out in the

More information

ALL LECTURES IN SB Introduction

ALL LECTURES IN SB Introduction 1. Introduction 2. Molecular Architecture I 3. Molecular Architecture II 4. Molecular Simulation I 5. Molecular Simulation II 6. Bioinformatics I 7. Bioinformatics II 8. Prediction I 9. Prediction II ALL

More information

PROTEIN'STRUCTURE'DETERMINATION'

PROTEIN'STRUCTURE'DETERMINATION' PROTEIN'STRUCTURE'DETERMINATION' USING'NMR'RESTRAINTS' BCMB/CHEM'8190' Programs for NMR Based Structure Determination CNS - Brünger, A. T.; Adams, P. D.; Clore, G. M.; DeLano, W. L.; Gros, P.; Grosse-Kunstleve,

More information

Protein folding. Today s Outline

Protein 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 information

Robotics. Path Planning. Marc Toussaint U Stuttgart

Robotics. Path Planning. Marc Toussaint U Stuttgart Robotics Path Planning Path finding vs. trajectory optimization, local vs. global, Dijkstra, Probabilistic Roadmaps, Rapidly Exploring Random Trees, non-holonomic systems, car system equation, path-finding

More information

Programming Robots in ROS Slides adapted from CMU, Harvard and TU Stuttgart

Programming Robots in ROS Slides adapted from CMU, Harvard and TU Stuttgart Programming Robots in ROS Slides adapted from CMU, Harvard and TU Stuttgart Path Planning Problem Given an initial configuration q_start and a goal configuration q_goal,, we must generate the best continuous

More information

Papers listed: Cell2. This weeks papers. Chapt 4. Protein structure and function. The importance of proteins

Papers listed: Cell2. This weeks papers. Chapt 4. Protein structure and function. The importance of proteins 1 Papers listed: Cell2 During the semester I will speak of information from several papers. For many of them you will not be required to read these papers, however, you can do so for the fun of it (and

More information

Bioinformatics. Macromolecular structure

Bioinformatics. Macromolecular structure Bioinformatics Macromolecular structure Contents Determination of protein structure Structure databases Secondary structure elements (SSE) Tertiary structure Structure analysis Structure alignment Domain

More information

Timescales of Protein Dynamics

Timescales 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 information

Using Dimensionality Reduction to Better Capture RNA and Protein Folding Motions

Using Dimensionality Reduction to Better Capture RNA and Protein Folding Motions Using Dimensionality Reduction to Better Capture RNA and Protein Folding Motions Lydia Tapia Shawna Thomas Nancy M. Amato ltapia@cs.tamu.edu sthomas@cs.tamu.edu amato@cs.tamu.edu Technical Report TR8-5

More information

April, The energy functions include:

April, The energy functions include: REDUX A collection of Python scripts for torsion angle Monte Carlo protein molecular simulations and analysis The program is based on unified residue peptide model and is designed for more efficient exploration

More information

Constraining Protein Structure Simulation with Aggregate Data Using Robotic Techniques

Constraining Protein Structure Simulation with Aggregate Data Using Robotic Techniques Grand Valley State University ScholarWorks@GVSU Masters Theses Graduate Research and Creative Practice 4-26-2007 Constraining Protein Structure Simulation with Aggregate Data Using Robotic Techniques Eric

More information

Robotics. Path Planning. University of Stuttgart Winter 2018/19

Robotics. Path Planning. University of Stuttgart Winter 2018/19 Robotics Path Planning Path finding vs. trajectory optimization, local vs. global, Dijkstra, Probabilistic Roadmaps, Rapidly Exploring Random Trees, non-holonomic systems, car system equation, path-finding

More information

Hyeyoung Shin a, Tod A. Pascal ab, William A. Goddard III abc*, and Hyungjun Kim a* Korea

Hyeyoung Shin a, Tod A. Pascal ab, William A. Goddard III abc*, and Hyungjun Kim a* Korea The Scaled Effective Solvent Method for Predicting the Equilibrium Ensemble of Structures with Analysis of Thermodynamic Properties of Amorphous Polyethylene Glycol-Water Mixtures Hyeyoung Shin a, Tod

More information

Monte Carlo simulations of polyalanine using a reduced model and statistics-based interaction potentials

Monte 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 information

BIMS 503 Exam I. Sign Pledge Here: Questions from Robert Nakamoto (40 pts. Total)

BIMS 503 Exam I. Sign Pledge Here: Questions from Robert Nakamoto (40 pts. Total) BIMS 503 Exam I September 24, 2007 _ /email: Sign Pledge Here: Questions from Robert Nakamoto (40 pts. Total) Questions 1-6 refer to this situation: You are able to partially purify an enzyme activity

More information

Timescales of Protein Dynamics

Timescales 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 information

Molecular Dynamics, Monte Carlo and Docking. Lecture 21. Introduction to Bioinformatics MNW2

Molecular Dynamics, Monte Carlo and Docking. Lecture 21. Introduction to Bioinformatics MNW2 Molecular Dynamics, Monte Carlo and Docking Lecture 21 Introduction to Bioinformatics MNW2 Allowed phi-psi angles Red areas are preferred, yellow areas are allowed, and white is avoided 2.3a Hamiltonian

More information

Computational Genomics

Computational Genomics Computational Genomics http://www.cs.cmu.edu/~02710 Introduction to probability, statistics and algorithms (brief) intro to probability Basic notations Random variable - referring to an element / event

More information

Introduction to" Protein Structure

Introduction 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 information

CHAPTER 29 HW: AMINO ACIDS + PROTEINS

CHAPTER 29 HW: AMINO ACIDS + PROTEINS CAPTER 29 W: AMI ACIDS + PRTEIS For all problems, consult the table of 20 Amino Acids provided in lecture if an amino acid structure is needed; these will be given on exams. Use natural amino acids (L)

More information

Protein structure (and biomolecular structure more generally) CS/CME/BioE/Biophys/BMI 279 Sept. 28 and Oct. 3, 2017 Ron Dror

Protein structure (and biomolecular structure more generally) CS/CME/BioE/Biophys/BMI 279 Sept. 28 and Oct. 3, 2017 Ron Dror Protein structure (and biomolecular structure more generally) CS/CME/BioE/Biophys/BMI 279 Sept. 28 and Oct. 3, 2017 Ron Dror Please interrupt if you have questions, and especially if you re confused! Assignment

More information

Free Radical-Initiated Unfolding of Peptide Secondary Structure Elements

Free Radical-Initiated Unfolding of Peptide Secondary Structure Elements Free Radical-Initiated Unfolding of Peptide Secondary Structure Elements Thesis of the Ph.D. Dissertation by Michael C. Owen, M.Sc. Department of Chemical Informatics Faculty of Education University of

More information

The Dominant Interaction Between Peptide and Urea is Electrostatic in Nature: A Molecular Dynamics Simulation Study

The 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 information

Large Macromolecular Assemblies (LMAs) Large Macromolecular Assemblies (LMAs) Large Macromolecular Assemblies (LMAs)

Large Macromolecular Assemblies (LMAs) Large Macromolecular Assemblies (LMAs) Large Macromolecular Assemblies (LMAs) Large Macromolecular Assemblies (LMAs) : an Automated Computational Tool for Identifying Conserved Domains in the Structures of Macromolecular Assemblies Large macromolecular assemblies (such as ribosomes,

More information

THE TANGO ALGORITHM: SECONDARY STRUCTURE PROPENSITIES, STATISTICAL MECHANICS APPROXIMATION

THE 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 information

Computer simulations of protein folding with a small number of distance restraints

Computer 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 information

1) NMR is a method of chemical analysis. (Who uses NMR in this way?) 2) NMR is used as a method for medical imaging. (called MRI )

1) NMR is a method of chemical analysis. (Who uses NMR in this way?) 2) NMR is used as a method for medical imaging. (called MRI ) Uses of NMR: 1) NMR is a method of chemical analysis. (Who uses NMR in this way?) 2) NMR is used as a method for medical imaging. (called MRI ) 3) NMR is used as a method for determining of protein, DNA,

More information

NMR-Structure determination with the program CNS

NMR-Structure determination with the program CNS NMR-Structure determination with the program CNS Blockkurs 2013 Exercise 11.10.2013, room Mango? 1 NMR-Structure determination - Overview Amino acid sequence Topology file nef_seq.mtf loop cns_mtf_atom.id

More information

Protein dynamics. Folding/unfolding dynamics. Fluctuations near the folded state

Protein dynamics. Folding/unfolding dynamics. Fluctuations near the folded state Protein dynamics Folding/unfolding dynamics Passage over one or more energy barriers Transitions between infinitely many conformations B. Ozkan, K.A. Dill & I. Bahar, Protein Sci. 11, 1958-1970, 2002 Fluctuations

More information

CS 273 Prof. Serafim Batzoglou Prof. Jean-Claude Latombe Spring Lecture 12 : Energy maintenance (1) Lecturer: Prof. J.C.

CS 273 Prof. Serafim Batzoglou Prof. Jean-Claude Latombe Spring Lecture 12 : Energy maintenance (1) Lecturer: Prof. J.C. CS 273 Prof. Serafim Batzoglou Prof. Jean-Claude Latombe Spring 2006 Lecture 12 : Energy maintenance (1) Lecturer: Prof. J.C. Latombe Scribe: Neda Nategh How do you update the energy function during the

More information

Quiz 2 Morphology of Complex Materials

Quiz 2 Morphology of Complex Materials 071003 Quiz 2 Morphology of Complex Materials 1) Explain the following terms: (for states comment on biological activity and relative size of the structure) a) Native State b) Unfolded State c) Denatured

More information

BCMP 201 Protein biochemistry

BCMP 201 Protein biochemistry BCMP 201 Protein biochemistry BCMP 201 Protein biochemistry with emphasis on the interrelated roles of protein structure, catalytic activity, and macromolecular interactions in biological processes. The

More information

Discrimination of Near-Native Protein Structures From Misfolded Models by Empirical Free Energy Functions

Discrimination of Near-Native Protein Structures From Misfolded Models by Empirical Free Energy Functions PROTEINS: Structure, Function, and Genetics 41:518 534 (2000) Discrimination of Near-Native Protein Structures From Misfolded Models by Empirical Free Energy Functions David W. Gatchell, Sheldon Dennis,

More information

Protein 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. 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 information

Protein Structures. 11/19/2002 Lecture 24 1

Protein Structures. 11/19/2002 Lecture 24 1 Protein Structures 11/19/2002 Lecture 24 1 All 3 figures are cartoons of an amino acid residue. 11/19/2002 Lecture 24 2 Peptide bonds in chains of residues 11/19/2002 Lecture 24 3 Angles φ and ψ in the

More information

Cecilia Clementi s research group.

Cecilia Clementi s research group. Cecilia Clementi s research group http://leonardo.rice.edu/~cecilia/research/ Proteins don t have a folding problem it s we humans that do! Cartoons by Larry Gonick In principle, the laws of physics completely

More information

What is Protein Design?

What is Protein Design? Protein Design What is Protein Design? Given a fixed backbone, find the optimal sequence. Given a fixed backbone and native sequence, redesign a subset of positions (e.g. in the active site). What does

More information

Modeling Biological Systems Opportunities for Computer Scientists

Modeling Biological Systems Opportunities for Computer Scientists Modeling Biological Systems Opportunities for Computer Scientists Filip Jagodzinski RBO Tutorial Series 25 June 2007 Computer Science Robotics & Biology Laboratory Protein: πρώτα, "prota, of Primary Importance

More information

Lecture 26: Polymers: DNA Packing and Protein folding 26.1 Problem Set 4 due today. Reading for Lectures 22 24: PKT Chapter 8 [ ].

Lecture 26: Polymers: DNA Packing and Protein folding 26.1 Problem Set 4 due today. Reading for Lectures 22 24: PKT Chapter 8 [ ]. Lecture 26: Polymers: DA Packing and Protein folding 26.1 Problem Set 4 due today. eading for Lectures 22 24: PKT hapter 8 DA Packing for Eukaryotes: The packing problem for the larger eukaryotic genomes

More information

STRUCTURAL BIOINFORMATICS I. Fall 2015

STRUCTURAL BIOINFORMATICS I. Fall 2015 STRUCTURAL BIOINFORMATICS I Fall 2015 Info Course Number - Classification: Biology 5411 Class Schedule: Monday 5:30-7:50 PM, SERC Room 456 (4 th floor) Instructors: Vincenzo Carnevale - SERC, Room 704C;

More information