Kirkwood-Buff Analysis of Aqueous N-Methylacetamide and Acetamide Solutions. Modeled by the CHARMM Additive and Drude Polarizable Force Fields

Size: px
Start display at page:

Download "Kirkwood-Buff Analysis of Aqueous N-Methylacetamide and Acetamide Solutions. Modeled by the CHARMM Additive and Drude Polarizable Force Fields"

Transcription

1 Supporting Information for Kirkwood-Buff Analysis of Aqueous N-Methylacetamide and Acetamide Solutions Modeled by the CHARMM Additive and Drude Polarizable Force Fields Bin Lin, 1 Pedro E. M. Lopes, 1 Benoît Roux, 2 Alexander D. MacKerell, Jr. 1 1 Department of Pharmaceutical Sciences, University of Maryland, School of Pharmacy, 20 Penn Street HSFII, Baltimore, Maryland 21201, USA 2 Department of Biochemistry and Molecular Biology, University of Chicago, Chicago, Illinois 60637, USA

2 Figure S1. Dihedral scan for the O6=C5- N7- C9 torsion of NMA. The energy at 0 is set to zero.

3 N cc Excess Coordination Number acetamide (ACEM) additive drude exp N cw N ww x c Figure S2. Excess coordination numbers as a function of mole fraction for the ACEM solutions.

4 Table S1. Vibrational data for ACEM a. Ab initio (MP2/6-31G*) Drude Mode frequency assignment % cont. frequency assignment % cont tcch tcch wnh bccn 100 tcn bccn wnh wco tcn tcn 96 wnh rco rco ncn wco wco wnh ncc ncc rch rnh 53 ncn rch 3 ' wco dch 3 as rch rnh dch 3 as rch 3 ' dch 3 s ncn bhnh dch 3 s rco' dch 3 s dch 3 as 36 ' rch dch 3 as dch 3 as rch 3 ' rch dch 3 as bhnh 76 ' rch bhnh dch 3 s ncc nco nco nch 3 s nch 3 s nch 3 as' nch 3 as nch 3 as' nch 3 as nch 3 as' nch 3 as nnht nnhc nnht nnhc nnhc nnht nnhc nnht rnh

5 a Frequencies in cm -1. Normal mode contributions determined with the MOLVIB module in CHARMM. Only modes contributing more than 20% are shown. n indicates a stretching mode, b indicates a bending mode, d indicates a deformation mode, t is used for torsional rotations and r is used for rocking modes. Ab initio data was obtained at the MP2/6-31G* level of theory and frequencies were scaled by (Scott and Radom, J. Phys. Chem. 1996, 100, ).

6 Table S2. Vibrational data for NMAa. Mode Ab initio (MP2/6-31G*) frequency assignment % cont tnch tcch tcn bcnc 56 bccn wnh Drude frequency assignment 51.7 tcch tnch tcn bccn rco bccn wco wco vcn rch3' vcc rch3 vnch3 rch3' rch3 rnh vcn dch3s dch3s dch3as' dch3as dch3as' dch3as rnh vco vch3s vch3s vch3as' vch3as vch3as' vch3as vnh % cont wnh wco bccn ncn rco wco wnh vcn vcc vcn rch3 rch3' rch3 rch3' rco dch3as dch3as' rch3' dch3as' rch3 dch3as rch3 dch3s dch3s vcc vco rnh vch3s vch3s vch3as vch3as vch3as' vch3as' vnh

7 a Frequencies in cm -1. Normal mode contributions determined with the MOLVIB module in CHARMM. Notations are the same as in Table S1. Ab initio data was obtained at the MP2/6-31G* level of theory and frequencies were scaled by (Scott and Radom, J. Phys. Chem. 1996, 100, ).

8 Table S3. Topology entries for NMA in CHARMM format. RESI NMA 0.000! N-methylacetamide!! Hl1 H! \ /! Hl2--CL--C---N! / \! Hl3 O Cc--Hc3 (typically cis to O)! / \! Hc1 Hc2! GROUP ATOM CL CD33C ALPHA THOLE ATOM HL1 HDA3A ATOM HL2 HDA3A ATOM HL3 HDA3A ATOM C CD2O1A ALPHA THOLE ATOM O OD2C1A ALPHA THOLE ATOM N ND2A ALPHA THOLE ATOM H HDP1A ATOM CC CD33C ALPHA THOLE ATOM HC1 HDA3A ATOM HC2 HDA3A ATOM HC3 HDA3A ATOM LPA LP ATOM LPB LP BOND HL1 CL HL2 CL HL3 CL BOND CL C C N N CC BOND C O N H BOND HC1 CC HC2 CC HC3 CC BOND O LPA O LPB IMPR C CL N O N C CC H LONEPAIR relative LPA O C CL distance 0.30 angle 91.0 dihe 0.0 LONEPAIR relative LPB O C CL distance 0.30 angle 91.0 dihe ANISOTROPY O C LPA LPB A A

9 Table S4. Topology entries for ACEM in CHARMM format. RESI ACEM 0.000! acetamide!! the amide charges listed below are used in asparagine and glutamine!! HL1 Ht (trans to O)! \ /! HL2--CL--C---N! / \! HL3 O Hc (cis to O)! GROUP ATOM CL CD33C ALPHA THOLE ATOM HL1 HDA3A ATOM HL2 HDA3A ATOM HL3 HDA3A ATOM C CD2O1A ALPHA THOLE ATOM O OD2C1A ALPHA THOLE ATOM N ND2A ALPHA THOLE ATOM HC HDP1A ATOM HT HDP1A ATOM LPA LP ATOM LPB LP BOND C N N Hc N Ht C O BOND C CL CL HL1 CL HL2 CL HL3 BOND O LPA O LPB IMPR C CL N O IMPR N HC HT C!from NMA LONEPAIR relative LPA O C CL distance 0.30 angle 91.0 dihe 0.0 LONEPAIR relative LPB O C CL distance 0.30 angle 91.0 dihe ANISOTROPY O C LPA LPB A A

10 Table S5. Internal parameters for NMA and ACEM in CHARMM format. BONDS K b : kcal/mol/å 2 B 0 : Å K b B 0 ND2A2 CD2O1A ND2A2 CD33C ND2A2 HDP1A OD2C1A CD2O1A CD2O1A CD33C CD33C HDA3A ND2A1 CD2O1A ND2A1 HDP1A ANGLES K theta : kcal/mol/rad 2 Theta 0 : degree K UB : kcal/mol/å 2 S 0 : Å K theta Theta 0 K UB S 0 CD2O1A ND2A2 CD33C CD2O1A ND2A2 HDP1A CD33C ND2A2 HDP1A ND2A2 CD2O1A OD2C1A ND2A2 CD2O1A CD33C OD2C1A CD2O1A CD33C ND2A2 CD33C HDA3A CD2O1A CD33C HDA3A CD2O1A ND2A1 HDP1A HDP1A ND2A1 HDP1A ND2A1 CD2O1A OD2C1A ND2A1 CD2O1A CD33C

11 DIHEDRALS K chi : kcal/mol n: multiplicity delta: degrees K chi n delta CD33C ND2A2 CD2O1A CD33C CD33C ND2A2 CD2O1A CD33C HDP1A ND2A2 CD2O1A CD33C CD2O1A ND2A2 CD33C HDA3A HDP1A ND2A2 CD33C HDA3A ND2A2 CD2O1A CD33C HDA3A OD2C1A CD2O1A CD33C HDA3A HDP1A ND2A1 CD2O1A OD2C1A HDP1A ND2A1 CD2O1A CD33C ND2A1 CD2O1A CD33C HDA3A IMPROPERS K psi : kcal/mol/å 2 PSI 0 : Å K psi PSI 0 CD2O1A CD33C ND2A2 OD2C1A ND2A2 CD2O1A CD33C HDP1A OD2C1A ND2A1 CD33C CD2O1A ND2A1 HDP1A HDP1A CD2O1A

12 Table S6. Lennard-Jones parameters. Epsilon in kcal/mol and Rmin in Å. Lennard-Jones parameters epsilon Rmin/2 CD2O1A CD33C ND2A OD2C1A CD2O1A HDA3A ND2A Lennard-Jones off diagonal terms to be used in combination with the SWM4-NDP water model oxygen LJ parameters epsilon Rmin ND2A2 ODW OD2C1A ODW

CHEM 498Q / CHEM 630Q: Molecular Modeling of Proteins TUTORIAL #3a: Empirical force fields

CHEM 498Q / CHEM 630Q: Molecular Modeling of Proteins TUTORIAL #3a: Empirical force fields Department of Chemistry and Biochemistry, Concordia University! page 1 of 6 CHEM 498Q / CHEM 630Q: Molecular Modeling of Proteins TUTORIAL #3a: Empirical force fields INTRODUCTION The goal of this tutorial

More information

CHEM 498Q / CHEM 630Q: Molecular Modeling of Proteins TUTORIAL #3a: Empirical force fields

CHEM 498Q / CHEM 630Q: Molecular Modeling of Proteins TUTORIAL #3a: Empirical force fields Department of Chemistry and Biochemistry, Concordia University page 1 of 6 CHEM 498Q / CHEM 630Q: Molecular Modeling of Proteins TUTORIAL #3a: Empirical force fields INTRODUCTION The goal of this tutorial

More information

Figure 1. Molecules geometries of 5021 and Each neutral group in CHARMM topology was grouped in dash circle.

Figure 1. Molecules geometries of 5021 and Each neutral group in CHARMM topology was grouped in dash circle. Project I Chemistry 8021, Spring 2005/2/23 This document was turned in by a student as a homework paper. 1. Methods First, the cartesian coordinates of 5021 and 8021 molecules (Fig. 1) are generated, in

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

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

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

Classical Force Fields

Classical Force Fields Classical Force Fields Coarse-grained, Full atom, Hybrid Knowledge-based and Physics-based Problems: structure prediction, protein folding kinetics membrane, protein insertion kinetics mechanism, protein

More information

Overview of Classical Force Fields and Parameter Optimization Strategy

Overview of Classical Force Fields and Parameter Optimization Strategy Overview of Classical Force Fields and Parameter Optimization Strategy Part I - Overview of CHARMM FF and Parameter Optimization Part II - Introduction to Quantum Chemistry Calculations (SPARTAN) Application

More information

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

Supplemental Material: Solute Adsorption at Air-Water Interfaces and Induced. Interface Fluctuations: The Hydrophobic Nature of Ions?

Supplemental 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

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

Overview of Classical Force Fields and Parameter Optimization Strategy

Overview of Classical Force Fields and Parameter Optimization Strategy Overview of Classical Force Fields and Parameter Optimization Strategy Part I - Overview of CHARMM FF and Parameter Optimization Part II - Introduction to Quantum Chemistry Calculations (SPARTAN) Application

More information

Understanding the determinants of selectivity in drug metabolism through modeling of dextromethorphan oxidation by cytochrome P450

Understanding the determinants of selectivity in drug metabolism through modeling of dextromethorphan oxidation by cytochrome P450 Understanding the determinants of selectivity in drug metabolism through modeling of dextromethorphan oxidation by cytochrome P450 Julianna Oláh, 1 Adrian J. Mulholland* and Jeremy N. Harvey* School of

More information

SUPPLEMENTARY INFORMATION

SUPPLEMENTARY INFORMATION DOI: 10.1038/NCHEM.1299 Protein fold determined by paramagnetic magic-angle spinning solid-state NMR spectroscopy Ishita Sengupta 1, Philippe S. Nadaud 1, Jonathan J. Helmus 1, Charles D. Schwieters 2

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

Nanotube AFM Probe Resolution

Nanotube AFM Probe Resolution Influence of Elastic Deformation on Single-Wall Carbon Nanotube AFM Probe Resolution Ian R. Shapiro, Santiago D. Solares, Maria J. Esplandiu, Lawrence A. Wade, William A. Goddard,* and C. Patrick Collier*

More information

Kirkwood-Buff Integrals for Aqueous Urea Solutions Based upon the Quantum Chemical Electrostatic Potential and Interaction Energies

Kirkwood-Buff Integrals for Aqueous Urea Solutions Based upon the Quantum Chemical Electrostatic Potential and Interaction Energies Supporting Information for Kirkwood-Buff Integrals for Aqueous Urea Solutions Based upon the Quantum Chemical Electrostatic Potential and Interaction Energies Shuntaro Chiba, 1* Tadaomi Furuta, 2 and Seishi

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

Pc2h = 22A~-}-lTB~+18Au + 21B~.

Pc2h = 22A~-}-lTB~+18Au + 21B~. INTERRETATION OF THE VIBRATIONAL SECTRA OF TERT-BUTYL AND DIMETHYLETHYNILMETHYL EROXIDES I.. Zyat'kov, G, A. itsevich, A.. Yuvchenko, Yu. L. Ol'dekop, V. I. Gogolinskii, V. L. Antonovskii, and D. I. Sagaidak

More information

Semi Empirical Force Fields and Their Limitations. Potential Energy Surface (PES)

Semi 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 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

TOPOLOGIES AND FORCE FIELD PARAMETERS FOR NITROXIDE SPIN LABELS

TOPOLOGIES AND FORCE FIELD PARAMETERS FOR NITROXIDE SPIN LABELS TOPOLOGIES AND FORCE FIELD PARAMETERS FOR NITROXIDE SPIN LABELS 1. Introduction The force field used in these studies is the CHARMM19 (Chemistry at HARvard Molecular Mechanics) extended atom force field

More information

Practical 1: Structure and electronic properties of organic molecules. B/ Structure, electronic and vibrational properties of the water molecule

Practical 1: Structure and electronic properties of organic molecules. B/ Structure, electronic and vibrational properties of the water molecule D1CH9116 - MMCO Molecular Modelling applied to organic chemistry Practical 1: tructure and electronic properties of organic molecules B/ tructure, electronic and vibrational properties of the water molecule

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

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

Pancake Bonding in 1,2,3,5-dithiadiazolyl and 1,2,3,5-diselenadiazolyl Radical Dimers and their Derivatives

Pancake Bonding in 1,2,3,5-dithiadiazolyl and 1,2,3,5-diselenadiazolyl Radical Dimers and their Derivatives This journal is The Owner Societies 12 Bonds or Not Bonds? Electronic Supplementary Information Pancake Bonding in 1,2,3,-dithiadiazolyl and 1,2,3,-diselenadiazolyl Radical Dimers and their Derivatives

More information

Advanced Electronic Structure Theory Density functional theory. Dr Fred Manby

Advanced Electronic Structure Theory Density functional theory. Dr Fred Manby Advanced Electronic Structure Theory Density functional theory Dr Fred Manby fred.manby@bris.ac.uk http://www.chm.bris.ac.uk/pt/manby/ 6 Strengths of DFT DFT is one of many theories used by (computational)

More information

Polarizable Empirical Force Field for Hexopyranose Monosaccharides Based on the Classical Drude Oscillator

Polarizable Empirical Force Field for Hexopyranose Monosaccharides Based on the Classical Drude Oscillator This is an open access article published under an ACS AuthorChoice License, which permits copying and redistribution of the article or any adaptations for non-commercial purposes. pubs.acs.org/jpcb Polarizable

More information

Supporting Information

Supporting Information Supporting Information Remote Stereoinductive Intramolecular Nitrile Oxide Cycloaddition: Asymmetric Total Synthesis and Structure Revision of ( )-11 -Hydroxycurvularin Hyeonjeong Choe, Thuy Trang Pham,

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

Example questions for Molecular modelling (Level 4) Dr. Adrian Mulholland

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

Supplementary Information

Supplementary Information Supplementary Information Supplementary Figure 1. Photographs show the titration experiments by dropwise adding ~5 times number of moles of (a) LiOH and LiOH+H 2 O, (b) H 2 O 2 and H 2 O 2 +LiOH, (c) Li

More information

Supplementary Information

Supplementary Information Supplementary Information Ballistic Thermal Transport in Carbyne and Cumulene with Micron-Scale Spectral Acoustic Phonon Mean Free Path Mingchao Wang and Shangchao Lin * Department of Mechanical Engineering,

More information

The nonplanarity of the peptide group: Molecular dynamics simulations

The nonplanarity of the peptide group: Molecular dynamics simulations University of New Orleans ScholarWorks@UNO Chemistry Faculty Publications Department of Chemistry 3-1-2000 The nonplanarity of the peptide group: Molecular dynamics simulations Steven W. Rick University

More information

Why the Sulfinyl Group is special in DMSO? Chao Lv June 4, 2014

Why the Sulfinyl Group is special in DMSO? Chao Lv June 4, 2014 Why the Sulfinyl Group is special in DMSO? Chao Lv June 4, 2014 The Parameterization of Dimethyl Sulfoxide (DMSO) Nucleic Acids are known to be difficult to be parameterized because: 1. The interac

More information

Available online at

Available online at Available online at www.sciencedirect.com Chemical Physics Letters 452 (2008) 78 83 www.elsevier.com/locate/cplett Geometry and excitation energy fluctuations of NMA in aqueous solution with CHARMM, AMBER,

More information

Rotamers in the CHARMM19 Force Field

Rotamers in the CHARMM19 Force Field Appendix A Rotamers in the CHARMM19 Force Field The people may be made to follow a path of action, but they may not be made to understand it. Confucius (551 BC - 479 BC) ( ) V r 1 (j),r 2 (j),r 3 (j),...,r

More information

SUPPLEMENTAL MATERIAL

SUPPLEMENTAL MATERIAL SUPPLEMENTAL MATERIAL Systematic Coarse-Grained Modeling of Complexation between Small Interfering RNA and Polycations Zonghui Wei 1 and Erik Luijten 1,2,3,4,a) 1 Graduate Program in Applied Physics, Northwestern

More information

Supporting Information. Substitutent Rate Effects

Supporting Information. Substitutent Rate Effects Supporting Information Gosteli Claisen Rearrangement: DFT Study of Substitutent Rate Effects Julia Rehbein* and Martin Hiersemann* Fakultät Chemie, Technische Universität Dortmund, 44227 Dortmund, Germany

More information

Types of Molecular Vibrations

Types of Molecular Vibrations Important concepts in IR spectroscopy Vibrations that result in change of dipole moment give rise to IR absorptions. The oscillating electric field of the radiation couples with the molecular vibration

More information

Development and Application of Combined Quantum Mechanical and Molecular Mechanical Methods

Development and Application of Combined Quantum Mechanical and Molecular Mechanical Methods University of Nebraska - Lincoln DigitalCommons@University of Nebraska - Lincoln Student Research Projects, Dissertations, and Theses - Chemistry Department Chemistry, Department of Fall 12-2014 Development

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

3rd Advanced in silico Drug Design KFC/ADD Molecular mechanics intro Karel Berka, Ph.D. Martin Lepšík, Ph.D. Pavel Polishchuk, Ph.D.

3rd Advanced in silico Drug Design KFC/ADD Molecular mechanics intro Karel Berka, Ph.D. Martin Lepšík, Ph.D. Pavel Polishchuk, Ph.D. 3rd Advanced in silico Drug Design KFC/ADD Molecular mechanics intro Karel Berka, Ph.D. Martin Lepšík, Ph.D. Pavel Polishchuk, Ph.D. Thierry Langer, Ph.D. Jana Vrbková, Ph.D. UP Olomouc, 23.1.-26.1. 2018

More information

Parameter Choice Matters: Validating Probe Parameters for Use in Mixed-Solvent Simulations

Parameter Choice Matters: Validating Probe Parameters for Use in Mixed-Solvent Simulations Parameter Choice Matters: Validating Probe Parameters for Use in Mixed-Solvent Simulations Katrina W. Lexa, Garrett B. Goh, Heather A. Carlson, Department of Medicinal Chemistry, College of Pharmacy, University

More information

Reactive Empirical Force Fields

Reactive Empirical Force Fields Reactive Empirical Force Fields Jason Quenneville jasonq@lanl.gov X-1: Solid Mechanics, EOS and Materials Properties Applied Physics Division Los Alamos National Laboratory Timothy C. Germann, Los Alamos

More information

NIH Public Access Author Manuscript J Comput Chem. Author manuscript; available in PMC 2010 September 1.

NIH Public Access Author Manuscript J Comput Chem. Author manuscript; available in PMC 2010 September 1. NIH Public Access Author Manuscript Published in final edited form as: J Comput Chem. 2009 September ; 30(12): 1821 1838. doi:10.1002/jcc.21183. Polarizable Empirical Force Field for Nitrogen-containing

More information

Biochemistry,530:,, Introduc5on,to,Structural,Biology, Autumn,Quarter,2015,

Biochemistry,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 information

Supporting Information

Supporting Information Supporting Information Copyright Wiley-VCH Verlag GmbH & Co. KGaA, 69451 Weinheim, 2013 Photochemistry of N-Methylformamide: Matrix Isolation and Nonadiabatic Dynamics Rachel Crespo-Otero, [a] Artur Mardyukov,

More information

Infrared Spectroscopy. Provides information about the vibraions of functional groups in a molecule

Infrared Spectroscopy. Provides information about the vibraions of functional groups in a molecule Infrared Spectroscopy Provides information about the vibraions of functional groups in a molecule Therefore, the functional groups present in a molecule can be deduced from an IR spectrum Two important

More information

CHEMISTRY 4021/8021 MIDTERM EXAM 1 SPRING 2014

CHEMISTRY 4021/8021 MIDTERM EXAM 1 SPRING 2014 CHEMISTRY 4021/8021 Q1) Propose a simple, united-atom molecular mechanics force-field needed to generate a potential energy surface for an isolated molecule of acetone (Me 2 CO). I.e., provide an energy

More information

Revised CHARMM Force Field Parameters for Iron- Containing Cofactors of Photosystem II

Revised CHARMM Force Field Parameters for Iron- Containing Cofactors of Photosystem II WWW.C-CHEM.ORG FULL PAPER Revised CHARMM Force Field Parameters for Iron- Containing Cofactors of Photosystem II Suliman Adam, [a] Michaela Knapp-Mohammady, [b] Jun Yi, [c] and Ana-Nicoleta Bondar [a]

More information

k θ (θ θ 0 ) 2 angles r i j r i j

k θ (θ θ 0 ) 2 angles r i j r i j 1 Force fields 1.1 Introduction The term force field is slightly misleading, since it refers to the parameters of the potential used to calculate the forces (via gradient) in molecular dynamics simulations.

More 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

Specific Ion Solvtion in Ethylene Carbonate and Propylene Carbonate

Specific Ion Solvtion in Ethylene Carbonate and Propylene Carbonate Specific Ion Solvtion in Ethylene Carbonate and Propylene Carbonate A. Arslanargin, A. Powers, S. Rick, T. Pollard, T. Beck Univ Cincinnati Chemistry Support: NSF, OSC TSRC 2016 November 2, 2016 A. Arslanargin,

More information

Theoretical study of unusual Bis(amino) (2,4,6-tri-t-butylphenyl)borane B(NH 2 ) 2 NHAr)

Theoretical study of unusual Bis(amino) (2,4,6-tri-t-butylphenyl)borane B(NH 2 ) 2 NHAr) Theoretical study of unusual Bis(amino) (2,4,6-tri-t-butylphenyl)borane B(NH 2 ) 2 NHAr) Hatam A. Jasim 1* Hanan M. Ali 2 1* Depart. Of Pharmaceutical Chemistry, College of Pharmacy, University of Basra

More information

read rtf card append * Initial topology guesses generated by * CHARMM General Force Field (CGenFF) program version 0.9.

read rtf card append * Initial topology guesses generated by * CHARMM General Force Field (CGenFF) program version 0.9. Electronic Supplementary Material (ESI) for RSC Advances. This journal is The Royal Society of Chemistry 2014 * Stream file for dehydro-amino residues * If using these parameters please cite: * Turpin,

More information

Electronic Supplementary Information for A molecular dynamics investigation on the cross-linking and physical properties of epoxy-based materials

Electronic Supplementary Information for A molecular dynamics investigation on the cross-linking and physical properties of epoxy-based materials Electronic Supplementary Material (ESI) for RSC Advances. This journal is The Royal Society of Chemistry 2014 Electronic Supplementary Information for A molecular dynamics investigation on the cross-linking

More information

HOMOLOGY MODELING. The sequence alignment and template structure are then used to produce a structural model of the target.

HOMOLOGY MODELING. The sequence alignment and template structure are then used to produce a structural model of the target. HOMOLOGY MODELING Homology modeling, also known as comparative modeling of protein refers to constructing an atomic-resolution model of the "target" protein from its amino acid sequence and an experimental

More information

Water models in classical simulations

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

The Effect of Model Internal Flexibility Upon NEMD Simulations of Viscosity

The Effect of Model Internal Flexibility Upon NEMD Simulations of Viscosity Draft: September 29, 1999 The Effect of Model Internal Flexibility Upon NEMD Simulations of Viscosity N. G. Fuller 1 and R. L. Rowley 1,2 Abstract The influence of model flexibility upon simulated viscosity

More information

Table 8.2 Detailed Table of Characteristic Infrared Absorption Frequencies

Table 8.2 Detailed Table of Characteristic Infrared Absorption Frequencies Table 8.2 Detailed Table of Characteristic Infrared Absorption Frequencies The hydrogen stretch region (3600 2500 cm 1 ). Absorption in this region is associated with the stretching vibration of hydrogen

More information

AB INITIO MODELING OF THE STRUCTURAL DEFECTS IN AMIDES

AB INITIO MODELING OF THE STRUCTURAL DEFECTS IN AMIDES Int. J. Chem. Sci.: 9(4), 2011, 1564-1568 ISSN 0972-768X www.sadgurupublications.com AB INITIO MODELING OF THE STRUCTURAL DEFECTS IN AMIDES M. FATHIMA BEGUM, HEMA TRESA VARGHESE a, Y. SHEENA MARY a, C.

More information

Supporting Information

Supporting Information Projection of atomistic simulation data for the dynamics of entangled polymers onto the tube theory: Calculation of the segment survival probability function and comparison with modern tube models Pavlos

More information

Algebraic Study of Stretching and Bending Modes in Linear Tetra-atomic Molecules: HCCCl

Algebraic Study of Stretching and Bending Modes in Linear Tetra-atomic Molecules: HCCCl The African Review of Physics (2013) 8:0016 99 Algebraic Study of Stretching and Bending Modes in Linear Tetra-atomic Molecules: HCCCl Kamal Ziadi * Department of Chemistry, Faculty of Science, University

More information

Chapter 11 Molecular Mechanics

Chapter 11 Molecular Mechanics Chapter 11 Molecular Mechanics Molecular Mechanics uses an analytical, differentiable, and relatively simple potential energy function, (R), for describing the interactions between a set of atoms specified

More information

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

On the accurate calculation of the dielectric constant and the diffusion coefficient from molecular dynamics simulations: the case of SPC/E water

On the accurate calculation of the dielectric constant and the diffusion coefficient from molecular dynamics simulations: the case of SPC/E water On the accurate calculation of the dielectric constant and the diffusion coefficient from molecular dynamics simulations: the case of SPC/E water Orsolya Gereben and László Pusztai Research Institute for

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

Extended examples and description of various jcompoundmapper fingerprints in string format

Extended examples and description of various jcompoundmapper fingerprints in string format Extended examples and description of various jcompoundmapper fingerprints in string format Figure 1: The geometry and topology of Oxaceprol. Pharmacophore types shown in the 3D structure are 1=[L], 3=[A

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

Structural, Electronic, Topological and Vibrational Properties of Isomers of Thymidine - A Study Combining Spectroscopic Data with DFT Calculations

Structural, Electronic, Topological and Vibrational Properties of Isomers of Thymidine - A Study Combining Spectroscopic Data with DFT Calculations Human Journals Research Article December 2017 Vol.:8, Issue:2 All rights are reserved by Silvia Antonia Brandán et al. Structural, Electronic, Topological and Vibrational Properties of Isomers of Thymidine

More information

Supporting Information: Improved Parametrization. of Lithium, Sodium, Potassium, and Magnesium ions. for All-Atom Molecular Dynamics Simulations of

Supporting Information: Improved Parametrization. of Lithium, Sodium, Potassium, and Magnesium ions. for All-Atom Molecular Dynamics Simulations of Supporting Information: Improved Parametrization of Lithium, Sodium, Potassium, and Magnesium ions for All-Atom Molecular Dynamics Simulations of Nucleic Acid Systems. Jejoong Yoo and Aleksei Aksimentiev,,

More information

The Long and Short of Weak Hydrogen Bonds

The Long and Short of Weak Hydrogen Bonds The Long and Short of Weak Hydrogen Bonds Eluvathingal D. Jemmis Department of Inorganic and Physical Chemistry Indian Institute of Science Bangalore 560012, India jemmis@ipc.iisc.ernet.in http://ipc.iisc.ernet.in/~jemmis

More information

Dissecting the stabilization of iodide at the air water interface into components: A free energy analysis

Dissecting the stabilization of iodide at the air water interface into components: A free energy analysis Chemical Physics Letters 42 (26) 199 23 www.elsevier.com/locate/cplett Dissecting the stabilization of iodide at the air water interface into components: A free energy analysis Georgios Archontis a, *,

More information

Dynamical Fluctuating Charge Force Fields: The Aqueous Solvation of Amides

Dynamical Fluctuating Charge Force Fields: The Aqueous Solvation of Amides 672 J. Am. Chem. Soc. 1996, 118, 672-679 Dynamical Fluctuating Charge Force Fields: The Aqueous Solvation of Amides Steven W. ick and B. J. Berne* Contribution from the Department of Chemistry and Center

More information

Photoelectron Spectroscopy of the Hydroxymethoxide Anion, H 2 C(OH)O

Photoelectron Spectroscopy of the Hydroxymethoxide Anion, H 2 C(OH)O Supplementary Material for: Photoelectron Spectroscopy of the Hydroxymethoxide Anion, H 2 C(OH)O Allan M. Oliveira, Julia H. Lehman, Anne B. McCoy 2 and W. Carl Lineberger JILA and Department of Chemistry

More information

POTENTIAL FUNCTIONS FOR HYDROGEN BONDS IN PROTEIN STRUCTURE PREDICTION AND DESIGN. By ALEXANDRE V. MOROZOV* AND TANJA KORTEMME {

POTENTIAL FUNCTIONS FOR HYDROGEN BONDS IN PROTEIN STRUCTURE PREDICTION AND DESIGN. By ALEXANDRE V. MOROZOV* AND TANJA KORTEMME { Comp. by:jaykumar Date:7/11/05 Time:15:40:53 Stage:First Proof File Path:// POTENTIAL FUNCTIONS FOR HYDROGEN BONDS IN PROTEIN STRUCTURE PREDICTION AND DESIGN By ALEXANDRE V. MOROZOV* AND TANJA KORTEMME

More information

Supporting Information

Supporting Information Supporting Information Ionic Liquid Designed for PEDOT:PSS Conductivity Enhancement Ambroise de Izarra, a,b,1 Seongjin Park, a,1 Jinhee Lee, a Yves Lansac, b,c, * Yun Hee Jang a, * a Department of Energy

More information

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

Density Functional Theory: from theory to Applications

Density Functional Theory: from theory to Applications Density Functional Theory: from theory to Applications Uni Mainz January 17, 2011 CP-PAW/COSMO Interface Mixed quantum/ classical molecular dynamics (QM/MM) Quantum mechanics is computationally expensive

More information

Infrared spectra of small biomolecules from first-principle molecular dynamics simulations and effective normal mode analysis

Infrared spectra of small biomolecules from first-principle molecular dynamics simulations and effective normal mode analysis Infrared spectra of small biomolecules from first-principle molecular dynamics simulations and effective normal mode analysis R. Vuilleumier, M.-P. Gaigeot and D. Borgis Département de chimie, Ecole Normale

More information

NUMERICAL SIMULATION OF THE IR SPECTRA OF DNA BASES

NUMERICAL SIMULATION OF THE IR SPECTRA OF DNA BASES BIOPHYSICS NUMERICAL SIMULATION OF THE IR SPECTRA OF DNA BASES C. I. MORARI, CRISTINA MUNTEAN National Institute of Research and Development for Isotopic and Molecular Technologies P.O. Box 700, R-400293

More information

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

4 th Advanced in silico Drug Design KFC/ADD Molecular Modelling Intro. Karel Berka, Ph.D.

4 th Advanced in silico Drug Design KFC/ADD Molecular Modelling Intro. Karel Berka, Ph.D. 4 th Advanced in silico Drug Design KFC/ADD Molecular Modelling Intro Karel Berka, Ph.D. UP Olomouc, 21.1.-25.1. 2019 Motto A theory is something nobody believes, except the person who made it An experiment

More information

Ion-Gated Gas Separation through Porous Graphene

Ion-Gated Gas Separation through Porous Graphene Online Supporting Information for: Ion-Gated Gas Separation through Porous Graphene Ziqi Tian, Shannon M. Mahurin, Sheng Dai,*,, and De-en Jiang *, Department of Chemistry, University of California, Riverside,

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

An Empirical Polarizable Force Field Based on the Classical Drude Oscillator Model: Development History and Recent Applications

An Empirical Polarizable Force Field Based on the Classical Drude Oscillator Model: Development History and Recent Applications This is an open access article published under an ACS AuthorChoice License, which permits copying and redistribution of the article or any adaptations for non-commercial purposes. pubs.acs.org/cr An Empirical

More information

THE JOURNAL OF CHEMICAL PHYSICS 125,

THE JOURNAL OF CHEMICAL PHYSICS 125, THE JOURNAL OF CHEMICAL PHYSICS 125, 194510 2006 Correlated hydrogen bonding fluctuations and vibrational cross peaks in N-methyl acetamide: Simulation based on a complete electrostatic density functional

More information

Ab Initio Modeling of Amide Vibrational Bands in Solution

Ab Initio Modeling of Amide Vibrational Bands in Solution 10794 J. Phys. Chem. A 2004, 108, 10794-10800 Ab Initio Modeling of Amide Vibrational Bands in Solution Nicholas A. Besley* School of Chemistry, UniVersity of Nottingham, UniVersity Park, Nottingham NG7

More information

The Acidity Constant, K a

The Acidity Constant, K a The Acidity Constant, K a The strength of an acid is represented by its ionization constant (acidity constant), K a K a = product of concentrations of ionized species concentration of intact acid K a =

More information

Project 1 Report File: Chem4PB3_project_1_2017-solutions last changed: 02-Feb-2017

Project 1 Report File: Chem4PB3_project_1_2017-solutions last changed: 02-Feb-2017 Project 1 Report File: Chem4PB3_project_1_2017-solutions last changed: 02-Feb-2017 1. Formaldehyde comparison of results from different methods Yellow shaded boxes are closest to experimental Method #e-

More information

An Informal AMBER Small Molecule Force Field :

An Informal AMBER Small Molecule Force Field : An Informal AMBER Small Molecule Force Field : parm@frosst Credit Christopher Bayly (1992-2010) initiated, contributed and lead the efforts Daniel McKay (1997-2010) and Jean-François Truchon (2002-2010)

More information

Chem 341 Organic Chemistry I Lecture Summary 10 September 14, 2007

Chem 341 Organic Chemistry I Lecture Summary 10 September 14, 2007 Chem 34 Organic Chemistry I Lecture Summary 0 September 4, 007 Chapter 4 - Stereochemistry of Alkanes and Cycloalkanes Conformations of Cycloalkanes Cyclic compounds contain something we call Ring Strain.

More information

Biochemistry Prof. S. DasGupta Department of Chemistry Indian Institute of Technology Kharagpur. Lecture - 06 Protein Structure IV

Biochemistry Prof. S. DasGupta Department of Chemistry Indian Institute of Technology Kharagpur. Lecture - 06 Protein Structure IV Biochemistry Prof. S. DasGupta Department of Chemistry Indian Institute of Technology Kharagpur Lecture - 06 Protein Structure IV We complete our discussion on Protein Structures today. And just to recap

More information

Derivation of a molecular mechanics force field for cholesterol*

Derivation of a molecular mechanics force field for cholesterol* Pure Appl. Chem., Vol. 76, No. 1, pp. 189 196, 2004. 2004 IUPAC Derivation of a molecular mechanics force field for cholesterol* Zoe Cournia, Andrea C. Vaiana, G. Matthias Ullmann, and Jeremy C. Smith

More information

T6.2 Molecular Mechanics

T6.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 information

STRUCTURAL DEFECTS IN IMIDATES : AN AB INITIO STUDY

STRUCTURAL DEFECTS IN IMIDATES : AN AB INITIO STUDY Int. J. Chem. Sci.: 9(4), 2011, 1763-1767 ISSN 0972-768X www.sadgurupublications.com STRUCTURAL DEFECTS IN IMIDATES : AN AB INITIO STUDY M. FATHIMA BEGUM, HEMA TRESA VARGHESE a, Y. SHEENA MARY a, C. YOHANNAN

More information

Introduction to Classical Molecular Dynamics. Giovanni Chillemi HPC department, CINECA

Introduction to Classical Molecular Dynamics. Giovanni Chillemi HPC department, CINECA Introduction to Classical Molecular Dynamics Giovanni Chillemi g.chillemi@cineca.it HPC department, CINECA MD ingredients Coordinates Velocities Force field Topology MD Trajectories Input parameters Analysis

More information

Vibrations. Matti Hotokka

Vibrations. Matti Hotokka Vibrations Matti Hotokka Identify the stuff I ve seen this spectrum before. I know what the stuff is Identify the stuff Let s check the bands Film: Polymer Aromatic C-H Aliphatic C-H Group for monosubstituted

More information

CHARMM General Force Field: A Force Field for Drug-Like Molecules Compatible with the CHARMM All-Atom Additive Biological Force Fields

CHARMM General Force Field: A Force Field for Drug-Like Molecules Compatible with the CHARMM All-Atom Additive Biological Force Fields CHARMM General Force Field: A Force Field for Drug-Like Molecules Compatible with the CHARMM All-Atom Additive Biological Force Fields K. VANOMMESLAEGHE, E. HATCHER, C. ACHARYA, S. KUNDU, S. ZHONG, J.

More information