1
|
Li HQ, Xu JY, Wu SS, Jin L. Molecular Insights into Water-Chloride and Water-Water Interactions in the Supramolecular Architecture of Aconine Hydrochloride Dihydrate. ACS OMEGA 2024; 9:11925-11941. [PMID: 38496984 PMCID: PMC10938397 DOI: 10.1021/acsomega.3c09696] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Download PDF] [Figures] [Subscribe] [Scholar Register] [Received: 12/04/2023] [Revised: 01/23/2024] [Accepted: 02/14/2024] [Indexed: 03/19/2024]
Abstract
Despite the previous preparation of aconine hydrochloride monohydrate (AHM), accurate determination of the crystal's composition was hindered by severely disordered water molecules within the crystal. In this study, we successfully prepared a new dihydrate form of the aconine hydrochloride [C25H42NO9+Cl-·2(H2O), aconine hydrochloride dihydrate (AHD)] and accurately refined all water molecules within the AHD crystal. Our objective is to elucidate both water-chloride and water-water interactions in the AHD crystal. The crystal structure of AHD was determined at 136 K using X-ray diffraction and a multipolar atom model was constructed by transferring charge-density parameters to explore the topological features of key short contacts. By comparing the crystal structures of dihydrate and monohydrate forms, we have observed that both AHD and AHM exhibit identical aconine cations, except for variations in the number of water molecules present. In the AHD crystal, chloride anions and water molecules serve as pivotal connecting hubs to establish three-dimensional hydrogen bonding networks and one-dimensional hydrogen bonding chain; both water-chloride and water-water interactions assemble supramolecular architectures. The crystal packing of AHD exhibits a complete reversal in the stacking order compared to AHM, thereby emphasizing distinct disparities between them. Hirshfeld surface analysis reveals that H···Cl- and H···O contacts play a significant role in constructing the hydrogen bonding network and chain within these supramolecular architectures. Furthermore, topological analysis and electrostatic interaction energy confirm that both water-chloride and water-water interactions stabilize supramolecular architectures through electrostatic attraction facilitated by H···Cl- and H···O contacts. Importantly, these findings are strongly supported by the existing literature evidence. Consequently, navigating these water-chloride and water-water interactions is imperative for ensuring storage and safe processing of this pharmaceutical compound.
Collapse
Affiliation(s)
- Han-Qing Li
- State
Clinical Trial Institution of New Drugs, International Mongolian Hospital of Inner Mongolia, Hohhot, Inner Mongolia 010065, People’s Republic
of China
- Mongolian
Medicine Laboratory, International Mongolian
Hospital of Inner Mongolia, Hohhot, Inner Mongolia 010065, People’s Republic of China
| | - Jia-Yin Xu
- Mongolian
Pharmaceutical Preparation Center, International
Mongolian Hospital of Inner Mongolia, Hohhot, Inner Mongolia 010065, People’s Republic
of China
| | - Shan-Shan Wu
- State
Clinical Trial Institution of New Drugs, International Mongolian Hospital of Inner Mongolia, Hohhot, Inner Mongolia 010065, People’s Republic
of China
| | - Liang Jin
- State
Clinical Trial Institution of New Drugs, International Mongolian Hospital of Inner Mongolia, Hohhot, Inner Mongolia 010065, People’s Republic
of China
| |
Collapse
|
2
|
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
|
3
|
Klitou P, Parisi E, Bordignon S, Bravetti F, Rosbottom I, Dell’Aera M, Cuocci C, Chierotti MR, Altomare A, Simone E. Navigating the Complex Solid Form Landscape of the Quercetin Flavonoid Molecule. CRYSTAL GROWTH & DESIGN 2023; 23:6034-6045. [PMID: 37547879 PMCID: PMC10401642 DOI: 10.1021/acs.cgd.3c00584] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Figures] [Subscribe] [Scholar Register] [Received: 05/14/2023] [Revised: 07/01/2023] [Indexed: 08/08/2023]
Abstract
Quercetin, a naturally occurring bioflavonoid substance widely used in the nutraceutical and food industries, exists in various solid forms that can have different physicochemical properties, thus impacting this compound's performance in various applications. In this work, we will clarify the complex solid-form landscape of this molecule. Two elusive isostructural solvates of quercetin were obtained from ethanol and methanol. The obtained crystals were characterized experimentally, but the crystallographic structure could not be solved due to their high instability. Nevertheless, the desolvated structure resulting from a high-temperature treatment (or prolonged storage at ambient conditions) of both these two labile crystals was characterized and solved via powder X-ray diffraction and solid-state nuclear magnetic resonance (SSNMR). This anhydrous crystal structure was compared with another anhydrous quercetin form obtained in our previous work, indicating that, at least, two different anhydrous polymorphs of quercetin exist. Navigating the solid-form landscape of quercetin is essential to ensure accurate control of the functional properties of food, nutraceutical, or pharmaceutical products containing crystal forms of this substance.
Collapse
Affiliation(s)
- Panayiotis Klitou
- School
of Food Science and Nutrition, Food Colloids and Bioprocessing Group, University of Leeds, Leeds LS2 9JT, UK
| | - Emmanuele Parisi
- Department
of Applied Science and Technology (DISAT), Politecnico di Torino, Torino I-10129, Italy
| | - Simone Bordignon
- Dipartimento
di Chimica I.F.M, Università degli
Studi di Torino, Via P. Giuria 7, Torino I-10125, Italy
| | - Federica Bravetti
- Dipartimento
di Chimica I.F.M, Università degli
Studi di Torino, Via P. Giuria 7, Torino I-10125, Italy
| | - Ian Rosbottom
- School
of Chemical and Process Engineering, University
of Leeds, Woodhouse
Lane Leeds LS2 9JT, UK
| | - Marzia Dell’Aera
- Institute
of Crystallography IC − CNR, via Amendola 122/O, Bari I-70126, Italy
| | - Corrado Cuocci
- Institute
of Crystallography IC − CNR, via Amendola 122/O, Bari I-70126, Italy
| | - Michele R. Chierotti
- Dipartimento
di Chimica I.F.M, Università degli
Studi di Torino, Via P. Giuria 7, Torino I-10125, Italy
| | - Angela Altomare
- Institute
of Crystallography IC − CNR, via Amendola 122/O, Bari I-70126, Italy
| | - Elena Simone
- School
of Food Science and Nutrition, Food Colloids and Bioprocessing Group, University of Leeds, Leeds LS2 9JT, UK
- Department
of Applied Science and Technology (DISAT), Politecnico di Torino, Torino I-10129, Italy
| |
Collapse
|
4
|
Cook CJ, Li W, Lui BF, Gately TJ, Al-Kaysi RO, Mueller LJ, Bardeen CJ, Beran GJO. A theoretical framework for the design of molecular crystal engines. Chem Sci 2023; 14:937-949. [PMID: 36755715 PMCID: PMC9890974 DOI: 10.1039/d2sc05549j] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/06/2022] [Accepted: 12/19/2022] [Indexed: 12/24/2022] Open
Abstract
Photomechanical molecular crystals have garnered attention for their ability to transform light into mechanical work, but difficulties in characterizing the structural changes and mechanical responses experimentally have hindered the development of practical organic crystal engines. This study proposes a new computational framework for predicting the solid-state crystal-to-crystal photochemical transformations entirely from first principles, and it establishes a photomechanical engine cycle that quantifies the anisotropic mechanical performance resulting from the transformation. The approach relies on crystal structure prediction, solid-state topochemical principles, and high-quality electronic structure methods. After validating the framework on the well-studied [4 + 4] cycloadditions in 9-methyl anthracene and 9-tert-butyl anthracene ester, the experimentally-unknown solid-state transformation of 9-carboxylic acid anthracene is predicted for the first time. The results illustrate how the mechanical work is done by relaxation of the crystal lattice to accommodate the photoproduct, rather than by the photochemistry itself. The large ∼107 J m-3 work densities computed for all three systems highlight the promise of photomechanical crystal engines. This study demonstrates the importance of crystal packing in determining molecular crystal engine performance and provides tools and insights to design improved materials in silico.
Collapse
Affiliation(s)
- Cameron J. Cook
- Department of Chemistry, University of California RiversideRiverside CA 92521USA
| | - Wangxiang Li
- Department of Chemistry, University of California Riverside Riverside CA 92521 USA
| | - Brandon F. Lui
- Department of Chemistry, University of California RiversideRiverside CA 92521USA
| | - Thomas J. Gately
- Department of Chemistry, University of California RiversideRiverside CA 92521USA
| | - Rabih O. Al-Kaysi
- College of Science and Health Professions-3124, King Saud Bin Abdulaziz University for Health Sciences, and King Abdullah International Medical Research Center, Ministry of National Guard Health AffairsRiyadh 11426Kingdom of Saudi Arabia
| | - Leonard J. Mueller
- Department of Chemistry, University of California RiversideRiverside CA 92521USA
| | | | - Gregory J. O. Beran
- Department of Chemistry, University of California RiversideRiverside CA 92521USA
| |
Collapse
|
5
|
Mattei A, Hong RS, Dietrich H, Firaha D, Helfferich J, Liu YM, Sasikumar K, Abraham NS, Miglani Bhardwaj R, Neumann MA, Sheikh AY. Efficient Crystal Structure Prediction for Structurally Related Molecules with Accurate and Transferable Tailor-Made Force Fields. J Chem Theory Comput 2022; 18:5725-5738. [PMID: 35930763 PMCID: PMC9476662 DOI: 10.1021/acs.jctc.2c00451] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/06/2023]
Abstract
Crystal structure prediction (CSP) his generally used to complement experimental solid form screening and applied to individual molecules in drug development. The fast development of algorithms and computing resources offers the opportunity to use CSP earlier and for a broader range of applications in the drug design cycle. This study presents a novel paradigm of CSP specifically designed for structurally related molecules, referred to as Quick-CSP. The approach prioritizes more accurate physics through robust and transferable tailor-made force fields (TMFFs), such that significant efficiency gains are achieved through the reduction of expensive ab initio calculations. The accuracy of the TMFF is increased by the introduction of electrostatic multipoles, and the fragment-based force field parameterization scheme is demonstrated to be transferable for a family of chemically related molecules. The protocol is benchmarked with structurally related compounds from the Bromodomain and Extraterminal (BET) domain inhibitors series. A new convergence criterion is introduced that aims at performing only as many ab initio optimizations of crystal structures as required to locate the bottom of the crystal energy landscape within a user-defined accuracy. The overall approach provides significant cost savings ranging from three- to eight-fold less than the full-CSP workflow. The reported advancements expand the scope and utility of the underlying CSP building blocks as well as their novel reassembly to other applications earlier in the drug design cycle to guide molecule design and selection.
Collapse
Affiliation(s)
- Alessandra Mattei
- Solid State Chemistry, Research & Development, AbbVie Inc., 1 N Waukegan Road, North Chicago, Illinois 60064, United States
| | - Richard S Hong
- Solid State Chemistry, Research & Development, AbbVie Inc., 1 N Waukegan Road, North Chicago, Illinois 60064, United States
| | - Hanno Dietrich
- Avant-garde Materials Simulation, GmbH, Alte Str. 2, 79249 Merzhausen, Germany
| | - Dzmitry Firaha
- Avant-garde Materials Simulation, GmbH, Alte Str. 2, 79249 Merzhausen, Germany
| | - Julian Helfferich
- Avant-garde Materials Simulation, GmbH, Alte Str. 2, 79249 Merzhausen, Germany
| | - Yifei Michelle Liu
- Avant-garde Materials Simulation, GmbH, Alte Str. 2, 79249 Merzhausen, Germany
| | - Kiran Sasikumar
- Avant-garde Materials Simulation, GmbH, Alte Str. 2, 79249 Merzhausen, Germany
| | - Nathan S Abraham
- Solid State Chemistry, Research & Development, AbbVie Inc., 1 N Waukegan Road, North Chicago, Illinois 60064, United States
| | - Rajni Miglani Bhardwaj
- Solid State Chemistry, Research & Development, AbbVie Inc., 1 N Waukegan Road, North Chicago, Illinois 60064, United States
| | - Marcus A Neumann
- Avant-garde Materials Simulation, GmbH, Alte Str. 2, 79249 Merzhausen, Germany
| | - Ahmad Y Sheikh
- Solid State Chemistry, Research & Development, AbbVie Inc., 1 N Waukegan Road, North Chicago, Illinois 60064, United States
| |
Collapse
|
6
|
Sanii R, Andaloussi YH, Patyk-Kaźmierczak E, Zaworotko MJ. Polymorphism in Ionic Cocrystals Comprising Lithium Salts and l-Proline. CRYSTAL GROWTH & DESIGN 2022; 22:3786-3794. [PMID: 36160301 PMCID: PMC9490868 DOI: 10.1021/acs.cgd.2c00172] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Figures] [Subscribe] [Scholar Register] [Received: 02/09/2022] [Revised: 04/20/2022] [Indexed: 06/16/2023]
Abstract
The occurrence of polymorphism in ionic cocrystals formed by two lithium salts, lithium salicylate (LIS) and lithium 4-methoxybenzoate (L4M), and l-proline (PRO) has been investigated. The previously reported monoclinic form of the 1:1 cocrystal of LIS and PRO, LISPRO(α), and a new thermodynamically stable orthorhombic polymorph, LISPRO(β), were prepared and characterized. The two polymorphs form square grid, sql, topology coordination networks and differ mainly in the conformation of the salicylate ions and positioning of the sql nets. LISPRO(α) was observed to transform to LISPRO(β) under slurry conditions. The 1:1 ionic cocrystal of L4M and PRO (L4MPRO) was found to form three polymorphs. Apart from the previously reported orthorhombic crystal form, L4MPRO(α), two new monoclinic crystal forms, L4MPRO(β) and L4MPRO(γ), were obtained by modifying crystallization conditions. The new polymorphs were found to be metastable, undergoing transformations to L4MPRO(α) upon exposure to humidity. Experimental conditions that induce transformations between the polymorphs of LISPRO and L4MPRO are detailed, and the structural differences between the polymorphs are discussed in the broader context of polymorphism.
Collapse
Affiliation(s)
- Rana Sanii
- Department
of Chemical Sciences and Bernal Institute, University of Limerick, Co., Limerick V94T9PX, Ireland
| | - Yassin H. Andaloussi
- Department
of Chemical Sciences and Bernal Institute, University of Limerick, Co., Limerick V94T9PX, Ireland
| | - Ewa Patyk-Kaźmierczak
- Department
of Materials Chemistry, Faculty of Chemistry, Adam Mickiewicz University in Poznań, Uniwersytetu Poznańskiego 8, Poznań 61-614, Poland
| | - Michael J. Zaworotko
- Department
of Chemical Sciences and Bernal Institute, University of Limerick, Co., Limerick V94T9PX, Ireland
| |
Collapse
|
7
|
|
8
|
Abramov YA, Sun G, Zeng Q. Emerging Landscape of Computational Modeling in Pharmaceutical Development. J Chem Inf Model 2022; 62:1160-1171. [PMID: 35226809 DOI: 10.1021/acs.jcim.1c01580] [Citation(s) in RCA: 11] [Impact Index Per Article: 5.5] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/06/2023]
Abstract
Computational chemistry applications have become an integral part of the drug discovery workflow over the past 35 years. However, computational modeling in support of drug development has remained a relatively uncharted territory for a significant part of both academic and industrial communities. This review considers the computational modeling workflows for three key components of drug preclinical and clinical development, namely, process chemistry, analytical research and development, as well as drug product and formulation development. An overview of the computational support for each step of the respective workflows is presented. Additionally, in context of solid form design, special consideration is given to modern physics-based virtual screening methods. This covers rational approaches to polymorph, coformer, counterion, and solvent virtual screening in support of solid form selection and design.
Collapse
Affiliation(s)
- Yuriy A Abramov
- XtalPi, Inc., 245 Main St., Cambridge, Massachusetts 02142, United States.,Eshelman School of Pharmacy, University of North Carolina, Chapel Hill, North Carolina 27599, United States
| | - Guangxu Sun
- XtalPi, Inc., Shenzhen Jingtai Technology Co., Ltd., Floor 3, Sf Industrial Plant, No. 2 Hongliu road, Fubao Community, Fubao Street, Futian District, Shenzhen 518100, China
| | - Qun Zeng
- XtalPi, Inc., Shenzhen Jingtai Technology Co., Ltd., Floor 3, Sf Industrial Plant, No. 2 Hongliu road, Fubao Community, Fubao Street, Futian District, Shenzhen 518100, China
| |
Collapse
|
9
|
Beran GJO, Sugden IJ, Greenwell C, Bowskill DH, Pantelides CC, Adjiman CS. How many more polymorphs of ROY remain undiscovered. Chem Sci 2022; 13:1288-1297. [PMID: 35222912 PMCID: PMC8809489 DOI: 10.1039/d1sc06074k] [Citation(s) in RCA: 13] [Impact Index Per Article: 6.5] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/02/2021] [Accepted: 12/10/2021] [Indexed: 12/15/2022] Open
Abstract
With 12 crystal forms, 5-methyl-2-[(2-nitrophenyl)amino]-3-thiophenecabonitrile (a.k.a. ROY) holds the current record for the largest number of fully characterized organic crystal polymorphs. Four of these polymorph structures have been reported since 2019, raising the question of how many more ROY polymorphs await future discovery. Employing crystal structure prediction and accurate energy rankings derived from conformational energy-corrected density functional theory, this study presents the first crystal energy landscape for ROY that agrees well with experiment. The lattice energies suggest that the seven most stable ROY polymorphs (and nine of the twelve lowest-energy forms) on the Z' = 1 landscape have already been discovered experimentally. Discovering any new polymorphs at ambient pressure will likely require specialized crystallization techniques capable of trapping metastable forms. At pressures above 10 GPa, however, a new crystal form is predicted to become enthalpically more stable than all known polymorphs, suggesting that further high-pressure experiments on ROY may be warranted. This work highlights the value of high-accuracy crystal structure prediction for solid-form screening and demonstrates how pragmatic conformational energy corrections can overcome the limitations of conventional density functionals for conformational polymorphs.
Collapse
Affiliation(s)
- Gregory J O Beran
- Department of Chemistry, University of California Riverside Riverside CA 92521 USA
| | - Isaac J Sugden
- Department of Chemical Engineering, Sargent Centre for Process Systems Engineering, Imperial College London London SW7 2AZ UK
| | - Chandler Greenwell
- Department of Chemistry, University of California Riverside Riverside CA 92521 USA
| | - David H Bowskill
- Department of Chemical Engineering, Sargent Centre for Process Systems Engineering, Imperial College London London SW7 2AZ UK
| | - Constantinos C Pantelides
- Department of Chemical Engineering, Sargent Centre for Process Systems Engineering, Imperial College London London SW7 2AZ UK
| | - Claire S Adjiman
- Department of Chemical Engineering, Sargent Centre for Process Systems Engineering, Imperial College London London SW7 2AZ UK
| |
Collapse
|
10
|
Cook C, McKinley JL, Beran GJO. Modeling the α- and β-resorcinol phase boundary via combination of density functional theory and density functional tight-binding. J Chem Phys 2021; 154:134109. [PMID: 33832233 PMCID: PMC8019358 DOI: 10.1063/5.0044385] [Citation(s) in RCA: 3] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/15/2021] [Accepted: 03/10/2021] [Indexed: 02/06/2023] Open
Abstract
The ability to predict not only what organic crystal structures might occur but also the thermodynamic conditions under which they are the most stable would be extremely useful for discovering and designing new organic materials. The present study takes a step in that direction by predicting the temperature- and pressure-dependent phase boundary between the α and β polymorphs of resorcinol using density functional theory (DFT) and the quasi-harmonic approximation. To circumvent the major computational bottleneck associated with computing a well-converged phonon density of states via the supercell approach, a recently developed approximation is employed, which combines a supercell phonon density of states from dispersion-corrected third-order density functional tight binding [DFTB3-D3(BJ)] with frequency corrections derived from a smaller B86bPBE-XDM functional DFT phonon calculation on the crystallographic unit cell. This mixed DFT/DFTB quasi-harmonic approach predicts the lattice constants and unit cell volumes to within 1%-2% at lower pressures. It predicts the thermodynamic phase boundary in almost perfect agreement with the experiment, although this excellent agreement does reflect fortuitous cancellation of errors between the enthalpy and entropy of transition.
Collapse
Affiliation(s)
- Cameron Cook
- Department of Chemistry, University of California, Riverside, California 92521, USA
| | - Jessica L. McKinley
- Department of Chemistry, University of California, Riverside, California 92521, USA
| | - Gregory J. O. Beran
- Department of Chemistry, University of California, Riverside, California 92521, USA
| |
Collapse
|
11
|
Braun DE, Gelbrich T, Kahlenberg V, Griesser UJ. The Eight Hydrates of Strychnine Sulfate. CRYSTAL GROWTH & DESIGN 2020; 20:6069-6083. [PMID: 32913424 PMCID: PMC7472434 DOI: 10.1021/acs.cgd.0c00777] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Figures] [Subscribe] [Scholar Register] [Received: 06/06/2020] [Revised: 08/05/2020] [Indexed: 06/11/2023]
Abstract
Commercial samples of strychnine sulfate were used as the starting material in crystallization experiments accompanied by stability studies. Eight hydrate forms (HyA-HyG), including five novel hydrates, were verified. The crystal structures of HyA ("pentahydrate") and HyF ("hexahydrate") were determined from single-crystal X-ray diffraction data. HyF was identified as the most stable hydrate at high water activities at room temperature (RT), and HyA and HyC were also found to be stable at ambient conditions. Long-time storage experiments over nearly two decades confirm that these three hydrates are stable at ambient conditions (20-60% relative humidity). The other five hydrates, HyB ("dihydrate"), HyD, HyE, HyG, and HyH, are only observable at the low(est) relative humidity (RH) levels at RT. Some of these latter forms can only exist within a very narrow RH range and are therefore intermediate phases. By applying a range of complementary experimental techniques such as gravimetric moisture sorption analysis, thermal analysis, moisture controlled PXRD measurements, and variable temperature IR spectroscopy in combination with principal component analysis, it was possible to identify the distinct hydrate phases and elucidate their stability and dehydration pathways. The observed (de)hydration routes, HyA ↔ HyB, HyC ↔ HyD ↔ HyE, HyF ↔ HyG ↔ HyH and HyF → HyA ↔ HyB, depended on the initial hydrate form, particle size, and atmospheric conditions. In addition, a transformation from HyC/HyA to HyF occurs at high RH values at RT. The specific moisture and temperature conditions of none of the applied drying regimes yielded a crystalline water-free form, which highlights the essential role of water molecules for the formation and stability of the crystalline strychnine sulfate phases.
Collapse
Affiliation(s)
- Doris E. Braun
- Institute
of Pharmacy, University of Innsbruck, Innrain 52c, 6020 Innsbruck, Austria
| | - Thomas Gelbrich
- Institute
of Pharmacy, University of Innsbruck, Innrain 52c, 6020 Innsbruck, Austria
| | - Volker Kahlenberg
- Institute
of Mineralogy and Petrography, University
of Innsbruck, Innrain 52, 6020 Innsbruck, Austria
| | - Ulrich J. Griesser
- Institute
of Pharmacy, University of Innsbruck, Innrain 52c, 6020 Innsbruck, Austria
| |
Collapse
|
12
|
Qiu J, Patel A, Stevens JM. High-Throughput Salt Screening of Synthetic Intermediates: Effects of Solvents, Counterions, and Counterion Solubility. Org Process Res Dev 2020. [DOI: 10.1021/acs.oprd.0c00132] [Citation(s) in RCA: 5] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/13/2022]
Affiliation(s)
- Jun Qiu
- Product Development, Bristol-Myers Squibb Company, One Squibb Drive, New Brunswick, New Jersey 08903, United States
| | - Anisha Patel
- Product Development, Bristol-Myers Squibb Company, One Squibb Drive, New Brunswick, New Jersey 08903, United States
| | - Jason M. Stevens
- Product Development, Bristol-Myers Squibb Company, One Squibb Drive, New Brunswick, New Jersey 08903, United States
| |
Collapse
|
13
|
Kodama K, Yi M, Shitara H, Hirose T. Chirality switching in the enantioseparation of 2-hydroxy-4-phenylbutyric acid: Role of solvents in selective crystallization of the diastereomeric salt. Tetrahedron Lett 2020. [DOI: 10.1016/j.tetlet.2020.151773] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/24/2022]
|
14
|
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
|
15
|
Sacchi P, Lusi M, Cruz-Cabeza AJ, Nauha E, Bernstein J. Same or different – that is the question: identification of crystal forms from crystal structure data. CrystEngComm 2020. [DOI: 10.1039/d0ce00724b] [Citation(s) in RCA: 16] [Impact Index Per Article: 4.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/22/2022]
Abstract
An analysis of the CSD with structural comparison tools shows that differentiating between polymorphism and redeterminations is not always straight forward and requires of complementary tools at the hands of an expert practitioner.
Collapse
Affiliation(s)
- Pietro Sacchi
- Department of Chemical Engineering and Analytical Science
- University of Manchester
- Manchester
- UK
| | - Matteo Lusi
- SSPC
- The SFI Research Centre for Pharmaceuticals
- Bernal Institute
- Department of Chemical Sciences
- University of Limerick
| | - Aurora J. Cruz-Cabeza
- Department of Chemical Engineering and Analytical Science
- University of Manchester
- Manchester
- UK
| | - Elisa Nauha
- New York University Abu Dhabi
- Abu Dhabi
- United Arab Emirates
| | - Joel Bernstein
- New York University Abu Dhabi
- Abu Dhabi
- United Arab Emirates
| |
Collapse
|
16
|
Braun DE. Experimental and computational approaches to rationalise multicomponent supramolecular assemblies: dapsone monosolvates. Phys Chem Chem Phys 2019; 21:17288-17305. [PMID: 31348477 DOI: 10.1039/c9cp02572c] [Citation(s) in RCA: 5] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022]
Abstract
The monosolvate crystal energy landscapes of dapsone (DDS) including the solvents carbon tetrachloride, acetone, cyclohexanone, dimethyl formamide, tetrahydrofuran, methyl ethyl ketone, 1,2-dichloroethane, 1,4-dioxane, dichloromethane and chloroform were established using experimental and computational approaches. To rationalise and understand solvate formation, solvate stability and desolvation reactions a careful control of the experimental crystallisation and storage conditions, a range of thermoanalytical methods and crystal structure prediction were required. Six of the eight DDS monosolvates are reported and characterised for the first time. Structural similarity and diversity of the at ambient conditions unstable monosolvates were apparent from the computed crystal energy landscapes, which had the experimental packings as lowest energy structures. The computed structures were used as input for Rietveld refinements and isostructurality of four of the monosolvates was confirmed. Packing comparisons of the solvate structures and molecular properties of the solvent molecules indicated that both size/shape of the solvent molecule and the possible DDSsolvent interactions are the important factors for DDS solvate formation. Through the combination of experiment and theory solvate stability and structural features have been rationalised.
Collapse
Affiliation(s)
- Doris E Braun
- Institute of Pharmacy, University of Innsbruck, Innrain 52c, 6020 Innsbruck, Austria.
| |
Collapse
|
17
|
Dudek MK, Day GM. Explaining crystallization preferences of two polyphenolic diastereoisomers by crystal structure prediction. CrystEngComm 2019. [DOI: 10.1039/c8ce01783b] [Citation(s) in RCA: 17] [Impact Index Per Article: 3.4] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/07/2023]
Abstract
Crystal structure prediction is used to understand the differences in crystallization of catechin and epicatechin, and to explore the predictability of solvate formation.
Collapse
Affiliation(s)
- Marta K. Dudek
- Computational Systems Chemistry
- School of Chemistry
- University of Southampton
- UK
- Center of Molecular and Macromolecular Studies PAS
| | - Graeme M. Day
- Computational Systems Chemistry
- School of Chemistry
- University of Southampton
- UK
| |
Collapse
|
18
|
Stevenson EL, Lancaster RW, Buanz ABM, Price LS, Tocher DA, Price SL. The solid state forms of the sex hormone 17-β-estradiol. CrystEngComm 2019. [DOI: 10.1039/c8ce01874j] [Citation(s) in RCA: 10] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/06/2023]
Abstract
The crystal structure of the female sex hormone has been established despite its high affinity for water.
Collapse
Affiliation(s)
| | | | | | - Louise S. Price
- Department of Chemistry
- University College London
- London WC1H 0AJ
- UK
| | - Derek A. Tocher
- Department of Chemistry
- University College London
- London WC1H 0AJ
- UK
| | - Sarah L. Price
- Department of Chemistry
- University College London
- London WC1H 0AJ
- UK
| |
Collapse
|
19
|
Playford HY, Whale TF, Murray BJ, Tucker MG, Salzmann CG. Analysis of stacking disorder in ice I using pair distribution functions. J Appl Crystallogr 2018. [DOI: 10.1107/s1600576718009056] [Citation(s) in RCA: 14] [Impact Index Per Article: 2.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/10/2022] Open
Abstract
Stacking-disordered materials display crystalline order in two dimensions but are disordered along the direction in which layered structural motifs are stacked. Countless examples of stacking disorder exist, ranging from close-packed metals, ice I and diamond to open-framework materials and small-molecule pharmaceuticals. In general, the presence of stacking disorder can have profound consequences for the physical and chemical properties of a material. Traditional analyses of powder diffraction data are often complicated by the presence of memory effects in the stacking sequences. Here it is shown that experimental pair distribution functions of stacking-disordered ice I can be used to determine local information on the fractions of cubic and hexagonal stacking. Ice is a particularly challenging material in this respect, since both the stacking disorder and the orientational disorder of the water molecules need to be described. Memory effects are found to contribute very little to the pair distribution functions, and consequently, the analysis of pair distribution functions is the method of choice for characterizing stacking-disordered samples with complicated and high-order memory effects. In the context of this work, the limitations of current structure-reconstruction approaches are also discussed.
Collapse
|
20
|
Iuzzolino L, McCabe P, Price SL, Brandenburg JG. Crystal structure prediction of flexible pharmaceutical-like molecules: density functional tight-binding as an intermediate optimisation method and for free energy estimation. Faraday Discuss 2018; 211:275-296. [PMID: 30035288 DOI: 10.1039/c8fd00010g] [Citation(s) in RCA: 23] [Impact Index Per Article: 3.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022]
Abstract
Successful methodologies for theoretical crystal structure prediction (CSP) on flexible pharmaceutical-like organic molecules explore the lattice energy surface to find a set of plausible crystal structures. The initial search stages of CSP studies use relatively simple lattice energy approximations as hundreds of thousands of minima have to be considered. These generated crystal structures often have poor molecular geometries, as well as inaccurate lattice energy rankings, and performing reasonably accurate but computationally affordable optimisations of the crystal structures generated in a search would be highly desirable. Here, we seek to explore whether semi-empirical quantum-mechanical methods can perform this task. We employed the dispersion-corrected tight-binding Hamiltonian (DFTB3-D3) to relax all the inter- and intra-molecular degrees of freedom of several thousands of generated crystal structures of five pharmaceutical-like molecules, saving a large amount of computational effort compared to earlier studies. The computational cost scales better with molecular size and flexibility than other CSP methods, suggesting that it could be extended to even larger and more flexible molecules. On average, this optimisation improved the average reproduction of the eight experimental crystal structures (RMSD15) and experimental conformers (RMSD1) by 4% and 23%, respectively. The intermolecular interactions were then further optimised using distributed multipoles, derived from the molecular wave-functions, to accurately describe the electrostatic components of the intermolecular energies. In all cases, the experimental crystal structures are close to the top of the lattice energy ranking. Phonon calculations on some of the lowest energy structures were also performed with DFTB3-D3 methods to calculate the vibrational component of the Helmholtz free energy, providing further insights into the solid-state behaviour of the target molecules. We conclude that DFTB3-D3 is a cost-effective method for optimising flexible molecules, bridging the gap between the approximate methods used in CSP searches for generating crystal structures and more accurate methods required in the final energy ranking.
Collapse
Affiliation(s)
- Luca Iuzzolino
- Department of Chemistry, University College London, 20 Gordon Street, London WC1H 0AJ, UK.
| | | | | | | |
Collapse
|
21
|
Price SL. Is zeroth order crystal structure prediction (CSP_0) coming to maturity? What should we aim for in an ideal crystal structure prediction code? Faraday Discuss 2018; 211:9-30. [PMID: 30051901 DOI: 10.1039/c8fd00121a] [Citation(s) in RCA: 50] [Impact Index Per Article: 8.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022]
Abstract
Crystal structure prediction based on searching for the global minimum in the lattice energy (CSP_0) is growing in use for guiding the discovery of new materials, for example, new functional materials, new phases of interest to planetary scientists and new polymorphs relevant to pharmaceutical development. This Faraday Discussion can assess the progress of CSP_0 over the range of types of materials to which CSP is currently and could be applied, which depends on our ability to model the variety of interatomic forces in crystals. The basic hypothesis, that the outcome of crystallisation is determined by thermodynamics, needs examining by considering methods of modelling relative thermodynamic stability not only as a function of pressure and temperature, but also of size, solvent and the presence of heterogeneous templates or impurities (CSP_thd). Given that many important materials persist, and indeed may be formed, when they are not the most thermodynamically stable structure, we need to define what would be required of an ideal CSP code (CSP_aim).
Collapse
Affiliation(s)
- Sarah L Price
- Department of Chemistry, University College London, 20 Gordon Street, London WC1H 0AJ, UK.
| |
Collapse
|