Lecture 19: Free Energies in Modern Computational Statistical Thermodynamics: FEP and Related Methods
|
|
- Tamsyn Lawson
- 5 years ago
- Views:
Transcription
1 Statistical Thermodynamics Lecture 19: Free Energies in Modern Computational Statistical Thermodynamics: FEP and Related Methods Dr. Ronald M. Levy
2 Free energy calculations Free energy is the most important quantity that describes a thermodynamic process Two types of free energy calculations: 1. Path independent methods for calculation of free energies. e.g. Free Energy Perturbation (FEP), Thermodynamic Integration (TI) 2. Path dependent methods for calculations of free energies. e.g. Potential of Mean Force (PMF)
3 General perturbation theory Perturbation theory is one of the oldest and most useful, general techniques in applied mathematics. Its initial applications to physics were in celestial mechanics to explain how the presence of bodies other than the sun perturbed the elliptical orbits of planets. The main idea: One starts with an initial system, called the unperturbed or reference system. The system of interest, called the target system, is represented in terms of a perturbation to the reference system.
4 Statistical mechanics description of free energy In the canonical ensemble: A = k B T lnq N-particle reference system described by Hamiltonian H(x,p), function of 3N Cartesian coordinates x and conjugated momenta p. Partition function: Q = Free energy difference : ΔA = k B T ln 1 h 3N N! exp[ βh(x, p)]dx dp ΔA = A 1 A 0 = k B T ln Q 1 Q 0 exp[ βh 1 (x, p)] dxdp exp[ βh 0 (x, p)]dx dp
5 ΔA = k B T ln exp[ βh 1 (x, p)] dxdp exp[ βh 0 (x, p)]dx dp = k B T ln exp[ βh 1 (x, p)+ βh 0 (x, p)] exp[ βh 0 (x, p)]dxdp exp[ βh 0 (x, p)]dx dp = k B T ln exp[ βδh(x, p)] exp[ βh 0 (x, p)]dxdp exp[ βh 0 (x, p)]dx dp
6 Probability density function P 0 of finding the reference system in a state defined by x and p is: P 0 (x, p) = Substituting in expression of ΔA: ΔA = k B T ln exp[ βh 0 (x, p)] exp[ βh 0 (x, p)]dx dp exp[ βδh(x, p)] exp[ βh 0 (x, p)]dxdp exp[ βh 0 (x, p)]dx dp = k B T ln exp[ βδh(x, p)] P 0 (x, p)dxdp Fundamental FEP formula for the transformation 0 1: ΔA = k B T ln < exp[ βδh(x, p)] > 0 If the kinetic term cancels out : ΔA = k B T ln < exp[ βδu] > 0
7 Free energy perturbation method ΔA = k B T ln < exp[ βδu] > 0 In this formula, ΔU = U 1 (x) U 0 (x) is the difference in potential between target and reference system, and the average is over the ensemble of the initial state corresponding to reference system with potential U 0 (x). Similarly, the free energy difference can also be written in terms of an average over the ensemble of the final state corresponding to target system with potential U 1 (x): ΔA = A 1 A 0 = k B T ln Q 1 Q 0 = k B T ln Q 0 Q 1 = k B T ln < exp[βδu] > 1
8 Free energy perturbation calculations continue ΔA = k B T ln ΔA = k B T ln exp[ βδh(x, p)] P 0 (x, p)dxdp exp( βδu)p 0 (ΔU)dΔU According to the central limit theorem (CLT), ΔU would be Gaussian-distributed if U 1 and U 0 were the functions of a sufficient number of identically distributed random variables. Therefore: where: P 0 (ΔU) = 1 2πσ exp[ (ΔU < ΔU > 0 )2 2σ 2 ] σ 2 =< ΔU 2 > 0 < ΔU > 0 2 Analytical expression of ΔA: ΔA =< ΔU > βσ 2 Notice: use of this analytical expression can be successful only if P 0 (ΔU) is a narrow function of ΔU; this condition does not imply that the ΔA must be small.
9 Example: hydration free energy of benzene reference state: solute has high probability to overlap with solvent in the gas phase target state: solute has low probability to overlap with solvent in the aqueous phase ΔA = kcal/mol Although hydration free energy of benzene is small, this quantity cannot be successfully calculated directly using the FEP equation; because low-energy configurations in the target ensemble, which do not suffer from the overlap between the solute and solvent molecules, are not sampled in simulations of the reference state.
10 FEP calculation --- staging The difficulty in applying FEP formula can be circumvented through staging strategy. Construction of several intermediate states between reference and target state such that P(ΔU i,i+1 ) for two consecutive states i and i+1 sampled at state i is sufficiently narrow for the direct evaluation of the corresponding free energy difference ΔA i,i+1. With N-2 intermediate states, N 1 i=1 ΔA = ΔA i,i+1 = 1 β N 1 i=1 ln < exp( βδu i,i+1 ) > Intermediate states do not need to be physically meaningful, they do not have to correspond to systems that actually exist. i
11 FEP calculation --- formula The Hamiltonian can be considered to be a function of some parameter λ. λ can be defined between 0 and 1, such that λ = 0 for reference state and λ = 1 for target state. Dependence of hybrid Hamiltonian on λ : H(λ) = (1 λ)h 0 + λh 1 = H 0 + λδh ΔH is the perturbation Hamiltonian, equal to H 1 -H 0. If N-2 intermediate states are created to link the reference and target states, such that λ 1 = 0 and λ N = 1: ΔH i = H(λ i+1 ) H(λ i ) = (λ i+1 λ i )ΔH Total free energy difference: = Δλ i ΔH with Δλ i =λ i+1 -λ i ΔU i = Δλ i ΔU ΔA = 1 β N 1 i=1 ln < exp( βδλ i ΔU) > λ i
12 Important concept in FEP: Order parameter Order parameter: They are collective variables used to describe transformations from the reference system to the target one. order parameter: distance order parameter: dihedral An order parameter may or may not correspond to the path along which the transformation takes place in nature, and would be called the reaction coordinate if such were the case. There is more than one way to define an order parameter. The choice of order parameters may have a significant effect on the efficiency and accuracy of free energy calculations.
13 Application of FEP: binding free energy calculation The binding free energy of two molecules, ΔA binding, defined as the free energy difference between two molecules in the bound and free, unbound states, can be determined experimentally through the measurement of binding constants. For direct calculation of ΔA binding, it needs to define an order parameter that measures the separation between the ligand and the binding center of the protein. This may be difficult when the ligand is buried deep in the binding pocket. Alternative route: double annihilation of the ligand through a thermodynamic cycle.
14 Thermodynamic cycle and double annihilation ΔA binding = ΔA 0 annihilation ΔA1 annihilation + ΔA restrain Double annihilation involves the removal of the ligand molecule from aqueous solution in the unbound state and bound state, respectively. Through thermodynamic cycle, the double annihilation free energies ΔA 0 annihilation and ΔA 1 annihilation were related to ΔA binding The ligand in the binding pocket is annihilated from a strongly constrained position, whereas the unbound ligand can translate and rotate freely during annihilation. Proper corrections for the loss of translational, rotational and conformational entropies should be taken into account. ΔA binding protein + ligand protein ligand ΔA 0 annihilation ΔA1 annihilation ΔA restrain protein + nothing protein nothing
15 Example: Relative binding free energies for ligands In the case of seeking potential inhibitors of a target protein, determining relative binding free energies for a series of ligands is required. This can be handled by repeating the absolute binding free energy calculation for each ligand of interest. ΔA XLC = kcal/mol Ligand XLC Protein FXa ΔA XLD = kcal/mol Ligand XLD Protein FXa ΔΔA cal = ΔA XLD ΔA XLC = kcal/mol. Match well with ΔΔA exp = kcal/mol
16 Alchemical transformations ΔA A binding protein + ligand A protein ligand ΔA 0 mutation ΔA1 mutation ΔA B binding protein + ligand B protein ligand There is an alternate pathway for calculations of relative binding free energies. Mutation of ligand A into ligand B in both the bound and the free states, following a different thermodynamic cycle. Example: In the mutation of ethane into methanol, the former serves as the common topology. As the carbon atom is transformed into oxygen, two hydrogen atoms of the methyl moiety are turned into non-interacting, ghost particles by annihilating their point charges and van der Waals parameters.
17 Thermodynamic integration method If N-2 intermediate states are created to link the reference and target states, such that λ 1 = 0 and λ N = 1: ΔA = A 1 (λ N ) A 0 (λ 1 ) = λ da(λ) N dλ = dλ λ da(λ) dλ dλ Q = da(λ) dλ exp[ βu(x, λ)] dx A = kt lnq = kt 1 Q d( exp[ βu(x, λ)]dx) dλ = kt β = du dλ λ du dλ exp( βu)dx Q TI: ΔA = 1 0 du dλ λ dλ FEP: ΔA = k B T ln < exp[ βδu] > 0
18 Potential of mean force method The reaction coordinate ξ, describes the path along which the transformation takes place in nature. Potential of mean force (PMF) is a path dependent method for calculations of free energies based on the reaction coordinate ξ. ΔA = A 1 (ξ N ) A 0 (ξ 1 ) = da(ξ) dξ = dξ du dξ du dξ ξ = F ξ ξ = constraint force Constraint force: force required to constrain the system at a fix reaction coordinate ξ; if this force can be extracted, free energy can be evaluated from integration of the mean force along the reaction coordinate. ξ dξ
19 Example: Single ion at liquid/vapor interface Z direction ξ Zero Position Gibbs Dividing Surface The reaction coordinate is the separation between ion and zero position along the z direction. Iodide shows more stability than chloride around interface from PMF profile.
20 Summary Different method for calculation free energy: free energy perturbation (FEP), thermodynamic integration (TI) and potential of mean force (PMF). In the next lab, we will use FEP method, implemented in the Schrodinger package, to calculate the relative binding free energy between a pair of ligand XLC and XLD binding to the receptor protein Fxa.
Free 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 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 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 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 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 informationLecture 14: Advanced Conformational Sampling
Lecture 14: Advanced Conformational Sampling Dr. Ronald M. Levy ronlevy@temple.edu Multidimensional Rough Energy Landscapes MD ~ ns, conformational motion in macromolecules ~µs to sec Interconversions
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 informationFree Energy Calculations in Biological Systems. How Useful Are They in Practice?
Free Energy Calculations in Biological Systems. How Useful Are They in Practice? Christophe Chipot Equipe de dynamique des assemblages membranaires, UMR CNRS/UHP 7565, Institut nancéien de chimie moléculaire,
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 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 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 information1. Thermodynamics 1.1. A macroscopic view of matter
1. Thermodynamics 1.1. A macroscopic view of matter Intensive: independent of the amount of substance, e.g. temperature,pressure. Extensive: depends on the amount of substance, e.g. internal energy, enthalpy.
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 informationGases and the Virial Expansion
Gases and the irial Expansion February 7, 3 First task is to examine what ensemble theory tells us about simple systems via the thermodynamic connection Calculate thermodynamic quantities: average energy,
More informationG : Statistical Mechanics
G25.2651: Statistical Mechanics Notes for Lecture 9 I. OVERVIEW Our treatment of the classical ensembles makes clear that the free energy is a quantity of particular importance in statistical mechanics.
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, 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 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 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 informationCHEM-UA 652: Thermodynamics and Kinetics
1 CHEM-UA 652: Thermodynamics and Kinetics Notes for Lecture 11 I. PHYSICAL AND CHEMICAL RELEVANCE OF FREE ENERGY In this section, we will consider some examples showing the significance of free energies.
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 informationψ s a ˆn a s b ˆn b ψ Hint: Because the state is spherically symmetric the answer can depend only on the angle between the two directions.
1. Quantum Mechanics (Fall 2004) Two spin-half particles are in a state with total spin zero. Let ˆn a and ˆn b be unit vectors in two arbitrary directions. Calculate the expectation value of the product
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 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 information3. RATE LAW AND STOICHIOMETRY
Page 1 of 39 3. RATE LAW AND STOICHIOMETRY Professional Reference Shelf R3.2 Abbreviated Lecture Notes Full Lecture Notes I. Overview II. Introduction A. The Transition State B. Procedure to Calculate
More informationPHYS 352 Homework 2 Solutions
PHYS 352 Homework 2 Solutions Aaron Mowitz (, 2, and 3) and Nachi Stern (4 and 5) Problem The purpose of doing a Legendre transform is to change a function of one or more variables into a function of variables
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 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 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 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 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 informationPhysics 4230 Final Examination 10 May 2007
Physics 43 Final Examination May 7 In each problem, be sure to give the reasoning for your answer and define any variables you create. If you use a general formula, state that formula clearly before manipulating
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 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 informationHyeyoung Shin a, Tod A. Pascal ab, William A. Goddard III abc*, and Hyungjun Kim a* Korea
The Scaled Effective Solvent Method for Predicting the Equilibrium Ensemble of Structures with Analysis of Thermodynamic Properties of Amorphous Polyethylene Glycol-Water Mixtures Hyeyoung Shin a, Tod
More informationBiomolecular modeling III
2016, January 5 Déjà vu Enhanced sampling Biomolecular simulation Each atom x, y, z coordinates Déjà vu Enhanced sampling Expression for energy the force field = 1 2 + N i i k i (r i r 0 i ) 2 + 1 2 N
More informationStatistical physics models belonging to the generalised exponential family
Statistical physics models belonging to the generalised exponential family Jan Naudts Universiteit Antwerpen 1. The generalized exponential family 6. The porous media equation 2. Theorem 7. The microcanonical
More informationChapter 6- An Introduction to Metabolism*
Chapter 6- An Introduction to Metabolism* *Lecture notes are to be used as a study guide only and do not represent the comprehensive information you will need to know for the exams. The Energy of Life
More informationElectronic Supplementary Information Effective lead optimization targeted for displacing bridging water molecule
Electronic Supplementary Material (ESI) for Physical Chemistry Chemical Physics. This journal is the Owner Societies 2018 Electronic Supplementary Information Effective lead optimization targeted for displacing
More informationLecture 27. Transition States and Enzyme Catalysis
Lecture 27 Transition States and Enzyme Catalysis Reading for Today: Chapter 15 sections B and C Chapter 16 next two lectures 4/8/16 1 Pop Question 9 Binding data for your thesis protein (YTP), binding
More informationAdvanced sampling. fluids of strongly orientation-dependent interactions (e.g., dipoles, hydrogen bonds)
Advanced sampling ChE210D Today's lecture: methods for facilitating equilibration and sampling in complex, frustrated, or slow-evolving systems Difficult-to-simulate systems Practically speaking, one is
More informationComputational Predictions of 1-Octanol/Water Partition Coefficient for Imidazolium based Ionic Liquids.
Computational Predictions of 1-Octanol/Water Partition Coefficient for Imidazolium based Ionic Liquids. Ganesh Kamath,* a Navendu Bhatnagar b, Gary A. Baker a, Sheila N. Baker c and Jeffrey J. Potoff b
More informationThe Partition Function Statistical Thermodynamics. NC State University
Chemistry 431 Lecture 4 The Partition Function Statistical Thermodynamics NC State University Molecular Partition Functions In general, g j is the degeneracy, ε j is the energy: = j q g e βε j We assume
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 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 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 informationExploring the Free Energy Surface of Short Peptides by Using Metadynamics
John von Neumann Institute for Computing Exploring the Free Energy Surface of Short Peptides by Using Metadynamics C. Camilloni, A. De Simone published in From Computational Biophysics to Systems Biology
More information2.4 DNA structure. S(l) bl + c log l + d, with c 1.8k B. (2.72)
2.4 DNA structure DNA molecules come in a wide range of length scales, from roughly 50,000 monomers in a λ-phage, 6 0 9 for human, to 9 0 0 nucleotides in the lily. The latter would be around thirty meters
More informationGAS-SURFACE INTERACTIONS
Page 1 of 16 GAS-SURFACE INTERACTIONS In modern surface science, important technological processes often involve the adsorption of molecules in gaseous form onto a surface. We can treat this adsorption
More informationCHEM 344 Molecular Modeling
CHEM 344 Molecular Modeling The Use of Computational Chemistry to Support Experimental Organic Chemistry Part 1: Molecular Orbital Theory, Hybridization, & Formal Charge * all calculation data obtained
More informationUnderstanding Molecular Simulation 2009 Monte Carlo and Molecular Dynamics in different ensembles. Srikanth Sastry
JNCASR August 20, 21 2009 Understanding Molecular Simulation 2009 Monte Carlo and Molecular Dynamics in different ensembles Srikanth Sastry Jawaharlal Nehru Centre for Advanced Scientific Research, Bangalore
More informationWashington University in St. Louis Chemistry Tournament Sample Problems for Individual Round #3: Kinetics, Electrochemistry, and Thermodynamics
Individual Exam #3: Kinetics, Electrochemistry, and Thermodynamics 1) A galvanic cell is to be constructed using the following half reactions under acidic conditions: ε 1.49 V ε 1.66 V Given that the potentials
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 informationVirtual screening for drug discovery. Markus Lill Purdue University
Virtual screening for drug discovery Markus Lill Purdue University mlill@purdue.edu Lecture material http://people.pharmacy.purdue.edu/~mlill/teaching/eidelberg/ I.1 Drug discovery Cl N Disease I.1 Drug
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 information2m + U( q i), (IV.26) i=1
I.D The Ideal Gas As discussed in chapter II, micro-states of a gas of N particles correspond to points { p i, q i }, in the 6N-dimensional phase space. Ignoring the potential energy of interactions, the
More informationWhy Proteins Fold. How Proteins Fold? e - ΔG/kT. Protein Folding, Nonbonding Forces, and Free Energy
Why Proteins Fold Proteins are the action superheroes of the body. As enzymes, they make reactions go a million times faster. As versatile transport vehicles, they carry oxygen and antibodies to fight
More information3.320 Lecture 18 (4/12/05)
3.320 Lecture 18 (4/12/05) Monte Carlo Simulation II and free energies Figure by MIT OCW. General Statistical Mechanics References D. Chandler, Introduction to Modern Statistical Mechanics D.A. McQuarrie,
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 informationChemical Potential. Combining the First and Second Laws for a closed system, Considering (extensive properties)
Chemical Potential Combining the First and Second Laws for a closed system, Considering (extensive properties) du = TdS pdv Hence For an open system, that is, one that can gain or lose mass, U will also
More informationMCB100A/Chem130 MidTerm Exam 2 April 4, 2013
MCB1A/Chem13 MidTerm Exam 2 April 4, 213 Name Student ID True/False (2 points each). 1. The Boltzmann constant, k b T sets the energy scale for observing energy microstates 2. Atoms with favorable electronic
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 informationThe protein folding problem consists of two parts:
Energetics and kinetics of protein folding The protein folding problem consists of two parts: 1)Creating a stable, well-defined structure that is significantly more stable than all other possible structures.
More 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 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 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 informationMolecular Modeling Lecture 11 side chain modeling rotamers rotamer explorer buried cavities.
Molecular Modeling 218 Lecture 11 side chain modeling rotamers rotamer explorer buried cavities. Sidechain Rotamers Discrete approximation of the continuous space of backbone angles. Sidechain conformations
More informationrate of reaction forward conc. reverse time P time Chemical Equilibrium Introduction Dynamic Equilibrium Dynamic Equilibrium + RT ln f p
Chemical Equilibrium Chapter 9 of Atkins: Sections 9.1-9.2 Spontaneous Chemical Reactions The Gibbs Energy Minimum The reaction Gibbs energy Exergonic and endergonic reactions The Description of Equilibrium
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 informationMolecular Driving Forces
Molecular Driving Forces Statistical Thermodynamics in Chemistry and Biology SUBGfittingen 7 At 216 513 073 / / Ken A. Dill Sarina Bromberg With the assistance of Dirk Stigter on the Electrostatics chapters
More informationThus, the volume element remains the same as required. With this transformation, the amiltonian becomes = p i m i + U(r 1 ; :::; r N ) = and the canon
G5.651: Statistical Mechanics Notes for Lecture 5 From the classical virial theorem I. TEMPERATURE AND PRESSURE ESTIMATORS hx i x j i = kt ij we arrived at the equipartition theorem: * + p i = m i NkT
More informationA Brief Introduction to Statistical Mechanics
A Brief Introduction to Statistical Mechanics E. J. Maginn, J. K. Shah Department of Chemical and Biomolecular Engineering University of Notre Dame Notre Dame, IN 46556 USA Monte Carlo Workshop Universidade
More informationNon-degenerate Perturbation Theory. and where one knows the eigenfunctions and eigenvalues of
on-degenerate Perturbation Theory Suppose one wants to solve the eigenvalue problem ĤΦ = Φ where µ =,1,2,, E µ µ µ and where Ĥ can be written as the sum of two terms, ˆ ˆ ˆ ˆ ˆ ˆ H = H + ( H H ) = H +
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 informationEntropy enthalpy compensation: Fact or artifact?
FOR THE RECORD Entropy enthalpy compensation: Fact or artifact? KIM SHARP Johnson Research Foundation, Department of Biochemistry and Biophysics, University of Pennsylvania, Philadelphia, Pennsylvania
More informationObjectives INTRODUCTION TO METABOLISM. Metabolism. Catabolic Pathways. Anabolic Pathways 3/6/2011. How to Read a Chemical Equation
Objectives INTRODUCTION TO METABOLISM. Chapter 8 Metabolism, Energy, and Life Explain the role of catabolic and anabolic pathways in cell metabolism Distinguish between kinetic and potential energy Distinguish
More informationProtein Folding. I. Characteristics of proteins. C α
I. Characteristics of proteins Protein Folding 1. Proteins are one of the most important molecules of life. They perform numerous functions, from storing oxygen in tissues or transporting it in a blood
More informationORGANIC - EGE 5E CH. 2 - COVALENT BONDING AND CHEMICAL REACTIVITY
!! www.clutchprep.com CONCEPT: HYBRID ORBITAL THEORY The Aufbau Principle states that electrons fill orbitals in order of increasing energy. If carbon has only two unfilled orbitals, why does it like to
More informationFundamentals of Distribution Separations (III)
Fundamentals of Distribution Separations (III) (01/16/15) K = exp -Δμ 0 ext i - Δμ i RT distribution coefficient C i = exp -Δμ 0 RT - - Δμ i = ΔH i TΔS i 0 0 0 solubility q A---B A + B 0 0 0 ΔH i = ΔH
More informationSupplemental Material: Solute Adsorption at Air-Water Interfaces and Induced. Interface Fluctuations: The Hydrophobic Nature of Ions?
Electronic Supplementary Material (ESI) for Physical Chemistry Chemical Physics. This journal is the Owner Societies Supplemental Material: Solute Adsorption at Air-Water Interfaces and Induced Interface
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 informationClose agreement between the orientation dependence of hydrogen bonds observed in protein structures and quantum mechanical calculations
Close agreement between the orientation dependence of hydrogen bonds observed in protein structures and quantum mechanical calculations Alexandre V. Morozov, Tanja Kortemme, Kiril Tsemekhman, David Baker
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 informationLECTURE 2 STRUCTURE AND PROPERTIES OF ORGANIC MOLECULES
LECTURE 2 STRUCTURE AND PROPERTIES OF ORGANIC MOLECULES 1. Atomic wave functions and orbitals. LCAO. The important thing to know is that atomic orbitals are represented by wave functions, and they have
More informationChapter 2 Approximation Methods Can be Used When Exact Solutions to the Schrödinger Equation Can Not be Found.
Chapter 2 Approximation Methods Can be Used When Exact Solutions to the Schrödinger Equation Can Not be Found. In applying quantum mechanics to 'real' chemical problems, one is usually faced with a Schrödinger
More information4.1 LAWS OF MECHANICS - Review
4.1 LAWS OF MECHANICS - Review Ch4 9 SYSTEM System: Moving Fluid Definitions: System is defined as an arbitrary quantity of mass of fixed identity. Surrounding is everything external to this system. Boundary
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 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 informationLattice protein models
Lattice protein models Marc R. Roussel epartment of Chemistry and Biochemistry University of Lethbridge March 5, 2009 1 Model and assumptions The ideas developed in the last few lectures can be applied
More informationThermal and Statistical Physics Department Exam Last updated November 4, L π
Thermal and Statistical Physics Department Exam Last updated November 4, 013 1. a. Define the chemical potential µ. Show that two systems are in diffusive equilibrium if µ 1 =µ. You may start with F =
More informationChapter 04 Alcohols and Alkyl Halides part 01
hapter 04 Alcohols and Alkyl alides part 01 EM 341: Spring 2012 Prof. Greg ook Functional Groups A functional group is a structural feature in a molecule that has characteristic reactivity. A functional
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 information2. Thermodynamics. Introduction. Understanding Molecular Simulation
2. Thermodynamics Introduction Molecular Simulations Molecular dynamics: solve equations of motion r 1 r 2 r n Monte Carlo: importance sampling r 1 r 2 r n How do we know our simulation is correct? Molecular
More informationChap. 2. Polymers Introduction. - Polymers: synthetic materials <--> natural materials
Chap. 2. Polymers 2.1. Introduction - Polymers: synthetic materials natural materials no gas phase, not simple liquid (much more viscous), not perfectly crystalline, etc 2.3. Polymer Chain Conformation
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 informationEntropy online activity: accompanying handout I. Introduction
Entropy online activity: accompanying handout I. Introduction The goal of this activity is to provide insight into the ways modern science views the effects of temperature on chemical reactions, including
More informationFirst Law CML 100, IIT Delhi SS. The total energy of the system. Contribution from translation + rotation + vibrations.
Internal Energy he total energy of the system. Contribution from translation + rotation + vibrations. Equipartition theorem for the translation and rotational degrees of freedom. 1/ k B Work Path function,
More information