1
|
Das G, Ibrahim FA, Khalil ZA, Bazin P, Chandra F, AbdulHalim RG, Prakasam T, Das AK, Sharma SK, Varghese S, Kirmizialtin S, Jagannathan R, Saleh N, Benyettou F, Roz ME, Addicoat M, Olson MA, Rao DSS, Prasad SK, Trabolsi A. Ionic Covalent Organic Framework as a Dual Functional Sensor for Temperature and Humidity. SMALL (WEINHEIM AN DER BERGSTRASSE, GERMANY) 2024; 20:e2311064. [PMID: 38396219 DOI: 10.1002/smll.202311064] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Received: 11/30/2023] [Revised: 02/07/2024] [Indexed: 02/25/2024]
Abstract
Visual sensing of humidity and temperature by solids plays an important role in the everyday life and in industrial processes. Due to their hydrophobic nature, most covalent organic framework (COF) sensors often exhibit poor optical response when exposed to moisture. To overcome this challenge, the optical response is set out to improve, to moisture by incorporating H-bonding ionic functionalities into the COF network. A highly sensitive COF, consisting of guanidinium and diformylpyridine linkers (TG-DFP), capable of detecting changes in temperature and moisture content is fabricated. The hydrophilic nature of the framework enables enhanced water uptake, allowing the trapped water molecules to form a large number of hydrogen bonds. Despite the presence of non-emissive building blocks, the H-bonds restrict internal bond rotation within the COF, leading to reversible fluorescence and solid-state optical hydrochromism in response to relative humidity and temperature.
Collapse
Affiliation(s)
- Gobinda Das
- Chemistry Program, New York University Abu Dhabi (NYUAD), Saadiyat Island, Abu Dhabi, 129188, United Arab Emirates
| | - Fayrouz Abou Ibrahim
- Chemistry Program, New York University Abu Dhabi (NYUAD), Saadiyat Island, Abu Dhabi, 129188, United Arab Emirates
| | - Zahraa Abou Khalil
- Laboratoire Catalyse et Spectrochimie, CNRS, Ensicaen, Université de Caen, 6, Boulevard Maréchal Juin 14050, Caen, France
| | - Philippe Bazin
- Laboratoire Catalyse et Spectrochimie, CNRS, Ensicaen, Université de Caen, 6, Boulevard Maréchal Juin 14050, Caen, France
| | - Falguni Chandra
- Chemistry Department, College of Science, United Arab Emirates University, P.O. Box 15551, Al-Ain, United Arab Emirates
| | - Rasha G AbdulHalim
- Chemistry Program, New York University Abu Dhabi (NYUAD), Saadiyat Island, Abu Dhabi, 129188, United Arab Emirates
| | - Thirumurugan Prakasam
- Chemistry Program, New York University Abu Dhabi (NYUAD), Saadiyat Island, Abu Dhabi, 129188, United Arab Emirates
| | - Akshaya Kumar Das
- Chemistry Program, New York University Abu Dhabi (NYUAD), Saadiyat Island, Abu Dhabi, 129188, United Arab Emirates
| | - Sudhir Kumar Sharma
- Engineering Division, New York University Abu Dhabi (NYUAD), Abu Dhabi, 129188, United Arab Emirates
| | - Sabu Varghese
- New York University Abu Dhabi, Abu Dhabi, 129188, United Arab Emirates
| | - Serdal Kirmizialtin
- Chemistry Program, New York University Abu Dhabi (NYUAD), Saadiyat Island, Abu Dhabi, 129188, United Arab Emirates
| | - Ramesh Jagannathan
- Engineering Division, New York University Abu Dhabi (NYUAD), Abu Dhabi, 129188, United Arab Emirates
| | - Na'il Saleh
- Chemistry Department, College of Science, United Arab Emirates University, P.O. Box 15551, Al-Ain, United Arab Emirates
- National Water and Energy center, United Arab Emirates University, P.O. Box 15551, Al Ain, United Arab Emirates
| | - Farah Benyettou
- Chemistry Program, New York University Abu Dhabi (NYUAD), Saadiyat Island, Abu Dhabi, 129188, United Arab Emirates
| | - Mohamad El Roz
- Laboratoire Catalyse et Spectrochimie, CNRS, Ensicaen, Université de Caen, 6, Boulevard Maréchal Juin 14050, Caen, France
| | - Matthew Addicoat
- School of Science and Technology, Nottingham Trent University, Clifton Lane, NG11 8NS, Nottingham, NG118NS, UK
| | - Mark A Olson
- Department of Physical and Environmental Sciences, Texas A&M University Corpus Christi, 6300 Ocean Dr, Corpus Christi, TX, 78412, USA
| | - D S Shankar Rao
- Centre for Nano and Soft Matter Sciences(CeNS), Arkavathi, Survey No.7, Shivanapura, Dasanapura Hobli, Bengaluru, 562162, India
| | - S Krishna Prasad
- Centre for Nano and Soft Matter Sciences(CeNS), Arkavathi, Survey No.7, Shivanapura, Dasanapura Hobli, Bengaluru, 562162, India
| | - Ali Trabolsi
- Chemistry Program, New York University Abu Dhabi (NYUAD), Saadiyat Island, Abu Dhabi, 129188, United Arab Emirates
- NYUAD Water Research Center, New York University Abu Dhabi (NYUAD), Saadiyat Island, Abu Dhabi, 129188, United Arab Emirates
| |
Collapse
|
2
|
Junkermeier CE, Kobebel J, Lavarez K, Adra RM, Yang J, Diaz VA, Paupitz R, Psofogiannaki G. Density functional tight-binding derived data of gas capture in functionalized carbophenes. Data Brief 2024; 55:110652. [PMID: 39035841 PMCID: PMC11259943 DOI: 10.1016/j.dib.2024.110652] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/22/2023] [Revised: 05/28/2024] [Accepted: 06/11/2024] [Indexed: 07/23/2024] Open
Abstract
The presented data relates to the investigation of the adsorption properties of carbon dioxide (CO2), methane (CH4), and dihydrogen (H2) within pristine and functionalized carbophene pores. The carbophenes were functionalized with one of the groups carboxyl (COOH), amine (NH2), nitro (NO2), hydroxyl (OH), or an amide (CONH2, NHCOOH, and N(COOH)2) groups. The systems were optimized using the density functional tight-binding theory code DFTB+ (pre-compiled Version 19.1) with the matsci Slater-Koster files on the Mana high performance computing cluster at the University of Hawai'i at Mānoa. The dataset consists of the molecular geometries, lattice vectors, and the total energies for each specific system. One possible use of the data is for training or validating force fields for running molecular dynamics simulations.
Collapse
Affiliation(s)
- Chad E. Junkermeier
- Department of Physics and Astronomy, University of Hawai`i at Mānoa, Honolulu HI 96822, USA
- Science, Technology, Engineering, and Mathematics Department, University of Hawai`i Maui College, Kahului HI 96732, USA
| | - Jedediah Kobebel
- Science, Technology, Engineering, and Mathematics Department, University of Hawai`i Maui College, Kahului HI 96732, USA
| | - Kat Lavarez
- Department of Physics and Astronomy, University of Hawai`i at Mānoa, Honolulu HI 96822, USA
| | - R. Martin Adra
- Department of Physics and Astronomy, University of Hawai`i at Mānoa, Honolulu HI 96822, USA
| | - Jirui Yang
- Department of Physics and Astronomy, University of Hawai`i at Mānoa, Honolulu HI 96822, USA
| | - Valeria Aparicio Diaz
- Department of Physics and Astronomy, University of Hawai`i at Mānoa, Honolulu HI 96822, USA
| | - Ricardo Paupitz
- Departamento de Física, Universidade Estadual Paulista, UNESP, 13506-900, Rio Claro, SP, Brazil
| | - George Psofogiannaki
- Department of Chemical and Biological Engineering, University of Ottawa, Ottawa, Ontario, Canada
| |
Collapse
|
3
|
Fox R, Klug J, Thompson D, Reilly A. Computational predictions of cocrystal formation: A benchmark study of 28 assemblies comparing five methods from high-throughput to advanced models. J Comput Chem 2024. [PMID: 38958249 DOI: 10.1002/jcc.27454] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/18/2024] [Revised: 06/10/2024] [Accepted: 06/12/2024] [Indexed: 07/04/2024]
Abstract
Cocrystals are assemblies of more than one type of molecule stabilized through noncovalent interactions. They are promising materials for improved drug formulation in which the stability, solubility, or biocompatibility of the active pharmaceutical ingredient (API) is improved by including a coformer. In this work, a range of density functional theory (DFT) and density functional tight binding (DFTB) models are systematically compared for their ability to predict the lattice enthalpy of a broad range of existing pharmaceutically relevant cocrystals. These range from cocrystals containing model compounds 4,4'-bipyridine and oxalic acid to those with the well benchmarked APIs of aspirin and paracetamol, all tested with a large set of alternative coformers. For simple cocrystals, there is a general consensus in lattice enthalpy calculated by the different DFT models. For the cocrystals with API coformers the cocrystals, enthalpy predictions depend strongly on the DFT model. The significantly lighter DFTB models predict unrealistic values of lattice enthalpy even for simple cocrystals.
Collapse
Affiliation(s)
- Robert Fox
- School of Chemical Sciences, Dublin City University, Dublin, Ireland
| | - Joaquin Klug
- Department of Life Sciences, Faculty of Sciences, Atlantic Technological University, ATU Sligo, Sligo, Ireland
| | - Damien Thompson
- Department of Physics, Bernal Institute, University of Limerick, Limerick, Ireland
| | - Anthony Reilly
- School of Chemical Sciences, Dublin City University, Dublin, Ireland
| |
Collapse
|
4
|
Mun H, Lorpaiboon W, Ho J. In Search of the Best Low-Cost Methods for Efficient Screening of Conformers. J Phys Chem A 2024; 128:4391-4400. [PMID: 38754085 DOI: 10.1021/acs.jpca.4c01407] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 05/18/2024]
Abstract
Locating the lowest energy conformer is crucial for the accurate computation of equilibrium properties of molecular systems. This paper examines the performance of efficient low-cost methods in terms of the alignment and relative energies of their energy minima against the benchmark revDSD-PBEP86-D4/def2-TZVPP//MP2/cc-pVTZ potential energy surface. The low-cost methods considered include GFN-FF, GFN2-xTB, DFTB3, HF-3c, B97-3c, PBEh-3c, and r2SCAN-3c composite methods against a diverse test set of 20 compounds including alkanes, perfluoroalkyl molecules, peptides, open-shell radicals, and Zn(II) complexes of varying sizes. The "3c" composite methods are generally more accurate, but are at least 2-3 orders of magnitude more expensive than tight-binding methods which have energy minima that align well with the benchmark potential energy surface. The findings of this paper were further exploited to introduce a simple strategy involving Grimme's CENSO energy-sorting algorithm that resulted in up to an order of magnitude reduction in computational time for locating the lowest energy conformer on the revDSD-PBEP86-D4/def2-TZVPP//MP2/cc-pVTZ surface.
Collapse
Affiliation(s)
- Haedam Mun
- School of Chemistry, The University of New South Wales, Sydney, NSW 2052, Australia
| | - Wanutcha Lorpaiboon
- School of Chemistry, The University of New South Wales, Sydney, NSW 2052, Australia
| | - Junming Ho
- School of Chemistry, The University of New South Wales, Sydney, NSW 2052, Australia
| |
Collapse
|
5
|
Sahu R, Yamijala SSRKC, Rao KV, Reddy SK. Dispersion-Driven Cooperativity in Alkyl Perylene Diimide Oligomers: Insights from Density Functional Theory. Chemphyschem 2024:e202400235. [PMID: 38807431 DOI: 10.1002/cphc.202400235] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/03/2024] [Revised: 05/16/2024] [Accepted: 05/28/2024] [Indexed: 05/30/2024]
Abstract
The cooperative mechanism is of paramount importance in the synthesis of supramolecular polymers with desired characteristics, including molecular mass, polydispersity, and morphology. It is primarily driven by the presence of intermolecular interactions, which encompass strong hydrogen bonding, metal-ligand interactions, and dipole-dipole interactions. In this study, we utilize density functional theory and energy decomposition analysis to investigate the cooperative behavior of perylene diimide (PDI) oligomers with alkyl chains at their imide positions, which lack the previously mentioned interactions. Our systematic examination reveals that dispersion interactions originating from the alkyl side-chain substituents play an important role in promoting cooperativity within these PDIs. This influence becomes even more pronounced for alkyl chain lengths beyond hexyl groups. The energy decomposition analysis reveals that the delicate balance between dispersion energy and Pauli repulsion energy is the key driver of cooperative behavior in PDIs. Additionally, we have developed a mathematical model capable of predicting the saturated binding energies for PDI oligomers of varying sizes and alkyl chain lengths. Overall, our findings emphasize the previously undervalued significance of dispersion forces in cooperative supramolecular polymerization, enhancing our overall understanding of the cooperative mechanism.
Collapse
Affiliation(s)
- Rahul Sahu
- Centre for Computational and Data Sciences, Indian Institute of Technology Kharagpur, Kharagpur, West Bengal, Pin, 721302, India
| | - Sharma S R K C Yamijala
- Department of Chemistry, Centre for Atomistic Modelling and Materials Design, Centre for Quantum Information, Communication, and Computing, Centre for Molecular Materials and Functions, Indian Institute of Technology Madras, Chennai, Tamil Nadu, Pin, 600036, India
- Centre for Atomistic Modelling and Materials Design, Indian Institute of Technology Madras, Chennai, Tamil Nadu, Pin, 600036, India
- Centre for Quantum Information, Communication, and Computing, Indian Institute of Technology Madras, Chennai, Tamil Nadu, Pin, 600036, India
- Centre for Molecular Materials and Functions, Indian Institute of Technology Madras, Chennai, Tamil Nadu, Pin, 600036, India
| | - Kotagiri Venkata Rao
- Department of Chemistry, Indian Institute of Technology Hyderabad, Kandi, Sangareddy, Telangana, Pin, 502285, India
| | - Sandeep K Reddy
- Centre for Computational and Data Sciences, Indian Institute of Technology Kharagpur, Kharagpur, West Bengal, Pin, 721302, India
| |
Collapse
|
6
|
Higuchi Y, Saleh MA, Anada T, Tanaka M, Hishida M. Rotational Dynamics of Water near Osmolytes by Molecular Dynamics Simulations. J Phys Chem B 2024; 128:5008-5017. [PMID: 38728154 DOI: 10.1021/acs.jpcb.3c08470] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 05/12/2024]
Abstract
The behavior of water molecules around organic molecules has attracted considerable attention as a crucial factor influencing the properties and functions of soft matter and biomolecules. Recently, it has been suggested that the change in protein stability upon the addition of small organic molecules (osmolytes) is dominated by the change in the water dynamics caused by the osmolyte, where the dynamics of not only the directly interacting water molecules but also the long-range hydration layer affect the protein stability. However, the relation between the long-range structure of hydration water in various solutions and the water dynamics remains unclear at the molecular level. We performed density-functional tight-binding molecular dynamics simulations to elucidate the varying rotational dynamics of water molecules in 15 osmolyte solutions. A positive correlation was observed between the rotational relaxation time and our proposed normalized parameter obtained by dividing the number of hydrogen bonds between water molecules by the number of nearest-neighbor water molecules. For the 15 osmolyte solutions, an increase or a decrease in the value of the normalized parameter for the second hydration shell tended to result in water molecules with slow and fast rotational dynamics, respectively, thus illustrating the importance of the second hydration shell for the rotational dynamics of water molecules. Our simulation results are anticipated to advance the current understanding of water dynamics around organic molecules and the long-range structure of water molecules.
Collapse
Affiliation(s)
- Yuji Higuchi
- Research Institute for Information Technology, Kyushu University, 744 Motooka, Nishi-ku, Fukuoka 819-0395, Japan
| | - Md Abu Saleh
- Department of Applied Chemistry, Graduate School of Engineering, Kyushu University, Fukuoka 819-0395, Japan
| | - Takahisa Anada
- Department of Applied Chemistry, Graduate School of Engineering, Kyushu University, Fukuoka 819-0395, Japan
- Institute for Materials Chemistry and Engineering, Kyushu University, Fukuoka 819-0395, Japan
| | - Masaru Tanaka
- Department of Applied Chemistry, Graduate School of Engineering, Kyushu University, Fukuoka 819-0395, Japan
- Institute for Materials Chemistry and Engineering, Kyushu University, Fukuoka 819-0395, Japan
| | - Mafumi Hishida
- Department of Chemistry, Faculty of Science, Tokyo University of Science, 1-3 Kagurazaka, Shinjuku, Tokyo 162-8601, Japan
| |
Collapse
|
7
|
Díaz Mirón G, Lien-Medrano CR, Banerjee D, Morzan UN, Sentef MA, Gebauer R, Hassanali A. Exploring the Mechanisms behind Non-aromatic Fluorescence with the Density Functional Tight Binding Method. J Chem Theory Comput 2024; 20:3864-3878. [PMID: 38634760 DOI: 10.1021/acs.jctc.4c00125] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 04/19/2024]
Abstract
Recent experimental findings reveal nonconventional fluorescence emission in biological systems devoid of conjugated bonds or aromatic compounds, termed non-aromatic fluorescence (NAF). This phenomenon is exclusive to aggregated or solid states and remains absent in monomeric solutions. Previous studies focused on small model systems in vacuum show that the carbonyl stretching mode along with strong interaction of short hydrogen bonds (SHBs) remains the primary vibrational mode explaining NAF in these systems. In order to simulate larger model systems taking into account the effects of the surrounding environment, in this work we propose using the density functional tight-binding (DFTB) method in combination with non-adiabatic molecular dynamics (NAMD) and the mixed quantum/molecular mechanics (QM/MM) approach. We investigate the mechanism behind NAF in the crystal structure of l-pyroglutamine-ammonium, comparing it with the related nonfluorescent amino acid l-glutamine. Our results extend our previous findings to more realistic systems, demonstrating the efficiency and robustness of the proposed DFTB method in the context of NAMD in biological systems. Furthermore, due to its inherent low computational cost, this method allows for a better sampling of the nonradiative events at the conical intersection which is crucial for a complete understanding of this phenomenon. Beyond contributing to the ongoing exploration of NAF, this work paves the way for future application of this method in more complex biological systems such as amyloid aggregates, biomaterials, and non-aromatic proteins.
Collapse
Affiliation(s)
- Gonzalo Díaz Mirón
- Condensed Matter and Statistical Physics, The Abdus Salam International Centre for Theoretical Physics, 34151 Trieste, Italy
| | - Carlos R Lien-Medrano
- Institute for Theoretical Physics and Bremen Center for Computational Materials Science, University of Bremen, 28359 Bremen, Germany
| | - Debarshi Banerjee
- Condensed Matter and Statistical Physics, The Abdus Salam International Centre for Theoretical Physics, 34151 Trieste, Italy
- Scuola Internazionale Superiore di Studi Avanzati (SISSA), 34136 Trieste, Italy
| | - Uriel N Morzan
- Instituto de Fisica de Buenos Aires, Facultad de Ciencias Exactas y Naturales, Universidad de Buenos Aires, C1428EGA Buenos Aires, Argentina
| | - Michael A Sentef
- Institute for Theoretical Physics and Bremen Center for Computational Materials Science, University of Bremen, 28359 Bremen, Germany
- Center for Free-Electron Laser Science (CFEL), Max Planck Institute for the Structure and Dynamics of Matter, 22761 Hamburg, Germany
| | - Ralph Gebauer
- Condensed Matter and Statistical Physics, The Abdus Salam International Centre for Theoretical Physics, 34151 Trieste, Italy
| | - Ali Hassanali
- Condensed Matter and Statistical Physics, The Abdus Salam International Centre for Theoretical Physics, 34151 Trieste, Italy
| |
Collapse
|
8
|
Wu X, Hartmann P, Berne D, De Bruyn M, Cuminet F, Wang Z, Zechner JM, Boese AD, Placet V, Caillol S, Barta K. Closed-loop recyclability of a biomass-derived epoxy-amine thermoset by methanolysis. Science 2024; 384:eadj9989. [PMID: 38603486 DOI: 10.1126/science.adj9989] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/27/2023] [Accepted: 02/16/2024] [Indexed: 04/13/2024]
Abstract
Epoxy resin thermosets (ERTs) are an important class of polymeric materials. However, owing to their highly cross-linked nature, they suffer from poor recyclability, which contributes to an unacceptable level of environmental pollution. There is a clear need for the design of inherently recyclable ERTs that are based on renewable resources. We present the synthesis and closed-loop recycling of a fully lignocellulose-derivable epoxy resin (DGF/MBCA), prepared from dimethyl ester of 2,5-furandicarboxylic acid (DMFD), 4,4'-methylenebis(cyclohexylamine) (MBCA), and glycidol, which displays excellent thermomechanical properties (a glass transition temperature of 170°C, and a storage modulus at 25°C of 1.2 gigapascals). Notably, the material undergoes methanolysis in the absence of any catalyst, regenerating 90% of the original DMFD. The diamine MBCA and glycidol can subsequently be reformed by acetolysis. Application and recycling of DGF/MBCA in glass and plant fiber composites are demonstrated.
Collapse
Affiliation(s)
- Xianyuan Wu
- Stratingh Institute for Chemistry, University of Groningen, 9747AG Groningen, Netherlands
- Institute of Chemistry, Organic and Bioorganic Chemistry, University of Graz, 8010 Graz, Austria
| | - Peter Hartmann
- Institute of Chemistry, Organic and Bioorganic Chemistry, University of Graz, 8010 Graz, Austria
| | - Dimitri Berne
- ICGM, Univ Montpellier, CNRS, ENSCM, 34000 Montpellier, France
| | - Mario De Bruyn
- Institute of Chemistry, Organic and Bioorganic Chemistry, University of Graz, 8010 Graz, Austria
| | - Florian Cuminet
- ICGM, Univ Montpellier, CNRS, ENSCM, 34000 Montpellier, France
| | - Zhiwen Wang
- Institute of Chemistry, Organic and Bioorganic Chemistry, University of Graz, 8010 Graz, Austria
| | | | - Adrian Daniel Boese
- Institute of Chemistry, Organic and Bioorganic Chemistry, University of Graz, 8010 Graz, Austria
| | - Vincent Placet
- Université de Franche-Comté, CNRS, institut FEMTO-ST, 2500 Besançon, France
| | - Sylvain Caillol
- ICGM, Univ Montpellier, CNRS, ENSCM, 34000 Montpellier, France
| | - Katalin Barta
- Institute of Chemistry, Organic and Bioorganic Chemistry, University of Graz, 8010 Graz, Austria
| |
Collapse
|
9
|
Manikandan M, Nicolini P, Hapala P. Computational Design of Photosensitive Polymer Templates To Drive Molecular Nanofabrication. ACS NANO 2024; 18:9969-9979. [PMID: 38545921 PMCID: PMC11008366 DOI: 10.1021/acsnano.3c10575] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Received: 10/27/2023] [Revised: 03/15/2024] [Accepted: 03/21/2024] [Indexed: 04/10/2024]
Abstract
Molecular electronics promises the ultimate level of miniaturization of computers and other machines as organic molecules are the smallest known physical objects with nontrivial structure and function. But despite the plethora of molecular switches, memories, and motors developed during the almost 50-years long history of molecular electronics, mass production of molecular computers is still an elusive goal. This is mostly due to the lack of scalable nanofabrication methods capable of rapidly producing complex structures (similar to silicon chips or living cells) with atomic precision and a small number of defects. Living nature solves this problem by using linear polymer templates encoding large volumes of structural information into sequence of hydrogen bonded end groups which can be efficiently replicated and which can drive assembly of other molecular components into complex supramolecular structures. In this paper, we propose a nanofabrication method based on a class of photosensitive polymers inspired by these natural principles, which can operate in concert with UV photolithography used for fabrication of current microelectronic processors. We believe that such a method will enable a smooth transition from silicon toward molecular nanoelectronics and photonics. To demonstrate its feasibility, we performed a computational screening of candidate molecules that can selectively bind and therefore allow the deterministic assembly of molecular components. In the process, we unearthed trends and design principles applicable beyond the immediate scope of our proposed nanofabrication method, e.g., to biologically relevant DNA analogues and molecular recognition within hydrogen-bonded systems.
Collapse
Affiliation(s)
- Mithun Manikandan
- Institute of Physics (FZU), Czech
Academy of Sciences, Na Slovance 2, 182 00 Prague, Czech Republic
| | - Paolo Nicolini
- Institute of Physics (FZU), Czech
Academy of Sciences, Na Slovance 2, 182 00 Prague, Czech Republic
| | - Prokop Hapala
- Institute of Physics (FZU), Czech
Academy of Sciences, Na Slovance 2, 182 00 Prague, Czech Republic
| |
Collapse
|
10
|
Cordova M, Moutzouri P, Nilsson Lill SO, Cousen A, Kearns M, Norberg ST, Svensk Ankarberg A, McCabe J, Pinon AC, Schantz S, Emsley L. Atomic-level structure determination of amorphous molecular solids by NMR. Nat Commun 2023; 14:5138. [PMID: 37612269 PMCID: PMC10447443 DOI: 10.1038/s41467-023-40853-2] [Citation(s) in RCA: 4] [Impact Index Per Article: 4.0] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/08/2023] [Accepted: 08/10/2023] [Indexed: 08/25/2023] Open
Abstract
Structure determination of amorphous materials remains challenging, owing to the disorder inherent to these materials. Nuclear magnetic resonance (NMR) powder crystallography is a powerful method to determine the structure of molecular solids, but disorder leads to a high degree of overlap between measured signals, and prevents the unambiguous identification of a single modeled periodic structure as representative of the whole material. Here, we determine the atomic-level ensemble structure of the amorphous form of the drug AZD4625 by combining solid-state NMR experiments with molecular dynamics (MD) simulations and machine-learned chemical shifts. By considering the combined shifts of all 1H and 13C atomic sites in the molecule, we determine the structure of the amorphous form by identifying an ensemble of local molecular environments that are in agreement with experiment. We then extract and analyze preferred conformations and intermolecular interactions in the amorphous sample in terms of the stabilization of the amorphous form of the drug.
Collapse
Affiliation(s)
- Manuel Cordova
- Institut des Sciences et Ingénierie Chimiques, École Polytechnique Fédérale de Lausanne (EPFL), CH-1015, Lausanne, Switzerland
- National Centre for Computational Design and Discovery of Novel Materials MARVEL, École Polytechnique Fédérale de Lausanne (EPFL), Lausanne, Switzerland
| | - Pinelopi Moutzouri
- Institut des Sciences et Ingénierie Chimiques, École Polytechnique Fédérale de Lausanne (EPFL), CH-1015, Lausanne, Switzerland
| | - Sten O Nilsson Lill
- Data Science & Modelling, Pharmaceutical Sciences, R&D, AstraZeneca, Gothenburg, Sweden
| | - Alexander Cousen
- Early Chemical Development, Pharmaceutical Sciences, R&D, AstraZeneca, Macclesfield, UK
| | - Martin Kearns
- Early Product Development and Manufacturing, Pharmaceutical Sciences, R&D, AstraZeneca, Macclesfield, UK
| | - Stefan T Norberg
- Oral Product Development, Pharmaceutical Technology & Development, Operations, AstraZeneca, Gothenburg, Sweden
| | - Anna Svensk Ankarberg
- Oral Product Development, Pharmaceutical Technology & Development, Operations, AstraZeneca, Gothenburg, Sweden
| | - James McCabe
- Early Product Development and Manufacturing, Pharmaceutical Sciences, R&D, AstraZeneca, Macclesfield, UK
| | - Arthur C Pinon
- Swedish NMR Center, Department of Chemistry and Molecular Biology, University of Gothenburg, 41390, Gothenburg, Sweden
| | - Staffan Schantz
- Oral Product Development, Pharmaceutical Technology & Development, Operations, AstraZeneca, Gothenburg, Sweden.
| | - Lyndon Emsley
- Institut des Sciences et Ingénierie Chimiques, École Polytechnique Fédérale de Lausanne (EPFL), CH-1015, Lausanne, Switzerland.
- National Centre for Computational Design and Discovery of Novel Materials MARVEL, École Polytechnique Fédérale de Lausanne (EPFL), Lausanne, Switzerland.
| |
Collapse
|
11
|
Shahpasand-Kroner H, Siddique I, Malik R, Linares GR, Ivanova MI, Ichida J, Weil T, Münch J, Sanchez-Garcia E, Klärner FG, Schrader T, Bitan G. Molecular Tweezers: Supramolecular Hosts with Broad-Spectrum Biological Applications. Pharmacol Rev 2023; 75:263-308. [PMID: 36549866 PMCID: PMC9976797 DOI: 10.1124/pharmrev.122.000654] [Citation(s) in RCA: 4] [Impact Index Per Article: 4.0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/26/2022] [Revised: 10/14/2022] [Accepted: 10/19/2022] [Indexed: 12/24/2022] Open
Abstract
Lysine-selective molecular tweezers (MTs) are supramolecular host molecules displaying a remarkably broad spectrum of biologic activities. MTs act as inhibitors of the self-assembly and toxicity of amyloidogenic proteins using a unique mechanism. They destroy viral membranes and inhibit infection by enveloped viruses, such as HIV-1 and SARS-CoV-2, by mechanisms unrelated to their action on protein self-assembly. They also disrupt biofilm of Gram-positive bacteria. The efficacy and safety of MTs have been demonstrated in vitro, in cell culture, and in vivo, suggesting that these versatile compounds are attractive therapeutic candidates for various diseases, infections, and injuries. A lead compound called CLR01 has been shown to inhibit the aggregation of various amyloidogenic proteins, facilitate their clearance in vivo, prevent infection by multiple viruses, display potent anti-biofilm activity, and have a high safety margin in animal models. The inhibitory effect of CLR01 against amyloidogenic proteins is highly specific to abnormal self-assembly of amyloidogenic proteins with no disruption of normal mammalian biologic processes at the doses needed for inhibition. Therapeutic effects of CLR01 have been demonstrated in animal models of proteinopathies, lysosomal-storage diseases, and spinal-cord injury. Here we review the activity and mechanisms of action of these intriguing compounds and discuss future research directions. SIGNIFICANCE STATEMENT: Molecular tweezers are supramolecular host molecules with broad biological applications, including inhibition of abnormal protein aggregation, facilitation of lysosomal clearance of toxic aggregates, disruption of viral membranes, and interference of biofilm formation by Gram-positive bacteria. This review discusses the molecular and cellular mechanisms of action of the molecular tweezers, including the discovery of distinct mechanisms acting in vitro and in vivo, and the application of these compounds in multiple preclinical disease models.
Collapse
Affiliation(s)
- Hedieh Shahpasand-Kroner
- Department of Neurology, David Geffen School of Medicine (H.S.-K., I.S., R.M., G.B.), Brain Research Institute (G.B.), and Molecular Biology Institute (G.B.), University of California, Los Angeles, California; Department of Stem Cell Biology & Regenerative Medicine, Keck School of Medicine, University of Southern California, Los Angeles, California (G.R.L., J.I.); Department of Neurology, University of Michigan, Ann Arbor, Michigan (M.I.I.); Institute of Molecular Virology, Ulm University Medical Center, Ulm, Germany (T.W., J.M.); and Department of Computational Biochemistry (E.S.-G.) and Faculty of Chemistry (F-G.K., T.S.), University of Duisburg-Essen, Essen, Germany
| | - Ibrar Siddique
- Department of Neurology, David Geffen School of Medicine (H.S.-K., I.S., R.M., G.B.), Brain Research Institute (G.B.), and Molecular Biology Institute (G.B.), University of California, Los Angeles, California; Department of Stem Cell Biology & Regenerative Medicine, Keck School of Medicine, University of Southern California, Los Angeles, California (G.R.L., J.I.); Department of Neurology, University of Michigan, Ann Arbor, Michigan (M.I.I.); Institute of Molecular Virology, Ulm University Medical Center, Ulm, Germany (T.W., J.M.); and Department of Computational Biochemistry (E.S.-G.) and Faculty of Chemistry (F-G.K., T.S.), University of Duisburg-Essen, Essen, Germany
| | - Ravinder Malik
- Department of Neurology, David Geffen School of Medicine (H.S.-K., I.S., R.M., G.B.), Brain Research Institute (G.B.), and Molecular Biology Institute (G.B.), University of California, Los Angeles, California; Department of Stem Cell Biology & Regenerative Medicine, Keck School of Medicine, University of Southern California, Los Angeles, California (G.R.L., J.I.); Department of Neurology, University of Michigan, Ann Arbor, Michigan (M.I.I.); Institute of Molecular Virology, Ulm University Medical Center, Ulm, Germany (T.W., J.M.); and Department of Computational Biochemistry (E.S.-G.) and Faculty of Chemistry (F-G.K., T.S.), University of Duisburg-Essen, Essen, Germany
| | - Gabriel R Linares
- Department of Neurology, David Geffen School of Medicine (H.S.-K., I.S., R.M., G.B.), Brain Research Institute (G.B.), and Molecular Biology Institute (G.B.), University of California, Los Angeles, California; Department of Stem Cell Biology & Regenerative Medicine, Keck School of Medicine, University of Southern California, Los Angeles, California (G.R.L., J.I.); Department of Neurology, University of Michigan, Ann Arbor, Michigan (M.I.I.); Institute of Molecular Virology, Ulm University Medical Center, Ulm, Germany (T.W., J.M.); and Department of Computational Biochemistry (E.S.-G.) and Faculty of Chemistry (F-G.K., T.S.), University of Duisburg-Essen, Essen, Germany
| | - Magdalena I Ivanova
- Department of Neurology, David Geffen School of Medicine (H.S.-K., I.S., R.M., G.B.), Brain Research Institute (G.B.), and Molecular Biology Institute (G.B.), University of California, Los Angeles, California; Department of Stem Cell Biology & Regenerative Medicine, Keck School of Medicine, University of Southern California, Los Angeles, California (G.R.L., J.I.); Department of Neurology, University of Michigan, Ann Arbor, Michigan (M.I.I.); Institute of Molecular Virology, Ulm University Medical Center, Ulm, Germany (T.W., J.M.); and Department of Computational Biochemistry (E.S.-G.) and Faculty of Chemistry (F-G.K., T.S.), University of Duisburg-Essen, Essen, Germany
| | - Justin Ichida
- Department of Neurology, David Geffen School of Medicine (H.S.-K., I.S., R.M., G.B.), Brain Research Institute (G.B.), and Molecular Biology Institute (G.B.), University of California, Los Angeles, California; Department of Stem Cell Biology & Regenerative Medicine, Keck School of Medicine, University of Southern California, Los Angeles, California (G.R.L., J.I.); Department of Neurology, University of Michigan, Ann Arbor, Michigan (M.I.I.); Institute of Molecular Virology, Ulm University Medical Center, Ulm, Germany (T.W., J.M.); and Department of Computational Biochemistry (E.S.-G.) and Faculty of Chemistry (F-G.K., T.S.), University of Duisburg-Essen, Essen, Germany
| | - Tatjana Weil
- Department of Neurology, David Geffen School of Medicine (H.S.-K., I.S., R.M., G.B.), Brain Research Institute (G.B.), and Molecular Biology Institute (G.B.), University of California, Los Angeles, California; Department of Stem Cell Biology & Regenerative Medicine, Keck School of Medicine, University of Southern California, Los Angeles, California (G.R.L., J.I.); Department of Neurology, University of Michigan, Ann Arbor, Michigan (M.I.I.); Institute of Molecular Virology, Ulm University Medical Center, Ulm, Germany (T.W., J.M.); and Department of Computational Biochemistry (E.S.-G.) and Faculty of Chemistry (F-G.K., T.S.), University of Duisburg-Essen, Essen, Germany
| | - Jan Münch
- Department of Neurology, David Geffen School of Medicine (H.S.-K., I.S., R.M., G.B.), Brain Research Institute (G.B.), and Molecular Biology Institute (G.B.), University of California, Los Angeles, California; Department of Stem Cell Biology & Regenerative Medicine, Keck School of Medicine, University of Southern California, Los Angeles, California (G.R.L., J.I.); Department of Neurology, University of Michigan, Ann Arbor, Michigan (M.I.I.); Institute of Molecular Virology, Ulm University Medical Center, Ulm, Germany (T.W., J.M.); and Department of Computational Biochemistry (E.S.-G.) and Faculty of Chemistry (F-G.K., T.S.), University of Duisburg-Essen, Essen, Germany
| | - Elsa Sanchez-Garcia
- Department of Neurology, David Geffen School of Medicine (H.S.-K., I.S., R.M., G.B.), Brain Research Institute (G.B.), and Molecular Biology Institute (G.B.), University of California, Los Angeles, California; Department of Stem Cell Biology & Regenerative Medicine, Keck School of Medicine, University of Southern California, Los Angeles, California (G.R.L., J.I.); Department of Neurology, University of Michigan, Ann Arbor, Michigan (M.I.I.); Institute of Molecular Virology, Ulm University Medical Center, Ulm, Germany (T.W., J.M.); and Department of Computational Biochemistry (E.S.-G.) and Faculty of Chemistry (F-G.K., T.S.), University of Duisburg-Essen, Essen, Germany
| | - Frank-Gerrit Klärner
- Department of Neurology, David Geffen School of Medicine (H.S.-K., I.S., R.M., G.B.), Brain Research Institute (G.B.), and Molecular Biology Institute (G.B.), University of California, Los Angeles, California; Department of Stem Cell Biology & Regenerative Medicine, Keck School of Medicine, University of Southern California, Los Angeles, California (G.R.L., J.I.); Department of Neurology, University of Michigan, Ann Arbor, Michigan (M.I.I.); Institute of Molecular Virology, Ulm University Medical Center, Ulm, Germany (T.W., J.M.); and Department of Computational Biochemistry (E.S.-G.) and Faculty of Chemistry (F-G.K., T.S.), University of Duisburg-Essen, Essen, Germany
| | - Thomas Schrader
- Department of Neurology, David Geffen School of Medicine (H.S.-K., I.S., R.M., G.B.), Brain Research Institute (G.B.), and Molecular Biology Institute (G.B.), University of California, Los Angeles, California; Department of Stem Cell Biology & Regenerative Medicine, Keck School of Medicine, University of Southern California, Los Angeles, California (G.R.L., J.I.); Department of Neurology, University of Michigan, Ann Arbor, Michigan (M.I.I.); Institute of Molecular Virology, Ulm University Medical Center, Ulm, Germany (T.W., J.M.); and Department of Computational Biochemistry (E.S.-G.) and Faculty of Chemistry (F-G.K., T.S.), University of Duisburg-Essen, Essen, Germany
| | - Gal Bitan
- Department of Neurology, David Geffen School of Medicine (H.S.-K., I.S., R.M., G.B.), Brain Research Institute (G.B.), and Molecular Biology Institute (G.B.), University of California, Los Angeles, California; Department of Stem Cell Biology & Regenerative Medicine, Keck School of Medicine, University of Southern California, Los Angeles, California (G.R.L., J.I.); Department of Neurology, University of Michigan, Ann Arbor, Michigan (M.I.I.); Institute of Molecular Virology, Ulm University Medical Center, Ulm, Germany (T.W., J.M.); and Department of Computational Biochemistry (E.S.-G.) and Faculty of Chemistry (F-G.K., T.S.), University of Duisburg-Essen, Essen, Germany
| |
Collapse
|
12
|
Kong L, Bryce RA. Modeling pyranose ring pucker in carbohydrates using machine learning and semi-empirical quantum chemical methods. J Comput Chem 2022; 43:2009-2022. [PMID: 36165294 PMCID: PMC9828179 DOI: 10.1002/jcc.27000] [Citation(s) in RCA: 3] [Impact Index Per Article: 1.5] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/09/2022] [Revised: 08/17/2022] [Accepted: 08/22/2022] [Indexed: 01/12/2023]
Abstract
Pyranose ring pucker is a key coordinate governing the structure, interactions and reactivity of carbohydrates. We assess the ability of the machine learning potentials, ANI-1ccx and ANI-2x, and the GFN2-xTB semiempirical quantum chemical method, to model ring pucker conformers of five monosaccharides and oxane in the gas phase. Relative to coupled-cluster quantum mechanical calculations, we find that ANI-1ccx most accurately reproduces the ring pucker energy landscape for these molecules, with a correlation coefficient r2 of 0.83. This correlation in relative energies lowers to values of 0.70 for ANI-2x and 0.60 for GFN2-xTB. The ANI-1ccx also provides the most accurate estimate of the energetics of the 4 C1 -to-1 C4 minimum energy pathway for the six molecules. All three models reproduce chair more accurately than non-chair geometries. Analysis of small model molecules suggests that the ANI-1ccx model favors puckers with equatorial hydrogen bonding substituents; that ANI-2x and GFN2-xTB models overstabilize conformers with axially oriented groups; and that the endo-anomeric effect is overestimated by the machine learning models and underestimated via the GFN2-xTB method. While the pucker conformers considered in this study correspond to a gas phase environment, the accuracy and computational efficiency of the ANI-1ccx approach in modeling ring pucker in vacuo provides a promising basis for future evaluation and application to condensed phase environments.
Collapse
Affiliation(s)
- Linghan Kong
- Division of Pharmacy and Optometry, School of Health Sciences, Manchester Academic Health Sciences CentreUniversity of ManchesterManchesterUK
| | - Richard A. Bryce
- Division of Pharmacy and Optometry, School of Health Sciences, Manchester Academic Health Sciences CentreUniversity of ManchesterManchesterUK
| |
Collapse
|
13
|
Xu Y, Zhang S, Lindahl E, Friedman R, Wu W, Su P. A general tight-binding based energy decomposition analysis scheme for intermolecular interactions in large molecules. J Chem Phys 2022; 157:034104. [DOI: 10.1063/5.0091781] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
In this work, a general tight-binding based energy decomposition analysis (EDA) scheme for intermolecular interactions is proposed. Different from the earlier version [Xu et al., J. Chem. Phys. 154, 194106 (2021)], the current tight-binding based density functional theory (DFTB)-EDA is capable of performing interaction analysis with all the self-consistent charge (SCC) type DFTB methods, including SCC-DFTB2/3 and GFN1/2-xTB, despite their different formulas and parameterization schemes. In DFTB-EDA, the total interaction energy is divided into frozen, polarization, and dispersion terms. The performance of DFTB-EDA with SCC-DFTB2/3 and GFN1/2-xTB for various interaction systems is discussed and assessed.
Collapse
Affiliation(s)
- Yuan Xu
- The State Key Laboratory of Physical Chemistry of Solid Surfaces, Fujian Provincial Key Laboratory of Theoretical and Computational Chemistry, and College of Chemistry and Chemical Engineering, Xiamen University, Xiamen, Fujian 361005, China
| | - Shu Zhang
- The State Key Laboratory of Physical Chemistry of Solid Surfaces, Fujian Provincial Key Laboratory of Theoretical and Computational Chemistry, and College of Chemistry and Chemical Engineering, Xiamen University, Xiamen, Fujian 361005, China
| | - Erik Lindahl
- Department of Chemistry and Biomedical Sciences, Linnaeus University, 39182 Kalmar, Sweden
| | - Ran Friedman
- Department of Chemistry and Biomedical Sciences, Linnaeus University, 39182 Kalmar, Sweden
| | - Wei Wu
- The State Key Laboratory of Physical Chemistry of Solid Surfaces, Fujian Provincial Key Laboratory of Theoretical and Computational Chemistry, and College of Chemistry and Chemical Engineering, Xiamen University, Xiamen, Fujian 361005, China
| | - Peifeng Su
- The State Key Laboratory of Physical Chemistry of Solid Surfaces, Fujian Provincial Key Laboratory of Theoretical and Computational Chemistry, and College of Chemistry and Chemical Engineering, Xiamen University, Xiamen, Fujian 361005, China
| |
Collapse
|
14
|
Villot C, Ballesteros F, Wang D, Lao KU. Coupled Cluster Benchmarking of Large Noncovalent Complexes in L7 and S12L as Well as the C 60 Dimer, DNA-Ellipticine, and HIV-Indinavir. J Phys Chem A 2022; 126:4326-4341. [PMID: 35766331 DOI: 10.1021/acs.jpca.2c01421] [Citation(s) in RCA: 3] [Impact Index Per Article: 1.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
Abstract
In this work, we report the benchmark binding energies of the seven complexes within the L7 data set, six host-guest complexes from the S12L data set, a C60 dimer, the DNA-ellipticine intercalation complex, and the largest system of the study, the HIV-indinavir system, which contained 343 atoms or 139 heavy atoms. The high-quality values reported were obtained via a focal point method that relies on the canonical form of second-order Møller-Plesset theory and the domain-based local pair natural orbital scheme for the coupled cluster with single double and perturbative triple excitations [DLPNO-CCSD(T)] extrapolated to the complete basis set (CBS) limit. The results in this work not only corroborate but also improve upon some previous benchmark values for large noncovalent complexes albeit at a relatively steep cost. Although local CCSD(T) and the largely successful fixed-node diffusion Monte Carlo (FN-DMC) have been shown to generally agree for small- to medium-size systems, a discrepancy in their reported binding energy values arises for large complexes, where the magnitude of the disagreement is a definite cause for concern. For example, the largest deviation in the L7 data set was 2.8 kcal/mol (∼10%) on the low end in C3GC. Such a deviation only grows worse in the S12L set, which showed a difference of up to 10.4 kcal/mol (∼25%) by a conservative estimation in buckycatcher-C60. The DNA-ellipticine complex also generated a disagreement of 4.4 kcal/mol (∼10%) between both state-of-the-art methods. The disagreement between local CCSD(T) and FN-DMC in large noncovalent complexes shows that it is urgently needed to have the canonical CCSD(T), the Monte Carlo CCSD(T), or the full configuration interaction quantum Monte Carlo approaches available to large systems on the hundred-atom scale to solve this dilemma. In addition, the performances of cheaper popular computational methods were assessed for the studied complexes with respect to DLPNO-CCSD(T)/CBS. r2SCAN-3c, B97M-V, and PBE0+D4 work well in large noncovalent complexes in this work, and GFN2-xTB performs well in π-π stacking complexes. B97M-V is the most reliable computationally efficient approach to predicting noncovalent interactions for large complexes, being the only one to have binding errors within the so-called 1 kcal/mol "chemical accuracy". The benchmark interaction energies of these host-guest complexes, molecular materials, and biological systems with electronic and medicinal implications provide crucial reference data for the improvement of current and future lower-cost methods.
Collapse
Affiliation(s)
- Corentin Villot
- Department of Chemistry, Virginia Commonwealth University, Richmond, Virginia 23284 United States
| | - Francisco Ballesteros
- Department of Chemistry, Virginia Commonwealth University, Richmond, Virginia 23284 United States
| | - Danyang Wang
- Department of Chemistry, Virginia Commonwealth University, Richmond, Virginia 23284 United States
| | - Ka Un Lao
- Department of Chemistry, Virginia Commonwealth University, Richmond, Virginia 23284 United States
| |
Collapse
|
15
|
Jangrouei MR, Krzemińska A, Hapka M, Pastorczak E, Pernal K. Dispersion Interactions in Exciton-Localized States. Theory and Applications to π-π* and n-π* Excited States. J Chem Theory Comput 2022; 18:3497-3511. [PMID: 35587598 PMCID: PMC9202351 DOI: 10.1021/acs.jctc.2c00221] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Abstract
![]()
We address the problem
of intermolecular interaction energy calculations
in molecular complexes with localized excitons. Our focus is on the
correct representation of the dispersion energy. We derive an extended
Casimir-Polder formula for direct computation of this contribution
through second order in the intermolecular interaction operator V̂. An alternative formula, accurate to infinite order
in V̂, is derived within the framework of the
adiabatic connection (AC) theory. We also propose a new parametrization
of the VV10 nonlocal correlation density functional, so that it corrects
the CASSCF energy for the dispersion contribution and can be applied
to excited-state complexes. A numerical investigation is carried out
for benzene, pyridine, and peptide complexes with the local exciton
corresponding to the lowest π–π* or n– π*
states. The extended Casimir-Polder formula is implemented in the
framework of multiconfigurational symmetry-adapted perturbation theory,
SAPT(MC). A SAPT(MC) analysis shows that the creation of a localized
exciton affects mostly the electrostatic component of the interaction
energy of investigated complexes. Nevertheless, the changes in Pauli
repulsion and dispersion energies cannot be neglected. We verify the
performance of several perturbation- and AC-based methods. Best results
are obtained with a range-separated variant of an approximate AC approach
employing extended random phase approximation and CASSCF wave functions.
Collapse
Affiliation(s)
- Mohammad Reza Jangrouei
- Institute of Physics, Lodz University of Technology, ul. Wolczanska 217/221, 93-005, Lodz, Poland
| | - Agnieszka Krzemińska
- Institute of Physics, Lodz University of Technology, ul. Wolczanska 217/221, 93-005, Lodz, Poland
| | - Michał Hapka
- Faculty of Chemistry, University of Warsaw, ul. L. Pasteura 1, 02-093, Warsaw, Poland
| | - Ewa Pastorczak
- Institute of Physics, Lodz University of Technology, ul. Wolczanska 217/221, 93-005, Lodz, Poland
| | - Katarzyna Pernal
- Institute of Physics, Lodz University of Technology, ul. Wolczanska 217/221, 93-005, Lodz, Poland
| |
Collapse
|
16
|
Balodis M, Cordova M, Hofstetter A, Day GM, Emsley L. De Novo Crystal Structure Determination from Machine Learned Chemical Shifts. J Am Chem Soc 2022; 144:7215-7223. [PMID: 35416661 PMCID: PMC9052749 DOI: 10.1021/jacs.1c13733] [Citation(s) in RCA: 8] [Impact Index Per Article: 4.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/14/2022]
Abstract
Determination of the three-dimensional atomic-level structure of powdered solids is one of the key goals in current chemistry. Solid-state NMR chemical shifts can be used to solve this problem, but they are limited by the high computational cost associated with crystal structure prediction methods and density functional theory chemical shift calculations. Here, we successfully determine the crystal structures of ampicillin, piroxicam, cocaine, and two polymorphs of the drug molecule AZD8329 using on-the-fly generated machine-learned isotropic chemical shifts to directly guide a Monte Carlo-based structure determination process starting from a random gas-phase conformation.
Collapse
Affiliation(s)
- Martins Balodis
- Institut des Sciences et Ingénierie Chimiques, École Polytechnique Fédérale de Lausanne (EPFL), CH-1015 Lausanne, Switzerland
| | - Manuel Cordova
- Institut des Sciences et Ingénierie Chimiques, École Polytechnique Fédérale de Lausanne (EPFL), CH-1015 Lausanne, Switzerland.,National Centre for Computational Design and Discovery of Novel Materials MARVEL, École Polytechnique Fédérale de Lausanne (EPFL), Lausanne CH-1015, Switzerland
| | - Albert Hofstetter
- Institut des Sciences et Ingénierie Chimiques, École Polytechnique Fédérale de Lausanne (EPFL), CH-1015 Lausanne, Switzerland
| | - Graeme M Day
- School of Chemistry, University of Southampton, Highfield SO17 1BJ, Southampton, United Kingdom
| | - Lyndon Emsley
- Institut des Sciences et Ingénierie Chimiques, École Polytechnique Fédérale de Lausanne (EPFL), CH-1015 Lausanne, Switzerland.,National Centre for Computational Design and Discovery of Novel Materials MARVEL, École Polytechnique Fédérale de Lausanne (EPFL), Lausanne CH-1015, Switzerland
| |
Collapse
|
17
|
Fabregat R, Fabrizio A, Engel EA, Meyer B, Juraskova V, Ceriotti M, Corminboeuf C. Local Kernel Regression and Neural Network Approaches to the Conformational Landscapes of Oligopeptides. J Chem Theory Comput 2022; 18:1467-1479. [PMID: 35179897 PMCID: PMC8908737 DOI: 10.1021/acs.jctc.1c00813] [Citation(s) in RCA: 5] [Impact Index Per Article: 2.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
Abstract
![]()
The application of
machine learning to theoretical chemistry has
made it possible to combine the accuracy of quantum chemical energetics
with the thorough sampling of finite-temperature fluctuations. To
reach this goal, a diverse set of methods has been proposed, ranging
from simple linear models to kernel regression and highly nonlinear
neural networks. Here we apply two widely different approaches to
the same, challenging problem: the sampling of the conformational
landscape of polypeptides at finite temperature. We develop a local
kernel regression (LKR) coupled with a supervised sparsity method
and compare it with a more established approach based on Behler-Parrinello
type neural networks. In the context of the LKR, we discuss how the
supervised selection of the reference pool of environments is crucial
to achieve accurate potential energy surfaces at a competitive computational
cost and leverage the locality of the model to infer which chemical
environments are poorly described by the DFTB baseline. We then discuss
the relative merits of the two frameworks and perform Hamiltonian-reservoir
replica-exchange Monte Carlo sampling and metadynamics simulations,
respectively, to demonstrate that both frameworks can achieve converged
and transferable sampling of the conformational landscape of complex
and flexible biomolecules with comparable accuracy and computational
cost.
Collapse
Affiliation(s)
| | | | - Edgar A Engel
- Laboratory of Computational Science and Modeling, IMX, École Polytechnique Fédérale de Lausanne, 1015 Lausanne, Switzerland
| | | | | | - Michele Ceriotti
- Laboratory of Computational Science and Modeling, IMX, École Polytechnique Fédérale de Lausanne, 1015 Lausanne, Switzerland
| | | |
Collapse
|
18
|
Řezáč J. Non-Covalent Interactions Atlas Benchmark Data Sets 5: London Dispersion in an Extended Chemical Space. Phys Chem Chem Phys 2022; 24:14780-14793. [DOI: 10.1039/d2cp01602h] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022]
Abstract
The Non-Covalent Interactions Atlas (www.nciatlas.org) has been extended with two data sets of benchmark interaction energies in complexes dominated by London dispersion. The D1200 data set of equilibrium geometries provides...
Collapse
|
19
|
Cordova M, Balodis M, Simões de Almeida B, Ceriotti M, Emsley L. Bayesian probabilistic assignment of chemical shifts in organic solids. SCIENCE ADVANCES 2021; 7:eabk2341. [PMID: 34826232 PMCID: PMC8626066 DOI: 10.1126/sciadv.abk2341] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Download PDF] [Figures] [Subscribe] [Scholar Register] [Received: 06/30/2021] [Accepted: 10/08/2021] [Indexed: 06/13/2023]
Abstract
A prerequisite for NMR studies of organic materials is assigning each experimental chemical shift to a set of geometrically equivalent nuclei. Obtaining the assignment experimentally can be challenging and typically requires time-consuming multidimensional correlation experiments. An alternative solution for determining the assignment involves statistical analysis of experimental chemical shift databases, but no such database exists for molecular solids. Here, by combining the Cambridge Structural Database with a machine learning model of chemical shifts, we construct a statistical basis for probabilistic chemical shift assignment of organic crystals by calculating shifts for more than 200,000 compounds, enabling the probabilistic assignment of organic crystals directly from their two-dimensional chemical structure. The approach is demonstrated with the 13C and 1H assignment of 11 molecular solids with experimental shifts and benchmarked on 100 crystals using predicted shifts. The correct assignment was found among the two most probable assignments in more than 80% of cases.
Collapse
Affiliation(s)
- Manuel Cordova
- Laboratory of Magnetic Resonance, Institute of Chemical Sciences and Engineering, Ecole Polytechnique Fédérale de Lausanne, CH-1015 Lausanne, Switzerland
- National Centre for Computational Design and Discovery of Novel Materials MARVEL, Ecole Polytechnique Fédérale de Lausanne, CH-1015 Lausanne, Switzerland
| | - Martins Balodis
- Laboratory of Magnetic Resonance, Institute of Chemical Sciences and Engineering, Ecole Polytechnique Fédérale de Lausanne, CH-1015 Lausanne, Switzerland
| | - Bruno Simões de Almeida
- Laboratory of Magnetic Resonance, Institute of Chemical Sciences and Engineering, Ecole Polytechnique Fédérale de Lausanne, CH-1015 Lausanne, Switzerland
| | - Michele Ceriotti
- National Centre for Computational Design and Discovery of Novel Materials MARVEL, Ecole Polytechnique Fédérale de Lausanne, CH-1015 Lausanne, Switzerland
- Laboratory of Computational Science and Modelling, Institute of Materials, Ecole Polytechnique Fédérale de Lausanne, CH-1015 Lausanne, Switzerland
| | - Lyndon Emsley
- Laboratory of Magnetic Resonance, Institute of Chemical Sciences and Engineering, Ecole Polytechnique Fédérale de Lausanne, CH-1015 Lausanne, Switzerland
- National Centre for Computational Design and Discovery of Novel Materials MARVEL, Ecole Polytechnique Fédérale de Lausanne, CH-1015 Lausanne, Switzerland
| |
Collapse
|
20
|
Mondal P, Cazade PA, Das AK, Bereau T, Meuwly M. Multipolar Force Fields for Amide-I Spectroscopy from Conformational Dynamics of the Alanine Trimer. J Phys Chem B 2021; 125:10928-10938. [PMID: 34559531 DOI: 10.1021/acs.jpcb.1c05423] [Citation(s) in RCA: 6] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Abstract
The dynamics and spectroscopy of N-methyl-acetamide (NMA) and trialanine in solution are characterized from molecular dynamics simulations using different energy functions, including a conventional point charge (PC)-based force field, one based on a multipolar (MTP) representation of the electrostatics, and a semiempirical DFT method. For the 1D infrared spectra, the frequency splitting between the two amide-I groups is 10 cm-1 from the PC, 13 cm-1 from the MTP, and 47 cm-1 from self-consistent charge density functional tight-binding (SCC-DFTB) simulations, compared with 25 cm-1 from experiment. The frequency trajectory required for the frequency fluctuation correlation function (FFCF) is determined from individual normal mode (INM) and full normal mode (FNM) analyses of the amide-I vibrations. The spectroscopy, time-zero magnitude of the FFCF C(t = 0), and the static component Δ02 from simulations using MTP and analysis based on FNM are all consistent with experiments for (Ala)3. Contrary to this, for the analysis excluding mode-mode coupling (INM), the FFCF decays to zero too rapidly and for simulations with a PC-based force field, the Δ02 is too small by a factor of two compared with experiments. Simulations with SCC-DFTB agree better with experiment for these observables than those from PC-based simulations. The conformational ensemble sampled from simulations using PCs is consistent with the literature (including PII, β, αR, and αL), whereas that covered by the MTP-based simulations is dominated by PII with some contributions from β and αR. This agrees with and confirms recently reported Bayesian-refined populations based on 1D infrared experiments. FNM analysis together with a MTP representation provides a meaningful model to correctly describe the dynamics of hydrated trialanine.
Collapse
Affiliation(s)
- Padmabati Mondal
- Department of Chemistry, University of Basel, Klingelbergstrasse 80, Basel 4056, Switzerland
| | - Pierre-André Cazade
- Department of Chemistry, University of Basel, Klingelbergstrasse 80, Basel 4056, Switzerland
| | - Akshaya K Das
- Department of Chemistry, University of Basel, Klingelbergstrasse 80, Basel 4056, Switzerland
| | - Tristan Bereau
- Department of Chemistry, University of Basel, Klingelbergstrasse 80, Basel 4056, Switzerland
| | - Markus Meuwly
- Department of Chemistry, University of Basel, Klingelbergstrasse 80, Basel 4056, Switzerland.,Department of Chemistry, Brown University, Providence/RI 02912, United States
| |
Collapse
|
21
|
Charzewski Ł, Krzyśko KA, Lesyng B. Exploring Covalent Docking Mechanisms of Boron-Based Inhibitors to Class A, C and D β-Lactamases Using Time-dependent Hybrid QM/MM Simulations. Front Mol Biosci 2021; 8:633181. [PMID: 34434961 PMCID: PMC8380965 DOI: 10.3389/fmolb.2021.633181] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/24/2020] [Accepted: 07/19/2021] [Indexed: 11/13/2022] Open
Abstract
Recently, molecular covalent docking has been extensively developed to design new classes of inhibitors that form chemical bonds with their biological targets. This strategy for the design of such inhibitors, in particular boron-based inhibitors, holds great promise for the vast family of β-lactamases produced, inter alia, by Gram-negative antibiotic-resistant bacteria. However, the description of covalent docking processes requires a quantum-mechanical approach, and so far, only a few studies of this type have been presented. This study accurately describes the covalent docking process between two model inhibitors - representing two large families of inhibitors based on boronic-acid and bicyclic boronate scaffolds, and three β-lactamases which belong to the A, C, and D classes. Molecular fragments containing boron can be converted from a neutral, trigonal, planar state with sp2 hybridization to the anionic, tetrahedral sp3 state in a process sometimes referred to as morphing. This study applies multi-scale modeling methods, in particular, the hybrid QM/MM approach which has predictive power reaching well beyond conventional molecular modeling. Time-dependent QM/MM simulations indicated several structural changes and geometric preferences, ultimately leading to covalent docking processes. With current computing technologies, this approach is not computationally expensive, can be used in standard molecular modeling and molecular design works, and can effectively support experimental research which should allow for a detailed understanding of complex processes important to molecular medicine. In particular, it can support the rational design of covalent boron-based inhibitors for β-lactamases as well as for many other enzyme systems of clinical relevance, including SARS-CoV-2 proteins.
Collapse
Affiliation(s)
| | | | - Bogdan Lesyng
- Department of Biophysics, Faculty of Physics, University of Warsaw, Warsaw, Poland
| |
Collapse
|
22
|
Milovanović B, Etinski M, Popov I. Self-assembly of rylene-decorated guanine ribbons on graphene surface for optoelectronic applications: a theoretical study. NANOTECHNOLOGY 2021; 32:435405. [PMID: 34284367 DOI: 10.1088/1361-6528/ac162c] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Subscribe] [Scholar Register] [Received: 04/07/2021] [Accepted: 07/20/2021] [Indexed: 06/13/2023]
Abstract
We are witnessing a change of paradigm from the conventional top-down to the bottom-up fabrication of nanodevices and particularly optoelectronic devices. A promising example of the bottom-up approach is self-assembling of molecules into layers with predictable and reproducible structural, electronic and optical properties. Nucleobases possess extraordinary ability to self-assembly into one-, two-, and three-dimensional structures. Optical properties of nucleotides are not suitable for wider application to optoelectronics and photovoltaics due to their large optical band gap, which is in contrast to rylene-based dyes that have been intensively investigated in organic optoelectronics. However, these lack the self-assembly capability of nucleobases. Combinations of covalently decorated guanine molecules with rylene type chromophores present 'the best of the both worlds'. Due to the large size of such compounds and its flexible nature their self-assemblies have not been fully understood yet. Here, we use a theoretical approach to study the structural, energetic and optical properties of rylene-based dye decorated guanine (GPDI), as self-assembled on a graphene sheet. Particularly we utilize the density-functional based tight-binding method to study atomic structure of these systems including the potential energy surface of GPDI and stability and organization of single- and multilayered GPDIs on graphene sheet. Using density-functional theory (DFT) we employ the energy decomposition analysis to gain a deeper insight into the contributions of different moieties to stability of GPDI films. Using time dependent DFT we analyze optical properties of these systems. We find that atomically thin films consisting of only a few molecular layers with large surface areas are more favorable than isolated thick islands. Our study of excited states indicates existence of charge separated states similar to ones found in the well-studied hydrogen bonded organic frameworks. The self-assembly characterized with a large homogeneous coverage and long-living charge-separated states provide the great potential for optoelectronic applications.
Collapse
Affiliation(s)
- Branislav Milovanović
- University of Belgrade, Faculty of Physical Chemistry, Studentski trg 12-16, Belgrade, Serbia
| | - Mihajlo Etinski
- University of Belgrade, Faculty of Physical Chemistry, Studentski trg 12-16, Belgrade, Serbia
| | - Igor Popov
- Institute for Multidisciplinary Research, University of Belgrade, Kneza Višeslava 1, Belgrade, Serbia
- Institut of Physics, University of Belgrade, Pregrevica 118, Belgrade, Serbia
| |
Collapse
|
23
|
Xu Y, Friedman R, Wu W, Su P. Understanding intermolecular interactions of large systems in ground state and excited state by using density functional based tight binding methods. J Chem Phys 2021; 154:194106. [PMID: 34240911 DOI: 10.1063/5.0052060] [Citation(s) in RCA: 9] [Impact Index Per Article: 3.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/15/2022] Open
Abstract
A novel energy decomposition analysis scheme, named DFTB-EDA, is proposed based on the density functional based tight-binding method (DFTB/TD-DFTB), which is a semi-empirical quantum mechanical method based on Kohn-Sham-DFT for large-scale calculations. In DFTB-EDA, the total interaction energy is divided into three terms: frozen density, polarization, and dispersion. Owing to the small cost of DFTB/TD-DFTB, DFTB-EDA is capable of analyzing intermolecular interactions in large molecular systems containing several thousand atoms with high computational efficiency. It can be used not only for ground states but also for excited states. Test calculations, involving the S66 and L7 databases, several large molecules, and non-covalent bonding complexes in their lowest excited states, demonstrate the efficiency, usefulness, and capabilities of DFTB-EDA. Finally, the limits of DFTB-EDA are pointed out.
Collapse
Affiliation(s)
- Yuan Xu
- The State Key Laboratory of Physical Chemistry of Solid Surfaces, Fujian Provincial Key Laboratory of Theoretical and Computational Chemistry, and College of Chemistry and Chemical Engineering, Xiamen University, Xiamen, Fujian 361005, China
| | - Ran Friedman
- Department of Chemistry and Biomedical Sciences, Linnaeus University, 39182 Kalmar, Sweden
| | - Wei Wu
- The State Key Laboratory of Physical Chemistry of Solid Surfaces, Fujian Provincial Key Laboratory of Theoretical and Computational Chemistry, and College of Chemistry and Chemical Engineering, Xiamen University, Xiamen, Fujian 361005, China
| | - Peifeng Su
- The State Key Laboratory of Physical Chemistry of Solid Surfaces, Fujian Provincial Key Laboratory of Theoretical and Computational Chemistry, and College of Chemistry and Chemical Engineering, Xiamen University, Xiamen, Fujian 361005, China
| |
Collapse
|
24
|
Cordova M, Balodis M, Hofstetter A, Paruzzo F, Nilsson Lill SO, Eriksson ESE, Berruyer P, Simões de Almeida B, Quayle MJ, Norberg ST, Svensk Ankarberg A, Schantz S, Emsley L. Structure determination of an amorphous drug through large-scale NMR predictions. Nat Commun 2021; 12:2964. [PMID: 34016980 PMCID: PMC8137699 DOI: 10.1038/s41467-021-23208-7] [Citation(s) in RCA: 27] [Impact Index Per Article: 9.0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/04/2020] [Accepted: 03/25/2021] [Indexed: 11/17/2022] Open
Abstract
Knowledge of the structure of amorphous solids can direct, for example, the optimization of pharmaceutical formulations, but atomic-level structure determination in amorphous molecular solids has so far not been possible. Solid-state nuclear magnetic resonance (NMR) is among the most popular methods to characterize amorphous materials, and molecular dynamics (MD) simulations can help describe the structure of disordered materials. However, directly relating MD to NMR experiments in molecular solids has been out of reach until now because of the large size of these simulations. Here, using a machine learning model of chemical shifts, we determine the atomic-level structure of the hydrated amorphous drug AZD5718 by combining dynamic nuclear polarization-enhanced solid-state NMR experiments with predicted chemical shifts for MD simulations of large systems. From these amorphous structures we then identify H-bonding motifs and relate them to local intermolecular complex formation energies.
Collapse
Affiliation(s)
- Manuel Cordova
- Institut des Sciences et Ingénierie Chimiques, École Polytechnique Fédérale de Lausanne (EPFL), Lausanne, Switzerland
- National Centre for Computational Design and Discovery of Novel Materials MARVEL, École Polytechnique Fédérale de Lausanne (EPFL), Lausanne, Switzerland
| | - Martins Balodis
- Institut des Sciences et Ingénierie Chimiques, École Polytechnique Fédérale de Lausanne (EPFL), Lausanne, Switzerland
| | - Albert Hofstetter
- Institut des Sciences et Ingénierie Chimiques, École Polytechnique Fédérale de Lausanne (EPFL), Lausanne, Switzerland
| | - Federico Paruzzo
- Institut des Sciences et Ingénierie Chimiques, École Polytechnique Fédérale de Lausanne (EPFL), Lausanne, Switzerland
| | - Sten O Nilsson Lill
- Early Product Development and Manufacturing, Pharmaceutical Sciences, R&D, AstraZeneca, Gothenburg, Sweden
| | - Emma S E Eriksson
- Early Product Development and Manufacturing, Pharmaceutical Sciences, R&D, AstraZeneca, Gothenburg, Sweden
| | - Pierrick Berruyer
- Institut des Sciences et Ingénierie Chimiques, École Polytechnique Fédérale de Lausanne (EPFL), Lausanne, Switzerland
| | - Bruno Simões de Almeida
- Institut des Sciences et Ingénierie Chimiques, École Polytechnique Fédérale de Lausanne (EPFL), Lausanne, Switzerland
| | - Michael J Quayle
- New Modalities and Parenteral Development, Pharmaceutical Technology & Development, Operations, AstraZeneca, Gothenburg, Sweden
| | - Stefan T Norberg
- Oral Product Development, Pharmaceutical Technology & Development, Operations, AstraZeneca, Gothenburg, Sweden
| | - Anna Svensk Ankarberg
- Oral Product Development, Pharmaceutical Technology & Development, Operations, AstraZeneca, Gothenburg, Sweden
| | - Staffan Schantz
- Oral Product Development, Pharmaceutical Technology & Development, Operations, AstraZeneca, Gothenburg, Sweden.
| | - Lyndon Emsley
- Institut des Sciences et Ingénierie Chimiques, École Polytechnique Fédérale de Lausanne (EPFL), Lausanne, Switzerland.
- National Centre for Computational Design and Discovery of Novel Materials MARVEL, École Polytechnique Fédérale de Lausanne (EPFL), Lausanne, Switzerland.
| |
Collapse
|
25
|
Karimova NV, Alves MR, Luo M, Grassian VH, Gerber RB. Toward a microscopic model of light absorbing dissolved organic compounds in aqueous environments: theoretical and experimental study. Phys Chem Chem Phys 2021; 23:10487-10497. [PMID: 33899856 DOI: 10.1039/d0cp06554d] [Citation(s) in RCA: 4] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022]
Abstract
Water systems often contain complex macromolecular systems that absorb light. In marine environments, these light absorbing components are often at the air-water interface and can participate in the chemistry of the atmosphere in ways that are poorly understood. Understanding the photochemistry and photophysics of these systems represents a major challenge since their composition and structures are not unique. In this study, we present a successful microscopic model of this light absorbing macromolecular species termed "marine derived chromophoric dissolved organic matter" or "m-CDOM" in water. The approach taken involves molecular dynamics simulations in the ground state using on the fly Density Functional Tight-Binding (DFTB) electronic structure theory; Time Dependent DFTB (TD-DFTB) calculations of excited states, and experimental measurements of the optical absorption spectra in aqueous solution. The theoretical hydrated model shows key features seen in the experimental data for a collected m-CDOM sample. As will be discussed, insights from the model are: (i) the low-energy A-band (at 410 nm) is due to the carbon chains combined with the diol- and the oxy-groups present in the structure; (ii) the weak B-band (at 320-360 nm) appears due to the contribution of the ionized speciated form of m-CDOM; and (iii) the higher-energy C-band (at 280 nm) is due to the two fused ring system. Thus, this is a two-speciated formed model. Although a relatively simple system, these calculations represent an important step in understanding light absorbing compounds found in nature and the search for other microscopic models of related materials remains of major interest.
Collapse
Affiliation(s)
- Natalia V Karimova
- Department of Chemistry, University of California, Irvine, CA 92697, USA.
| | - Michael R Alves
- Department of Chemistry and Biochemistry, University of California, San Diego, CA 92093, USA.
| | - Man Luo
- Department of Chemistry and Biochemistry, University of California, San Diego, CA 92093, USA.
| | - Vicki H Grassian
- Department of Chemistry and Biochemistry, University of California, San Diego, CA 92093, USA. and Department of Nanoengineering and Scripps Institution of Oceanography, University of California, San Diego, CA 92093, USA
| | - R Benny Gerber
- Department of Chemistry, University of California, Irvine, CA 92697, USA. and Institute of Chemistry and Fritz Haber Research Center, Hebrew University of Jerusalem, Jerusalem 91904, Israel
| |
Collapse
|
26
|
Higuchi Y, Asano Y, Kuwahara T, Hishida M. Rotational Dynamics of Water at the Phospholipid Bilayer Depending on the Head Groups Studied by Molecular Dynamics Simulations. LANGMUIR : THE ACS JOURNAL OF SURFACES AND COLLOIDS 2021; 37:5329-5338. [PMID: 33881324 DOI: 10.1021/acs.langmuir.1c00417] [Citation(s) in RCA: 6] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 06/12/2023]
Abstract
Hydration states are a crucial factor that affect the self-assembly and properties of soft materials and biomolecules. Although previous experiments have revealed that the hydration state strongly depends on the chemical structure of lipid molecules, the mechanisms at the molecular level remain unknown. Classical and density-functional tight-binding (DFTB) molecular dynamics (MD) simulations are employed to determine the mechanisms underlying dissimilar water dynamics between lipid membranes with phosphatidylcholine (PC) and phosphatidylethanolamine (PE) head groups. The classical MD simulation shows that rotational relaxations of water are faster on the PE lipid than on the PC lipid, which is consistent with a previous experimental study using terahertz spectroscopy. Furthermore, DFTB-MD simulation of N(CH3)4+ and NH4+ ions, which correspond to the different head groups in PC and PE, respectively, shows qualitative agreement with the classical MD simulation. Remarkably, the PE lipids and the NH4+ ions break hydrogen bonds between water molecules in the secondary hydration shell. In contrast, the PC lipids and the N(CH3)4+ ions bind water molecules weakly in both the primary and secondary hydration shells and increase hydrogen bonds between water. Our atomistic simulations show that these changes in the hydrogen-bond network of water molecules cause the different rotational relaxation of water molecules between the two lipids.
Collapse
Affiliation(s)
- Yuji Higuchi
- Institute for Solid State Physics, The University of Tokyo, 5-1-5 Kashiwanoha, Kashiwa, Chiba 277-8581, Japan
| | - Yuta Asano
- Institute for Solid State Physics, The University of Tokyo, 5-1-5 Kashiwanoha, Kashiwa, Chiba 277-8581, Japan
| | - Takuya Kuwahara
- MicroTribology Center μTC, Fraunhofer IWM, Wöhlerstraße 11, Freiburg 79108, Germany
| | - Mafumi Hishida
- Department of Chemistry, Faculty of Pure and Applied Sciences, University of Tsukuba, Tsukuba 305-8571, Ibaraki, Japan
| |
Collapse
|
27
|
Sulimov VB, Kutov DC, Taschilova AS, Ilin IS, Tyrtyshnikov EE, Sulimov AV. Docking Paradigm in Drug Design. Curr Top Med Chem 2021; 21:507-546. [PMID: 33292135 DOI: 10.2174/1568026620666201207095626] [Citation(s) in RCA: 13] [Impact Index Per Article: 4.3] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/01/2020] [Revised: 09/28/2020] [Accepted: 10/16/2020] [Indexed: 11/22/2022]
Abstract
Docking is in demand for the rational computer aided structure based drug design. A review of docking methods and programs is presented. Different types of docking programs are described. They include docking of non-covalent small ligands, protein-protein docking, supercomputer docking, quantum docking, the new generation of docking programs and the application of docking for covalent inhibitors discovery. Taking into account the threat of COVID-19, we present here a short review of docking applications to the discovery of inhibitors of SARS-CoV and SARS-CoV-2 target proteins, including our own result of the search for inhibitors of SARS-CoV-2 main protease using docking and quantum chemical post-processing. The conclusion is made that docking is extremely important in the fight against COVID-19 during the process of development of antivirus drugs having a direct action on SARS-CoV-2 target proteins.
Collapse
Affiliation(s)
- Vladimir B Sulimov
- Research Computer Center of Lomonosov Moscow State University, Moscow, Russian Federation
| | - Danil C Kutov
- Research Computer Center of Lomonosov Moscow State University, Moscow, Russian Federation
| | - Anna S Taschilova
- Research Computer Center of Lomonosov Moscow State University, Moscow, Russian Federation
| | - Ivan S Ilin
- Research Computer Center of Lomonosov Moscow State University, Moscow, Russian Federation
| | - Eugene E Tyrtyshnikov
- Institute of Numerical Mathematics of Russian Academy of Sciences, Moscow, Russian Federation
| | - Alexey V Sulimov
- Research Computer Center of Lomonosov Moscow State University, Moscow, Russian Federation
| |
Collapse
|
28
|
Gutten O, Jurečka P, Aliakbar Tehrani Z, Buděšínský M, Řezáč J, Rulíšek L. Conformational energies and equilibria of cyclic dinucleotides in vacuo and in solution: computational chemistry vs. NMR experiments. Phys Chem Chem Phys 2021; 23:7280-7294. [PMID: 33876088 DOI: 10.1039/d0cp05993e] [Citation(s) in RCA: 5] [Impact Index Per Article: 1.7] [Reference Citation Analysis] [Abstract] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/22/2022]
Abstract
Performance of computational methods in modelling cyclic dinucleotides - an important and challenging class of compounds - has been evaluated by two different benchmarks: (1) gas-phase conformational energies and (2) qualitative agreement with NMR observations of the orientation of the χ-dihedral angle in solvent. In gas-phase benchmarks, where CCSD(T) and DLPNO-CCSD(T) methods have been used as the reference, most of the (dispersion corrected) density functional approximations are accurate enough to justify prioritizing computational cost and compatibility with other modelling options as the criterion of choice. NMR experiments of 3'3'-c-di-AMP, 3'3'-c-GAMP, and 3'3'-c-di-GMP show the overall prevalence of the anti-conformation of purine bases, but some population of syn-conformations is observed for guanines. Implicit solvation models combined with quantum-chemical methods struggle to reproduce this behaviour, probably due to a lack of dynamics and explicitly modelled solvent, leading to structures that are too compact. Molecular dynamics simulations overrepresent the syn-conformation of guanine due to the overestimation of an intramolecular hydrogen bond. Our combination of experimental and computational benchmarks provides "error bars" for modelling cyclic dinucleotides in solvent, where such information is generally difficult to obtain, and should help gauge the interpretability of studies dealing with binding of cyclic dinucleotides to important pharmaceutical targets. At the same time, the presented analysis calls for improvement in both implicit solvation models and force-field parameters.
Collapse
Affiliation(s)
- Ondrej Gutten
- Institute of Organic Chemistry and Biochemistry of the Czech Academy of Sciences, Flemingovo náměstí 2, 166 10, Praha 6, Czech Republic.
| | | | | | | | | | | |
Collapse
|
29
|
Tandiana R, Brun E, Sicard-Roselli C, Domin D, Van-Oanh NT, Clavaguéra C. Probing the structural properties of the water solvation shell around gold nanoparticles: A computational study. J Chem Phys 2021; 154:044706. [DOI: 10.1063/5.0037551] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/06/2023] Open
Affiliation(s)
- Rika Tandiana
- Université Paris-Saclay, CNRS, Institut de Chimie Physique, UMR8000, 91405 Orsay, France
| | - Emilie Brun
- Université Paris-Saclay, CNRS, Institut de Chimie Physique, UMR8000, 91405 Orsay, France
| | - Cécile Sicard-Roselli
- Université Paris-Saclay, CNRS, Institut de Chimie Physique, UMR8000, 91405 Orsay, France
| | - Dominik Domin
- Université Paris-Saclay, CNRS, Institut de Chimie Physique, UMR8000, 91405 Orsay, France
| | - Nguyen-Thi Van-Oanh
- Université Paris-Saclay, CNRS, Institut de Chimie Physique, UMR8000, 91405 Orsay, France
| | - Carine Clavaguéra
- Université Paris-Saclay, CNRS, Institut de Chimie Physique, UMR8000, 91405 Orsay, France
| |
Collapse
|
30
|
Cook C, Beran GJO. Reduced-cost supercell approach for computing accurate phonon density of states in organic crystals. J Chem Phys 2020; 153:224105. [PMID: 33317313 DOI: 10.1063/5.0032649] [Citation(s) in RCA: 6] [Impact Index Per Article: 1.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/10/2023] Open
Abstract
Phonon contributions to organic crystal structures and thermochemical properties can be significant, but computing a well-converged phonon density of states with lattice dynamics and periodic density functional theory (DFT) is often computationally expensive due to the need for large supercells. Using semi-empirical methods like density functional tight binding (DFTB) instead of DFT can reduce the computational costs dramatically, albeit with noticeable reductions in accuracy. This work proposes approximating the phonon density of states via a relatively inexpensive DFTB supercell treatment of the phonon dispersion that is then corrected by shifting the individual phonon modes according to the difference between the DFT and DFTB phonon frequencies at the Γ-point. The acoustic modes are then computed at the DFT level from the elastic constants. In several small-molecule crystal test cases, this combined approach reproduces DFT thermochemistry with kJ/mol accuracy and 1-2 orders of magnitude less computational effort. Finally, this approach is applied to computing the free energy differences between the five crystal polymorphs of oxalyl dihydrazide.
Collapse
Affiliation(s)
- Cameron Cook
- Department of Chemistry, University of California, Riverside, California 92521, USA
| | - Gregory J O Beran
- Department of Chemistry, University of California, Riverside, California 92521, USA
| |
Collapse
|
31
|
Chen Z, Liu Y, Liao S, Yi N, Hu Q. Treatment of adsorption of dioxane by using SiCNT toward efficient remediation of refractory organic contaminants from wastewater: DFT and DFTB-MD simulations. J Mol Liq 2020. [DOI: 10.1016/j.molliq.2020.113869] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/07/2023]
|
32
|
Řezáč J. Non-Covalent Interactions Atlas Benchmark Data Sets 2: Hydrogen Bonding in an Extended Chemical Space. J Chem Theory Comput 2020; 16:6305-6316. [DOI: 10.1021/acs.jctc.0c00715] [Citation(s) in RCA: 21] [Impact Index Per Article: 5.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Affiliation(s)
- Jan Řezáč
- Institute of Organic Chemistry and Biochemistry, Czech Academy of Sciences, 166 10 Prague, Czech Republic
| |
Collapse
|
33
|
Ono J, Imai M, Nishimura Y, Nakai H. Hydroxide Ion Carrier for Proton Pumps in Bacteriorhodopsin: Primary Proton Transfer. J Phys Chem B 2020; 124:8524-8539. [DOI: 10.1021/acs.jpcb.0c05507] [Citation(s) in RCA: 9] [Impact Index Per Article: 2.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Affiliation(s)
- Junichi Ono
- Waseda Research Institute for Science and Engineering (WISE), Waseda University, 3-4-1 Okubo, Shinjuku, Tokyo 169-8555, Japan
- Elements Strategy Initiative for Catalysts & Batteries (ESICB), Kyoto University, 1-30 Goryo-Ohara, Nishikyo-ku, Kyoto 615-8245, Japan
| | - Minori Imai
- Department of Chemistry and Biochemistry, School of Advanced Science and Engineering, Waseda University, 3-4-1 Okubo, Shinjuku, Tokyo 169-8555, Japan
| | - Yoshifumi Nishimura
- Waseda Research Institute for Science and Engineering (WISE), Waseda University, 3-4-1 Okubo, Shinjuku, Tokyo 169-8555, Japan
| | - Hiromi Nakai
- Waseda Research Institute for Science and Engineering (WISE), Waseda University, 3-4-1 Okubo, Shinjuku, Tokyo 169-8555, Japan
- Elements Strategy Initiative for Catalysts & Batteries (ESICB), Kyoto University, 1-30 Goryo-Ohara, Nishikyo-ku, Kyoto 615-8245, Japan
- Department of Chemistry and Biochemistry, School of Advanced Science and Engineering, Waseda University, 3-4-1 Okubo, Shinjuku, Tokyo 169-8555, Japan
| |
Collapse
|
34
|
Boyes M, Alieva A, Tong J, Nagyte V, Melle-Franco M, Vetter T, Casiraghi C. Exploiting the Surface Properties of Graphene for Polymorph Selectivity. ACS NANO 2020; 14:10394-10401. [PMID: 32692539 PMCID: PMC7605719 DOI: 10.1021/acsnano.0c04183] [Citation(s) in RCA: 8] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Subscribe] [Scholar Register] [Received: 05/19/2020] [Accepted: 07/21/2020] [Indexed: 05/29/2023]
Abstract
Producing crystals of the desired form (polymorph) is currently a challenge as nucleation is yet to be fully understood. Templated crystallization is an efficient approach to achieve polymorph selectivity; however, it is still unclear how to design the template to achieve selective crystallization of specific polymorphs. More insights into the nanoscale interactions happening during nucleation are needed. In this work, we investigate crystallization of glycine using graphene, with different surface chemistry, as a template. We show that graphene induces the preferential crystallization of the metastable α-polymorph compared to the unstable β-form at the contact region of an evaporating droplet. Computer modeling indicates the presence of a small amount of oxidized moieties on graphene to be responsible for the increased stabilization of the α-form. In conclusion, our work shows that graphene could become an attractive material for polymorph selectivity and screening by exploiting its tunable surface chemistry.
Collapse
Affiliation(s)
- Matthew Boyes
- Department
of Chemistry, University of Manchester, Manchester M13 9PL, United Kingdom
| | - Adriana Alieva
- Department
of Chemistry, University of Manchester, Manchester M13 9PL, United Kingdom
| | - Jincheng Tong
- Department
of Chemistry, University of Manchester, Manchester M13 9PL, United Kingdom
| | - Vaiva Nagyte
- Department
of Chemistry, University of Manchester, Manchester M13 9PL, United Kingdom
| | - Manuel Melle-Franco
- CICECO−Aveiro
Institute of Materials, Department of Chemistry, University of Aveiro, Aveiro 3810-193, Portugal
| | - Thomas Vetter
- Department
of Chemical Engineering and Analytical Science, University of Manchester, Manchester M13 9PL, United Kingdom
| | - Cinzia Casiraghi
- Department
of Chemistry, University of Manchester, Manchester M13 9PL, United Kingdom
| |
Collapse
|
35
|
Rossi K, Jurásková V, Wischert R, Garel L, Corminbœuf C, Ceriotti M. Simulating Solvation and Acidity in Complex Mixtures with First-Principles Accuracy: The Case of CH 3SO 3H and H 2O 2 in Phenol. J Chem Theory Comput 2020; 16:5139-5149. [PMID: 32567854 DOI: 10.1021/acs.jctc.0c00362] [Citation(s) in RCA: 19] [Impact Index Per Article: 4.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/18/2022]
Abstract
We present a generally applicable computational framework for the efficient and accurate characterization of molecular structural patterns and acid properties in an explicit solvent using H2O2 and CH3SO3H in phenol as an example. To address the challenges posed by the complexity of the problem, we resort to a set of data-driven methods and enhanced sampling algorithms. The synergistic application of these techniques makes the first-principle estimation of the chemical properties feasible without renouncing to the use of explicit solvation, involving extensive statistical sampling. Ensembles of neural network (NN) potentials are trained on a set of configurations carefully selected out of preliminary simulations performed at a low-cost density functional tight-binding (DFTB) level. The energy and forces of these configurations are then recomputed at the hybrid density functional theory (DFT) level and used to train the neural networks. The stability of the NN model is enhanced by using DFTB energetics as a baseline, but the efficiency of the direct NN (i.e., baseline-free) is exploited via a multiple-time-step integrator. The neural network potentials are combined with enhanced sampling techniques, such as replica exchange and metadynamics, and used to characterize the relevant protonated species and dominant noncovalent interactions in the mixture, also considering nuclear quantum effects.
Collapse
Affiliation(s)
- Kevin Rossi
- Laboratory of Computational Science and Modeling (COSMO), Institute of Materials, Ecole Polytechnique Fédérale de Lausanne (EPFL), Lausanne 1015, Switzerland
| | - Veronika Jurásková
- Laboratory for Computational Molecular Design (LCMD), Institute of Chemical Sciences and Engineering, Ecole Polytechnique Fédérale de Lausanne (EPFL), Lausanne 1015, Switzerland
| | - Raphael Wischert
- Eco-Efficient Products and Processes Laboratory, Solvay, RIC Shanghai, Shanghai 201108, China
| | - Laurent Garel
- Aroma Performance Laboratory, Solvay, RIC Lyon, 69190 Saint-Fons, France
| | - Clémence Corminbœuf
- Laboratory for Computational Molecular Design (LCMD), Institute of Chemical Sciences and Engineering, Ecole Polytechnique Fédérale de Lausanne (EPFL), Lausanne 1015, Switzerland.,National Centre for Computational Design and Discovery of Novel Materials (MARVEL), École Polytechnique Fédérale de Lausanne, 1015 Lausanne, Switzerland
| | - Michele Ceriotti
- Laboratory of Computational Science and Modeling (COSMO), Institute of Materials, Ecole Polytechnique Fédérale de Lausanne (EPFL), Lausanne 1015, Switzerland.,National Centre for Computational Design and Discovery of Novel Materials (MARVEL), École Polytechnique Fédérale de Lausanne, 1015 Lausanne, Switzerland
| |
Collapse
|
36
|
Pecina A, Eyrilmez SM, Köprülüoğlu C, Miriyala VM, Lepšík M, Fanfrlík J, Řezáč J, Hobza P. SQM/COSMO Scoring Function: Reliable Quantum-Mechanical Tool for Sampling and Ranking in Structure-Based Drug Design. Chempluschem 2020; 85:2362-2371. [PMID: 32609421 DOI: 10.1002/cplu.202000120] [Citation(s) in RCA: 6] [Impact Index Per Article: 1.5] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/20/2020] [Revised: 05/27/2020] [Indexed: 12/17/2022]
Abstract
Quantum mechanical (QM) methods have been gaining importance in structure-based drug design where a reliable description of protein-ligand interactions is of utmost significance. However, strategies i. e. QM/MM, fragmentation or semiempirical (SQM) methods had to be pursued to overcome the unfavorable scaling of QM methods. Various SQM-based approaches have significantly contributed to the accuracy of docking and improvement of lead compounds. Parametrizations of SQM and implicit solvent methods in our laboratory have been instrumental to obtain a reliable SQM-based scoring function. The experience gained in its application for activity ranking of ligands binding to tens of protein targets resulted in setting up a faster SQM/COSMO scoring approach, which outperforms standard scoring methods in native pose identification for two dozen protein targets with ten thousand poses. Recently, SQM/COSMO was effectively applied in a proof-of-concept study of enrichment in virtual screening. Due to its superior performance, feasibility and chemical generality, we propose the SQM/COSMO approach as an efficient tool in structure-based drug design.
Collapse
Affiliation(s)
- Adam Pecina
- Institute of Organic Chemistry, and Biochemistry of Czech Academy of Sciences, Flemingovo namesti 2, 166 10, Prague, Czech Republic
| | - Saltuk M Eyrilmez
- Institute of Organic Chemistry, and Biochemistry of Czech Academy of Sciences, Flemingovo namesti 2, 166 10, Prague, Czech Republic.,Regional Centre of Advanced Technologies and Materials, Department of Physical Chemistry, Palacky University, 771 46, Olomouc, Czech Republic
| | - Cemal Köprülüoğlu
- Institute of Organic Chemistry, and Biochemistry of Czech Academy of Sciences, Flemingovo namesti 2, 166 10, Prague, Czech Republic.,Regional Centre of Advanced Technologies and Materials, Department of Physical Chemistry, Palacky University, 771 46, Olomouc, Czech Republic
| | - Vijay Madhav Miriyala
- Institute of Organic Chemistry, and Biochemistry of Czech Academy of Sciences, Flemingovo namesti 2, 166 10, Prague, Czech Republic
| | - Martin Lepšík
- Institute of Organic Chemistry, and Biochemistry of Czech Academy of Sciences, Flemingovo namesti 2, 166 10, Prague, Czech Republic
| | - Jindřich Fanfrlík
- Institute of Organic Chemistry, and Biochemistry of Czech Academy of Sciences, Flemingovo namesti 2, 166 10, Prague, Czech Republic
| | - Jan Řezáč
- Institute of Organic Chemistry, and Biochemistry of Czech Academy of Sciences, Flemingovo namesti 2, 166 10, Prague, Czech Republic
| | - Pavel Hobza
- Institute of Organic Chemistry, and Biochemistry of Czech Academy of Sciences, Flemingovo namesti 2, 166 10, Prague, Czech Republic.,Regional Centre of Advanced Technologies and Materials, Department of Physical Chemistry, Palacky University, 771 46, Olomouc, Czech Republic
| |
Collapse
|
37
|
Hourahine B, Aradi B, Blum V, Bonafé F, Buccheri A, Camacho C, Cevallos C, Deshaye MY, Dumitrică T, Dominguez A, Ehlert S, Elstner M, van der Heide T, Hermann J, Irle S, Kranz JJ, Köhler C, Kowalczyk T, Kubař T, Lee IS, Lutsker V, Maurer RJ, Min SK, Mitchell I, Negre C, Niehaus TA, Niklasson AMN, Page AJ, Pecchia A, Penazzi G, Persson MP, Řezáč J, Sánchez CG, Sternberg M, Stöhr M, Stuckenberg F, Tkatchenko A, Yu VWZ, Frauenheim T. DFTB+, a software package for efficient approximate density functional theory based atomistic simulations. J Chem Phys 2020; 152:124101. [PMID: 32241125 DOI: 10.1063/1.5143190] [Citation(s) in RCA: 409] [Impact Index Per Article: 102.3] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/07/2023] Open
Abstract
DFTB+ is a versatile community developed open source software package offering fast and efficient methods for carrying out atomistic quantum mechanical simulations. By implementing various methods approximating density functional theory (DFT), such as the density functional based tight binding (DFTB) and the extended tight binding method, it enables simulations of large systems and long timescales with reasonable accuracy while being considerably faster for typical simulations than the respective ab initio methods. Based on the DFTB framework, it additionally offers approximated versions of various DFT extensions including hybrid functionals, time dependent formalism for treating excited systems, electron transport using non-equilibrium Green's functions, and many more. DFTB+ can be used as a user-friendly standalone application in addition to being embedded into other software packages as a library or acting as a calculation-server accessed by socket communication. We give an overview of the recently developed capabilities of the DFTB+ code, demonstrating with a few use case examples, discuss the strengths and weaknesses of the various features, and also discuss on-going developments and possible future perspectives.
Collapse
Affiliation(s)
- B Hourahine
- SUPA, Department of Physics, The University of Strathclyde, Glasgow G4 0NG, United Kingdom
| | - B Aradi
- Bremen Center for Computational Materials Science, University of Bremen, Bremen, Germany
| | - V Blum
- Department of Mechanical Engineering and Materials Science, Duke University, Durham, North Carolina 27708, USA
| | - F Bonafé
- Max Planck Institute for the Structure and Dynamics of Matter, Hamburg, Germany
| | - A Buccheri
- School of Chemistry, University of Bristol, Cantock's Close, Bristol BS8 1TS, United Kingdom
| | - C Camacho
- School of Chemistry, University of Costa Rica, San José 11501-2060, Costa Rica
| | - C Cevallos
- School of Chemistry, University of Costa Rica, San José 11501-2060, Costa Rica
| | - M Y Deshaye
- Department of Chemistry and Advanced Materials Science and Engineering Center, Western Washington University, Bellingham, Washington 98225, USA
| | - T Dumitrică
- Department of Mechanical Engineering, University of Minnesota, Minneapolis, Minnesota 55455, USA
| | - A Dominguez
- Bremen Center for Computational Materials Science, University of Bremen, Bremen, Germany
| | - S Ehlert
- University of Bonn, Bonn, Germany
| | - M Elstner
- Institute of Physical Chemistry, Karlsruhe Institute of Technology, Karlsruhe, Germany
| | - T van der Heide
- Bremen Center for Computational Materials Science, University of Bremen, Bremen, Germany
| | - J Hermann
- Freie Universität Berlin, Berlin, Germany
| | - S Irle
- Computational Sciences and Engineering Division, Oak Ridge National Laboratory, Oak Ridge, Tennessee 37831, USA
| | - J J Kranz
- Institute of Physical Chemistry, Karlsruhe Institute of Technology, Karlsruhe, Germany
| | - C Köhler
- Bremen Center for Computational Materials Science, University of Bremen, Bremen, Germany
| | - T Kowalczyk
- Department of Chemistry and Advanced Materials Science and Engineering Center, Western Washington University, Bellingham, Washington 98225, USA
| | - T Kubař
- Institute of Physical Chemistry, Karlsruhe Institute of Technology, Karlsruhe, Germany
| | - I S Lee
- Department of Chemistry, Ulsan National Institute of Science and Technology, Ulsan, South Korea
| | - V Lutsker
- Institut I - Theoretische Physik, University of Regensburg, Regensburg, Germany
| | - R J Maurer
- Department of Chemistry, University of Warwick, Gibbet Hill Road, Coventry CV4 7AL, United Kingdom
| | - S K Min
- Department of Chemistry, Ulsan National Institute of Science and Technology, Ulsan, South Korea
| | - I Mitchell
- Center for Multidimensional Carbon Materials, Institute for Basic Science (IBS), Ulsan 44919, South Korea
| | - C Negre
- Theoretical Division, Los Alamos National Laboratory, Los Alamos, New Mexico 87545, USA
| | - T A Niehaus
- Université de Lyon, Université Claude Bernard Lyon 1, CNRS, Institut Lumière Matière, F-69622 Villeurbanne, France
| | - A M N Niklasson
- Theoretical Division, Los Alamos National Laboratory, Los Alamos, New Mexico 87545, USA
| | - A J Page
- School of Environmental and Life Sciences, University of Newcastle, Callaghan, Australia
| | - A Pecchia
- CNR-ISMN, Via Salaria km 29.300, 00015 Monterotondo Stazione, Rome, Italy
| | - G Penazzi
- Bremen Center for Computational Materials Science, University of Bremen, Bremen, Germany
| | - M P Persson
- Dassault Systemes, Cambridge, United Kingdom
| | - J Řezáč
- Institute of Organic Chemistry and Biochemistry AS CR, Prague, Czech Republic
| | - C G Sánchez
- Instituto Interdisciplinario de Ciencias Básicas, Universidad Nacional de Cuyo, CONICET, Facultad de Ciencias Exactas y Naturales, Mendoza, Argentina
| | - M Sternberg
- Argonne National Laboratory, Lemont, Illinois 60439, USA
| | - M Stöhr
- Department of Physics and Materials Science, University of Luxembourg, L-1511 Luxembourg City, Luxembourg
| | - F Stuckenberg
- Bremen Center for Computational Materials Science, University of Bremen, Bremen, Germany
| | - A Tkatchenko
- Department of Physics and Materials Science, University of Luxembourg, L-1511 Luxembourg City, Luxembourg
| | - V W-Z Yu
- Department of Mechanical Engineering and Materials Science, Duke University, Durham, North Carolina 27708, USA
| | - T Frauenheim
- Bremen Center for Computational Materials Science, University of Bremen, Bremen, Germany
| |
Collapse
|
38
|
Structural Features of Triethylammonium Acetate through Molecular Dynamics. Molecules 2020; 25:molecules25061432. [PMID: 32245229 PMCID: PMC7146455 DOI: 10.3390/molecules25061432] [Citation(s) in RCA: 9] [Impact Index Per Article: 2.3] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/04/2020] [Accepted: 03/19/2020] [Indexed: 12/13/2022] Open
Abstract
I have explored the structural features and the dynamics of triethylammonium acetate by means of semi-empirical (density functional tight binding, DFTB) molecular dynamics. I find that the results from the present simulations agree with recent experimental determinations with only few minor differences in the structural interpretation. A mixture of triethylamine and acetic acid does not form an ionic liquid, but gives rise to a very complex system where ionization is only a partial process affecting only few molecules (1 over 4 experimentally). I have also found that the few ionic couples are stable and remain mainly embedded inside the AcOH neutral moiety.
Collapse
|
39
|
Řezáč J. Non-Covalent Interactions Atlas Benchmark Data Sets: Hydrogen Bonding. J Chem Theory Comput 2020; 16:2355-2368. [DOI: 10.1021/acs.jctc.9b01265] [Citation(s) in RCA: 30] [Impact Index Per Article: 7.5] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/14/2022]
Affiliation(s)
- Jan Řezáč
- Institute of Organic Chemistry and Biochemistry, Czech Academy of Sciences, 166 10 Prague, Czech Republic
| |
Collapse
|
40
|
Kříž K, Řezáč J. Benchmarking of Semiempirical Quantum-Mechanical Methods on Systems Relevant to Computer-Aided Drug Design. J Chem Inf Model 2020; 60:1453-1460. [PMID: 32062970 DOI: 10.1021/acs.jcim.9b01171] [Citation(s) in RCA: 35] [Impact Index Per Article: 8.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/21/2023]
Abstract
The semiempirical quantum mechanical (SQM) methods used in drug design are commonly parametrized and tested on data sets of systems that may not be representative models for drug-biomolecule interactions in terms of both size and chemical composition. This is addressed here with a new benchmark data set, PLF547, derived from protein-ligand complexes, consisting of complexes of ligands with protein fragments (such as amino-acid side chains), with interaction energies based on MP2-F12 and DLPNO-CCSD(T) calculations. From these, composite benchmark interaction energies are also built for complexes of the ligand with the complete active site of the protein (PLA15 data set). These data sets are used to test multiple SQM methods with corrections for noncovalent interactions; the role of the solvation model in the calculations is tested as well.
Collapse
Affiliation(s)
- Kristian Kříž
- Institute of Organic Chemistry and Biochemistry, Czech Academy of Sciences, 166 10 Prague, Czech Republic.,Department of Physical and Macromolecular Chemistry, Faculty of Science, Charles University, Hlavova 8, 128 40 Praha 2, Czech Republic
| | - Jan Řezáč
- Institute of Organic Chemistry and Biochemistry, Czech Academy of Sciences, 166 10 Prague, Czech Republic
| |
Collapse
|
41
|
Homochiral Supramolecular Thin Film from Self-Assembly of Achiral Triarylamine Molecules by Circularly Polarized Light. Molecules 2020; 25:molecules25020402. [PMID: 31963685 PMCID: PMC7024168 DOI: 10.3390/molecules25020402] [Citation(s) in RCA: 7] [Impact Index Per Article: 1.8] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/18/2019] [Revised: 01/12/2020] [Accepted: 01/16/2020] [Indexed: 01/19/2023] Open
Abstract
Here, we report the formation of homochiral supramolecular thin film from achiral molecules, by using circularly polarized light (CPL) only as a chiral source, on the condition that irradiation of CPL does not induce a photochemical change of the achiral molecules. Thin films of self-assembled structures consisting of chiral supramolecular fibrils was obtained from the triarylamine derivatives through evaporation of the self-assembled triarylamine solution. The homochiral supramolecular helices with the desired handedness was achieved by irradiation of circularly polarized visible light during the self-assembly process, and the chiral stability of supramolecular self-assembled product was achieved by photopolymerization of the diacetylene moieties at side chains of the building blocks, with irradiation of circularly polarized ultraviolet light. This work provides a novel methodology for the generation of homochiral supramolecular thin film from the corresponding achiral molecules.
Collapse
|
42
|
Sakti AW, Nishimura Y, Nakai H. Recent advances in quantum‐mechanical molecular dynamics simulations of proton transfer mechanism in various water‐based environments. WIRES COMPUTATIONAL MOLECULAR SCIENCE 2020. [DOI: 10.1002/wcms.1419] [Citation(s) in RCA: 5] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 11/07/2022]
Affiliation(s)
- Aditya W. Sakti
- Element Strategy Initiative for Catalysts and Batteries (ESICB) Kyoto University Kyoto Japan
| | - Yoshifumi Nishimura
- Waseda Research Institute for Science and Engineering (WISE) Waseda University Tokyo Japan
| | - Hiromi Nakai
- Element Strategy Initiative for Catalysts and Batteries (ESICB) Kyoto University Kyoto Japan
- Waseda Research Institute for Science and Engineering (WISE) Waseda University Tokyo Japan
- Department of Chemistry and Biochemistry, School of Advanced Science and Engineering Waseda University Tokyo Japan
| |
Collapse
|
43
|
Eyrilmez SM, Köprülüoğlu C, Řezáč J, Hobza P. Impressive Enrichment of Semiempirical Quantum Mechanics-Based Scoring Function: HSP90 Protein with 4541 Inhibitors and Decoys. Chemphyschem 2019; 20:2759-2766. [PMID: 31460692 DOI: 10.1002/cphc.201900628] [Citation(s) in RCA: 8] [Impact Index Per Article: 1.6] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/26/2019] [Revised: 08/21/2019] [Indexed: 12/11/2022]
Abstract
This paper describes the excellent performance of a newly developed scoring function (SF), based on the semiempirical QM (SQM) PM6-D3H4X method combined with the conductor-like screening implicit solvent model (COSMO). The SQM/COSMO, Amber/GB and nine widely used SFs have been evaluated in terms of ranking power on the HSP90 protein with 72 biologically active compounds and 4469 structurally similar decoys. Among conventional SFs, the highest early and overall enrichment measured by EF1 and AUC% obtained using single-scoring-function ranking has been found for Glide SP and Gold-ASP SFs, respectively (7, 75 % and 3, 76 %). The performance of other standard SFs has not been satisfactory, mostly even decreasing below random values. The SQM/COSMO SF, where P-L structures were optimised at the advanced Amber level, has resulted in a dramatic enrichment increase (47, 98 %), almost reaching the best possible receiver operator characteristic (ROC) curve. The best SQM frame thus inserts about seven times more active compounds into the selected dataset than the best standard SF.
Collapse
Affiliation(s)
- Saltuk M Eyrilmez
- Institute of Organic Chemistry and Biochemistry of the Czech Academy of Sciences, Flemingovo nam. 2, 16610, Prague 6, Czech Republic.,Regional Centre of Advanced Technologies and Materials, Palacký University, 77146, Olomouc, CzechRepublic
| | - Cemal Köprülüoğlu
- Institute of Organic Chemistry and Biochemistry of the Czech Academy of Sciences, Flemingovo nam. 2, 16610, Prague 6, Czech Republic.,Regional Centre of Advanced Technologies and Materials, Palacký University, 77146, Olomouc, CzechRepublic
| | - Jan Řezáč
- Institute of Organic Chemistry and Biochemistry of the Czech Academy of Sciences, Flemingovo nam. 2, 16610, Prague 6, Czech Republic
| | - Pavel Hobza
- Institute of Organic Chemistry and Biochemistry of the Czech Academy of Sciences, Flemingovo nam. 2, 16610, Prague 6, Czech Republic.,Regional Centre of Advanced Technologies and Materials, Palacký University, 77146, Olomouc, CzechRepublic
| |
Collapse
|
44
|
Lu X, Duchimaza-Heredia J, Cui Q. Analysis of Density Functional Tight Binding with Natural Bonding Orbitals. J Phys Chem A 2019; 123:7439-7453. [PMID: 31373822 DOI: 10.1021/acs.jpca.9b05072] [Citation(s) in RCA: 4] [Impact Index Per Article: 0.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/29/2022]
Abstract
The description of chemical bonding by the density functional tight binding (DFTB) model is analyzed using natural bonding orbitals (NBOs) and compared to results from density functional theory (B3LYP/aug-cc-pVTZ) calculations. Several molecular systems have been chosen to represent fairly diverse bonding scenarios that include standard covalent bonds, hypervalent interactions, multicenter bonds, metal-ligand interactions (with and without the pseudo-Jahn-Teller effect), and through-space donor-acceptor interactions. Overall, the results suggest that DFTB3/3OB provides physically sound descriptions for the different bonding scenarios analyzed here, as reflected by the general agreement between DFTB3 and B3LYP NBO properties, such as the nature of the NBOs, the magnitudes of natural charges and bond orders, and the dominant donor-acceptor interactions. The degree of ligand-to-metal charge transfer and the ionic nature of pentavalent phosphate are overestimated, likely reflecting the minimal-basis nature of DFTB3/3OB. Moreover, certain orbital interactions, such as geminal interactions, are observed to be grossly overestimated by DFTB3 for hypervalent phosphate and several transition metal compounds that involve copper and nickel. The study indicates that results from NBO analysis can be instructive for identifying electronic structure descriptions at the approximate quantum-mechanical level that require improvement and thus for guiding the systematic improvement of these methods.
Collapse
Affiliation(s)
- Xiya Lu
- Department of Chemistry and Theoretical Chemistry Institute , University of Wisconsin-Madison , 1101 University Avenue , Madison , Wisconsin 53706 , United States
| | | | | |
Collapse
|
45
|
Řezáč J. Description of halogen bonding in semiempirical quantum-mechanical and self-consistent charge density-functional tight-binding methods. J Comput Chem 2019; 40:1633-1642. [PMID: 30941801 DOI: 10.1002/jcc.25816] [Citation(s) in RCA: 6] [Impact Index Per Article: 1.2] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/22/2018] [Revised: 02/14/2019] [Accepted: 02/21/2019] [Indexed: 12/20/2022]
Abstract
This article analyzes the ability of semiempirical quantum-mechanical methods (PM6 and PM7) and self-consistent charge density-functional tight-binding (SCC-DFTB) method DFTB3 to describe halogen bonds. Calculations of the electrostatic potential on the surface of molecules containing halogens show that the σ-hole could be described well in modified neglect of diatomic overlap-based methods. The situation is more complex in the case of DFTB3 where a simpler model is used for the electrostatics, but short-ranged effects are covered in the Hamiltonian. All these methods can thus capture the effects that, for example, define the geometry of halogen bonds. The interaction energies are, however, affected by generally underestimated repulsion, which has been addressed earlier by standalone empirical corrections. Another approach to correcting this issue in DFTB3 is presented here-a modification of the energies of d-orbitals on halogens yields better results than the empirical correction in DFTB3-D3X, although it remains difficult to describe halogen and hydrogen bonds simultaneously. © 2019 Wiley Periodicals, Inc.
Collapse
Affiliation(s)
- Jan Řezáč
- Institute of Organic Chemistry and Biochemistry, Czech Academy of Sciences, 166 10, Prague, Czech Republic
| |
Collapse
|
46
|
Nishimura Y, Nakai H. D
cdftbmd
: Divide‐and‐Conquer Density Functional Tight‐Binding Program for Huge‐System Quantum Mechanical Molecular Dynamics Simulations. J Comput Chem 2019; 40:1538-1549. [DOI: 10.1002/jcc.25804] [Citation(s) in RCA: 36] [Impact Index Per Article: 7.2] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/04/2018] [Revised: 01/31/2019] [Accepted: 02/05/2019] [Indexed: 12/17/2022]
Affiliation(s)
- Yoshifumi Nishimura
- Waseda Research Institute for Science and Engineering, Waseda University 3‐4‐1 Okubo, Shinjuku‐ku, Tokyo 169‐8555 Japan
| | - Hiromi Nakai
- Waseda Research Institute for Science and Engineering, Waseda University 3‐4‐1 Okubo, Shinjuku‐ku, Tokyo 169‐8555 Japan
- Department of Chemistry and BiochemistrySchool of Advanced Science and Engineering, Waseda University 3‐4‐1 Okubo, Shinjuku‐ku, Tokyo 169‐8555 Japan
- ESICB, Kyoto University Kyotodaigaku‐Katsura, Kyoto 615‐8520 Japan
| |
Collapse
|
47
|
Pokorná P, Kruse H, Krepl M, Šponer J. QM/MM Calculations on Protein-RNA Complexes: Understanding Limitations of Classical MD Simulations and Search for Reliable Cost-Effective QM Methods. J Chem Theory Comput 2018; 14:5419-5433. [PMID: 30199638 DOI: 10.1021/acs.jctc.8b00670] [Citation(s) in RCA: 28] [Impact Index Per Article: 4.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/15/2022]
Abstract
Although atomistic explicit-solvent Molecular Dynamics (MD) is a popular tool to study protein-RNA recognition, satisfactory MD description of protein-RNA complexes is not always achieved. Unfortunately, it is often difficult to separate MD simulation instabilities primarily caused by the simple point-charge molecular mechanics (MM) force fields from problems related to the notorious uncertainties in the starting structures. Herein, we report a series of large-scale QM/MM calculations on the U1A protein-RNA complex. This experimentally well-characterized system has an intricate protein-RNA interface, which is very unstable in MD simulations. The QM/MM calculations identify several H-bonds poorly described by the MM method and thus indicate the sources of instabilities of the U1A interface in MD simulations. The results suggest that advanced QM/MM computations could be used to indirectly rationalize problems seen in MM-based MD simulations of protein-RNA complexes. As the most accurate QM method, we employ the computationally demanding meta-GGA density functional TPSS-D3(BJ)/def2-TZVP level of theory. Because considerably faster methods would be needed to extend sampling and to study even larger protein-RNA interfaces, a set of low-cost QM/MM methods is compared to the TPSS-D3(BJ)/def2-TZVP data. The PBEh-3c and B97-3c density functional composite methods appear to be suitable for protein-RNA interfaces. In contrast, HF-3c and the tight-binding Hamiltonians DFTB3-D3 and GFN-xTB perform unsatisfactorily and do not provide any advantage over the MM description. These conclusions are supported also by similar analysis of a simple HutP protein-RNA interface, which is well-described by MD with the exception of just one H-bond. Some other methodological aspects of QM/MM calculations on protein-RNA interfaces are discussed.
Collapse
Affiliation(s)
- Pavlína Pokorná
- Institute of Biophysics of the Czech Academy of Sciences , Královopolská 135 , Brno 612 65 , Czech Republic
| | - Holger Kruse
- Institute of Biophysics of the Czech Academy of Sciences , Královopolská 135 , Brno 612 65 , Czech Republic
| | - Miroslav Krepl
- Institute of Biophysics of the Czech Academy of Sciences , Královopolská 135 , Brno 612 65 , Czech Republic
| | - Jiří Šponer
- Institute of Biophysics of the Czech Academy of Sciences , Královopolská 135 , Brno 612 65 , Czech Republic.,CEITEC - Central European Institute of Technology, Masaryk University , Campus Bohunice, Kamenice 5 , Brno 625 00 , Czech Republic
| |
Collapse
|
48
|
Miriyala VM, Řezáč J. Testing Semiempirical Quantum Mechanical Methods on a Data Set of Interaction Energies Mapping Repulsive Contacts in Organic Molecules. J Phys Chem A 2018; 122:2801-2808. [PMID: 29473742 DOI: 10.1021/acs.jpca.8b00260] [Citation(s) in RCA: 15] [Impact Index Per Article: 2.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Abstract
Semiempirical quantum mechanical (QM) methods with corrections for noncovalent interactions provide a favorable combination of accuracy and computational efficiency that makes them a useful tool for a study of large molecular systems. It was, however, noted that the accuracy of these methods deteriorates at intermolecular distances shorter than equilibrium. In this work, we explore this issue systematically using a newly developed data set of benchmark interaction energies named R160×6. This data set maps repulsive contacts in organic molecules, and it consists of 160 model complexes for which six points along the dissociation curve are provided. Testing a wide range of semiempirical QM methods against the CCSD(T)/CBS benchmark revealed that most methods, and all the dispersion-corrected ones, underestimate the repulsion systematically. The worst cases are usually hydrogen-hydrogen contacts. The best results were obtained with PM6-D3H4 and DFTB3-D3H4, as these methods already contain a correction for the H-H repulsion, but the errors are still about twice as large as in equilibrium geometries.
Collapse
Affiliation(s)
- V M Miriyala
- Department of Computational Chemistry , Institute of Organic Chemistry and Biochemistry, Czech Academy of Sciences , Flemingovo Náměstí 542/2 , 16610 Prague , Czech Republic
| | - J Řezáč
- Department of Computational Chemistry , Institute of Organic Chemistry and Biochemistry, Czech Academy of Sciences , Flemingovo Náměstí 542/2 , 16610 Prague , Czech Republic
| |
Collapse
|
49
|
Řezáč J, Bím D, Gutten O, Rulíšek L. Toward Accurate Conformational Energies of Smaller Peptides and Medium-Sized Macrocycles: MPCONF196 Benchmark Energy Data Set. J Chem Theory Comput 2018; 14:1254-1266. [DOI: 10.1021/acs.jctc.7b01074] [Citation(s) in RCA: 55] [Impact Index Per Article: 9.2] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/23/2023]
Affiliation(s)
- Jan Řezáč
- Gilead Sciences Research Center and The Institute of Organic Chemistry and Biochemistry, Czech Academy of Sciences, Flemingovo náměstí 2, 166 10 Praha 6, Czech Republic
| | - Daniel Bím
- Gilead Sciences Research Center and The Institute of Organic Chemistry and Biochemistry, Czech Academy of Sciences, Flemingovo náměstí 2, 166 10 Praha 6, Czech Republic
| | - Ondrej Gutten
- Gilead Sciences Research Center and The Institute of Organic Chemistry and Biochemistry, Czech Academy of Sciences, Flemingovo náměstí 2, 166 10 Praha 6, Czech Republic
| | - Lubomír Rulíšek
- Gilead Sciences Research Center and The Institute of Organic Chemistry and Biochemistry, Czech Academy of Sciences, Flemingovo náměstí 2, 166 10 Praha 6, Czech Republic
| |
Collapse
|
50
|
Mortazavi M, Brandenburg JG, Maurer RJ, Tkatchenko A. Structure and Stability of Molecular Crystals with Many-Body Dispersion-Inclusive Density Functional Tight Binding. J Phys Chem Lett 2018; 9:399-405. [PMID: 29298075 DOI: 10.1021/acs.jpclett.7b03234] [Citation(s) in RCA: 28] [Impact Index Per Article: 4.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 06/07/2023]
Abstract
Accurate prediction of structure and stability of molecular crystals is crucial in materials science and requires reliable modeling of long-range dispersion interactions. Semiempirical electronic structure methods are computationally more efficient than their ab initio counterparts, allowing structure sampling with significant speedups. We combine the Tkatchenko-Scheffler van der Waals method (TS) and the many-body dispersion method (MBD) with third-order density functional tight-binding (DFTB3) via a charge population-based method. We find an overall good performance for the X23 benchmark database of molecular crystals, despite an underestimation of crystal volume that can be traced to the DFTB parametrization. We achieve accurate lattice energy predictions with DFT+MBD energetics on top of vdW-inclusive DFTB3 structures, resulting in a speedup of up to 3000 times compared with a full DFT treatment. This suggests that vdW-inclusive DFTB3 can serve as a viable structural prescreening tool in crystal structure prediction.
Collapse
Affiliation(s)
- Majid Mortazavi
- Fritz-Haber-Institut der Max-Planck-Gesellschaft , Faradayweg 4-6, 14195 Berlin, Germany
| | - Jan Gerit Brandenburg
- Department of Chemistry, University College London , 20 Gordon Street, WC1H 0AJ London, United Kingdom
- London Centre for Nanotechnology, University College London , 17-19 Gordon Street, WC1H 0AJ London, United Kingdom
- Thomas Young Centre, University College London , Gower Street, WC1E 6BT London, United Kingdom
| | - Reinhard J Maurer
- Department of Chemistry and Centre for Scientific Computing, University of Warwick , Gibbet Hill Road, Coventry CV4 7AL, United Kingdom
| | - Alexandre Tkatchenko
- Fritz-Haber-Institut der Max-Planck-Gesellschaft , Faradayweg 4-6, 14195 Berlin, Germany
- Physics and Materials Science Research Unit, University of Luxembourg , L-1511, Luxembourg
| |
Collapse
|