1
|
Thürlemann M, Riniker S. Hybrid classical/machine-learning force fields for the accurate description of molecular condensed-phase systems. Chem Sci 2023; 14:12661-12675. [PMID: 38020395 PMCID: PMC10646964 DOI: 10.1039/d3sc04317g] [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] [Received: 08/17/2023] [Accepted: 10/24/2023] [Indexed: 12/01/2023] Open
Abstract
Electronic structure methods offer in principle accurate predictions of molecular properties, however, their applicability is limited by computational costs. Empirical methods are cheaper, but come with inherent approximations and are dependent on the quality and quantity of training data. The rise of machine learning (ML) force fields (FFs) exacerbates limitations related to training data even further, especially for condensed-phase systems for which the generation of large and high-quality training datasets is difficult. Here, we propose a hybrid ML/classical FF model that is parametrized exclusively on high-quality ab initio data of dimers and monomers in vacuum but is transferable to condensed-phase systems. The proposed hybrid model combines our previous ML-parametrized classical model with ML corrections for situations where classical approximations break down, thus combining the robustness and efficiency of classical FFs with the flexibility of ML. Extensive validation on benchmarking datasets and experimental condensed-phase data, including organic liquids and small-molecule crystal structures, showcases how the proposed approach may promote FF development and unlock the full potential of classical FFs.
Collapse
Affiliation(s)
- Moritz Thürlemann
- Department of Chemistry and Applied Biosciences, ETH Zürich Vladimir-Prelog-Weg 2 Zürich 8093 Switzerland
| | - Sereina Riniker
- Department of Chemistry and Applied Biosciences, ETH Zürich Vladimir-Prelog-Weg 2 Zürich 8093 Switzerland
| |
Collapse
|
2
|
Han Y, Luo H, Lu Q, Liu Z, Liu J, Zhang J, Wei Z, Li J. Quantum Mechanical-Based Stability Evaluation of Crystal Structures for HIV-Targeted Drug Cabotegravir. Molecules 2021; 26:molecules26237178. [PMID: 34885762 PMCID: PMC8659202 DOI: 10.3390/molecules26237178] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/09/2021] [Revised: 11/21/2021] [Accepted: 11/23/2021] [Indexed: 11/23/2022] Open
Abstract
The long-acting parenteral formulation of the HIV integrase inhibitor cabotegravir (GSK744) is currently being developed to prevent HIV infections, benefiting from infrequent dosing and high efficacy. The crystal structure can affect the bioavailability and efficacy of cabotegravir. However, the stability determination of crystal structures of GSK744 have remained a challenge. Here, we introduced an ab initio protocol to determine the stability of the crystal structures of pharmaceutical molecules, which were obtained from crystal structure prediction process starting from the molecular diagram. Using GSK744 as a case study, the ab initio predicted that Gibbs free energy provides reliable further refinement of the predicted crystal structures and presents its capability for becoming a crystal stability determination approach in the future. The proposed work can assist in the comprehensive screening of pharmaceutical design and can provide structural predictions and stability evaluation for pharmaceutical crystals.
Collapse
Affiliation(s)
- Yanqiang Han
- Shanghai Key Laboratory of Maternal Fetal Medicine, Shanghai First Maternity and Infant Hospital, School of Medicine, Tongji University, Shanghai 200092, China; (Y.H.); (Z.L.)
- Key Laboratory for Thin Film and Microfabrication of Ministry of Education, Department of Micro/Nano-electronics, Shanghai Jiao Tong University, Shanghai 200240, China; (H.L.); (Q.L.)
| | - Hongyuan Luo
- Key Laboratory for Thin Film and Microfabrication of Ministry of Education, Department of Micro/Nano-electronics, Shanghai Jiao Tong University, Shanghai 200240, China; (H.L.); (Q.L.)
| | - Qianqian Lu
- Key Laboratory for Thin Film and Microfabrication of Ministry of Education, Department of Micro/Nano-electronics, Shanghai Jiao Tong University, Shanghai 200240, China; (H.L.); (Q.L.)
| | - Zeying Liu
- Shanghai Key Laboratory of Maternal Fetal Medicine, Shanghai First Maternity and Infant Hospital, School of Medicine, Tongji University, Shanghai 200092, China; (Y.H.); (Z.L.)
| | - Jinyun Liu
- Key Laboratory of Functional Molecular Solids of the Ministry of Education, Anhui Provincial Engineering Laboratory for New-Energy Vehicle Battery Energy-Storage Materials, School of Chemistry and Materials Science, Anhui Normal University, Wuhu 241002, China
- Correspondence: (J.L.); (Z.W.); (J.L.)
| | - Jiarui Zhang
- Division of Computational Biomedicine, Boston University School of Medicine, Boston, MA 02118, USA;
| | - Zhiyun Wei
- Shanghai Key Laboratory of Maternal Fetal Medicine, Shanghai First Maternity and Infant Hospital, School of Medicine, Tongji University, Shanghai 200092, China; (Y.H.); (Z.L.)
- Correspondence: (J.L.); (Z.W.); (J.L.)
| | - Jinjin Li
- Shanghai Key Laboratory of Maternal Fetal Medicine, Shanghai First Maternity and Infant Hospital, School of Medicine, Tongji University, Shanghai 200092, China; (Y.H.); (Z.L.)
- Key Laboratory for Thin Film and Microfabrication of Ministry of Education, Department of Micro/Nano-electronics, Shanghai Jiao Tong University, Shanghai 200240, China; (H.L.); (Q.L.)
- Correspondence: (J.L.); (Z.W.); (J.L.)
| |
Collapse
|
3
|
Bowskill DH, Sugden IJ, Konstantinopoulos S, Adjiman CS, Pantelides CC. Crystal Structure Prediction Methods for Organic Molecules: State of the Art. Annu Rev Chem Biomol Eng 2021; 12:593-623. [PMID: 33770462 DOI: 10.1146/annurev-chembioeng-060718-030256] [Citation(s) in RCA: 20] [Impact Index Per Article: 6.7] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/09/2022]
Abstract
The prediction of the crystal structures that a given organic molecule is likely to form is an important theoretical problem of significant interest for the pharmaceutical and agrochemical industries, among others. As evidenced by a series of six blind tests organized over the past 2 decades, methodologies for crystal structure prediction (CSP) have witnessed substantial progress and have now reached a stage of development where they can begin to be applied to systems of practical significance. This article reviews the state of the art in general-purpose methodologies for CSP, placing them within a common framework that highlights both their similarities and their differences. The review discusses specific areas that constitute the main focus of current research efforts toward improving the reliability and widening applicability of these methodologies, and offers some perspectives for the evolution of this technology over the next decade.
Collapse
Affiliation(s)
- David H Bowskill
- Molecular Systems Engineering Group, Centre for Process Systems Engineering, Department of Chemical Engineering, and Institute for Molecular Science and Engineering, Imperial College London, South Kensington Campus, London SW7 2AZ, United Kingdom;
| | - Isaac J Sugden
- Molecular Systems Engineering Group, Centre for Process Systems Engineering, Department of Chemical Engineering, and Institute for Molecular Science and Engineering, Imperial College London, South Kensington Campus, London SW7 2AZ, United Kingdom;
| | - Stefanos Konstantinopoulos
- Molecular Systems Engineering Group, Centre for Process Systems Engineering, Department of Chemical Engineering, and Institute for Molecular Science and Engineering, Imperial College London, South Kensington Campus, London SW7 2AZ, United Kingdom;
| | - Claire S Adjiman
- Molecular Systems Engineering Group, Centre for Process Systems Engineering, Department of Chemical Engineering, and Institute for Molecular Science and Engineering, Imperial College London, South Kensington Campus, London SW7 2AZ, United Kingdom;
| | - Constantinos C Pantelides
- Molecular Systems Engineering Group, Centre for Process Systems Engineering, Department of Chemical Engineering, and Institute for Molecular Science and Engineering, Imperial College London, South Kensington Campus, London SW7 2AZ, United Kingdom;
| |
Collapse
|
4
|
Greenwell C, McKinley JL, Zhang P, Zeng Q, Sun G, Li B, Wen S, Beran GJO. Overcoming the difficulties of predicting conformational polymorph energetics in molecular crystals via correlated wavefunction methods. Chem Sci 2020; 11:2200-2214. [PMID: 32190277 PMCID: PMC7059316 DOI: 10.1039/c9sc05689k] [Citation(s) in RCA: 38] [Impact Index Per Article: 9.5] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/09/2019] [Accepted: 01/13/2020] [Indexed: 11/21/2022] Open
Abstract
Molecular crystal structure prediction is increasingly being applied to study the solid form landscapes of larger, more flexible pharmaceutical molecules. Despite many successes in crystal structure prediction, van der Waals-inclusive density functional theory (DFT) methods exhibit serious failures predicting the polymorph stabilities for a number of systems exhibiting conformational polymorphism, where changes in intramolecular conformation lead to different intermolecular crystal packings. Here, the stabilities of the conformational polymorphs of o-acetamidobenzamide, ROY, and oxalyl dihydrazide are examined in detail. DFT functionals that have previously been very successful in crystal structure prediction perform poorly in all three systems, due primarily to the poor intramolecular conformational energies, but also due to the intermolecular description in oxalyl dihydrazide. In all three cases, a fragment-based dispersion-corrected second-order Møller-Plesset perturbation theory (MP2D) treatment of the crystals overcomes these difficulties and predicts conformational polymorph stabilities in good agreement with experiment. These results highlight the need for methods which go beyond current-generation DFT functionals to make crystal polymorph stability predictions truly reliable.
Collapse
Affiliation(s)
- Chandler Greenwell
- Department of Chemistry , University of California , Riverside , California 92521 , USA . ; Tel: +1-951-827-7869
| | - Jessica L McKinley
- Department of Chemistry , University of California , Riverside , California 92521 , USA . ; Tel: +1-951-827-7869
| | - Peiyu Zhang
- Xtalpi, Inc. , 245 Main St, 12th Floor , Cambridge , MA 02142 , USA
| | - Qun Zeng
- Xtalpi, Inc. , 245 Main St, 12th Floor , Cambridge , MA 02142 , USA
| | - Guangxu Sun
- Xtalpi, Inc. , 245 Main St, 12th Floor , Cambridge , MA 02142 , USA
| | - Bochen Li
- Xtalpi, Inc. , 245 Main St, 12th Floor , Cambridge , MA 02142 , USA
| | - Shuhao Wen
- Xtalpi, Inc. , 245 Main St, 12th Floor , Cambridge , MA 02142 , USA
| | - Gregory J O Beran
- Department of Chemistry , University of California , Riverside , California 92521 , USA . ; Tel: +1-951-827-7869
| |
Collapse
|
5
|
Toader AM, Zarić SD, Zalaru CM, Ferbinteanu M. The Structural Details of Aspirin Molecules and Crystals. Curr Med Chem 2020; 27:99-120. [PMID: 30381068 DOI: 10.2174/0929867325666181031132823] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/17/2017] [Revised: 01/29/2018] [Accepted: 02/06/2018] [Indexed: 11/22/2022]
Abstract
We revisit, in the key of structural chemistry, one of the most known and important drugs: the aspirin. Although apparently simple, the factors determining the molecular structure and supramolecular association in crystals are not trivial. We addressed the problem from experimental and theoretical sides, considering issues from X-ray measurements and results of first-principle reconstruction of molecule and lattices by ab initio calculations. Some puzzling problems can give headaches to specialists and intrigue the general public. Thus, the reported polymorphism of aspirin is disputed, a so-called form II being alleged as a result of misinterpretation. At the same time, were presented evidences that the structure of common form I can be disrupted by domains where the regular packing is changed to the pattern of form II. The problems appear even at the level of independent molecule: the most stable conformation computed by various techniques of electronic structure differs from those encountered in crystals. Because the energy difference between the related conformational isomers (computed as most stable vs. the experimental structure) is small, about 1 kcal/mol, comprised in the error bars of used methods, the unresting question is whether the modelling is imprecise, or the supramolecular factors are mutating the conformational preferences. By a detective following of the issue, the intermolecular effects were made responsible for the conformation of the molecule in crystal. The presented problems were gathered from literature results, debates, glued with modelling and analysis redone by ourselves, in order to secure the unitary view of the considered prototypic topic.
Collapse
Affiliation(s)
| | - Snezana D Zarić
- Department of Chemistry, University of Belgrade, Belgrade, Serbia.,Department of Chemistry, Texas A&M University at Qatar, Doha, Qatar
| | - Christina M Zalaru
- Department of Organic Chemistry, Biochemistry and Catalysis, Faculty of Chemistry, University of Bucharest, Bucharest, Romania
| | - Marilena Ferbinteanu
- Department of Inorganic Chemistry, Faculty of Chemistry, University of Bucharest, Bucharest, Romania
| |
Collapse
|
6
|
Olsen R, Kvamme B, Kuznetsova T. Hydrogen bond lifetimes and statistics of aqueous mono-, di- and tri-ethylene glycol. AIChE J 2016. [DOI: 10.1002/aic.15539] [Citation(s) in RCA: 13] [Impact Index Per Article: 1.6] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/01/2023]
Affiliation(s)
- Richard Olsen
- Dept. of Physics and Technology; University of Bergen; Allégaten 55 Bergen 5007 Norway
| | - Bjørn Kvamme
- Dept. of Physics and Technology; University of Bergen; Allégaten 55 Bergen 5007 Norway
| | - Tatiana Kuznetsova
- Dept. of Physics and Technology; University of Bergen; Allégaten 55 Bergen 5007 Norway
| |
Collapse
|
7
|
Elking DM, Fusti-Molnar L, Nichols A. Crystal structure prediction of rigid molecules. ACTA CRYSTALLOGRAPHICA SECTION B-STRUCTURAL SCIENCE CRYSTAL ENGINEERING AND MATERIALS 2016; 72:488-501. [DOI: 10.1107/s2052520616010118] [Citation(s) in RCA: 5] [Impact Index Per Article: 0.6] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Received: 02/15/2016] [Accepted: 06/21/2016] [Indexed: 11/11/2022]
Abstract
A non-polarizable force field based on atomic multipoles fit to reproduce experimental crystal properties andab initiogas-phase dimers is described. The Ewald method is used to calculate both long-range electrostatic and 1/r6dispersion energies of crystals. The dispersion energy of a crystal calculated by a cutoff method is shown to converge slowly to the exact Ewald result. A method for constraining space-group symmetry during unit-cell optimization is derived. Results for locally optimizing 4427 unit cells including volume, cell parameters, unit-cell r.m.s.d. and CPU timings are given for both flexible and rigid molecule optimization. An algorithm for randomly generating rigid molecule crystals is described. Using the correct experimentally determined space group, the average and maximum number of random crystals needed to find the correct experimental structure is given for 2440 rigid single component crystals. The force field energy rank of the correct experimental structure is presented for the same set of 2440 rigid single component crystals assuming the correct space group. A complete crystal prediction is performed for two rigid molecules by searching over the 32 most probable space groups.
Collapse
|
8
|
Reilly AM, Cooper RI, Adjiman CS, Bhattacharya S, Boese AD, Brandenburg JG, Bygrave PJ, Bylsma R, Campbell JE, Car R, Case DH, Chadha R, Cole JC, Cosburn K, Cuppen HM, Curtis F, Day GM, DiStasio Jr RA, Dzyabchenko A, van Eijck BP, Elking DM, van den Ende JA, Facelli JC, Ferraro MB, Fusti-Molnar L, Gatsiou CA, Gee TS, de Gelder R, Ghiringhelli LM, Goto H, Grimme S, Guo R, Hofmann DWM, Hoja J, Hylton RK, Iuzzolino L, Jankiewicz W, de Jong DT, Kendrick J, de Klerk NJJ, Ko HY, Kuleshova LN, Li X, Lohani S, Leusen FJJ, Lund AM, Lv J, Ma Y, Marom N, Masunov AE, McCabe P, McMahon DP, Meekes H, Metz MP, Misquitta AJ, Mohamed S, Monserrat B, Needs RJ, Neumann MA, Nyman J, Obata S, Oberhofer H, Oganov AR, Orendt AM, Pagola GI, Pantelides CC, Pickard CJ, Podeszwa R, Price LS, Price SL, Pulido A, Read MG, Reuter K, Schneider E, Schober C, Shields GP, Singh P, Sugden IJ, Szalewicz K, Taylor CR, Tkatchenko A, Tuckerman ME, Vacarro F, Vasileiadis M, Vazquez-Mayagoitia A, Vogt L, Wang Y, Watson RE, de Wijs GA, Yang J, Zhu Q, Groom CR. Report on the sixth blind test of organic crystal structure prediction methods. ACTA CRYSTALLOGRAPHICA SECTION B, STRUCTURAL SCIENCE, CRYSTAL ENGINEERING AND MATERIALS 2016; 72:439-59. [PMID: 27484368 PMCID: PMC4971545 DOI: 10.1107/s2052520616007447] [Citation(s) in RCA: 308] [Impact Index Per Article: 38.5] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Subscribe] [Scholar Register] [Received: 02/15/2016] [Accepted: 05/04/2016] [Indexed: 05/05/2023]
Abstract
The sixth blind test of organic crystal structure prediction (CSP) methods has been held, with five target systems: a small nearly rigid molecule, a polymorphic former drug candidate, a chloride salt hydrate, a co-crystal and a bulky flexible molecule. This blind test has seen substantial growth in the number of participants, with the broad range of prediction methods giving a unique insight into the state of the art in the field. Significant progress has been seen in treating flexible molecules, usage of hierarchical approaches to ranking structures, the application of density-functional approximations, and the establishment of new workflows and `best practices' for performing CSP calculations. All of the targets, apart from a single potentially disordered Z' = 2 polymorph of the drug candidate, were predicted by at least one submission. Despite many remaining challenges, it is clear that CSP methods are becoming more applicable to a wider range of real systems, including salts, hydrates and larger flexible molecules. The results also highlight the potential for CSP calculations to complement and augment experimental studies of organic solid forms.
Collapse
Affiliation(s)
- Anthony M. Reilly
- The Cambridge Crystallographic Data Centre, 12 Union Road, Cambridge CB2 1EZ, England
| | - Richard I. Cooper
- Chemical Crystallography, Chemistry Research Laboratory, Mansfield Road, Oxford OX1 3TA, England
| | - Claire S. Adjiman
- Department of Chemical Engineering, Centre for Process Systems Engineering, Imperial College London, London SW7 2AZ, England
| | - Saswata Bhattacharya
- Fritz-Haber-Institut der Max-Planck-Gesellschaft, Faradayweg 4-6, 14195, Berlin, Germany
| | - A. Daniel Boese
- Department of Chemistry, Institute of Physical and Theoretical Chemistry, University of Graz, Heinrichstraße 28/IV, 8010 Graz, Austria
| | - Jan Gerit Brandenburg
- Mulliken Center for Theoretical Chemistry, Institut für Physikalische und Theoretische Chemie, Rheinische Friedrich-Wilhelms Universität Bonn, Beringstraße 4, 53115 Bonn, Germany
| | - Peter J. Bygrave
- School of Chemistry, University of Southampton, Southampton SO17 1BJ, England
| | - Rita Bylsma
- Radboud University, Institute for Molecules and Materials, Heyendaalseweg 135, 6525 AJ Nijmegen, The Netherlands
| | - Josh E. Campbell
- School of Chemistry, University of Southampton, Southampton SO17 1BJ, England
| | - Roberto Car
- Department of Chemistry, Princeton University, Princeton, NJ 08544, USA
| | - David H. Case
- School of Chemistry, University of Southampton, Southampton SO17 1BJ, England
| | - Renu Chadha
- University Institute of Pharmaceutical Sciences, Panjab University, Chandigarh, India
| | - Jason C. Cole
- The Cambridge Crystallographic Data Centre, 12 Union Road, Cambridge CB2 1EZ, England
| | - Katherine Cosburn
- Department of Physics and Engineering Physics, Tulane University, New Orleans, LA 70118, USA
- Department of Physics, University of Toronto, Toronto, Canada M5S 1A7
| | - Herma M. Cuppen
- Radboud University, Institute for Molecules and Materials, Heyendaalseweg 135, 6525 AJ Nijmegen, The Netherlands
| | - Farren Curtis
- Department of Physics and Engineering Physics, Tulane University, New Orleans, LA 70118, USA
- Department of Physics, Carnegie Mellon University, Pittsburgh, PA 15213, USA
| | - Graeme M. Day
- School of Chemistry, University of Southampton, Southampton SO17 1BJ, England
| | - Robert A. DiStasio Jr
- Department of Chemistry, Princeton University, Princeton, NJ 08544, USA
- Department of Chemistry and Chemical Biology, Cornell University, Ithaca, NY 14853, USA
| | | | | | - Dennis M. Elking
- OpenEye Scientific Software, 9 Bisbee Court, Suite D, Santa Fe, NM 87508, USA
| | - Joost A. van den Ende
- Radboud University, Institute for Molecules and Materials, Heyendaalseweg 135, 6525 AJ Nijmegen, The Netherlands
| | - Julio C. Facelli
- Center for High Performance Computing, University of Utah, 155 South 1452 East Room 405, Salt Lake City, UT 84112-0190, USA
- Department of Biomedical Informatics, University of Utah, 155 South 1452 East Room 405, Salt Lake City, UT 84112-0190, USA
| | - Marta B. Ferraro
- Departamento de Física and Ifiba (CONICET) Facultad de Ciencias Exactas y Naturales, Universidad de Buenos Aires, Ciudad Universitaria, Pab. I (1428), Buenos Aires, Argentina
| | - Laszlo Fusti-Molnar
- OpenEye Scientific Software, 9 Bisbee Court, Suite D, Santa Fe, NM 87508, USA
| | - Christina-Anna Gatsiou
- Department of Chemical Engineering, Centre for Process Systems Engineering, Imperial College London, London SW7 2AZ, England
| | - Thomas S. Gee
- School of Chemistry, University of Southampton, Southampton SO17 1BJ, England
| | - René de Gelder
- Radboud University, Institute for Molecules and Materials, Heyendaalseweg 135, 6525 AJ Nijmegen, The Netherlands
| | - Luca M. Ghiringhelli
- Fritz-Haber-Institut der Max-Planck-Gesellschaft, Faradayweg 4-6, 14195, Berlin, Germany
| | - Hitoshi Goto
- Educational Programs on Advanced Simulation Engineering, Toyohashi University of Technology, 1-1 Hibarigaoka, Tempaku-cho, Toyohashi, Aichi 441-8580, Japan
- Department of Computer Science and Engineering, Graduate School of Engineering, Toyohashi University of Technology, 1-1 Hibarigaoka, Tempaku-cho, Toyohashi, Aichi 441-8580, Japan
| | - Stefan Grimme
- Mulliken Center for Theoretical Chemistry, Institut für Physikalische und Theoretische Chemie, Rheinische Friedrich-Wilhelms Universität Bonn, Beringstraße 4, 53115 Bonn, Germany
| | - Rui Guo
- Department of Chemistry, University College London, 20 Gordon Street, London WC1H 0AJ, England
| | - Detlef W. M. Hofmann
- CRS4, Parco Scientifico e Tecnologico, POLARIS, Edificio 1, 09010 PULA, Italy
- FlexCryst, Schleifweg 23, 91080 Uttenreuth, Germany
| | - Johannes Hoja
- Fritz-Haber-Institut der Max-Planck-Gesellschaft, Faradayweg 4-6, 14195, Berlin, Germany
| | - Rebecca K. Hylton
- Department of Chemistry, University College London, 20 Gordon Street, London WC1H 0AJ, England
| | - Luca Iuzzolino
- Department of Chemistry, University College London, 20 Gordon Street, London WC1H 0AJ, England
| | - Wojciech Jankiewicz
- Institute of Chemistry, University of Silesia, Szkolna 9, 40-006 Katowice, Poland
| | - Daniël T. de Jong
- Radboud University, Institute for Molecules and Materials, Heyendaalseweg 135, 6525 AJ Nijmegen, The Netherlands
| | - John Kendrick
- Faculty of Life Sciences, University of Bradford, Richmond Road, Bradford BD7 1DP, England
| | - Niek J. J. de Klerk
- Radboud University, Institute for Molecules and Materials, Heyendaalseweg 135, 6525 AJ Nijmegen, The Netherlands
| | - Hsin-Yu Ko
- Department of Chemistry, Princeton University, Princeton, NJ 08544, USA
| | | | - Xiayue Li
- Department of Physics and Engineering Physics, Tulane University, New Orleans, LA 70118, USA
- Argonne Leadership Computing Facility, Argonne National Laboratory, Lemont, IL 60439, USA
| | - Sanjaya Lohani
- Department of Physics and Engineering Physics, Tulane University, New Orleans, LA 70118, USA
| | - Frank J. J. Leusen
- Faculty of Life Sciences, University of Bradford, Richmond Road, Bradford BD7 1DP, England
| | - Albert M. Lund
- OpenEye Scientific Software, 9 Bisbee Court, Suite D, Santa Fe, NM 87508, USA
- Department of Chemistry, University of Utah, 155 South 1452 East Room 405, Salt Lake City, UT 84112-0190, USA
| | - Jian Lv
- State Key Laboratory of Superhard Materials, Jilin University, Changchun 130012, People’s Republic of China
| | - Yanming Ma
- State Key Laboratory of Superhard Materials, Jilin University, Changchun 130012, People’s Republic of China
| | - Noa Marom
- Department of Physics and Engineering Physics, Tulane University, New Orleans, LA 70118, USA
- Department of Materials Science and Engineering and Department of Physics, Carnegie Mellon University, Pittsburgh, PA 15213, USA
| | - Artëm E. Masunov
- NanoScience Technology Center, University of Central Florida, 12424 Research Parkway, PAV400, Orlando, FL 32826, USA
- Department of Chemistry, University of Central Florida, 4111 Libra Drive PSB225, Orlando, FL 32816, USA
- Department of Physics, University of Central Florida, 4111 Libra Drive PSB430, Orlando, FL 32816, USA
- Department of Condensed Matter Physics, National Research Nuclear University MEPhI, Kashirskoye shosse 31, Moscow 115409, Russia
| | - Patrick McCabe
- The Cambridge Crystallographic Data Centre, 12 Union Road, Cambridge CB2 1EZ, England
| | - David P. McMahon
- School of Chemistry, University of Southampton, Southampton SO17 1BJ, England
| | - Hugo Meekes
- Radboud University, Institute for Molecules and Materials, Heyendaalseweg 135, 6525 AJ Nijmegen, The Netherlands
| | - Michael P. Metz
- Department of Physics and Astronomy, University of Delaware, Newark, DE 19716, USA
| | - Alston J. Misquitta
- School of Physics and Astronomy, Queen Mary University of London, London E1 4NS, England
| | | | - Bartomeu Monserrat
- Cavendish Laboratory, 19, J. J. Thomson Avenue, Cambridge CB3 0HE, England
- Department of Physics and Astronomy, Rutgers University, Piscataway, NJ 08854-8019, USA
| | - Richard J. Needs
- Cavendish Laboratory, 19, J. J. Thomson Avenue, Cambridge CB3 0HE, England
| | | | - Jonas Nyman
- School of Chemistry, University of Southampton, Southampton SO17 1BJ, England
| | - Shigeaki Obata
- Educational Programs on Advanced Simulation Engineering, Toyohashi University of Technology, 1-1 Hibarigaoka, Tempaku-cho, Toyohashi, Aichi 441-8580, Japan
| | - Harald Oberhofer
- Chair for Theoretical Chemistry and Catalysis Research Center, Technische Universität München, Lichtenbergstr. 4, D-85747 Garching, Germany
| | - Artem R. Oganov
- Department of Geosciences, Center for Materials by Design, and Institute for Advanced Computational Science, SUNY Stony Brook, NY 11794-2100, USA
- Skolkovo Institute of Science and Technology, Skolkovo Innovation Centers, Bldg. 3, Moscow Region, 143026, Russia
- Moscow Institute of Physics and Technology, 9 Institutskiy Lane, Dolgoprudny City, Moscow Region 141700, Russia
- International Center for Materials Discovery, School of Materials Science and Engineering, Northwestern Polytechnical University, Xi’an 710072, China
| | - Anita M. Orendt
- Center for High Performance Computing, University of Utah, 155 South 1452 East Room 405, Salt Lake City, UT 84112-0190, USA
| | - Gabriel I. Pagola
- Departamento de Física and Ifiba (CONICET) Facultad de Ciencias Exactas y Naturales, Universidad de Buenos Aires, Ciudad Universitaria, Pab. I (1428), Buenos Aires, Argentina
| | - Constantinos C. Pantelides
- Department of Chemical Engineering, Centre for Process Systems Engineering, Imperial College London, London SW7 2AZ, England
| | - Chris J. Pickard
- Department of Materials Science and Metallurgy, University of Cambridge, 27 Charles Babbage Road, Cambridge CB3 0FS, England
- Department of Physics and Astronomy, University College London, Gower St., London WC1E 6BT, England
| | - Rafal Podeszwa
- Institute of Chemistry, University of Silesia, Szkolna 9, 40-006 Katowice, Poland
| | - Louise S. Price
- Department of Chemistry, University College London, 20 Gordon Street, London WC1H 0AJ, England
| | - Sarah L. Price
- Department of Chemistry, University College London, 20 Gordon Street, London WC1H 0AJ, England
| | - Angeles Pulido
- School of Chemistry, University of Southampton, Southampton SO17 1BJ, England
| | - Murray G. Read
- The Cambridge Crystallographic Data Centre, 12 Union Road, Cambridge CB2 1EZ, England
| | - Karsten Reuter
- Chair for Theoretical Chemistry and Catalysis Research Center, Technische Universität München, Lichtenbergstr. 4, D-85747 Garching, Germany
| | - Elia Schneider
- Department of Chemistry, New York University, New York, NY 10003, USA
| | - Christoph Schober
- Chair for Theoretical Chemistry and Catalysis Research Center, Technische Universität München, Lichtenbergstr. 4, D-85747 Garching, Germany
| | - Gregory P. Shields
- The Cambridge Crystallographic Data Centre, 12 Union Road, Cambridge CB2 1EZ, England
| | - Pawanpreet Singh
- University Institute of Pharmaceutical Sciences, Panjab University, Chandigarh, India
| | - Isaac J. Sugden
- Department of Chemical Engineering, Centre for Process Systems Engineering, Imperial College London, London SW7 2AZ, England
| | - Krzysztof Szalewicz
- Department of Physics and Astronomy, University of Delaware, Newark, DE 19716, USA
| | | | - 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
| | - Mark E. Tuckerman
- Department of Chemistry, New York University, New York, NY 10003, USA
- Courant Institute of Mathematical Sciences, New York University, New York, NY 10012, USA
- NYU-ECNU Center for Computational Chemistry at NYU Shanghai, 3663 Zhongshan Road North, Shanghai 200062, China
| | - Francesca Vacarro
- Department of Physics and Engineering Physics, Tulane University, New Orleans, LA 70118, USA
- Department of Chemistry, Loyola University, New Orleans, LA 70118, USA
| | - Manolis Vasileiadis
- Department of Chemical Engineering, Centre for Process Systems Engineering, Imperial College London, London SW7 2AZ, England
| | | | - Leslie Vogt
- Department of Chemistry, New York University, New York, NY 10003, USA
| | - Yanchao Wang
- State Key Laboratory of Superhard Materials, Jilin University, Changchun 130012, People’s Republic of China
| | - Rona E. Watson
- Department of Chemistry, University College London, 20 Gordon Street, London WC1H 0AJ, England
| | - Gilles A. de Wijs
- Radboud University, Institute for Molecules and Materials, Heyendaalseweg 135, 6525 AJ Nijmegen, The Netherlands
| | - Jack Yang
- School of Chemistry, University of Southampton, Southampton SO17 1BJ, England
| | - Qiang Zhu
- Department of Geosciences, Center for Materials by Design, and Institute for Advanced Computational Science, SUNY Stony Brook, NY 11794-2100, USA
| | - Colin R. Groom
- The Cambridge Crystallographic Data Centre, 12 Union Road, Cambridge CB2 1EZ, England
| |
Collapse
|
9
|
Abstract
Interest in molecular crystals has grown thanks to their relevance to pharmaceuticals, organic semiconductor materials, foods, and many other applications. Electronic structure methods have become an increasingly important tool for modeling molecular crystals and polymorphism. This article reviews electronic structure techniques used to model molecular crystals, including periodic density functional theory, periodic second-order Møller-Plesset perturbation theory, fragment-based electronic structure methods, and diffusion Monte Carlo. It also discusses the use of these models for predicting a variety of crystal properties that are relevant to the study of polymorphism, including lattice energies, structures, crystal structure prediction, polymorphism, phase diagrams, vibrational spectroscopies, and nuclear magnetic resonance spectroscopy. Finally, tools for analyzing crystal structures and intermolecular interactions are briefly discussed.
Collapse
Affiliation(s)
- Gregory J O Beran
- Department of Chemistry, University of California , Riverside, California 92521, United States
| |
Collapse
|
10
|
Habgood M, Sugden IJ, Kazantsev AV, Adjiman CS, Pantelides CC. Efficient Handling of Molecular Flexibility in Ab Initio Generation of Crystal Structures. J Chem Theory Comput 2015; 11:1957-69. [DOI: 10.1021/ct500621v] [Citation(s) in RCA: 50] [Impact Index Per Article: 5.6] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/01/2023]
Affiliation(s)
- Matthew Habgood
- Molecular Systems Engineering
Group, Centre for Process Systems Engineering, Department of Chemical
Engineering, Imperial College London, London SW7 2AZ, United Kingdom
| | - Isaac J. Sugden
- Molecular Systems Engineering
Group, Centre for Process Systems Engineering, Department of Chemical
Engineering, Imperial College London, London SW7 2AZ, United Kingdom
| | - Andrei V. Kazantsev
- Molecular Systems Engineering
Group, Centre for Process Systems Engineering, Department of Chemical
Engineering, Imperial College London, London SW7 2AZ, United Kingdom
| | - Claire S. Adjiman
- Molecular Systems Engineering
Group, Centre for Process Systems Engineering, Department of Chemical
Engineering, Imperial College London, London SW7 2AZ, United Kingdom
| | - Constantinos C. Pantelides
- Molecular Systems Engineering
Group, Centre for Process Systems Engineering, Department of Chemical
Engineering, Imperial College London, London SW7 2AZ, United Kingdom
| |
Collapse
|
11
|
Abramov YA. Current Computational Approaches to Support Pharmaceutical Solid Form Selection. Org Process Res Dev 2012. [DOI: 10.1021/op300274s] [Citation(s) in RCA: 67] [Impact Index Per Article: 5.6] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Affiliation(s)
- Yuriy A. Abramov
- Pfizer Global Research and Development, Groton, Connecticut 06340, USA
| |
Collapse
|
12
|
|
13
|
Bardwell DA, Adjiman CS, Arnautova YA, Bartashevich E, Boerrigter SXM, Braun DE, Cruz-Cabeza AJ, Day GM, Della Valle RG, Desiraju GR, van Eijck BP, Facelli JC, Ferraro MB, Grillo D, Habgood M, Hofmann DWM, Hofmann F, Jose KVJ, Karamertzanis PG, Kazantsev AV, Kendrick J, Kuleshova LN, Leusen FJJ, Maleev AV, Misquitta AJ, Mohamed S, Needs RJ, Neumann MA, Nikylov D, Orendt AM, Pal R, Pantelides CC, Pickard CJ, Price LS, Price SL, Scheraga HA, van de Streek J, Thakur TS, Tiwari S, Venuti E, Zhitkov IK. Towards crystal structure prediction of complex organic compounds--a report on the fifth blind test. ACTA CRYSTALLOGRAPHICA. SECTION B, STRUCTURAL SCIENCE 2011; 67:535-51. [PMID: 22101543 PMCID: PMC3222142 DOI: 10.1107/s0108768111042868] [Citation(s) in RCA: 247] [Impact Index Per Article: 19.0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Subscribe] [Scholar Register] [Received: 08/01/2011] [Accepted: 10/16/2011] [Indexed: 12/04/2022]
Abstract
Following on from the success of the previous crystal structure prediction blind tests (CSP1999, CSP2001, CSP2004 and CSP2007), a fifth such collaborative project (CSP2010) was organized at the Cambridge Crystallographic Data Centre. A range of methodologies was used by the participating groups in order to evaluate the ability of the current computational methods to predict the crystal structures of the six organic molecules chosen as targets for this blind test. The first four targets, two rigid molecules, one semi-flexible molecule and a 1:1 salt, matched the criteria for the targets from CSP2007, while the last two targets belonged to two new challenging categories - a larger, much more flexible molecule and a hydrate with more than one polymorph. Each group submitted three predictions for each target it attempted. There was at least one successful prediction for each target, and two groups were able to successfully predict the structure of the large flexible molecule as their first place submission. The results show that while not as many groups successfully predicted the structures of the three smallest molecules as in CSP2007, there is now evidence that methodologies such as dispersion-corrected density functional theory (DFT-D) are able to reliably do so. The results also highlight the many challenges posed by more complex systems and show that there are still issues to be overcome.
Collapse
Affiliation(s)
- David A Bardwell
- Cambridge Crystallographic Data Centre, 12 Union Road, Cambridge CB2 1EZ, England.
| | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | |
Collapse
|
14
|
Kazantsev AV, Karamertzanis PG, Adjiman CS, Pantelides CC. Efficient Handling of Molecular Flexibility in Lattice Energy Minimization of Organic Crystals. J Chem Theory Comput 2011; 7:1998-2016. [DOI: 10.1021/ct100597e] [Citation(s) in RCA: 97] [Impact Index Per Article: 7.5] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Affiliation(s)
- A. V. Kazantsev
- Centre for Process Systems Engineering, Department of Chemical Engineering, Imperial College London, South Kensington Campus, London SW7 2AZ, United Kingdom
| | - P. G. Karamertzanis
- Centre for Process Systems Engineering, Department of Chemical Engineering, Imperial College London, South Kensington Campus, London SW7 2AZ, United Kingdom
| | - C. S. Adjiman
- Centre for Process Systems Engineering, Department of Chemical Engineering, Imperial College London, South Kensington Campus, London SW7 2AZ, United Kingdom
| | - C. C. Pantelides
- Centre for Process Systems Engineering, Department of Chemical Engineering, Imperial College London, South Kensington Campus, London SW7 2AZ, United Kingdom
| |
Collapse
|
15
|
|
16
|
Price SL, Leslie M, Welch GWA, Habgood M, Price LS, Karamertzanis PG, Day GM. Modelling organic crystal structures using distributed multipole and polarizability-based model intermolecular potentials. Phys Chem Chem Phys 2010; 12:8478-90. [PMID: 20607186 DOI: 10.1039/c004164e] [Citation(s) in RCA: 202] [Impact Index Per Article: 14.4] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/09/2023]
Abstract
Crystal structure prediction for organic molecules requires both the fast assessment of thousands to millions of crystal structures and the greatest possible accuracy in their relative energies. We describe a crystal lattice simulation program, DMACRYS, emphasizing the features that make it suitable for use in crystal structure prediction for pharmaceutical molecules using accurate anisotropic atom-atom model intermolecular potentials based on the theory of intermolecular forces. DMACRYS can optimize the lattice energy of a crystal, calculate the second derivative properties, and reduce the symmetry of the spacegroup to move away from a transition state. The calculated terahertz frequency k = 0 rigid-body lattice modes and elastic tensor can be used to estimate free energies. The program uses a distributed multipole electrostatic model (Q, t = 00,...,44s) for the electrostatic fields, and can use anisotropic atom-atom repulsion models, damped isotropic dispersion up to R(-10), as well as a range of empirically fitted isotropic exp-6 atom-atom models with different definitions of atomic types. A new feature is that an accurate model for the induction energy contribution to the lattice energy has been implemented that uses atomic anisotropic dipole polarizability models (alpha, t = (10,10)...(11c,11s)) to evaluate the changes in the molecular charge density induced by the electrostatic field within the crystal. It is demonstrated, using the four polymorphs of the pharmaceutical carbamazepine C(15)H(12)N(2)O, that whilst reproducing crystal structures is relatively easy, calculating the polymorphic energy differences to the accuracy of a few kJ mol(-1) required for applications is very demanding of assumptions made in the modelling. Thus DMACRYS enables the comparison of both known and hypothetical crystal structures as an aid to the development of pharmaceuticals and other speciality organic materials, and provides a tool to develop the modelling of the intermolecular forces involved in molecular recognition processes.
Collapse
Affiliation(s)
- Sarah L Price
- Department of Chemistry, 20 Gordon Street, London WC1H 0AJ, UK.
| | | | | | | | | | | | | |
Collapse
|
17
|
Day GM, Cooper TG. Crystal packing predictions of the alpha-amino acids: methods assessment and structural observations. CrystEngComm 2010. [DOI: 10.1039/c002213f] [Citation(s) in RCA: 31] [Impact Index Per Article: 2.2] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022]
|
18
|
Kim S, Orendt AM, Ferraro MB, Facelli JC. Crystal structure prediction of flexible molecules using parallel genetic algorithms with a standard force field. J Comput Chem 2009; 30:1973-85. [PMID: 19130496 DOI: 10.1002/jcc.21189] [Citation(s) in RCA: 21] [Impact Index Per Article: 1.4] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/06/2022]
Abstract
This article describes the application of our distributed computing framework for crystal structure prediction (CSP) the modified genetic algorithms for crystal and cluster prediction (MGAC), to predict the crystal structure of flexible molecules using the general Amber force field (GAFF) and the CHARMM program. The MGAC distributed computing framework includes a series of tightly integrated computer programs for generating the molecule's force field, sampling crystal structures using a distributed parallel genetic algorithm and local energy minimization of the structures followed by the classifying, sorting, and archiving of the most relevant structures. Our results indicate that the method can consistently find the experimentally known crystal structures of flexible molecules, but the number of missing structures and poor ranking observed in some crystals show the need for further improvement of the potential.
Collapse
Affiliation(s)
- Seonah Kim
- Center for High Performance Computing, University of Utah, 155 South 1452 East Room 405, Salt Lake City, Utah 84112-0190, USA
| | | | | | | |
Collapse
|
19
|
Karamertzanis PG, Kazantsev AV, Issa N, Welch GW, Adjiman CS, Pantelides CC, Price SL. Can the Formation of Pharmaceutical Cocrystals Be Computationally Predicted? 2. Crystal Structure Prediction. J Chem Theory Comput 2009; 5:1432-48. [DOI: 10.1021/ct8004326] [Citation(s) in RCA: 106] [Impact Index Per Article: 7.1] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Affiliation(s)
- Panagiotis G. Karamertzanis
- Centre for Process Systems Engineering, Department of Chemical Engineering, Imperial College London, SW7 2AZ, United Kingdom, and Department of Chemistry, University College London, 20 Gordon Street, London, WC1H 0AJ, United Kingdom
| | - Andrei V. Kazantsev
- Centre for Process Systems Engineering, Department of Chemical Engineering, Imperial College London, SW7 2AZ, United Kingdom, and Department of Chemistry, University College London, 20 Gordon Street, London, WC1H 0AJ, United Kingdom
| | - Nizar Issa
- Centre for Process Systems Engineering, Department of Chemical Engineering, Imperial College London, SW7 2AZ, United Kingdom, and Department of Chemistry, University College London, 20 Gordon Street, London, WC1H 0AJ, United Kingdom
| | - Gareth W.A. Welch
- Centre for Process Systems Engineering, Department of Chemical Engineering, Imperial College London, SW7 2AZ, United Kingdom, and Department of Chemistry, University College London, 20 Gordon Street, London, WC1H 0AJ, United Kingdom
| | - Claire S. Adjiman
- Centre for Process Systems Engineering, Department of Chemical Engineering, Imperial College London, SW7 2AZ, United Kingdom, and Department of Chemistry, University College London, 20 Gordon Street, London, WC1H 0AJ, United Kingdom
| | - Constantinos C. Pantelides
- Centre for Process Systems Engineering, Department of Chemical Engineering, Imperial College London, SW7 2AZ, United Kingdom, and Department of Chemistry, University College London, 20 Gordon Street, London, WC1H 0AJ, United Kingdom
| | - Sarah L. Price
- Centre for Process Systems Engineering, Department of Chemical Engineering, Imperial College London, SW7 2AZ, United Kingdom, and Department of Chemistry, University College London, 20 Gordon Street, London, WC1H 0AJ, United Kingdom
| |
Collapse
|
20
|
|
21
|
Day GM, S Motherwell WD, Jones W. A strategy for predicting the crystal structures of flexible molecules: the polymorphism of phenobarbital. Phys Chem Chem Phys 2007; 9:1693-704. [PMID: 17396181 DOI: 10.1039/b612190j] [Citation(s) in RCA: 68] [Impact Index Per Article: 4.0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022]
Abstract
A computational exploration of the low energy crystal structures of the pharmaceutical molecule phenobarbital is presented as a test of an approach for the crystal structure prediction of flexible molecules. Traditional transferable force field methods of modelling flexible molecules are unreliable for the level of accuracy required in crystal structure prediction and we outline a strategy for improving the evaluation of relative energies of large sets of crystal structures. The approach involves treating the molecule as a set of linked rigid units, whose conformational energy is expressed as a function of the relative orientations of the rigid groups. The conformational energy is calculated by electronic structure methods and the intermolecular interactions using an atomic multipole description of electrostatics. A key consideration in our approach is the scalability to more typical pharmaceutical molecules of higher molecular weight with many more atoms and degrees of flexibility. Based on our calculations, crystal structures are proposed for the as-yet uncharacterised forms IV and V, as well as further polymorphs of phenobarbital.
Collapse
Affiliation(s)
- G M Day
- Department of Chemistry, University of Cambridge, Lensfield Road, Cambridge, UK.
| | | | | |
Collapse
|
22
|
Karamertzanis† PG, Pantelides CC. Ab initiocrystal structure prediction. II. Flexible molecules. Mol Phys 2007. [DOI: 10.1080/00268970601143317] [Citation(s) in RCA: 61] [Impact Index Per Article: 3.6] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/23/2022]
|
23
|
|
24
|
Bazterra VE, Thorley M, Ferraro MB, Facelli JC. A Distributed Computing Method for Crystal Structure Prediction of Flexible Molecules: An Application to N-(2-Dimethyl-4,5-dinitrophenyl) Acetamide. J Chem Theory Comput 2006; 3:201-9. [DOI: 10.1021/ct6002115] [Citation(s) in RCA: 14] [Impact Index Per Article: 0.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
Affiliation(s)
- Victor E. Bazterra
- Departamento de Física, Facultad de Ciencias Exactas y Naturales, Universidad de Buenos Aires, Ciudad Universitaria, Pab. I (1428), Buenos Aires, Argentina, and Center for High Performance Computing, University of Utah, 155 South 1452 East Rm 405, Salt Lake City, Utah 84112-0190
| | - Matthew Thorley
- Departamento de Física, Facultad de Ciencias Exactas y Naturales, Universidad de Buenos Aires, Ciudad Universitaria, Pab. I (1428), Buenos Aires, Argentina, and Center for High Performance Computing, University of Utah, 155 South 1452 East Rm 405, Salt Lake City, Utah 84112-0190
| | - Marta B. Ferraro
- Departamento de Física, Facultad de Ciencias Exactas y Naturales, Universidad de Buenos Aires, Ciudad Universitaria, Pab. I (1428), Buenos Aires, Argentina, and Center for High Performance Computing, University of Utah, 155 South 1452 East Rm 405, Salt Lake City, Utah 84112-0190
| | - Julio C. Facelli
- Departamento de Física, Facultad de Ciencias Exactas y Naturales, Universidad de Buenos Aires, Ciudad Universitaria, Pab. I (1428), Buenos Aires, Argentina, and Center for High Performance Computing, University of Utah, 155 South 1452 East Rm 405, Salt Lake City, Utah 84112-0190
| |
Collapse
|
25
|
van Mourik T, Karamertzanis PG, Price SL. Molecular Conformations and Relative Stabilities Can Be as Demanding of the Electronic Structure Method as Intermolecular Calculations. J Phys Chem A 2005; 110:8-12. [PMID: 16392833 DOI: 10.1021/jp0563181] [Citation(s) in RCA: 108] [Impact Index Per Article: 5.7] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/18/2022]
Abstract
We have performed a variety of high-level electronic structure calculations on two moderately sized organic molecules and found considerable sensitivity of the intramolecular potential energy surface to the method employed. The gas-phase structure of tyrosine-glycine varies qualitatively between B3LYP and MP2 optimizations, producing different close contacts between the tyrosine ring and the glycine moiety. The relative energies of the 2-(acetylamino)benzamide conformations found in its two polymorphs can vary by over 20 kJ mol-1 between MP2 and B3LYP calculations, using the same basis set. It is shown by a novel analysis that the intramolecular equivalent of basis set superposition error competes with the errors in the intramolecular dispersion in causing this sensitivity.
Collapse
Affiliation(s)
- Tanja van Mourik
- Department of Chemistry, University College London, 20 Gordon Street, London WC1H 0AJ, UK.
| | | | | |
Collapse
|
26
|
Modelling Intermolecular Forces for Organic Crystal Structure Prediction. STRUCTURE AND BONDING 2005. [DOI: 10.1007/b135616] [Citation(s) in RCA: 7] [Impact Index Per Article: 0.4] [Reference Citation Analysis] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 03/07/2023]
|
27
|
Karamertzanis PG, Pantelides CC. Ab initio crystal structure prediction?I. Rigid molecules. J Comput Chem 2004; 26:304-24. [PMID: 15622548 DOI: 10.1002/jcc.20165] [Citation(s) in RCA: 108] [Impact Index Per Article: 5.4] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/11/2022]
Abstract
A new methodology for the prediction of molecular crystal structures using only the atomic connectivity of the molecule under consideration is presented. The approach is based on the global minimization of the lattice enthalpy of the crystal. The modeling of the electrostatic interactions is accomplished through a set of distributed charges that are optimally and automatically selected and positioned based on results of quantum mechanical calculations. A four-step global optimization algorithm is used for the identification of the local minima of the lattice enthalpy surface. A parallelized implementation of the algorithm permits a much more extensive search of the solution space than has hitherto been possible, allowing the identification of crystal structures in less frequently occurring space groups and with more than one molecule in the asymmetric unit. The algorithm has been applied successfully to the prediction of the crystal structures of 3-aza-bicyclo(3.3.1)nonane-2,4-dione (P2(1)/a, Z' = 1), allopurinol (P2(1)/c, Z' = 1), 1,3,4,6,7,9-hexa-azacycl(3.3.3)azine (Pbca, Z' = 2), and triethylenediamine (P6(3)/m, Z' = 1). In all cases, the experimentally known structure is among the most stable predicted structures, but not necessarily the global minimum.
Collapse
Affiliation(s)
- Panagiotis G Karamertzanis
- Centre for Process Systems Engineering, Department of Chemical Engineering and Chemical Technology, Imperial College London, South Kensington Campus, London SW7 2AZ, United Kingdom
| | | |
Collapse
|
28
|
Klein RA. Ab initio conformational studies on diols and binary diol-water systems using DFT methods. Intramolecular hydrogen bonding and 1:1 complex formation with water. J Comput Chem 2002; 23:585-99. [PMID: 11939593 DOI: 10.1002/jcc.10053] [Citation(s) in RCA: 95] [Impact Index Per Article: 4.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/10/2022]
Abstract
Studies on the conformational equilibrium for the following diols, ethane-1,2-diol (12EG, CAS 107-21-1), 2R-D-(-)-propane-1,2-diol (12PG, CAS 4254-14-2), (2S,3S)-L-(+)-butane-2,3-diol (L23BD, CAS 19132-06-0), and (2S,3R)-meso-butane-2,3-diol (m23BD, CAS 5341-95-7), are described using Gaussian ab initio calculations involving density functional theory (DFT) methods. We also report in this article results on the stability and conformation for the 1:1 water-diol complex formed by ethane-1,2-diol, propane-1,2-diol, and L- and meso-butane-2,3-diol. The relative stability of the intramolecular (internal) hydrogen bond in a range of diols (n = 2 to 6), based on ab initio geometry optimization and determination of the -O...H- distance, dOH, and -O-H...O- angle, theta, increases through the sequence 1,2 approximately equals 2,3 < 1,3 < 1,4 approximately equals 1,5 approximately equals 1,6, as judged from the bond linearity and -O...H- separation. Quantum mechanical and topological analysis of possible intramolecular hydrogen bonding in this complete series of diols provides convincing evidence for this in diols in which the hydroxyl groups are separated by three or more carbon atoms, that is, in (n, n+m) diols for m > or = 2, but not for ethane-1,2-diol or other vicinal diols, which do not satisfy Popelier's topological and electron density criteria based on the AIM theory of Bader. Based on these criteria it is unlikely that vicinal diols are in fact capable of forming an intramolecular hydrogen bond, in spite of geometric and spectroscopic data in the literature suggesting otherwise.
Collapse
Affiliation(s)
- Roger A Klein
- Institute for Physiological Chemistry, University of Bonn, Federal Republic of Germany.
| |
Collapse
|
29
|
Swerts B, Van Droogenbroeck J, Peeters A, Van Alsenoy C. A Linearly Scaling QM/MM Method to Study Molecular Crystals Using BRABO/CHARMM: Application to 2-(2-Methyl-3-chloroanilino) Nicotinic Acid. J Phys Chem A 2002. [DOI: 10.1021/jp014513j] [Citation(s) in RCA: 10] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Affiliation(s)
- Ben Swerts
- Structural Chemistry Group, Department of Chemistry, University of Antwerp, Universiteitsplein 1, B-2610 Wilrijk, Belgium
| | - Joris Van Droogenbroeck
- Structural Chemistry Group, Department of Chemistry, University of Antwerp, Universiteitsplein 1, B-2610 Wilrijk, Belgium
| | - Anik Peeters
- Structural Chemistry Group, Department of Chemistry, University of Antwerp, Universiteitsplein 1, B-2610 Wilrijk, Belgium
| | - Christian Van Alsenoy
- Structural Chemistry Group, Department of Chemistry, University of Antwerp, Universiteitsplein 1, B-2610 Wilrijk, Belgium
| |
Collapse
|
30
|
Gavezzotti A. Calculation of Intermolecular Interaction Energies by Direct Numerical Integration over Electron Densities. I. Electrostatic and Polarization Energies in Molecular Crystals. J Phys Chem B 2002. [DOI: 10.1021/jp0144202] [Citation(s) in RCA: 244] [Impact Index Per Article: 11.1] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
Affiliation(s)
- A. Gavezzotti
- Dipartimento di Chimica Strutturale e Stereochimica inorganica, Università di Milano, via Venezian 21, 20133 Milano, Italy
| |
Collapse
|
31
|
Van Eijck BP. Crystal structure predictions using five space groups with two independent molecules. The case of small organic acids. J Comput Chem 2002; 23:456-62. [PMID: 11908081 DOI: 10.1002/jcc.10042] [Citation(s) in RCA: 14] [Impact Index Per Article: 0.6] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/05/2022]
Abstract
Crystal structure generations with two independent molecules have been performed for a series of carboxylic acids, using a slightly modified version of the OPLS force field. It was found that in this way the experimental structures with one independent molecule were produced as special cases, except for the molecules with four or more internal degrees of freedom. This work shows that a search with two independent molecules in only five space groups, although costly in computer power, can automatically also find structures with one independent molecule in many supergroups. Considering the observed abundances of structural classes, such a search should cover more than 95% of the possible homomolecular crystal structures.
Collapse
Affiliation(s)
- Bouke P Van Eijck
- Department of Crystal and Structural Chemistry, Bijvoet Center for Biomolecular Research, Utrecht University, Padualaan 8, 3584 CH Utrecht, The Netherlands.
| |
Collapse
|
32
|
van Eijck BP, Mooij WTM, Kroon J. Crystal Structure Prediction for Six Monosaccharides Revisited. J Phys Chem B 2001. [DOI: 10.1021/jp012366j] [Citation(s) in RCA: 25] [Impact Index Per Article: 1.1] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Affiliation(s)
- Bouke P. van Eijck
- Department of Crystal and Structural Chemistry, Bijvoet Center for Biomolecular Research, Utrecht University, Padualaan 8, 3584 CH Utrecht, The Netherlands, and Accelrys Ltd., 230/250 The Quorum, Barnwell Road, Cambridge CB5 8RE, England
| | - Wijnand T. M. Mooij
- Department of Crystal and Structural Chemistry, Bijvoet Center for Biomolecular Research, Utrecht University, Padualaan 8, 3584 CH Utrecht, The Netherlands, and Accelrys Ltd., 230/250 The Quorum, Barnwell Road, Cambridge CB5 8RE, England
| | - Jan Kroon
- Department of Crystal and Structural Chemistry, Bijvoet Center for Biomolecular Research, Utrecht University, Padualaan 8, 3584 CH Utrecht, The Netherlands, and Accelrys Ltd., 230/250 The Quorum, Barnwell Road, Cambridge CB5 8RE, England
| |
Collapse
|
33
|
Beyer T, Lewis T, Price SL. Which organic crystal structures are predictable by lattice energy minimisation? CrystEngComm 2001. [DOI: 10.1039/b108135g] [Citation(s) in RCA: 52] [Impact Index Per Article: 2.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022]
|
34
|
van Eijck BP. Ab initio crystal structure predictions for flexible hydrogen-bonded molecules. Part III. Effect of lattice vibrations. J Comput Chem 2001. [DOI: 10.1002/jcc.1047] [Citation(s) in RCA: 54] [Impact Index Per Article: 2.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/09/2022]
|