Statistical Mechanics for Proteins
|
|
- Gervais Ramsey
- 5 years ago
- Views:
Transcription
1 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 ),, ( ),, ( ln ),, ( T V Q kt T V F = T V F p, = V T F S, = T V F, = μ Statistical Mechanics for Proteins
2 Exact evaluation of Q only possible for idealized systems (free particles, harmonic oscillator) From this: 2 / 3 3 free 2! 1 ),, ( h mkt T V Q = π Statistical Mechanics for Proteins Λ + = 3 ln 1 ),, ( V kt T V F V kt p = Λ + = 3 ln 2 5 V k S Λ = 3 kt ln V μ kt TS F U 2 3 = + =
3 Statistical Mechanics for Proteins Since force field is not entirely harmonic, Q for a protein can not be evaluated in closed form: Q(, V, T ) = Q = free (, V, T )* Q 1 2πmkT 3! h (, V, T ) The second term is called configurational integral and describes the contribution of the interaction part to Q. Free energies can therefore be decomposed into 3 / 2 inter * ( q) V dq exp kt ( ln Q ln Q ) F = kt ln Q = kt + free inter
4 Statistical Mechanics for Proteins Two primary strategies to evaluate the excess (nonideal, interaction-related) contribution to Q: Molecular Dynamics Simulations Sampling of the phase space (p,q) Monte Carlo Sampling of the integral
5 Free energy: classical definition + The free energy is the energy left for once you paid the tax to entropy:!g =!H"T!S Enthalpic! Hydrogen bonds! Polar interactions! Van der Waals interactions!... Entropic! Loss of degrees of freedom! Gain of vibrational modes! Loss of solvent/protein structure!... Theoretical Predictions:! Approximate: empirical formula for all contributions! Exact: using statistical physics definition of G G = -K B T ln(z)
6 Examples of factors determining the binding free energy Electrostatic interactions - Strength depends on microscopic environment (!) - Case of hydrogen bonds H H eutral : Charge assisted : H H H O O O H H H O H O H O H solvent -1.2 ± 0.6 kcal/mol -2.4 to -4.8 kcal/mol E H-bond (solv.) - E H-bond (comp.) H O H O S complex determines if H-bonds contributes to affinity or not Unpaired polar groups upon binding are detrimental Strong directional nature Specificity of molecular recognition
7 Free energy: statistical mechanics definition G =!k B Tln(Z) where Z = # e!"ei i is the partition function Free energy differences between 2 states (bound/unbound, É) are, therefore, ratios of partition functions #!G =G A "G B = "k B Tln Z A $ % Z B & ' ( Free energy simulations aim at computing this ratio using various techniques.
8 Relation with chemical equilibrium A + B " # AÕBÕ K A = K b = [ A'B' ] A [ ] [ B] AÕBÕ "# A + B [ ] [ B] [ ] K D = K i = A A'B' K b : binding constant, K d : dissociation constant, K i : inhibition constant!g binding = "RTlnK A =RTlnK D =!H"T!S "G binding (kcal/mol) Weak asso. K D (mol/l) Strong asso.
9 Connection micro/macroscopic: thermodynamics and kinetics e - RT!G = Free Energy K A Association Constant Absolute binding free energies:!g " K A Relative binding free energies:!!g " K AÕ / K A Microscopic Structure Biological function Binding free energy profiles:!g(#) " K A, K on, K off
10 The free energy is the main function behind all process A) Chemical equilibrium!g binding =RTlnK A + K A = [ AB] A [ ] [ B] A B K D =1/K A AB B) Conformational changes!g conf =k B Tln P Conf1 P Conf2 R=k B A C) Ligand binding D) É!G binding =k B Tln P Unbound P Bound
11 Free energy: computational approaches #!G =G A "G B = "k B Tln Z A $ % Z B & ' ( Free energy simulations techniques aim at computing ratios of partition functions using various techniques. Z = # i e!"ei Sampling of important microstates of the system (MD, MC, GA, É) Computation of energy of each microstate (force fields, QM, CP)
12 Connection micro/macroscopic: intuitive view E 1, P 1 ~ e -$E1 Expectation value E 2, P 2 ~ e -$E2 O = 1 Z # i O i e!"ei E 3, P 3 ~ e -$E3 # i Where Z = e!"ei is the partition function E 4, P 4 ~ e -$E4 E 5, P 5 ~ e -$E5
13 Central Role of the Partition Function Z = # i e!"ei O = 1 Z # i O i e!"ei... Expectation Value E =!!" ln(z)=u " p =k B T!ln(Z) % # $!V & ',T G = -k B T ln(z) Internal Energy Pressure Gibbs free energy
14 Molecular Modeling Principles 1) Modeling of molecular interactions Electrostatics Van der Waals Covalent bonds Solvent 2) Simulation of time evolution (ewton) 3) Computation of average values O = < O > Ensemble = < O > Temps (Ergodicity) Free energy landscape Connection microscopic/ macroscopic Macroscopic value Average simulation value
15
16
17 Ergodic Hypothesis MD Trajectory E VT simulation $ # VE simulation 3 Spatial coordinates Dialanine Protein O Ensemble = 1 % Z O(!,")e#$E(!,") d!d" = 1 % O(t)dt = O & Time? & 0
18 Free energy calculation: Main approaches Sampling, Exact Free Energy Perturbation (FEP) on Equilibrium Statistical Mechanics (Jarzynski) Thermodynamical Integration (TI) Sampling, Approx. Linear Interaction Energy (LIE) Molecular Mechanics/Poisson- Boltzmann/Surface area (MM-PBSA) CPU Time Approx. Quantitative Structure Activity Relationship (QSAR) G k 0 k i X i!g =F(X) (X is a descriptor)
19 Linear interaction energy (LIE) J. qvist, J. Phys. Chem., 1994, 98, 8253 Two MD runs : free state and bound state Free state ÒsolventÓ = water Bound state ÒsolventÓ = water + protein vdw!g bind =" E l#s ( ) + $ ( elec E ) l#s # E vdw bound l#s free # E elec bound l#s free J. qvist, J. Phys. Chem., 1994, 98, 8253 '=0.165 and (=0.5 T. Hansson et al., J. Comp.-Aided Molec. Design, 1998, 12, 27 '=0.181 and (=0.5, 0.43, 0.37, 0.33 W. Wang, Proteins, 1999, 34, 395
20 Linear interaction energy (LIE) Advantages : Drawbacks : - Faster than free energy simulation - More structurally different ligands than for free energy simulation. But generally restricted to rather similar ligands. - Slower than scores based on a single conformation (LUDI, PMF,...) - ot really universal (' and ( system dependent) - eed experimental binding affinities of known complexes Modifications : - Additional term proportional to buried surface upon complexation D.K. Jones-Hertzog and W.L. Jorgensen, J. Med. Chem., 1997, 40, Use of continuum solvent model instead of explicit solvent R. Zhou and W.L. Jorgensen et al., J. Phys. Chem., 2001, 105, 10388
21 Binding free energy decomposition: MM-PBSA, MM-GBSA Gaz lig!g solv Lig + Prot prot!g solv!e gaz comp!g solv Lig:Prot Averaged over an MD simulation trajectory of the complex (and isolated parts)!g bind =!E gaz +!G desolv "T!S Sol Lig + Prot!G bind Lig:Prot E gaz =E elec +E vdw +!E intra!g desolv =!G comp solv "!G solv lig prot ( +!G solv )!T"S =!T(S comp!(s prot +S lig )) S =S trans +S rot +S vib B. Tidor and M. Karplus, J. Mol. Biol., 1994, 238, 405!G solv =!G solv,elec +!G solv,np comp!g desolv =!G solv,elec Depending on the way!g solv,elec is calculated: prot ( +!G solv,elec ) +# SASA comp " SASA lig +SASA prot lig "!G solv,elec Molecular mechanics Ð Poisson-Boltzmann Surface Area (MM- PBSA) J. Srinivasan, P.A. Kollmann et al., J. Am. Chem. Soc., 1998, 120, 9401 Molecular mechanics Ð Generalized Born Surface Area (MM- GBSA) H. Gohlke, C. Kiel and D.A. Case, J. Mol. Biol., 2003, 330, 891 ( ) ( )
22 Binding free energy decomposition MM- PBSA, MM-GBSA Advantages : Drawbacks : - Used for ligand:protein and protein:protein complexes - Could be applied to structurally different ligands (but in fact applied to similar ones) - ÒUniversalÓ (no parameter to be fitted) - MM-GBSA allows a per-atom decomposition of "G bind (e.g. contribution of side chains) - Rather time consuming - In some cases, found unable to rank ligands -T"S is necessary to find the order of magnitude of the absolute binding free energies but, in some cases, it is not necessary to estimate relative binding free energies W. Wang and P.A. Kollman, J. Mol. Biol., 2000, 303, 567 H. Gohlke, C. Kiel and D.A. Case, J. Mol. Biol., 2003, 330, 891
23
24 Use of thermodynamical cycles L1 + Prot "G 1 L1:Prot Thermodynamic cycle perturbation approach: "G 3 "G 4 ""G bind ="G 2 -"G 1 ="G 4 -"G 3 L2 + Prot "G 2 L2:Prot "G 4 -"G 3 is computationally accessible ÒAlchemicalÓ reaction. MD or MC at different *. Coupling parameter * H! = H 0 +! H L1 + ( 1"! ) H L2 O O H H O CH 2 O OH CH 2 OH Cl O H O CH 2 Cl *=0 * *=1 Free energy perturbation (FEP) Thermodynamic integration (TI) $$ G bind =# RT n " # 1 i= 0 ln exp (# ( H # H ) RT)! i! i+ i! i ## G bind = $! = 1 " H $ = 0 " $ $ d$ $
25 Alchemical free energy formalism
26 ÒAlchemicalÓ Free Energy Calculations
27 Hybrid Side Chain for P "A Mutation
28 Results of the Free Energy Simulations Total (Path Independent) Experimental: Theoretical: 2.9 (0.2) kcal/mol 2.9 (1.1) kcal/mol Components (Path Dependent) Experimental: - Theoretical: TCR 25% Solvent 20% HLA A2 40% Peptide 15% } } 45% 55%
29 Free energy formalism From the statistical definition of the free energy,
30 Concluding Remarks 1) Free energy simulations can reproduce accurately experimental changes in association constant between to closely related protein systems if detailed structural knowledge is available (X-ray, MR or model) 2) The formalism is exact from a statistical physics stand point and accurate treatment of entropic terms, solvent effect or conformational changes can be obtained 3) Convergence of the free energy derivative is still problematic. The situation should improve with new methodological enhancements as well as longer simulation time 4) Absolute free energies can also be computed but the convergence is even more difficult 5) Much details about the specificity of the association can be gained using component analysis, opening the door to rational peptide or protein design
Molecular 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 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 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 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 informationMMGBSA: Thermodynamics of Biomolecular Systems
MMGBSA: Thermodynamics of Biomolecular Systems The MMGBSA approach employs molecular mechanics, the generalized Born model and solvent accessibility method to elicit free energies from structural information
More informationAlchemical 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 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 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 informationInsights into Protein Protein Binding by Binding Free Energy Calculation and Free Energy Decomposition for the Ras Raf and Ras RalGDS Complexes
doi:10.1016/s0022-2836(03)00610-7 J. Mol. Biol. (2003) 330, 891 913 Insights into Protein Protein Binding by Binding Free Energy Calculation and Free Energy Decomposition for the Ras Raf and Ras RalGDS
More informationGherman Group Meeting. Thermodynamics and Kinetics and Applications. June 25, 2009
Gherman Group Meeting Thermodynamics and Kinetics and Applications June 25, 2009 Outline Calculating H f, S, G f Components which contribute to H f, S, G f Calculating ΔH, ΔS, ΔG Calculating rate constants
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 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 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 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 informationk θ (θ θ 0 ) 2 angles r i j r i j
1 Force fields 1.1 Introduction The term force field is slightly misleading, since it refers to the parameters of the potential used to calculate the forces (via gradient) in molecular dynamics simulations.
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 informationMM-GBSA for Calculating Binding Affinity A rank-ordering study for the lead optimization of Fxa and COX-2 inhibitors
MM-GBSA for Calculating Binding Affinity A rank-ordering study for the lead optimization of Fxa and COX-2 inhibitors Thomas Steinbrecher Senior Application Scientist Typical Docking Workflow Databases
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 informationMonte Carlo (MC) Simulation Methods. Elisa Fadda
Monte Carlo (MC) Simulation Methods Elisa Fadda 1011-CH328, Molecular Modelling & Drug Design 2011 Experimental Observables A system observable is a property of the system state. The system state i is
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 informationPhase Equilibria and Molecular Solutions Jan G. Korvink and Evgenii Rudnyi IMTEK Albert Ludwig University Freiburg, Germany
Phase Equilibria and Molecular Solutions Jan G. Korvink and Evgenii Rudnyi IMTEK Albert Ludwig University Freiburg, Germany Preliminaries Learning Goals Phase Equilibria Phase diagrams and classical thermodynamics
More informationOther Cells. Hormones. Viruses. Toxins. Cell. Bacteria
Other Cells Hormones Viruses Toxins Cell Bacteria ΔH < 0 reaction is exothermic, tells us nothing about the spontaneity of the reaction Δ H > 0 reaction is endothermic, tells us nothing about the spontaneity
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: Potential Energy Functions
Lecture 11: Potential Energy Functions Dr. Ronald M. Levy ronlevy@temple.edu Originally contributed by Lauren Wickstrom (2011) Microscopic/Macroscopic Connection The connection between microscopic interactions
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 informationProtein-Ligand Docking Methods
Review Goal: Given a protein structure, predict its ligand bindings Protein-Ligand Docking Methods Applications: Function prediction Drug discovery etc. Thomas Funkhouser Princeton University S597A, Fall
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 information7/19/2011. Models of Solution. State of Equilibrium. State of Equilibrium Chemical Reaction
Models of Solution Chemistry- I State of Equilibrium A covered cup of coffee will not be colder than or warmer than the room temperature Heat is defined as a form of energy that flows from a high temperature
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 informationThermodynamics and Kinetics
Thermodynamics and Kinetics C. Paolucci University of Notre Dame Department of Chemical & Biomolecular Engineering What is the energy we calculated? You used GAMESS to calculate the internal (ground state)
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 informationSolvent & geometric effects on non-covalent interactions
Solvent & geometric effects on non-covalent interactions Scott L. Cockroft PhysChem Forum 10, Syngenta, Jealott s Hill, 23 rd March 11 QSAR & Physical Organic Chemistry Quantifiable Physicochemical Properties
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 informationProtein-Ligand Interactions* and Energy Evaluation Methods
Protein-Ligand Interactions* and Energy Evaluation Methods *with a revealing look at roles of water Glen E. Kellogg Department of Medicinal Chemistry Institute for Structural Biology, Drug Discovery &
More informationDocking. GBCB 5874: Problem Solving in GBCB
Docking Benzamidine Docking to Trypsin Relationship to Drug Design Ligand-based design QSAR Pharmacophore modeling Can be done without 3-D structure of protein Receptor/Structure-based design Molecular
More informationKd = koff/kon = [R][L]/[RL]
Taller de docking y cribado virtual: Uso de herramientas computacionales en el diseño de fármacos Docking program GLIDE El programa de docking GLIDE Sonsoles Martín-Santamaría Shrödinger is a scientific
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 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 informationStatistical thermodynamics for MD and MC simulations
Statistical thermodynamics for MD and MC simulations knowing 2 atoms and wishing to know 10 23 of them Marcus Elstner and Tomáš Kubař 22 June 2016 Introduction Thermodynamic properties of molecular systems
More informationStatistical Thermodynamics and Monte-Carlo Evgenii B. Rudnyi and Jan G. Korvink IMTEK Albert Ludwig University Freiburg, Germany
Statistical Thermodynamics and Monte-Carlo Evgenii B. Rudnyi and Jan G. Korvink IMTEK Albert Ludwig University Freiburg, Germany Preliminaries Learning Goals From Micro to Macro Statistical Mechanics (Statistical
More informationProtein-Ligand Docking Methods
Review Goal: Given a protein structure, predict its ligand bindings Protein-Ligand Docking Methods Applications: Function prediction Drug discovery etc. Thomas Funkhouser Princeton University S597A, Fall
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 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 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 informationCOSMO-RS Theory. The Basics
Theory The Basics From µ to properties Property µ 1 µ 2 activity coefficient vapor pressure Infinite dilution Gas phase Pure compound Pure bulk compound Partition coefficient Phase 1 Phase 2 Liquid-liquid
More information3.320: Lecture 19 (4/14/05) Free Energies and physical Coarse-graining. ,T) + < σ > dµ
3.320: Lecture 19 (4/14/05) F(µ,T) = F(µ ref,t) + < σ > dµ µ µ ref Free Energies and physical Coarse-graining T S(T) = S(T ref ) + T T ref C V T dt Non-Boltzmann sampling and Umbrella sampling Simple
More informationSupplementary Methods
Supplementary Methods MMPBSA Free energy calculation Molecular Mechanics/Poisson Boltzmann Surface Area (MM/PBSA) has been widely used to calculate binding free energy for protein-ligand systems (1-7).
More informationLecture 5: Electrostatic Interactions & Screening
Lecture 5: Electrostatic Interactions & Screening Lecturer: Prof. Brigita Urbanc (brigita@drexel.edu) PHYS 461 & 561, Fall 2009-2010 1 A charged particle (q=+1) in water, at the interface between water
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 informationTransition Theory Abbreviated Derivation [ A - B - C] # E o. Reaction Coordinate. [ ] # æ Æ
Transition Theory Abbreviated Derivation A + BC æ Æ AB + C [ A - B - C] # E A BC D E o AB, C Reaction Coordinate A + BC æ æ Æ æ A - B - C [ ] # æ Æ æ A - B + C The rate of reaction is the frequency of
More informationCE 530 Molecular Simulation
1 CE 530 Molecular Simulation Lecture 14 Molecular Models David A. Kofke Department of Chemical Engineering SUNY Buffalo kofke@eng.buffalo.edu 2 Review Monte Carlo ensemble averaging, no dynamics easy
More informationStatistical mechanics of biological processes
Statistical mechanics of biological processes 1 Modeling biological processes Describing biological processes requires models. If reaction occurs on timescales much faster than that of connected processes
More informationIntro/Review of Quantum
Intro/Review of Quantum QM-1 So you might be thinking I thought I could avoid Quantum Mechanics?!? Well we will focus on thermodynamics and kinetics, but we will consider this topic with reference to the
More information5th CCPN Matt Crump. Thermodynamic quantities derived from protein dynamics
5th CCPN 2005 -Matt Crump Thermodynamic quantities derived from protein dynamics Relaxation in Liquids (briefly!) The fluctuations of each bond vector can be described in terms of an angular correlation
More informationIntro/Review of Quantum
Intro/Review of Quantum QM-1 So you might be thinking I thought I could avoid Quantum Mechanics?!? Well we will focus on thermodynamics and kinetics, but we will consider this topic with reference to the
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 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 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 informationMonte Carlo. Lecture 15 4/9/18. Harvard SEAS AP 275 Atomistic Modeling of Materials Boris Kozinsky
Monte Carlo Lecture 15 4/9/18 1 Sampling with dynamics In Molecular Dynamics we simulate evolution of a system over time according to Newton s equations, conserving energy Averages (thermodynamic properties)
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 informationMD Simulation in Pose Refinement and Scoring Using AMBER Workflows
MD Simulation in Pose Refinement and Scoring Using AMBER Workflows Yuan Hu (On behalf of Merck D3R Team) D3R Grand Challenge 2 Webinar Department of Chemistry, Modeling & Informatics Merck Research Laboratories,
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 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 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 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 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 informationBasics of Statistical Mechanics
Basics of Statistical Mechanics Review of ensembles Microcanonical, canonical, Maxwell-Boltzmann Constant pressure, temperature, volume, Thermodynamic limit Ergodicity (see online notes also) Reading assignment:
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 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 information510 Subject Index. Hamiltonian 33, 86, 88, 89 Hamilton operator 34, 164, 166
Subject Index Ab-initio calculation 24, 122, 161. 165 Acentric factor 279, 338 Activity absolute 258, 295 coefficient 7 definition 7 Atom 23 Atomic units 93 Avogadro number 5, 92 Axilrod-Teller-forces
More informationProtein-Ligand Docking
Protein-Ligand Docking Matthias Rarey GMD - German National Research Center for Information Technology Institute for Algorithms and Scientific Computing (SCAI) 53754Sankt Augustin, Germany rarey@gmd.de
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 If you throw up a stone, it is Physics. If you throw up a stone, it is Physics. If it lands on your head, it is
More informationChemical properties that affect binding of enzyme-inhibiting drugs to enzymes
Chemical properties that affect binding of enzyme-inhibiting drugs to enzymes Introduction The production of new drugs requires time for development and testing, and can result in large prohibitive costs
More informationClassical Monte-Carlo simulations
Classical Monte-Carlo simulations Graduate Summer Institute Complex Plasmas at the Stevens Insitute of Technology Henning Baumgartner, A. Filinov, H. Kählert, P. Ludwig and M. Bonitz Christian-Albrechts-University
More informationExample questions for Molecular modelling (Level 4) Dr. Adrian Mulholland
Example questions for Molecular modelling (Level 4) Dr. Adrian Mulholland 1) Question. Two methods which are widely used for the optimization of molecular geometies are the Steepest descents and Newton-Raphson
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 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 informationMM-PBSA Validation Study. Trent E. Balius Department of Applied Mathematics and Statistics AMS
MM-PBSA Validation Study Trent. Balius Department of Applied Mathematics and Statistics AMS 535 11-26-2008 Overview MM-PBSA Introduction MD ensembles one snap-shots relaxed structures nrichment Computational
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 informationUnderstanding Chemical Reactions through Computer Modeling. Tyler R. Josephson University of Delaware 4/21/16
Understanding Chemical Reactions through Computer Modeling Tyler R. Josephson University of Delaware 4/21/16 A little about me B.S. in Chem E from U of M, 2011 Currently, Ph.D. student at University of
More informationThermodynamics of Three-phase Equilibrium in Lennard Jones System with a Simplified Equation of State
23 Bulletin of Research Center for Computing and Multimedia Studies, Hosei University, 28 (2014) Thermodynamics of Three-phase Equilibrium in Lennard Jones System with a Simplified Equation of State Yosuke
More informationEntropy and Free Energy in Biology
Entropy and Free Energy in Biology Energy vs. length from Phillips, Quake. Physics Today. 59:38-43, 2006. kt = 0.6 kcal/mol = 2.5 kj/mol = 25 mev typical protein typical cell Thermal effects = deterministic
More informationSolubility Properties
Solubility Properties X-ray crystal structure of : dicyclohexyl[18]crown-6 and potassium complex of [18]crown-6. Structures are different: In the solid state In polar and apolar solvent As a complex 1
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 informationResearch Statement. Shenggao Zhou. November 3, 2014
Shenggao Zhou November 3, My research focuses on: () Scientific computing and numerical analysis (numerical PDEs, numerical optimization, computational fluid dynamics, and level-set method for interface
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 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 informationMOLECULAR RECOGNITION DOCKING ALGORITHMS. Natasja Brooijmans 1 and Irwin D. Kuntz 2
Annu. Rev. Biophys. Biomol. Struct. 2003. 32:335 73 doi: 10.1146/annurev.biophys.32.110601.142532 Copyright c 2003 by Annual Reviews. All rights reserved First published online as a Review in Advance on
More informationENZYME SCIENCE AND ENGINEERING PROF. SUBHASH CHAND DEPARTMENT OF BIOCHEMICAL ENGINEERING AND BIOTECHNOLOGY IIT DELHI LECTURE 3
ENZYME SCIENCE AND ENGINEERING PROF. SUBHASH CHAND DEPARTMENT OF BIOCHEMICAL ENGINEERING AND BIOTECHNOLOGY IIT DELHI LECTURE 3 ENZYMES AS BIOCATALYSTS * CATALYTIC EFFICIENCY *SPECIFICITY Having discussed
More informationBioengineering 215. An Introduction to Molecular Dynamics for Biomolecules
Bioengineering 215 An Introduction to Molecular Dynamics for Biomolecules David Parker May 18, 2007 ntroduction A principal tool to study biological molecules is molecular dynamics simulations (MD). MD
More informationWater models in classical simulations
Water models in classical simulations Maria Fyta Institut für Computerphysik, Universität Stuttgart Stuttgart, Germany Water transparent, odorless, tasteless and ubiquitous really simple: two H atoms attached
More informationDFT modeling of novel materials for hydrogen storage
DFT modeling of novel materials for hydrogen storage Tejs Vegge 1, J Voss 1,2, Q Shi 1, HS Jacobsen 1, JS Hummelshøj 1,2, AS Pedersen 1, JK Nørskov 2 1 Materials Research Department, Risø National Laboratory,
More informationFree energy calculations with alchemlyb
Free energy calculations with alchemlyb Oliver Beckstein Arizona State University SPIDAL Teleconference 2019-02-01 Binding free energy measure of how strong a protein P and a ligand X stick together key
More information= (-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 informationBIOC : Homework 1 Due 10/10
Contact information: Name: Student # BIOC530 2012: Homework 1 Due 10/10 Department Email address The following problems are based on David Baker s lectures of forces and protein folding. When numerical
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 information1 What is energy?
http://www.intothecool.com/ 1 What is energy? the capacity to do work? (Greek: en-, in; + ergon, work) the capacity to cause change to produce an effect? a certain quantity that does not change in the
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 informationIntroduction Statistical Thermodynamics. Monday, January 6, 14
Introduction Statistical Thermodynamics 1 Molecular Simulations Molecular dynamics: solve equations of motion Monte Carlo: importance sampling r 1 r 2 r n MD MC r 1 r 2 2 r n 2 3 3 4 4 Questions How can
More information