Macrocycle Conformational Sampling by DFT-
|
|
- Tyrone Ernest Tate
- 5 years ago
- Views:
Transcription
1 Supporting Information Macrocycle Conformational Sampling by DFT- D3/COSMO-RS Methodology. Ondrej Gutten, Daniel Bím, Jan Řezáč, Lubomír Rulíšek* The Institute of Organic Chemistry and Biochemistry of the Czech Academy of Sciences, Gilead Sciences Research Center & IOCB, Flemingovo náměstí 2, Praha 6, Czech Republic * 1
2 Individual components of G DFT/COSMO-RS : detailed analysis (Section S3.1) Importance of individual components of the free energy. From the computational perspective, the total free energy is composed of three components, discussed individually below. Knowledge of relative importance of these components allows for designing a protocol with cost and accuracy tailored to available resources and specific requirements. Figure S3 (below) presents the range of values of the components for two of the systems in the test set. The three components, in order of importance, are the gas-phase energy, E gp, the solvation energy, G solv, and the sum of zero-point vibrational energy, entropic, and thermal contributions, µ. The energy spans for these components are loosely related to the size of the system. The total free energy, G DFT/COSMO-RS, which is a sum of the three contributions, exhibits ranges that are significantly smaller, because the individual components are strongly correlated, as discussed below (Figure 4, main text). Gas-Phase Energies: Performance of DFT Functionals. The same conclusion is obtained from the pairwise average absolute differences for several popular DFT functionals (vide supra). The data are reported in the Table SI_test.data. It can be seen that the MAE for any two of the D3-corrected functionals is within 5 kj.mol -1. The same holds true for any pair from the D3-uncorrected set whereas the mutual comparison (D3-corrected vs. D3- uncorrected) yields much larger errors (> 10 kj.mol -1 ). With respect to our calibration work 1 and similar test sets reported by Grimme and coworkers, 2 the D3-corrected functionals should provide considerably more accurate values. Therefore, our choice of functional is driven by computational cost towards pure GGA functionals and BLYP-D3 is used in the remaining part of this work. The BP86-D3 is also a reasonable alternative, since the COMSO-RS protocol has been parametrized for this functional and we have demonstrated in our previous work concerning thermodynamics of metal ion complexation 3,4 that G BP86-D3/COSMO-RS values were 2
3 in the best agreement with experimental stability constants despite the poor performance of BP86 in the in vacuo component of the overall free energy change. The apparent contradiction had been attributed to cancellation of errors between gas-phase thermodynamics and the COSMO-RS solvation energies. Solvation Energies. Previous studies carried out in our group identified COSMO-RS as one of the most robust implicit solvation methods and therefore, it was used as the solvation method of choice in this study. Importantly, we have observed that it is critical to use the so-called FINE parametrizations and the $cosmo_isorad keyword in the COSMO input which prevents rare numerical failures of the protocol for conformers with very small ring cavity. As can be seen in Table SI_test.data, the computed COSMO-RS solvation energies of the conformers of the single compound spanned a seemingly wide range of energies (~100 kj.mol -1 ). Somewhat surprisingly, size of a system is not a decisive factor. A careful inspection of several extreme cases revealed that differences originate, inter alia, in large differences in dipole moments between such conformations, which appear as a result of an alignment of local dipoles corresponding to polar groups in the given conformer. As discussed below, such changes are largely compensated by an opposing change in gas-phase electronic energies. Low-lying vibrational modes and the treatment of imaginary frequencies. In our test set of 120 structures, there are 22 such cases (DZ gas-phase frequency analysis, Table SI_test.data). Upon re-optimizing all structures with finer DFT grid and weighting of the derivatives in calculations of the analytical second derivatives (denoted FINE below), only one system with imaginary frequency remained. In systems with no imaginary frequencies in REGULAR setup (non-problematic cases), the thermal contributions differ almost negligibly from their FINE counterparts (MAE of ~1.3 kj.mol -1 ). Thus, in well converged systems the increased numerical accuracy does not provide 3
4 an attractive tradeoff since it comes at the expense of approximately one order of magnitude greater computational price (REGULAR vs. FINE setup). For the 21 structures with vanishing imaginary frequencies we observe only negligible structural changes after re-optimization and re-calculation in FINE setup. Moreover, the vibrational spectrum is almost identical, the only change being the imaginary frequency turning into a real low-lying frequency (0-200 cm -1 ). A quick and dirty trick that circumvents the need for costly re-optimization is to simply flip an imaginary frequency into its positive value. This approximate solution is also implemented in the Grimme s program thermo which has been employed throughout for the gas-phase thermodynamics. 5 In this small subset of 21 structures, MAE increased to 2.3 kj.mol -1 (again, using FINE setup as reference). 4
5 Table S1: Number of low-energy conformers similar to a crystal structure, as found by a composite protocol and by standalone protocols (the remaining numerical columns). The number of similar structures is judged by energy (<20kJ/mol above global minimum), geometry (<1.0 Å RMSD) or both. True positives represent number of structures found by both the composite protocol and the standalone protocol. False positives represents number of structures found by standalone protocol that are not found to be similar by the composite protocol. False negatives represent number of structures found by the composite protocol but not found by the standalone protocol. Even using just the geometry or just the energy criterium, the number of true positivies is significantly lower for standalone protocols than for the composite protocol. Energy & geometry Geometry only Energy only System Composite True False False True False True False protocol positives positives negatives positives negatives positives negatives CAMVES CHPSAR COHVAW defllmod GS POXTRD SANGLI YIVNOG CAMVES CHPSAR COHVAW extllmod GS
6 POXTRD SANGLI YIVNOG CAMVES CHPSAR COHVAW PM6MD GS POXTRD SANGLI YIVNOG
7 Table S2: Number of unique low-energy conformers found simultaneously by two protocols (intersection). The second row of this table was used to construct Figure 8 corresponds to the size of the overlapping region of the full and dashed lines. 1 2 CAMVES CHPSAR COHVAW GS GS POXTRD SANGLI YIVNOG defllmod PM6MD defllmod extllmod PM6MD extllmod Table S3: Number of unique low-energy conformers found by one or more protocols (union). The first, third, and fifth row of this table was used to construct Figure 8 corresponds to the height of the full and dashed lines and their total height, respectively. 1 2 CAMVES CHPSAR COHVAW GS GS POXTRD SANGLI YIVNOG defllmod PM6MD extllmod defllmod PM6MD defllmod extllmod PM6MD extllmod all three
8 Table S4: Chapman estimate, E(N), of the total number of low-energy conformers based on hypergeometric distribution H(N,M,n,m). N is the total number of low-energy conformers, M is the number of low-energy conformers found by one of the protocols, n is the number of low-energy conformers found by the other protocol, and m is the number of low-energy conformers found by both of these protocols simultaneously. Chapman estimate is calculated as EE(NN) = (MM+1)(nn+1) 1. The (mm+1) second row of this table was used to construct Figure 8 corresponds to the coloured horizontal bar. 1 2 CAMVES CHPSAR COHVAW GS GS POXTRD SANGLI YIVNOG defllmod PM6MD defllmod extllmod PM6MD extllmod Table S5: Lower boundary of a 95%-confidence interval of the total number of low-energy conformers based on hypergeometric distribution H(N,M,n,m). N is the total number of low-energy conformers, M is the number of low-energy conformers found by one of the protocols, n is the number of low-energy conformers found by the other protocol, and m is the number of low-energy conformers found by both of these protocols simultaneously. For details on calculation of confidence intervals see Ref (VI). The second row of this table was used to construct Figure 8 corresponds to the lower end of the coloured bar. 1 2 CAMVES CHPSAR COHVAW GS GS POXTRD SANGLI YIVNOG defllmod PM6MD defllmod extllmod PM6MD extllmod
9 Table S6: Upper boundary of a 95%-confidence interval of the total number of low-energy conformers based on hypergeometric distribution H(N,M,n,m). N is the total number of low-energy conformers, M is the number of low-energy conformers found by one of the protocols, n is the number of low-energy conformers found by the other protocol, and m is the number of low-energy conformers found by both of these protocols simultaneously. For details on calculation of confidence intervals see Ref (VI). The second row of this table was used to construct Figure 8 corresponds to the upper end of the coloured bar. 1 2 CAMVES CHPSAR COHVAW GS GS POXTRD SANGLI YIVNOG defllmod PM6MD infinity 472 infinity 459 defllmod extllmod infinity 2722 PM6MD extllmod infinity infinity infinity infinity 9
10 Figure S1: A relationship between mutual RMSD of OPLS2005-optimized CAMVES conformers in the unpruned set (see Section 3.4 for details) and their final G DFT/COSMO-RS difference. The fraction of pairs with significant G DFT/COSMO-RS differences rises significantly above RMSD ~ 0.4 Å. 10
11 (a) (b) Figure S2: The number of total (a) and low-energy (b) unique conformers found in the unpruned MD/LLMOD sampling (see Section 3.4 for details) as a function of RMSD threshold used for pruning of conformers. An RMSD threshold of 0.4 Å was used for definitino of a unique conformer. The energy criterium of kj.mol -1 is used in all cases. The number of total unique conformers decreases approximately linearly as the threshold is loosened (a), however, a significant loss of the number of unique low-energy conformers only appears at thresholds above 0.6 Å. 11
12 Figure S3: Range of individual components ( X max X min, in kj.mol -1 ) of the overall free energy as defined in (1) equation in main text for two selected systems in the test set. The ranges of individual components sum up to total free energy, G DFT/COSMO-RS. The substantial decrease in the total range is due to correlation of individual components. 12
13 REFERENCES (1) Řezáč, J.; Bím, D.; Gutten, O.; Rulíšek, L.: Accuracy of Standard DFT Functionals and DLPNO-CCSD(T) Method in Conformational Sampling of Smaller Peptides and Medium- Sized Macrocycles: MP-CONF196 Data Set. Submitted. (2) Goerigk, L.; Kruse, H.; Grimme, S. Benchmarking density functional methods against the S66 and S66x8 datasets for non-covalent interactions. ChemPhysChem 2011, 12, (3) Gutten, O.; Beššeová, I.; Rulíšek, L.: Interaction of Metal Ions with Biomolecular Ligands: How Accurate Are Calculated Free Energies Associated with Metal Ion Complexation? J. Phys. Chem. A 2011, 115, (4) Gutten, O.; Rulíšek, L.: Predicting the Stability Constants of Metal-Ion Complexes from First Principles. Inorg. Chem. 2013, 52, (5) Grimme, S.: Supramolecular Binding Thermodynamics by Dispersion-Corrected Density Functional Theory. Chem. Eur. J. 2012, 18,
Intermolecular Forces in Density Functional Theory
Intermolecular Forces in Density Functional Theory Problems of DFT Peter Pulay at WATOC2005: There are 3 problems with DFT 1. Accuracy does not converge 2. Spin states of open shell systems often incorrect
More informationAdvanced in silico drug design
Advanced in silico drug design RNDr. Martin Lepšík, Ph.D. Lecture: Advanced scoring Palacky University, Olomouc 2016 1 Outline 1. Scoring Definition, Types 2. Physics-based Scoring: Master Equation Terms
More informationInteraction 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 informationIncrease in Solubility of Poorly-Ionizable Pharmaceuticals by Salt Formation. A Case of Agomelatine Sulfonates.
Supporting Information Increase in Solubility of Poorly-Ionizable Pharmaceuticals by Salt Formation. A Case of Agomelatine Sulfonates. Eliška Skořepová a,b, Daniel Bím e, Michal Hušák b, Jiří Klimeš d,
More informationPeptide folding in non-aqueous environments investigated with molecular dynamics simulations Soto Becerra, Patricia
University of Groningen Peptide folding in non-aqueous environments investigated with molecular dynamics simulations Soto Becerra, Patricia IMPORTANT NOTE: You are advised to consult the publisher's version
More informationSupporting Information: Predicting the Ionic Product of Water
Supporting Information: Predicting the Ionic Product of Water Eva Perlt 1,+, Michael von Domaros 1,+, Barbara Kirchner 1, Ralf Ludwig 2, and Frank Weinhold 3,* 1 Mulliken Center for Theoretical Chemistry,
More informationQuantum Chemistry for Supramolecular Complexes
Quantum Chemistry for Supramolecular Complexes Stefan Grimme Mulliken Center for Theoretical Chemistry, University of Bonn COSMO Symposium, March 2015 Supramolecular chemistry: chemistry beyond the covalent
More informationThis is a simple input file for the calculation of NMR chemical shieldings for a given molecule using the B3LYP functional and def2-tzvpp basis set:
Computing NMR parameters using ORCA This practical comes with a short lecture on the basics of the computation of NMR parameters using standard electronic structure theory methods. By now you should have
More informationTailoring 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 informationSUPPORTING INFORMATION. Table S1: Use of different functionals and variation of HF exchange on IS/HS splitting
SUPPORTING INFORMATION List of Contents Table S1: Use of different functionals and variation of HF exchange on IS/HS splitting S2 Table S2: Structural parameters using the B3LYP** functional for the IS
More informationLecture 2 and 3: Review of forces (ctd.) and elementary statistical mechanics. Contributions to protein stability
Lecture 2 and 3: Review of forces (ctd.) and elementary statistical mechanics. Contributions to protein stability Part I. Review of forces Covalent bonds Non-covalent Interactions: Van der Waals Interactions
More informationMolecular tailoring: a possible synthetic route to hexasilabenzene
SUPPORTING INFORMATION Molecular tailoring: a possible synthetic route to hexasilabenzene Zsolt Benedek, Tibor Szilvási and Tamás Veszprémi* Department of Inorganic and Analytical Chemistry, Budapest University
More informationFree-Energy Perturbation, Thermodynamic Integration, Potential of Mean Force (+working examples on chemical reactivity) Lubomír Rulíšek, Martin Srnec
Free-Energy Perturbation, Thermodynamic Integration, Potential of Mean Force (+working examples on chemical reactivity) Lubomír Rulíšek, Martin Srnec Institute of Organic Chemistry and Biochemistry AS
More informationBiochemistry,530:,, Introduc5on,to,Structural,Biology, Autumn,Quarter,2015,
Biochemistry,530:,, Introduc5on,to,Structural,Biology, Autumn,Quarter,2015, Course,Informa5on, BIOC%530% GraduateAlevel,discussion,of,the,structure,,func5on,,and,chemistry,of,proteins,and, nucleic,acids,,control,of,enzyma5c,reac5ons.,please,see,the,course,syllabus,and,
More informationRelative binding enthalpies from molecular dynamics simulations. using a direct method
Relative binding enthalpies from molecular dynamics simulations using a direct method Amitava Roy, 1 Duy P. Hua, 1 Joshua M. Ward, 1 and Carol Beth Post* Department of Medicinal Chemistry, Markey Center
More informationConformational 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 informationKeto-Enol Thermodynamics of Breslow Intermediates
Keto-Enol Thermodynamics of Breslow Intermediates Mathias Paul, Martin Breugst, Jörg-Martin Neudörfl, Raghavan B. Sunoj, and Albrecht Berkessel Computational Section of the Supporting Information 1 Computational
More informationSolutions 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 informationReviewers' comments: Reviewer #1 (Remarks to the Author):
Reviewers' comments: Reviewer #1 (Remarks to the Author): Review for Nature Communications manuscript NCOMMS-17-03154-T "Polyaromatic Molecular Peanuts" by Yoshizawa and co-workers reports the synthesis
More informationEur. J. Inorg. Chem WILEY-VCH Verlag GmbH & Co. KGaA, Weinheim, 2012 ISSN SUPPORTING INFORMATION
Eur. J. Inorg. Chem. 2012 WILEY-VCH Verlag GmbH & Co. KGaA, 69451 Weinheim, 2012 ISSN 1434 1948 SUPPORTING INFORMATION DOI: 10.1002/ejic.201200391 Title: Structural Insight into the Prolyl Hydroxylase
More informationSUPPLEMENTARY INFORMATION
Calculations predict a stable molecular crystal of N 8 : Barak Hirshberg a, R. Benny Gerber a,b, and Anna I. Krylov c a Institute of Chemistry and The Fritz Haber Center for Molecular Dynamics, The Hebrew
More informationUptake of OH radical to aqueous aerosol: a computational study
Uptake of OH radical to aqueous aerosol: a computational study Grigory Andreev Karpov Institute of Physical Chemistry 10 Vorontsovo pole, Moscow, 105064, Russia Institute of Physical Chemistry and Electrochemistry
More informationThe broad topic of physical metallurgy provides a basis that links the structure of materials with their properties, focusing primarily on metals.
Physical Metallurgy The broad topic of physical metallurgy provides a basis that links the structure of materials with their properties, focusing primarily on metals. Crystal Binding In our discussions
More informationSupplementary Information
Supplementary Information Supplementary Figure 1: Electronic Kohn-Sham potential profile of a charged monolayer MoTe 2 calculated using PBE-DFT. Plotted is the averaged electronic Kohn- Sham potential
More informationQUANTUM CHEMISTRY FOR TRANSITION METALS
QUANTUM CHEMISTRY FOR TRANSITION METALS Outline I Introduction II Correlation Static correlation effects MC methods DFT III Relativity Generalities From 4 to 1 components Effective core potential Outline
More informationCCSD(T) benchmarks of non-equilibrium water clusters: the importance of monomer deformation
CCSD(T) benchmarks of non-equilibrium water clusters: the importance of monomer deformation Biswajit Santra 1, Angelos Michaelides 1,2, and Matthias Scheffler 1 1 Fritz-Haber-Institut der MPG, Berlin,
More informationExcited States Calculations for Protonated PAHs
52 Chapter 3 Excited States Calculations for Protonated PAHs 3.1 Introduction Protonated PAHs are closed shell ions. Their electronic structure should therefore be similar to that of neutral PAHs, but
More informationAdvanced 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 informationRelative energetics of acetyl-histidine. protomers with and without Zn 2+ and a. benchmark of energy methods
Relative energetics of acetyl-histidine protomers with and without Zn 2+ and a benchmark of energy methods Markus Schneider and Carsten Baldauf,, Fritz-aber-Institut der Max-Planck-Gesellschaft, Theory
More informationMultivalent interactions in human biology
Cooperativity Multivalent interactions in human biology Multivalent interactions in supramolecular chemistry Additivity (?) Multivalent interactions in supramolecular chemistry In order to obtain a
More informationPaper No. 1: ORGANIC CHEMISTRY- I (Nature of Bonding and Stereochemistry)
Subject Chemistry Paper No and Title Paper 1: ORGANIC - I (Nature of Bonding Module No and Title Module Tag CHE_P1_M10 TABLE OF CONTENTS 1. Learning Outcomes 2. Introduction 3. Non-Covalent Interactions
More informationFrequency Estimation of Rare Events by Adaptive Thresholding
Frequency Estimation of Rare Events by Adaptive Thresholding J. R. M. Hosking IBM Research Division 2009 IBM Corporation Motivation IBM Research When managing IT systems, there is a need to identify transactions
More informationPerforming a Pharmacophore Search using CSD-CrossMiner
Table of Contents Introduction... 2 CSD-CrossMiner Terminology... 2 Overview of CSD-CrossMiner... 3 Searching with a Pharmacophore... 4 Performing a Pharmacophore Search using CSD-CrossMiner Version 2.0
More informationNucleotides 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 informationAb-initio modeling of opto-electronic properties of molecules in solvents and in proximity to a semiconductor nanoparticle
Ab-initio modeling of opto-electronic properties of molecules in solvents and in proximity to a semiconductor nanoparticle Alain Delgado (a,b), Stefano Corni (b), Carlo Andrea Rozzi (b) Stefano Pittalis
More informationSupporting Information
Supporting Information Wiley-VCH 2007 69451 Weinheim, Germany Synthesis and Properties of the THF Solvates of Extremely Soluble Bis(2,4,6-trimethylphenyl)calcium and Tris(2,6-dimethoxyphenyl)dicalcium
More informationMethods 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 informationConformational Searching using MacroModel and ConfGen. John Shelley Schrödinger Fellow
Conformational Searching using MacroModel and ConfGen John Shelley Schrödinger Fellow Overview Types of conformational searching applications MacroModel s conformation generation procedure General features
More informationType of file: PDF Title of file for HTML: Peer Review File Description:
Type of file: PDF Title of file for HTML: Supplementary Information Description: Supplementary Figures, Supplementary Table, Supplementary Notes and Supplementary References. Type of file: PDF Title of
More informationUsing Phase for Pharmacophore Modelling. 5th European Life Science Bootcamp March, 2017
Using Phase for Pharmacophore Modelling 5th European Life Science Bootcamp March, 2017 Phase: Our Pharmacohore generation tool Significant improvements to Phase methods in 2016 New highly interactive interface
More informationLecture 2-3: Review of forces (ctd.) and elementary statistical mechanics. Contributions to protein stability
Lecture 2-3: Review of forces (ctd.) and elementary statistical mechanics. Contributions to protein stability Part I. Review of forces Covalent bonds Non-covalent Interactions Van der Waals Interactions
More informationMOPAC Manual. ADF Modeling Suite
MOPAC Manual ADF Modeling Suite 2018 www.scm.com Sep 18, 2018 CONTENTS 1 Introduction 1 2 Input keywords 3 3 References 5 i ii CHAPTER ONE INTRODUCTION MOPAC [1 (page 5)] is a general-purpose semiempirical
More informationOcta-Acid Binding Affinities with QM Methods
Octa-Acid Binding Affinities with QM Methods Ulf Ryde Theoretical Chemistry Lund University Sweden 8 H B 7 Distance (Å) 6 5 4 3 2 1 0 0 1 2 3 4 5 6 Time (ns) 7 8 9 10 1. Intro 2. DFT-D3 3. CCSD(T) 4. MD
More informationMany-body electronic structure theory: Is Quantum Monte-Carlo affordable for molecular crystals?
Many-body electronic structure theory: Is Quantum Monte-Carlo affordable for molecular crystals? Gerit Brandenburg 23 rd of November 2017 CHEMISTRY DEPARTMENT, UNIVERSITY OF TURIN,
More informationElectronic structure theory: Fundamentals to frontiers. 2. Density functional theory
Electronic structure theory: Fundamentals to frontiers. 2. Density functional theory MARTIN HEAD-GORDON, Department of Chemistry, University of California, and Chemical Sciences Division, Lawrence Berkeley
More informationJoana Pereira Lamzin Group EMBL Hamburg, Germany. Small molecules How to identify and build them (with ARP/wARP)
Joana Pereira Lamzin Group EMBL Hamburg, Germany Small molecules How to identify and build them (with ARP/wARP) The task at hand To find ligand density and build it! Fitting a ligand We have: electron
More informationLecture 7: Numerical Tools
Lecture 7: Numerical Tools Fatih Guvenen January 10, 2016 Fatih Guvenen Lecture 7: Numerical Tools January 10, 2016 1 / 18 Overview Three Steps: V (k, z) =max c,k 0 apple u(c)+ Z c + k 0 =(1 + r)k + z
More informationScoring functions for of protein-ligand docking: New routes towards old goals
3nd Strasbourg Summer School on Chemoinformatics Strasbourg, June 25-29, 2012 Scoring functions for of protein-ligand docking: New routes towards old goals Christoph Sotriffer Institute of Pharmacy and
More informationAccurate and efficient modeling of inter- and intramolecular
pubs.acs.org/jpcl Accurate Modeling of Organic Molecular Crystals by Dispersion- Corrected Density Functional Tight Binding (DFTB) Jan Gerit Brandenburg and Stefan Grimme* Mulliken Center for Theoretical
More information6 Hydrophobic interactions
The Physics and Chemistry of Water 6 Hydrophobic interactions A non-polar molecule in water disrupts the H- bond structure by forcing some water molecules to give up their hydrogen bonds. As a result,
More informationAnion-π and π-π cooperative interactions
Chapter 5 Anion-π and π-π cooperative interactions 5.1 Introduction The design of selective receptors of anionic species is a very active area of research within supramolecular chemistry due to the potential
More informationProbing the Origins of Intermolecular Vibrational and Relaxational Dynamics in Organic Solids with CP2K
Probing the Origins of Intermolecular Vibrational and Relaxational Dynamics in Organic Solids with CP2K Michael Ruggiero Department of Chemical Engineering and Biotechnology, University of Cambridge CP2K
More informationElectronic, magnetic and spectroscopic properties of free Fe clusters
Electronic, magnetic and spectroscopic properties of free Fe clusters O. Šipr 1, M. Košuth 2, J. Minár 2, S. Polesya 2 and H. Ebert 2 1 Institute of Physics, Academy of Sciences of the Czech Republic,
More informationLecture 4: Band theory
Lecture 4: Band theory Very short introduction to modern computational solid state chemistry Band theory of solids Molecules vs. solids Band structures Analysis of chemical bonding in Reciprocal space
More informationCovalent Bonding. a. O b. Mg c. Ar d. C. a. K b. N c. Cl d. B
Covalent Bonding 1. Obtain the number of valence electrons for each of the following atoms from its group number and draw the correct Electron Dot Notation (a.k.a. Lewis Dot Structures). a. K b. N c. Cl
More informationMolecular 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 informationPhonon frequency (THz) R-3m: 60GPa. Supplementary Fig. 1. Phonon dispersion curves of R-3m-Ca 5 C 2 at 60 GPa.
Phonon frequency (THz) 2 15 R-3m: 6GPa 1 5 F Z Supplementary Fig. 1. Phonon dispersion curves of R-3m-Ca 5 C 2 at 6 GPa. 1 Phonon frequency (THz) 16 Pnma: 8 GPa 12 8 4 Z T Y S X U R Supplementary Fig.
More information1. Introduction to Clusters
1. Introduction to Clusters 1.1 The Field of Clusters Atomic clusters are aggregates of atoms containing from few to a few thousand atoms. Due to their small size, the properties of the clusters are, in
More informationSparse regularization of NIR spectra using implicit spiking and derivative spectroscopy
20th International Congress on Modelling and Simulation, Adelaide, Australia, 1 6 December 2013 www.mssanz.org.au/modsim2013 Sparse regularization of NIR spectra using implicit spiking and derivative spectroscopy
More informationAccurate and Efficient Model Energies for. Exploring Intermolecular Interactions in Molecular
Accurate and Efficient Model Energies for Exploring Intermolecular Interactions in Molecular Crystals Michael J. Turner, Simon Grabowsky, Dylan Jayatilaka and Mark A. Spackman* School of Chemistry & Biochemistry,
More informationProcess Performance and Quality
Chapter 5 Process Performance and Quality Evaluating Process Performance Identify opportunity 1 Define scope 2 Document process 3 Figure 5.1 Implement changes 6 Redesign process 5 Evaluate performance
More informationThe use of solvation models and the ONIOM layered approach in Gaussian.
The use of solvation models and the NIM layered approach in Gaussian. In this lab we will consider two techniques that are very useful to model larger systems: the use of solvation models to mimic systems
More informationOn the Chemical Free Energy of the Electrical Double Layer
1114 Langmuir 23, 19, 1114-112 On the Chemical Free Energy of the Electrical Double Layer Marian Manciu and Eli Ruckenstein* Department of Chemical Engineering, State University of New York at Buffalo,
More informationDispersion Interactions from the Exchange-Hole Dipole Moment
Dispersion Interactions from the Exchange-Hole Dipole Moment Erin R. Johnson and Alberto Otero-de-la-Roza Chemistry and Chemical Biology, University of California, Merced E. R. Johnson (UC Merced) Dispersion
More informationThe Schrödinger equation for many-electron systems
The Schrödinger equation for many-electron systems Ĥ!( x,, x ) = E!( x,, x ) 1 N 1 1 Z 1 Ĥ = " $ # " $ + $ 2 r 2 A j j A, j RAj i, j < i a linear differential equation in 4N variables (atomic units) (3
More informationTitle Theory of solutions in the energy r of the molecular flexibility Author(s) Matubayasi, N; Nakahara, M Citation JOURNAL OF CHEMICAL PHYSICS (2003), 9702 Issue Date 2003-11-08 URL http://hdl.handle.net/2433/50354
More informationIntermolecular Forces & Condensed Phases
Intermolecular Forces & Condensed Phases CHEM 107 T. Hughbanks READING We will discuss some of Chapter 5 that we skipped earlier (Van der Waals equation, pp. 145-8), but this is just a segue into intermolecular
More informationINTRODUCTION TO GIS. Dr. Ori Gudes
INTRODUCTION TO GIS Dr. Ori Gudes Outline of the Presentation What is GIS? What s the rational for using GIS, and how GIS can be used to solve problems? Explore a GIS map and get information about map
More informationAtom-molecule molecule collisions in spin-polarized polarized alkalis: potential energy surfaces and quantum dynamics
Atom-molecule molecule collisions in spin-polarized polarized alkalis: potential energy surfaces and quantum dynamics Pavel Soldán, Marko T. Cvitaš and Jeremy M. Hutson University of Durham with Jean-Michel
More informationUnified approach to the classical statistical analysis of small signals
PHYSICAL REVIEW D VOLUME 57, NUMBER 7 1 APRIL 1998 Unified approach to the classical statistical analysis of small signals Gary J. Feldman * Department of Physics, Harvard University, Cambridge, Massachusetts
More informationWalter Kohn was awarded with the Nobel Prize in Chemistry in 1998 for his development of the density functional theory.
Walter Kohn was awarded with the Nobel Prize in Chemistry in 1998 for his development of the density functional theory. Walter Kohn receiving his Nobel Prize from His Majesty the King at the Stockholm
More informationLecture 11: Potential Energy Functions
Lecture 11: Potential Energy Functions Dr. Ronald M. Levy ronlevy@temple.edu Originally contributed by Lauren Wickstrom (2011) Microscopic/Macroscopic Connection The connection between microscopic interactions
More informationHigh-level Quantum Chemistry Methods and Benchmark Datasets for Molecules
High-level Quantum Chemistry Methods and Benchmark Datasets for Molecules Markus Schneider Fritz Haber Institute of the MPS, Berlin, Germany École Polytechnique Fédérale de Lausanne, Switzerland دانشگاه
More informationResearch Article Polymorphs of Tolfenamic Acids: Stability Analysis Using Cluster Method
Physics Research International Volume 2016, Article ID 3537842, 5 pages http://dx.doi.org/10.1155/2016/3537842 Research Article Polymorphs of Tolfenamic Acids: Stability Analysis Using Cluster Method Lee
More informationRF System Calibration Using Beam Orbits at LEP
EUROPEAN ORGANIZATION FOR NUCLEAR RESEARCH CERN SL DIVISION CERN-SL-22-28 OP LEP Energy Working Group 2/1 RF System Calibration Using Beam Orbits at LEP J. Wenninger Abstract The target for beam energy
More informationBinding free energies in the SAMPL5 octa-acid host guest challenge calculated with DFT-D3 and CCSD(T)
J Comput Aided Mol Des (2017) 31:87 106 DOI 10.1007/s10822-016-9957-5 Binding free energies in the SAMPL5 octa-acid host guest challenge calculated with DFT-D3 and CCSD(T) Octav Caldararu 1 Martin A. Olsson
More informationAccounting for Solvation in Quantum Chemistry. Comp chem spring school 2014 CSC, Finland
Accounting for Solvation in Quantum Chemistry Comp chem spring school 2014 CSC, Finland In solution or in a vacuum? Solvent description is important when: Polar solvent: electrostatic stabilization Solvent
More informationFormation of complexes: thermodynamics
LN08-1 Formation of complexes: thermodynamics Or, How to think about making and using transition metal complexes Its important to always consider, and distinguish between, stability and inertness of complexes.
More informationSearch for an electric dipole moment with polarized beams in storage rings
Search for an electric dipole moment with polarized beams in storage rings Università di Ferrara and INFN 44 Ferrara, Italy E-mail: lenisa@fe.infn.it E. J. Stephenson Center for Exploration of Energy and
More informationNoncovalent 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 informationplatinum and iridium complexes
SUPPORTING INFORMATION A relativistic DFT methodology for calculating the structures and NMR chemical shifts of octahedral platinum and iridium complexes Jan Vícha, Michael Patzschke, Radek Marek Figure
More informationMachine learning the Born-Oppenheimer potential energy surface: from molecules to materials. Gábor Csányi Engineering Laboratory
Machine learning the Born-Oppenheimer potential energy surface: from molecules to materials Gábor Csányi Engineering Laboratory Interatomic potentials for molecular dynamics Transferability biomolecular
More informationS66: 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 informationAccurate van der Waals interactions from ground state electron density
Accurate van der Waals interactions from ground state electron density Alexandre Tkatchenko Theory Department, Fritz Haber Institut der MPG Berlin, Germany tkatchen@fhi berlin.mpg.de Haber Institute EXCITCM09,
More informationCHAPTER 9: LIQUIDS AND SOLIDS
CHAPTER 9: LIQUIDS AND SOLIDS Section 9.1 Liquid/Vapor Equilibrium Vaporization process in which a liquid vapor open container - evaporation continues until all liquid evaporates closed container 1) Liquid
More informationLecture 17 - The Secrets we have Swept Under the Rug
1.0 0.5 0.0-0.5-0.5 0.0 0.5 1.0 Lecture 17 - The Secrets we have Swept Under the Rug A Puzzle... What makes 3D Special? Example (1D charge distribution) A stick with uniform charge density λ lies between
More informationBonding in solids The interaction of electrons in neighboring atoms of a solid serves the very important function of holding the crystal together.
Bonding in solids The interaction of electrons in neighboring atoms of a solid serves the very important function of holding the crystal together. For example Nacl In the Nacl lattice, each Na atom is
More informationConcepts of Chemical Bonding and Molecular Geometry Part 1: Ionic and Covalent Bonds. David A. Katz Pima Community College Tucson, AZ
Concepts of Chemical Bonding and Molecular Geometry Part 1: Ionic and Covalent Bonds David A. Katz Pima Community College Tucson, AZ Chemical Bonds Three basic types of bonds: Ionic Electrostatic attraction
More informationWater complexes in the keto-enol tautomeric equilibrium
Universidad Autónoma Metropolitana Water complexes in the keto-enol tautomeric equilibrium Nelly González-Rivas, Mariano Méndez-Chávez and Andrés Cedillo Departamento de Química, UAM-Iztapalapa, San Rafael
More informationBonding. Computer Lab: Ionic Bonds. Important Notes 3/22/18
Bonding What are ionic bonds, and how are they formed? Computer Lab: Ionic Bonds Go to http://www.pbslearningmedia.org/asset/ lsps07_int_ionicbonding/ Read each screen and follow the directions where appropriate.
More informationCreating a Pharmacophore Query from a Reference Molecule & Scaffold Hopping in CSD-CrossMiner
Table of Contents Creating a Pharmacophore Query from a Reference Molecule & Scaffold Hopping in CSD-CrossMiner Introduction... 2 CSD-CrossMiner Terminology... 2 Overview of CSD-CrossMiner... 3 Features
More informationSensitive NMR Approach for Determining the Binding Mode of Tightly Binding Ligand Molecules to Protein Targets
Supporting information Sensitive NMR Approach for Determining the Binding Mode of Tightly Binding Ligand Molecules to Protein Targets Wan-Na Chen, Christoph Nitsche, Kala Bharath Pilla, Bim Graham, Thomas
More informationChem 673, Problem Sets 4 & 5 Due Tuesday, December 3, Problems from Carter: Chapter 6: 6.1, 6.3, 6.7, 6.9 Chapter 7: 7.2a,b,e,g,i,j, 7.
Chem 673, Problem Sets 4 & 5 Due Tuesday, December 3, 2013 Problems from Carter: Chapter 6: 6.1, 6.3, 6.7, 6.9 Chapter 7: 7.2a,b,e,g,i,j, 7.6, (1) Use the angular overlap table given on the back page of
More informationWhy Is Molecular Interaction Important in Our Life
Why Is Molecular Interaction Important in ur Life QuLiS and Graduate School of Science iroshima University http://www.nabit.hiroshima-u.ac.jp/iwatasue/indexe.htm Suehiro Iwata Sept. 29, 2007 Department
More informationElectronic 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 informationMultiwalled nanotube faceting unravelled
Multiwalled nanotube faceting unravelled Itai Leven, Roberto Guerra, Andrea Vanossi, Erio Tosatti, and Oded Hod NATURE NANOTECHNOLOGY www.nature.com/naturenanotechnology 1 1. Achiral DWCNTs optimized using
More informationTowards Physics-based Models for ADME/Tox. Tyler Day
Towards Physics-based Models for ADME/Tox Tyler Day Overview Motivation Application: P450 Site of Metabolism Application: Membrane Permeability Future Directions and Applications Motivation Advantages
More informationINSPIRE Monitoring and Reporting Implementing Rule Draft v2.1
INSPIRE Infrastructure for Spatial Information in Europe INSPIRE Monitoring and Reporting Implementing Rule Draft v2.1 Title INSPIRE Monitoring and Reporting Implementing Rule v2.1 Creator DT Monitoring
More informationwith the larger dimerization energy also exhibits the larger structural changes.
A7. Looking at the image and table provided below, it is apparent that the monomer and dimer are structurally almost identical. Although angular and dihedral data were not included, these data are also
More informationJulian Hassinger Madeleine Jensen
UNIVERSITY OF CALIFORNIA, BERKELEY CHEM C130/MCB C100A MIDTERM EXAMINATION #2 October 22, 2014 INSTRUCTORS: John Kuriyan and David Savage THE TIME LIMIT FOR THIS EXAMINATION: 1 HOUR 50 MINUTES SIGNATURE:
More information