Supporting Information

Size: px
Start display at page:

Download "Supporting Information"

Transcription

1 Supporting Information Energies and 2 -hydroxyl group orientations of RNA backbone conformations. Benchmark CCSD(T)/CBS database, electronic analysis and assessment of DFT methods and MD simulations Arnošt Mládek,* 1 Pavel Banáš, 2 Petr Jurečka, 2 Michal Otyepka, 2 Marie Zgarbová, 2 and Jiří Šponer* 1,3 1 Institute of Biophysics, Academy of Sciences of the Czech Republic, Královopolská 135, Brno, Czech Republic. 2 Regional Centre of Advanced Technologies and Materials, Department of Physical Chemistry, Faculty of Science, Palacky University, tr. 17. listopadu 12, , Olomouc, Czech Republic. 3 CEITEC Central European Institute of Technology, Campus Bohunice, Kamenice 5, Brno, Czech Republic. * Corresponding authors: arnost.mladek@gmail.com, sponer@ncbr.muni.cz

2 Table S1 List of solvent relative energies (kcal.mol -1 ) referenced against the 1a conformation. Applied methods are tagged as follows: M1: B-LYP, M2: PBE, M3: revpbe, M4: TPSS, M5: otpss, M6: B3-LYP, M7: PW6B95, M8: MPW1B95, M9: PWPB95, M10: DSD-B- LYP, M11: MP2(D,T), and M12: MP2(T,Q). For further details see Methods. MAD and MAX (both in kcal.mol -1 ) represent the mean absolute deviation and the maximum absolute deviation from the reference CBS(T) energies. MM labels force field calculations. Listed DFT energies include D3 correction; MAD and MAX values are quoted for DFs with and without D3 term. Label GGA meta-gga Hybrid Double-Hyb. MP2/CBS M1 M2 M3 M4 M5 M6 M7 M8 M9 M10 M11 M12 CBS(T) MM 1a m L &a a a a g d d d e c f j b [ b z z p t q o r a a a #a g g d d d h n i n j [ b b p p s o MAD- D MAX- D N/A N/A N/A MAD MAX

3 Table S2 List of O2 H group acceptors and respective H-bond strengths estimated by AIM analysis for RNA backbone conformations. The (i+1) designates that respective H-bond acceptor belongs to the i+1-th (i.e. 5-3 downstream) nucleotide. Electronic densities (ρ) and magnitudes of electronic density Laplacians ( ρ) at the H-bond critical point are given in a.u. δ, δ + 1, γ + 1 Group Label ρ ρ O2 H Acceptor 1a O3 1m O3 1L O3 &a O5 (i+1) O4 (i+1) 3, 3, p 7a O3 3a O3 Group I 9a O3 1g O3 7d O3 3d O3 5d O3 1e O4 (i+1) 3, 3, t 1c O5 (i+1) 1f O5 (i+1) 5j O3 1b O3 1[ O3 3, 2, p 3b O3 1z O3 Group II 5z O2P (i+1) 7p O3 3, 2, t 1t O5 (i+1) 5q O3 3, 2, m 1o O3 7r O3 2a O3 4a O3 0a O1P (i+1) O3 2, 3, p #a O1P (i+1) 4g O3 Group III 6g O3 8d O3 4d O3 6d O3 2h O3 2, 3, t 4n O1P (i+1) 0i O3 6n O3

4 6j O3 2[ O3 4b O3 2, 2, p 0b O3 Group IV 4p O3 6p O3 2, 2, t 4s O3 2, 2, m 2o O3 Comment on MD simulation protocol The simulation protocol was as follows both systems (SRL and glms, for more details see main text) were minimized, by first optimizing the waters and ions, while the positions of the RNA molecule remained constrained. Subsequently, all RNA atoms were frozen and the solvent molecules with counter-ions were allowed to move during a 1000-step minimization, followed by 10 ps long MD runs under [NpT] conditions (p = 1 atm., T = K) in order to relax the total density. After this, the solute and nucleobases were relaxed by several minimization runs, with decreasing force constants applied to the sugar-phosphate backbone atoms. After the relaxation, each system was heated from K to K for 100 ps. The particle-mesh Ewald (PME) method for treating electrostatic interactions was used, and all simulations were performed under periodic boundary conditions in the [NpT] ensemble at K using weak-coupling Berendsen thermostat with coupling time of 1 ps. The SHAKE algorithm, with a tolerance of 10 5 Å, was used to fix the positions of all hydrogen atoms, and a 10.0 Å cutoff was applied to non-bonding interactions to allow a 2 fs integration step. The simulation lengths of production phases of the simulations were 1 µs.

5 Figure S1 Hydroxyl group rotation energy profiles (kcal.mol -1 ). The zero energy refers to C3 -C2 -O2 -O2 H torsion of 180. The sequence of conformations corresponds to that in Table 1 in the main text. For methodology details, see Methods in the main text.

6

7

8

9

10

11

12

13

14

15

16

17 Cartesian coordinates (xyz format) of the conformations optimized in gas phase with respective relative CBS(T) energies (kcal.mol -1 ). System: 1a, CBS(T): 0.00 C C O C C O C O P O O O C C O C C C O H H H H H H H H H H H H H H H H H H System: 1m, CBS(T): C C O C C O C O P O O

18 O C C O C C C O H H H H H H H H H H H H H H H H H H System: 1L, CBS(T): 2.60 C C O C C O C O P O O O C C O C C C O H H H H H H H H H

19 H H H H H H H H H System: &a, CBS(T): C C O C C O C O P O O O C C O C C C O H H H H H H H H H H H H H H H H H H System: 7a, CBS(T): 1.20 C C O C C

20 O C O P O O O C C O C C C O H H H H H H H H H H H H H H H H H H System: 3a, CBS(T): 0.64 C C O C C O C O P O O O C C O C C C O H H H

21 H H H H H H H H H H H H H H H System: 9a, CBS(T): 0.56 C C O C C O C O P O O O C C O C C C O H H H H H H H H H H H H H H H H H H

22 System: 1g, CBS(T): 0.89 C C O C C O C O P O O O C C O C C C O H H H H H H H H H H H H H H H H H H System: 7d, CBS(T): C C O C C O C O P O O O C C O C

23 C C O H H H H H H H H H H H H H H H H H H System: 3d, CBS(T): 1.04 C C O C C O C O P O O O C C O C C C O H H H H H H H H H H H H H H

24 H H H H System: 5d, CBS(T): 0.60 C C O C C O C O P O O O C C O C C C O H H H H H H H H H H H H H H H H H H System: 1e, CBS(T): 1.92 C C O C C O C O P O

25 O O C C O C C C O H H H H H H H H H H H H H H H H H H System: 1c, CBS(T): 1.22 C C O C C O C O P O O O C C O C C C O H H H H H H H H

26 H H H H H H H H H H System: 1f, CBS(T): 2.15 C C O C C O C O P O O O C C O C C C O H H H H H H H H H H H H H H H H H H System: 5j, CBS(T): 4.75 C C O C

27 C O C O P O O O C C O C C C O H H H H H H H H H H H H H H H H H H System: 1b, CBS(T): C C O C C O C O P O O O C C O C C C O H H

28 H H H H H H H H H H H H H H H H System: 1[, CBS(T): C C O C C O C O P O O O C C O C C C O H H H H H H H H H H H H H H H H H H

29 System: 3b, CBS(T): C C O C C O C O P O O O C C O C C C O H H H H H H H H H H H H H H H H H H System: 1z, CBS(T): C C O C C O C O P O O O C C O

30 C C C O H H H H H H H H H H H H H H H H H H System: 5z, CBS(T): C C O C C O C O P O O O C C O C C C O H H H H H H H H H H H H H

31 H H H H H System: 7p, CBS(T): C C O C C O C O P O O O C C O C C C O H H H H H H H H H H H H H H H H H H System: 1t, CBS(T): 2.45 C C O C C O C O P

32 O O O C C O C C C O H H H H H H H H H H H H H H H H H H System: 5q, CBS(T): 3.19 C C O C C O C O P O O O C C O C C C O H H H H H H H

33 H H H H H H H H H H H System: 1o, CBS(T): 0.11 C C O C C O C O P O O O C C O C C C O H H H H H H H H H H H H H H H H H H System: 7r, CBS(T): 0.77 C C O

34 C C O C O P O O O C C O C C C O H H H H H H H H H H H H H H H H H H System: 2a, CBS(T): C C O C C O C O P O O O C C O C C C O H

35 H H H H H H H H H H H H H H H H H System: 4a, CBS(T): 1.06 C C O C C O C O P O O O C C O C C C O H H H H H H H H H H H H H H H H H H

36 System: 0a, CBS(T): 0.52 C C O C C O C O P O O O C C O C C C O H H H H H H H H H H H H H H H H H H System: #a, CBS(T): 0.05 C C O C C O C O P O O O C C

37 O C C C O H H H H H H H H H H H H H H H H H H System: 4g, CBS(T): 3.62 C C O C C O C O P O O O C C O C C C O H H H H H H H H H H H H

38 H H H H H H System: 6g, CBS(T): C C O C C O C O P O O O C C O C C C O H H H H H H H H H H H H H H H H H H System: 8d, CBS(T): 1.69 C C O C C O C O

39 P O O O C C O C C C O H H H H H H H H H H H H H H H H H H System: 4d, CBS(T): 1.63 C C O C C O C O P O O O C C O C C C O H H H H H H

40 H H H H H H H H H H H H System: 6d, CBS(T): 1.99 C C O C C O C O P O O O C C O C C C O H H H H H H H H H H H H H H H H H H System: 2h, CBS(T): 3.10 C C

41 O C C O C O P O O O C C O C C C O H H H H H H H H H H H H H H H H H H System: 4n, CBS(T): 6.51 C C O C C O C O P O O O C C O C C C O

Supporting Information. Chemical Feasibility of the General Acid/Base Mechanism of glms Ribozyme Self-Cleavage

Supporting Information. Chemical Feasibility of the General Acid/Base Mechanism of glms Ribozyme Self-Cleavage Supporting Information to Chemical Feasibility of the General Acid/Base Mechanism of glms Ribozyme Self-Cleavage Matúš Dubecký, 1 Nils G. Walter, 2 Jiří Šponer, 3,4 Michal Otyepka 1 and Pavel Banáš *,1,3

More information

Tailoring the Properties of Quadruplex Nucleobases for Biological and Nanomaterial Applications

Tailoring the Properties of Quadruplex Nucleobases for Biological and Nanomaterial Applications Electronic Supplementary Material (ESI) for Physical Chemistry Chemical Physics. This journal is the Owner Societies 2014 Supporting Information for: Tailoring the Properties of Quadruplex Nucleobases

More information

Nucleotides containing variously modified sugars: energetics, structure, and mechanical properties

Nucleotides containing variously modified sugars: energetics, structure, and mechanical properties Electronic Supplementary Material (ESI) for Physical Chemistry Chemical Physics. This journal is the Owner Societies 2015 ELECTRONIC SUPPLEMENTARY INFORMATION Nucleotides containing variously modified

More information

Structure and Dynamics of Chitin Nanofibrils in Aqueous Environment Revealed by Molecular Dynamics Simulations

Structure and Dynamics of Chitin Nanofibrils in Aqueous Environment Revealed by Molecular Dynamics Simulations Electronic Supplementary Material (ESI) for RSC Advances. This journal is The Royal Society of Chemistry 2016 Structure and Dynamics of Chitin Nanofibrils in Aqueous Environment Revealed by Molecular Dynamics

More information

Performance of Molecular Mechanics Force Fields for RNA Simulations: Stability of UUCG and GNRA Hairpins

Performance of Molecular Mechanics Force Fields for RNA Simulations: Stability of UUCG and GNRA Hairpins 3836 J. Chem. Theory Comput. 2010, 6, 3836 3849 Performance of Molecular Mechanics Force Fields for RNA Simulations: Stability of UUCG and GNRA Hairpins Pavel Banáš,, Daniel Hollas, Marie Zgarbová, Petr

More information

SANDRO BOTTARO, PAVEL BANÁŠ, JIŘÍ ŠPONER, AND GIOVANNI BUSSI

SANDRO BOTTARO, PAVEL BANÁŠ, JIŘÍ ŠPONER, AND GIOVANNI BUSSI FREE ENERGY LANDSCAPE OF GAGA AND UUCG RNA TETRALOOPS. SANDRO BOTTARO, PAVEL BANÁŠ, JIŘÍ ŠPONER, AND GIOVANNI BUSSI Contents 1. Supplementary Text 1 2 Sample PLUMED input file. 2. Supplementary Figure

More information

With the advance of fast computers, explicit solvent

With the advance of fast computers, explicit solvent pubs.acs.org/jpcl Terms of Use Molecular Dynamics Simulations of Nucleic Acids. From Tetranucleotides to the Ribosome Jir í S poner,*,, Pavel Banaś, Petr Jurecǩa, Marie Zgarbova, Petra Kuḧrova, Marek Havrila,,

More information

Universal Repulsive Contribution to the. Solvent-Induced Interaction Between Sizable, Curved Hydrophobes: Supporting Information

Universal Repulsive Contribution to the. Solvent-Induced Interaction Between Sizable, Curved Hydrophobes: Supporting Information Universal Repulsive Contribution to the Solvent-Induced Interaction Between Sizable, Curved Hydrophobes: Supporting Information B. Shadrack Jabes, Dusan Bratko, and Alenka Luzar Department of Chemistry,

More information

Potential Energy Surface and Binding Energy in External Electric Field: Modulation of Anion π Interactions for Graphene Based Receptors

Potential Energy Surface and Binding Energy in External Electric Field: Modulation of Anion π Interactions for Graphene Based Receptors Potential Energy Surface and Binding Energy in External Electric Field: Modulation of Anion π Interactions for Graphene Based Receptors Cina Foroutan Nejad a and Radek Marek a,b a National Center for Biomolecular

More information

Interaction energy of DNA base pairs and amino acid pairs: DFT, DFTB and WFT calculations

Interaction energy of DNA base pairs and amino acid pairs: DFT, DFTB and WFT calculations Interaction energy of DNA base pairs and amino acid pairs: DFT, DFTB and WFT calculations Tomáš Kubař, Petr Jurečka, Jirka Černý, Honza Řezáč, Michal Otyepka, Haydée Valdés and Pavel Hobza* Institute of

More information

Thiourea Derivatives as Brønsted Acid Organocatalysts

Thiourea Derivatives as Brønsted Acid Organocatalysts Supporting Information Thiourea Derivatives as Brønsted Acid Organocatalysts Ádám Madarász, Zsolt Dósa, Szilárd Varga, * Tibor Soós, Antal Csámpai, Imre Pápai * Institute of Organic Chemistry, Research

More information

Hands-on : Model Potential Molecular Dynamics

Hands-on : Model Potential Molecular Dynamics Hands-on : Model Potential Molecular Dynamics OUTLINE 0. DL_POLY code introduction 0.a Input files 1. THF solvent molecule 1.a Geometry optimization 1.b NVE/NVT dynamics 2. Liquid THF 2.a Equilibration

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

Advanced in silico drug design

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

Introduction. Claudio A. Morgado,, Petr Jurečka,, Daniel Svozil, Pavel Hobza,, and

Introduction. Claudio A. Morgado,, Petr Jurečka,, Daniel Svozil, Pavel Hobza,, and 1524 J. Chem. Theory Comput. 2009, 5, 1524 1544 Balance of Attraction and Repulsion in Nucleic-Acid Base Stacking: CCSD(T)/Complete-Basis-Set-Limit Calculations on Uracil Dimer and a Comparison with the

More information

ARTICLE IN PRESS. Methods xxx (2009) xxx xxx. Contents lists available at ScienceDirect. Methods. journal homepage:

ARTICLE IN PRESS. Methods xxx (2009) xxx xxx. Contents lists available at ScienceDirect. Methods. journal homepage: Methods xxx (2009) xxx xxx Contents lists available at ScienceDirect Methods journal homepage: www.elsevier.com/locate/ymeth Review Article Theoretical studies of RNA catalysis: Hybrid QM/MM methods and

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

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

Medical Research, Medicinal Chemistry, University of Leuven, Leuven, Belgium.

Medical Research, Medicinal Chemistry, University of Leuven, Leuven, Belgium. Supporting Information Towards peptide vaccines against Zika virus: Immunoinformatics combined with molecular dynamics simulations to predict antigenic epitopes of Zika viral proteins Muhammad Usman Mirza

More information

Folding of guanine quadruplex molecules funnel-like mechanism or kinetic partitioning? An overview from MD simulation studies.

Folding of guanine quadruplex molecules funnel-like mechanism or kinetic partitioning? An overview from MD simulation studies. Folding of guanine quadruplex molecules funnel-like mechanism or kinetic partitioning? An overview from MD simulation studies. Jiří Šponer, a,b * Giovanni Bussi, c Petr Stadlbauer, a,b Petra Kührová, b

More information

Routine access to millisecond timescale events with accelerated molecular dynamics

Routine access to millisecond timescale events with accelerated molecular dynamics Routine access to millisecond timescale events with accelerated molecular dynamics Levi C.T. Pierce, Romelia Salomon-Ferrer, Cesar Augusto F. de Oliveira #, J. Andrew McCammon #, Ross C. Walker * SUPPORTING

More information

PCCP Accepted Manuscript

PCCP Accepted Manuscript PCCP Accepted Manuscript This is an Accepted Manuscript, which has been through the Royal Society of Chemistry peer review process and has been accepted for publication. Accepted Manuscripts are published

More information

- 1 - Institute of Organic Chemistry and Biochemistry AS CR, v.v.i., Flemingovo náměstí 2, CZ, Praha, Czech Republic

- 1 - Institute of Organic Chemistry and Biochemistry AS CR, v.v.i., Flemingovo náměstí 2, CZ, Praha, Czech Republic Electronic Supplementary Material (ESI) for RSC Advances. This journal is The Royal Society of Chemistry 2014-1 - The activation of N-glycosidic bond cleavage operated by base-excision repair enzyme hogg1;

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

Computational Modeling of Protein Kinase A and Comparison with Nuclear Magnetic Resonance Data

Computational Modeling of Protein Kinase A and Comparison with Nuclear Magnetic Resonance Data Computational Modeling of Protein Kinase A and Comparison with Nuclear Magnetic Resonance Data ABSTRACT Keyword Lei Shi 1 Advisor: Gianluigi Veglia 1,2 Department of Chemistry 1, & Biochemistry, Molecular

More information

Force Field for Water Based on Neural Network

Force Field for Water Based on Neural Network Force Field for Water Based on Neural Network Hao Wang Department of Chemistry, Duke University, Durham, NC 27708, USA Weitao Yang* Department of Chemistry, Duke University, Durham, NC 27708, USA Department

More information

Hairpins Participating in Folding of Human Telomeric Sequence Quadruplexes Studied by Standard and T-REMD Simulations

Hairpins Participating in Folding of Human Telomeric Sequence Quadruplexes Studied by Standard and T-REMD Simulations Hairpins Participating in Folding of Human Telomeric Sequence Quadruplexes Studied by Standard and T-REMD Simulations Petr Stadlbauer 1, Petra Kührová 2, Pavel Banáš 2, Jaroslav Koča 3,4, Giovanni Bussi

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

Type II Kinase Inhibitors Show an Unexpected Inhibition Mode against Parkinson s Disease-Linked LRRK2 Mutant G2019S.

Type II Kinase Inhibitors Show an Unexpected Inhibition Mode against Parkinson s Disease-Linked LRRK2 Mutant G2019S. Type II Kinase Inhibitors Show an Unexpected Inhibition Mode against Parkinson s Disease-Linked LRRK2 Mutant G219S. Min Liu@&*, Samantha A. Bender%*, Gregory D Cuny@, Woody Sherman, Marcie Glicksman@ Soumya

More information

Supplementary Information for Electronic signature of the instantaneous asymmetry in the first coordination shell in liquid water

Supplementary Information for Electronic signature of the instantaneous asymmetry in the first coordination shell in liquid water Supplementary Information for Electronic signature of the instantaneous asymmetry in the first coordination shell in liquid water Thomas D. Kühne 1, 2 and Rustam Z. Khaliullin 1, 1 Institute of Physical

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

Molecular dynamics simulations of a single stranded (ss) DNA

Molecular dynamics simulations of a single stranded (ss) DNA Molecular dynamics simulations of a single stranded (ss) DNA Subhasish Chatterjee 1, Bonnie Gersten 1, Siddarth Thakur 2, Alexander Burin 2 1 Department of Chemistry, Queens College and the Graduate Center

More information

Applications of the van der Waals density functional: DNA and metal organic framework materials. Chalmers

Applications of the van der Waals density functional: DNA and metal organic framework materials. Chalmers Applications of the van der Waals density functional: DNA and metal organic framework materials Valentino R. Cooper Chalmers Svetla Chakarova-Käck Per Hyldgaard Jesper Kleis Henrik Rydberg Elsebeth Schröder

More information

Journal of Pharmacology and Experimental Therapy-JPET#172536

Journal of Pharmacology and Experimental Therapy-JPET#172536 A NEW NON-PEPTIDIC INHIBITOR OF THE 14-3-3 DOCKING SITE INDUCES APOPTOTIC CELL DEATH IN CHRONIC MYELOID LEUKEMIA SENSITIVE OR RESISTANT TO IMATINIB Manuela Mancini, Valentina Corradi, Sara Petta, Enza

More information

XYZ file format Protein Data Bank (pdb) file format Z - matrix

XYZ file format Protein Data Bank (pdb) file format Z - matrix Chemistry block (exercise 1) In this exercise, students will be introduced how to preform simple quantum chemical calculations. Input files for Gaussian09. Output file structure. Geometry optimization,

More information

Doping with Graphitic Nitrogen Triggers Ferromagnetism in Graphene

Doping with Graphitic Nitrogen Triggers Ferromagnetism in Graphene Supporting Information for Doping with Graphitic Nitrogen Triggers Ferromagnetism in Graphene Piotr Błoński, Jiří Tuček, Zdeněk Sofer, Vlastimil Mazánek, Martin Petr, Martin Pumera, Michal Otyepka, and

More information

SUPPLEMENTARY MATERIAL. Supplementary material and methods:

SUPPLEMENTARY MATERIAL. Supplementary material and methods: Electronic Supplementary Material (ESI) for Catalysis Science & Technology. This journal is The Royal Society of Chemistry 2015 SUPPLEMENTARY MATERIAL Supplementary material and methods: - Computational

More information

Towards accurate calculations of Zn 2+ binding free energies in zinc finger proteins

Towards accurate calculations of Zn 2+ binding free energies in zinc finger proteins Towards accurate calculations of Zn 2+ binding free energies in zinc finger proteins UNDERGRADUATE HONORS RESEARCH THESIS Presented in Partial Fulfillment of the Requirements for the Bachelor of Science

More information

Supplementary material. From cellulose to kerogen: molecular simulation. of a geological process

Supplementary material. From cellulose to kerogen: molecular simulation. of a geological process Electronic Supplementary Material (ESI) for Chemical Science. This journal is The Royal Society of Chemistry 2017 Supplementary material From cellulose to kerogen: molecular simulation of a geological

More information

Hairpins participating in folding of human telomeric sequence quadruplexes studied by standard and T-REMD simulations

Hairpins participating in folding of human telomeric sequence quadruplexes studied by standard and T-REMD simulations Nucleic Acids Research Advance Access published October 3, 2015 Nucleic Acids Research, 2015 1 doi: 10.1093/nar/gkv994 Hairpins participating in folding of human telomeric sequence quadruplexes studied

More information

Design of Efficient Catalysts with Double Transition Metal. Atoms on C 2 N Layer

Design of Efficient Catalysts with Double Transition Metal. Atoms on C 2 N Layer Supporting Information Design of Efficient Catalysts with Double Transition Metal Atoms on C 2 N Layer Xiyu Li, 1, Wenhui Zhong, 2, Peng Cui, 1 Jun Li, 1 Jun Jiang 1, * 1 Hefei National Laboratory for

More information

Molecular Dynamics Simulations. Dr. Noelia Faginas Lago Dipartimento di Chimica,Biologia e Biotecnologie Università di Perugia

Molecular Dynamics Simulations. Dr. Noelia Faginas Lago Dipartimento di Chimica,Biologia e Biotecnologie Università di Perugia Molecular Dynamics Simulations Dr. Noelia Faginas Lago Dipartimento di Chimica,Biologia e Biotecnologie Università di Perugia 1 An Introduction to Molecular Dynamics Simulations Macroscopic properties

More information

Supplementary Materials

Supplementary Materials Supplementary Materials Atomistic Origin of Brittle Failure of Boron Carbide from Large Scale Reactive Dynamics Simulations; Suggestions toward Improved Ductility Qi An and William A. Goddard III * Materials

More information

2008 Biowerkzeug Ltd.

2008 Biowerkzeug Ltd. 2008 Biowerkzeug Ltd. 1 Contents Summary...3 1 Simulation...4 1.1 Setup...4 1.2 Output...4 2 Settings...5 3 Analysis...9 3.1 Setup...9 3.2 Input options...9 3.3 Descriptions...10 Please note that we cannot

More information

DFT STUDY OF Co 2+ AND Fe 2+ - URACIL COMPLEXES IN THE GAS PHASE

DFT STUDY OF Co 2+ AND Fe 2+ - URACIL COMPLEXES IN THE GAS PHASE INTERNATIONAL JOURNAL OF RESEARCH IN PHARMACY AND CHEMISTRY Available online at www.ijrpc.com Research Article DFT STUDY OF Co 2+ AND Fe 2+ - URACIL COMPLEXES IN THE GAS PHASE Raghab Parajuli* Department

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

WATER PERMEATION THROUGH GRAPHENE NANOSLIT BY MOLECULAR DYNAMICS SIMULATION

WATER PERMEATION THROUGH GRAPHENE NANOSLIT BY MOLECULAR DYNAMICS SIMULATION WATER PERMEATION THROUGH GRAPHENE NANOSLIT BY MOLECULAR DYNAMICS SIMULATION Taro Yamada 1 and Ryosuke Matsuzaki 2 1 Department of Mechanical Engineering, Tokyo University of Science, 2641 Yamazaki, Noda,

More information

Supporting Information Synthesis of (D)-erythrose from glycolaldehyde aqueous solutions under electric field

Supporting Information Synthesis of (D)-erythrose from glycolaldehyde aqueous solutions under electric field Electronic Supplementary Material (ESI) for ChemComm. This journal is The Royal Society of Chemistry 2018 Supporting Information Synthesis of (D)-erythrose from glycolaldehyde aqueous solutions under electric

More information

Synergy between NMR measurements and MD simulations of protein/rna complexes: application to the RRMs, the most common RNA recognition motifs

Synergy between NMR measurements and MD simulations of protein/rna complexes: application to the RRMs, the most common RNA recognition motifs 6452 6470 Nucleic Acids Research, 2016, Vol. 44, No. 13 Published online 18 May 2016 doi: 10.1093/nar/gkw438 Synergy between NMR measurements and MD simulations of protein/rna complexes: application to

More information

Classification and Energetics of the Base-Phosphate Interactions in RNA

Classification and Energetics of the Base-Phosphate Interactions in RNA Bowling Green State University ScholarWorks@BGSU Chemistry Faculty Publications Chemistry 8-2009 Classification and Energetics of the Base-Phosphate Interactions in RNA Craig L. Zirbel Bowling Green State

More information

DISCRETE TUTORIAL. Agustí Emperador. Institute for Research in Biomedicine, Barcelona APPLICATION OF DISCRETE TO FLEXIBLE PROTEIN-PROTEIN DOCKING:

DISCRETE TUTORIAL. Agustí Emperador. Institute for Research in Biomedicine, Barcelona APPLICATION OF DISCRETE TO FLEXIBLE PROTEIN-PROTEIN DOCKING: DISCRETE TUTORIAL Agustí Emperador Institute for Research in Biomedicine, Barcelona APPLICATION OF DISCRETE TO FLEXIBLE PROTEIN-PROTEIN DOCKING: STRUCTURAL REFINEMENT OF DOCKING CONFORMATIONS Emperador

More information

A conserved P-loop anchor limits the structural dynamics that mediate. nucleotide dissociation in EF-Tu.

A conserved P-loop anchor limits the structural dynamics that mediate. nucleotide dissociation in EF-Tu. Supplemental Material for A conserved P-loop anchor limits the structural dynamics that mediate nucleotide dissociation in EF-Tu. Evan Mercier 1,2, Dylan Girodat 1, and Hans-Joachim Wieden 1 * 1 Alberta

More information

MARTINI simulation details

MARTINI simulation details S1 Appendix MARTINI simulation details MARTINI simulation initialization and equilibration In this section, we describe the initialization of simulations from Main Text section Residue-based coarsegrained

More information

Peptide folding in non-aqueous environments investigated with molecular dynamics simulations Soto Becerra, Patricia

Peptide folding in non-aqueous environments investigated with molecular dynamics simulations Soto Becerra, Patricia University of Groningen Peptide folding in non-aqueous environments investigated with molecular dynamics simulations Soto Becerra, Patricia IMPORTANT NOTE: You are advised to consult the publisher's version

More information

S66: A Well-balanced Database of Benchmark Interaction Energies Relevant to Biomolecular Structures

S66: A Well-balanced Database of Benchmark Interaction Energies Relevant to Biomolecular Structures pubs.acs.org/jctc S66: A Well-balanced Database of Benchmark Interaction Energies Relevant to Biomolecular Structures Jan Rezac,*, Kevin E. Riley, and Pavel Hobza, Institute of Organic Chemistry and Biochemistry,

More information

Methods for van der Waals Interactions

Methods for van der Waals Interactions Methods for van der Waals Interactions Alexandre Tkatchenko Theory Department, Fritz Haber Institut der MPG Berlin, Germany tkatchen@fhi berlin.mpg.de Haber Institute FHI DFT and Beyond Workshop, Jul.

More information

Simulation of molecular systems by molecular dynamics

Simulation of molecular systems by molecular dynamics Simulation of molecular systems by molecular dynamics Yohann Moreau yohann.moreau@ujf-grenoble.fr November 26, 2015 Yohann Moreau (UJF) Molecular Dynamics, Label RFCT 2015 November 26, 2015 1 / 35 Introduction

More information

SUPPLEMENTARY INFORMATION

SUPPLEMENTARY INFORMATION DOI: 10.1038/NCHEM.1680 On the nature and origin of dicationic, charge-separated species formed in liquid water on X-ray irradiation Stephan Thürmer, 1 Milan Ončák, 2 Niklas Ottosson, 3 Robert Seidel,

More information

Supplementary Information

Supplementary Information Electronic Supplementary Material (ESI) for New Journal of Chemistry. This journal is The Royal Society of Chemistry and the Centre National de la Recherche Scientifique 2018 Supplementary Information

More information

Supporting Material for. Microscopic origin of gating current fluctuations in a potassium channel voltage sensor

Supporting Material for. Microscopic origin of gating current fluctuations in a potassium channel voltage sensor Supporting Material for Microscopic origin of gating current fluctuations in a potassium channel voltage sensor J. Alfredo Freites, * Eric V. Schow, * Stephen H. White, and Douglas J. Tobias * * Department

More information

Canadian Journal of Chemistry. The INV24 Test Set: How Well do Quantum-Chemical Methods Describe Inversion and Racemization Barriers?

Canadian Journal of Chemistry. The INV24 Test Set: How Well do Quantum-Chemical Methods Describe Inversion and Racemization Barriers? The INV24 Test Set: How Well do Quantum-Chemical Methods Describe Inversion and Racemization Barriers? Journal: Manuscript ID cjc-2016-0290.r1 Manuscript Type: Article Date Submitted by the Author: 15-Jul-2016

More information

Unit Cell-Level Thickness Control of Single-Crystalline Zinc Oxide Nanosheets Enabled by Electrical Double Layer Confinement

Unit Cell-Level Thickness Control of Single-Crystalline Zinc Oxide Nanosheets Enabled by Electrical Double Layer Confinement Unit Cell-Level Thickness Control of Single-Crystalline Zinc Oxide Nanosheets Enabled by Electrical Double Layer Confinement Xin Yin, Yeqi Shi, Yanbing Wei, Yongho Joo, Padma Gopalan, Izabela Szlufarska,

More information

Molecular Dynamics. A very brief introduction

Molecular Dynamics. A very brief introduction Molecular Dynamics A very brief introduction Sander Pronk Dept. of Theoretical Physics KTH Royal Institute of Technology & Science For Life Laboratory Stockholm, Sweden Why computer simulations? Two primary

More information

Van der Waals interactions in DFT

Van der Waals interactions in DFT Van der Waals interactions in DFT Maxime Dion*, Aaron Puzder*, T. Thonhauser,* Valentino R. Cooper*, Shen Li*, Eamonn Murray, Lingzhu Kong, Kyuho Lee, and David C. Langreth Department of Physics and Astronomy,

More information

Problem Set 1

Problem Set 1 2006 7.012 Problem Set 1 Due before 5 PM on FRIDAY, September 15, 2006. Turn answers in to the box outside of 68-120. PLEASE WRITE YOUR ANSWERS ON THIS PRINTOUT. 1. For each of the following parts, pick

More information

Conformations of Flanking Bases in HIV-1 RNA DIS Kissing Complexes Studied by Molecular Dynamics

Conformations of Flanking Bases in HIV-1 RNA DIS Kissing Complexes Studied by Molecular Dynamics 3932 Biophysical Journal Volume 93 December 2007 3932 3949 Conformations of Flanking Bases in HIV-1 RNA DIS Kissing Complexes Studied by Molecular Dynamics Kamila Réblová,* Eva Fadrná, y Joanna Sarzynska,

More information

Small Details Matter: The 2 -Hydroxyl as a Conformational Switch in RNA

Small Details Matter: The 2 -Hydroxyl as a Conformational Switch in RNA Small Details Matter: The 2 -Hydroxyl as a Conformational Switch in RNA Leonardo Darré 1,2, Ivan Ivani 1,2, Pablo D. Dans 1,2, Hansel Gómez 1,2, Adam Hospital 1,2 and Modesto Orozco 1,2,3 * 1 Institute

More information

Explicit Ion, Implicit Water Solvation for Molecular Dynamics of Nucleic Acids and Highly Charged Molecules

Explicit Ion, Implicit Water Solvation for Molecular Dynamics of Nucleic Acids and Highly Charged Molecules Explicit Ion, Implicit Water Solvation for Molecular Dynamics of Nucleic Acids and Highly Charged Molecules NINAD V. PRABHU, 1 MANORANJAN PANDA, 2 QINGYI YANG, 3 KIM A. SHARP 3 1 Schering-Plough, Kenilworth,

More information

Introductory Biochemistry

Introductory Biochemistry Introductory Biochemistry Instructors Dr. Nafez Abu Tarboush Dr. Mamoun Ahram Recommended textbooks Biochemistry; Mary K. Campbell and Shawn O. Farrell, Brooks Cole; 6 th edition Recommended electronic

More information

Supplementary Information for Atomistic Simulation of Spinodal Phase Separation Preceding Polymer Crystallization

Supplementary Information for Atomistic Simulation of Spinodal Phase Separation Preceding Polymer Crystallization Supplementary Information for Atomistic Simulation of Spinodal Phase Separation Preceding Polymer Crystallization Richard H. Gee * Naida Lacevic and Laurence E. Fried University of California Lawrence

More information

Structures of the Molecular Components in DNA and RNA with Bond Lengths Interpreted as Sums of Atomic Covalent Radii

Structures of the Molecular Components in DNA and RNA with Bond Lengths Interpreted as Sums of Atomic Covalent Radii The Open Structural Biology Journal, 2008, 2, 1-7 1 Structures of the Molecular Components in DNA and RNA with Bond Lengths Interpreted as Sums of Atomic Covalent Radii Raji Heyrovska * Institute of Biophysics

More information

Pseudopotentials for hybrid density functionals and SCAN

Pseudopotentials for hybrid density functionals and SCAN Pseudopotentials for hybrid density functionals and SCAN Jing Yang, Liang Z. Tan, Julian Gebhardt, and Andrew M. Rappe Department of Chemistry University of Pennsylvania Why do we need pseudopotentials?

More information

ACMAC s PrePrint Repository

ACMAC s PrePrint Repository ACMAC s PrePrint Repository Properties of short polystyrene chains confined between two Gold surfaces through a combined Density Functional Theory and classical Molecular Dynamics approach Karen Johnston

More information

MM-PBSA Validation Study. Trent E. Balius Department of Applied Mathematics and Statistics AMS

MM-PBSA Validation Study. Trent E. Balius Department of Applied Mathematics and Statistics AMS MM-PBSA Validation Study Trent. Balius Department of Applied Mathematics and Statistics AMS 535 11-26-2008 Overview MM-PBSA Introduction MD ensembles one snap-shots relaxed structures nrichment Computational

More information

Supporting Information

Supporting Information Electronic Supplementary Material (ESI) for Physical Chemistry Chemical Physics. This journal is the Owner Societies 2016 Supporting Information Lipid molecules can induce an opening of membrane-facing

More information

Role of van der Waals Interactions in Physics, Chemistry, and Biology

Role of van der Waals Interactions in Physics, Chemistry, and Biology Role of van der Waals Interactions in Physics, Chemistry, and Biology How can we describe vdw forces in materials accurately? Failure of DFT Approximations for (Long-Range) Van der Waals Interactions 1

More information

Evolution of Translation: Dynamics of Recognition in RNA:Protein Complexes

Evolution of Translation: Dynamics of Recognition in RNA:Protein Complexes Evolution of Translation: Dynamics of Recognition in RNA:Protein Complexes Zaida (Zan) Luthey-Schulten Dept. Chemistry, Beckman Institute, Biophysics, Institute of Genomics Biology, & Physics NIH Workshop

More information

Ranjit P. Bahadur Assistant Professor Department of Biotechnology Indian Institute of Technology Kharagpur, India. 1 st November, 2013

Ranjit P. Bahadur Assistant Professor Department of Biotechnology Indian Institute of Technology Kharagpur, India. 1 st November, 2013 Hydration of protein-rna recognition sites Ranjit P. Bahadur Assistant Professor Department of Biotechnology Indian Institute of Technology Kharagpur, India 1 st November, 2013 Central Dogma of life DNA

More information

The micro-properties of [hmpy+] [Tf 2 N-] Ionic liquid: a simulation. study. 1. Introduction

The micro-properties of [hmpy+] [Tf 2 N-] Ionic liquid: a simulation. study. 1. Introduction ISBN 978-1-84626-081-0 Proceedings of the 2010 International Conference on Application of Mathematics and Physics Volume 1: Advances on Space Weather, Meteorology and Applied Physics Nanjing, P. R. China,

More information

List of Figures Page Figure No. Figure Caption No. Figure 1.1.

List of Figures Page Figure No. Figure Caption No. Figure 1.1. List of Figures Figure No. Figure Caption Page No. Figure 1.1. Cation- interactions and their modulations. 4 Figure 1.2. Three conformations of benzene dimer, S is not a minimum on the potential energy

More information

LAMMPS Performance Benchmark on VSC-1 and VSC-2

LAMMPS Performance Benchmark on VSC-1 and VSC-2 LAMMPS Performance Benchmark on VSC-1 and VSC-2 Daniel Tunega and Roland Šolc Institute of Soil Research, University of Natural Resources and Life Sciences VSC meeting, Neusiedl am See, February 27-28,

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

Noncovalent interactions in biochemistry

Noncovalent interactions in biochemistry Noncovalent interactions in biochemistry Kevin E. Riley 1 and Pavel Hobza 2,3 Noncovalent interactions are known to play a key role in biochemistry. The knowledge of stabilization (relative) energies and

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

Improved Resolution of Tertiary Structure Elasticity in Muscle Protein

Improved Resolution of Tertiary Structure Elasticity in Muscle Protein Improved Resolution of Tertiary Structure Elasticity in Muscle Protein Jen Hsin and Klaus Schulten* Department of Physics and Beckman Institute, University of Illinois at Urbana-Champaign, Urbana, Illinois

More information

Ab initio molecular dynamics

Ab initio molecular dynamics Ab initio molecular dynamics Kari Laasonen, Physical Chemistry, Aalto University, Espoo, Finland (Atte Sillanpää, Jaakko Saukkoriipi, Giorgio Lanzani, University of Oulu) Computational chemistry is a field

More information

Molecular Dynamics Simulations of Guanine Quadruplex Loops: Advances and Force Field Limitations

Molecular Dynamics Simulations of Guanine Quadruplex Loops: Advances and Force Field Limitations Biophysical Journal Volume 87 July 2004 227 242 227 Molecular Dynamics Simulations of Guanine Quadruplex Loops: Advances and Force Field Limitations Eva Fadrná,* Nad a Špačková, y Richard Štefl,* Jaroslav

More information

NWChem: Molecular Dynamics and QM/MM

NWChem: Molecular Dynamics and QM/MM NWChem: Molecular Dynamics and QM/MM Molecular Dynamics Functionality Target systems: General features: biomolecules (proteins, DNA/RNA, biomembranes) energy evaluation (SP) energy minimization (EM) molecular

More information

Molecular Dynamics. Coherency of Molecular Motions. Local Minima & Global Minimum. Theoretical Basis of Molecular Dynamic Calculations

Molecular Dynamics. Coherency of Molecular Motions. Local Minima & Global Minimum. Theoretical Basis of Molecular Dynamic Calculations Molecular Dynamics Molecular Dynamics provide a good approach to determine the preferred conformers and global minimum of a molecule Theoretical Basis of Molecular Dynamic Calculations Molecular dynamics

More information

Supplementary Information

Supplementary Information Supplementary Information Resveratrol Serves as a Protein-Substrate Interaction Stabilizer in Human SIRT1 Activation Xuben Hou,, David Rooklin, Hao Fang *,,, Yingkai Zhang Department of Medicinal Chemistry

More information

Jorge Kohanoff. DNA damage from first principles

Jorge Kohanoff. DNA damage from first principles DNA damage from first principles Jorge Kohanoff Atomistic Simulation Centre School of Maths and Physics Queen s University Belfast Northern Ireland UK CP2K UK Users Meeting London, 29 th January 2014 Queen

More information

Supplementary Figure 2 FF at concentration of 8 mg/ml dissolved in 100% DMSO. The molecules assemble into plate-like structures. Scale bar is 50 µm.

Supplementary Figure 2 FF at concentration of 8 mg/ml dissolved in 100% DMSO. The molecules assemble into plate-like structures. Scale bar is 50 µm. 1 2 3 4 5 6 7 8 9 10 11 Supplementary Figure 1 Microfluidics device mixing efficiency. Yellow and blue food coloring solutions are introduced to inlets 1 and 2, respectively. After approximately one tenth

More information

I690/B680 Structural Bioinformatics Spring Protein Structure Determination by NMR Spectroscopy

I690/B680 Structural Bioinformatics Spring Protein Structure Determination by NMR Spectroscopy I690/B680 Structural Bioinformatics Spring 2006 Protein Structure Determination by NMR Spectroscopy Suggested Reading (1) Van Holde, Johnson, Ho. Principles of Physical Biochemistry, 2 nd Ed., Prentice

More information

Statistical thermodynamics for MD and MC simulations

Statistical thermodynamics for MD and MC simulations Statistical thermodynamics for MD and MC simulations knowing 2 atoms and wishing to know 10 23 of them Marcus Elstner and Tomáš Kubař 22 June 2016 Introduction Thermodynamic properties of molecular systems

More 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

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

Modeling Ultrafast Deactivation in Oligothiophenes via Nonadiabatic Dynamics

Modeling Ultrafast Deactivation in Oligothiophenes via Nonadiabatic Dynamics Electronic Supplementary Material (ESI) for Physical Chemistry Chemical Physics. This journal is the Owner Societies 2015 Supplementary Data for Modeling Ultrafast Deactivation in Oligothiophenes via Nonadiabatic

More information

Supporting Information

Supporting Information Supporting Information Solid-State 17 O NMR Reveals Hydrogen-Bonding Energetics: Not All Low-Barrier Hydrogen Bonds Are Strong Jiasheng Lu, Ivan Hung, Andreas Brinkmann, Zhehong Gan, Xianqi Kong, and Gang

More information

Diffusion of Water and Diatomic Oxygen in Poly(3-hexylthiophene) Melt: A Molecular Dynamics Simulation Study

Diffusion of Water and Diatomic Oxygen in Poly(3-hexylthiophene) Melt: A Molecular Dynamics Simulation Study Diffusion of Water and Diatomic Oxygen in Poly(3-hexylthiophene) Melt: A Molecular Dynamics Simulation Study Julia Deitz, Yeneneh Yimer, and Mesfin Tsige Department of Polymer Science University of Akron

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