Lecture 12: Solvation Models: Molecular Mechanics Modeling of Hydration Effects
|
|
- Leslie Bradford
- 5 years ago
- Views:
Transcription
1 Statistical Thermodynamics Lecture 12: Solvation Models: Molecular Mechanics Modeling of Hydration Effects Dr. Ronald M. Levy
2 Bare Molecular Mechanics Atomistic Force Fields: torsion stretching E stretch = bonds K r r r 0 2 bending E bend = angles K q q q 0 2 non-bonded E torsion = torsions V n 2 [ 1±cos n n ] E non bonded = ij [ q i q j r ij 4 ij 12 ij r ij ] 6 ij 12 6 r ij
3 Hydration has a large effect on the conformations of macromolecules MD simulation, RMSD from native native:
4 and on ligand binding Distribution of complex decoy binding energies:
5 Challenge Model hydration conveniently, as accurately as needed by the application and with the least computational cost. Explicit solvation models. Implicit solvation models. Two main approaches
6 Explicit Solvation Each solvent molecule is represented with a set of atomic interaction centers (just as for the solute). Most accurate/detailed. Computationally expensive. Requires averaging over solvent coordinates. Difficult to obtain relative free energies of solute conformations.
7 Implicit Solvation The solvent is represented by a continuum described by macroscopic parameters such as the dielectric constant, density, surface tension, etc. Theoretical framework based on solvent PMF. Not as accurate, especially for short-range solute-solvent interactions. Reduced dimensionality. Relative solvation free energies from single point effective potential energy calculations.
8 Many recent developments in molecular modeling have focused on solvation models. Explicit: long-range electrostatic models Implicit: improve coverage, accuracy and efficiency. Sometimes it is hard to keep track of all the choices. From a 2007 publication: The peptide is simulated in TIP3P and several variations of the GB implicit solvent model: GBHCT, GBOBC, and GBNeck (igb = 1, 5, and 7, respectively, in Amber 9). [ ] For consistency, MBondi radii were used in both the GB REMD simulations and subsequent GB and PE energy calculations described below. For TIP3P simulations, Ala10 was solvated in a truncated octahedral box with 983 solvent molecules [ ] long range electrostatic interactions were calculated using periodic boundary conditions via the particle mesh Ewald (PME). Roe, Wickstrom, et al. JPC B 111:1846 (2007)
9 Explicit Hydration: Water Models A water model specifies the interaction sites of a water molecule and their partial charges and LJ parameters. There is a large variety of models: rigid/flexible, 3-points/4-points or more, polarizable, dissociable, The 3 rigid models most frequently used are SPC, TIP3P and TIP4P: SPC TIP3P TIP4P r(oh),å HOH, deg σ(o), Å ε(o),kcal/mol q(o) q(h) q(m) -1.04
10 Periodic Boundary Conditions (PBC s) PBC s are often used to simulate a bulk solution with a relatively small number of molecules. The central simulation box is replicated in 3D space so that the system has no edges. The minimum image convention says that the actual distance between two particles is the smallest distance between all of the distances between their images: Δx min =Δx L nint Δx L where L is the box size and x is the distance (along x) between any two images of a pair of particles. nint()=nearest integer. A variety of space filling box shapes can be used. The truncated octahedron has particularly nice properties.
11 Treatment of long-range electrostatic interactions Electrostatic energy (per cell) of a periodic system: U c = 1 2 i q i φ r i φ r i = j,n q j r ij nl The sum for the electrostatic potential φ includes all of the particles in the unit cell plus all of their images. Straight sum is conditionally convergent ; for example it would diverge if the positive and negative charges are summed separately. Unlike LJ interactions, in this case truncation of the sum (with a distance cutoff) can be inaccurate. In explicit solvent molecular simulation programs the infinite sum problem is circumvented using the Ewald summation class of models.
12 Ewald summation methods The Ewald sum is based on a decomposition of the 1/r interaction into a fast decaying component (handled with distance cutoffs) and a smooth component (summed using a Fourier transformation in reciprocal space). Think of surrounding each charge by a canceling Gaussian charge distribution: c, i = q i / 3 /2 e r 2 ij 1 r ij erfc r r ρ k is the Fourier transform of the charge distribution: k ρ k 2 exp k 2 /4α k 2 k = i q i e i k r i measures the softness of the canceling distribution - the larger the more convergent is the direct sum (but less convergent is the reciprocal sum).
13 Particle mesh Ewald methods Straightforward implementations of the Ewald sum scale at best as N 3/2 More efficient implementations of complexity NlogN exploit FFT techniques by expressing the charge density on a regular grid. The original idea of Hockeny and Eastwood (1981) consisted of assigning fractional charges to grid points: PPPM or P 3 M Darden, York and Pedersen (1993) presented an (equivalent) approach based on Lagrange interpolation of the structure factors exp(ikr) on a grid: PME PME was later refined using B-spline interpolation by Essmann et al. (1995) to obtain gradients for MD: SPME In most recent papers, PME is SPME
14 Fast Multipole Methods Based on representing distant charge distributions by multipole expansions. Very favorable asymptotic scaling of order N but with high initial overhead. For values of N typical for biomolecular simulations FMM tends not to be as efficient as PME-based methods. Still the only efficient method to treat longrange interactions without periodic boundary conditions. Used routinely in simulations of million-body simulations of galaxies, etc. A future for FMM in the context of implicit solvent modeling? Figuerido, Levy, Zhou, Berne. JCP 106:9835 (1997)
15 Implicit Solvation: the solvent potential of mean force (SPMF) Write the canonical partition function Q xy of the solution in terms of solute coordinates x, and solvent coordinates y. The average of any observable that depends only on x can be written in terms of an effective solvation potential W(x). where 1 O x =Q xy dxdyo x exp { β [U x U x, y U y ] } O( x) = 1 Q x dxo( x)exp { β [U ( x)+w ( x )] } exp [ βw ( x)]= dy exp { β [U ( x, y)+u ( y) ] } dy exp { βu ( y) } defines the solvent potential of mean force So called because the gradients of W(x) with respect to the solute coordinates are equal to the average forces of the solvent on the solute atoms at fixed solute conformation. W(x) implicitly contains all of the effects of solvation.
16 The solvent potential of mean force is a free energy Conformational equilibria: x 2 x 1 ΔG 21 = kt ln P( x 2 ) P( x 1 ) = ΔU 21 + ΔW 21 Solvation (fixed solute conformation): y x Q sep =exp [ βu x ]Q y y x Q sol =exp [ βu ( x)] dy exp { β [U ( x, y)+u ( y )] } =exp { β [U ( x)+w ( x)] }Q y ΔG solv ( x )= kt ln Q sol Q sep = kt ln exp [ βw ( x)]=w ( x)
17 Polar/Non-polar decomposition of the solvent PMF ( v) U ch Δ G solv ( x)=w ( x) G np G cav ( w ) +U ch Δ G vdw G solv = G elec G np G np = G cav G vdw ( Δ G elec =U w) ( v ) ch U ch
18 Typical Modern Implicit Solvent Model Electrostatic Component: Continuum Dielectric Poisson-Boltzmann solvers (accurate but numerical and slow). Generalized Born models (faster, can be expressed as analytic function). Non-Polar Component: Solute surface area models Cavity + van der Waals NP models.
19 The Poisson Equation (PE) and the Poisson-Boltzmann (PB) equations PE for a non-homogeneous dielectric: [ ε x ϕ x ]= ρ x Electrostatic free energy of solute: G elec = 1 2 q i i ϕ( x i ) Electrostatic solvation free energy: ε=1 ε= ε slv + - W elec =G elec ( ε in =1, ε out =ε slv ) G elec (ε in =1, ε out =1) W elec = 1 2 i q i ϕ ( x i ) 1 2 i q i ϕ q ( x i )= 1 2 i q i ϕ rf ( x i ) direct Coulomb The reaction field potential φ rf is due to the polarization of the medium.
20 The PB equation takes into account the effect of free ions in solution [ x x ]= x = i q i x x i ions c j z j e z j x Debye-Hückel ion distribution - ε=1 + - When βz j ϕ x << 1 (otherwise ion adsorption): ε= ε slv + linearized form of the PB equation. [ ε x ϕ x ] k 2 x ε x ϕ x = q i δ x x i atoms where k 2 x = 2I kt ε x with I = 1 2 j c j z j 2 (ionic strength).
21 Numerical solutions of the PB equation The PB equation is solved on a grid in both surface and volume formulations. Finite difference: solves the PB equation on a volume grid (APBS, Delphi, UHBD) Finite element: solves integral form of the equation on a volume grid (PBF) Boundary element: surface grid. PB solvers often available in molecular simulation packages: Amber, CHARMM, IMPACT, etc. Main drawback: continuum dielectric models are not suitable for specific short-range solute-solvent interactions, finite size effects, non-linear effects, high ionization states. Other limitations are dependence on atomic radii parameters, speed, lack of analytical derivatives, dependence on frame of reference.
22 Approximate continuum dielectric models The basic idea is that a dielectric model of hydration should describe these two basic effects: 1.Dielectric polarization around polar groups Favorable interaction between exposed charged atoms and the polarized dielectric. Born model of ion hydration: ΔG elec = ε w z2 R Dielectric screening of electrostatic interactions The dielectric weakens the interactions between charges Distance-dependent dielectric models u ij = q i q j ε r ij r ij ε r r
23 Generalized Born Model W elec = 1 2 ( 1 1 ε w ) ij q i q j f ij (r ij ) f ij =[r 2 ij B i B j exp r 2 ij /4 B i B j ] 1/2 B i is the Born radius of atom i defined by: i W single = 1 ( 1 1 ) q 2 i 1 2 ε w B i 8π ( 1 1 ε w ) V q i 2 r r i 4 d3 r Satisfies the Born model in the two limits of infinite separation and complete overlap of solute atoms. Basically it s an interpolation formula.
24 Overall Features of Generalized Born Models W elec = 1 2 ( 1 1 ε w ) ij q i q j f ij (r ij ) The GB model works because it describes both dielectric polarization and dielectric screening effects. Polarization i=j ( self energy): i W single = 1 ( 1 1 ) q 2 i 2 ε w B i Dielectric screening i j (pair energies): u ij = q i q j ε r ij r ij ε ij r =S S x = [ ε 1 Favors the solvent exposure (small B i ) of polar groups (large q). r B i B j 1 x 2 exp x 2 /4 ] 1 S x
25 Born Radii: Pairwise Descreening Scheme Born radius calculated as a sum of pairwise atomic contributions: 1 = 1 B i 4 dielectric 1 r i r d 3 r= R i 4 solute-i 1 r i r 4 d 3 r Solute i i i Solvent Dielectric Summing over all j s: 1 B i 1 R i 1 4 j 1 = 1 1 = 1 1 B i R i B i R i 4 atom j Q ij j 1 r i r 4 d 3 r Q ij : pairwise descreening function. Forms the basis for most analytical GB models used in molecular simulations.
26 Pairwise Descreening Gotchas Q ij is not simply the integral of 1/ r i - r 4 over atom j because: 1. Overlaps with atom j and other atoms would be over-counted j i Reduce descreening contribution from atom j 2. Atoms i and j may overlap j i Integrate only over portion of j not overlapping with i 3. Both situations can and do occur simultaneously.
27 GB implementations Most major biomolecular simulation packages (CHARMM, Amber, IMPACT, Gromacs, etc.) include pairwise descreening GB implementations suitable for MD calculations. Key ingredients are the atomic radii and the description of the solute volume. The atoms overlap problem is generally addressed by empirical scaling coefficients parameterized with respect to higher level calculations that is the geometric model is parameterized in addition to the energetic model (ACE, GB/SA, GBHCT, GBSW) Work on the AGBNP series of models shows that geometric parameterization is unnecessary. Some implementations (GBMV, SGB) perform numerical integration on a grid (volume or surface) non-analytic, higher computational cost, difficulties with derivatives, dependence on coordinate frame. Some implementations differ in the choice of the GB distance function f(r) Many of the models include continuum dielectric correction terms. Recent developments have focused on the interstitial volumes problem (GBneck, GBMV, AGBNP2).
28 Non-Polar Hydration Free Energy Defined as the work of introducing the uncharged solute into solution. In some ways a harder problem for modeling than electrostatics. Often (and, nearly as often, incorrectly) ignored because appears small in magnitude. Traditionally modeled in terms of the solvent-exposed surface of the solute by means of a surface tension parameter (SA models) ΔG np =γa Motivated by macroscopic interfacial models and theories of hydration of cavities (probability of spontaneous occurrence of voids) in water. Recent developments have moved beyond surface-only models adopting distinct geometric models for the cavity and van der Waals components of non-polar hydration (NP models): G np = G cav G vdw
29 Example of an analytical NP model (the NP in AGBNP) DG np = G cav G vdw Δ G np = i [ γ i A i +α i ω( B i )] A i ω( B i ) i,α i : Surface area of atom i : Geometrical predictor based on Born radius : Surface tension and van der Waals adjustable parameters 6-4ϵ ω i ρ w i σ i slv. r r i =-16 π ρ 6 wϵ i σ i 6 3 3C i C i= 3 4 slv. 1 r r i 6 1 /3 B i
30 Advantages of NP formulation The free energy of cavity formation is defined by the solute geometry (by the surface area, say). Van der Walls interactions are longer-ranged and also depend on the properties of individual atoms (C vs. O, say). It makes sense to model these two components independently. PMF of dimerization of uncharged alanine dipeptide. Solute-solvent van der Waals energy changes upon binding. NP explicit solvent Non-monotonic behavior can not be reproduced with a surface area model. Large scatter relative to surface area model due to residual interactions of buried atoms.
Lecture 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 informationNon-bonded interactions
speeding up the number-crunching continued Marcus Elstner and Tomáš Kubař December 3, 2013 why care? number of individual pair-wise interactions bonded interactions proportional to N: O(N) non-bonded interactions
More informationNon-bonded interactions
speeding up the number-crunching Marcus Elstner and Tomáš Kubař May 8, 2015 why care? key to understand biomolecular structure and function binding of a ligand efficiency of a reaction color of a chromophore
More informationComputer simulation methods (2) Dr. Vania Calandrini
Computer simulation methods (2) Dr. Vania Calandrini in the previous lecture: time average versus ensemble average MC versus MD simulations equipartition theorem (=> computing T) virial theorem (=> computing
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 informationAll-atom Molecular Mechanics. Trent E. Balius AMS 535 / CHE /27/2010
All-atom Molecular Mechanics Trent E. Balius AMS 535 / CHE 535 09/27/2010 Outline Molecular models Molecular mechanics Force Fields Potential energy function functional form parameters and parameterization
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 informationComputational Molecular Biophysics. Computational Biophysics, GRS Jülich SS 2013
Computational Molecular Biophysics Computational Biophysics, GRS Jülich SS 2013 Computational Considerations Number of terms for different energy contributions / E el ~N ~N ~N ~N 2 Computational cost for
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 informationThe generalized Born model: its foundation, applications, and limitations
The generalized Born model: its foundation, applications, and limitations Alexey Onufriev September 8, 2010 Abstract The generalized Born (GB) approximation is introduced in the context of the implicit
More informationAdvanced Molecular Molecular Dynamics
Advanced Molecular Molecular Dynamics Technical details May 12, 2014 Integration of harmonic oscillator r m period = 2 k k and the temperature T determine the sampling of x (here T is related with v 0
More informationIntroduction to molecular dynamics
1 Introduction to molecular dynamics Yves Lansac Université François Rabelais, Tours, France Visiting MSE, GIST for the summer Molecular Simulation 2 Molecular simulation is a computational experiment.
More informationSolvation and Macromolecular Structure. The structure and dynamics of biological macromolecules are strongly influenced by water:
Overview Solvation and Macromolecular Structure The structure and dynamics of biological macromolecules are strongly influenced by water: Electrostatic effects: charges are screened by water molecules
More informationPerformance Comparison of Generalized Born and Poisson Methods in the Calculation of Electrostatic Solvation Energies for Protein Structures
Performance Comparison of Generalized Born and Poisson Methods in the Calculation of Electrostatic Solvation Energies for Protein Structures MICHAEL FEIG, 1 ALEXEY ONUFRIEV, 1 MICHAEL S. LEE, 2 WONPIL
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 informationThe AGBNP2 Implicit Solvation Model
2544 J. Chem. Theory Comput. 2009, 5, 2544 2564 The AGBNP2 Implicit Solvation Model Emilio Gallicchio,* Kristina Paris, and Ronald M. Levy Department of Chemistry and Chemical Biology and BioMaPS Institute
More informationMolecular Simulation III
Molecular Simulation III Quantum Chemistry Classical Mechanics E = Ψ H Ψ ΨΨ U = E bond +E angle +E torsion +E non-bond Molecular Dynamics Jeffry D. Madura Department of Chemistry & Biochemistry Center
More informationPeptide folding in non-aqueous environments investigated with molecular dynamics simulations Soto Becerra, Patricia
University of Groningen Peptide folding in non-aqueous environments investigated with molecular dynamics simulations Soto Becerra, Patricia IMPORTANT NOTE: You are advised to consult the publisher's version
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 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 informationImplicit modeling of nonpolar solvation for simulating protein folding and conformational transitions
PERSPECTIVE www.rsc.org/pccp Physical Chemistry Chemical Physics Implicit modeling of nonpolar solvation for simulating protein folding and conformational transitions Jianhan Chenw and Charles L. Brooks
More 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 informationComputer simulation of Biological Systems
UNIVERSITY OF TRENTO Faculty of Mathematical, Physical and Natural Sciences Physics Department Ph.D. Thesis in Physics Computer simulation of Biological Systems Supervisors: Dr. Maurizio Dapor Dr. Giovanni
More informationINTRODUCTION TO FINITE ELEMENT METHODS ON ELLIPTIC EQUATIONS LONG CHEN
INTROUCTION TO FINITE ELEMENT METHOS ON ELLIPTIC EQUATIONS LONG CHEN CONTENTS 1. Poisson Equation 1 2. Outline of Topics 3 2.1. Finite ifference Method 3 2.2. Finite Element Method 3 2.3. Finite Volume
More informationMolecular Modeling -- Lecture 15 Surfaces and electrostatics
Molecular Modeling -- Lecture 15 Surfaces and electrostatics Molecular surfaces The Hydrophobic Effect Electrostatics Poisson-Boltzmann Equation Electrostatic maps Electrostatic surfaces in MOE 15.1 The
More informationAn FPGA Implementation of Reciprocal Sums for SPME
An FPGA Implementation of Reciprocal Sums for SPME Sam Lee and Paul Chow Edward S. Rogers Sr. Department of Electrical and Computer Engineering University of Toronto Objectives Accelerate part of Molecular
More informationScuola di Chimica Computazionale
Societa Chimica Italiana Gruppo Interdivisionale di Chimica Computazionale Scuola di Chimica Computazionale Introduzione, per Esercizi, all Uso del Calcolatore in Chimica Organica e Biologica Modellistica
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 informationGENERALIZED BORN MODELS OF MACROMOLECULARSOLVATION EFFECTS
Annu. Rev. Phys. Chem. 2000. 51:129 52 Copyright c 2000 by Annual Reviews. All rights reserved GENERALIZED BORN MODELS OF MACROMOLECULARSOLVATION EFFECTS Donald Bashford and David A. Case Department of
More informationG pol = 1 2 where ρ( x) is the charge density at position x and the reaction electrostatic potential
EFFICIENT AND ACCURATE HIGHER-ORDER FAST MULTIPOLE BOUNDARY ELEMENT METHOD FOR POISSON BOLTZMANN ELECTROSTATICS CHANDRAJIT BAJAJ SHUN-CHUAN CHEN Abstract. The Poisson-Boltzmann equation is a partial differential
More informationOptimizing GROMACS for parallel performance
Optimizing GROMACS for parallel performance Outline 1. Why optimize? Performance status quo 2. GROMACS as a black box. (PME) 3. How does GROMACS spend its time? (MPE) 4. What you can do What I want to
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 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 informationCoupling the Level-Set Method with Variational Implicit Solvent Modeling of Molecular Solvation
Coupling the Level-Set Method with Variational Implicit Solvent Modeling of Molecular Solvation Bo Li Math Dept & CTBP, UCSD Li-Tien Cheng (Math, UCSD) Zhongming Wang (Math & Biochem, UCSD) Yang Xie (MAE,
More informationNeighbor Tables Long-Range Potentials
Neighbor Tables Long-Range Potentials Today we learn how we can handle long range potentials. Neighbor tables Long-range potential Ewald sums MSE485/PHY466/CSE485 1 Periodic distances Minimum Image Convention:
More informationBiochemistry,530:,, Introduc5on,to,Structural,Biology, Autumn,Quarter,2015,
Biochemistry,530:,, Introduc5on,to,Structural,Biology, Autumn,Quarter,2015, Course,Informa5on, BIOC%530% GraduateAlevel,discussion,of,the,structure,,func5on,,and,chemistry,of,proteins,and, nucleic,acids,,control,of,enzyma5c,reac5ons.,please,see,the,course,syllabus,and,
More informationMOLECULAR DYNAMICS SIMULATIONS OF BIOMOLECULES: Long-Range Electrostatic Effects
Annu. Rev. Biophys. Biomol. Struct. 1999. 28:155 79 Copyright c 1999 by Annual Reviews. All rights reserved MOLECULAR DYNAMICS SIMULATIONS OF BIOMOLECULES: Long-Range Electrostatic Effects Celeste Sagui
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 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 informationThe AGBNP2 Implicit Solvation Model
The AGBNP2 Implicit Solvation Model arxiv:0903.3058v1 [q-bio.bm] 17 Mar 2009 Emilio Gallicchio, Kristina Paris, and Ronald M. Levy Department of Chemistry and Chemical Biology and BioMaPS Institute for
More informationTitle Theory of solutions in the energy r of the molecular flexibility Author(s) Matubayasi, N; Nakahara, M Citation JOURNAL OF CHEMICAL PHYSICS (2003), 9702 Issue Date 2003-11-08 URL http://hdl.handle.net/2433/50354
More informationGeneralized Born Model with a Simple, Robust Molecular Volume Correction
156 J. Chem. Theory Comput. 007, 3, 156-169 Generalized Born Model with a Simple, Robust Molecular Volume Correction John Mongan,,, Carlos Simmerling, # J. Andrew McCammon,,,,,,3 David A. Case,,f and Alexey
More 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 informationInvestigation of physiochemical interactions in
Investigation of physiochemical interactions in Bulk and interfacial water Aqueous salt solutions (new endeavor) Polypeptides exhibiting a helix-coil transition Aqueous globular proteins Protein-solvent
More informationTutorial on the smooth particle-mesh Ewald algorithm
Tutorial on the smooth particle-mesh Ewald algorithm These notes explain the smooth particle-mesh Ewald algorithm (SPME) in detail. The accompanying library libpme6 is a fully-featured implementation of
More informationImproving of the accuracy and efficiency of implicit solvent models in Biomolecular Modeling
Improving of the accuracy and efficiency of implicit solvent models in Biomolecular Modeling Boris Aguilar Huacan Dissertation submitted to the Faculty of the Virginia Polytechnic Institute and State University
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 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 informationMolecular Dynamics Simulations
Molecular Dynamics Simulations Dr. Kasra Momeni www.knanosys.com Outline Long-range Interactions Ewald Sum Fast Multipole Method Spherically Truncated Coulombic Potential Speeding up Calculations SPaSM
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 informationA generalized Born formalism for heterogeneous dielectric environments: Application to the implicit modeling of biological membranes
THE JOURNAL OF CHEMICAL PHYSICS 122, 124706 2005 A generalized Born formalism for heterogeneous dielectric environments: Application to the implicit modeling of biological membranes Seiichiro Tanizaki
More informationForce fields in computer simulation of soft nanomaterials
Lecture 2 Part B Force fields in computer simulation of soft nanomaterials Recommended reading: Leach Chapter 4 1 Force Field Methods Force field methods are simulation methods that use classical force
More informationModeling in Solution. Important Observations. Energy Concepts
Modeling in Solution 1 Important Observations Various electrostatic effects for a solvated molecule are often less important than for an isolated gaseous molecule when the molecule is dissolved in a solvent
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 informationLevel-Set Variational Solvation Coupling Solute Molecular Mechanics with Continuum Solvent
Level-Set Variational Solvation Coupling Solute Molecular Mechanics with Continuum Solvent Bo Li Department of Mathematics and Center for Theoretical Biological Physics (CTBP) University of California,
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 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 informationEfficient multiple time step method for use with Ewald and particle mesh Ewald for large biomolecular systems
JOURNAL OF CHEMICAL PHYSICS VOLUME 115, NUMBER 5 1 AUGUST 2001 Efficient multiple time step method for use with Ewald and particle mesh Ewald for large biomolecular systems Ruhong Zhou IBM Thomas J. Watson
More informationWhat is Classical Molecular Dynamics?
What is Classical Molecular Dynamics? Simulation of explicit particles (atoms, ions,... ) Particles interact via relatively simple analytical potential functions Newton s equations of motion are integrated
More informationMultiple time step Monte Carlo simulations: Application to charged systems with Ewald summation
JOURNAL OF CHEMICAL PHYSICS VOLUME 11, NUMBER 1 1 JULY 004 Multiple time step Monte Carlo simulations: Application to charged systems with Ewald summation Katarzyna Bernacki a) Department of Chemistry
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 informationChallenges in generation of conformational ensembles for peptides and small proteins
Challenges in generation of conformational ensembles for peptides and small proteins Carlos Simmerling Stony Brook University What could (and does) go wrong? 1. Sampling: difficult to obtain converged
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 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 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 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 informationMolecular Dynamics Simulations. Dr. Noelia Faginas Lago Dipartimento di Chimica,Biologia e Biotecnologie Università di Perugia
Molecular Dynamics Simulations Dr. Noelia Faginas Lago Dipartimento di Chimica,Biologia e Biotecnologie Università di Perugia 1 An Introduction to Molecular Dynamics Simulations Macroscopic properties
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 informationStatistical Theory and Learning from Molecular Simulations
Statistical Theory and Learning from Molecular Simulations Lawrence R. Pratt 1 and Susan B. Rempe 2 1 Department of Chemical & Biomolecular Engineering 2 Center for Biological and Material Sciences, Sandia
More informationAlek Marenich Cramer -Truhlar Group. University of Minnesota
and Computational Electrochemistry Alek Marenich Cramer -Truhlar Group University of Minnesota Outline Solvation and charge models Energetics and dynamics in solution Computational electrochemistry Future
More informationLevel-Set Variational Implicit-Solvent Modeling of Biomolecular Solvation
Level-Set Variational Implicit-Solvent Modeling of Biomolecular Solvation Bo Li Department of Mathematics and Quantitative Biology Graduate Program UC San Diego The 7 th International Congress of Chinese
More informationT6.2 Molecular Mechanics
T6.2 Molecular Mechanics We have seen that Benson group additivities are capable of giving heats of formation of molecules with accuracies comparable to those of the best ab initio procedures. However,
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 informationFast Analytical Methods for Macroscopic Electrostatic Models in Biomolecular Simulations
SIAM REVIEW Vol. 53, No. 4, pp. 683 720 c 20 Society for Industrial and Applied Mathematics Fast Analytical Methods for Macroscopic Electrostatic Models in Biomolecular Simulations Zhenli Xu Wei Cai Abstract.
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 informationThe change in free energy on transferring an ion from a medium of low dielectric constantε1 to one of high dielectric constant ε2:
The Born Energy of an Ion The free energy density of an electric field E arising from a charge is ½(ε 0 ε E 2 ) per unit volume Integrating the energy density of an ion over all of space = Born energy:
More informationSolvation and reorganization energies in polarizable molecular and continuum solvents
Solvation and reorganization energies in polarizable molecular and continuum solvents Joel S. Bader CuraGen Corporation, 322 East Main Street, Branford, Connecticut 06405 Christian M. Cortis Department
More informationSimulations of bulk phases. Periodic boundaries. Cubic boxes
Simulations of bulk phases ChE210D Today's lecture: considerations for setting up and running simulations of bulk, isotropic phases (e.g., liquids and gases) Periodic boundaries Cubic boxes In simulations
More informationA Fast N-Body Solver for the Poisson(-Boltzmann) Equation
A Fast N-Body Solver for the Poisson(-Boltzmann) Equation Robert D. Skeel Departments of Computer Science (and Mathematics) Purdue University http://bionum.cs.purdue.edu/2008december.pdf 1 Thesis Poisson(-Boltzmann)
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 informationSemi Empirical Force Fields and Their Limitations. Potential Energy Surface (PES)
Semi Empirical Force Fields and Their Limitations Ioan Kosztin Beckman Institute University of Illinois at Urbana-Champaign Potential Energy Surface (PES) Schrödinger equation: H T Ψ( r, = E Ψ( r, H =
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 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 informationLong-range interactions: P 3 M, MMMxD, ELC, MEMD and ICC
Long-range interactions: P 3 M, MMMxD, ELC, MEMD and ICC Axel Arnold Institute for Computational Physics Universität Stuttgart October 10, 2012 Long-range interactions gravity Coulomb interaction dipolar
More informationAccounting for Solvation in Quantum Chemistry. Comp chem spring school 2014 CSC, Finland
Accounting for Solvation in Quantum Chemistry Comp chem spring school 2014 CSC, Finland In solution or in a vacuum? Solvent description is important when: Polar solvent: electrostatic stabilization Solvent
More informationBiomolecules are dynamic no single structure is a perfect model
Molecular Dynamics Simulations of Biomolecules References: A. R. Leach Molecular Modeling Principles and Applications Prentice Hall, 2001. M. P. Allen and D. J. Tildesley "Computer Simulation of Liquids",
More informationFree Energy Landscape of Protein Folding in Water: Explicit vs. Implicit Solvent
PROTEINS: Structure, Function, and Genetics 53:148 161 (2003) Free Energy Landscape of Protein Folding in Water: Explicit vs. Implicit Solvent Ruhong Zhou* IBM T.J. Watson Research Center, Yorktown Heights,
More informationElectrolytes. Chapter Basics = = 131 2[ ]. (c) From both of the above = = 120 8[
Chapter 1 Electrolytes 1.1 Basics Here we consider species that dissociate into positively and negatively charged species in solution. 1. Consider: 1 H (g) + 1 Cl (g) + ()+ () = { } = (+ )+ ( ) = 167[
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 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 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 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 informationTheory and Modeling of Specific Ion Hydration
Theory and Modeling of Specific Ion Hydration Thomas L. Beck, University of Cincinnati November 21, 2013 Thanks: Lawrence Pratt, Hank Ashbaugh, Dilip Asthagiri, PNNL group, Kevin Leung, John Weeks, Barry
More informationCoarse-Grained Models!
Coarse-Grained Models! Large and complex molecules (e.g. long polymers) can not be simulated on the all-atom level! Requires coarse-graining of the model! Coarse-grained models are usually also particles
More informationHYBRID SOLVATION MODELS FOR ELECTROSTATIC INTERACTIONS IN MOLECULAR DYNAMICS SIMULATIONS OF IONIC SOLVENT. Ming Xiang
HYBRID SOLVATION MODELS FOR ELECTROSTATIC INTERACTIONS IN MOLECULAR DYNAMICS SIMULATIONS OF IONIC SOLVENT by Ming Xiang A dissertation submitted to the faculty of The University of North Carolina at Charlotte
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 informationMultiscale Materials Modeling
Multiscale Materials Modeling Lecture 09 Quantum Mechanics/Molecular Mechanics (QM/MM) Techniques Fundamentals of Sustainable Technology These notes created by David Keffer, University of Tennessee, Knoxville,
More informationMCB100A/Chem130 MidTerm Exam 2 April 4, 2013
MCBA/Chem Miderm Exam 2 April 4, 2 Name Student ID rue/false (2 points each).. he Boltzmann constant, k b sets the energy scale for observing energy microstates 2. Atoms with favorable electronic configurations
More informationExercise 2: Solvating the Structure Before you continue, follow these steps: Setting up Periodic Boundary Conditions
Exercise 2: Solvating the Structure HyperChem lets you place a molecular system in a periodic box of water molecules to simulate behavior in aqueous solution, as in a biological system. In this exercise,
More informationRapid evaluation of electrostatic interactions in multi-phase media
Rapid evaluation of electrostatic interactions in multi-phase media P.G. Martinsson, The University of Colorado at Boulder Acknowledgements: Some of the work presented is joint work with Mark Tygert and
More information