The computational Core of Molecular Modeling. (The What? Why? And How? )
|
|
- Maria Marsha Watkins
- 5 years ago
- Views:
Transcription
1 The computational Core of Molecular Modeling. (The What? Why? And How? ) CS Seminar. Feb 09. Alexey Onufriev, Dept of Computer Science; Dept. of Physics and the GBCB program.
2 Thanks to: Faculty co-authors: T.M. Murali, M. Prisant, L. Heath, C. Simmerling Student co-authors: J. Gordon, J. Myers, A. Fenley, J. Ruscio, R. Anandakrishnan, D. Kumar, M. Shukla, V. Sojia Sponsors: NIH, VT. Support: System-X team, N. Polys (myoglobin graphics)
3 Will focus on the computational side of solving problems in the following areas: a. Rational Drug Design b. From atomic motion to biological function c. The grand challenge of computational science: the Protein Folding Problem.
4 The emergence of in virtuo Science. in vivo in vitro in virtuo
5 Biological function = f( 3D molecular structure ) A T G C DNA sequence Protein structure Bilogical function How can we predict/understand/modify function if we know the structure? Can we predict the structure? Key challenges: Biomolecular structures are complex (e.g. compared to crystal solids). Biology works on many time scales. Experiments can only go so far. A solution: Computational methods. Model movements of individual atoms.
6 The paradigm: All things are made of atoms, and everything that living things do can be understood in terms of the wigglings and jigglings of atoms R. Feynmann Suggests the approach: model what nature does, i.e. let the molecule evolve with time according to underlying physics laws.
7 A protein on a surface. Atomic resolution
8 Theme 1. Rational (structurebased) design of new medicines: Picture: Design of a HIV protease inhibitor. Hornak, V.; Okur, A., Rizzo, R. and Simmerling, C., HIV-1 protease flaps spontaneously open and reclose in molecular dynamics simulations, Proc. Nat. Acad. Sci. USA, 103: (2006)
9 Example: rational drug design. If you block the enzymes function you kill the virus. e.g:viral protease (chops up proteins) Drug agent
10 Example of successful computer-aided (rational) drug design: One of the drugs that helped slow down the AIDS epidemic (part of anti-retro viral cocktail). The drug blocks the function of a key viral protein. To design the drug, one needs a precise 3D structure of that protein.
11 A computational challenge: Need high quality, complete protein structures. But the experiment (Xray) does not ``see hydrogen atoms. These have to be placed computationally.
12 Combinatorial explosion problem: A molecule with N sites, each of which can be occupied by an H+ atom (or be empty), has 2 N possible states. X 1 X W 23 X 3 All of the possible charge (protonation) states must be taken into account. For a typical protein with 100 ionizable amino-acids this means ~ possible variants to consider! And we need to find the minimum energy state among them (for the experts: what we really need is, of course, the partition function Z, see below. Which is an even more computationally demanding job) N ΔG k (ph) = Σ x i k (kt ln10 ph - Δ E i calc ) + 1/2ΣΣx ik x j k W ij Free energy i X i = 1 or 0 (occupied or empty in state k) k protonation state (out of 2 N ) N N i j trouble Matrix of site-site interactions. Z = Σ all states exp(-δg k /kt)
13 Beating the combinatorics: the clustering approach. Example: a protein with N = 6 ionizable groups. Total number of protonation states = 2 N = 64 Every site interacts with every other site, a complete graph Cluster 1, N=3 Cluster 2, N=3 Solution: neglect the ``weak edges. Cluster the strong ones. After clustering: Total number of states = = 16 < 64.
14 A web-server that adds hydrogens to molecular structures. Launched by Onufriev s group in June ~1000 registered users since.
15 Intrigued? Suggested reading: The Many roles of computation in drug discovery, W. Jorgensen, Sceince 303, 1813 (2004). + references therein.
16 THEME II The protein folding challenge. Nature does it all the time. Can we? Amino-acid sequence translated genetic code. MET ALA ALA ASP GLU GLU--. How? Experiment: amino acid sequence uniquely determines proteinʼs 3D shape (ground state). Why bother: proteinʼs shape determines its biological function.
17 Protein Structure in 3 steps. Step 1. Two amino-acids together (di-peptide) Peptide bond Amino-acid #1 Amino-acid #2
18 Protein Structure in 3 steps. Step 2: Most flexible degrees of freedom:
19 A protein is simply a chain of amino-acids: φ 1 φ 2 φ 3 φ 4 Each configuration {Φ1, Φ2, Φκ} has some energy. The folded (biologically functional) protein has the lowest possible energy - global minimum. So just find this conformation by some kind of a minimization algorithm what s he big deal?
20 The magnitude of the protein folding challenge: Enormous number of the possible conformations of the polypeptide chain φ 1 φ 2 φ 3 φ 4 A typical protein is a chain of ~ 100 mino acids. Assume that each amino acid can take up only 10 conformations (vast underestimation) Total number of possible conformations: Suppose each energy estimate is just 1 float point operation. Suppose you have a Penta-Flop supercomputer. An exhaustive search for the global minimum would take seconds ~ 3*10 78 years. Age of the Universe ~ 2*10 10 years.
21 2 Free energy 1 3 Folding coordinate Adopted from Ken Dill s web site at UCSF Finding a global minimum in a multidimensional case is easy only when the landscape is smooth. No matter where you start (1, 2 or 3), you quickly end up at the bottom - - the Native (N), functional state of the protein.
22 Adopted from Ken Dill s web site at UCSF Realistic landscapes are much more complex, with multiple local minima folding traps. Proteins trapped in those minima may lead to disease, such as Altzheimer s
23 Adopted from Dobson, NATURE 426,
24 Since minimization won t work, choose an alternative. Do what Nature does: just let it fold on its own, at normal temperature. Method: Molecular Dynamics
25 Principles of Molecular Dynamics (MD): Y Each atom moves by Newton s 2 nd Law: F = ma F = de/dr System s energy - Bond spring + x E = Kr 2 Bond stretching + A/r 12 B/r 6 VDW interaction + Q 1 Q 2 /r Electrostatic forces +
26 Now we have positions of all atoms as a function of time. Can compute statistical averages, fluctuations; Analyze side chain movements, Cavity dynamics, Domain motion, Etc.
27 Molecular Dynamics: PRICIPLE: Given positions of each atom x(t) at time t, its position at next time-step t + Δt is given by: force x(t + Δt) x(t) + v(t) Δt + ½*F/m * (Δt) 2 Key parameter: integration time step Δt. Controls accuracy and speed of numerical integration routines. Smaller Δt more accurate, but need more steps. How many needed to simulate biology? How many can one afford?
28 As a result, we can not quite get into the biological time scales. Currently accessible times biology H-C bond vibration Protein folding Characteristic time scales [sec] Time-step, Δt For stability, Δt must be at least an order of magnitude less than the fastest motion, i.e Δt ~ s. Example: to simulate folding of the fastest folding protein, at least 10-6 /10-15 = 10 9 steps will be needed.
29 The bottleneck of the methodology: computation of long-range interactions. Electrostatic interactions fall of as inverse distance between atoms. Too strong to neglect. Need to account for all of them. Very expensive. Up to 99% of total cost for a protein.
30 Massive parallel machines help. Virginia Tech s supercomputer, System-X
31 The worst problem for parallel computations: Processor #1 Processor #2 X 1, F(on X 1 ) X 2, F(on X 2 ) Force acting on each atom depends upon positions of every other atom in the system. Processor #4 Processor #3 Computed coordinates have to be communicated between all processors at each step
32 Without approximations, computation of long-range electrostatic forces will cost you O(N 2 ), where the number of atoms N may be as large as N ~ Too expensive for large systems. Every atom interacts with every other atom, a complete graph A solution: combine charges (vertices) into groups, that is use coarse-graining. +3 After coarse-graining costs can be as low as Nlog(N) Works because macromolecules are naturally partitioned into hierarchical levels: atoms -> amino-acids -> proteins -> complexes..
33
34
35 Simulated Refolding pathway of the 46-residue protein. Molecular dynamics based on AMBER-7 Movie available at: NB: due to the absence of viscosity, folding occurs on much shorter time-scale than in an experiment.
36 Intrigued? Suggested articles: 1. Protein Folding and Misfolding, C. Dobson, Nature 426, 884 (2003). 2. Design of a Novel Globular Protein Fold with Atomic-level Accuracy, Kuhlman et al., Science, 302, 1364, (2003) + references therein.
37 Theme III: how does oxygen get inside the protein myoglobin? Trivia: it is myoglobin (with oxygen bound to it) that gives fresh meat its red color. The protein stores oxygen.
38 Myoglobin the hydrogen atom of structural biology Exactly what are the migration pathways of small non-polar ligands (O2, CO, NO) inside the protein? A 50 year old question. O 2
39 The heme group the heart of myoglobin. O2
40 Myoglobin: protein responsible for oxygen storage in all breathing organisms. O 2? How does Oxygen, or other small non-polar ligands get in to bind to Fe atom in its core? Fe
41 After 50 years of research many questions still remain. Example: one or many pathways? Single dominant pathway, HIS gate e.g. Olson et al. Multiple pathays, e.g. Elber and Karplus
42 Two complementary approaches: #1 Simulate explicit ligand (CO) trajectories via microseconds of room temperature, explicit solvent MD. #2 Compute dynamic volume fluctuations in the protein.
43 Approach #1: Straightforward MD (explicit solvent, PME, NPT, 300K, 7 µs total) 1. Simulation type #1. Model photodetachement of carbon monoxyde (CO) ligand. 20 trajectories, 90 ns each. The ligand starts at the docking site (Fe) and comes out, sometimes. 2. Simulation type #2. Model diffusion of the ligand from the outside. The ligand starts in the solvent, and sometimes diffuses all the way to the docking site inside. 48 trajectories, each 90 ns long.
44 Raw Result #1: a sample of simulated trajectories that show ligand migration paths multiple individual trajectories combined in one slide Photodetachement experiment. CO Initially at the binding site CO initially in the water
45 Summary of the MD trajectories: Pathway #1 Pathway #2 Concern: Have we explored ALL the pathways? Challenge: can we get the same pathways cheaper?
46 Approach #2. PATHFINDER dynamic free volume analysis. Run a separate, much shorter (15 ns) MD in which the ligand does not move. In each snapshot, identify cavities large enough to hold the ligand. Purely steric criterion. Combine the cavities over the entire trajectory, determine connectivity. timeline pathway
47 The challenge: protein shape is very complex, internal cavities are hiding within, also very complex shape. Connectivity to the outside is a non-trivia issue. In this case, it is critical to have a SIMPLE and robust algorithm, to have just ONE complexity to deal with, not a product of the two: (problem complexity) x (algorithm complexity).
48 Free volume analysis (PATHFINDER)
49 Pathways from volume fluctuations (PATHFINDER):
50 Combining the two approaches: Explicit trajectories (MD) Free volume fluctuations J. Ruscio et al. Atomic level computational identification of ligand migration pathways between solvent and binding site in myoglobin, Proc. Natl. Acad. Sci. USA, 105, 9204 (2008)
51 The end result: oxygen diffusion pathways in myoglobin: The pathways (blue) occur inside the rigid scaffold structures (helices, white).
52 The law of (bio/phys/chem)- computational science: speed easy hard genius 0 accuracy
53 Lunch time
Computational structural biology and bioinformatics
Computational structural biology and bioinformatics What is it all about? Why take it? What are we going to be doing? Organizational notes. Grades etc. Books. CS6104. Spring CS6104. 04. Spring Alexey 04.
More informationWhy Proteins Fold. How Proteins Fold? e - ΔG/kT. Protein Folding, Nonbonding Forces, and Free Energy
Why Proteins Fold Proteins are the action superheroes of the body. As enzymes, they make reactions go a million times faster. As versatile transport vehicles, they carry oxygen and antibodies to fight
More informationLecture 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 informationProtein 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 informationFondamenti di Chimica Farmaceutica. Computer Chemistry in Drug Research: Introduction
Fondamenti di Chimica Farmaceutica Computer Chemistry in Drug Research: Introduction Introduction Introduction Introduction Computer Chemistry in Drug Design Drug Discovery: Target identification Lead
More informationProteins 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 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 informationMolecular 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 informationLecture 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 informationProtein 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 informationThermodynamics. Entropy and its Applications. Lecture 11. NC State University
Thermodynamics Entropy and its Applications Lecture 11 NC State University System and surroundings Up to this point we have considered the system, but we have not concerned ourselves with the relationship
More informationAn introduction to Molecular Dynamics. EMBO, June 2016
An introduction to Molecular Dynamics EMBO, June 2016 What is MD? everything that living things do can be understood in terms of the jiggling and wiggling of atoms. The Feynman Lectures in Physics vol.
More informationThe 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 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 informationLimitations 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 informationEarly Stages of Drug Discovery in the Pharmaceutical Industry
Early Stages of Drug Discovery in the Pharmaceutical Industry Daniel Seeliger / Jan Kriegl, Discovery Research, Boehringer Ingelheim September 29, 2016 Historical Drug Discovery From Accidential Discovery
More informationMolecular dynamics simulation. CS/CME/BioE/Biophys/BMI 279 Oct. 5 and 10, 2017 Ron Dror
Molecular dynamics simulation CS/CME/BioE/Biophys/BMI 279 Oct. 5 and 10, 2017 Ron Dror 1 Outline Molecular dynamics (MD): The basic idea Equations of motion Key properties of MD simulations Sample applications
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 informationModeling Biomolecular Systems II. BME 540 David Sept
Modeling Biomolecular Systems II BME 540 David Sept Introduction Why do we perform simulations? What makes simulations possible? How do we perform simulations? What types of things/systems do we simulate?
More informationMolecular Interactions F14NMI. Lecture 4: worked answers to practice questions
Molecular Interactions F14NMI Lecture 4: worked answers to practice questions http://comp.chem.nottingham.ac.uk/teaching/f14nmi jonathan.hirst@nottingham.ac.uk (1) (a) Describe the Monte Carlo algorithm
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 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 informationSAM Teacher s Guide Protein Partnering and Function
SAM Teacher s Guide Protein Partnering and Function Overview Students explore protein molecules physical and chemical characteristics and learn that these unique characteristics enable other molecules
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 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 informationHomology modeling. Dinesh Gupta ICGEB, New Delhi 1/27/2010 5:59 PM
Homology modeling Dinesh Gupta ICGEB, New Delhi Protein structure prediction Methods: Homology (comparative) modelling Threading Ab-initio Protein Homology modeling Homology modeling is an extrapolation
More informationDISCRETE 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 informationProtein Structure. W. M. Grogan, Ph.D. OBJECTIVES
Protein Structure W. M. Grogan, Ph.D. OBJECTIVES 1. Describe the structure and characteristic properties of typical proteins. 2. List and describe the four levels of structure found in proteins. 3. Relate
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 informationLecture 34 Protein Unfolding Thermodynamics
Physical Principles in Biology Biology 3550 Fall 2018 Lecture 34 Protein Unfolding Thermodynamics Wednesday, 21 November c David P. Goldenberg University of Utah goldenberg@biology.utah.edu Clicker Question
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 informationThe Molecular Dynamics Method
The Molecular Dynamics Method Thermal motion of a lipid bilayer Water permeation through channels Selective sugar transport Potential Energy (hyper)surface What is Force? Energy U(x) F = d dx U(x) Conformation
More informationDr. Sander B. Nabuurs. Computational Drug Discovery group Center for Molecular and Biomolecular Informatics Radboud University Medical Centre
Dr. Sander B. Nabuurs Computational Drug Discovery group Center for Molecular and Biomolecular Informatics Radboud University Medical Centre The road to new drugs. How to find new hits? High Throughput
More informationFree energy simulations
Free energy simulations Marcus Elstner and Tomáš Kubař January 14, 2013 Motivation a physical quantity that is of most interest in chemistry? free energies Helmholtz F or Gibbs G holy grail of computational
More informationWhy Proteins Fold? (Parts of this presentation are based on work of Ashok Kolaskar) CS490B: Introduction to Bioinformatics Mar.
Why Proteins Fold? (Parts of this presentation are based on work of Ashok Kolaskar) CS490B: Introduction to Bioinformatics Mar. 25, 2002 Molecular Dynamics: Introduction At physiological conditions, the
More informationEnergy Landscapes and Accelerated Molecular- Dynamical Techniques for the Study of Protein Folding
Energy Landscapes and Accelerated Molecular- Dynamical Techniques for the Study of Protein Folding John K. Prentice Boulder, CO BioMed Seminar University of New Mexico Physics and Astronomy Department
More informationPROTEIN-PROTEIN DOCKING REFINEMENT USING RESTRAINT MOLECULAR DYNAMICS SIMULATIONS
TASKQUARTERLYvol.20,No4,2016,pp.353 360 PROTEIN-PROTEIN DOCKING REFINEMENT USING RESTRAINT MOLECULAR DYNAMICS SIMULATIONS MARTIN ZACHARIAS Physics Department T38, Technical University of Munich James-Franck-Str.
More informationEvolutionary design of energy functions for protein structure prediction
Evolutionary design of energy functions for protein structure prediction Natalio Krasnogor nxk@ cs. nott. ac. uk Paweł Widera, Jonathan Garibaldi 7th Annual HUMIES Awards 2010-07-09 Protein structure prediction
More informationProtein Folding. I. Characteristics of proteins. C α
I. Characteristics of proteins Protein Folding 1. Proteins are one of the most important molecules of life. They perform numerous functions, from storing oxygen in tissues or transporting it in a blood
More informationWhy study protein dynamics?
Why study protein dynamics? Protein flexibility is crucial for function. One average structure is not enough. Proteins constantly sample configurational space. Transport - binding and moving molecules
More informationShort 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 informationALL 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 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 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 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 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 informationStructural biology and drug design: An overview
Structural biology and drug design: An overview livier Taboureau Assitant professor Chemoinformatics group-cbs-dtu otab@cbs.dtu.dk Drug discovery Drug and drug design A drug is a key molecule involved
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 informationEnhancing Specificity in the Janus Kinases: A Study on the Thienopyridine. JAK2 Selective Mechanism Combined Molecular Dynamics Simulation
Electronic Supplementary Material (ESI) for Molecular BioSystems. This journal is The Royal Society of Chemistry 2015 Supporting Information Enhancing Specificity in the Janus Kinases: A Study on the Thienopyridine
More informationStructure Investigation of Fam20C, a Golgi Casein Kinase
Structure Investigation of Fam20C, a Golgi Casein Kinase Sharon Grubner National Taiwan University, Dr. Jung-Hsin Lin University of California San Diego, Dr. Rommie Amaro Abstract This research project
More informationMolecular Dynamics. A very brief introduction
Molecular Dynamics A very brief introduction Sander Pronk Dept. of Theoretical Physics KTH Royal Institute of Technology & Science For Life Laboratory Stockholm, Sweden Why computer simulations? Two primary
More informationProtein Structure Analysis and Verification. Course S Basics for Biosystems of the Cell exercise work. Maija Nevala, BIO, 67485U 16.1.
Protein Structure Analysis and Verification Course S-114.2500 Basics for Biosystems of the Cell exercise work Maija Nevala, BIO, 67485U 16.1.2008 1. Preface When faced with an unknown protein, scientists
More informationLecture 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 informationDevelopment of Software Package for Determining Protein Titration Properties
Development of Software Package for Determining Protein Titration Properties By Kaila Bennett, Amitoj Chopra, Jesse Johnson, Enrico Sagullo Bioengineering 175-Senior Design University of California Riverside
More informationBIOCHEMISTRY 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 informationBiophysics II. Key points to be covered. Molecule and chemical bonding. Life: based on materials. Molecule and chemical bonding
Biophysics II Life: based on materials By A/Prof. Xiang Yang Liu Biophysics & Micro/nanostructures Lab Department of Physics, NUS Organelle, formed by a variety of molecules: protein, DNA, RNA, polysaccharides
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 informationSupplementary Figure S1. Urea-mediated buffering mechanism of H. pylori. Gastric urea is funneled to a cytoplasmic urease that is presumably attached
Supplementary Figure S1. Urea-mediated buffering mechanism of H. pylori. Gastric urea is funneled to a cytoplasmic urease that is presumably attached to HpUreI. Urea hydrolysis products 2NH 3 and 1CO 2
More informationProtein Simulations in Confined Environments
Critical Review Lecture Protein Simulations in Confined Environments Murat Cetinkaya 1, Jorge Sofo 2, Melik C. Demirel 1 1. 2. College of Engineering, Pennsylvania State University, University Park, 16802,
More informationIn silico pharmacology for drug discovery
In silico pharmacology for drug discovery In silico drug design In silico methods can contribute to drug targets identification through application of bionformatics tools. Currently, the application of
More informationFree Energy. because H is negative doesn't mean that G will be negative and just because S is positive doesn't mean that G will be negative.
Biochemistry 462a Bioenergetics Reading - Lehninger Principles, Chapter 14, pp. 485-512 Practice problems - Chapter 14: 2-8, 10, 12, 13; Physical Chemistry extra problems, free energy problems Free Energy
More informationMolecular Dynamics. Molecules in motion
Molecular Dynamics Molecules in motion 1 Molecules in mo1on Molecules are not sta1c, but move all the 1me Source: h9p://en.wikipedia.org/wiki/kine1c_theory 2 Gasses, liquids and solids Gasses, liquids
More informationNB-DNJ/GCase-pH 7.4 NB-DNJ+/GCase-pH 7.4 NB-DNJ+/GCase-pH 4.5
SUPPLEMENTARY TABLES Suppl. Table 1. Protonation states at ph 7.4 and 4.5. Protonation states of titratable residues in GCase at ph 7.4 and 4.5. Histidine: HID, H at δ-nitrogen; HIE, H at ε-nitrogen; HIP,
More informationAdvanced Molecular Dynamics
Advanced Molecular Dynamics Introduction May 2, 2017 Who am I? I am an associate professor at Theoretical Physics Topics I work on: Algorithms for (parallel) molecular simulations including GPU acceleration
More informationDana Alsulaibi. Jaleel G.Sweis. Mamoon Ahram
15 Dana Alsulaibi Jaleel G.Sweis Mamoon Ahram Revision of last lectures: Proteins have four levels of structures. Primary,secondary, tertiary and quaternary. Primary structure is the order of amino acids
More informationEnzyme Catalysis & Biotechnology
L28-1 Enzyme Catalysis & Biotechnology Bovine Pancreatic RNase A Biochemistry, Life, and all that L28-2 A brief word about biochemistry traditionally, chemical engineers used organic and inorganic chemistry
More informationIntroduction 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 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 informationLecture 11: Enzymes: Kinetics [PDF] Reading: Berg, Tymoczko & Stryer, Chapter 8, pp
Lecture 11: Enzymes: Kinetics [PDF] Reading: Berg, Tymoczko & Stryer, Chapter 8, pp. 216-225 Updated on: 2/4/07 at 9:00 pm Key Concepts Kinetics is the study of reaction rates. Study of enzyme kinetics
More informationF. Piazza Center for Molecular Biophysics and University of Orléans, France. Selected topic in Physical Biology. Lecture 1
Zhou Pei-Yuan Centre for Applied Mathematics, Tsinghua University November 2013 F. Piazza Center for Molecular Biophysics and University of Orléans, France Selected topic in Physical Biology Lecture 1
More informationThe Computational Microscope
The Computational Microscope Computational microscope views at atomic resolution... Rs SER RER C E M RER N GA L... how living cells maintain health and battle disease John Stone Our Microscope is Made
More informationThe Molecular Dynamics Simulation Process
The Molecular Dynamics Simulation Process For textbooks see: M.P. Allen and D.J. Tildesley. Computer Simulation of Liquids.Oxford University Press, New York, 1987. D. Frenkel and B. Smit. Understanding
More information2-4 Chemical Reactions and Enzymes
2-4 Chemical Reactions and Enzymes Living things, as you have seen, are made up of chemical compounds-some simple and some complex. But chemistry isn t just what life is made of-chemistry is also what
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 informationPrinciples and Applications of Molecular Dynamics Simulations with NAMD
Principles and Applications of Molecular Dynamics Simulations with NAMD Nov. 14, 2016 Computational Microscope NCSA supercomputer JC Gumbart Assistant Professor of Physics Georgia Institute of Technology
More informationOverview & Applications. T. Lezon Hands-on Workshop in Computational Biophysics Pittsburgh Supercomputing Center 04 June, 2015
Overview & Applications T. Lezon Hands-on Workshop in Computational Biophysics Pittsburgh Supercomputing Center 4 June, 215 Simulations still take time Bakan et al. Bioinformatics 211. Coarse-grained Elastic
More informationCrystal Structure Prediction using CRYSTALG program
Crystal Structure Prediction using CRYSTALG program Yelena Arnautova Baker Laboratory of Chemistry and Chemical Biology, Cornell University Problem of crystal structure prediction: - theoretical importance
More informationBiological Thermodynamics
Biological Thermodynamics Classical thermodynamics is the only physical theory of universal content concerning which I am convinced that, within the framework of applicability of its basic contents, will
More informationImproving Protein Function Prediction with Molecular Dynamics Simulations. Dariya Glazer Russ Altman
Improving Protein Function Prediction with Molecular Dynamics Simulations Dariya Glazer Russ Altman Motivation Sometimes the 3D structure doesn t score well for a known function. The experimental structure
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 informationIdentification of core amino acids stabilizing rhodopsin
Identification of core amino acids stabilizing rhodopsin A. J. Rader, Gulsum Anderson, Basak Isin, H. Gobind Khorana, Ivet Bahar, and Judith Klein- Seetharaman Anes Aref BBSI 2004 Outline Introduction
More informationModeling 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 informationFree energy, electrostatics, and the hydrophobic effect
Protein Physics 2016 Lecture 3, January 26 Free energy, electrostatics, and the hydrophobic effect Magnus Andersson magnus.andersson@scilifelab.se Theoretical & Computational Biophysics Recap Protein structure
More informationStructural Bioinformatics (C3210) Molecular Mechanics
Structural Bioinformatics (C3210) Molecular Mechanics How to Calculate Energies Calculation of molecular energies is of key importance in protein folding, molecular modelling etc. There are two main computational
More informationVMD: Visual Molecular Dynamics. Our Microscope is Made of...
VMD: Visual Molecular Dynamics Computational Microscope / Tool to Think amino acid tyrosine enzymatic control BPTI VMD tutorial traficking Ubiquitin case study http://www.ks.uiuc.edu/training/casestudies/
More informationMolecular Structure Prediction by Global Optimization
Molecular Structure Prediction by Global Optimization K.A. DILL Department of Pharmaceutical Chemistry, University of California at San Francisco, San Francisco, CA 94118 A.T. PHILLIPS Computer Science
More informationBiochemistry Prof. S. DasGupta Department of Chemistry Indian Institute of Technology Kharagpur. Lecture - 06 Protein Structure IV
Biochemistry Prof. S. DasGupta Department of Chemistry Indian Institute of Technology Kharagpur Lecture - 06 Protein Structure IV We complete our discussion on Protein Structures today. And just to recap
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 informationThe Molecular Dynamics Method
H-bond energy (kcal/mol) - 4.0 The Molecular Dynamics Method Fibronectin III_1, a mechanical protein that glues cells together in wound healing and in preventing tumor metastasis 0 ATPase, a molecular
More informationEnergy functions and their relationship to molecular conformation. CS/CME/BioE/Biophys/BMI 279 Oct. 3 and 5, 2017 Ron Dror
Energy functions and their relationship to molecular conformation CS/CME/BioE/Biophys/BMI 279 Oct. 3 and 5, 2017 Ron Dror Yesterday s Nobel Prize: single-particle cryoelectron microscopy 2 Outline Energy
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 informationMolecular 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 informationDistance Constraint Model; Donald J. Jacobs, University of North Carolina at Charlotte Page 1 of 11
Distance Constraint Model; Donald J. Jacobs, University of North Carolina at Charlotte Page 1 of 11 Taking the advice of Lord Kelvin, the Father of Thermodynamics, I describe the protein molecule and other
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 informationPhysical Models of Allostery: Allosteric Regulation in Capsid Assembly
Physical Models of Allostery: Allosteric Regulation in Capsid Assembly QCB Journal Club Prof. Sima Setayeshgar JB Holmes Nov. 2, 2017 Mechanisms of Allosteric Regulation From R.A. Laskowski, FEBS Letters,
More informationProtein 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 informationMolecular 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 informationCase Study: Protein Folding using Homotopy Methods
Case Study: Protein Folding using Homotopy Methods james robert white Since the late 1990s, the world of genomics has exploded. Beginning with the bacterium Haemophilus influenzae in 1995, researchers
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 informationBiochemistry 3100 Sample Problems Binding proteins, Kinetics & Catalysis
(1) Draw an approximate denaturation curve for a typical blood protein (eg myoglobin) as a function of ph. (2) Myoglobin is a simple, single subunit binding protein that has an oxygen storage function
More information