Computational Studies of the Photoreceptor Rhodopsin. Scott E. Feller Wabash College

Size: px
Start display at page:

Download "Computational Studies of the Photoreceptor Rhodopsin. Scott E. Feller Wabash College"

Transcription

1 Computational Studies of the Photoreceptor Rhodopsin Scott E. Feller Wabash College

2 Rhodopsin Photocycle Dark-adapted Rhodopsin hn Isomerize retinal Photorhodopsin ~200 fs Bathorhodopsin Meta-II ms timescale Meta-I s timescale ns timescale Lumirhodopsin

3 Challenging Problem Requires Multiple Complimentary Techniques Changes in conformation and protonation of protein Changes in conformation and protonation of retinal ligand Changes in conformation of solvating phospholipids?

4 Dark state retinal orientation: Combining 2 H NMR and MD simulation NMR on reconstituted oriented bilayers with labeled retinal methyl groups MD simulations on BlueGene

5 MD simulation reveals fluctuations in retinal torsions Carried out ns simulations of rhodopsin in an explicit lipid bilayer using different random lipid conformations and positions as starting points As expected the conjugated polyene chain is stiff with all dihedrals in cis or trans conformations in all simulations

6 Dynamic Structure of Retinylidene Ligand of Rhodopsin Probed by Molecular Simulations, J. Mol. Biol. 372, MD simulation reveals fluctuations in b-ionone ring orientation Isomerization of C5-C6-C7-C8 dihedral is slow on MD timescale but multiple conformers observed Pooled distribution of conformers from 2.6 s gives excellent agreement with NMR

7 Modeling the Meta-I state MI formation involves isomerization of retinal and may involve changes in protein protonation state in the binding pocket Carried out simulations (> 1 s each) of competing mechanisms suggested by vibrational spectroscopy experiments Counterion switch Yan, E. et al PNAS 2003, 100, 9262 Complex counterion - Lüdeke, S. et al. J. Mol. Biol. 2005, 353, 345

8 Counterion switch Compare MD with 2 H NMR Complex counterion Retinal Counterion Switch Mechanism in Vision Evaluated by Molecular Simulations, J. Am. Chem. Soc 28;

9 b-ionone ring conformation changes upon isomerization of C11=C12 After isomerization of the 11-cis position the C6-C7 dihedral rotates from s-cis to s-trans state Fluctuations to (-) twisted s-cis state occur on the 100 ns scale

10 Summary of Dark State/MI MD simulations of dark state (11-cis, PSB) reproduce NMR observables After forcing isomerization of C11=C12, the MI state (11- trans, PSB) is produced on the s time scale. MD reproduces NMR provided the correct side chain protonation state is chosen. The b-ionone ring has significant mobility, sampling twisted s-cis conformations (+/- 60) in the dark state and rotating to s-trans in MI

11 Challenges Modelling the MII state MII formation occurs on the millisecond time scale, inaccessible to atomistic MD simulation Fortunately crystal structures of the MII state have become available CHARMM style force fields were published only for the PSB, inappropriate for the deprotonated retinal of MII We set out to refine force fields for protonated and deprotonated retinal consistent with each other and with the CHARMM force field for lipids and proteins

12 Retinal force field parameters depend on protonation state Bond lengths are intermediate between single and double bonds in protonated retinal

13 Retinal force field parameters depend on protonation state Torsional bariers are intermediate between single and double bonds in protonated retinal

14 Examples of torsional energy surfaces Large differences are observed near the Schiff base Small differences are observed near the b-ionone ring

15 An surprising result The C6-C7 bond, that defines the orientation of the ring with respect to the polyene chain, goes against the trend Deprotonated surface shows two minima corresponding to +/- twisted s-cis conformers Protonated surface has minima shifted and the appearance of a local minima at s-trans relative energy (kcal/mol) Deprotonated Protonated Lower barriers to rotation in protonated (dark, MI) states C5=C6-C7=C8

16 Proton affinity depends on b-ionone ring orientation Proton affinity (PA) is the difference in absolute gas-phase energies between protonated and neutral species p r o to n a ffin ity (k c a l/m o l) PA = E A- - E HA PA is highest for planar conformations, i.e. planar angles favor the protonated form raising pk A while deviations from planarity make retinal more acidic C 5 = C 6 -C 7 = C 8

17 pk A of Schiff based depends on b-ionone ring orientation pk A can be computed in a thermodynamic cycle that combines gas phase PA with QM calculation of DG solv in implicit solvent Accurate pka calculations for carboxylic acids using complete basis set and Gaussian-n models combined with CPCM continuum solvation methods MD Liptak & GC Shields, JACS Absolute values of pk A likely shift due to environment but conformation dependent shifts of ~2 units is predicted among thermally accessible states Is this relevant for MI MII transition? pk A C5=C6-C7=C8

18 Simulation suggests fluctuations in MI state correspond to large pk A changes Thermally accessible transitions between s-cis and s-trans involve pk A change of ~ 2 units, i.e. a 100 fold change in equilibrium constant

19 Potential implications for Rhodopsin activation Potential for reciprocal control where ligand conformational change (11-cis to all-trans) leads to protein conformational change which subsequently shifts the ligand conformation to one that favors deprotonation Implies that ligand-protein interactions change as a function of the Rhodopsin photocycle

20 Preliminary results on ligand-protein interactions Begin by examining dark state and MII where crystal structures are available Carry out Umbrella Sampling with Hamiltonian Replica exchange using the C5=C6-C7=C8 dihedral as the reaction coordinate First compute potential of mean force in vacuum, repeat in full membrane protein environment Compute DDG between vacuum and protein. Does it change with photocycle?

21 Preliminary results on ligand-protein interactions The effect of the protein is to change the free energy profile by ~1 kcal/mol in the deprotonated (MII) case

22 Preliminary results on ligand-protein interactions The s-trans state is destabilized in MII (but not in dark state) consistent with the protein selecting the ligand conformation that favors deprotonation

23 Preliminary results on ligand-protein interactions We have not computed the PMF for MI but the unbiased simulation suggests that s-trans is stabilized in MI

24 Summary MD simulations, after validation by experimental observables, can provide an atomic level interpretation of the structures and fluctuations in membrane proteins The retinal ligand is a complicated species requiring careful parameterization to obtain accurate MD simulations of rhodopsin Fluctuations in b-ionone ring conformation may provide a mechanism to tune the Schiff base pk A

25 Acknowledgements Mike Lu, Shenghuang Zhu, Lu Hong (Wabash College) Michael Brown (University of Arizona) Blake Mertz (West Virginia University) Alan Grossfield (University of Rochester) NSF MCB

Retinal Proteins (Rhodopsins) Vision, Bioenergetics, Phototaxis. Bacteriorhodopsin s Photocycle. Bacteriorhodopsin s Photocycle

Retinal Proteins (Rhodopsins) Vision, Bioenergetics, Phototaxis. Bacteriorhodopsin s Photocycle. Bacteriorhodopsin s Photocycle Molecular chanisms of Photoactivation and Spectral Tuning in Retinal Proteins Emad Tajkhorshid Theoretical and Computational Biophysics Group Beckman Institute University of Illinois at Urbana-Champaign

More information

Biochimica et Biophysica Acta

Biochimica et Biophysica Acta Biochimica et Biophysica Acta 1818 (2012) 241 251 Contents lists available at ScienceDirect Biochimica et Biophysica Acta journal homepage: www.elsevier.com/locate/bbamem Review Molecular simulations and

More information

The Molecular Dynamics Method

The 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 information

TDDFT as a tool in biophysics

TDDFT as a tool in biophysics TDDFT as a tool in biophysics The primary event in vision Robert Send Universität Karlsruhe 09.09.08 Robert Send TDDFT as a tool in biophysics 09.09.08 1 / 28 Outline 1 Human vision 2 The methods 3 The

More information

A molecular-dynamics study of the rhodopsin chromophore using ultrasoft pseudopotentials

A molecular-dynamics study of the rhodopsin chromophore using ultrasoft pseudopotentials Progress of Theoretical Physics Supplement No. 138, 2000 107 A molecular-dynamics study of the rhodopsin chromophore using ultrasoft pseudopotentials Minoru Sugihara, 1 ) Peter Entel, 1 Hendrik Meyer,

More information

Comparison between Bacteriorhodopsin and Halorhodopsin. Halorhodopsin (HR) and Bacteriorhodopsin (BR) belong to a subfamily of

Comparison between Bacteriorhodopsin and Halorhodopsin. Halorhodopsin (HR) and Bacteriorhodopsin (BR) belong to a subfamily of Comparison between Bacteriorhodopsin and Halorhodopsin Halorhodopsin (HR) and Bacteriorhodopsin (BR) belong to a subfamily of heptahelical membrane proteins, the archaeal rhodopsins. They are found in

More information

An introduction to Molecular Dynamics. EMBO, June 2016

An 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 information

Force Fields for MD simulations

Force Fields for MD simulations Force Fields for MD simulations Topology/parameter files Where do the numbers an MD code uses come from? ow to make topology files for ligands, cofactors, special amino acids, ow to obtain/develop missing

More information

Structural Analysis and Dynamics of Retinal Chromophore in Dark and Meta I States of Rhodopsin from 2 H NMR of Aligned Membranes

Structural Analysis and Dynamics of Retinal Chromophore in Dark and Meta I States of Rhodopsin from 2 H NMR of Aligned Membranes doi:10.1016/j.jmb.2007.03.046 J. Mol. Biol. (2007) 372, 50 66 Structural Analysis and Dynamics of Retinal Chromophore in Dark and Meta I States of Rhodopsin from 2 H NMR of Aligned Membranes Andrey V.

More information

Three-layer ONIOM studies of the dark state of Rhodopsin: the protonation state of Glu181

Three-layer ONIOM studies of the dark state of Rhodopsin: the protonation state of Glu181 Accepted Manuscript Three-layer ONIOM studies of the dark state of Rhodopsin: the protonation state of Glu181 Katherine. F. Hall, Thom Vreven, Michael J. Frisch, Michael J. Bearpark PII: S0022-2836(08)00987-X

More information

Structural Bioinformatics (C3210) Molecular Mechanics

Structural 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 information

Force Fields for Classical Molecular Dynamics simulations of Biomolecules. Emad Tajkhorshid

Force Fields for Classical Molecular Dynamics simulations of Biomolecules. Emad Tajkhorshid Force Fields for Classical Molecular Dynamics simulations of Biomolecules Emad Tajkhorshid Theoretical and Computational Biophysics Group, Beckman Institute Departments of Biochemistry and Pharmacology,

More information

ONETEP PB/SA: Application to G-Quadruplex DNA Stability. Danny Cole

ONETEP 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 information

Molecular Dynamics Study of the. Proton Pump Cycle of. Bacteriorhodopsin

Molecular Dynamics Study of the. Proton Pump Cycle of. Bacteriorhodopsin Molecular Dynamics Study of the Proton Pump Cycle of Bacteriorhodopsin Feng Zhou, Andreas Windemuth, Klaus Schulten Beckman Institute and Departments of Biophysics and Physics, University of Illinois at

More information

Free energy simulations

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 information

MOLECULAR DYNAMICS SIMULATIONS ON STRUCTURAL CONFORMATIONS OF RHODOPSIN AND PRION PROTEINS

MOLECULAR DYNAMICS SIMULATIONS ON STRUCTURAL CONFORMATIONS OF RHODOPSIN AND PRION PROTEINS Ó³ Ÿ. 2005.. 2, º 6(129).. 87Ä97 Š [51-72:530.145]+[51-72:541.1] MOLECULAR DYNAMICS SIMULATIONS ON STRUCTURAL CONFORMATIONS OF RHODOPSIN AND PRION PROTEINS Kh. T. Kholmurodov 1 Joint Institute for Nuclear

More information

Dynamic Structure of Retinylidene Ligand of Rhodopsin Probed by Molecular Simulations

Dynamic Structure of Retinylidene Ligand of Rhodopsin Probed by Molecular Simulations doi:10.1016/j.jmb.2007.06.047 J. Mol. Biol. (2007) 372, 906 917 Dynamic Structure of Retinylidene Ligand of Rhodopsin Probed by Molecular Simulations Pick-Wei Lau 1, Alan Grossfield 2, Scott E. Feller

More information

Potential Energy (hyper)surface

Potential 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 information

TD!DFT applied to biological problems. Marcus Elstner TU Braunschweig

TD!DFT applied to biological problems. Marcus Elstner TU Braunschweig TD!DFT applied to biological problems Marcus Elstner TU Braunschweig 0) biological structures Biological structures: proteins, DNA, lipids N lipids Amino acids: 20 DNA bases: 4 (A,C,G,T) Proteins - contain

More information

T H E J O U R N A L O F G E N E R A L P H Y S I O L O G Y. jgp

T H E J O U R N A L O F G E N E R A L P H Y S I O L O G Y. jgp S u p p l e m e n ta l m at e r i a l jgp Lee et al., http://www.jgp.org/cgi/content/full/jgp.201411219/dc1 T H E J O U R N A L O F G E N E R A L P H Y S I O L O G Y S u p p l e m e n ta l D I S C U S

More information

Coupling of Retinal, Protein, and Water Dynamics in Squid Rhodopsin

Coupling of Retinal, Protein, and Water Dynamics in Squid Rhodopsin 2200 Biophysical Journal Volume 99 October 2010 2200 2207 Coupling of Retinal, Protein, and Water Dynamics in Squid Rhodopsin Eduardo Jardón-Valadez, Ana-Nicoleta Bondar, and Douglas J. Tobias * Departments

More information

Supplementary Figure 1

Supplementary Figure 1 Supplementary Figure 1 a Tyr191 Trp265 4-5 Å Tyr268 ~14 Å 4.3 Å 16 17 19 20 4.4 Å 11 15 Lys296 N 8 10 12 14 7.3 Å 18 Tyr191 4-5 Å 18 3.8-4 Å 8 10 12 11 17 16 19 ~ Å ~4 Å Tyr268 5.5-6.3 Å 14 N Lys296 15

More information

Free energy calculations using molecular dynamics simulations. Anna Johansson

Free 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 information

Force Fields for Classical Molecular Dynamics simulations of Biomolecules. Emad Tajkhorshid

Force Fields for Classical Molecular Dynamics simulations of Biomolecules. Emad Tajkhorshid Force Fields for Classical Molecular Dynamics simulations of Biomolecules Emad Tajkhorshid Beckman Institute Departments of Biochemistry Center for Biophysics and Computational Biology University of Illinois

More information

Can 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? 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 information

Computational and spectroscopic investigation of 7-azaindole: Solvation and intermolecular interactions

Computational and spectroscopic investigation of 7-azaindole: Solvation and intermolecular interactions Computational and spectroscopic investigation of 7-azaindole: Solvation and intermolecular interactions Michael Kamrath, Krista Cruse, Nathan Erickson, Molly Beernink Abstract We report results of an experimental

More information

Electronic Supplementary Information

Electronic Supplementary Information Electronic Supplementary Material (ESI) for Chemical Science. This journal is The Royal Society of Chemistry 2016 Electronic Supplementary Information Active Site Structure and Absorption Spectrum of Channelrhodopsin-2

More information

Ab initio molecular dynamics of rhodopsin

Ab initio molecular dynamics of rhodopsin Pure&App/. Chem.,Vol. 69, No. 10, pp. 2105-2110, 1997. Printed in Great Britain. Q 1997 IUPAC Ab initio molecular dynamics of rhodopsin Angelo Bifone", Huub J.M. de Groot" and Francesco Budab Leiden Institute

More information

Why study protein dynamics?

Why 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 information

Supporting Information

Supporting Information Supporting Information Formation of Ruthenium Carbenes by gem-hydrogen Transfer to Internal Alkynes: Implications for Alkyne trans-hydrogenation Markus Leutzsch, Larry M. Wolf, Puneet Gupta, Michael Fuchs,

More information

What is the course about?

What is the course about? What is the course about? Part 1 Organic Chemistry: organocatalysis 10-10 (1 Å) - 10-9 m Part 2 Structural Biology: peptide conformations Part 3 10-9 - 10-8 m Photochemistry and Photobiology: olefins &

More information

q 2 Da This looks very similar to Coulomb doesn t it? The only difference to Coulomb is the factor of 1/2.

q 2 Da This looks very similar to Coulomb doesn t it? The only difference to Coulomb is the factor of 1/2. Born Lets now think about a somewhat different, but related problem that is associated with the name Born. Using a name (usually the name of the guy who first discovered or solved the problem) to describe

More information

Dye Sensitized Solar Cells

Dye Sensitized Solar Cells Large Scale simulations with CP2K Dye Sensitized Solar Cells Florian Schiffmann 2nd CP2K tutorial f.schiffmann@ucl.ac.uk University College London Dye sensitized solar cells Grätzel, Nature (1991,2001)

More information

Free energy calculations and the potential of mean force

Free 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 information

Proton Acidity. (b) For the following reaction, draw the arrowhead properly to indicate the position of the equilibrium: HA + K + B -

Proton Acidity. (b) For the following reaction, draw the arrowhead properly to indicate the position of the equilibrium: HA + K + B - Proton Acidity A01 Given that acid A has a pk a of 15 and acid B has a pk a of 10, then: (a) Which of the two acids is stronger? (b) For the following reaction, draw the arrowhead properly to indicate

More information

5th CCPN Matt Crump. Thermodynamic quantities derived from protein dynamics

5th 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 information

Lecture 11: Potential Energy Functions

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 information

Conformational sampling of macrocycles in solution and in the solid state

Conformational sampling of macrocycles in solution and in the solid state Conformational sampling of macrocycles in solution and in the solid state Paul Hawkins, Ph.D. Head of Scientific Solutions Stanislaw Wlodek, Ph.D. Senior Scientific Developer 6/6/2018 https://berkonomics.com/?p=2437

More information

Interpreting and evaluating biological NMR in the literature. Worksheet 1

Interpreting and evaluating biological NMR in the literature. Worksheet 1 Interpreting and evaluating biological NMR in the literature Worksheet 1 1D NMR spectra Application of RF pulses of specified lengths and frequencies can make certain nuclei detectable We can selectively

More information

Proton-Pump Mechanism in Retinal Schiff Base: On the molecular structure of the M-state

Proton-Pump Mechanism in Retinal Schiff Base: On the molecular structure of the M-state Proton-Pump Mechanism in Retinal Schiff Base: On the molecular structure of the M-state Ayan Datta and Swapan K. Pati * Theoretical Sciences Unit and Chemistry and Physics of Materials Unit, Jawaharlal

More information

- 13-cis photoisomerization: one surface hindering a Cl4-CI5 single-bond corotation and the other surface

- 13-cis photoisomerization: one surface hindering a Cl4-CI5 single-bond corotation and the other surface J. Phys. Chem. 1995,99, 14549-14560 14549 Molecular Dynamics Study of the Early Intermediates in the Bacteriorhodopsin Photocycle William Humphrey, Dong Xu, Mordechai Sheves, and Klaus Schulten*J Beckman

More information

Supporting Information

Supporting Information Supporting Information Constant ph molecular dynamics reveals ph-modulated binding of two small-molecule BACE1 inhibitors Christopher R. Ellis 1,, Cheng-Chieh Tsai 1,, Xinjun Hou 2, and Jana Shen 1, 1

More information

Ion Translocation Across Biological Membranes. Janos K. Lanyi University of California, Irvine

Ion Translocation Across Biological Membranes. Janos K. Lanyi University of California, Irvine Ion Translocation Across Biological Membranes Janos K. Lanyi University of California, Irvine Examples of transmembrane ion pumps Protein Cofactor, substrate, etc. MW Subunits Mitoch. cytochrome oxidase

More information

Electronic excitations in conjugated. Many-Body Green's Functions. Behnaz Bagheri Varnousfaderani. Behnaz Bagheri Varnousfaderani

Electronic excitations in conjugated. Many-Body Green's Functions. Behnaz Bagheri Varnousfaderani. Behnaz Bagheri Varnousfaderani Technische Universiteit Eindhoven University of Technology Electronic excitations in conjugated polymers Electronic from excitations QM/MM simulations in conjugated pol based from on Many-Body QM/MM simulations

More information

Applications of Molecular Dynamics

Applications 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 information

UV-vis (Electronic) Spectra Ch.13 Atkins, Ch.19 Engel

UV-vis (Electronic) Spectra Ch.13 Atkins, Ch.19 Engel XV 74 UV-vis (Electronic) Spectra-2014 -Ch.13 Atkins, Ch.19 Engel Most broadly used analytical tech / especially bio-applic. inexpensive optics / solvent & cell usually not problem intense transitions

More information

Quantifying Uncertainty and Sampling Quality in Biomolecular Simulations

Quantifying Uncertainty and Sampling Quality in Biomolecular Simulations CHAPTER 2 Quantifying Uncertainty and Sampling Quality in Biomolecular Simulations Alan Grossfield 1 and Daniel M. Zuckerman 2 Contents 1. Introduction 24 1.1 Examples big and small: butane and rhodopsin

More information

The Mechanistic Studies of the Wacker Oxidation. Tyler W. Wilson SED Group Meeting

The Mechanistic Studies of the Wacker Oxidation. Tyler W. Wilson SED Group Meeting The Mechanistic Studies of the Wacker xidation Tyler W. Wilson SE Group Meeting 11.27.2007 Introduction xidation of ethene by (II) chloride solutions (Phillips, 1894) -First used as a test for alkenes

More information

Supporting Infromation

Supporting Infromation Supporting Infromation Aromatic Triazole Foldamers Induced by C H X (X = F, Cl) Intramolecular Hydrogen Bonding Jie Shang,, Nolan M. Gallagher, Fusheng Bie,, Qiaolian Li,, Yanke Che, Ying Wang,*,, and

More information

CD Basis Set of Spectra that is used is that derived from comparing the spectra of globular proteins whose secondary structures are known from X-ray

CD Basis Set of Spectra that is used is that derived from comparing the spectra of globular proteins whose secondary structures are known from X-ray CD Basis Set of Spectra that is used is that derived from comparing the spectra of globular proteins whose secondary structures are known from X-ray crystallography An example of the use of CD Modeling

More information

Hyeyoung Shin a, Tod A. Pascal ab, William A. Goddard III abc*, and Hyungjun Kim a* Korea

Hyeyoung 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 information

Investigating excited state dynamics in 7 azaindole. Nathan Erickson, Molly Beernink, and Nathaniel Swenson

Investigating excited state dynamics in 7 azaindole. Nathan Erickson, Molly Beernink, and Nathaniel Swenson Investigating excited state dynamics in 7 azaindole Nathan Erickson, Molly Beernink, and Nathaniel Swenson 1 Background I 7AI Dimer Previous studies have shown that 7 azaindole (7AI) readily forms H bonded

More information

Supplementary Materials for

Supplementary Materials for advances.sciencemag.org/cgi/content/full/3/4/e1600663/dc1 Supplementary Materials for A dynamic hydrophobic core orchestrates allostery in protein kinases Jonggul Kim, Lalima G. Ahuja, Fa-An Chao, Youlin

More information

Computational and Spectroscopic Investigation of Solution Phase Excited State Dynamics in 7 azaindole

Computational and Spectroscopic Investigation of Solution Phase Excited State Dynamics in 7 azaindole Computational and Spectroscopic Investigation of Solution Phase Excited State Dynamics in 7 azaindole Nathan Erickson, Molly Beernink, and Nathaniel Swenson Midwest Undergraduate Computational Chemistry

More information

Replica exchange methodology. John Karanicolas June 2003

Replica exchange methodology. John Karanicolas June 2003 Replica exchange methodology John Karanicolas MMTSB@PSC, June 2003 Outline o Motivation o Theory o Practical considerations o MMTSB Tool Set Why? o Sampling on rugged potential energy surfaces is difficult

More information

The Dominant Interaction Between Peptide and Urea is Electrostatic in Nature: A Molecular Dynamics Simulation Study

The Dominant Interaction Between Peptide and Urea is Electrostatic in Nature: A Molecular Dynamics Simulation Study Dror Tobi 1 Ron Elber 1,2 Devarajan Thirumalai 3 1 Department of Biological Chemistry, The Hebrew University, Jerusalem 91904, Israel 2 Department of Computer Science, Cornell University, Ithaca, NY 14853

More information

Molecular 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 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 information

All-atom Molecular Mechanics. Trent E. Balius AMS 535 / CHE /27/2010

All-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 information

Multi-Ensemble Markov Models and TRAM. Fabian Paul 21-Feb-2018

Multi-Ensemble Markov Models and TRAM. Fabian Paul 21-Feb-2018 Multi-Ensemble Markov Models and TRAM Fabian Paul 21-Feb-2018 Outline Free energies Simulation types Boltzmann reweighting Umbrella sampling multi-temperature simulation accelerated MD Analysis methods

More information

Chemistry 4021/8021 Computational Chemistry 3/4 Credits Spring Semester 2011 ( Due 2 / 21 / 11 )

Chemistry 4021/8021 Computational Chemistry 3/4 Credits Spring Semester 2011 ( Due 2 / 21 / 11 ) Chemistry 4021/8021 Computational Chemistry 3/4 Credits Spring Semester 2011 ( Due 2 / 21 / 11 ) Using PC Model, answer the questions below. 1. Consider the water molecule, 2. Using the data for the force

More information

RHODOPSIN: Insights from Recent Structural Studies

RHODOPSIN: Insights from Recent Structural Studies Annu. Rev. Biophys. Biomol. Struct. 2002. 31:443 84 DOI: 10.1146/annurev.biophys.31.082901.134348 Copyright c 2002 by Annual Reviews. All rights reserved RHODOPSIN: Insights from Recent Structural Studies

More information

4 Examples of enzymes

4 Examples of enzymes Catalysis 1 4 Examples of enzymes Adding water to a substrate: Serine proteases. Carbonic anhydrase. Restrictions Endonuclease. Transfer of a Phosphoryl group from ATP to a nucleotide. Nucleoside monophosphate

More information

CHAPTER-IV. FT-IR and FT-Raman investigation on m-xylol using ab-initio HF and DFT calculations

CHAPTER-IV. FT-IR and FT-Raman investigation on m-xylol using ab-initio HF and DFT calculations 4.1. Introduction CHAPTER-IV FT-IR and FT-Raman investigation on m-xylol using ab-initio HF and DFT calculations m-xylol is a material for thermally stable aramid fibers or alkyd resins [1]. In recent

More information

Chap. 8 Substitution Reactions

Chap. 8 Substitution Reactions Chap. 8 Substitution Reactions Y + R X R' Y + X Nucleophilic not necessarily the same as R Electrophilic S N 1 slow (C 3 ) 3 CCl (C + Cl - 3 ) 3 C + (C 3 ) 3 C + OC 2 C 3 C 3 C 2 O C 3 C 2 O d[( C ) 3CCl]

More information

Solvent & geometric effects on non-covalent interactions

Solvent & 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 information

Molecular dynamics simulations of anti-aggregation effect of ibuprofen. Wenling E. Chang, Takako Takeda, E. Prabhu Raman, and Dmitri Klimov

Molecular 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 information

LS1a Midterm Exam 1 Review Session Problems

LS1a Midterm Exam 1 Review Session Problems LS1a Midterm Exam 1 Review Session Problems 1. n aqueous mixture of a weak acid and its conjugate base is often used in the laboratory to prepare solutions referred to as buffers. ne commonly used acid

More information

sensors ISSN

sensors ISSN Sensors 2009, 9, 9741-9804; doi:10.3390/s91209741 OPEN ACCESS sensors ISSN 1424-8220 www.mdpi.com/journal/sensors Review Photoreactions and Structural Changes of Anabaena Sensory Rhodopsin Akira Kawanabe

More information

Basic Organic Chemistry Course code : CHEM (Pre-requisites : CHEM 11122)

Basic Organic Chemistry Course code : CHEM (Pre-requisites : CHEM 11122) Basic Organic Chemistry Course code : CHEM 12162 (Pre-requisites : CHEM 11122) Chapter 01 Mechanistic Aspects of S N2,S N1, E 2 & E 1 Reactions Dr. Dinesh R. Pandithavidana Office: B1 222/3 Phone: (+94)777-745-720

More information

Molecular dynamics simulation of Aquaporin-1. 4 nm

Molecular 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 information

Molecular Mechanics, Dynamics & Docking

Molecular 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 information

Calculation of Proton Transfers in Bacteriorhodopsin br and M Intermediates

Calculation of Proton Transfers in Bacteriorhodopsin br and M Intermediates Biochemistry 2003, 42, 9875-9888 9875 Calculation of Proton Transfers in Bacteriorhodopsin br and M Intermediates Yifan Song, Junjun Mao, and M. R. Gunner* Physics Department J-419, City College of New

More information

Characterizing Structural Transitions of Membrane Transport Proteins at Atomic Detail Mahmoud Moradi

Characterizing Structural Transitions of Membrane Transport Proteins at Atomic Detail Mahmoud Moradi Characterizing Structural Transitions of Membrane Transport Proteins at Atomic Detail Mahmoud Moradi NCSA Blue Waters Symposium for Petascale Science and Beyond Sunriver, Oregon May 11, 2015 Outline Introduction

More information

Biomolecular modeling. Theoretical Chemistry, TU Braunschweig (Dated: December 10, 2010)

Biomolecular 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 information

Lectures Spectroscopy. Fall 2012

Lectures Spectroscopy. Fall 2012 Lectures 19-20 Spectroscopy Fall 2012 1 spectroscopic principles (Chem 1M/1N exps. #6 and #11) 4 spectroscopic excitations ( E = h = hc/ ; = c ) (nm) (sec -1 ) radiation technique molecular excitation

More information

Protein Simulations in Confined Environments

Protein 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 information

Solutions and Non-Covalent Binding Forces

Solutions and Non-Covalent Binding Forces Chapter 3 Solutions and Non-Covalent Binding Forces 3.1 Solvent and solution properties Molecules stick together using the following forces: dipole-dipole, dipole-induced dipole, hydrogen bond, van der

More information

Improving Protein Function Prediction with Molecular Dynamics Simulations. Dariya Glazer Russ Altman

Improving Protein Function Prediction with Molecular Dynamics Simulations. Dariya Glazer Russ Altman Improving Protein Function Prediction with Molecular Dynamics Simulations Dariya Glazer Russ Altman Motivation Sometimes the 3D structure doesn t score well for a known function. The experimental structure

More information

The Signaling Pathway of Rhodopsin

The Signaling Pathway of Rhodopsin Article The Yifei Kong 1 and Martin Karplus 1,2, * 1 Department of Chemistry and Chemical Biology, Harvard University, Cambridge, MA 02138, USA 2 Laboratoire de Chimie Biophysique, ISIS, Université Louis

More information

The Molecular Dynamics Method

The Molecular Dynamics Method H-bond energy (kcal/mol) - 4.0 The Molecular Dynamics Method Fibronectin III_1, a mechanical protein that glues cells together in wound healing and in preventing tumor metastasis 0 ATPase, a molecular

More information

Molecular Mechanics. I. Quantum mechanical treatment of molecular systems

Molecular 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 information

Computing free energy: Thermodynamic perturbation and beyond

Computing free energy: Thermodynamic perturbation and beyond Computing free energy: Thermodynamic perturbation and beyond Extending the scale Length (m) 1 10 3 Potential Energy Surface: {Ri} 10 6 (3N+1) dimensional 10 9 E Thermodynamics: p, T, V, N continuum ls

More information

Molecular 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 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 information

Solid-State Deuterium NMR Spectroscopy of Rhodopsin

Solid-State Deuterium NMR Spectroscopy of Rhodopsin Solid-State Deuterium NMR Spectroscopy of Rhodopsin Suchithranga M.D.C. Perera, Xiaolin Xu, Trivikram R. Molugu, Andrey V. Struts, and Michael F. Brown Abstract Rhodopsin is a prototype for the large Family

More information

The Molecular Dynamics Simulation Process

The Molecular Dynamics Simulation Process The Molecular Dynamics Simulation Process For textbooks see: M.P. Allen and D.J. Tildesley. Computer Simulation of Liquids.Oxford University Press, New York, 1987. D. Frenkel and B. Smit. Understanding

More information

Model-Free Approach to Internal Motions in Proteins

Model-Free Approach to Internal Motions in Proteins Model-Free Approach to Internal Motions in Proteins Lipari & Szabo, JACS 104, 4546 (1982) Palmer AG. Ann. Rev. Biophys. Biomol. Struc., 30, 129-155 (2001) Palmer AG, Kroenke CD, Loria JP, Meth. Enzymol.

More information

Some Accelerated/Enhanced Sampling Techniques

Some Accelerated/Enhanced Sampling Techniques Some Accelerated/Enhanced Sampling Techniques The Sampling Problem Accelerated sampling can be achieved by addressing one or more of the MD ingredients Back to free energy differences... Examples of accelerated

More information

Supporting Information

Supporting Information Supporting Information ph-responsive self-assembly of polysaccharide through a rugged energy landscape Brian H. Morrow, Gregory F. Payne, and Jana Shen Department of Pharmaceutical Sciences, School of

More information

Hamiltonian Replica Exchange Molecular Dynamics Using Soft-Core Interactions to Enhance Conformational Sampling

Hamiltonian Replica Exchange Molecular Dynamics Using Soft-Core Interactions to Enhance Conformational Sampling John von Neumann Institute for Computing Hamiltonian Replica Exchange Molecular Dynamics Using Soft-Core Interactions to Enhance Conformational Sampling J. Hritz, Ch. Oostenbrink published in From Computational

More information

Fondamenti di Chimica Farmaceutica. Computer Chemistry in Drug Research: Introduction

Fondamenti 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 information

Toward An Understanding Of The Retinal Chromophore In Rhodopsin Mimics

Toward An Understanding Of The Retinal Chromophore In Rhodopsin Mimics Bowling Green State University ScholarWorks@BGSU Chemistry Faculty Publications Chemistry 9-2013 Toward An Understanding Of The Retinal Chromophore In Rhodopsin Mimics Mark M. Huntress Samer Gozem Konstantin

More information

Effects of Chemical Exchange on NMR Spectra

Effects of Chemical Exchange on NMR Spectra Effects of Chemical Exchange on NMR Spectra Chemical exchange refers to any process in which a nucleus exchanges between two or more environments in which its NMR parameters (e.g. chemical shift, scalar

More information

The intricate dance of electrons, photons, and ions when matter responds to light. Yusheng Dou and Roland E. Allen Texas A&M University

The intricate dance of electrons, photons, and ions when matter responds to light. Yusheng Dou and Roland E. Allen Texas A&M University The intricate dance of electrons, photons, and ions when matter responds to light Yusheng Dou and Roland E. Allen Texas A&M University This talk is based on papers at our web site, http://lightandmatter.net.

More information

Molecular Modelling for Medicinal Chemistry (F13MMM) Room A36

Molecular Modelling for Medicinal Chemistry (F13MMM) Room A36 Molecular Modelling for Medicinal Chemistry (F13MMM) jonathan.hirst@nottingham.ac.uk Room A36 http://comp.chem.nottingham.ac.uk Assisted reading Molecular Modelling: Principles and Applications. Andrew

More information

The Molecular Dynamics Method

The Molecular Dynamics Method H-bond energy (kcal/mol) - 4.0 The Molecular Dynamics Method Fibronectin III_1, a mechanical protein that glues cells together in wound healing and in preventing tumor metastasis 0 ATPase, a molecular

More information

Experimental Methods in Kinetics (2014)

Experimental Methods in Kinetics (2014) IIIa- 35 Experimental Methods in Kinetics (2014) Initial examples assume start R = A 0 + B 0 etc., P = 0 follow reaction forward initial rate, etc. measure [A] vs. t Chemical methods: Could mix and take

More information

Homework Problem Set 4 Solutions

Homework Problem Set 4 Solutions Chemistry 380.37 Dr. Jean M. Standard omework Problem Set 4 Solutions 1. A conformation search is carried out on a system and four low energy stable conformers are obtained. Using the MMFF force field,

More information

Protein Folding & Stability. Lecture 11: Margaret A. Daugherty. Fall How do we go from an unfolded polypeptide chain to a

Protein 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 information

(1) Solvent interactions Water is a ligand itself so it can become involved in the chemistry: Overall reaction: [L 5 MX] + Y [L 5 MY] + X

(1) Solvent interactions Water is a ligand itself so it can become involved in the chemistry: Overall reaction: [L 5 MX] + Y [L 5 MY] + X (1) Solvent interactions Water is a ligand itself so it can become involved in the chemistry: Overall reaction: [L 5 MX] + Y [L 5 MY] + X 270 1.(slow) [L 5 MX] + H 2 O [L 5 M(H 2 O)] + X 2. (fast) [L 5

More information

Molecular Mechanics. Yohann Moreau. November 26, 2015

Molecular 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 information