1
|
Wu EJ, Kelly AW, Iuzzolino L, Lee AY, Zhu X. Unprecedented Packing Polymorphism of Oxindole: An Exploration Inspired by Crystal Structure Prediction. Angew Chem Int Ed Engl 2024; 63:e202406214. [PMID: 38825853 DOI: 10.1002/anie.202406214] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/01/2024] [Revised: 05/13/2024] [Accepted: 05/29/2024] [Indexed: 06/04/2024]
Abstract
Crystal polymorphism, characterized by different packing arrangements of the same compound, strongly ties to the physical properties of a molecule. Determining the polymorphic landscape is complex and time-consuming, with the number of experimentally observed polymorphs varying widely from molecule to molecule. Furthermore, disappearing polymorphs, the phenomenon whereby experimentally observed forms cannot be reproduced, pose a significant challenge for the pharmaceutical industry. Herein, we focused on oxindole (OX), a small rigid molecule with four known polymorphs, including a reported disappearing form. Using crystal structure prediction (CSP), we assessed OX solid-state landscape and thermodynamic stability by comparing predicted structures with experimentally known forms. We then performed melt and solution crystallization in bulk and nanoconfinement to validate our predictions. These experiments successfully reproduced the known forms and led to the discovery of four novel polymorphs. Our approach provided insights into reconstructing disappearing polymorphs and building more comprehensive polymorph landscapes. These results also establish a new record of packing polymorphism for rigid molecules.
Collapse
Affiliation(s)
- Emily J Wu
- Analytical Research & Development, Merck & Co., Inc., Rahway, New Jersey, 07065, United States
| | - Andrew W Kelly
- Analytical Research & Development, Merck & Co., Inc., Rahway, New Jersey, 07065, United States
| | - Luca Iuzzolino
- Modeling & Informatics, Discovery Chemistry, Merck & Co., Inc., Rahway, New Jersey, 07065, United States
| | - Alfred Y Lee
- Analytical Research & Development, Merck & Co., Inc., Rahway, New Jersey, 07065, United States
| | - Xiaolong Zhu
- Analytical Research & Development, Merck & Co., Inc., Rahway, New Jersey, 07065, United States
| |
Collapse
|
2
|
Nessler A, Okada O, Kinoshita Y, Nishimura K, Nagata H, Fukuzawa K, Yonemochi E, Schnieders MJ. Crystal Polymorph Search in the NPT Ensemble via a Deposition/Sublimation Alchemical Path. CRYSTAL GROWTH & DESIGN 2024; 24:3205-3217. [PMID: 38659664 PMCID: PMC11036363 DOI: 10.1021/acs.cgd.3c01358] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Figures] [Subscribe] [Scholar Register] [Received: 11/14/2023] [Revised: 02/22/2024] [Accepted: 02/23/2024] [Indexed: 04/26/2024]
Abstract
The formulation of active pharmaceutical ingredients involves discovering stable crystal packing arrangements or polymorphs, each of which has distinct pharmaceutically relevant properties. Traditional experimental screening techniques utilizing various conditions are commonly supplemented with in silico crystal structure prediction (CSP) to inform the crystallization process and mitigate risk. Predictions are often based on advanced classical force fields or quantum mechanical calculations that model the crystal potential energy landscape but do not fully incorporate temperature, pressure, or solution conditions during the search procedure. This study proposes an innovative alchemical path that utilizes an advanced polarizable atomic multipole force field to predict crystal structures based on direct sampling of the NPT ensemble. The use of alchemical (i.e., nonphysical) intermediates, a novel Monte Carlo barostat, and an orthogonal space tempering bias combine to enhance the sampling efficiency of the deposition/sublimation phase transition. The proposed algorithm was applied to 2-((4-(2-(3,4-dichlorophenyl)ethyl)phenyl)amino)benzoic acid (Cambridge Crystallography Database Centre ID: XAFPAY) as a case study to showcase the algorithm. Each experimentally determined polymorph with one molecule in the asymmetric unit was successfully reproduced via approximately 1000 short 1 ns simulations per space group where each simulation was initiated from random rigid body coordinates and unit cell parameters. Utilizing two threads of a recent Intel CPU (a Xeon Gold 6330 CPU at 2.00 GHz), 1 ns of sampling using the polarizable AMOEBA force field can be acquired in 4 h (equating to more than 300 ns/day using all 112 threads/56 cores of a dual CPU node) within the Force Field X software (https://ffx.biochem.uiowa.edu). These results demonstrate a step forward in the rigorous use of the NPT ensemble during the CSP search process and open the door to future algorithms that incorporate solution conditions using continuum solvation methods.
Collapse
Affiliation(s)
- Aaron
J. Nessler
- Department
of Biomedical Engineering, University of
Iowa, 103 South Capitol
Street, 5601 Seamans Center for the Engineering Arts and Sciences, Iowa City, Iowa 52242, United States
| | - Okimasa Okada
- Sohyaku
Innovative Research Division, Mitsubishi
Tanabe Pharma Corporation, 1000 Kamoshida-cho, Aoba-ku, Yokohama, Kanagawa 227-0033, Japan
| | - Yuya Kinoshita
- Analytical
Development, Pharmaceutical Sciences, Takeda
Pharmaceutical Company Limited, 2-26-1, Muraoka-Higashi, Fujisawa 251-8555, Kanagawa, Japan
| | - Koki Nishimura
- Analytical
Development, Pharmaceutical Sciences, Takeda
Pharmaceutical Company Limited, 2-26-1, Muraoka-Higashi, Fujisawa 251-8555, Kanagawa, Japan
| | - Hiroomi Nagata
- CMC
Modality Technology Laboratories, Production Technology and Supply
Chain Management Division, Mitsubishi Tanabe
Pharma Corporation, Osaka 541-8505, Japan
| | - Kaori Fukuzawa
- Graduate
School of Pharmaceutical Sciences, Osaka
University, 1-6 Yamadaoka, Suita, Osaka 565-0871, Japan
| | - Etsuo Yonemochi
- Department
of Physical Chemistry, School of Pharmacy and Pharmaceutical Sciences, Hoshi University, 2-4-41 Ebara, Shinagawa-ku, Tokyo 142-8501, Japan
| | - Michael J. Schnieders
- Department
of Biomedical Engineering, University of
Iowa, 103 South Capitol
Street, 5601 Seamans Center for the Engineering Arts and Sciences, Iowa City, Iowa 52242, United States
- Department
of Biochemistry, University of Iowa, 51 Newton Road, 4-403 Bowen Science
Building, Iowa City, Iowa 52242, United States
| |
Collapse
|
3
|
Shishkina SV, Shaposhnyk AM, Konovalova IS, Dyakonenko VV, Vaksler YO. Concomitant polymorphs of 2-imino-2H-chromene-3-carboxylic acid amide: experimental and quantum chemical study. ACTA CRYSTALLOGRAPHICA SECTION B, STRUCTURAL SCIENCE, CRYSTAL ENGINEERING AND MATERIALS 2024; 80:27-37. [PMID: 38205837 DOI: 10.1107/s2052520623010193] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Received: 09/29/2023] [Accepted: 11/24/2023] [Indexed: 01/12/2024]
Abstract
2-Iminocoumarin-3-carboxamide (2-imino-2 H-chromene-3-carboxylic acid) is a perspective compound for use in the pharmaceutical industry. This compound crystallized from several solvents as two concomitant polymorphic forms. The monoclinic polymorph, crystallized initially, is formed due to strong N-H...O hydrogen bonds, weak C-H...O and C-H...N(π) hydrogen bonds, and stacking interactions of `head-to-head' type. The triclinic polymorphic form obtained due to slow evaporation of the same solution is formed due to only strong intermolecular interactions, N-H...O hydrogen bonds of two types, and stacking interactions of two types. Analysis of pairwise interaction energies showed that the monoclinic structure is columnar while the triclinic one is layered. Calculations in a periodic approximation of their lattice energies confirmed that the monoclinic polymorphic crystals are metastable as compared to the stable triclinic polymorph. Further quantum chemical modeling of possible structure deformations proved that both concomitant polymorphs can not be transformed into a new polymorphic form under external influence.
Collapse
Affiliation(s)
- Svitlana V Shishkina
- SSI "Institute for Single Crystals" NAS of Ukraine, 60 Nauky ave., Kharkiv, 61001, Ukraine
| | - Anna M Shaposhnyk
- SSI "Institute for Single Crystals" NAS of Ukraine, 60 Nauky ave., Kharkiv, 61001, Ukraine
| | - Irina S Konovalova
- SSI "Institute for Single Crystals" NAS of Ukraine, 60 Nauky ave., Kharkiv, 61001, Ukraine
| | - Victoriya V Dyakonenko
- SSI "Institute for Single Crystals" NAS of Ukraine, 60 Nauky ave., Kharkiv, 61001, Ukraine
| | - Yevhenii O Vaksler
- SSI "Institute for Single Crystals" NAS of Ukraine, 60 Nauky ave., Kharkiv, 61001, Ukraine
| |
Collapse
|
4
|
Hoja J, List A, Boese AD. Multimer Embedding Approach for Molecular Crystals up to Harmonic Vibrational Properties. J Chem Theory Comput 2024; 20:357-367. [PMID: 38109226 PMCID: PMC10782452 DOI: 10.1021/acs.jctc.3c01082] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/29/2023] [Revised: 11/16/2023] [Accepted: 11/20/2023] [Indexed: 12/20/2023]
Abstract
Accurate calculations of molecular crystals are crucial for drug design and crystal engineering. However, periodic high-level density functional calculations using hybrid functionals are often prohibitively expensive for the relevant systems. These expensive periodic calculations can be circumvented by the usage of embedding methods in which, for instance, the periodic calculation is only performed at a lower-cost level and then monomer energies and dimer interactions are replaced by those of the higher-level method. Herein, we extend such a multimer embedding approach to enable energy corrections for trimer interactions and the calculation of harmonic vibrational properties up to the dimer level. We evaluate this approach for the X23 benchmark set of molecular crystals by approximating a periodic hybrid density functional (PBE0+MBD) by embedding multimers into less expensive calculations using a generalized-gradient approximation functional (PBE+MBD). We show that trimer interactions are crucial for accurately approximating lattice energies within 1 kJ/mol and might also be needed for further improvement of lattice constants and hence cell volumes. Finally, the vibrational properties are already very well captured at the monomer and dimer level, making it possible to approximate vibrational free energies at room temperature within 1 kJ/mol.
Collapse
Affiliation(s)
- Johannes Hoja
- Department of Chemistry, University
of Graz, Heinrichstraße 28/IV, Graz 8010, Austria
| | - Alexander List
- Department of Chemistry, University
of Graz, Heinrichstraße 28/IV, Graz 8010, Austria
| | - A. Daniel Boese
- Department of Chemistry, University
of Graz, Heinrichstraße 28/IV, Graz 8010, Austria
| |
Collapse
|
5
|
Mathew R, Mazumder A, Kumar P, Matula J, Mohamed S, Brazda P, Hariharan M, Thomas B. Unveiling the topology of partially disordered micro-crystalline nitro-perylenediimide with X-aggregate stacking: an integrated approach. Chem Sci 2024; 15:490-499. [PMID: 38179523 PMCID: PMC10762722 DOI: 10.1039/d3sc05514k] [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: 10/17/2023] [Accepted: 11/23/2023] [Indexed: 01/06/2024] Open
Abstract
Profound knowledge of the molecular structure and supramolecular organization of organic molecules is essential to understand their structure-property relationships. Herein we demonstrate the packing arrangement of partially disordered nitro-perylenediimide (NO2-PDI), revealing that the perylenediimide units exhibit an X-shaped packing pattern. The packing of NO2-PDI is derived using a complementary approach that utilises solid-state NMR (ssNMR) and 3D electron diffraction (3D ED) techniques. Perylenediimide (PDI) molecules are captivating due to their high luminescence efficiency and optoelectronic properties, which are related to supramolecular self-assembly. Increasing the alkyl chain length on the imide substituent poses a more significant challenge in crystallizing the resulting molecule. In addition to the alkyl tails, other functional groups, like the nitro group attached as a bay substituent, can also cause disorder. Such heterogeneity could lead to diffuse scattering, which then complicates the interpretation of diffraction experiment data, where perfect periodicity is expected. As a result, there is an unmet need to develop a methodology for solving the structures of difficult-to-crystallize materials. A synergistic approach is utilised in this manuscript to understand the packing arrangement of the disordered material NO2-PDI by making use of 3D ED, ssNMR and density functional theory calculations (DFT). The combination of these experimental and theoretical approaches provides great promise in enabling the structural investigation of novel materials with customized properties across various applications, which are, due to the internal disorder, very difficult to study by diffraction techniques. By effectively addressing these challenges, our methodology opens up new avenues for material characterization, thereby driving exciting advancements in the field.
Collapse
Affiliation(s)
- Renny Mathew
- Science Division, New York University Abu Dhabi P.O. Box 129188 Abu Dhabi United Arab Emirates
| | - Aniruddha Mazumder
- School of Chemistry, Indian Institute of Science Education and Research Thiruvananthapuram (IISER TVM) Maruthamala P.O., Vithura Thiruvananthapuram 695551 Kerala India
| | - Praveen Kumar
- Science Division, New York University Abu Dhabi P.O. Box 129188 Abu Dhabi United Arab Emirates
| | - Julie Matula
- Science Division, New York University Abu Dhabi P.O. Box 129188 Abu Dhabi United Arab Emirates
| | - Sharmarke Mohamed
- Department of Chemistry, Green Chemistry & Materials Modelling Laboratory, Khalifa University of Science and Technology P.O. Box 127788 Abu Dhabi United Arab Emirates
- Advanced Materials Chemistry Center (AMCC), Khalifa University of Science and Technology P.O. Box 127788 Abu Dhabi United Arab Emirates
| | - Petr Brazda
- Institute of Physics of the Czech Academy of Sciences Na Slovance 2/1999 18200 Prague 8 Czech Republic
| | - Mahesh Hariharan
- School of Chemistry, Indian Institute of Science Education and Research Thiruvananthapuram (IISER TVM) Maruthamala P.O., Vithura Thiruvananthapuram 695551 Kerala India
| | - Brijith Thomas
- Science Division, New York University Abu Dhabi P.O. Box 129188 Abu Dhabi United Arab Emirates
| |
Collapse
|
6
|
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
|
7
|
Hermann J, Stöhr M, Góger S, Chaudhuri S, Aradi B, Maurer RJ, Tkatchenko A. libMBD: A general-purpose package for scalable quantum many-body dispersion calculations. J Chem Phys 2023; 159:174802. [PMID: 37933783 DOI: 10.1063/5.0170972] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/05/2023] [Accepted: 10/17/2023] [Indexed: 11/08/2023] Open
Abstract
Many-body dispersion (MBD) is a powerful framework to treat van der Waals (vdW) dispersion interactions in density-functional theory and related atomistic modeling methods. Several independent implementations of MBD with varying degree of functionality exist across a number of electronic structure codes, which both limits the current users of those codes and complicates dissemination of new variants of MBD. Here, we develop and document libMBD, a library implementation of MBD that is functionally complete, efficient, easy to integrate with any electronic structure code, and already integrated in FHI-aims, DFTB+, VASP, Q-Chem, CASTEP, and Quantum ESPRESSO. libMBD is written in modern Fortran with bindings to C and Python, uses MPI/ScaLAPACK for parallelization, and implements MBD for both finite and periodic systems, with analytical gradients with respect to all input parameters. The computational cost has asymptotic cubic scaling with system size, and evaluation of gradients only changes the prefactor of the scaling law, with libMBD exhibiting strong scaling up to 256 processor cores. Other MBD properties beyond energy and gradients can be calculated with libMBD, such as the charge-density polarization, first-order Coulomb correction, the dielectric function, or the order-by-order expansion of the energy in the dipole interaction. Calculations on supramolecular complexes with MBD-corrected electronic structure methods and a meta-review of previous applications of MBD demonstrate the broad applicability of the libMBD package to treat vdW interactions.
Collapse
Affiliation(s)
- Jan Hermann
- Department of Mathematics and Computer Science, FU Berlin, 14195 Berlin, Germany
| | - Martin Stöhr
- Department of Physics and Materials Science, University of Luxembourg, L-1511 Luxembourg City, Luxembourg
| | - Szabolcs Góger
- Department of Physics and Materials Science, University of Luxembourg, L-1511 Luxembourg City, Luxembourg
| | - Shayantan Chaudhuri
- Department of Chemistry, University of Warwick, Coventry CV4 7AL, United Kingdom
| | - Bálint Aradi
- Bremen Center for Computational Materials Science, University of Bremen, 28359 Bremen, Germany
| | - Reinhard J Maurer
- Department of Chemistry, University of Warwick, Coventry CV4 7AL, United Kingdom
- Department of Physics, University of Warwick, Coventry CV4 7AL, United Kingdom
| | - Alexandre Tkatchenko
- Department of Physics and Materials Science, University of Luxembourg, L-1511 Luxembourg City, Luxembourg
| |
Collapse
|
8
|
Li S, Liu B, Chen Z, Ou X, Rong H, Lu M. Ritonavir Revisited: Melt Crystallization Can Easily Find the Late-Appearing Polymorph II and Unexpectedly Discover a New Polymorph III. Mol Pharm 2023; 20:3854-3863. [PMID: 37450774 DOI: 10.1021/acs.molpharmaceut.2c00994] [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] [Indexed: 07/18/2023]
Abstract
Identification of a thermodynamically stable polymorph is an important step in the early stage of drug development. Ritonavir (RIT) is a well-known case where the most stable polymorph II emerged after being marketed, leading to a loss of $250 million. Herein, we report the findings that routine melt crystallization can reveal the late-appearing polymorph II of RIT at small supercooling, but the probability of nucleation is very low. The addition of 30-50% polyethylene glycol (PEG) promotes the crystallization of Form II as the only phase at low supercooling, making it easier to detect in polymorphism screening. During the course of our research, a new polymorph, denoted Form III, was unexpectedly discovered, crystallizing as the major phase from neat RIT melts. Single crystals of Form III were grown from melt microdroplets. Benefiting from the ability of synchrotron radiation to detect weak diffraction signals that cannot be accessible by a laboratory diffractometer, a reasonable structure of Form III was solved with slight disorder relative to thiazole groups (P1 space group and Z' = 4). The thermodynamic stability ranking of the three true polymorphs is Form II > Form I > Form III, as opposed to the order of solubility. The capacity to efficiently reveal rich polymorphs, especially the kinetically hindered polymorph, and rapidly grow single crystals of a new phase for structure determination together highlights the necessity of incorporating melt crystallization into routine methods for pharmaceutical polymorphism screening.
Collapse
Affiliation(s)
- Shuting Li
- School of Pharmaceutical Sciences, Sun Yat-sen University, Guangzhou 510006, China
| | - Binbin Liu
- School of Pharmaceutical Sciences, Sun Yat-sen University, Guangzhou 510006, China
| | - Ziqiao Chen
- School of Pharmaceutical Sciences, Sun Yat-sen University, Guangzhou 510006, China
| | - Xiao Ou
- School of Pharmaceutical Sciences, Sun Yat-sen University, Guangzhou 510006, China
| | - Haowei Rong
- School of Pharmaceutical Sciences, Sun Yat-sen University, Guangzhou 510006, China
| | - Ming Lu
- School of Pharmaceutical Sciences, Sun Yat-sen University, Guangzhou 510006, China
| |
Collapse
|
9
|
Abramov YA, Iuzzolino L, Jin Y, York G, Chen CH, Shultz CS, Yang Z, Chang C, Shi B, Zhou T, Greenwell C, Sekharan S, Lee AY. Cocrystal Synthesis through Crystal Structure Prediction. Mol Pharm 2023. [PMID: 37279175 DOI: 10.1021/acs.molpharmaceut.2c01098] [Citation(s) in RCA: 2] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 06/08/2023]
Abstract
Crystal structure prediction (CSP) is an invaluable tool in the pharmaceutical industry because it allows to predict all the possible crystalline solid forms of small-molecule active pharmaceutical ingredients. We have used a CSP-based cocrystal prediction method to rank ten potential cocrystal coformers by the energy of the cocrystallization reaction with an antiviral drug candidate, MK-8876, and a triol process intermediate, 2-ethynylglyclerol. For MK-8876, the CSP-based cocrystal prediction was performed retrospectively and successfully predicted the maleic acid cocrystal as the most likely cocrystal to be observed. The triol is known to form two different cocrystals with 1,4-diazabicyclo[2.2.2]octane (DABCO), but a larger solid form landscape was desired. CSP-based cocrystal screening predicted the triol-DABCO cocrystal as rank one, while a triol-l-proline cocrystal was predicted as rank two. Computational finite-temperature corrections enabled determination of relative crystallization propensities of the triol-DABCO cocrystals with different stoichiometries and prediction of the triol-l-proline polymorphs in the free-energy landscape. The triol-l-proline cocrystal was obtained during subsequent targeted cocrystallization experiments and was found to exhibit an improved melting point and deliquescence behavior over the triol-free acid, which could be considered as an alternative solid form in the synthesis of islatravir.
Collapse
Affiliation(s)
- Yuriy A Abramov
- XtalPi Inc., 245 Main Street, Cambridge, Massachusetts 02142, United States
- Eshelman School of Pharmacy, University of North Carolina, Chapel Hill, North Carolina 27599, United States
| | - Luca Iuzzolino
- Computational and Structural Chemistry, Merck & Co., Inc., 126 East Lincoln Avenue, Rahway, New Jersey 07065, United States
| | - Yingdi Jin
- XtalPi, Inc., Shenzhen Jingtai Technology Co., Ltd., Futian District, Shenzhen 518100, China
| | - Gregory York
- Analytical Research and Development, Merck & Co., Inc., 126 East Lincoln Avenue, Rahway, New Jersey 07065, United States
| | - Chien-Hung Chen
- Analytical Research and Development, Merck & Co., Inc., 126 East Lincoln Avenue, Rahway, New Jersey 07065, United States
| | - C Scott Shultz
- Analytical Research and Development, Merck & Co., Inc., 126 East Lincoln Avenue, Rahway, New Jersey 07065, United States
| | - Zhuocen Yang
- XtalPi, Inc., Shenzhen Jingtai Technology Co., Ltd., Futian District, Shenzhen 518100, China
| | - Chao Chang
- XtalPi, Inc., Shenzhen Jingtai Technology Co., Ltd., Futian District, Shenzhen 518100, China
| | - Baimei Shi
- XtalPi, Inc., Shenzhen Jingtai Technology Co., Ltd., Futian District, Shenzhen 518100, China
| | - Tian Zhou
- XtalPi, Inc., Shenzhen Jingtai Technology Co., Ltd., Futian District, Shenzhen 518100, China
| | - Chandler Greenwell
- XtalPi Inc., 245 Main Street, Cambridge, Massachusetts 02142, United States
| | - Sivakumar Sekharan
- XtalPi Inc., 245 Main Street, Cambridge, Massachusetts 02142, United States
| | - Alfred Y Lee
- Analytical Research and Development, Merck & Co., Inc., 126 East Lincoln Avenue, Rahway, New Jersey 07065, 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
|
Chen B, Xu X. Discriminating and understanding molecular crystal polymorphism. J Comput Chem 2023; 44:969-979. [PMID: 36585855 DOI: 10.1002/jcc.27057] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/04/2022] [Revised: 11/01/2022] [Accepted: 11/30/2022] [Indexed: 01/01/2023]
Abstract
Polymorph discrimination for a molecular crystal has long been a challenging task, which, nonetheless, is a major concern in the pharmaceutical industry. In this work, we have investigated polymorph discrimination on three different molecular crystals, tetrolic acid, oxalic acid, and oxalyl dihydrazide, covering both packing polymorphism and conformational polymorphism. To gain more understanding, we have performed energy decomposition analysis based on many-body expansion, and have compared the results from the XO-PBC method, that is, the eXtended ONIOM method (XO) with the periodic boundary condition (PBC), with those from some commonly used dispersion corrected density functional theory (DFT-D) methods. It is shown here that, with the XYG3 doubly hybrid functional chosen as the target high level to capture the intra- and short-range intermolecular interactions, and the periodic PBE as the basic low level to take long range interactions into account, the XO-PBC(XYG3:PBE) method not only obtains the correct experimental stability orderings, but also predicts reasonable polymorph energy ranges for all three cases. Our results have demonstrated the usefulness of the present theoretical methods, in particular XO-PBC, while highlighted the importance of a better treatment of different kinds of interactions to be beneficial to polymorph control.
Collapse
Affiliation(s)
- Bozhu Chen
- Department of Chemistry, Collaborative Innovation Center of Chemistry for Energy Materials, Shanghai Key Laboratory of Molecular Catalysis and Innovative Materials, Ministry of Education Key Laboratory of Computational Physical Sciences, Fudan University, Shanghai, China
| | - Xin Xu
- Department of Chemistry, Collaborative Innovation Center of Chemistry for Energy Materials, Shanghai Key Laboratory of Molecular Catalysis and Innovative Materials, Ministry of Education Key Laboratory of Computational Physical Sciences, Fudan University, Shanghai, China.,Hefei National Laboratory, Hefei, China
| |
Collapse
|
12
|
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
|
13
|
Jung DŠ, Lukin S, Halasz I. Improving the Accuracy of Small‐Molecule Crystal Structures Solved from Powder X‐Ray Diffraction Data by Using External Sources. Helv Chim Acta 2023. [DOI: 10.1002/hlca.202200087] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/09/2023]
Affiliation(s)
| | - Stipe Lukin
- Ruđer Bošković Institute Bijenička 54 HR-10000 Zagreb Croatia
| | - Ivan Halasz
- Ruđer Bošković Institute Bijenička 54 HR-10000 Zagreb Croatia
| |
Collapse
|
14
|
Price AJA, Otero-de-la-Roza A, Johnson ER. XDM-corrected hybrid DFT with numerical atomic orbitals predicts molecular crystal lattice energies with unprecedented accuracy. Chem Sci 2023; 14:1252-1262. [PMID: 36756332 PMCID: PMC9891363 DOI: 10.1039/d2sc05997e] [Citation(s) in RCA: 10] [Impact Index Per Article: 10.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/29/2022] [Accepted: 12/13/2022] [Indexed: 12/23/2022] Open
Abstract
Molecular crystals are important for many applications, including energetic materials, organic semiconductors, and the development and commercialization of pharmaceuticals. The exchange-hole dipole moment (XDM) dispersion model has shown good performance in the calculation of relative and absolute lattice energies of molecular crystals, although it has traditionally been applied in combination with plane-wave/pseudopotential approaches. This has limited XDM to use with semilocal functional approximations, which suffer from delocalization error and poor quality conformational energies, and to systems with a few hundreds of atoms at most due to unfavorable scaling. In this work, we combine XDM with numerical atomic orbitals, which enable the efficient use of XDM-corrected hybrid functionals for molecular crystals. We test the new XDM-corrected functionals for their ability to predict the lattice energies of molecular crystals for the X23 set and 13 ice phases, the latter being a particularly stringent test. A composite approach using a XDM-corrected, 25% hybrid functional based on B86bPBE achieves a mean absolute error of 0.48 kcal mol-1 per molecule for the X23 set and 0.19 kcal mol-1 for the total lattice energies of the ice phases, compared to recent diffusion Monte-Carlo data. These results make the new XDM-corrected hybrids not only far more computationally efficient than previous XDM implementations, but also the most accurate density-functional methods for molecular crystal lattice energies to date.
Collapse
Affiliation(s)
- Alastair J. A. Price
- Department of Chemistry, Dalhousie University6274 Coburg RdHalifaxB3H 4R2Nova ScotiaCanada
| | - Alberto Otero-de-la-Roza
- Departamento de Química Física y Analítica and MALTA-Consolider Team, Facultad de Química, Universidad de Oviedo Oviedo 33006 Spain
| | - Erin R. Johnson
- Department of Chemistry, Dalhousie University6274 Coburg RdHalifaxB3H 4R2Nova ScotiaCanada
| |
Collapse
|
15
|
Li S, Lightowler M, Ou X, Huang S, Jiang Y, Li X, Zou X, Xu H, Lu M. Direct structure determination of vemurafenib polymorphism from compact spherulites using 3D electron diffraction. Commun Chem 2023; 6:18. [PMID: 36697943 PMCID: PMC9871043 DOI: 10.1038/s42004-022-00804-2] [Citation(s) in RCA: 6] [Impact Index Per Article: 6.0] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/05/2022] [Accepted: 12/22/2022] [Indexed: 01/25/2023] Open
Abstract
The spherulitic morphology is considered to be the most common morphology of crystalline materials and is particularly apparent in melt-crystallized products. Yet, historically, the polycrystalline nature of spherulites has hindered successful crystal structure determination. Here, we report the direct structure determination of a clinical drug, vemurafenib (VMN), in compact spherulite form using 3D electron diffraction (3D ED). VMN has four known polymorphs. We first solved the crystal structures of α-, β-, and γ-VMN from compact spherulites using 3D ED, and the resulting structures were highly consistent with those obtained by single-crystal X-ray diffraction. We then determined the crystal structure of δ-VMN-the least stable polymorph which cannot be cultivated as a single crystal-directly from the compact spherulite sample. We unexpectedly discovered a new polymorph during our studies, denoted as ε-VMN. Single crystals of ε-VMN are extremely thin and not suitable for study by X-ray diffraction. Again, we determined the structure of ε-VMN in a compact spherulite form. This successful structure elucidation of all five VMN polymorphs demonstrates the possibility of directly determining structures from melt-grown compact spherulite samples. Thereby, this discovery will improve the efficiency and broaden the scope of polymorphism research, especially within the field of melt crystallization.
Collapse
Affiliation(s)
- Shuting Li
- School of Pharmaceutical Sciences, Sun Yat-sen University, Guangzhou, China
| | - Molly Lightowler
- Department of Materials and Environmental Chemistry, Stockholm University, Stockholm, Sweden
| | - Xiao Ou
- School of Pharmaceutical Sciences, Sun Yat-sen University, Guangzhou, China
| | - Siyong Huang
- School of Pharmaceutical Sciences, Sun Yat-sen University, Guangzhou, China
| | - Yifan Jiang
- School of Pharmaceutical Sciences, Sun Yat-sen University, Guangzhou, China
| | - Xizhen Li
- School of Pharmaceutical Sciences, Sun Yat-sen University, Guangzhou, China
| | - Xiaodong Zou
- Department of Materials and Environmental Chemistry, Stockholm University, Stockholm, Sweden
| | - Hongyi Xu
- Department of Materials and Environmental Chemistry, Stockholm University, Stockholm, Sweden.
| | - Ming Lu
- School of Pharmaceutical Sciences, Sun Yat-sen University, Guangzhou, China.
| |
Collapse
|
16
|
Saha BK, Nath NK, Thakuria R. Polymorphs with Remarkably Distinct Physical and/or Chemical Properties. CHEM REC 2023; 23:e202200173. [PMID: 36166697 DOI: 10.1002/tcr.202200173] [Citation(s) in RCA: 1] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/06/2022] [Revised: 08/30/2022] [Indexed: 01/21/2023]
Abstract
Polymorphism in crystals is known since 1822 and the credit goes to Mitscherlich who realized the existence of different crystal structures of the same compound while working with some arsenate and phosphate salts. Later on, this phenomenon was observed also in organic crystals. With the advent of different technologies, especially the easy availability of single crystal XRD instruments, polymorphism in crystals has become a common phenomenon. Almost 37 % of compounds (single component) are polymorphic to date. As the energies of the different polymorphic forms are very close to each other, small changes in crystallization conditions might lead to different polymorphic structures. As a result, sometimes it is difficult to control polymorphism. For this reason, it is considered to be a nuisance to crystal engineering. It has been realized that the property of a material depends not only on the molecular structure but also on its crystal structure. Therefore, it is not only of interest to academia but also has widespread applications in the materials science as well as pharmaceutical industries. In this review, we have discussed polymorphism which causes significant changes in materials properties in different fields of solid-state science, such as electrical, magnetic, SHG, thermal expansion, mechanical, luminescence, color, and pharmaceutical. Therefore, this review will interest researchers from supramolecular chemistry, materials science as well as medicinal chemistry.
Collapse
Affiliation(s)
- Binoy K Saha
- Department of Chemistry, Pondicherry University, Puducherry, 605014, India
| | - Naba K Nath
- Department of Chemistry, National Institute of Technology Meghalaya, Shillong, Meghalaya 793003, India
| | - Ranjit Thakuria
- Department of Chemistry, Gauhati University, Guwahati, 781014, India
| |
Collapse
|
17
|
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: 2] [Impact Index Per Article: 1.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
|
18
|
Application of functionalized magnetic silica nanoparticles for selective induction of three coumarin metastable polymorphs. Chin J Chem Eng 2022. [DOI: 10.1016/j.cjche.2022.05.008] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/17/2022]
|
19
|
Smalley CJH, Hoskyns HE, Hughes CE, Johnstone DN, Willhammar T, Young MT, Pickard CJ, Logsdail AJ, Midgley PA, Harris KDM. A structure determination protocol based on combined analysis of 3D-ED data, powder XRD data, solid-state NMR data and DFT-D calculations reveals the structure of a new polymorph of l-tyrosine. Chem Sci 2022; 13:5277-5288. [PMID: 35655549 PMCID: PMC9093151 DOI: 10.1039/d1sc06467c] [Citation(s) in RCA: 9] [Impact Index Per Article: 4.5] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/19/2021] [Accepted: 03/29/2022] [Indexed: 11/21/2022] Open
Abstract
We report the crystal structure of a new polymorph of l-tyrosine (denoted the β polymorph), prepared by crystallization from the gas phase following vacuum sublimation. Structure determination was carried out by combined analysis of three-dimensional electron diffraction (3D-ED) data and powder X-ray diffraction (XRD) data. Specifically, 3D-ED data were required for reliable unit cell determination and space group assignment, with structure solution carried out independently from both 3D-ED data and powder XRD data, using the direct-space strategy for structure solution implemented using a genetic algorithm. Structure refinement was carried out both from powder XRD data, using the Rietveld profile refinement technique, and from 3D-ED data. The final refined structure was validated both by periodic DFT-D calculations, which confirm that the structure corresponds to an energy minimum on the energy landscape, and by the fact that the values of isotropic 13C NMR chemical shifts calculated for the crystal structure using DFT-D methodology are in good agreement with the experimental high-resolution solid-state 13C NMR spectrum. Based on DFT-D calculations using the PBE0-MBD method, the β polymorph is meta-stable with respect to the previously reported crystal structure of l-tyrosine (now denoted the α polymorph). Crystal structure prediction calculations using the AIRSS approach suggest that there are three other plausible crystalline polymorphs of l-tyrosine, with higher energy than the α and β polymorphs.
Collapse
Affiliation(s)
| | - Harriet E Hoskyns
- School of Chemistry, Cardiff University Park Place Cardiff CF10 3AT Wales UK
| | - Colan E Hughes
- School of Chemistry, Cardiff University Park Place Cardiff CF10 3AT Wales UK
| | - Duncan N Johnstone
- Department of Materials Science, University of Cambridge 27 Charles Babbage Road Cambridge CB3 0FS England UK
| | - Tom Willhammar
- Department of Materials and Environmental Chemistry, Stockholm University Svante Arrhenius väg 16C 106 91 Stockholm Sweden
| | - Mark T Young
- School of Biosciences, Cardiff University Cardiff CF10 3AX Wales UK
| | - Christopher J Pickard
- Department of Materials Science, University of Cambridge 27 Charles Babbage Road Cambridge CB3 0FS England UK
- Advanced Institute for Materials Research, Tohoku University 2-1-1 Katahira Aoba Sendai 980-8577 Japan
| | - Andrew J Logsdail
- Cardiff Catalysis Institute, School of Chemistry, Cardiff University Park Place Cardiff CF10 3AT Wales UK
| | - Paul A Midgley
- Department of Materials Science, University of Cambridge 27 Charles Babbage Road Cambridge CB3 0FS England UK
| | - Kenneth D M Harris
- School of Chemistry, Cardiff University Park Place Cardiff CF10 3AT Wales UK
| |
Collapse
|
20
|
Pham CH, Lindsey RK, Fried LE, Goldman N. High-Accuracy Semiempirical Quantum Models Based on a Minimal Training Set. J Phys Chem Lett 2022; 13:2934-2942. [PMID: 35343698 DOI: 10.1021/acs.jpclett.2c00453] [Citation(s) in RCA: 5] [Impact Index Per Article: 2.5] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 06/14/2023]
Abstract
A great need exists for computationally efficient quantum simulation approaches that can achieve an accuracy similar to high-level theories at a fraction of the computational cost. In this regard, we have leveraged a machine-learned interaction potential based on Chebyshev polynomials to improve density functional tight binding (DFTB) models for organic materials. The benefit of our approach is two-fold: (1) many-body interactions can be corrected for in a systematic and rapidly tunable process, and (2) high-level quantum accuracy for a broad range of compounds can be achieved with ∼0.3% of data required for one advanced deep learning potential. Our model exhibits both transferability and extensibility through comparison to quantum chemical results for organic clusters, solid carbon phases, and molecular crystal phase stability rankings. Our efforts thus allow for high-throughput physical and chemical predictions with up to coupled-cluster accuracy for systems that are computationally intractable with standard approaches.
Collapse
Affiliation(s)
- Cong Huy Pham
- Physical and Life Sciences Directorate, Lawrence Livermore National Laboratory, Livermore, California 94550, United States
| | - Rebecca K Lindsey
- Physical and Life Sciences Directorate, Lawrence Livermore National Laboratory, Livermore, California 94550, United States
| | - Laurence E Fried
- Physical and Life Sciences Directorate, Lawrence Livermore National Laboratory, Livermore, California 94550, United States
| | - Nir Goldman
- Physical and Life Sciences Directorate, Lawrence Livermore National Laboratory, Livermore, California 94550, United States
- Department of Chemical Engineering, University of California, Davis, California 95616, United States
| |
Collapse
|
21
|
Ramuthai M, Jeyavijayan S, Premkumar R, Uma Priya M, Jayram ND. Structure, Spectroscopic Investigation, Molecular Docking and In vitro Cytotoxicity Studies on 4,7-dihydroxycoumarin: A Breast Cancer Drug. JOURNAL OF COMPUTATIONAL BIOPHYSICS AND CHEMISTRY 2022. [DOI: 10.1142/s2737416522500119] [Citation(s) in RCA: 2] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 11/18/2022]
Abstract
Coumarin derivatives are broadly used as anti-inflammatory, antioxidants, anticancer, and antiviral drugs in recent years. In particular, hydroxy coumarins have great importance because of their various biological and pharmacological purposes. The quantum chemical studies of 4,7-dihydroxycoumarin (DHC) have been performed using the cc-pVTZ level of basis set. The DHC molecular structure has been optimized and the computed frequency assignments have been correlated well with the experimental results. The experimental [Formula: see text]C NMR shifts of DHC have been compared with the computed [Formula: see text]C NMR in the dimethyl sulfoxide (DMSO) solution using the Gauge-invariant atomic orbital (GIAO) method. The electron delocalization within the DHC is shown by highest occupied molecular orbitals (HOMO)-lowest unoccupied molecular orbitals (LUMO) energy analysis, and the resulting small energy gap value reveal the molecule’s bioactive characteristics. The natural bond orbital (NBO) analysis approves the bioactive property of the DHC molecule. The DHC compound has a cytotoxic impact on the MCF-7 breast cancer cell line, according to in vitro cytotoxicity studies. The docking study approves that the DHC works as a new inhibitor of breast cancer targeted proteins such as epidermal growth factor receptor (EGFR), estrogen receptor (ER), and progesterone receptor (PR). Thus, this work covers the approach for the evolution of new drugs against breast cancer.
Collapse
Affiliation(s)
- M. Ramuthai
- Department of Physics, Kalasalingam Academy of Research and Education Krishnankoil, 626126 Tamil Nadu, India
| | - S. Jeyavijayan
- Department of Physics, Kalasalingam Academy of Research and Education Krishnankoil, 626126 Tamil Nadu, India
| | - R. Premkumar
- PG and Research Department of Physics, N.M.S.S.V.N. College Madurai, 625019 Tamil Nadu, India
| | - M. Uma Priya
- Department of Biotechnology, Kalasalingam Academy of Research and Education Krishnankoil 626126 Tamil Nadu, India
| | - Naidu Dhanpal Jayram
- Department of Physics, Kalasalingam Academy of Research and Education Krishnankoil, 626126 Tamil Nadu, India
| |
Collapse
|
22
|
Thandra DR, Allikayala R. Synthesis, characterization, molecular structure determination by single crystal X-ray diffraction, and Hirshfeld surface analysis of 7-fluoro-6-morpholino-3-phenylquinolin-1-ium chloride salt and computational studies of its cation. J Mol Struct 2022. [DOI: 10.1016/j.molstruc.2021.131701] [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]
|
23
|
Smalley CH, Logsdail AJ, Hughes CE, Iuga D, Young MT, Harris KDM. Solid-State Structural Properties of Alloxazine Determined from Powder XRD Data in Conjunction with DFT-D Calculations and Solid-State NMR Spectroscopy: Unraveling the Tautomeric Identity and Pathways for Tautomeric Interconversion. CRYSTAL GROWTH & DESIGN 2022; 22:524-534. [PMID: 35024003 PMCID: PMC8739831 DOI: 10.1021/acs.cgd.1c01114] [Citation(s) in RCA: 4] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Received: 09/25/2021] [Revised: 10/25/2021] [Indexed: 06/02/2023]
Abstract
We report the solid-state structural properties of alloxazine, a tricyclic ring system found in many biologically important molecules, with structure determination carried out directly from powder X-ray diffraction (XRD) data. As the crystal structures containing the alloxazine and isoalloxazine tautomers both give a high-quality fit to the powder XRD data in Rietveld refinement, other techniques are required to establish the tautomeric form in the solid state. In particular, high-resolution solid-state 15N NMR data support the presence of the alloxazine tautomer, based on comparison between isotropic chemical shifts in the experimental 15N NMR spectrum and the corresponding values calculated for the crystal structures containing the alloxazine and isoalloxazine tautomers. Furthermore, periodic DFT-D calculations at the PBE0-MBD level indicate that the crystal structure containing the alloxazine tautomer has significantly lower energy. We also report computational investigations of the interconversion between the tautomeric forms in the crystal structure via proton transfer along two intermolecular N-H···N hydrogen bonds; DFT-D calculations at the PBE0-MBD level indicate that the tautomeric interconversion is associated with a lower energy transition state for a mechanism involving concerted (rather than sequential) proton transfer along the two hydrogen bonds. However, based on the relative energies of the crystal structures containing the alloxazine and isoalloxazine tautomers, it is estimated that under conditions of thermal equilibrium at ambient temperature, more than 99.9% of the molecules in the crystal structure will exist as the alloxazine tautomer.
Collapse
Affiliation(s)
| | - Andrew J. Logsdail
- Cardiff
Catalysis Institute, School of Chemistry, Cardiff University, Park Place, Cardiff CF10
3AT, Wales, United Kingdom
| | - Colan E. Hughes
- School
of Chemistry, Cardiff University, Park Place, Cardiff CF10 3AT, Wales,
United Kingdom
| | - Dinu Iuga
- Department
of Physics, University of Warwick, Coventry CV4 7AL, England, United Kingdom
| | - Mark T. Young
- School
of Biosciences, Cardiff University, Cardiff CF10 3AX, Wales, United Kingdom
| | - Kenneth D. M. Harris
- School
of Chemistry, Cardiff University, Park Place, Cardiff CF10 3AT, Wales,
United Kingdom
| |
Collapse
|
24
|
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
|
25
|
Zhu X, Hu CT, Erriah B, Vogt-Maranto L, Yang J, Yang Y, Qiu M, Fellah N, Tuckerman ME, Ward MD, Kahr B. Imidacloprid Crystal Polymorphs for Disease Vector Control and Pollinator Protection. J Am Chem Soc 2021; 143:17144-17152. [PMID: 34634905 DOI: 10.1021/jacs.1c07610] [Citation(s) in RCA: 11] [Impact Index Per Article: 3.7] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Abstract
Imidacloprid, the world's leading insecticide, has been approved recently for controlling infectious disease vectors; yet, in agricultural settings, it has been implicated in the frightening decline of pollinators. This argues for strategies that sharply reduce the environmental impact of imidacloprid. When used as a contact insecticide, the effectiveness of imidacloprid relies on physical contact between its crystal surfaces and insect tarsi. Herein, seven new imidacloprid crystal polymorphs are reported, adding to two known forms. Anticipating that insect uptake of imidacloprid molecules would depend on the respective free energies of crystal polymorph surfaces, measurements of insect knockdown times for the metastable crystal forms were as much as nine times faster acting than the commercial form against Aedes, Anopheles, and Culex mosquitoes as well as Drosophila (fruit flies). These results suggest that replacement of commercially available imidacloprid crystals (a.k.a. Form I) in space-spraying with any one of three new polymorphs, Forms IV, VI, IX, would suppress vector-borne disease transmission while reducing environmental exposure and harm to nontarget organisms.
Collapse
Affiliation(s)
- Xiaolong Zhu
- Department of Chemistry and Molecular Design Institute, New York University, New York, New York 10003 United States
| | - Chunhua T Hu
- Department of Chemistry and Molecular Design Institute, New York University, New York, New York 10003 United States
| | - Bryan Erriah
- Department of Chemistry and Molecular Design Institute, 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
| | - Jingxiang Yang
- Department of Chemistry and Molecular Design Institute, New York University, New York, New York 10003 United States
| | - Yongfan Yang
- Department of Chemistry and Molecular Design Institute, New York University, New York, New York 10003 United States
| | - Mengdi Qiu
- Department of Chemistry and Molecular Design Institute, New York University, New York, New York 10003 United States
| | - Noalle Fellah
- Department of Chemistry and Molecular Design Institute, New York University, New York, New York 10003 United States
| | - Mark E Tuckerman
- Department of Chemistry and Molecular Design Institute, 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, New York University Shanghai, Shanghai 200062, China
| | - Michael D Ward
- Department of Chemistry and Molecular Design Institute, New York University, New York, New York 10003 United States
| | - Bart Kahr
- Department of Chemistry and Molecular Design Institute, New York University, New York, New York 10003 United States
| |
Collapse
|
26
|
Wang H, Schrage BR, Takematsu K, Ziegler CJ. Photophysical properties of a boron analogue of coumarin. Phys Chem Chem Phys 2021; 23:18855-18862. [PMID: 34612423 DOI: 10.1039/d1cp02386a] [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
In this report, we present a study into the structure and electronic properties of difluoroboronsalicylaldoxime (DFBS), a boron-based structural analog of coumarin. The modification of the heterocyclic ring of coumarin with boron results in a compound with similar structural parameters and molecular orbitals to coumarin. DFT and TDDFT calculations reveal a significant stabilization of the LUMO in DFBS; this is supported by a ∼40 nm red shift of the lowest electronic transition in the absorption spectrum. Interestingly, DFBS is emissive, while unmodified coumarin is effectively non-radiative. Comparisons between DFBS, emissive coumarin variants, and unmodified coumarin suggest that the charge transfer character of the transition contributes to the fluorescence.
Collapse
Affiliation(s)
- Huayi Wang
- Department of Chemistry, University of Akron, Akron, Ohio 44325-3601, USA.
| | | | | | | |
Collapse
|
27
|
Surov AO, Voronin AP, Drozd KV, Gruzdev MS, Perlovich GL, Prashanth J, Balasubramanian S. Polymorphic forms of antiandrogenic drug nilutamide: structural and thermodynamic aspects. Phys Chem Chem Phys 2021; 23:9695-9708. [PMID: 33908506 DOI: 10.1039/d1cp00793a] [Citation(s) in RCA: 5] [Impact Index Per Article: 1.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022]
Abstract
Attempts to obtain new cocrystals of nonsteroidal antiandrogenic drug nilutamide produced alternative polymorphic forms of the compound (Form II and Form III) and their crystal structures were elucidated by single-crystal X-ray diffraction. Apart from the cocrystallization technique, lyophilization was found to be an effective strategy for achieving polymorph control of nilutamide, which was difficult to obtain by other methods. The physicochemical properties and relative stability of the commercial Form I and newly obtained Form II were comprehensively investigated by a variety of analytical methods (thermal analysis, solution calorimetry, solubility, and sublimation), whereas for Form III, only a handful of experimental parameters were obtained due to the elusive nature of the polymorph. Form I and Form II were found to be monotropically related, with Form I being confirmed as the thermodynamically most stable solid phase. In addition, the performance of different DFT-D and semi-empirical schemes for lattice energy calculation and polymorph energy ranking was compared and analysed. Lattice energy calculations using periodic DFT at B3LYP-D3/6-31(F+)G(d,p) and PBEh-3c/def2-mSVP levels of theory were found to provide the most accurate lattice energy values for Form I against experimental data, while PIXEL and PBEh-3c/def2-mSVP were the only methods that predicted the correct order of stability of Forms I and II.
Collapse
Affiliation(s)
- Artem O Surov
- G. A. Krestov Institute of Solution Chemistry of the Russian Academy of Sciences, 1 Akademicheskaya St., 153045 Ivanovo, Russia.
| | | | | | | | | | | | | |
Collapse
|
28
|
Demir S, Tekin A. FFCASP: A Massively Parallel Crystal Structure Prediction Algorithm. J Chem Theory Comput 2021; 17:2586-2598. [PMID: 33798330 DOI: 10.1021/acs.jctc.0c01197] [Citation(s) in RCA: 5] [Impact Index Per Article: 1.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/23/2022]
Abstract
A new algorithm called Fast and Flexible CrystAl Structure Predictor (FFCASP) was developed to predict the structure of covalent and molecular crystals. FFCASP is massively parallel and able to handle more than 200 atoms in the unit cell (in other terms, it allows global optimization around 100 individual parameters). It uses a global optimizer specialized for Crystal Structure Prediction (CSP) which combines particle swarm and simulated annealing optimizers. Three different molecular crystals, including diverse intermolecular interactions, namely, cytosine, coumarin, and pyrazinamide, have been selected to evaluate the performance of FFCASP. While cytosine polymorphs have been searched by employing two different force fields (a DFT-SAPT based intermolecular potential and generalized amber force field (GAFF)) up to Z = 16, only GAFF has been used both in coumarin and pyrazinamide polymorph searches up to Z = 4. For these three molecular crystals, FFCASP generated more than 20 000 crystal structures, and the unique ones have been further treated by DFT-D3. A combination of data mining and a machine learning approach was introduced to determine the unique structures and their distribution into different clusters, which ultimately gives an opportunity to retrieve the common features and relations between the resulting structures. There are two known experimental crystal structures of cytosine, and both were successfully located with FFCASP. Two of the reported crystal structures of coumarin have been reproduced. Similarly, in pyrazinamide, three known experimental structures have been rediscovered. In addition to finding the experimentally known structures, FFCASP also located other low-energy structures for each considered molecular crystals. These successes of FFCASP offer the possibility to discover the polymorphic nature of other important molecular crystals (e.g., drugs) as well.
Collapse
Affiliation(s)
- Samet Demir
- Informatics Institute, Istanbul Technical University, 34469 Maslak, Istanbul, Turkey.,TÜBİTAK Research Institute for Fundamental Sciences, 41470 Gebze, Kocaeli, Turkey
| | - Adem Tekin
- Informatics Institute, Istanbul Technical University, 34469 Maslak, Istanbul, Turkey.,TÜBİTAK Research Institute for Fundamental Sciences, 41470 Gebze, Kocaeli, Turkey
| |
Collapse
|
29
|
Kostritskiy AY, Dmitriev MV, Grinev VS, Fedotova OV. CRYSTAL STRUCTURE AND PACKING FEATURES
OF 3-(5-METHYL-1H-PYRAZOL-3-YL)-2H-CHROMEN-2-ONE
AND 3-(3-METHYL-1H-PYRAZOL-3-YL)-2H-CHROMEN-2-ONE. J STRUCT CHEM+ 2021. [DOI: 10.1134/s0022476621030112] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/23/2022]
|
30
|
Wengert S, Csányi G, Reuter K, Margraf JT. Data-efficient machine learning for molecular crystal structure prediction. Chem Sci 2021; 12:4536-4546. [PMID: 34163719 PMCID: PMC8179468 DOI: 10.1039/d0sc05765g] [Citation(s) in RCA: 28] [Impact Index Per Article: 9.3] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/19/2020] [Accepted: 02/05/2021] [Indexed: 12/16/2022] Open
Abstract
The combination of modern machine learning (ML) approaches with high-quality data from quantum mechanical (QM) calculations can yield models with an unrivalled accuracy/cost ratio. However, such methods are ultimately limited by the computational effort required to produce the reference data. In particular, reference calculations for periodic systems with many atoms can become prohibitively expensive for higher levels of theory. This trade-off is critical in the context of organic crystal structure prediction (CSP). Here, a data-efficient ML approach would be highly desirable, since screening a huge space of possible polymorphs in a narrow energy range requires the assessment of a large number of trial structures with high accuracy. In this contribution, we present tailored Δ-ML models that allow screening a wide range of crystal candidates while adequately describing the subtle interplay between intermolecular interactions such as H-bonding and many-body dispersion effects. This is achieved by enhancing a physics-based description of long-range interactions at the density functional tight binding (DFTB) level-for which an efficient implementation is available-with a short-range ML model trained on high-quality first-principles reference data. The presented workflow is broadly applicable to different molecular materials, without the need for a single periodic calculation at the reference level of theory. We show that this even allows the use of wavefunction methods in CSP.
Collapse
Affiliation(s)
- Simon Wengert
- Chair of Theoretical Chemistry, Technische Universität München 85747 Garching Germany
| | - Gábor Csányi
- Engineering Laboratory, University of Cambridge Cambridge CB2 1PZ UK
| | - Karsten Reuter
- Chair of Theoretical Chemistry, Technische Universität München 85747 Garching Germany
- Fritz-Haber-Institut der Max-Planck-Gesellschaft Faradayweg 4-6 14195 Berlin Germany
| | - Johannes T Margraf
- Chair of Theoretical Chemistry, Technische Universität München 85747 Garching Germany
- Fritz-Haber-Institut der Max-Planck-Gesellschaft Faradayweg 4-6 14195 Berlin Germany
| |
Collapse
|
31
|
Shishkina SV, Baumer VN, Kovalenko SM, Trostianko PV, Bunyatyan ND. Usage of Quantum Chemical Methods to Understand the Formation of Concomitant Polymorphs of Acetyl 2-( N-(2-Fluorophenyl)imino)coumarin-3-carboxamide. ACS OMEGA 2021; 6:3120-3129. [PMID: 33553928 PMCID: PMC7860053 DOI: 10.1021/acsomega.0c05516] [Citation(s) in RCA: 5] [Impact Index Per Article: 1.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Figures] [Subscribe] [Scholar Register] [Received: 11/11/2020] [Accepted: 01/12/2021] [Indexed: 06/12/2023]
Abstract
Crystallization of concomitant polymorphs is a very intriguing process that is difficult to be studied experimentally. A comprehensive study of two polymorphic modifications of acetyl 2-(N-(2-fluorophenyl)imino)coumarin-3-carboxamide using quantum chemical methods has revealed molecular and crystal structure dependence on crystallization conditions. Fast crystallization associated with a kinetically controlled process results in the formation of a columnar structure with a nonequilibrium molecular conformation and more isotropic topology of interaction energies between molecules. Slow crystallization may be considered as a thermodynamically controlled process and leads to the formation of a layered crystal structure with the conformation of the molecule corresponding to local minima and anisotropic topology of interaction energies. Fast crystallization results in the formation of a lot of weak intermolecular interactions, while slow crystallization leads to the formation of small amounts of stronger interactions.
Collapse
Affiliation(s)
- Svitlana V. Shishkina
- Department
of X-ray Diffraction Study and Quantum Chemistry, SSI Institute for Single Crystals NAS of Ukraine, 60 Nauky Avenue, Kharkiv 61001, Ukraine
- Department
of Inorganic Chemistry, V. N. Karazin Kharkiv
National University, 4 Svobodi sq., Kharkiv 61022, Ukraine
| | - Vyacheslav N. Baumer
- Department
of X-ray Diffraction Study and Quantum Chemistry, SSI Institute for Single Crystals NAS of Ukraine, 60 Nauky Avenue, Kharkiv 61001, Ukraine
| | - Sergiy M. Kovalenko
- Department
of Organic Chemistry, V. N. Karazin Kharkiv
National University, 4 Svobodi sq., Kharkiv 61022, Ukraine
- I.M.
Sechenov First Moscow State Medical University, 8 Trubeckaya Str., Moscow 119991, Russia
| | - Pavel V. Trostianko
- Department
of Organic Chemistry, V. N. Karazin Kharkiv
National University, 4 Svobodi sq., Kharkiv 61022, Ukraine
| | - Natalya D. Bunyatyan
- I.M.
Sechenov First Moscow State Medical University, 8 Trubeckaya Str., Moscow 119991, Russia
- Federal
State Budgetary Institution Scientific Centre for Expert Evaluation
of Medicinal Products of the Ministry of Health of the Russian Federation, 8 Petrovsky blv., Moscow 127051, Russia
| |
Collapse
|
32
|
Rodríguez I, Gautam R, Tinoco AD. Using X-ray Diffraction Techniques for Biomimetic Drug Development, Formulation, and Polymorphic Characterization. Biomimetics (Basel) 2020; 6:1. [PMID: 33396786 PMCID: PMC7838816 DOI: 10.3390/biomimetics6010001] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.8] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/19/2020] [Revised: 12/19/2020] [Accepted: 12/23/2020] [Indexed: 12/31/2022] Open
Abstract
Drug development is a decades-long, multibillion dollar investment that often limits itself. To decrease the time to drug approval, efforts are focused on drug targets and drug formulation for optimal biocompatibility and efficacy. X-ray structural characterization approaches have catalyzed the drug discovery and design process. Single crystal X-ray diffraction (SCXRD) reveals important structural details and molecular interactions for the manifestation of a disease or for therapeutic effect. Powder X-ray diffraction (PXRD) has provided a method to determine the different phases, purity, and stability of biological drug compounds that possess crystallinity. Recently, synchrotron sources have enabled wider access to the study of noncrystalline or amorphous solids. One valuable technique employed to determine atomic arrangements and local atom ordering of amorphous materials is the pair distribution function (PDF). PDF has been used in the study of amorphous solid dispersions (ASDs). ASDs are made up of an active pharmaceutical ingredient (API) within a drug dispersed at the molecular level in an amorphous polymeric carrier. This information is vital for appropriate formulation of a drug for stability, administration, and efficacy purposes. Natural or biomimetic products are often used as the API or the formulation agent. This review profiles the deep insights that X-ray structural techniques and associated analytical methods can offer in the development of a drug.
Collapse
Affiliation(s)
- Israel Rodríguez
- Department of Chemistry, University of Puerto Rico Río Piedras, San Juan, PR 00925, USA
| | - Ritika Gautam
- Department of Chemistry, Indian Institute of Technology Kanpur, Kanpur 208016, India
| | - Arthur D. Tinoco
- Department of Chemistry, University of Puerto Rico Río Piedras, San Juan, PR 00925, USA
| |
Collapse
|
33
|
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] [Grants] [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
- School of Pharmaceutical Sciences, Sun Yat-sen University, Guangzhou, China
| | - Xiao Ou
- School of Pharmaceutical Sciences, Sun Yat-sen University, Guangzhou, China
| | - Bingquan Wang
- School of Pharmaceutical Sciences, Sun Yat-sen University, Guangzhou, China
| | - Haowei Rong
- School 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
- School of Pharmacy, University of Wisconsin - Madison, Madison, WI, USA
| | - Ming Lu
- School of Pharmaceutical Sciences, Sun Yat-sen University, Guangzhou, China.
- Guangdong Provincial Key Laboratory of New Drug Design and Evaluation, Sun Yat-sen University, Guangzhou, China.
| |
Collapse
|
34
|
Abraham NS, Shirts MR. Statistical Mechanical Approximations to More Efficiently Determine Polymorph Free Energy Differences for Small Organic Molecules. J Chem Theory Comput 2020; 16:6503-6512. [PMID: 32877183 DOI: 10.1021/acs.jctc.0c00570] [Citation(s) in RCA: 8] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
Abstract
Methods to efficiently determine the relative stability of polymorphs of organic crystals are highly desired in crystal structure predictions (CSPs). Current methodologies include calculating the free energy of static lattice phonons, quasi-harmonic approximations (QHA), and computing the full thermodynamic cycle using replica exchange molecular dynamics (REMD). We found that 13 out of the 29 systems minimized from experimental crystal structures restructured to a lower energy minimum when heated and annealed using REMD, a phenomenon that QHA alone cannot capture. Here, we present a series of methods that are intermediate in accuracy and expense between QHA and computing the full thermodynamic cycle, which can save 42-80% of the computational cost and introduces, on this benchmark, a relatively small (0.16 ± 0.04 kcal/mol) error relative to the full thermodynamic cycle. In particular, a method that Boltzmann weights harmonic free energies from along the trajectory of REMD replica exchange appears to be an appropriate intermediate between QHA and the full thermodynamic cycle using MD when screening crystal polymorph stability.
Collapse
Affiliation(s)
- Nathan S Abraham
- Department of Chemical and Biological Engineering, University of Colorado Boulder, Boulder, Colorado 80309, United States
| | - Michael R Shirts
- Department of Chemical and Biological Engineering, University of Colorado Boulder, Boulder, Colorado 80309, United States
| |
Collapse
|
35
|
Sun X, Pradeepkumar P, Rajendran NK, Shakila H, Houreld NN, Al Farraj DA, Elnahas YM, Elumalai N, Rajan M. Natural deep eutectic solvent supported targeted solid-liquid polymer carrier for breast cancer therapy. RSC Adv 2020; 10:36989-37004. [PMID: 35521273 PMCID: PMC9057073 DOI: 10.1039/d0ra03790g] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/27/2020] [Accepted: 08/11/2020] [Indexed: 11/21/2022] Open
Abstract
Solid–liquid nanocarriers (SLNs) are at the front of the rapidly emerging field of medicinal applications with a potential role in the delivery of bioactive agents. Here, we report a new SLN of natural deep eutectic solvent (NADES) and biotin-conjugated lysine–polyethylene glycol copolymer. The SLN system was analyzed for its functional groups, thermal stability, crystalline nature, particle size, and surface morphology through the instrumental analysis of FT-IR, TGA, XRD, DLS, SEM, and TEM. Encapsulation of PTX (paclitaxel) and 7-HC (7-hydroxycoumarin) with the SLN was carried out by dialysis, and UV-visible spectra evidenced the drug loading capacity and higher encapsulation efficiency obtained. The enhanced anticancer potential of PTX- and 7-HC-loaded SLN was assessed in vitro, and the system reduces the cell viability of MDA-MB-231 cells. The PTX- and 7-HC-loaded SLN system was investigated in a breast cancer-induced rat model via in vivo studies. It shows decreased lysosomal enzymes and increased levels of caspase to cure breast tumors. It very well may be reasoned that the designed PTX- and 7-HC-loaded SLN system has strong anticancer properties and exhibits potential for delivery of drug molecules in cancer treatment. Solid–liquid nanocarriers (SLNs) are at the front of the rapidly emerging field of medicinal applications with a potential role in the delivery of bioactive agents.![]()
Collapse
Affiliation(s)
- Xianfu Sun
- Department of Breast, The Affiliated Tumor Hospital of Zhengzhou University Zhengzhou Henan 450008 China
| | - Periyakaruppan Pradeepkumar
- Biomaterials in Medicinal Chemistry Laboratory, Department of Natural Products Chemistry, School of Chemistry, Madurai Kamaraj University Madurai-625021 Tamil Nadu India
| | - Naresh Kumar Rajendran
- Laser Research Centre, Faculty of Health Sciences, University of Johannesburg PO Box 17011 Doornfontein 2028 South Africa
| | - Harshavardhan Shakila
- Department of Molecular Microbiology, School of Biotechnology, Madurai Kamaraj University Madurai-625021 India
| | - Nicolette Nadene Houreld
- Laser Research Centre, Faculty of Health Sciences, University of Johannesburg PO Box 17011 Doornfontein 2028 South Africa
| | - Dunia A Al Farraj
- Department of Botany and Microbiology, College of Science, King Saud University Riyadh 11451 Saudi Arabia
| | - Yousif M Elnahas
- Department of Surgery, College of Medicine, King Saud University Medical City, King Saud University Riyadh Saudi Arabia
| | - Nandhakumar Elumalai
- Department of Biochemistry, Sri Muthukumaran Medical College and Research Institute Chennai-600069 Tamil Nadu India
| | - Mariappan Rajan
- Biomaterials in Medicinal Chemistry Laboratory, Department of Natural Products Chemistry, School of Chemistry, Madurai Kamaraj University Madurai-625021 Tamil Nadu India
| |
Collapse
|
36
|
Chen B, Xu X. XO-PBC: An Accurate and Efficient Method for Molecular Crystals. J Chem Theory Comput 2020; 16:4271-4285. [PMID: 32456429 DOI: 10.1021/acs.jctc.0c00232] [Citation(s) in RCA: 8] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Abstract
In this work, we propose the XO-PBC method, which combines the eXtended ONIOM method (XO) with the periodic boundary condition (PBC) for the description of molecular crystals. XO-PBC tries to embed a finite cluster cut out from the solid into the periodic environment, making it feasible to employ advanced molecular quantum chemistry methods, which are usually prohibitively expensive for direct PBC calculations. In particular, XO-PBC utilizes the results from force calculations to design the scheme to fragment the molecule when crystals are made of large molecules and to select cluster model systems automatically consisting of dimer up to tetramer interactions for embedding. By applying an appropriate theory to each model, a satisfactory accuracy for the system under study is ensured, while a high efficiency is achieved with massively parallel computing by distributing model systems onto different processors. A comparison of the XO-PBC calculations with the conventional direct PBC calculations at the B3LYP level demonstrates its accuracy at substantially low cost for the description of molecular crystals. The usefulness of the XO-PBC method is further exemplified, showing that XO-PBC is able to predict the lattice energies of various types of molecular crystals within chemical accuracy (<4 kJ/mol) when the doubly hybrid density functional XYG3 is used as the target high level and the periodic PBE as the basic low level. The XO-PBC method provides a general protocol that brings the great predictive power of advanced electronic structure methods from molecular systems to the extended solids.
Collapse
Affiliation(s)
- Bozhu Chen
- Collaborative Innovation Center of Chemistry for Energy Materials, Shanghai Key Laboratory of Molecular Catalysis and Innovative Materials, Ministry of Education Key Laboratory of Computational Physical Sciences, Department of Chemistry, Fudan University, Shanghai, 200433, China
| | - Xin Xu
- Collaborative Innovation Center of Chemistry for Energy Materials, Shanghai Key Laboratory of Molecular Catalysis and Innovative Materials, Ministry of Education Key Laboratory of Computational Physical Sciences, Department of Chemistry, Fudan University, Shanghai, 200433, China
| |
Collapse
|
37
|
Greenwell C, McKinley JL, Zhang P, Zeng Q, Sun G, Li B, Wen S, Beran GJO. Overcoming the difficulties of predicting conformational polymorph energetics in molecular crystals via correlated wavefunction methods. Chem Sci 2020; 11:2200-2214. [PMID: 32190277 PMCID: PMC7059316 DOI: 10.1039/c9sc05689k] [Citation(s) in RCA: 38] [Impact Index Per Article: 9.5] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/09/2019] [Accepted: 01/13/2020] [Indexed: 11/21/2022] Open
Abstract
Molecular crystal structure prediction is increasingly being applied to study the solid form landscapes of larger, more flexible pharmaceutical molecules. Despite many successes in crystal structure prediction, van der Waals-inclusive density functional theory (DFT) methods exhibit serious failures predicting the polymorph stabilities for a number of systems exhibiting conformational polymorphism, where changes in intramolecular conformation lead to different intermolecular crystal packings. Here, the stabilities of the conformational polymorphs of o-acetamidobenzamide, ROY, and oxalyl dihydrazide are examined in detail. DFT functionals that have previously been very successful in crystal structure prediction perform poorly in all three systems, due primarily to the poor intramolecular conformational energies, but also due to the intermolecular description in oxalyl dihydrazide. In all three cases, a fragment-based dispersion-corrected second-order Møller-Plesset perturbation theory (MP2D) treatment of the crystals overcomes these difficulties and predicts conformational polymorph stabilities in good agreement with experiment. These results highlight the need for methods which go beyond current-generation DFT functionals to make crystal polymorph stability predictions truly reliable.
Collapse
Affiliation(s)
- Chandler Greenwell
- Department of Chemistry , University of California , Riverside , California 92521 , USA . ; Tel: +1-951-827-7869
| | - Jessica L McKinley
- Department of Chemistry , University of California , Riverside , California 92521 , USA . ; Tel: +1-951-827-7869
| | - Peiyu Zhang
- Xtalpi, Inc. , 245 Main St, 12th Floor , Cambridge , MA 02142 , USA
| | - Qun Zeng
- Xtalpi, Inc. , 245 Main St, 12th Floor , Cambridge , MA 02142 , USA
| | - Guangxu Sun
- Xtalpi, Inc. , 245 Main St, 12th Floor , Cambridge , MA 02142 , USA
| | - Bochen Li
- Xtalpi, Inc. , 245 Main St, 12th Floor , Cambridge , MA 02142 , USA
| | - Shuhao Wen
- Xtalpi, Inc. , 245 Main St, 12th Floor , Cambridge , MA 02142 , USA
| | - Gregory J O Beran
- Department of Chemistry , University of California , Riverside , California 92521 , USA . ; Tel: +1-951-827-7869
| |
Collapse
|
38
|
Ishii H, Obata S, Niitsu N, Watanabe S, Goto H, Hirose K, Kobayashi N, Okamoto T, Takeya J. Charge mobility calculation of organic semiconductors without use of experimental single-crystal data. Sci Rep 2020; 10:2524. [PMID: 32066751 PMCID: PMC7026405 DOI: 10.1038/s41598-020-59238-2] [Citation(s) in RCA: 8] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/29/2019] [Accepted: 01/24/2020] [Indexed: 01/23/2023] Open
Abstract
Prediction of material properties of newly designed molecules is a long-term goal in organic electronics. In general, it is a difficult problem, because the material properties are dominated by the unknown packing structure. We present a practical method to obtain charge transport properties of organic single crystals, without use of experimental single-crystal data. As a demonstration, we employ the promising molecule C10-DNBDT. We succeeded in quantitative evaluation of charge mobility of the single crystal using our quantum wave-packet dynamical simulation method. Here, the single-crystal data is computationally obtained by searching possible packing structures from structural formula of the molecule. We increase accuracy in identifying the actual crystal structure from suggested ones by using not only crystal energy but also similarity between calculated and experimental powder X-ray diffraction patterns. The proposed methodology can be a theoretical design technique for efficiently developing new high-performance organic semiconductors, since it can estimate the charge transport properties at early stage in the process of material development.
Collapse
Affiliation(s)
- Hiroyuki Ishii
- Department of Applied Physics, Faculty of Pure and Applied Sciences, University of Tsukuba, 1-1-1 Tennodai, Tsukuba, Ibaraki, 305-8573, Japan.
| | - Shigeaki Obata
- Educational Programs on Advanced Simulation Engineering, Toyohashi University of Technology, 1-1 Hibarigaoka, Tempaku-cho, Toyohashi, Aichi, 441-8580, Japan.
- CONFLEX Corporation, Shinagawa Center Bldg. 6F, 3-23-17 Takanawa, Minato-ku, Tokyo, 108-0074, Japan.
| | - Naoyuki Niitsu
- Material Innovation Research Center (MIRC) and Department of Advanced Materials Science, Graduate School of Frontier Sciences, The University of Tokyo, 5-1-5 Kashiwanoha, Kashiwa, Chiba, 277-8561, Japan
| | - Shun Watanabe
- Material Innovation Research Center (MIRC) and Department of Advanced Materials Science, Graduate School of Frontier Sciences, The University of Tokyo, 5-1-5 Kashiwanoha, Kashiwa, Chiba, 277-8561, Japan
- JST, PRESTO, 4-1-8 Honcho, Kawaguchi, Saitama, 332-0012, Japan
| | - Hitoshi Goto
- Educational Programs on Advanced Simulation Engineering, Toyohashi University of Technology, 1-1 Hibarigaoka, Tempaku-cho, Toyohashi, Aichi, 441-8580, Japan
- CONFLEX Corporation, Shinagawa Center Bldg. 6F, 3-23-17 Takanawa, Minato-ku, Tokyo, 108-0074, Japan
- Department of Computer Science and Engineering, Toyohashi University of Technology, 1-1 Hibarigaoka, Tempaku-cho, Toyohashi, Aichi, 441-8580, Japan
| | - Kenji Hirose
- Department of Applied Physics, Faculty of Pure and Applied Sciences, University of Tsukuba, 1-1-1 Tennodai, Tsukuba, Ibaraki, 305-8573, Japan
| | - Nobuhiko Kobayashi
- Department of Applied Physics, Faculty of Pure and Applied Sciences, University of Tsukuba, 1-1-1 Tennodai, Tsukuba, Ibaraki, 305-8573, Japan
| | - Toshihiro Okamoto
- Material Innovation Research Center (MIRC) and Department of Advanced Materials Science, Graduate School of Frontier Sciences, The University of Tokyo, 5-1-5 Kashiwanoha, Kashiwa, Chiba, 277-8561, Japan
- JST, PRESTO, 4-1-8 Honcho, Kawaguchi, Saitama, 332-0012, Japan
| | - Jun Takeya
- Material Innovation Research Center (MIRC) and Department of Advanced Materials Science, Graduate School of Frontier Sciences, The University of Tokyo, 5-1-5 Kashiwanoha, Kashiwa, Chiba, 277-8561, Japan
- International Center of Materials Nanoarchitectonics, National Institute for Materials Science (NIMS), 1-1 Namiki, Tsukuba, 305-0044, Japan
| |
Collapse
|
39
|
Ortiz KS, Hernández Espinell JR, Ortiz Torres D, Lopéz-Mejías V, Stelzer T. Polymorphism in Solid Dispersions. CRYSTAL GROWTH & DESIGN 2020; 20:713-722. [PMID: 38107251 PMCID: PMC10723824 DOI: 10.1021/acs.cgd.9b01138] [Citation(s) in RCA: 4] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 12/19/2023]
Abstract
Solid dispersions embed active pharmaceutical ingredients in polymeric carriers to improve their solubility. Three solid dispersion preparation techniques are typically employed: solvent evaporation, solvent-fusion, and fusion methods. Although these are also widely recommended as preparative methods for phase diagram determination, few examples exist concerning their effect on the resulting polymorph, once the solid dispersion is produced. In this study, the influence of these methods on the polymorphic form obtained in crystalline solid dispersions (CSDs) composed of flufenamic acid (FFA) and poly(ethylene glycol) was investigated. The physical mixtures and CSDs were characterized by powder X-ray diffraction, infrared spectroscopy, and differential scanning calorimetry. The results reveal that the fusion method leads to concomitant polymorphs (mainly FFA I and III) in the CSDs. In contrast, the solvent evaporation and solvent-fusion methods lead to FFA III. Collectively, these results demonstrate that preparative methods have a significant influence on the phase diagrams determined (average relative deviation ≤8%), which are often used to justify the design space of manufacturing processes, including those deemed "continuous." Consequently, choosing a preparation method that results in the desired polymorph is crucial to ensure accurate determination of phase diagrams and critical quality attributes of formulations.
Collapse
Affiliation(s)
- Karina Sanabria Ortiz
- Department of Chemistry, University of Puerto Rico, Río Piedras Campus, San Juan, Puerto Rico 00931, United States
- Crystallization Design Institute, Molecular Sciences Research Center, University of Puerto Rico, San Juan, Puerto Rico 00926, United States
| | - José R. Hernández Espinell
- Department of Chemistry, University of Puerto Rico, Río Piedras Campus, San Juan, Puerto Rico 00931, United States
- Crystallization Design Institute, Molecular Sciences Research Center, University of Puerto Rico, San Juan, Puerto Rico 00926, United States
| | - Desire Ortiz Torres
- Department of Chemistry, University of Puerto Rico, Río Piedras Campus, San Juan, Puerto Rico 00931, United States
- Crystallization Design Institute, Molecular Sciences Research Center, University of Puerto Rico, San Juan, Puerto Rico 00926, United States
| | - Vilmalí Lopéz-Mejías
- Department of Chemistry, University of Puerto Rico, Río Piedras Campus, San Juan, Puerto Rico 00931, United States
- Crystallization Design Institute, Molecular Sciences Research Center, University of Puerto Rico, San Juan, Puerto Rico 00926, United States
| | - Torsten Stelzer
- Crystallization Design Institute, Molecular Sciences Research Center, University of Puerto Rico, San Juan, Puerto Rico 00926, United States
- Department of Pharmaceutical Sciences, University of Puerto Rico, Medical Sciences Campus, San Juan, Puerto Rico 00936, United States
| |
Collapse
|
40
|
Zhang K, Fellah N, Shtukenberg AG, Fu X, Hu C, Ward MD. Discovery of new polymorphs of the tuberculosis drug isoniazid. CrystEngComm 2020. [DOI: 10.1039/d0ce00440e] [Citation(s) in RCA: 18] [Impact Index Per Article: 4.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022]
Abstract
Two new metastable polymorphs of the tuberculosis drug isoniazid, considered monomorphic for sixty years, were discovered using melt crystallization and nanoscale confinement.
Collapse
Affiliation(s)
- Keke Zhang
- Department of Chemistry and Molecular Design Institute
- New York University
- New York City
- USA
- School of Chemical Engineering and Technology
| | - Noalle Fellah
- Department of Chemistry and Molecular Design Institute
- New York University
- New York City
- USA
| | | | - Xiaoyan Fu
- State Key Laboratory of Catalysis
- Dalian Institute of Chemical Physics
- Chinese Academy of Sciences
- Dalian 116023
- People's Republic of China
| | - Chunhua Hu
- Department of Chemistry and Molecular Design Institute
- New York University
- New York City
- USA
| | - Michael D. Ward
- Department of Chemistry and Molecular Design Institute
- New York University
- New York City
- USA
| |
Collapse
|
41
|
Al Rahal O, Hughes CE, Williams PA, Logsdail AJ, Diskin‐Posner Y, Harris KDM. Polymorphism of
l
‐Tryptophan. Angew Chem Int Ed Engl 2019; 58:18788-18792. [DOI: 10.1002/anie.201908247] [Citation(s) in RCA: 16] [Impact Index Per Article: 3.2] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/03/2019] [Revised: 10/11/2019] [Indexed: 11/10/2022]
Affiliation(s)
- Okba Al Rahal
- School of ChemistryCardiff University Park Place Cardiff Wales CF10 3AT UK
| | - Colan E. Hughes
- School of ChemistryCardiff University Park Place Cardiff Wales CF10 3AT UK
| | - P. Andrew Williams
- School of ChemistryCardiff University Park Place Cardiff Wales CF10 3AT UK
| | - Andrew J. Logsdail
- School of ChemistryCardiff University Park Place Cardiff Wales CF10 3AT UK
| | - Yael Diskin‐Posner
- Department of Chemical Research SupportWeizmann Institute of Science Rehovot 76100 Israel
| | | |
Collapse
|
42
|
Al Rahal O, Hughes CE, Williams PA, Logsdail AJ, Diskin‐Posner Y, Harris KDM. Polymorphism of
l
‐Tryptophan. Angew Chem Int Ed Engl 2019. [DOI: 10.1002/ange.201908247] [Citation(s) in RCA: 5] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/12/2022]
Affiliation(s)
- Okba Al Rahal
- School of ChemistryCardiff University Park Place Cardiff Wales CF10 3AT UK
| | - Colan E. Hughes
- School of ChemistryCardiff University Park Place Cardiff Wales CF10 3AT UK
| | - P. Andrew Williams
- School of ChemistryCardiff University Park Place Cardiff Wales CF10 3AT UK
| | - Andrew J. Logsdail
- School of ChemistryCardiff University Park Place Cardiff Wales CF10 3AT UK
| | - Yael Diskin‐Posner
- Department of Chemical Research SupportWeizmann Institute of Science Rehovot 76100 Israel
| | | |
Collapse
|
43
|
Shishkina SV, Konovalova IS, Trostianko PV, Geleverya AO, Kovalenko SM, Bunyatyan ND. Polymorphism of 3-(5-phenyl-1,3,4-oxadiazol-2-yl)- and 3-[5-(pyridin-4-yl)-1,3,4-oxadiazol-2-yl]-2H-chromen-2-ones. ACTA CRYSTALLOGRAPHICA SECTION C-STRUCTURAL CHEMISTRY 2019; 75:1541-1553. [PMID: 31686666 DOI: 10.1107/s2053229619014256] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.2] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Subscribe] [Scholar Register] [Received: 09/03/2019] [Accepted: 10/18/2019] [Indexed: 02/08/2023]
Abstract
This study of 3-(5-phenyl-1,3,4-oxadiazol-2-yl)-2H-chromen-2-one, C17H10N2O3, 1, and 3-[5-(pyridin-4-yl)-1,3,4-oxadiazol-2-yl]-2H-chromen-2-one, C16H9N3O3, 2, was performed on the assumption of the potential anticancer activity of the compounds. Three polymorphic structures for 1 and two polymorphic structures for 2 have been studied thoroughly. The strongest intermolecular interaction is stacking of the `head-to-head' type in all the studied crystals. The polymorphic structures of 1 differ with respect to the intermolecular interactions between stacked columns. Two of the polymorphs have a columnar or double columnar type of crystal organization, while the third polymorphic structure can be classified as columnar-layered. The difference between the two structures of 2 is less pronounced. Both crystals can be considered as having very similar arrangements of neighbouring columns. The formation of polymorphic modifications is caused by a subtle balance of very weak intermolecular interactions and packing differences can be identified only using an analysis based on a study of the pairwise interaction energies.
Collapse
Affiliation(s)
- Svitlana V Shishkina
- SSI `Institute for Single Crystals' NAS of Ukraine, 60 Nauky ave., Kharkiv 61001, Ukraine
| | - Irina S Konovalova
- SSI `Institute for Single Crystals' NAS of Ukraine, 60 Nauky ave., Kharkiv 61001, Ukraine
| | - Pavlo V Trostianko
- V.N. Karazin Kharkiv National University, 4 Svobody sq., Kharkiv 61077, Ukraine
| | - Anna O Geleverya
- V.N. Karazin Kharkiv National University, 4 Svobody sq., Kharkiv 61077, Ukraine
| | - Sergiy M Kovalenko
- V.N. Karazin Kharkiv National University, 4 Svobody sq., Kharkiv 61077, Ukraine
| | - Natalya D Bunyatyan
- Federal State Autonomous Educational Institution of Higher Education, I.M. Sechenov First Moscow State Medical University, 8 Trubeckaya, Moscow 119991, Russian Federation
| |
Collapse
|
44
|
Bhardwaj RM, McMahon JA, Nyman J, Price LS, Konar S, Oswald IDH, Pulham CR, Price SL, Reutzel-Edens SM. A Prolific Solvate Former, Galunisertib, under the Pressure of Crystal Structure Prediction, Produces Ten Diverse Polymorphs. J Am Chem Soc 2019; 141:13887-13897. [DOI: 10.1021/jacs.9b06634] [Citation(s) in RCA: 78] [Impact Index Per Article: 15.6] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Affiliation(s)
- Rajni M. Bhardwaj
- Small Molecule Design & Development, Eli Lilly and Company, Indianapolis, Indiana 46285, United States
| | - Jennifer A. McMahon
- Small Molecule Design & Development, Eli Lilly and Company, Indianapolis, Indiana 46285, United States
| | - Jonas Nyman
- Small Molecule Design & Development, Eli Lilly and Company, Indianapolis, Indiana 46285, United States
- School of Pharmacy, University of Wisconsin—Madison, 777 Highland Avenue, Madison, Wisconsin 53705, United States
| | - Louise S. Price
- Department of Chemistry, University College London, 20 Gordon Street, London WC1H 0AJ, U.K
| | - Sumit Konar
- School of Chemistry, University of Edinburgh, David Brewster Road, Edinburgh EH9 3FJ, U.K
| | - Iain D. H. Oswald
- Strathclyde Institute of Pharmacy and Biomedical Sciences, University of Strathclyde, 161 Cathedral St, Glasgow G4 0RE, U.K
| | - Colin R. Pulham
- School of Chemistry, University of Edinburgh, David Brewster Road, Edinburgh EH9 3FJ, U.K
| | - Sarah L. Price
- Department of Chemistry, University College London, 20 Gordon Street, London WC1H 0AJ, U.K
| | - Susan M. Reutzel-Edens
- Small Molecule Design & Development, Eli Lilly and Company, Indianapolis, Indiana 46285, United States
| |
Collapse
|
45
|
Mortazavi M, Hoja J, Aerts L, Quéré L, van de Streek J, Neumann MA, Tkatchenko A. Computational polymorph screening reveals late-appearing and poorly-soluble form of rotigotine. Commun Chem 2019. [DOI: 10.1038/s42004-019-0171-y] [Citation(s) in RCA: 30] [Impact Index Per Article: 6.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/09/2022] Open
|
46
|
Crystal Structure Optimization and Gibbs Free Energy Comparison of Five Sulfathiazole Polymorphs by the Embedded Fragment QM Method at the DFT Level. CRYSTALS 2019. [DOI: 10.3390/cryst9050256] [Citation(s) in RCA: 6] [Impact Index Per Article: 1.2] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 11/16/2022]
Abstract
Molecular crystal plays an important role in many fields of science and technology, but it often crystallizes in different polymorphs with different physical properties. To guide the experimental synthesis of candidate materials, the atomic-scale model is frequently used to predict the most stable polymorph and its structural properties. Here, we show how an ab initio method can be used to achieve a rapid and accurate prediction of sulfathiazole crystal polymorphs (an antibiotic drug), based on the Gibbs free energy calculation and Raman spectra analysis. At the atmospheric pressure and the temperature of 300 K, we demonstrate that form III (FIII) is the most stable structure of sulfathiazole. The agreement between the predicted and experimental crystal structures corresponds to the order of stability for five sulfathiazole polymorphs as FI < FV < FIV < FII < FIII, which is achieved by employing the density functional theory (DFT) calculations.
Collapse
|
47
|
Veit M, Jain SK, Bonakala S, Rudra I, Hohl D, Csányi G. Equation of State of Fluid Methane from First Principles with Machine Learning Potentials. J Chem Theory Comput 2019; 15:2574-2586. [DOI: 10.1021/acs.jctc.8b01242] [Citation(s) in RCA: 31] [Impact Index Per Article: 6.2] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/06/2023]
Affiliation(s)
- Max Veit
- Engineering Laboratory, University of Cambridge, Trumpington Street, Cambridge CB2 1PZ, United Kingdom
| | | | | | - Indranil Rudra
- Shell India Markets
Pvt. Ltd., Bengaluru 562149, Karnataka, India
| | - Detlef Hohl
- Shell Global Solutions
International BV, Grasweg 31, 1031 HW Amsterdam, The Netherlands
| | - Gábor Csányi
- Engineering Laboratory, University of Cambridge, Trumpington Street, Cambridge CB2 1PZ, United Kingdom
| |
Collapse
|
48
|
Hoja J, Ko HY, Neumann MA, Car R, DiStasio RA, Tkatchenko A. Reliable and practical computational description of molecular crystal polymorphs. SCIENCE ADVANCES 2019; 5:eaau3338. [PMID: 30746448 PMCID: PMC6357866 DOI: 10.1126/sciadv.aau3338] [Citation(s) in RCA: 91] [Impact Index Per Article: 18.2] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Received: 05/30/2018] [Accepted: 11/28/2018] [Indexed: 05/12/2023]
Abstract
Reliable prediction of the polymorphic energy landscape of a molecular crystal would yield profound insight into drug development in terms of the existence and likelihood of late-appearing polymorphs. However, the computational prediction of molecular crystal polymorphs is highly challenging due to the high dimensionality of conformational and crystallographic space accompanied by the need for relative free energies to within 1 kJ/mol per molecule. In this study, we combine the most successful crystal structure sampling strategy with the most successful first-principles energy ranking strategy of the latest blind test of organic crystal structure prediction methods. Specifically, we present a hierarchical energy ranking approach intended for the refinement of relative stabilities in the final stage of a crystal structure prediction procedure. Such a combined approach provides excellent stability rankings for all studied systems and can be applied to molecular crystals of pharmaceutical importance.
Collapse
Affiliation(s)
- Johannes Hoja
- Physics and Materials Science Research Unit, University of Luxembourg, L-1511 Luxembourg, Luxembourg
| | - Hsin-Yu Ko
- Department of Chemistry, Princeton University, Princeton, NJ 08544, USA
| | - Marcus A. Neumann
- Avant-garde Materials Simulation Deutschland GmbH, Alte Str. 2, 79249 Merzhausen, Germany
| | - Roberto Car
- Department of Chemistry, Princeton University, Princeton, NJ 08544, USA
| | - Robert A. DiStasio
- Department of Chemistry and Chemical Biology, Cornell University, Ithaca, NY 14853, USA
| | - Alexandre Tkatchenko
- Physics and Materials Science Research Unit, University of Luxembourg, L-1511 Luxembourg, Luxembourg
- Corresponding author.
| |
Collapse
|
49
|
Shishkina SV, Konovalova IS, Kovalenko SM, Trostianko PV, Geleverya AO, Bunyatyan ND. Hydrogen bonding vs. stacking interaction in the crystals of the simplest coumarin derivatives: a study from the energetic viewpoint. CrystEngComm 2019. [DOI: 10.1039/c9ce01344j] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.6] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022]
Abstract
Hydrogen bonded and stacked dimers have very close interaction energies in the crystals of the simplest coumarin derivatives according to the data of the crystal structure analysis based on the comparison of pairwise interaction energies.
Collapse
Affiliation(s)
- Svitlana V. Shishkina
- SSI “Institute for Single Crystals” National Academy of Science of Ukraine
- Kharkiv
- Ukraine
- V. N. Karazin Kharkiv National University
- Kharkiv 61077
| | - Irina S. Konovalova
- SSI “Institute for Single Crystals” National Academy of Science of Ukraine
- Kharkiv
- Ukraine
| | - Sergiy M. Kovalenko
- V. N. Karazin Kharkiv National University
- Kharkiv 61077
- Ukraine
- Federal State Autonomous Educational Institution of Higher Education I.M. Sechenov First Moscow State Medical University of the Ministry of Healthcare of the Russian Federation (Sechenovskiy University)
- Moscow
| | | | | | - Natalya D. Bunyatyan
- Federal State Autonomous Educational Institution of Higher Education I.M. Sechenov First Moscow State Medical University of the Ministry of Healthcare of the Russian Federation (Sechenovskiy University)
- Moscow
- Russia
- Federal State Budgetary Institution “Scientific Centre for Expert Evaluation of Medicinal Products” of the Ministry of Health of the Russian Federation
- Moscow
| |
Collapse
|
50
|
Stöhr M, Van Voorhis T, Tkatchenko A. Theory and practice of modeling van der Waals interactions in electronic-structure calculations. Chem Soc Rev 2019; 48:4118-4154. [PMID: 31190037 DOI: 10.1039/c9cs00060g] [Citation(s) in RCA: 68] [Impact Index Per Article: 13.6] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/21/2022]
Abstract
The accurate description of long-range electron correlation, most prominently including van der Waals (vdW) dispersion interactions, represents a particularly challenging task in the modeling of molecules and materials. vdW forces arise from the interaction of quantum-mechanical fluctuations in the electronic charge density. Within (semi-)local density functional approximations or Hartree-Fock theory such interactions are neglected altogether. Non-covalent vdW interactions, however, are ubiquitous in nature and play a key role for the understanding and accurate description of the stability, dynamics, structure, and response properties in a plethora of systems. During the last decade, many promising methods have been developed for modeling vdW interactions in electronic-structure calculations. These methods include vdW-inclusive Density Functional Theory and correlated post-Hartree-Fock approaches. Here, we focus on the methods within the framework of Density Functional Theory, including non-local van der Waals density functionals, interatomic dispersion models within many-body and pairwise formulation, and random phase approximation-based approaches. This review aims to guide the reader through the theoretical foundations of these methods in a tutorial-style manner and, in particular, highlight practical aspects such as the applicability and the advantages and shortcomings of current vdW-inclusive approaches. In addition, we give an overview of complementary experimental approaches, and discuss tools for the qualitative understanding of non-covalent interactions as well as energy decomposition techniques. Besides representing a reference for the current state-of-the-art, this work is thus also designed as a concise and detailed introduction to vdW-inclusive electronic structure calculations for a general and broad audience.
Collapse
Affiliation(s)
- Martin Stöhr
- Physics and Materials Science Research Unit, University of Luxembourg, L-1511 Luxembourg, Luxembourg.
| | | | | |
Collapse
|