Free energy calculations with alchemlyb
|
|
- Jacob Tate
- 5 years ago
- Views:
Transcription
1 Free energy calculations with alchemlyb Oliver Beckstein Arizona State University SPIDAL Teleconference
2 Binding free energy measure of how strong a protein P and a ligand X stick together key quantity in quantitative mechanistic explanations of biological processes and in drug discovery free energy (i.e., averaged over all other degrees of freedom (such as solvent, protein motions, )) free energy difference exp[ A(T, V, N)/kT ] = dp 3N dx 3N e H(p,x)/kT ΔA = A bound A unbound ΔA = ln P + X ΔF PX dx3n exp[ U(x)/kT ] χ bound (x) dx 3N exp[ U(x)/kT ] χ unbound (x) H(p, x) = N i=1 p 2 i 2m i + U(x)
3 Alchemical free energy calculations force fields for H, molecular dynamics (MD) simulations for sampling free energy is a state function: generate non-physical paths between physical end states (bound/unbound) use stratification ( windows ) of path with parameter λ H(λ) = (1 λ)h bound + λh unbound, 0 λ 1 U Coulomb (x 1, x 2 ; λ) = 1 (1 λ)q 1 q 2 4πϵ 0 x 1 x 2 λ 3 λ 4 λ 2 λ 5 λ 6 λ 1 = 0 λ n = 1
4 Methods Free Energy Perturbation (FEP): Zwanzig FEP, BAR, MBAR (overlaps of distributions) ΔA = kt ln exp[ ΔU(x)/kT] 1, Thermodynamic Integration (TI): TI ΔA = 1 0 dλ H(λ) λ with ΔU(x) = U 1 (x) U 0 (x) exp( ΔA/kT ) = 1 + exp[(δu C)/kT ] exp[(δu C)/kT ] 1 1 exp( C/kT ) ΔA/kT = C/kT ln n 1 n 2 (MBAR is more complicated: uses overlaps between all windows) ΔU(x) = U λn+1 (x) U λn (x) λ FEP BAR TI C. Chipot and A. Pohorille, editors. Free energy calculations. Number 86 in Springer Series in Chemical Physics. Springer, Berlin, 2007.
5 Thermodynamic cycle for absolute binding free energy calculations appearance in solvent (no protein) decoupling from protein P + X ΔF PX MD simulations with multiple λ 5 50 λ-windows per free energy component 10 ns 250 ns per window D. L. Mobley, J. D. Chodera, and K. A. Dill. On the use of orientational restraints and symmetry corrections in alchemical free energy calculations. 125:084902, 2006.
6 Simulation output Per timestep FEP: ΔU ΔU λi,λ j j for window i TI: U/ λ But: different file formats for different codes (Gromacs, Amber, NAMD, ) solution: common interface via alchemlyb U(x; λ) λ λ i
7 alchemlyb O Beckstein (ASU), D Mobley (UC Irvine), M Shirts (U Colorado, Boulder) David Dotson (ASU), Dominik Wille (Freie Univ. Berlin) Silicon Therapeutics (STX) (Bryce Allen, Shuai Liu)
8 alchemlyb: basic idea data data data parse standard form (pandas DataFrame) preprocess standard form (pandas DataFrame) estimate ΔA ΔU λi,λ j j U(x; λ) λ λ i preprocessors slicing statistical inefficiency equilibrium detection estimators TI BAR MBAR
9 TI example
10
11
12 NapA elevator mechanism K305 D156 TM 5 T126:O TM 4 D157 Nature Struct. Mol. Biol., 23 (2016): ion binding in different states? dimer domain e. TM 5 D156 Inward-facing S4 A153:O [ need coordination image ] T126:O K305 TM 4 D157 Nature, 501 (2013): Na+ ion core domain
13 Absolute binding free energies: alchemistry Ion-in-solvent i. j. + alchemical binding pathway h. g. (analytic) f. equivalent e a. b. c. d. NapA complex Windowed alchemical free energy calculations (TI or MBAR) 150 ns 250 ns per lambda window (Coulomb/vdW decoupling) for windows ~8 µs (!) Position restraints (and analytical removal) Additional repulsive ion-ion potential to enforce one-ion occupancy (rigorously removed in calculation)
14 Amount of raw free energy data NapA FEP simulations (one free energy) windows time µs size GB total time µs total size (GB) VDW Coulomb repulsion restraint conformations: IF and OF (2) protonation states: 3 repeats: x2 (some) ~12 sets of simulations: ~2 TB (in ~130,000 files)
15 <latexit sha1_base64="0i1nre9o2xwd6b+g814wnznx2wk=">aaaczhicbvdrshtbfj2strvaa1t6jmhgkai2ysmf7ymgktrhbanckoa7k5tkchzmmlkrjct+l9/ig6/6fujnyx6memdgzlnnzp17equkpzi+q0qzsx8+fpr7pl/wzfhrunv55dybzalscqomu0zao5iamyrj4av1cgmi8ck5oizrf9fovdt6jeywoykmtoxlarskbvw0fyskgp/5g3cl3+mvrnmb8b/l1hlcqbecvfjbzcf/vmmajnpu/gjrdku1ub6pwd+sxotu2aqn3eptu2delqimocd7vio21mnbkrqki/l25tgcuiibtglvkklv5opvc/49kd3eny4ctxysvuziifv+lcbbmqin/etakb5xa2xu3+3kutumuivnqf1mctk8zjh3penbahqiccfdx7kyggnbie2pkrzuieabqu3ycqkfg1pm1xidzlvs3k7/rsenv2r7b5pg5tga22cbrmf22d47ziesyqs7zffsgt1whqpfadx69mynkpoevtafap0/fjy7ig==</latexit> Absolute binding free energies: alchemistry DG 0 i = DG 0 protein+ion,i DG 0 hydration NapA IF D156 D157 K305 ΔG 0 (kj/mol) 103±1 H 44.6±0.9 H H 24.2±13.0 (not binding ignore state) Δ G (kj/mol) Convergence of ion hydration calculation coulomb vdw repulsion simulation time used (ns) Δ G (kj/mol) simulation time used (ns) Δ G (kj/mol) ion-ion simulation time used (ns) Na + ion in CHARMM TIP3P water ΔGhydration (kj/mol) Coulomb 382.7±0.2 VdW 8.81±0.05 repulsion 0.258±0.005 restraint total 392.0±0.2
16 Convergence of protein ion calculations ion-ion Coulomb H χ 180º χ 90º all kj/mol restraint kj/mol D156 D157 K305 repulsion Very slow convergence of Coulomb Slow degrees of freedom (e.g., D157 χ1 dihedral) Time (ns) χ1 Time (ns) Na+ D157
17 Challenge data analysis is cumbersome and slow(ish), even with dask on a 6 core workstation (hours) on-demand analysis with all current data/while new data is coming in? run on HPC system (e.g. XSEDE PSC Bridges)?
Alchemical free energy calculations in OpenMM
Alchemical free energy calculations in OpenMM Lee-Ping Wang Stanford Department of Chemistry OpenMM Workshop, Stanford University September 7, 2012 Special thanks to: John Chodera, Morgan Lawrenz Outline
More informationLecture 19: Free Energies in Modern Computational Statistical Thermodynamics: FEP and Related Methods
Statistical Thermodynamics Lecture 19: Free Energies in Modern Computational Statistical Thermodynamics: FEP and Related Methods Dr. Ronald M. Levy ronlevy@temple.edu Free energy calculations Free energy
More informationFree energy calculations
Free energy calculations Berk Hess May 5, 2017 Why do free energy calculations? The free energy G gives the population of states: ( ) P 1 G = exp, G = G 2 G 1 P 2 k B T Since we mostly simulate in the
More informationFree energy calculations
Free energy calculations Jochen Hub & David van der Spoel Overview Free energies and Probabilities Thermodynamic cycles (Free energy perturbation (FEP)) Thermodynamic integration (TI) (Jarzynski equality
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 informationFree energy calculations using molecular dynamics simulations. Anna Johansson
Free energy calculations using molecular dynamics simulations Anna Johansson 2007-03-13 Outline Introduction to concepts Why is free energy important? Calculating free energy using MD Thermodynamical Integration
More informationBiomolecular modeling. Theoretical Chemistry, TU Braunschweig (Dated: December 10, 2010)
Biomolecular modeling Marcus Elstner and Tomáš Kubař Theoretical Chemistry, TU Braunschweig (Dated: December 10, 2010) IX. FREE ENERGY SIMULATIONS When searching for a physical quantity that is of most
More informationw REXAMD: A Hamiltonian Replica Exchange Approach to Improve Free Energy Calculations for Systems with Kinetically Trapped Conformations
pubs.acs.org/jctc Downloaded via 148.251.232.83 on March 8, 2019 at 14:33:02 (UTC). See https://pubs.acs.org/sharingguidelines for options on how to legitimately share published articles. w REXAMD: A Hamiltonian
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 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 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 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 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 informationConformational Free-Energy Differences by Confinement Simulations
Conformational Free-Energy Differences by Confinement Simulations Marco Cecchini Laboratoire d Ingeniérie des Fonctions Moléculaires (ISIS) UMR 76 - Université de Strasbourg mcecchini@unistra.fr HPC days
More informationT H E J O U R N A L O F G E N E R A L P H Y S I O L O G Y. jgp
S u p p l e m e n ta l m at e r i a l jgp Lee et al., http://www.jgp.org/cgi/content/full/jgp.201411219/dc1 T H E J O U R N A L O F G E N E R A L P H Y S I O L O G Y S u p p l e m e n ta l D I S C U S
More informationStatistical Mechanics for Proteins
The Partition Function From Q all relevant thermodynamic properties can be obtained by differentiation of the free energy F: = kt q p E q pd d h T V Q ), ( exp 1! 1 ),, ( 3 3 3 ),, ( ln ),, ( T V Q kt
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 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 informationEstimating Atomic Contributions to Hydration and Binding Using Free Energy Perturbation
This is an open access article published under a Creative Commons Attribution (CC-BY) License, which permits unrestricted use, distribution and reproduction in any medium, provided the author and source
More informationNWChem: Molecular Dynamics and QM/MM
NWChem: Molecular Dynamics and QM/MM Molecular Dynamics Functionality Target systems: General features: biomolecules (proteins, DNA/RNA, biomembranes) energy evaluation (SP) energy minimization (EM) molecular
More informationComputing free energies with PLUMED 2.0
Computing free energies with PLUMED 2.0 Davide Branduardi Formerly at MPI Biophysics, Frankfurt a.m. TyOutline of the talk Relevance of free energy computation Free-energy from histograms: issues Tackling
More informationThermodynamic Integration with Enhanced Sampling (TIES)
Thermodynamic Integration with Enhanced Sampling (TIES) A. P. Bhati, S. Wan, D. W. Wright and P. V. Coveney agastya.bhati.14@ucl.ac.uk Centre for Computational Science Department of Chemistry University
More informationVariational Implicit Solvation of Biomolecules: From Theory to Numerical Computations
Variational Implicit Solvation of Biomolecules: From Theory to Numerical Computations Bo Li Department of Mathematics and Center for Theoretical Biological Physics UC San Diego CECAM Workshop: New Perspectives
More informationHow is molecular dynamics being used in life sciences? Davide Branduardi
How is molecular dynamics being used in life sciences? Davide Branduardi davide.branduardi@schrodinger.com Exploring molecular processes with MD Drug discovery and design Protein-protein interactions Protein-DNA
More informationMolecular simulation and structure prediction using CHARMM and the MMTSB Tool Set Free Energy Methods
Molecular simulation and structure prediction using CHARMM and the MMTSB Tool Set Free Energy Methods Charles L. Brooks III MMTSB/CTBP 2006 Summer Workshop CHARMM Simulations The flow of data and information
More informationFree-energy calculations
Measuring free-energy differences using computer simulations Theoretical and Computational Biophysics Group University of Illinois, Urbana Champaign and Équipe de dynamique des assemblages membranaires,
More informationFree Energy Simulation Methods
Free Energy Simulation Methods Free energy simulation methods Many methods have been developed to compute (relative) free energies on the basis of statistical mechanics Free energy perturbation Thermodynamic
More informationLecture 20: Binding Free Energy Calculations
Statistical Thermodynamics ecture 20: Binding Free Energy Calculations Dr. onald M. evy ronlevy@temple.edu Outline focus on the binding problem of protein-ligand systems. Statistical theory of non-covalent
More informationBinding Free Energy Calculations for Lead Optimization: Assessment of Their Accuracy in an Industrial Drug Design Context
pubs.acs.org/jctc Binding Free Energy Calculations for Lead Optimization: Assessment of Their Accuracy in an Industrial Drug Design Context Nadine Homeyer, Friederike Stoll, Alexander Hillisch, and Holger
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 informationComputational Studies of the Photoreceptor Rhodopsin. Scott E. Feller Wabash College
Computational Studies of the Photoreceptor Rhodopsin Scott E. Feller Wabash College Rhodopsin Photocycle Dark-adapted Rhodopsin hn Isomerize retinal Photorhodopsin ~200 fs Bathorhodopsin Meta-II ms timescale
More informationA general overview of. Free energy methods. Comer et al. (2014) JPCB 119:1129. James C. (JC) Gumbart Georgia Institute of Technology, Atlanta
A general overview of Free energy methods Comer et al. (2014) JPCB 119:1129. James C. (JC) Gumbart Georgia Institute of Technology, Atlanta Computational Biophysics Workshop Georgia Tech Nov. 16 2016 Outline
More information4 th Advanced in silico Drug Design KFC/ADD Molecular Modelling Intro. Karel Berka, Ph.D.
4 th Advanced in silico Drug Design KFC/ADD Molecular Modelling Intro Karel Berka, Ph.D. UP Olomouc, 21.1.-25.1. 2019 Motto A theory is something nobody believes, except the person who made it An experiment
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 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 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 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 informationPotential Energy (hyper)surface
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 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 informationLarge Scale FEP on Congeneric Ligand Series Have Practical Free Energy Calculations arrived at Last?
Large Scale FEP on Congeneric Ligand Series Have Practical Free Energy Calculations arrived at Last? Thomas Steinbrecher, Teng Lin, Lingle Wang, Goran Krilov, Robert Abel, Woody Sherman, Richard Friesner
More informationEffective dynamics for the (overdamped) Langevin equation
Effective dynamics for the (overdamped) Langevin equation Frédéric Legoll ENPC and INRIA joint work with T. Lelièvre (ENPC and INRIA) Enumath conference, MS Numerical methods for molecular dynamics EnuMath
More informationAvailable online at ScienceDirect. Physics Procedia 57 (2014 ) 7 15
Available online at www.sciencedirect.com ScienceDirect Physics Procedia 57 (2014 ) 7 15 27th Annual CSP Workshops on Recent Developments in Computer Simulation Studies in Condensed Matter Physics, CSP
More informationProblem solving steps
Problem solving steps Determine the reaction Write the (balanced) equation ΔG K v Write the equilibrium constant v Find the equilibrium constant using v If necessary, solve for components K K = [ p ] ν
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 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 informationExam 3 Solutions. ClO g. At 200 K and a total pressure of 1.0 bar, the partial pressure ratio for the chlorine-containing compounds is p ClO2
Chemistry 360 Dr. Jean M. Standard Fall 2016 Name KEY Exam 3 Solutions 1.) (14 points) Consider the gas phase decomposition of chlorine dioxide, ClO 2, ClO 2 ( g) ClO ( g) + O ( g). At 200 K and a total
More informationDon t forget to bring your MD tutorial. Potential Energy (hyper)surface
Don t forget to bring your MD tutorial Lab session starts at 1pm You will have to finish an MD/SMD exercise on α-conotoxin in oxidized and reduced forms Potential Energy (hyper)surface What is Force? Energy
More informationThe Second Law of Thermodynamics (Chapter 4)
The Second Law of Thermodynamics (Chapter 4) First Law: Energy of universe is constant: ΔE system = - ΔE surroundings Second Law: New variable, S, entropy. Changes in S, ΔS, tell us which processes made
More information3. Solutions W = N!/(N A!N B!) (3.1) Using Stirling s approximation ln(n!) = NlnN N: ΔS mix = k (N A lnn + N B lnn N A lnn A N B lnn B ) (3.
3. Solutions Many biological processes occur between molecules in aqueous solution. In addition, many protein and nucleic acid molecules adopt three-dimensional structure ( fold ) in aqueous solution.
More informationComputing Relative Free Energies of Solvation Using Single Reference Thermodynamic Integration Augmented with Hamiltonian Replica Exchange
J. Chem. Theory Comput. 2010, 6, 3427 3441 3427 Computing Relative Free Energies of Solvation Using Single Reference Thermodynamic Integration Augmented with Hamiltonian Replica Exchange Ilja V. Khavrutskii*
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 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 informationLots of thanks for many reasons!
Lots of thanks for many reasons! Structure and function of glutamate transporters from free energy simulations Turgut Baştuğ TOBB ETU Transporter structures Transporters have larger structures, which
More informationProtein:ligand standard binding free energies: A tutorial for alchemical and geometrical transformations
School of Physics Georgia Institute of Technology Department of Biochemistry and Molecular Biology Gordon Center for Integrative Science The University of Chicago Centre National de la Recherche Scientifique
More informationBasic Ingredients of Free Energy Calculations: A Review
Feature Article Basic Ingredients of Free Energy Calculations: A Review CLARA D. CHRIST, 1 ALAN E. MARK, 2 WILFRED F. van GUNSTEREN 1 1 Laboratory of Physical Chemistry, Swiss Federal Institute of Technology,
More informationSupporting Information Quantitative characterization of the binding and unbinding of millimolar drug fragments with molecular dynamics simulations
Supporting Information Quantitative characterization of the binding and unbinding of millimolar drug fragments with molecular dynamics simulations Simulation Details All binding simulations were run on
More informationTowards accurate calculations of Zn 2+ binding free energies in zinc finger proteins
Towards accurate calculations of Zn 2+ binding free energies in zinc finger proteins UNDERGRADUATE HONORS RESEARCH THESIS Presented in Partial Fulfillment of the Requirements for the Bachelor of Science
More informationAdvanced in silico drug design
Advanced in silico drug design RNDr. Martin Lepšík, Ph.D. Lecture: Advanced scoring Palacky University, Olomouc 2016 1 Outline 1. Scoring Definition, Types 2. Physics-based Scoring: Master Equation Terms
More informationGromacs Workshop Spring CSC
Gromacs Workshop Spring 2007 @ CSC Erik Lindahl Center for Biomembrane Research Stockholm University, Sweden David van der Spoel Dept. Cell & Molecular Biology Uppsala University, Sweden Berk Hess Max-Planck-Institut
More informationComputational Study of Gleevec and G6G Reveals Molecular Determinants of Kinase Inhibitor Selectivity
This is an open access article published under an ACS AuthorChoice License, which permits copying and redistribution of the article or any adaptations for non-commercial purposes. pubs.acs.org/jacs Computational
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 informationNew Method for Calculating the Absolute Free Energy of Binding: The Effect of a Mobile Loop on the Avidin/Biotin Complex
168 J. Phys. Chem. B 2011, 115, 168 175 New Method for Calculating the Absolute Free Energy of Binding: The Effect of a Mobile Loop on the Avidin/Biotin Complex Ignacio J. General, Ralitsa Dragomirova,
More informationSampling the free energy surfaces of collective variables
Sampling the free energy surfaces of collective variables Jérôme Hénin Enhanced Sampling and Free-Energy Calculations Urbana, 12 September 2018 Please interrupt! struct bioinform phys chem theoretical
More informationSimulating biomolecular function from motions across multiple scales (I) Peter J. Bond (BII)
Simulating biomolecular function from motions across multiple scales (I) Peter J. Bond (BII) peterjb@bii.a-star.edu.sg Structural Biology: Why the Need for Simulation? 2017 year 1972 0 RCSB PDB: RCSB Protein
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 informationA Cavity Corrected 3D-RISM Functional for Accurate Solvation Free Energies
pubs.acs.org/jctc Terms of Use A Cavity Corrected 3D-RISM Functional for Accurate Solvation Free Energies Jean-Franc ois Truchon,*, B. Montgomery Pettitt, and Paul Labute Chemical Computing Group Inc.,
More informationMolecular modeling with InsightII
Molecular modeling with InsightII Yuk Sham Computational Biology/Biochemistry Consultant Phone: (612) 624 7427 (Walter Library) Phone: (612) 624 0783 (VWL) Email: shamy@msi.umn.edu How to run InsightII
More informationBiomolecular modeling I
2015, December 15 Biomolecular simulation Elementary body atom Each atom x, y, z coordinates A protein is a set of coordinates. (Gromacs, A. P. Heiner) Usually one molecule/complex of interest (e.g. protein,
More informationMolecular Dynamics Simulation of a Nanoconfined Water Film
Molecular Dynamics Simulation of a Nanoconfined Water Film Kyle Lindquist, Shu-Han Chao May 7, 2013 1 Introduction The behavior of water confined in nano-scale environment is of interest in many applications.
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 informationA Large-Scale Test of Free-Energy Simulation Estimates of Protein-Ligand Binding Affinities.
A Large-Scale Test of Free-Energy Simulation Estimates of Protein-Ligand Binding Affinities. Mikulskis, Paulius; Genheden, Samuel; Ryde, Ulf Published in: DOI: 10.1021/ci5004027 2014 Link to publication
More informationschematic diagram; EGF binding, dimerization, phosphorylation, Grb2 binding, etc.
Lecture 1: Noncovalent Biomolecular Interactions Bioengineering and Modeling of biological processes -e.g. tissue engineering, cancer, autoimmune disease Example: RTK signaling, e.g. EGFR Growth responses
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 informationTITAN: Two-dimensional lineshape analysis
TITAN: Two-dimensional lineshape analysis Chris Waudby Christodoulou Group c.waudby@ucl.ac.uk Andres Ramos Lisa Cabrita John Christodoulou Inhibition of fatty acid synthesis for treatment of tularemia
More informationAssessing the potential of atomistic molecular dynamics simulations to probe reversible. protein-protein recognition and binding
Supporting information for Assessing the potential of atomistic molecular dynamics simulations to probe reversible protein-protein recognition and binding Luciano A. Abriata and Matteo Dal Peraro Laboratory
More informationD3RGC2: free energy scoring by alchemical free energy implementation in SOMD
D3RGC2: free energy scoring by alchemical free energy implementation in SOMD Julien Michel EaStCHEM School of Chemistry, University of Edinburgh, United Kingdom D3R webinar 27/03/17 Acknowledgements: group
More informationMolecular Dynamics Simulations of the Mammalian Glutamate Transporter EAAT3
Molecular Dynamics Simulations of the Mammalian Glutamate Transporter EAAT3 Germano Heinzelmann, Serdar Kuyucak* School of Physics, University of Sydney, NSW, Australia Abstract Excitatory amino acid transporters
More informationChemical properties that affect binding of enzyme-inhibiting drugs to enzymes
Introduction Chemical properties that affect binding of enzyme-inhibiting drugs to enzymes The production of new drugs requires time for development and testing, and can result in large prohibitive costs
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 informationOrthogonal Space Sampling of Slow Environment Responses
IMA University of Minnesota 2015 Orthogonal Space Sampling of Slow Environment Responses Lianqing Zheng, Chao Lv, Dongsheng Wu, William Harris, Xubin Li, Erick Aitchison, and Wei Yang Institute of Molecular
More informationRate Law Summary. Rate Laws vary as a function of time
Rate Law Summary Measure the instantaneous rate of a reaction: this is a number with units of M/s! Measure the rate of loss of a reactant r... the rate of appearance of a product Repeat the experiment
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 informationCoupling Protein Dynamics with Enzyme Catalysis in Human Carbonic Anhydrase II
Supporting Information Coupling Protein Dynamics with Enzyme Catalysis in Human Carbonic Anhydrase II Srabani Taraphder,*, C. Mark Maupin, Jessica M. J. Swanson and Gregory A. Voth,* Department of Chemistry,
More informationTowards fast and accurate binding affinity. prediction with pmemdgti: an efficient. implementation of GPU-accelerated. Thermodynamic Integration
Towards fast and accurate binding affinity prediction with pmemdgti: an efficient implementation of GPU-accelerated Thermodynamic Integration Tai-Sung Lee,, Yuan Hu, Brad Sherborne, Zhuyan Guo, and Darrin
More informationEnergy Barriers and Rates - Transition State Theory for Physicists
Energy Barriers and Rates - Transition State Theory for Physicists Daniel C. Elton October 12, 2013 Useful relations 1 cal = 4.184 J 1 kcal mole 1 = 0.0434 ev per particle 1 kj mole 1 = 0.0104 ev per particle
More informationDeveloping Monovalent Ion Parameters for the Optimal Point Charge (OPC) Water Model. John Dood Hope College
Developing Monovalent Ion Parameters for the Optimal Point Charge (OPC) Water Model John Dood Hope College What are MD simulations? Model and predict the structure and dynamics of large macromolecules.
More informationSupporting Information
Supporting Information Constant ph molecular dynamics reveals ph-modulated binding of two small-molecule BACE1 inhibitors Christopher R. Ellis 1,, Cheng-Chieh Tsai 1,, Xinjun Hou 2, and Jana Shen 1, 1
More informationLarge scale free energy calculations for blind predictions of protein ligand binding: the D3R Grand Challenge 2015
J Comput Aided Mol Des (2016) 30:743 751 DOI 10.1007/s10822-016-9952-x Large scale free energy calculations for blind predictions of protein ligand binding: the D3R Grand Challenge 2015 Nanjie Deng 1 William
More information3rd Advanced in silico Drug Design KFC/ADD Molecular mechanics intro Karel Berka, Ph.D. Martin Lepšík, Ph.D. Pavel Polishchuk, Ph.D.
3rd Advanced in silico Drug Design KFC/ADD Molecular mechanics intro Karel Berka, Ph.D. Martin Lepšík, Ph.D. Pavel Polishchuk, Ph.D. Thierry Langer, Ph.D. Jana Vrbková, Ph.D. UP Olomouc, 23.1.-26.1. 2018
More informationPhotosynthetic and Protein Electron Transfer: Is Biology (Thermo) Dynamic?
Photosynthetic and Protein Electron Transfer: Is Biology (Thermo) Dynamic? University of Calgary, December 13, 2014 Dmitry Matyushov Department of Physics/Chemistry Center for Biological Physics Arizona
More informationSupporting Information. The hinge region strengthens the nonspecific interaction. between lac-repressor and DNA: a computer simulation.
Supporting Information The hinge region strengthens the nonspecific interaction between lac-repressor and DNA: a computer simulation study Lili Sun 1, Marcin Tabaka 1, Sen Hou 1, Lin Li 2, Krzysztof Burdzy
More informationComputational Chemistry - MD Simulations
Computational Chemistry - MD Simulations P. Ojeda-May pedro.ojeda-may@umu.se Department of Chemistry/HPC2N, Umeå University, 901 87, Sweden. May 2, 2017 Table of contents 1 Basics on MD simulations Accelerated
More informationEfficient Parallelization of Molecular Dynamics Simulations on Hybrid CPU/GPU Supercoputers
Efficient Parallelization of Molecular Dynamics Simulations on Hybrid CPU/GPU Supercoputers Jaewoon Jung (RIKEN, RIKEN AICS) Yuji Sugita (RIKEN, RIKEN AICS, RIKEN QBiC, RIKEN ithes) Molecular Dynamics
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 information10/26/2010. An Example of a Polar Reaction: Addition of H 2 O to Ethylene. to Ethylene
6.5 An Example of a Polar Reaction: Addition of H 2 O to Ethylene Addition of water to ethylene Typical polar process Acid catalyzed addition reaction (Electophilic addition reaction) Polar Reaction All
More informationSupplementary Information Intrinsic Localized Modes in Proteins
Supplementary Information Intrinsic Localized Modes in Proteins Adrien Nicolaï 1,, Patrice Delarue and Patrick Senet, 1 Department of Physics, Applied Physics and Astronomy, Rensselaer Polytechnic Institute,
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 informationWhen does TMAO fold a polymer chain and urea unfold it?
When does TMAO fold a polymer chain and urea unfold it? Jagannath Mondal, Guillaume Stirnemann, and B. J. Berne Department of Chemistry, Columbia University, New York, NY 7 (Dated: March 1, 1) Other studies
More informationMolecular Mechanics. Yohann Moreau. November 26, 2015
Molecular Mechanics Yohann Moreau yohann.moreau@ujf-grenoble.fr November 26, 2015 Yohann Moreau (UJF) Molecular Mechanics, Label RFCT 2015 November 26, 2015 1 / 29 Introduction A so-called Force-Field
More informationComparison of the Efficiency of the LIE and MM/GBSA Methods to Calculate LigandBinding Energies
Comparison of the Efficiency of the LIE and MM/GBSA Methods to Calculate LigandBinding Energies Genheden, Samuel; Ryde, Ulf Published in: Journal of Chemical Theory and Computation DOI: 10.1021/ct200163c
More information