1
|
Grelska J, Temleitner L, Park C, Jurkiewicz K, Pawlus S. High-Pressure and Temperature Effects on the Clustering Ability of Monohydroxy Alcohols. J Phys Chem Lett 2024; 15:3118-3126. [PMID: 38471115 PMCID: PMC10961836 DOI: 10.1021/acs.jpclett.4c00085] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/09/2024] [Revised: 02/20/2024] [Accepted: 02/22/2024] [Indexed: 03/14/2024]
Abstract
This study examined the clustering behavior of monohydroxy alcohols, where hydrogen-bonded clusters of up to a hundred molecules on the nanoscale can form. By performing X-ray diffraction experiments at different temperatures and under high pressure, we investigated how these conditions affect the ability of alcohols to form clusters. The pioneering high-pressure experiment performed on liquid alcohols contributes to the emerging knowledge in this field. Implementation of molecular dynamics simulations yielded excellent agreement with the experimental results, enabling the analysis of theoretical models. Here we show that at the same global density achieved either by alteration of pressure or temperature, the local aggregation of molecules at the nanoscale may significantly differ. Surprisingly, high pressure not only promotes the formation of hydrogen-bonded clusters but also induces the serious reorganization of molecules. This research represents a milestone in understanding association under extreme thermodynamic conditions in other hydrogen bonding systems such as water.
Collapse
Affiliation(s)
- Joanna Grelska
- A.
Chełkowski Institute of Physics, University
of Silesia in Katowice, 75 Pułku Piechoty 1, 41-500 Chorzów, Poland
| | - László Temleitner
- HUN-REN
Wigner Research Centre for Physics, Konkoly Thege út 29-33, H-1121 Budapest, Hungary
| | - Changyong Park
- High
Pressure Collaborative Access Team (HPCAT), X-ray Science Division, Argonne National Laboratory, Lemont, Illinois 60439, United States
| | - Karolina Jurkiewicz
- A.
Chełkowski Institute of Physics, University
of Silesia in Katowice, 75 Pułku Piechoty 1, 41-500 Chorzów, Poland
| | - Sebastian Pawlus
- A.
Chełkowski Institute of Physics, University
of Silesia in Katowice, 75 Pułku Piechoty 1, 41-500 Chorzów, Poland
| |
Collapse
|
2
|
Rahman M, Dannatt HRW, Blundell CD, Hughes LP, Blade H, Carson J, Tatman BP, Johnston ST, Brown SP. Polymorph Identification for Flexible Molecules: Linear Regression Analysis of Experimental and Calculated Solution- and Solid-State NMR Data. J Phys Chem A 2024; 128:1793-1816. [PMID: 38427685 PMCID: PMC10945485 DOI: 10.1021/acs.jpca.3c07732] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/24/2023] [Revised: 02/06/2024] [Accepted: 02/07/2024] [Indexed: 03/03/2024]
Abstract
The Δδ regression approach of Blade et al. [ J. Phys. Chem. A 2020, 124(43), 8959-8977] for accurately discriminating between solid forms using a combination of experimental solution- and solid-state NMR data with density functional theory (DFT) calculation is here extended to molecules with multiple conformational degrees of freedom, using furosemide polymorphs as an exemplar. As before, the differences in measured 1H and 13C chemical shifts between solution-state NMR and solid-state magic-angle spinning (MAS) NMR (Δδexperimental) are compared to those determined by gauge-including projector augmented wave (GIPAW) calculations (Δδcalculated) by regression analysis and a t-test, allowing the correct furosemide polymorph to be precisely identified. Monte Carlo random sampling is used to calculate solution-state NMR chemical shifts, reducing computation times by avoiding the need to systematically sample the multidimensional conformational landscape that furosemide occupies in solution. The solvent conditions should be chosen to match the molecule's charge state between the solution and solid states. The Δδ regression approach indicates whether or not correlations between Δδexperimental and Δδcalculated are statistically significant; the approach is differently sensitive to the popular root mean squared error (RMSE) method, being shown to exhibit a much greater dynamic range. An alternative method for estimating solution-state NMR chemical shifts by approximating the measured solution-state dynamic 3D behavior with an ensemble of 54 furosemide crystal structures (polymorphs and cocrystals) from the Cambridge Structural Database (CSD) was also successful in this case, suggesting new avenues for this method that may overcome its current dependency on the prior determination of solution dynamic 3D structures.
Collapse
Affiliation(s)
- Mohammed Rahman
- Department
of Physics, University of Warwick, Coventry CV4 7AL, U.K.
- Department
of Chemistry, University of Warwick, Coventry CV4 7AL, U.K.
| | | | | | - Leslie P. Hughes
- Oral
Product Development, Pharmaceutical Technology & Development, Operations, AstraZeneca, Macclesfield SK10 2NA, U.K.
| | - Helen Blade
- Oral
Product Development, Pharmaceutical Technology & Development, Operations, AstraZeneca, Macclesfield SK10 2NA, U.K.
| | - Jake Carson
- Mathematics
Institute at Warwick, University of Warwick, Coventry CV4 7AL, U.K.
| | - Ben P. Tatman
- Department
of Physics, University of Warwick, Coventry CV4 7AL, U.K.
- Department
of Chemistry, University of Warwick, Coventry CV4 7AL, U.K.
| | | | - Steven P. Brown
- Department
of Physics, University of Warwick, Coventry CV4 7AL, U.K.
| |
Collapse
|
3
|
Ye Z, Wang N, Zhou J, Ouyang D. Organic crystal structure prediction via coupled generative adversarial networks and graph convolutional networks. Innovation (N Y) 2024; 5:100562. [PMID: 38379785 PMCID: PMC10878116 DOI: 10.1016/j.xinn.2023.100562] [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: 07/19/2023] [Accepted: 12/29/2023] [Indexed: 02/22/2024] Open
Abstract
Organic crystal structures exert a profound impact on the physicochemical properties and biological effects of organic compounds. Quantum mechanics (QM)-based crystal structure predictions (CSPs) have somewhat alleviated the dilemma that experimental crystal structure investigations struggle to conduct complete polymorphism studies, but the high computing cost poses a challenge to its widespread application. The present study aims to construct DeepCSP, a feasible pure machine learning framework for minute-scale rapid organic CSP. Initially, based on 177,746 data entries from the Cambridge Crystal Structure Database, a generative adversarial network was built to conditionally generate trial crystal structures under selected feature constraints for the given molecule. Simultaneously, a graph convolutional attention network was used to predict the density of stable crystal structures for the input molecule. Subsequently, the distances between the predicted density and the definition-based calculated density would be considered to be the crystal structure screening and ranking basis, and finally, the density-based crystal structure ranking would be output. Two such distinct algorithms, performing the generation and ranking functionalities, respectively, collectively constitute the DeepCSP, which has demonstrated compelling performance in marketed drug validations, achieving an accuracy rate exceeding 80% and a hit rate surpassing 85%. Inspiringly, the computing speed of the pure machine learning methodology demonstrates the potential of artificial intelligence in advancing CSP research.
Collapse
Affiliation(s)
- Zhuyifan Ye
- State Key Laboratory of Quality Research in Chinese Medicine, Institute of Chinese Medical Sciences, University of Macau, Macau 999078, China
- Faculty of Applied Sciences, Macao Polytechnic University, Macao 999078, China
| | - Nannan Wang
- State Key Laboratory of Quality Research in Chinese Medicine, Institute of Chinese Medical Sciences, University of Macau, Macau 999078, China
| | - Jiantao Zhou
- State Key Laboratory of Internet of Things for Smart City, University of Macau, Macau 999078, China
- Department of Computer and Information Science, Faculty of Science and Technology, University of Macau, Macau 999078, China
| | - Defang Ouyang
- State Key Laboratory of Quality Research in Chinese Medicine, Institute of Chinese Medical Sciences, University of Macau, Macau 999078, China
- Department of Public Health and Medicinal Administration, Faculty of Health Sciences, University of Macau, Macau 999078, China
| |
Collapse
|
4
|
Beran GJO. Frontiers of molecular crystal structure prediction for pharmaceuticals and functional organic materials. Chem Sci 2023; 14:13290-13312. [PMID: 38033897 PMCID: PMC10685338 DOI: 10.1039/d3sc03903j] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/27/2023] [Accepted: 11/02/2023] [Indexed: 12/02/2023] Open
Abstract
The reliability of organic molecular crystal structure prediction has improved tremendously in recent years. Crystal structure predictions for small, mostly rigid molecules are quickly becoming routine. Structure predictions for larger, highly flexible molecules are more challenging, but their crystal structures can also now be predicted with increasing rates of success. These advances are ushering in a new era where crystal structure prediction drives the experimental discovery of new solid forms. After briefly discussing the computational methods that enable successful crystal structure prediction, this perspective presents case studies from the literature that demonstrate how state-of-the-art crystal structure prediction can transform how scientists approach problems involving the organic solid state. Applications to pharmaceuticals, porous organic materials, photomechanical crystals, organic semi-conductors, and nuclear magnetic resonance crystallography are included. Finally, efforts to improve our understanding of which predicted crystal structures can actually be produced experimentally and other outstanding challenges are discussed.
Collapse
Affiliation(s)
- Gregory J O Beran
- Department of Chemistry, University of California Riverside Riverside CA 92521 USA
| |
Collapse
|
5
|
Firaha D, Liu YM, van de Streek J, Sasikumar K, Dietrich H, Helfferich J, Aerts L, Braun DE, Broo A, DiPasquale AG, Lee AY, Le Meur S, Nilsson Lill SO, Lunsmann WJ, Mattei A, Muglia P, Putra OD, Raoui M, Reutzel-Edens SM, Rome S, Sheikh AY, Tkatchenko A, Woollam GR, Neumann MA. Predicting crystal form stability under real-world conditions. Nature 2023; 623:324-328. [PMID: 37938708 PMCID: PMC10632141 DOI: 10.1038/s41586-023-06587-3] [Citation(s) in RCA: 5] [Impact Index Per Article: 5.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/14/2022] [Accepted: 08/30/2023] [Indexed: 11/09/2023]
Abstract
The physicochemical properties of molecular crystals, such as solubility, stability, compactability, melting behaviour and bioavailability, depend on their crystal form1. In silico crystal form selection has recently come much closer to realization because of the development of accurate and affordable free-energy calculations2-4. Here we redefine the state of the art, primarily by improving the accuracy of free-energy calculations, constructing a reliable experimental benchmark for solid-solid free-energy differences, quantifying statistical errors for the computed free energies and placing both hydrate crystal structures of different stoichiometries and anhydrate crystal structures on the same energy landscape, with defined error bars, as a function of temperature and relative humidity. The calculated free energies have standard errors of 1-2 kJ mol-1 for industrially relevant compounds, and the method to place crystal structures with different hydrate stoichiometries on the same energy landscape can be extended to other multi-component systems, including solvates. These contributions reduce the gap between the needs of the experimentalist and the capabilities of modern computational tools, transforming crystal structure prediction into a more reliable and actionable procedure that can be used in combination with experimental evidence to direct crystal form selection and establish control5.
Collapse
Affiliation(s)
| | | | | | | | | | - Julian Helfferich
- Avant-garde Materials Simulation, Merzhausen, Germany
- JobRad, Freiburg, Germany
| | - Luc Aerts
- UCB Pharma SA, Chemin du Foriest, Braine-l'Alleud, Belgium
| | - Doris E Braun
- Institute of Pharmacy, University of Innsbruck, Innsbruck, Austria
| | - Anders Broo
- Data Science and Modelling, Pharmaceutical Sciences, R&D, AstraZeneca Gothenburg, Mölndal, Sweden
| | | | - Alfred Y Lee
- Merck, Analytical Research & Development, Rahway, NJ, USA
| | - Sarah Le Meur
- UCB Pharma SA, Chemin du Foriest, Braine-l'Alleud, Belgium
| | - Sten O Nilsson Lill
- Data Science and Modelling, Pharmaceutical Sciences, R&D, AstraZeneca Gothenburg, Mölndal, Sweden
| | | | - Alessandra Mattei
- Solid State Chemistry, Research & Development, AbbVie, North Chicago, IL, USA
| | | | - Okky Dwichandra Putra
- Early Product Development and Manufacturing, Pharmaceutical Sciences R&D, AstraZeneca Gothenburg, Mölndal, Sweden
| | | | - Susan M Reutzel-Edens
- Cambridge Crystallographic Data Centre, Cambridge, UK
- SuRE Pharma Consulting, Zionsville, IN, USA
| | - Sandrine Rome
- UCB Pharma SA, Chemin du Foriest, Braine-l'Alleud, Belgium
| | - Ahmad Y Sheikh
- Solid State Chemistry, Research & Development, AbbVie, North Chicago, IL, USA
| | - Alexandre Tkatchenko
- Department of Physics and Materials Science, University of Luxembourg, Luxembourg City, Luxembourg
| | | | | |
Collapse
|
6
|
Henaff C, Siepmann J, Siepmann F, Danède F, Avettand-Fènoël MN, Vérin J, Fadel A, Willart JF. Exploration of the physical states of riboflavin (free base) by mechanical milling. Int J Pharm 2023; 645:123416. [PMID: 37716485 DOI: 10.1016/j.ijpharm.2023.123416] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/06/2023] [Revised: 09/11/2023] [Accepted: 09/12/2023] [Indexed: 09/18/2023]
Abstract
Amorphous riboflavin (free base) could be produced for the first time via high energy ball milling of a commercial crystalline form (Form I). Importantly, this solid state amorphization process allowed to circumvent chemical degradation occurring during melting as well as the lack of suitable solvents, which are required for amorphization via spray- or freeze-drying. The amorphous state of riboflavin was thoroughly characterized, revealing a complex recrystallization pattern upon heating, involving two enantiotropic polymorphic forms (II and III) and a dihydrate. The glass transition temperature (Tg) and heat capacity (Cp) jump of the amorphous form were determined as 144 °C and 0.68 J/g/°C. Moreover, the relative physical stability of the different physical states has been elucidated, e.g., at room temperature: I > II > III. The following rank order was observed for the dissolution rates in water at 37 °C during the first 4 h: amorphous > III ≈ II > I. Afterwards, a dihydrate crystallized from the solutions of amorphous and metastable crystalline riboflavin forms, the solubility of which was well above the solubility of the stable FormI.
Collapse
Affiliation(s)
- C Henaff
- Univ. Lille, CNRS, INRAE, Centrale Lille, UMR 8207 UMET - Unité Matériaux et Transformations F-59000 Lille, France; Univ. Lille, Inserm, CHU Lille, U1008, F-59000 Lille, France
| | - J Siepmann
- Univ. Lille, Inserm, CHU Lille, U1008, F-59000 Lille, France
| | - F Siepmann
- Univ. Lille, Inserm, CHU Lille, U1008, F-59000 Lille, France
| | - F Danède
- Univ. Lille, CNRS, INRAE, Centrale Lille, UMR 8207 UMET - Unité Matériaux et Transformations F-59000 Lille, France
| | - M-N Avettand-Fènoël
- Univ. Lille, CNRS, INRAE, Centrale Lille, UMR 8207 UMET - Unité Matériaux et Transformations F-59000 Lille, France
| | - J Vérin
- Univ. Lille, Inserm, CHU Lille, U1008, F-59000 Lille, France
| | - A Fadel
- Univ. Lille, CNRS, INRAE, Centrale Lille, Univ. Artois, FR 2638 - IMEC - Institut Michel-Eugène Chevreul, F-59000 Lille, France
| | - J-F Willart
- Univ. Lille, CNRS, INRAE, Centrale Lille, UMR 8207 UMET - Unité Matériaux et Transformations F-59000 Lille, France.
| |
Collapse
|
7
|
Ward M, Taylor CR, Mulvee MT, Lampronti GI, Belenguer AM, Steed JW, Day GM, Oswald IDH. Pushing Technique Boundaries to Probe Conformational Polymorphism. CRYSTAL GROWTH & DESIGN 2023; 23:7217-7230. [PMID: 37808905 PMCID: PMC10557047 DOI: 10.1021/acs.cgd.3c00641] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Figures] [Subscribe] [Scholar Register] [Received: 05/26/2023] [Revised: 08/11/2023] [Indexed: 10/10/2023]
Abstract
We present an extensive exploration of the solid-form landscape of chlorpropamide (CPA) using a combined experimental-computational approach at the frontiers of both fields. We have obtained new conformational polymorphs of CPA, placing them into context with known forms using flexible-molecule crystal structure prediction. We highlight the formation of a new polymorph (ζ-CPA) via spray-drying experiments despite its notable metastability (14 kJ/mol) relative to the thermodynamic α-form, and we identify and resolve the ball-milled η-form isolated in 2019. Additionally, we employ impurity- and gel-assisted crystallization to control polymorphism and the formation of novel multicomponent forms. We, thus, demonstrate the power of this collaborative screening approach to observe, rationalize, and control the formation of new metastable forms.
Collapse
Affiliation(s)
- Martin
R. Ward
- Strathclyde
Institute of Pharmacy and Biomedical Sciences, University of Strathclyde, 161 Cathedral Street, Glasgow G4 0RE, U.K.
| | - Christopher R. Taylor
- Computational
Systems Chemistry, School of Chemistry, University of Southampton, Southampton SO17 1BJ, U.K.
| | - Matthew T. Mulvee
- Department
of Chemistry, Durham University, South Road, Durham DH1 3LE, U.K.
| | - Giulio I. Lampronti
- Department
of Materials Science & Metallurgy, University
of Cambridge, 27 Charles Babbage Rd, Cambridge CB3 0FS, U.K.
| | - Ana M. Belenguer
- Yusuf
Hamied Department of Chemistry, University
of Cambridge, Lensfield
Road, Cambridge CB2 1EW, U.K.
| | - Jonathan W. Steed
- Department
of Chemistry, Durham University, South Road, Durham DH1 3LE, U.K.
| | - Graeme M. Day
- Computational
Systems Chemistry, School of Chemistry, University of Southampton, Southampton SO17 1BJ, U.K.
| | - Iain D. H. Oswald
- Strathclyde
Institute of Pharmacy and Biomedical Sciences, University of Strathclyde, 161 Cathedral Street, Glasgow G4 0RE, U.K.
| |
Collapse
|
8
|
Napiórkowska E, Milcarz K, Szeleszczuk Ł. Review of Applications of Density Functional Theory (DFT) Quantum Mechanical Calculations to Study the High-Pressure Polymorphs of Organic Crystalline Materials. Int J Mol Sci 2023; 24:14155. [PMID: 37762459 PMCID: PMC10532210 DOI: 10.3390/ijms241814155] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/20/2023] [Revised: 09/05/2023] [Accepted: 09/11/2023] [Indexed: 09/29/2023] Open
Abstract
Since its inception, chemistry has been predominated by the use of temperature to generate or change materials, but applications of pressure of more than a few tens of atmospheres for such purposes have been rarely observed. However, pressure is a very effective thermodynamic variable that is increasingly used to generate new materials or alter the properties of existing ones. As computational approaches designed to simulate the solid state are normally tuned using structural data at ambient pressure, applying them to high-pressure issues is a highly challenging test of their validity from a computational standpoint. However, the use of quantum chemical calculations, typically at the level of density functional theory (DFT), has repeatedly been shown to be a great tool that can be used to both predict properties that can be later confirmed by experimenters and to explain, at the molecular level, the observations of high-pressure experiments. This article's main goal is to compile, analyze, and synthesize the findings of works addressing the use of DFT in the context of molecular crystals subjected to high-pressure conditions in order to give a general overview of the possibilities offered by these state-of-the-art calculations.
Collapse
Affiliation(s)
| | | | - Łukasz Szeleszczuk
- Department of Organic and Physical Chemistry, Faculty of Pharmacy, Medical University of Warsaw, Banacha 1, 02-093 Warsaw, Poland; (E.N.); (K.M.)
| |
Collapse
|
9
|
O’Connor D, Bier I, Tom R, Hiszpanski AM, Steele BA, Marom N. Ab Initio Crystal Structure Prediction of the Energetic Materials LLM-105, RDX, and HMX. CRYSTAL GROWTH & DESIGN 2023; 23:6275-6289. [PMID: 38173900 PMCID: PMC10763925 DOI: 10.1021/acs.cgd.3c00027] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Figures] [Subscribe] [Scholar Register] [Received: 01/07/2023] [Revised: 08/02/2023] [Indexed: 01/05/2024]
Abstract
Crystal structure prediction (CSP) is performed for the energetic materials (EMs) LLM-105 and α-RDX, as well as the α and β conformational polymorphs of 1,3,5,7-tetranitro-1,3,5,7-tetraazacyclooctane (HMX), using the genetic algorithm (GA) code, GAtor, and its associated random structure generator, Genarris. Genarris and GAtor successfully generate the experimental structures of all targets. GAtor's symmetric crossover scheme, where the space group symmetries of parent structures are treated as genes inherited by offspring, is found to be particularly effective. However, conducting several GA runs with different settings is still important for achieving diverse samplings of the potential energy surface. For LLM-105 and α-RDX, the experimental structure is ranked as the most stable, with all of the dispersion-inclusive density functional theory (DFT) methods used here. For HMX, the α form was persistently ranked as more stable than the β form, in contrast to experimental observations, even when correcting for vibrational contributions and thermal expansion. This may be attributed to insufficient accuracy of dispersion-inclusive DFT methods or to kinetic effects not considered here. In general, the ranking of some putative structures is found to be sensitive to the choice of the DFT functional and the dispersion method. For LLM-105, GAtor generates a putative structure with a layered packing motif, which is desirable thanks to its correlation with low sensitivity. Our results demonstrate that CSP is a useful tool for studying the ubiquitous polymorphism of EMs and shows promise of becoming an integral part of the EM development pipeline.
Collapse
Affiliation(s)
- Dana O’Connor
- Department
of Materials Science and Engineering, Carnegie
Mellon University, Pittsburgh, Pennsylvania 15213, United States
| | - Imanuel Bier
- Department
of Materials Science and Engineering, Carnegie
Mellon University, Pittsburgh, Pennsylvania 15213, United States
| | - Rithwik Tom
- Department
of Physics, Carnegie Mellon University, Pittsburgh, Pennsylvania 15213, United States
| | - Anna M. Hiszpanski
- Materials
Science Division, Lawrence Livermore National
Laboratory, Livermore, California 94550, United States
| | - Brad A. Steele
- Materials
Science Division, Lawrence Livermore National
Laboratory, Livermore, California 94550, United States
| | - Noa Marom
- Department
of Materials Science and Engineering, Carnegie
Mellon University, Pittsburgh, Pennsylvania 15213, United States
- Department
of Physics, Carnegie Mellon University, Pittsburgh, Pennsylvania 15213, United States
- Department
of Chemistry, Carnegie Mellon University, Pittsburgh, Pennsylvania 15213, United States
| |
Collapse
|
10
|
Butler PWV, Day GM. Reducing overprediction of molecular crystal structures via threshold clustering. Proc Natl Acad Sci U S A 2023; 120:e2300516120. [PMID: 37252993 PMCID: PMC10266058 DOI: 10.1073/pnas.2300516120] [Citation(s) in RCA: 3] [Impact Index Per Article: 3.0] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/11/2023] [Accepted: 05/01/2023] [Indexed: 06/01/2023] Open
Abstract
Crystal structure prediction is becoming an increasingly valuable tool for assessing polymorphism of crystalline molecular compounds, yet invariably, it overpredicts the number of polymorphs. One of the causes for this overprediction is in neglecting the coalescence of potential energy minima, separated by relatively small energy barriers, into a single basin at finite temperature. Considering this, we demonstrate a method underpinned by the threshold algorithm for clustering potential energy minima into basins, thereby identifying kinetically stable polymorphs and reducing overprediction.
Collapse
Affiliation(s)
- Patrick W. V. Butler
- School of Chemistry, University of Southampton, SouthamptonSO17 1BJ, United Kingdom
| | - Graeme M. Day
- School of Chemistry, University of Southampton, SouthamptonSO17 1BJ, United Kingdom
| |
Collapse
|
11
|
Rietveld IB, Akiba H, Yamamuro O, Barrio M, Céolin R, Tamarit JL. The Phase Diagram of the API Benzocaine and Its Highly Persistent, Metastable Crystalline Polymorphs. Pharmaceutics 2023; 15:pharmaceutics15051549. [PMID: 37242790 DOI: 10.3390/pharmaceutics15051549] [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: 04/06/2023] [Revised: 05/05/2023] [Accepted: 05/14/2023] [Indexed: 05/28/2023] Open
Abstract
The availability of sufficient amounts of form I of benzocaine has led to the investigation of its phase relationships with the other two existing forms, II and III, using adiabatic calorimetry, powder X-ray diffraction, and high-pressure differential thermal analysis. The latter two forms were known to have an enantiotropic phase relationship in which form III is stable at low-temperatures and high-pressures, while form II is stable at room temperature with respect to form III. Using adiabatic calorimetry data, it can be concluded, that form I is the stable low-temperature, high-pressure form, which also happens to be the most stable form at room temperature; however, due to its persistence at room temperature, form II is still the most convenient polymorph to use in formulations. Form III presents a case of overall monotropy and does not possess any stability domain in the pressure-temperature phase diagram. Heat capacity data for benzocaine have been obtained by adiabatic calorimetry from 11 K to 369 K above its melting point, which can be used to compare to results from in silico crystal structure prediction.
Collapse
Affiliation(s)
- Ivo B Rietveld
- Institute for Solid State Physics, University of Tokyo, 5-1-5 Kashiwanoha, Kashiwa 217-8581, Chiba, Japan
- Université Rouen Normandie, SMS, UR 3233, F-76000 Rouen, France
- Faculté de Pharmacie, Université Paris Cité, F-75006 Paris, France
| | - Hiroshi Akiba
- Institute for Solid State Physics, University of Tokyo, 5-1-5 Kashiwanoha, Kashiwa 217-8581, Chiba, Japan
| | - Osamu Yamamuro
- Institute for Solid State Physics, University of Tokyo, 5-1-5 Kashiwanoha, Kashiwa 217-8581, Chiba, Japan
| | - Maria Barrio
- Group de Caracterizació de Materials, Departament de Fisica, EEBE, Universitat Politècnica de Catalunya, Eduard Maristany, 10-14, 08019 Barcelona, Catalonia, Spain
- Barcelona Research Center in Multiscale Science and Engineering, Universitat Politècnica de Catalunya, Eduard Maristany, 10-14, 08019 Barcelona, Catalonia, Spain
| | - René Céolin
- Group de Caracterizació de Materials, Departament de Fisica, EEBE, Universitat Politècnica de Catalunya, Eduard Maristany, 10-14, 08019 Barcelona, Catalonia, Spain
| | - Josep-Lluís Tamarit
- Group de Caracterizació de Materials, Departament de Fisica, EEBE, Universitat Politècnica de Catalunya, Eduard Maristany, 10-14, 08019 Barcelona, Catalonia, Spain
- Barcelona Research Center in Multiscale Science and Engineering, Universitat Politècnica de Catalunya, Eduard Maristany, 10-14, 08019 Barcelona, Catalonia, Spain
| |
Collapse
|
12
|
Athiyarath V, Mathew LA, Zhao Y, Khazeber R, Ramamurty U, Sureshan KM. Rational design and topochemical synthesis of polymorphs of a polymer. Chem Sci 2023; 14:5132-5140. [PMID: 37206383 PMCID: PMC10189859 DOI: 10.1039/d3sc00053b] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/04/2023] [Accepted: 04/13/2023] [Indexed: 05/21/2023] Open
Abstract
Packing a polymer in different ways can give polymorphs of the polymer having different properties. β-Turn forming peptides such as 2-aminoisobutyric acid (Aib)-rich peptides adopt several conformations by varying the dihedral angles. Aiming at this, a β-turn-forming peptide monomer would give different polymorphs and these polymorphs upon topochemical polymerization would yield polymorphs of the polymer, we designed an Aib-rich monomer N3-(Aib)3-NHCH2-C[triple bond, length as m-dash]CH. This monomer crystallizes as two polymorphs and one hydrate. In all forms, the peptide adopts β-turn conformations and arranges in a head-to-tail manner with their azide and alkyne units proximally placed in a ready-to-react alignment. On heating, both the polymorphs undergo topochemical azide-alkyne cycloaddition polymerization. Polymorph I polymerized in a single-crystal-to-single-crystal (SCSC) fashion and the single-crystal X-ray diffraction analysis of the polymer revealed its screw-sense reversing helical structure. Polymorph II maintains its crystallinity during polymerization but gradually becomes amorphous upon storage. The hydrate III undergoes a dehydrative transition to polymorph II. Nanoindentation studies revealed that different polymorphs of the monomer and the corresponding polymers exhibited different mechanical properties, in accordance with their crystal packing. This work demonstrates the promising future of the marriage of polymorphism and topochemistry for obtaining polymorphs of polymers.
Collapse
Affiliation(s)
- Vignesh Athiyarath
- School of Chemistry, Indian Institute of Science Education and Research Thiruvananthapuram Kerala 695551 India
| | - Liby Ann Mathew
- School of Chemistry, Indian Institute of Science Education and Research Thiruvananthapuram Kerala 695551 India
| | - Yakai Zhao
- School of Mechanical and Aerospace Engineering, Nanyang Technological University 639798 Singapore
| | - Ravichandran Khazeber
- School of Chemistry, Indian Institute of Science Education and Research Thiruvananthapuram Kerala 695551 India
| | - Upadrasta Ramamurty
- School of Mechanical and Aerospace Engineering, Nanyang Technological University 639798 Singapore
| | - Kana M Sureshan
- School of Chemistry, Indian Institute of Science Education and Research Thiruvananthapuram Kerala 695551 India
| |
Collapse
|
13
|
Andrews J, Yufit DS, McCabe JF, Fox MA, Steed JW. Vapor Sorption and Halogen-Bond-Induced Solid-Form Rearrangement of a Porous Pharmaceutical. CRYSTAL GROWTH & DESIGN 2023; 23:2628-2633. [PMID: 37038401 PMCID: PMC10080649 DOI: 10.1021/acs.cgd.2c01464] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Figures] [Subscribe] [Scholar Register] [Received: 12/09/2022] [Revised: 02/08/2023] [Indexed: 06/19/2023]
Abstract
A porous, nonsolvated polymorph of the voltage-gated sodium channel blocker mexiletine hydrochloride absorbs iodine vapor to give a pharmaceutical cocrystal incorporating an I2Cl- anion that forms a halogen-π interaction with the mexiletine cations. The most thermodynamically stable form of the compound does not absorb iodine. This example shows that vapor sorption is a potentially useful and underused tool for bringing about changes in pharmaceutical solid form as part of a solid form screening protocol.
Collapse
Affiliation(s)
- Jessica
L. Andrews
- Department
of Chemistry, Durham University, South Road, Durham DH1 3LE, U.K.
| | - Dmitry S. Yufit
- Department
of Chemistry, Durham University, South Road, Durham DH1 3LE, U.K.
| | - James F. McCabe
- Pharmaceutical
Sciences, R&D, AstraZeneca, Charter Way, Silk Road Business
Park, Macclesfield SK10
2NA, U.K.
| | - Mark A. Fox
- Department
of Chemistry, Durham University, South Road, Durham DH1 3LE, U.K.
| | - Jonathan W. Steed
- Department
of Chemistry, Durham University, South Road, Durham DH1 3LE, U.K.
| |
Collapse
|
14
|
Tom R, Gao S, Yang Y, Zhao K, Bier I, Buchanan EA, Zaykov A, Havlas Z, Michl J, Marom N. Inverse Design of Tetracene Polymorphs with Enhanced Singlet Fission Performance by Property-Based Genetic Algorithm Optimization. CHEMISTRY OF MATERIALS : A PUBLICATION OF THE AMERICAN CHEMICAL SOCIETY 2023; 35:1373-1386. [PMID: 36999121 PMCID: PMC10042130 DOI: 10.1021/acs.chemmater.2c03444] [Citation(s) in RCA: 3] [Impact Index Per Article: 3.0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Figures] [Subscribe] [Scholar Register] [Received: 11/17/2022] [Revised: 01/06/2023] [Indexed: 06/19/2023]
Abstract
The efficiency of solar cells may be improved by using singlet fission (SF), in which one singlet exciton splits into two triplet excitons. SF occurs in molecular crystals. A molecule may crystallize in more than one form, a phenomenon known as polymorphism. Crystal structure may affect SF performance. In the common form of tetracene, SF is experimentally known to be slightly endoergic. A second, metastable polymorph of tetracene has been found to exhibit better SF performance. Here, we conduct inverse design of the crystal packing of tetracene using a genetic algorithm (GA) with a fitness function tailored to simultaneously optimize the SF rate and the lattice energy. The property-based GA successfully generates more structures predicted to have higher SF rates and provides insight into packing motifs associated with improved SF performance. We find a putative polymorph predicted to have superior SF performance to the two forms of tetracene, whose structures have been determined experimentally. The putative structure has a lattice energy within 1.5 kJ/mol of the most stable common form of tetracene.
Collapse
Affiliation(s)
- Rithwik Tom
- Department
of Physics, Carnegie Mellon University, Pittsburgh, Pennsylvania15213, United States
| | - Siyu Gao
- Department
of Materials Science and Engineering, Carnegie
Mellon University, Pittsburgh, Pennsylvania15213, United States
| | - Yi Yang
- Department
of Materials Science and Engineering, Carnegie
Mellon University, Pittsburgh, Pennsylvania15213, United States
| | - Kaiji Zhao
- Department
of Materials Science and Engineering, Carnegie
Mellon University, Pittsburgh, Pennsylvania15213, United States
| | - Imanuel Bier
- Department
of Materials Science and Engineering, Carnegie
Mellon University, Pittsburgh, Pennsylvania15213, United States
| | - Eric A. Buchanan
- Department
of Chemistry, University of Colorado, Boulder, Colorado80309, United States
| | - Alexandr Zaykov
- Institute
of Organic Chemistry and Biochemistry, Czech
Academy of Sciences, 16610Prague 6, Czech
Republic
- Department
of Physical Chemistry, University of Chemistry
and Technology, 166 28Prague 6, Czech Republic
| | - Zdeněk Havlas
- Institute
of Organic Chemistry and Biochemistry, Czech
Academy of Sciences, 16610Prague 6, Czech
Republic
| | - Josef Michl
- Department
of Chemistry, University of Colorado, Boulder, Colorado80309, United States
- Institute
of Organic Chemistry and Biochemistry, Czech
Academy of Sciences, 16610Prague 6, Czech
Republic
| | - Noa Marom
- Department
of Physics, Carnegie Mellon University, Pittsburgh, Pennsylvania15213, United States
- Department
of Materials Science and Engineering, Carnegie
Mellon University, Pittsburgh, Pennsylvania15213, United States
- Department
of Chemistry, Carnegie Mellon University, Pittsburgh, Pennsylvania15213, United States
| |
Collapse
|
15
|
Abe H, Kishimura H, Uruichi M. A phase variety of fluorinated ionic liquids: Molecular conformational and crystal polymorph. SPECTROCHIMICA ACTA. PART A, MOLECULAR AND BIOMOLECULAR SPECTROSCOPY 2023; 286:121948. [PMID: 36252301 DOI: 10.1016/j.saa.2022.121948] [Citation(s) in RCA: 1] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Received: 08/21/2022] [Revised: 09/26/2022] [Accepted: 10/03/2022] [Indexed: 06/16/2023]
Abstract
Crystal polymorphs of fluorinated ionic liquids (fILs) were examined at low-temperature (LT) by Raman spectroscopy. The fILs were 1-alkyl-3-methylimidazolium perfluorobutanesulfonate, [Cnmim][PFBS] (n = 4, 6, and 8). The cations and anion possess conformational degrees of freedom. Various LT phases were derived from the conformational polymorphs of the cations and the anion. Conformational flexibility depended on alkyl chain length. The crystal polymorphs in the fILs were sensitive to molecular conformations and flexibility.
Collapse
Affiliation(s)
- Hiroshi Abe
- Department of Materials Science and Engineering, National Defense Academy, Yokosuka 239-8686, Japan
| | - Hiroaki Kishimura
- Department of Materials Science and Engineering, National Defense Academy, Yokosuka 239-8686, Japan
| | - Mikio Uruichi
- Institute for Molecular Science, Myodaiji, Okazaki 444-8585, Japan
| |
Collapse
|
16
|
Cook CJ, Li W, Lui BF, Gately TJ, Al-Kaysi RO, Mueller LJ, Bardeen CJ, Beran GJO. A theoretical framework for the design of molecular crystal engines. Chem Sci 2023; 14:937-949. [PMID: 36755715 PMCID: PMC9890974 DOI: 10.1039/d2sc05549j] [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: 10/06/2022] [Accepted: 12/19/2022] [Indexed: 12/24/2022] Open
Abstract
Photomechanical molecular crystals have garnered attention for their ability to transform light into mechanical work, but difficulties in characterizing the structural changes and mechanical responses experimentally have hindered the development of practical organic crystal engines. This study proposes a new computational framework for predicting the solid-state crystal-to-crystal photochemical transformations entirely from first principles, and it establishes a photomechanical engine cycle that quantifies the anisotropic mechanical performance resulting from the transformation. The approach relies on crystal structure prediction, solid-state topochemical principles, and high-quality electronic structure methods. After validating the framework on the well-studied [4 + 4] cycloadditions in 9-methyl anthracene and 9-tert-butyl anthracene ester, the experimentally-unknown solid-state transformation of 9-carboxylic acid anthracene is predicted for the first time. The results illustrate how the mechanical work is done by relaxation of the crystal lattice to accommodate the photoproduct, rather than by the photochemistry itself. The large ∼107 J m-3 work densities computed for all three systems highlight the promise of photomechanical crystal engines. This study demonstrates the importance of crystal packing in determining molecular crystal engine performance and provides tools and insights to design improved materials in silico.
Collapse
Affiliation(s)
- Cameron J. Cook
- Department of Chemistry, University of California RiversideRiverside CA 92521USA
| | - Wangxiang Li
- Department of Chemistry, University of California Riverside Riverside CA 92521 USA
| | - Brandon F. Lui
- Department of Chemistry, University of California RiversideRiverside CA 92521USA
| | - Thomas J. Gately
- Department of Chemistry, University of California RiversideRiverside CA 92521USA
| | - Rabih O. Al-Kaysi
- College of Science and Health Professions-3124, King Saud Bin Abdulaziz University for Health Sciences, and King Abdullah International Medical Research Center, Ministry of National Guard Health AffairsRiyadh 11426Kingdom of Saudi Arabia
| | - Leonard J. Mueller
- Department of Chemistry, University of California RiversideRiverside CA 92521USA
| | | | - Gregory J. O. Beran
- Department of Chemistry, University of California RiversideRiverside CA 92521USA
| |
Collapse
|
17
|
Phase Transitions and Stabilities among Three Phases of Di-p-tolyl Disulfides. Molecules 2022; 27:molecules27238342. [PMID: 36500435 PMCID: PMC9739323 DOI: 10.3390/molecules27238342] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/25/2022] [Revised: 11/22/2022] [Accepted: 11/28/2022] [Indexed: 12/05/2022] Open
Abstract
Di-p-tolyl disulfides (p-Tol2S2) are employed as load-carrying additives because of their anti-wear and extreme load-bearing qualities. External pressure triggers conformational up-conversion (leads to phase transition) in the molecules of p-Tol2S2, by compensating for the stress and absorbing its energy. These features make p-Tol2S2 a potential candidate for next-generation energy storage devices. Upon lithiation, MoS2 expands up to 103% which causes stress and affects battery stability and performance. Therefore, it is essential to study these materials under different physical conditions. In this work, we used density functional theory (DFT) at ωB97XD/6-31G* functional level, to calculate lattice parameters, Gibbs free energies, and vibrational spectra of three phases (i.e., α, β, and γ) of p-Tol2S2 under different pressure and temperature conditions. The phase transition between phases α and β occurred at a pressure and temperature of 0.65 GPa and 463 K, respectively. Furthermore, phase transition between phases α and γ was found at a pressure and temperature of 0.35 GPa and 400 K, respectively. Moreover, no phase transition was observed between phases β and γ under the pressure range studied (0 GPa to 5.5 GPa). We also computed and compared the FT-IR spectra of the three phases. These results can guide scientists and chemists in designing more stable battery materials.
Collapse
|
18
|
Mattei A, Hong RS, Dietrich H, Firaha D, Helfferich J, Liu YM, Sasikumar K, Abraham NS, Miglani Bhardwaj R, Neumann MA, Sheikh AY. Efficient Crystal Structure Prediction for Structurally Related Molecules with Accurate and Transferable Tailor-Made Force Fields. J Chem Theory Comput 2022; 18:5725-5738. [PMID: 35930763 PMCID: PMC9476662 DOI: 10.1021/acs.jctc.2c00451] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/06/2023]
Abstract
Crystal structure prediction (CSP) his generally used to complement experimental solid form screening and applied to individual molecules in drug development. The fast development of algorithms and computing resources offers the opportunity to use CSP earlier and for a broader range of applications in the drug design cycle. This study presents a novel paradigm of CSP specifically designed for structurally related molecules, referred to as Quick-CSP. The approach prioritizes more accurate physics through robust and transferable tailor-made force fields (TMFFs), such that significant efficiency gains are achieved through the reduction of expensive ab initio calculations. The accuracy of the TMFF is increased by the introduction of electrostatic multipoles, and the fragment-based force field parameterization scheme is demonstrated to be transferable for a family of chemically related molecules. The protocol is benchmarked with structurally related compounds from the Bromodomain and Extraterminal (BET) domain inhibitors series. A new convergence criterion is introduced that aims at performing only as many ab initio optimizations of crystal structures as required to locate the bottom of the crystal energy landscape within a user-defined accuracy. The overall approach provides significant cost savings ranging from three- to eight-fold less than the full-CSP workflow. The reported advancements expand the scope and utility of the underlying CSP building blocks as well as their novel reassembly to other applications earlier in the drug design cycle to guide molecule design and selection.
Collapse
Affiliation(s)
- Alessandra Mattei
- Solid State Chemistry, Research & Development, AbbVie Inc., 1 N Waukegan Road, North Chicago, Illinois 60064, United States
| | - Richard S Hong
- Solid State Chemistry, Research & Development, AbbVie Inc., 1 N Waukegan Road, North Chicago, Illinois 60064, United States
| | - Hanno Dietrich
- Avant-garde Materials Simulation, GmbH, Alte Str. 2, 79249 Merzhausen, Germany
| | - Dzmitry Firaha
- Avant-garde Materials Simulation, GmbH, Alte Str. 2, 79249 Merzhausen, Germany
| | - Julian Helfferich
- Avant-garde Materials Simulation, GmbH, Alte Str. 2, 79249 Merzhausen, Germany
| | - Yifei Michelle Liu
- Avant-garde Materials Simulation, GmbH, Alte Str. 2, 79249 Merzhausen, Germany
| | - Kiran Sasikumar
- Avant-garde Materials Simulation, GmbH, Alte Str. 2, 79249 Merzhausen, Germany
| | - Nathan S Abraham
- Solid State Chemistry, Research & Development, AbbVie Inc., 1 N Waukegan Road, North Chicago, Illinois 60064, United States
| | - Rajni Miglani Bhardwaj
- Solid State Chemistry, Research & Development, AbbVie Inc., 1 N Waukegan Road, North Chicago, Illinois 60064, United States
| | - Marcus A Neumann
- Avant-garde Materials Simulation, GmbH, Alte Str. 2, 79249 Merzhausen, Germany
| | - Ahmad Y Sheikh
- Solid State Chemistry, Research & Development, AbbVie Inc., 1 N Waukegan Road, North Chicago, Illinois 60064, United States
| |
Collapse
|
19
|
Newman JA, Iuzzolino L, Tan M, Orth P, Bruhn J, Lee AY. From Powders to Single Crystals: A Crystallographer's Toolbox for Small-Molecule Structure Determination. Mol Pharm 2022; 19:2133-2141. [PMID: 35576503 PMCID: PMC10152450 DOI: 10.1021/acs.molpharmaceut.2c00020] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Abstract
Although the crystal structures of small-molecule compounds are often determined from single-crystal X-ray diffraction (scXRD), recent advances in three-dimensional electron diffraction (3DED) and crystal structure prediction (CSP) methods promise to expand the structure elucidation toolbox available to the crystallographer. Herein, a comparative assessment of scXRD, 3DED, and CSP in combination with powder X-ray diffraction is carried out on two former drug candidate compounds and a multicomponent crystal of a key building block in the synthesis of gefapixant citrate.
Collapse
Affiliation(s)
- Justin A. Newman
- Department
of Analytical Research and Development, Merck & Co., Inc., 126 East Lincoln Avenue, Rahway, New Jersey 07065, United States
| | - Luca Iuzzolino
- Department
of Computational and Structural Chemistry, Merck & Co., Inc., Rahway, 126 East Lincoln Avenue, Rahway, New Jersey 07065, United States
| | - Melissa Tan
- Department
of Analytical Research and Development, Merck & Co., Inc., 126 East Lincoln Avenue, Rahway, New Jersey 07065, United States
| | - Peter Orth
- Department
of Computational and Structural Chemistry, Merck & Co., Inc., 2000 Galloping Hill Road, Kenilworth, New Jersey 07033, United States
| | - Jessica Bruhn
- Nanoimaging
Services, San Diego, California 92121, United States
| | - Alfred Y. Lee
- Department
of Analytical Research and Development, Merck & Co., Inc., 126 East Lincoln Avenue, Rahway, New Jersey 07065, United States
| |
Collapse
|
20
|
Mathew R, Sergeyev IV, Aussenac F, Gkoura L, Rosay M, Baias M. Complete resonance assignment of a pharmaceutical drug at natural isotopic abundance from DNP-Enhanced solid-state NMR. SOLID STATE NUCLEAR MAGNETIC RESONANCE 2022; 119:101794. [PMID: 35462269 DOI: 10.1016/j.ssnmr.2022.101794] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Received: 10/20/2021] [Revised: 03/31/2022] [Accepted: 04/01/2022] [Indexed: 06/14/2023]
Abstract
Solid-state dynamic nuclear polarization enhanced magic angle spinning (DNP-MAS) NMR measurements coupled with density functional theory (DFT) calculations enable the full resonance assignment of a complex pharmaceutical drug molecule without the need for isotopic enrichment. DNP dramatically enhances the NMR signals, thereby making possible previously intractable two-dimensional correlation NMR spectra at natural abundance. Using inputs from DFT calculations, herein we describe a significant improvement to the structure elucidation process for complex organic molecules. Further, we demonstrate that a series of two-dimensional correlation experiments, including 15N-13C TEDOR, 13C-13C INADEQUATE/SARCOSY, 19F-13C HETCOR, and 1H-13C HETCOR, can be obtained at natural isotopic abundance within reasonable experiment times, thus enabling a complete resonance assignment of sitagliptin, a pharmaceutical used for the treatment of type 2 diabetes.
Collapse
Affiliation(s)
- Renny Mathew
- Division of Science, New York University Abu Dhabi, P.O. Box 129188, Abu Dhabi, United Arab Emirates
| | - Ivan V Sergeyev
- Bruker Biospin Corporation, 15 Fortune Drive, Billerica, MA, USA
| | - Fabien Aussenac
- Bruker France, 34 rue de l'industrie, 67166, Wissembourg, France.
| | - Lydia Gkoura
- Division of Science, New York University Abu Dhabi, P.O. Box 129188, Abu Dhabi, United Arab Emirates.
| | - Melanie Rosay
- Bruker Biospin Corporation, 15 Fortune Drive, Billerica, MA, USA
| | - Maria Baias
- Division of Science, New York University Abu Dhabi, P.O. Box 129188, Abu Dhabi, United Arab Emirates
| |
Collapse
|
21
|
Molecular Simulation Analyses of Polymorphism Control Factors by the example of Carbamazepine Forms I-IV: a Blueprint for Industrial Drug Formulation? J Pharm Sci 2022; 111:2898-2906. [DOI: 10.1016/j.xphs.2022.06.001] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/29/2022] [Revised: 06/02/2022] [Accepted: 06/02/2022] [Indexed: 11/18/2022]
|
22
|
High-Pressure Polymorphism in Hydrogen-Bonded Crystals: A Concise Review. CRYSTALS 2022. [DOI: 10.3390/cryst12050739] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 11/16/2022]
Abstract
High-pressure polymorphism is a developing interdisciplinary field. Pressure up to 20 GPa is a powerful thermodynamic parameter for the study and fabrication of hydrogen-bonded polymorphic systems. This review describes how pressure can be used to explore polymorphism and surveys the reports on examples of compounds that our group has studied at high pressures. Such studies have provided insight into the nature of structure–property relationships, which will enable crystal engineering to design crystals with desired architectures through hydrogen-bonded networks. Experimental methods are also briefly surveyed, along with two methods that have proven to be very helpful in the analysis of high-pressure polymorphs, namely, the ab initio pseudopotential plane–wave density functional method and using Hirshfeld surfaces to construct a graphical overview of intermolecular interactions.
Collapse
|
23
|
Abe H, Koyama Y, Shimono S, Kishimura H, Matsuishi K. High-pressure crystal polymorphs and multiple pathways in 1-hexyl-3-methylimidazolium perfluorobutanesulfonate ionic liquid. Chem Phys 2022. [DOI: 10.1016/j.chemphys.2022.111479] [Citation(s) in RCA: 2] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/27/2022]
|
24
|
Beran GJO, Wright SE, Greenwell C, Cruz-Cabeza AJ. The interplay of intra- and intermolecular errors in modeling conformational polymorphs. J Chem Phys 2022; 156:104112. [DOI: 10.1063/5.0088027] [Citation(s) in RCA: 4] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/26/2023] Open
Abstract
Conformational polymorphs of organic molecular crystals represent a challenging test for quantum chemistry because they require careful balancing of the intra- and intermolecular interactions. This study examines 54 molecular conformations from 20 sets of conformational polymorphs, along with the relative lattice energies and 173 dimer interactions taken from six of the polymorph sets. These systems are studied with a variety of van der Waals-inclusive density functionals theory models; dispersion-corrected spin-component-scaled second-order Møller–Plesset perturbation theory (SCS-MP2D); and domain local pair natural orbital coupled cluster singles, doubles, and perturbative triples [DLPNO-CCSD(T)]. We investigate how delocalization error in conventional density functionals impacts monomer conformational energies, systematic errors in the intermolecular interactions, and the nature of error cancellation that occurs in the overall crystal. The density functionals B86bPBE-XDM, PBE-D4, PBE-MBD, PBE0-D4, and PBE0-MBD are found to exhibit sizable one-body and two-body errors vs DLPNO-CCSD(T) benchmarks, and the level of success in predicting the relative polymorph energies relies heavily on error cancellation between different types of intermolecular interactions or between intra- and intermolecular interactions. The SCS-MP2D and, to a lesser extent, ωB97M-V models exhibit smaller errors and rely less on error cancellation. Implications for crystal structure prediction of flexible compounds are discussed. Finally, the one-body and two-body DLPNO-CCSD(T) energies taken from these conformational polymorphs establish the CP1b and CP2b benchmark datasets that could be useful for testing quantum chemistry models in challenging real-world systems with complex interplay between intra- and intermolecular interactions, a number of which are significantly impacted by delocalization error.
Collapse
Affiliation(s)
- Gregory J. O. Beran
- Department of Chemistry, University of California, Riverside, California 92521, USA
| | - Sarah E. Wright
- Department of Chemical Engineering and Analytical Science, University of Manchester, Manchester, United Kingdom
| | - Chandler Greenwell
- Department of Chemistry, University of California, Riverside, California 92521, USA
| | - Aurora J. Cruz-Cabeza
- Department of Chemical Engineering and Analytical Science, University of Manchester, Manchester, United Kingdom
| |
Collapse
|
25
|
Kargbo RB, Sherwood AM, Meisenheimer P, Lenoch K, Abebe S. Psilocybin: Characterization of the Metastable Zone Width (MSZW), Control of Anhydrous Polymorphs, and Particle Size Distribution (PSD). ACS OMEGA 2022; 7:5429-5436. [PMID: 35187358 PMCID: PMC8851641 DOI: 10.1021/acsomega.1c06708] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Subscribe] [Scholar Register] [Received: 11/29/2021] [Accepted: 01/17/2022] [Indexed: 06/14/2023]
Abstract
Psilocybin, a serotonergic agonist, was granted a "breakthrough therapy" status by the Food and Drug Administration for clinical trials involving major depressive disorder and treatment-resistant depression. The direct phosphorylation of psilocin to psilocybin that uses a fast crystallization associated with a kinetically controlled process resulted in a smaller particle size distribution. Herein, the measurement of the metastable zone width (MSZW) and nucleation induction enabled a thermodynamically controlled crystallization process, which leads to the formation of a crystal structure with stronger interactions, controlled particle size distribution (PSD), and improved impurity profile. Employing a high-resolution inline microscopy viewer allowed the real-time monitoring of the crystallization process and the measurement of the particle size. We also present a comprehensive study of the formation of polymorph B (trihydrate), polymorph A (anhydrate), and polymorph H (anhydrate) using water recrystallization, which indicates that the formation of polymorph B (trihydrate) is independent of the crystallization method. However, polymorphs A and H are dependent on the mode of drying: drying at room temperature under vacuum gives rise to mainly polymorph A, and when heated even at relatively low temperatures, a mixture of polymorphs A and H beings to form.
Collapse
Affiliation(s)
- Robert B. Kargbo
- Usona
Institute, 2780 Woods
Hollow Road, Madison, Wisconsin 53711, United States
| | | | - Poncho Meisenheimer
- Usona
Institute, 2780 Woods
Hollow Road, Madison, Wisconsin 53711, United States
| | - Kelsey Lenoch
- Usona
Institute, 2780 Woods
Hollow Road, Madison, Wisconsin 53711, United States
| | - Solomon Abebe
- Almac
Sciences, 20 Seagoe Industrial
Estate, Craigavon BT63
5QD, United Kingdom
| |
Collapse
|
26
|
Momenzadeh Abardeh Z, Salimi A, Oganov AR. Crystal structure prediction of N-halide phthalimide compounds: Halogen bond synthon as a touchstone. CrystEngComm 2022. [DOI: 10.1039/d2ce00476c] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022]
Abstract
We address crystal structure prediction problem by combining evolutionary algorithm USPEX (used to predict sets of low-energy crystal structures) and synthon approach (extracting preferable supramolecular synthons from Cambridge Structural Database,...
Collapse
|
27
|
Jia S, Zhou Q, Li F, Hu Y, Wang C, Wang X, He S, Li X, Li L, Cui T. High-pressure Bandgap Engineering and Amorphization in TiNb2O7 Single Crystals. CrystEngComm 2022. [DOI: 10.1039/d2ce00168c] [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
Titanium niobate (TiNb2O7) possesses excellent photocatalytic, dielectric properties, and lithium-insertion capacity. And high-pressure (HP) is a powerful tool for bandgap engineering aiming at widening their applications. Herein, we report the...
Collapse
|
28
|
Dudek MK, Druzbicki K. Along the road to Crystal Structure Prediction (CSP) of pharmaceutical-like molecules. CrystEngComm 2022. [DOI: 10.1039/d1ce01564h] [Citation(s) in RCA: 2] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022]
Abstract
Computational methods used for predicting crystal structures of organic compounds are mature enough to be routinely used with many rigid and semi-rigid organic molecules. The usefulness of Crystal Structure Prediction...
Collapse
|
29
|
Botes DS, Khorasani S, Levendis DC, Fernandes MA. Accessing a regiospecific isomer and a metastable polymorph through crystal engineering and solid-state reaction. CrystEngComm 2022. [DOI: 10.1039/d2ce01094a] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022]
Abstract
We describe a solid-state Diels–Alder reaction where crystal engineering was used to design a reaction site yielding one regioisomer. Reaction was followed with SCXRD, compared to solution synthesis and rationalised using computational modelling.
Collapse
Affiliation(s)
- Delbert S. Botes
- Molecular Sciences Institute, School of Chemistry, University of the Witwatersrand, PO Wits 2050, Johannesburg, South Africa
| | - Sanaz Khorasani
- Molecular Sciences Institute, School of Chemistry, University of the Witwatersrand, PO Wits 2050, Johannesburg, South Africa
| | - Demetrius C. Levendis
- Molecular Sciences Institute, School of Chemistry, University of the Witwatersrand, PO Wits 2050, Johannesburg, South Africa
| | - Manuel A. Fernandes
- Molecular Sciences Institute, School of Chemistry, University of the Witwatersrand, PO Wits 2050, Johannesburg, South Africa
| |
Collapse
|
30
|
Beran GJO, Sugden IJ, Greenwell C, Bowskill DH, Pantelides CC, Adjiman CS. How many more polymorphs of ROY remain undiscovered. Chem Sci 2022; 13:1288-1297. [PMID: 35222912 PMCID: PMC8809489 DOI: 10.1039/d1sc06074k] [Citation(s) in RCA: 13] [Impact Index Per Article: 6.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/02/2021] [Accepted: 12/10/2021] [Indexed: 12/15/2022] Open
Abstract
With 12 crystal forms, 5-methyl-2-[(2-nitrophenyl)amino]-3-thiophenecabonitrile (a.k.a. ROY) holds the current record for the largest number of fully characterized organic crystal polymorphs. Four of these polymorph structures have been reported since 2019, raising the question of how many more ROY polymorphs await future discovery. Employing crystal structure prediction and accurate energy rankings derived from conformational energy-corrected density functional theory, this study presents the first crystal energy landscape for ROY that agrees well with experiment. The lattice energies suggest that the seven most stable ROY polymorphs (and nine of the twelve lowest-energy forms) on the Z′ = 1 landscape have already been discovered experimentally. Discovering any new polymorphs at ambient pressure will likely require specialized crystallization techniques capable of trapping metastable forms. At pressures above 10 GPa, however, a new crystal form is predicted to become enthalpically more stable than all known polymorphs, suggesting that further high-pressure experiments on ROY may be warranted. This work highlights the value of high-accuracy crystal structure prediction for solid-form screening and demonstrates how pragmatic conformational energy corrections can overcome the limitations of conventional density functionals for conformational polymorphs. Crystal structure prediction suggests that the low-energy polymorphs of ROY have already been found, but a new high-pressure form is predicted.![]()
Collapse
Affiliation(s)
- Gregory J. O. Beran
- Department of Chemistry, University of California Riverside, Riverside, CA 92521, USA
| | - Isaac J. Sugden
- Department of Chemical Engineering, Sargent Centre for Process Systems Engineering, Imperial College London, London, SW7 2AZ, UK
| | - Chandler Greenwell
- Department of Chemistry, University of California Riverside, Riverside, CA 92521, USA
| | - David H. Bowskill
- Department of Chemical Engineering, Sargent Centre for Process Systems Engineering, Imperial College London, London, SW7 2AZ, UK
| | - Constantinos C. Pantelides
- Department of Chemical Engineering, Sargent Centre for Process Systems Engineering, Imperial College London, London, SW7 2AZ, UK
| | - Claire S. Adjiman
- Department of Chemical Engineering, Sargent Centre for Process Systems Engineering, Imperial College London, London, SW7 2AZ, UK
| |
Collapse
|
31
|
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
|
32
|
|
33
|
Schlesinger C, Habermehl S, Prill D. Structure determination of organic compounds by a fit to the pair distribution function from scratch without prior indexing. J Appl Crystallogr 2021; 54:776-786. [PMID: 34188612 PMCID: PMC8202035 DOI: 10.1107/s1600576721002569] [Citation(s) in RCA: 6] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/04/2020] [Accepted: 03/08/2021] [Indexed: 11/20/2022] Open
Abstract
A method for the ab initio crystal structure determination of organic compounds by a fit to the pair distribution function (PDF), without prior knowledge of lattice parameters and space group, has been developed. The method is called 'PDF-Global-Fit' and is implemented by extension of the program FIDEL (fit with deviating lattice parameters). The structure solution is based on a global optimization approach starting from random structural models in selected space groups. No prior indexing of the powder data is needed. The new method requires only the molecular geometry and a carefully determined PDF. The generated random structures are compared with the experimental PDF and ranked by a similarity measure based on cross-correlation functions. The most promising structure candidates are fitted to the experimental PDF data using a restricted simulated annealing structure solution approach within the program TOPAS, followed by a structure refinement against the PDF to identify the correct crystal structure. With the PDF-Global-Fit it is possible to determine the local structure of crystalline and disordered organic materials, as well as to determine the local structure of unindexable powder patterns, such as nanocrystalline samples, by a fit to the PDF. The success of the method is demonstrated using barbituric acid as an example. The crystal structure of barbituric acid form IV solved and refined by the PDF-Global-Fit is in excellent agreement with the published crystal structure data.
Collapse
Affiliation(s)
- Carina Schlesinger
- Institut für Anorganische und Analytische Chemie, Goethe Universität, Max-von-Laue-Strasse 7, Frankfurt am Main, 60437, Germany
| | - Stefan Habermehl
- Institut für Anorganische und Analytische Chemie, Goethe Universität, Max-von-Laue-Strasse 7, Frankfurt am Main, 60437, Germany
| | - Dragica Prill
- Institut für Anorganische und Analytische Chemie, Goethe Universität, Max-von-Laue-Strasse 7, Frankfurt am Main, 60437, Germany
| |
Collapse
|
34
|
|
35
|
Cook C, McKinley JL, Beran GJO. Modeling the α- and β-resorcinol phase boundary via combination of density functional theory and density functional tight-binding. J Chem Phys 2021; 154:134109. [PMID: 33832233 PMCID: PMC8019358 DOI: 10.1063/5.0044385] [Citation(s) in RCA: 3] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/15/2021] [Accepted: 03/10/2021] [Indexed: 02/06/2023] Open
Abstract
The ability to predict not only what organic crystal structures might occur but also the thermodynamic conditions under which they are the most stable would be extremely useful for discovering and designing new organic materials. The present study takes a step in that direction by predicting the temperature- and pressure-dependent phase boundary between the α and β polymorphs of resorcinol using density functional theory (DFT) and the quasi-harmonic approximation. To circumvent the major computational bottleneck associated with computing a well-converged phonon density of states via the supercell approach, a recently developed approximation is employed, which combines a supercell phonon density of states from dispersion-corrected third-order density functional tight binding [DFTB3-D3(BJ)] with frequency corrections derived from a smaller B86bPBE-XDM functional DFT phonon calculation on the crystallographic unit cell. This mixed DFT/DFTB quasi-harmonic approach predicts the lattice constants and unit cell volumes to within 1%-2% at lower pressures. It predicts the thermodynamic phase boundary in almost perfect agreement with the experiment, although this excellent agreement does reflect fortuitous cancellation of errors between the enthalpy and entropy of transition.
Collapse
Affiliation(s)
- Cameron Cook
- Department of Chemistry, University of California, Riverside, California 92521, USA
| | - Jessica L. McKinley
- 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
|
36
|
Hao X, Liu J, Ali I, Luo H, Han Y, Hu W, Liu J, He X, Li J. Ab initio determination of crystal stability of di-p-tolyl disulfide. Sci Rep 2021; 11:7076. [PMID: 33782489 PMCID: PMC8007795 DOI: 10.1038/s41598-021-86519-1] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/01/2020] [Accepted: 03/16/2021] [Indexed: 11/09/2022] Open
Abstract
With the rapid growth of energy demand and the depletion of existing energy resources, the new materials with superior performances, low costs and environmental friendliness for energy production and storage are explored. Di-p-tolyl disulfide (p-Tol2S2) is a typical lubricating material, which has been applied in the field of energy storage. The conformational properties and phase transformations of p-Tol2S2 have been studied by pioneers, but their polymorphs and the polymorphism induced crystal structure changes require further analysis. In this study, we perform the crystal structural screening, prediction and optimization of p-Tol2S2 crystal with quantum mechanical calculations, i.e., density functional theory (DFT) and second-order Møller–Plesset perturbation (MP2) methods. A series of crystal structures with different molecular arrangements are generated based on the crystal structure screening. As compared to long-established lattice energy calculation, we take an advantage of using more accurate technique, which is Gibbs free energy calculation. It considers the effects of entropy and temperature to predict the crystal structures and energy landscape. By comparing the Gibbs free energies between predicted and experimental structures, we found that phase α is the most stable structure for p-Tol2S2 crystal at ambient temperature and standard atmospheric pressure. Furthermore, we provide an efficient method to discriminate different polymorphs that are otherwise difficult to be identified based on the Raman/IR spectra. The proposed work enable us to evaluate the quality of various crystal polymorphs rapidly.
Collapse
Affiliation(s)
- Xuan Hao
- School of Chemistry and Molecular Engineering, Shanghai Engineering Research Center of Molecular Therapeutics and New Drug Development, East China Normal University, Shanghai, 200062, China.,Key Laboratory for Thin Film and Microfabrication of Ministry of Education, Department of Micro/Nano-Electronics, Shanghai Jiao Tong University, Shanghai, 200240, China
| | - Jinfeng Liu
- State Key Laboratory of Natural Medicines, Department of Basic Medicine and Clinical Pharmacy, China Pharmaceutical University, Nanjing, 210009, China
| | - Imran Ali
- Key Laboratory for Thin Film and Microfabrication of Ministry of Education, Department of Micro/Nano-Electronics, Shanghai Jiao Tong University, Shanghai, 200240, China
| | - 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
| | - Yanqiang Han
- Key Laboratory for Thin Film and Microfabrication of Ministry of Education, Department of Micro/Nano-Electronics, Shanghai Jiao Tong University, Shanghai, 200240, China
| | - Wenxin Hu
- School of Computer Science and Software Engineering, The Computer Center, East China Normal University, Shanghai, 200062, China
| | - Jinyun Liu
- Anhui Provincial Engineering Laboratory for New-Energy Vehicle Battery Energy-Storage Materials, Anhui Laboratory of Molecule-Based Materials, College of Chemistry and Materials Science, Anhui Normal University, Wuhu, 241002, Anhui, China.
| | - Xiao He
- School of Chemistry and Molecular Engineering, Shanghai Engineering Research Center of Molecular Therapeutics and New Drug Development, East China Normal University, Shanghai, 200062, China. .,NYU-ECNU Center for Computational Chemistry at NYU Shanghai, Shanghai, 200062, China.
| | - Jinjin Li
- Key Laboratory for Thin Film and Microfabrication of Ministry of Education, Department of Micro/Nano-Electronics, Shanghai Jiao Tong University, Shanghai, 200240, China.
| |
Collapse
|
37
|
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
|
38
|
Tang J, Han Y, Ali I, Luo H, Nowak A, Li J. Stability and phase transition investigation of olanzapine polymorphs. Chem Phys Lett 2021. [DOI: 10.1016/j.cplett.2021.138384] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/22/2022]
|
39
|
Dittrich B. On modelling disordered crystal structures through restraints from molecule-in-cluster computations, and distinguishing static and dynamic disorder. IUCRJ 2021; 8:305-318. [PMID: 33708406 PMCID: PMC7924241 DOI: 10.1107/s2052252521000531] [Citation(s) in RCA: 4] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Figures] [Subscribe] [Scholar Register] [Received: 10/02/2020] [Accepted: 01/14/2021] [Indexed: 06/12/2023]
Abstract
Distinguishing disorder into static and dynamic based on multi-temperature X-ray or neutron diffraction experiments is the current state of the art, but is only descriptive, not predictive. Here, several disordered structures are revisited from the Cambridge Crystallographic Data Center 'drug subset', the Cambridge Structural Database and own earlier work, where experimental intensities of Bragg diffraction data were available. Using the molecule-in-cluster approach, structures with distinguishable conformations were optimized separately, as extracted from available or generated disorder models of the respective disordered crystal structures. Re-combining these 'archetype structures' by restraining positional and constraining displacement parameters for conventional least-squares refinement, based on the optimized geometries, then often achieves a superior fit to the experimental diffraction data compared with relying on experimental information alone. It also simplifies and standardizes disorder refinement. Ten example structures were analysed. It is observed that energy differences between separate disorder conformations are usually within a small energy window of RT (T = crystallization temperature). Further computations classify disorder into static or dynamic, using single experiments performed at one single temperature, and this was achieved for propionamide.
Collapse
Affiliation(s)
- Birger Dittrich
- Novartis Campus, Novartis Pharma AG, Postfach, Basel, CH-4002, Switzerland
- Mathematisch-Naturwissenschaftliche Fakultät, Universität Zürich, Winterthurerstrasse 190, Zürich, CH-8057, Switzerland
| |
Collapse
|
40
|
Mayo RA, Johnson ER. Improved quantitative crystal-structure comparison using powder diffractograms via anisotropic volume correction. CrystEngComm 2021. [DOI: 10.1039/d1ce01058a] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/18/2022]
Abstract
A new anisotropic volume correction improves quantitative crystal structure comparison. Benchmarking against the 6th crystal structure prediction blind test data results in identification of two previously uncredited matching structures.
Collapse
Affiliation(s)
- R. Alex Mayo
- Department of Chemistry, Dalhousie University, 6274 Coburg Road, PO Box 15000, Halifax, Nova Scotia, B3H 4R2, Canada
| | - Erin R. Johnson
- Department of Chemistry, Dalhousie University, 6274 Coburg Road, PO Box 15000, Halifax, Nova Scotia, B3H 4R2, Canada
| |
Collapse
|
41
|
Woodley SM, Day GM, Catlow R. Structure prediction of crystals, surfaces and nanoparticles. PHILOSOPHICAL TRANSACTIONS. SERIES A, MATHEMATICAL, PHYSICAL, AND ENGINEERING SCIENCES 2020; 378:20190600. [PMID: 33100162 DOI: 10.1098/rsta.2019.0600] [Citation(s) in RCA: 15] [Impact Index Per Article: 3.8] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 05/12/2023]
Abstract
We review the current techniques used in the prediction of crystal structures and their surfaces and of the structures of nanoparticles. The main classes of search algorithm and energy function are summarized, and we discuss the growing role of methods based on machine learning. We illustrate the current status of the field with examples taken from metallic, inorganic and organic systems. This article is part of a discussion meeting issue 'Dynamic in situ microscopy relating structure and function'.
Collapse
Affiliation(s)
- Scott M Woodley
- Department of Chemistry, University College London, 20 Gordon Street, London WC1H 0AJ, UK
| | - Graeme M Day
- Computational Systems Chemistry, School of Chemistry, University of Southampton, Southampton SO17 1BJ, UK
| | - R Catlow
- Department of Chemistry, University College London, 20 Gordon Street, London WC1H 0AJ, UK
- School of Chemistry, Cardiff University, Park Place, Cardiff CF10 3AT, UK
| |
Collapse
|
42
|
Song H, Vogt-Maranto L, Wiscons R, Matzger AJ, Tuckerman ME. Generating Cocrystal Polymorphs with Information Entropy Driven by Molecular Dynamics-Based Enhanced Sampling. J Phys Chem Lett 2020; 11:9751-9758. [PMID: 33141590 DOI: 10.1021/acs.jpclett.0c02647] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 06/11/2023]
Abstract
Predicting structures of organic molecular cocrystals is a challenging task when considering the immense number of possible intermolecular orientations. Use of the Shannon information entropy, constructed from an intermolecular orientational spatial distribution function, to drive a search for crystal structures via enhanced molecular dynamics can be an efficient way to map out a landscape of putative polymorphs. Here, the Shannon entropy is used to generate a set of collective variables for differentiating polymorphs of a 1:1 cocrystal of resorcinol and urea. We show that driven adiabatic free energy dynamics, a particular enhanced-sampling approach, combined with these entropy variables, can transform the stable phase into alternate polymorphs. Density functional theory calculations confirm that a structure obtained from the enhanced molecular dynamics is stable at pressures above 1 GPa. We thus show that enhanced sampling should be considered an integral component of crystal structure searching protocols for systems with multiple independent molecules.
Collapse
Affiliation(s)
- Hongxing Song
- Department of Chemistry, New York University, New York, New York 10003, United States
| | - Leslie Vogt-Maranto
- Department of Chemistry, New York University, New York, New York 10003, United States
| | - Ren Wiscons
- Department of Chemistry and Macromolecular Science and Engineering Program, University of Michigan, 930 North University Avenue, Ann Arbor, Michigan 48109-1055, United States
| | - Adam J Matzger
- Department of Chemistry and Macromolecular Science and Engineering Program, University of Michigan, 930 North University Avenue, Ann Arbor, Michigan 48109-1055, United States
| | - Mark E Tuckerman
- Department of Chemistry, New York University, New York, New York 10003, United States
- Courant Institute of Mathematical Sciences, New York University, New York, New York 10012, United States
- NYU-ECNU Center for Computational Chemistry at NYU Shanghai, 3663 Zhongshan Road North, Shanghai 200062, China
| |
Collapse
|
43
|
Montis R, Davey RJ, Wright SE, Woollam GR, Cruz‐Cabeza AJ. Transforming Computed Energy Landscapes into Experimental Realities: The Role of Structural Rugosity. Angew Chem Int Ed Engl 2020; 59:20357-20360. [PMID: 32730676 PMCID: PMC7693212 DOI: 10.1002/anie.202006939] [Citation(s) in RCA: 13] [Impact Index Per Article: 3.3] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/13/2020] [Indexed: 11/14/2022]
Abstract
We exploit the possible link between structural surface roughness and difficulty of crystallisation. Polymorphs with smooth surfaces may nucleate and crystallise more readily than polymorphs with rough surfaces. The concept is applied to crystal structure prediction landscapes and reveals a promising complementary way of ranking putative crystal structures.
Collapse
Affiliation(s)
- Riccardo Montis
- Department of Chemical Engineering and Analytical ScienceThe University of ManchesterOxford RoadManchesterM13 9PLUK
| | - Roger J. Davey
- Department of Chemical Engineering and Analytical ScienceThe University of ManchesterOxford RoadManchesterM13 9PLUK
| | - Sarah E. Wright
- Department of Chemical Engineering and Analytical ScienceThe University of ManchesterOxford RoadManchesterM13 9PLUK
| | | | - Aurora J. Cruz‐Cabeza
- Department of Chemical Engineering and Analytical ScienceThe University of ManchesterOxford RoadManchesterM13 9PLUK
| |
Collapse
|
44
|
Liu C, Brandenburg JG, Valsson O, Kremer K, Bereau T. Free-energy landscape of polymer-crystal polymorphism. SOFT MATTER 2020; 16:9683-9692. [PMID: 33000842 DOI: 10.1039/d0sm01342k] [Citation(s) in RCA: 6] [Impact Index Per Article: 1.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 06/11/2023]
Abstract
Polymorphism rationalizes how processing can control the final structure of a material. The rugged free-energy landscape and exceedingly slow kinetics in the solid state have so far hampered computational investigations. We report for the first time the free-energy landscape of a polymorphic crystalline polymer, syndiotactic polystyrene. Coarse-grained metadynamics simulations allow us to efficiently sample the landscape at large. The free-energy difference between the two main polymorphs, α and β, is further investigated by quantum-chemical calculations. The results of the two methods are in line with experimental observations: they predict β as the more stable polymorph under standard conditions. Critically, the free-energy landscape suggests how the α polymorph may lead to experimentally observed kinetic traps. The combination of multiscale modeling, enhanced sampling, and quantum-chemical calculations offers an appealing strategy to uncover complex free-energy landscapes with polymorphic behavior.
Collapse
Affiliation(s)
- Chan Liu
- Max Planck Institute for Polymer Research, 55128 Mainz, Germany
| | | | | | | | | |
Collapse
|
45
|
Li X, Ou X, Wang B, Rong H, Wang B, Chang C, Shi B, Yu L, Lu M. Rich polymorphism in nicotinamide revealed by melt crystallization and crystal structure prediction. Commun Chem 2020; 3:152. [PMID: 36703331 PMCID: PMC9814109 DOI: 10.1038/s42004-020-00401-1] [Citation(s) in RCA: 19] [Impact Index Per Article: 4.8] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/15/2020] [Accepted: 10/09/2020] [Indexed: 01/29/2023] Open
Abstract
Overprediction is a major limitation of current crystal structure prediction (CSP) methods. It is difficult to determine whether computer-predicted polymorphic structures are artefacts of the calculation model or are polymorphs that have not yet been found. Here, we reported the well-known vitamin nicotinamide (NIC) to be a highly polymorphic compound with nine solved single-crystal structures determined by performing melt crystallization. A CSP calculation successfully identifies all six Z' = 1 and 2 experimental structures, five of which defy 66 years of attempts at being explored using solution crystallization. Our study demonstrates that when combined with our strategy for cultivating single crystals from melt microdroplets, melt crystallization has turned out to be an efficient tool for exploring polymorphic landscapes to better understand polymorphic crystallization and to more effectively test the accuracy of theoretical predictions, especially in regions inaccessible by solution crystallization.
Collapse
Affiliation(s)
- Xizhen Li
- grid.12981.330000 0001 2360 039XSchool of Pharmaceutical Sciences, Sun Yat-sen University, Guangzhou, China
| | - Xiao Ou
- grid.12981.330000 0001 2360 039XSchool of Pharmaceutical Sciences, Sun Yat-sen University, Guangzhou, China
| | - Bingquan Wang
- grid.12981.330000 0001 2360 039XSchool of Pharmaceutical Sciences, Sun Yat-sen University, Guangzhou, China
| | - Haowei Rong
- grid.12981.330000 0001 2360 039XSchool of Pharmaceutical Sciences, Sun Yat-sen University, Guangzhou, China
| | - Bing Wang
- Shenzhen Jingtai Technology Co., Ltd. (XtalPi Inc.), Shenzhen, China
| | - Chao Chang
- Shenzhen Jingtai Technology Co., Ltd. (XtalPi Inc.), Shenzhen, China
| | - Baimei Shi
- Shenzhen Jingtai Technology Co., Ltd. (XtalPi Inc.), Shenzhen, China
| | - Lian Yu
- grid.14003.360000 0001 2167 3675School of Pharmacy, University of Wisconsin – Madison, Madison, WI USA
| | - Ming Lu
- grid.12981.330000 0001 2360 039XSchool of Pharmaceutical Sciences, Sun Yat-sen University, Guangzhou, China ,grid.12981.330000 0001 2360 039XGuangdong Provincial Key Laboratory of New Drug Design and Evaluation, Sun Yat-sen University, Guangzhou, China
| |
Collapse
|
46
|
Abe H, Koyama Y, Kishimura H, Matsuishi K. High-pressure crystal polymorph of the protic ionic liquid: Ethylammonium nitrate. J Mol Liq 2020. [DOI: 10.1016/j.molliq.2020.113959] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/23/2022]
|
47
|
Reutzel-Edens SM, Bhardwaj RM. Crystal forms in pharmaceutical applications: olanzapine, a gift to crystal chemistry that keeps on giving. IUCRJ 2020; 7:955-964. [PMID: 33209310 PMCID: PMC7642794 DOI: 10.1107/s2052252520012683] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Figures] [Subscribe] [Scholar Register] [Received: 08/26/2020] [Accepted: 09/17/2020] [Indexed: 06/11/2023]
Abstract
This contribution reviews the efforts of many scientists around the world to discover and structurally characterize olanzapine crystal forms, clearing up inconsistencies in the scientific and patent literature and highlighting the challenges in identifying new forms amidst 60+ known polymorphs and solvates. Owing to its remarkable solid-state chemistry, olanzapine has emerged over the last three decades as a popular tool compound for developing new experimental and computational methods for enhanced molecular level understanding of solid-state structure, form diversity and crystallization outcomes. This article highlights the role of olanzapine in advancing the fundamental understanding of crystal forms, interactions within crystal structures, and growth units in molecular crystallization, as well as influencing the way in which drugs are developed today.
Collapse
Affiliation(s)
- Susan M. Reutzel-Edens
- Synthetic Molecule Design and Development, Eli Lilly and Company, Indianapolis, IN 46285, USA
| | - Rajni M. Bhardwaj
- Synthetic Molecule Design and Development, Eli Lilly and Company, Indianapolis, IN 46285, USA
| |
Collapse
|
48
|
Blade H, Blundell CD, Brown SP, Carson J, Dannatt HRW, Hughes LP, Menakath AK. Conformations in Solution and in Solid-State Polymorphs: Correlating Experimental and Calculated Nuclear Magnetic Resonance Chemical Shifts for Tolfenamic Acid. J Phys Chem A 2020; 124:8959-8977. [DOI: 10.1021/acs.jpca.0c07000] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/24/2022]
Affiliation(s)
- Helen Blade
- Oral Product Development, Pharmaceutical Technology & Development, Operations, AstraZeneca, Macclesfield SK10 2NA, U.K
| | | | - Steven P. Brown
- Department of Physics, University of Warwick, Coventry CV4 7AL, U.K
| | - Jake Carson
- Department of Statistics, University of Warwick, Coventry CV4 7AL, U.K
| | | | - Leslie P. Hughes
- Oral Product Development, Pharmaceutical Technology & Development, Operations, AstraZeneca, Macclesfield SK10 2NA, U.K
| | | |
Collapse
|
49
|
Fedorov AY, Rychkov DA. COMPARISON OF DIFFERENT COMPUTATIONAL APPROACHES FOR UNVEILING THE HIGH-PRESSURE BEHAVIOR OF ORGANIC CRYSTALS AT A MOLECULAR LEVEL. CASE STUDY OF TOLAZAMIDE POLYMORPHS. J STRUCT CHEM+ 2020. [DOI: 10.1134/s0022476620090024] [Citation(s) in RCA: 5] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/13/2022]
|
50
|
Taylor CR, Mulvee MT, Perenyi DS, Probert MR, Day GM, Steed JW. Minimizing Polymorphic Risk through Cooperative Computational and Experimental Exploration. J Am Chem Soc 2020; 142:16668-16680. [PMID: 32897065 PMCID: PMC7586337 DOI: 10.1021/jacs.0c06749] [Citation(s) in RCA: 20] [Impact Index Per Article: 5.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/16/2022]
Abstract
![]()
We
combine state-of-the-art computational crystal structure prediction
(CSP) techniques with a wide range of experimental crystallization
methods to understand and explore crystal structure in pharmaceuticals
and minimize the risk of unanticipated late-appearing polymorphs.
Initially, we demonstrate the power of CSP to rationalize the difficulty
in obtaining polymorphs of the well-known pharmaceutical isoniazid
and show that CSP provides the structure of the recently obtained,
but unsolved, Form III of this drug despite there being only a single
resolved form for almost 70 years. More dramatically, our blind CSP
study predicts a significant risk of polymorphism for the related
iproniazid. Employing a wide variety of experimental techniques, including
high-pressure experiments, we experimentally obtained the first three
known nonsolvated crystal forms of iproniazid, all of which were successfully
predicted in the CSP procedure. We demonstrate the power of CSP methods
and free energy calculations to rationalize the observed elusiveness
of the third form of iproniazid, the success of high-pressure experiments
in obtaining it, and the ability of our synergistic computational-experimental
approach to “de-risk” solid form landscapes.
Collapse
Affiliation(s)
- Christopher R Taylor
- Computational Systems Chemistry, School of Chemistry, University of Southampton, Southampton SO17 1NX, U.K
| | - Matthew T Mulvee
- Department of Chemistry, Durham University, South Road, Durham DH1 3LE, U.K
| | - Domonkos S Perenyi
- Department of Chemistry, Durham University, South Road, Durham DH1 3LE, U.K
| | - Michael R Probert
- Chemistry, School of Natural and Environmental Sciences, Newcastle University, Newcastle Upon Tyne NE1 7RU, U.K
| | - Graeme M Day
- Computational Systems Chemistry, School of Chemistry, University of Southampton, Southampton SO17 1NX, U.K
| | - Jonathan W Steed
- Department of Chemistry, Durham University, South Road, Durham DH1 3LE, U.K
| |
Collapse
|