1
|
Rahman M, Dannatt HRW, Blundell CD, Hughes LP, Blade H, Carson J, Tatman BP, Johnston ST, Brown SP. Polymorph Identification for Flexible Molecules: Linear Regression Analysis of Experimental and Calculated Solution- and Solid-State NMR Data. J Phys Chem A 2024; 128:1793-1816. [PMID: 38427685 PMCID: PMC10945485 DOI: 10.1021/acs.jpca.3c07732] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/24/2023] [Revised: 02/06/2024] [Accepted: 02/07/2024] [Indexed: 03/03/2024]
Abstract
The Δδ regression approach of Blade et al. [ J. Phys. Chem. A 2020, 124(43), 8959-8977] for accurately discriminating between solid forms using a combination of experimental solution- and solid-state NMR data with density functional theory (DFT) calculation is here extended to molecules with multiple conformational degrees of freedom, using furosemide polymorphs as an exemplar. As before, the differences in measured 1H and 13C chemical shifts between solution-state NMR and solid-state magic-angle spinning (MAS) NMR (Δδexperimental) are compared to those determined by gauge-including projector augmented wave (GIPAW) calculations (Δδcalculated) by regression analysis and a t-test, allowing the correct furosemide polymorph to be precisely identified. Monte Carlo random sampling is used to calculate solution-state NMR chemical shifts, reducing computation times by avoiding the need to systematically sample the multidimensional conformational landscape that furosemide occupies in solution. The solvent conditions should be chosen to match the molecule's charge state between the solution and solid states. The Δδ regression approach indicates whether or not correlations between Δδexperimental and Δδcalculated are statistically significant; the approach is differently sensitive to the popular root mean squared error (RMSE) method, being shown to exhibit a much greater dynamic range. An alternative method for estimating solution-state NMR chemical shifts by approximating the measured solution-state dynamic 3D behavior with an ensemble of 54 furosemide crystal structures (polymorphs and cocrystals) from the Cambridge Structural Database (CSD) was also successful in this case, suggesting new avenues for this method that may overcome its current dependency on the prior determination of solution dynamic 3D structures.
Collapse
Affiliation(s)
- Mohammed Rahman
- Department
of Physics, University of Warwick, Coventry CV4 7AL, U.K.
- Department
of Chemistry, University of Warwick, Coventry CV4 7AL, U.K.
| | | | | | - Leslie P. Hughes
- Oral
Product Development, Pharmaceutical Technology & Development, Operations, AstraZeneca, Macclesfield SK10 2NA, U.K.
| | - Helen Blade
- Oral
Product Development, Pharmaceutical Technology & Development, Operations, AstraZeneca, Macclesfield SK10 2NA, U.K.
| | - Jake Carson
- Mathematics
Institute at Warwick, University of Warwick, Coventry CV4 7AL, U.K.
| | - Ben P. Tatman
- Department
of Physics, University of Warwick, Coventry CV4 7AL, U.K.
- Department
of Chemistry, University of Warwick, Coventry CV4 7AL, U.K.
| | | | - Steven P. Brown
- Department
of Physics, University of Warwick, Coventry CV4 7AL, U.K.
| |
Collapse
|
2
|
Wang R, Mehdi S, Zou Z, Tiwary P. Is the Local Ion Density Sufficient to Drive NaCl Nucleation from the Melt and Aqueous Solution? J Phys Chem B 2024; 128:1012-1021. [PMID: 38262436 DOI: 10.1021/acs.jpcb.3c06735] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/25/2024]
Abstract
Even though nucleation is ubiquitous in different science and engineering problems, investigating nucleation is extremely difficult due to the complicated ranges of time and length scales involved. In this work, we simulate NaCl nucleation in both molten and aqueous environments using enhanced sampling of all-atom molecular dynamics with deep-learning-based estimation of reaction coordinates. By incorporating various structural order parameters and learning the reaction coordinate as a function thereof, we achieve significantly improved sampling relative to traditional ad hoc descriptions of what drives nucleation, particularly in an aqueous medium. Our results reveal a one-step nucleation mechanism in both environments, with reaction coordinate analysis highlighting the importance of local ion density in distinguishing solid and liquid states. However, although fluctuations in the local ion density are necessary to drive nucleation, they are not sufficient. Our analysis shows that near the transition states, descriptors such as enthalpy and local structure become crucial. Our protocol proposed here enables robust nucleation analysis and phase sampling and could offer insights into nucleation mechanisms for generic small molecules in different environments.
Collapse
Affiliation(s)
- Ruiyu Wang
- Institute for Physical Science and Technology, University of Maryland, College Park, Maryland 20742, United States
| | - Shams Mehdi
- Institute for Physical Science and Technology, University of Maryland, College Park, Maryland 20742, United States
- Biophysics Program, University of Maryland, College Park, Maryland 20742, United States
| | - Ziyue Zou
- Department of Chemistry and Biochemistry, University of Maryland, College Park, Maryland 20742, United States
| | - Pratyush Tiwary
- Institute for Physical Science and Technology, University of Maryland, College Park, Maryland 20742, United States
- Department of Chemistry and Biochemistry, University of Maryland, College Park, Maryland 20742, United States
| |
Collapse
|
3
|
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
|
4
|
Firaha D, Liu YM, van de Streek J, Sasikumar K, Dietrich H, Helfferich J, Aerts L, Braun DE, Broo A, DiPasquale AG, Lee AY, Le Meur S, Nilsson Lill SO, Lunsmann WJ, Mattei A, Muglia P, Putra OD, Raoui M, Reutzel-Edens SM, Rome S, Sheikh AY, Tkatchenko A, Woollam GR, Neumann MA. Predicting crystal form stability under real-world conditions. Nature 2023; 623:324-328. [PMID: 37938708 PMCID: PMC10632141 DOI: 10.1038/s41586-023-06587-3] [Citation(s) in RCA: 5] [Impact Index Per Article: 5.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/14/2022] [Accepted: 08/30/2023] [Indexed: 11/09/2023]
Abstract
The physicochemical properties of molecular crystals, such as solubility, stability, compactability, melting behaviour and bioavailability, depend on their crystal form1. In silico crystal form selection has recently come much closer to realization because of the development of accurate and affordable free-energy calculations2-4. Here we redefine the state of the art, primarily by improving the accuracy of free-energy calculations, constructing a reliable experimental benchmark for solid-solid free-energy differences, quantifying statistical errors for the computed free energies and placing both hydrate crystal structures of different stoichiometries and anhydrate crystal structures on the same energy landscape, with defined error bars, as a function of temperature and relative humidity. The calculated free energies have standard errors of 1-2 kJ mol-1 for industrially relevant compounds, and the method to place crystal structures with different hydrate stoichiometries on the same energy landscape can be extended to other multi-component systems, including solvates. These contributions reduce the gap between the needs of the experimentalist and the capabilities of modern computational tools, transforming crystal structure prediction into a more reliable and actionable procedure that can be used in combination with experimental evidence to direct crystal form selection and establish control5.
Collapse
Affiliation(s)
| | | | | | | | | | - Julian Helfferich
- Avant-garde Materials Simulation, Merzhausen, Germany
- JobRad, Freiburg, Germany
| | - Luc Aerts
- UCB Pharma SA, Chemin du Foriest, Braine-l'Alleud, Belgium
| | - Doris E Braun
- Institute of Pharmacy, University of Innsbruck, Innsbruck, Austria
| | - Anders Broo
- Data Science and Modelling, Pharmaceutical Sciences, R&D, AstraZeneca Gothenburg, Mölndal, Sweden
| | | | - Alfred Y Lee
- Merck, Analytical Research & Development, Rahway, NJ, USA
| | - Sarah Le Meur
- UCB Pharma SA, Chemin du Foriest, Braine-l'Alleud, Belgium
| | - Sten O Nilsson Lill
- Data Science and Modelling, Pharmaceutical Sciences, R&D, AstraZeneca Gothenburg, Mölndal, Sweden
| | | | - Alessandra Mattei
- Solid State Chemistry, Research & Development, AbbVie, North Chicago, IL, USA
| | | | - Okky Dwichandra Putra
- Early Product Development and Manufacturing, Pharmaceutical Sciences R&D, AstraZeneca Gothenburg, Mölndal, Sweden
| | | | - Susan M Reutzel-Edens
- Cambridge Crystallographic Data Centre, Cambridge, UK
- SuRE Pharma Consulting, Zionsville, IN, USA
| | - Sandrine Rome
- UCB Pharma SA, Chemin du Foriest, Braine-l'Alleud, Belgium
| | - Ahmad Y Sheikh
- Solid State Chemistry, Research & Development, AbbVie, North Chicago, IL, USA
| | - Alexandre Tkatchenko
- Department of Physics and Materials Science, University of Luxembourg, Luxembourg City, Luxembourg
| | | | | |
Collapse
|
5
|
Shamshina JL, Rogers RD. Ionic Liquids: New Forms of Active Pharmaceutical Ingredients with Unique, Tunable Properties. Chem Rev 2023; 123:11894-11953. [PMID: 37797342 DOI: 10.1021/acs.chemrev.3c00384] [Citation(s) in RCA: 4] [Impact Index Per Article: 4.0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/07/2023]
Abstract
This Review aims to summarize advances over the last 15 years in the development of active pharmaceutical ingredient ionic liquids (API-ILs), which make up a prospective game-changing strategy to overcome multiple problems with conventional solid-state drugs, for example, polymorphism. A critical part of the present Review is the collection of API-ILs and deep eutectic solvents (DESs) prepared to date. The Review covers rules for rational design of API-ILs and tools for API-IL formation, syntheses, and characterization. Nomenclature and ionic speciation, and the confusion that these may cause, are highlighted, particularly for speciation in both ILs and DESs of intermediate ionicity. We also highlight in vivo and in vitro pharmaceutical activity studies, with differences in pharmacokinetic/pharmacodynamic depending on ionicity of API-ILs. A brief overview is provided for the ILs used to deliver drugs, and the Review concludes with key prospects and roadblocks in translating API-ILs into pharmaceutical manufacturing.
Collapse
Affiliation(s)
- Julia L Shamshina
- Fiber and Biopolymer Research Institute (FBRI), Texas Tech University, Lubbock, Texas 79409, United States
| | - Robin D Rogers
- 525 Solutions, Inc., P.O. Box 2206, Tuscaloosa, Alabama 35403, United States
| |
Collapse
|
6
|
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
|
7
|
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
|
8
|
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
|
9
|
Newman JA, Iuzzolino L, Tan M, Orth P, Bruhn J, Lee AY. From Powders to Single Crystals: A Crystallographer's Toolbox for Small-Molecule Structure Determination. Mol Pharm 2022; 19:2133-2141. [PMID: 35576503 PMCID: PMC10152450 DOI: 10.1021/acs.molpharmaceut.2c00020] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Abstract
Although the crystal structures of small-molecule compounds are often determined from single-crystal X-ray diffraction (scXRD), recent advances in three-dimensional electron diffraction (3DED) and crystal structure prediction (CSP) methods promise to expand the structure elucidation toolbox available to the crystallographer. Herein, a comparative assessment of scXRD, 3DED, and CSP in combination with powder X-ray diffraction is carried out on two former drug candidate compounds and a multicomponent crystal of a key building block in the synthesis of gefapixant citrate.
Collapse
Affiliation(s)
- Justin A. Newman
- Department
of Analytical Research and Development, Merck & Co., Inc., 126 East Lincoln Avenue, Rahway, New Jersey 07065, United States
| | - Luca Iuzzolino
- Department
of Computational and Structural Chemistry, Merck & Co., Inc., Rahway, 126 East Lincoln Avenue, Rahway, New Jersey 07065, United States
| | - Melissa Tan
- Department
of Analytical Research and Development, Merck & Co., Inc., 126 East Lincoln Avenue, Rahway, New Jersey 07065, United States
| | - Peter Orth
- Department
of Computational and Structural Chemistry, Merck & Co., Inc., 2000 Galloping Hill Road, Kenilworth, New Jersey 07033, United States
| | - Jessica Bruhn
- Nanoimaging
Services, San Diego, California 92121, United States
| | - Alfred Y. Lee
- Department
of Analytical Research and Development, Merck & Co., Inc., 126 East Lincoln Avenue, Rahway, New Jersey 07065, United States
| |
Collapse
|
10
|
Reliable crystal structure predictions from first principles. Nat Commun 2022; 13:3095. [PMID: 35654882 PMCID: PMC9163189 DOI: 10.1038/s41467-022-30692-y] [Citation(s) in RCA: 5] [Impact Index Per Article: 2.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/17/2021] [Accepted: 05/10/2022] [Indexed: 11/28/2022] Open
Abstract
An inexpensive and reliable method for molecular crystal structure predictions (CSPs) has been developed. The new CSP protocol starts from a two-dimensional graph of crystal’s monomer(s) and utilizes no experimental information. Using results of quantum mechanical calculations for molecular dimers, an accurate two-body, rigid-monomer ab initio-based force field (aiFF) for the crystal is developed. Since CSPs with aiFFs are essentially as expensive as with empirical FFs, tens of thousands of plausible polymorphs generated by the crystal packing procedures can be optimized. Here we show the robustness of this protocol which found the experimental crystal within the 20 most stable predicted polymorphs for each of the 15 investigated molecules. The ranking was further refined by performing periodic density-functional theory (DFT) plus dispersion correction (pDFT+D) calculations for these 20 top-ranked polymorphs, resulting in the experimental crystal ranked as number one for all the systems studied (and the second polymorph, if known, ranked in the top few). Alternatively, the polymorphs generated can be used to improve aiFFs, which also leads to rank one predictions. The proposed CSP protocol should result in aiFFs replacing empirical FFs in CSP research. Developing theoretical frameworks to predict new polymorphs is highly desirable. Here the authors present an ab initio based force-field approach for crystal structure prediction offering a dramatic computational speed-up over fully ab initio schemes.
Collapse
|
11
|
Weatherby JA, Rumson AF, Price AJA, Otero de la Roza A, Johnson ER. A density-functional benchmark of vibrational free-energy corrections for molecular crystal polymorphism. J Chem Phys 2022; 156:114108. [DOI: 10.1063/5.0083082] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/11/2022] Open
Abstract
Many crystal structure prediction protocols only concern themselves with the electronic energy of molecular crystals. However, vibrational contributions to the free energy ( Fvib) can be significant in determining accurate stability rankings for crystal candidates. While force-field studies have been conducted to gauge the magnitude of these free-energy corrections, highly accurate results from quantum mechanical methods, such as density-functional theory (DFT), are desirable. Here, we introduce the PV17 set of 17 polymorphic pairs of organic molecular crystals, for which plane wave DFT is used to calculate the vibrational free energies and free-energy differences (Δ Fvib) between each pair. Our DFT results confirm that the vibrational free-energy corrections are small, having a mean value of 1.0 kJ/mol and a maximum value of 2.3 kJ/mol for the PV17 set. Furthermore, we assess the accuracy of a series of lower-cost DFT, semi-empirical, and force-field models for computing Δ Fvib that have been proposed in the literature. It is found that calculating Fvib using the Γ-point frequencies does not provide Δ Fvib values of sufficiently high quality. In addition, Δ Fvib values calculated using various approximate methods have mean absolute errors relative to our converged DFT results of equivalent or larger magnitude than the vibrational free-energy corrections themselves. Thus, we conclude that, in a crystal structure prediction protocol, it is preferable to forego the inclusion of vibrational free-energy corrections than to estimate them with any of the approximate methods considered here.
Collapse
Affiliation(s)
- Joseph A. Weatherby
- Department of Chemistry, Dalhousie University, 6274 Coburg Rd, Halifax, Nova Scotia B3H 4R2, Canada
| | - Adrian F. Rumson
- Department of Chemistry, Dalhousie University, 6274 Coburg Rd, Halifax, Nova Scotia B3H 4R2, Canada
| | - Alastair J. A. Price
- Department of Chemistry, Dalhousie University, 6274 Coburg Rd, Halifax, Nova Scotia B3H 4R2, Canada
| | - 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, 33006 Oviedo, Spain
| | - Erin R. Johnson
- Department of Chemistry, Dalhousie University, 6274 Coburg Rd, Halifax, Nova Scotia B3H 4R2, Canada
- Department of Physics and Atmospheric Science, Dalhousie University, 6310 Coburg Rd, Halifax, Nova Scotia B3H 4R2, Canada
| |
Collapse
|
12
|
Beran GJO, Wright SE, Greenwell C, Cruz-Cabeza AJ. The interplay of intra- and intermolecular errors in modeling conformational polymorphs. J Chem Phys 2022; 156:104112. [DOI: 10.1063/5.0088027] [Citation(s) in RCA: 4] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/26/2023] Open
Abstract
Conformational polymorphs of organic molecular crystals represent a challenging test for quantum chemistry because they require careful balancing of the intra- and intermolecular interactions. This study examines 54 molecular conformations from 20 sets of conformational polymorphs, along with the relative lattice energies and 173 dimer interactions taken from six of the polymorph sets. These systems are studied with a variety of van der Waals-inclusive density functionals theory models; dispersion-corrected spin-component-scaled second-order Møller–Plesset perturbation theory (SCS-MP2D); and domain local pair natural orbital coupled cluster singles, doubles, and perturbative triples [DLPNO-CCSD(T)]. We investigate how delocalization error in conventional density functionals impacts monomer conformational energies, systematic errors in the intermolecular interactions, and the nature of error cancellation that occurs in the overall crystal. The density functionals B86bPBE-XDM, PBE-D4, PBE-MBD, PBE0-D4, and PBE0-MBD are found to exhibit sizable one-body and two-body errors vs DLPNO-CCSD(T) benchmarks, and the level of success in predicting the relative polymorph energies relies heavily on error cancellation between different types of intermolecular interactions or between intra- and intermolecular interactions. The SCS-MP2D and, to a lesser extent, ωB97M-V models exhibit smaller errors and rely less on error cancellation. Implications for crystal structure prediction of flexible compounds are discussed. Finally, the one-body and two-body DLPNO-CCSD(T) energies taken from these conformational polymorphs establish the CP1b and CP2b benchmark datasets that could be useful for testing quantum chemistry models in challenging real-world systems with complex interplay between intra- and intermolecular interactions, a number of which are significantly impacted by delocalization error.
Collapse
Affiliation(s)
- Gregory J. O. Beran
- Department of Chemistry, University of California, Riverside, California 92521, USA
| | - Sarah E. Wright
- Department of Chemical Engineering and Analytical Science, University of Manchester, Manchester, United Kingdom
| | - Chandler Greenwell
- Department of Chemistry, University of California, Riverside, California 92521, USA
| | - Aurora J. Cruz-Cabeza
- Department of Chemical Engineering and Analytical Science, University of Manchester, Manchester, United Kingdom
| |
Collapse
|
13
|
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
|
14
|
Price SL. Progress in understanding crystallisation: a personal perspective. Faraday Discuss 2022; 235:569-581. [DOI: 10.1039/d2fd00077f] [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
After this discussion meeting, most participants felt that we do not understand crystallisation. However, in the 1980s I believe that most scientists would have considered that crystallisation was adequately understood....
Collapse
|
15
|
Dudek MK, Druzbicki K. Along the road to Crystal Structure Prediction (CSP) of pharmaceutical-like molecules. CrystEngComm 2022. [DOI: 10.1039/d1ce01564h] [Citation(s) in RCA: 2] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022]
Abstract
Computational methods used for predicting crystal structures of organic compounds are mature enough to be routinely used with many rigid and semi-rigid organic molecules. The usefulness of Crystal Structure Prediction...
Collapse
|
16
|
Beran GJO, Sugden IJ, Greenwell C, Bowskill DH, Pantelides CC, Adjiman CS. How many more polymorphs of ROY remain undiscovered. Chem Sci 2022; 13:1288-1297. [PMID: 35222912 PMCID: PMC8809489 DOI: 10.1039/d1sc06074k] [Citation(s) in RCA: 13] [Impact Index Per Article: 6.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/02/2021] [Accepted: 12/10/2021] [Indexed: 12/15/2022] Open
Abstract
With 12 crystal forms, 5-methyl-2-[(2-nitrophenyl)amino]-3-thiophenecabonitrile (a.k.a. ROY) holds the current record for the largest number of fully characterized organic crystal polymorphs. Four of these polymorph structures have been reported since 2019, raising the question of how many more ROY polymorphs await future discovery. Employing crystal structure prediction and accurate energy rankings derived from conformational energy-corrected density functional theory, this study presents the first crystal energy landscape for ROY that agrees well with experiment. The lattice energies suggest that the seven most stable ROY polymorphs (and nine of the twelve lowest-energy forms) on the Z′ = 1 landscape have already been discovered experimentally. Discovering any new polymorphs at ambient pressure will likely require specialized crystallization techniques capable of trapping metastable forms. At pressures above 10 GPa, however, a new crystal form is predicted to become enthalpically more stable than all known polymorphs, suggesting that further high-pressure experiments on ROY may be warranted. This work highlights the value of high-accuracy crystal structure prediction for solid-form screening and demonstrates how pragmatic conformational energy corrections can overcome the limitations of conventional density functionals for conformational polymorphs. Crystal structure prediction suggests that the low-energy polymorphs of ROY have already been found, but a new high-pressure form is predicted.![]()
Collapse
Affiliation(s)
- Gregory J. O. Beran
- Department of Chemistry, University of California Riverside, Riverside, CA 92521, USA
| | - Isaac J. Sugden
- Department of Chemical Engineering, Sargent Centre for Process Systems Engineering, Imperial College London, London, SW7 2AZ, UK
| | - Chandler Greenwell
- Department of Chemistry, University of California Riverside, Riverside, CA 92521, USA
| | - David H. Bowskill
- Department of Chemical Engineering, Sargent Centre for Process Systems Engineering, Imperial College London, London, SW7 2AZ, UK
| | - Constantinos C. Pantelides
- Department of Chemical Engineering, Sargent Centre for Process Systems Engineering, Imperial College London, London, SW7 2AZ, UK
| | - Claire S. Adjiman
- Department of Chemical Engineering, Sargent Centre for Process Systems Engineering, Imperial College London, London, SW7 2AZ, UK
| |
Collapse
|
17
|
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
|
18
|
Hong RS, Mattei A, Sheikh AY, Bhardwaj RM, Bellucci MA, McDaniel KF, Pierce MO, Sun G, Li S, Wang L, Mondal S, Ji J, Borchardt TB. Novel Physics-Based Ensemble Modeling Approach That Utilizes 3D Molecular Conformation and Packing to Access Aqueous Thermodynamic Solubility: A Case Study of Orally Available Bromodomain and Extraterminal Domain Inhibitor Lead Optimization Series. J Chem Inf Model 2021; 61:1412-1426. [PMID: 33661005 DOI: 10.1021/acs.jcim.0c01410] [Citation(s) in RCA: 6] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Abstract
Drug design with patient centricity for ease of administration and pill burden requires robust understanding of the impact of chemical modifications on relevant physicochemical properties early in lead optimization. To this end, we have developed a physics-based ensemble approach to predict aqueous thermodynamic crystalline solubility, with a 2D chemical structure as the input. Predictions for the bromodomain and extraterminal domain (BET) inhibitor series show very close match (0.5 log unit) with measured thermodynamic solubility for cases with low crystal anisotropy and good match (1 log unit) for high anisotropy structures. The importance of thermodynamic solubility is clearly demonstrated by up to a 4 log unit drop in solubility compared to kinetic (amorphous) solubility in some cases and implications thereof, for instance on human dose. We have also demonstrated that incorporating predicted crystal structures in thermodynamic solubility prediction is necessary to differentiate (up to 4 log unit) between solubility of molecules within the series. Finally, our physics-based ensemble approach provides valuable structural insights into the origins of 3-D conformational landscapes, crystal polymorphism, and anisotropy that can be leveraged for both drug design and development.
Collapse
Affiliation(s)
- Richard S Hong
- Research & Development, AbbVie Inc., 1 N Waukegan Road, North Chicago, Illinois 60064, United States
| | - Alessandra Mattei
- Research & Development, AbbVie Inc., 1 N Waukegan Road, North Chicago, Illinois 60064, United States
| | - Ahmad Y Sheikh
- Research & Development, AbbVie Inc., 1 N Waukegan Road, North Chicago, Illinois 60064, United States
| | - Rajni Miglani Bhardwaj
- Research & Development, AbbVie Inc., 1 N Waukegan Road, North Chicago, Illinois 60064, United States
| | - Michael A Bellucci
- XtalPi, Inc., 245 Main Street, Cambridge, Massachusetts 02142, United States
| | - Keith F McDaniel
- Research & Development, AbbVie Inc., 1 N Waukegan Road, North Chicago, Illinois 60064, United States
| | - M Olivia Pierce
- Schrödinger Inc., 120 W 45th Street, New York, New York 10036, United States
| | - Guangxu Sun
- XtalPi, Inc., 245 Main Street, Cambridge, Massachusetts 02142, United States
| | - Sizhu Li
- XtalPi, Inc., 245 Main Street, Cambridge, Massachusetts 02142, United States
| | - Lingle Wang
- Schrödinger Inc., 120 W 45th Street, New York, New York 10036, United States
| | - Sayan Mondal
- Schrödinger Inc., 120 W 45th Street, New York, New York 10036, United States
| | - Jianguo Ji
- Research & Development, AbbVie Inc., 1 N Waukegan Road, North Chicago, Illinois 60064, United States
| | - Thomas B Borchardt
- Research & Development, AbbVie Inc., 1 N Waukegan Road, North Chicago, Illinois 60064, United States
| |
Collapse
|
19
|
Mayo RA, Johnson ER. Improved quantitative crystal-structure comparison using powder diffractograms via anisotropic volume correction. CrystEngComm 2021. [DOI: 10.1039/d1ce01058a] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/18/2022]
Abstract
A new anisotropic volume correction improves quantitative crystal structure comparison. Benchmarking against the 6th crystal structure prediction blind test data results in identification of two previously uncredited matching structures.
Collapse
Affiliation(s)
- R. Alex Mayo
- Department of Chemistry, Dalhousie University, 6274 Coburg Road, PO Box 15000, Halifax, Nova Scotia, B3H 4R2, Canada
| | - Erin R. Johnson
- Department of Chemistry, Dalhousie University, 6274 Coburg Road, PO Box 15000, Halifax, Nova Scotia, B3H 4R2, Canada
| |
Collapse
|
20
|
Bier I, Marom N. Machine Learned Model for Solid Form Volume Estimation Based on Packing-Accessible Surface and Molecular Topological Fragments. J Phys Chem A 2020; 124:10330-10345. [DOI: 10.1021/acs.jpca.0c06791] [Citation(s) in RCA: 4] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/31/2022]
Affiliation(s)
- Imanuel Bier
- Department of Materials Science and Engineering, Carnegie Mellon University, Pittsburgh, Pennsylvania 15213, United States
| | - Noa Marom
- Department of Materials Science and Engineering, Carnegie Mellon University, Pittsburgh, Pennsylvania 15213, United States
- Department of Chemistry, Carnegie Mellon University, Pittsburgh, Pennsylvania 15213, United States
- Department of Physics, Carnegie Mellon University, Pittsburgh, Pennsylvania 15213, United States
| |
Collapse
|
21
|
Montis R, Davey RJ, Wright SE, Woollam GR, Cruz‐Cabeza AJ. Transforming Computed Energy Landscapes into Experimental Realities: The Role of Structural Rugosity. Angew Chem Int Ed Engl 2020; 59:20357-20360. [PMID: 32730676 PMCID: PMC7693212 DOI: 10.1002/anie.202006939] [Citation(s) in RCA: 13] [Impact Index Per Article: 3.3] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/13/2020] [Indexed: 11/14/2022]
Abstract
We exploit the possible link between structural surface roughness and difficulty of crystallisation. Polymorphs with smooth surfaces may nucleate and crystallise more readily than polymorphs with rough surfaces. The concept is applied to crystal structure prediction landscapes and reveals a promising complementary way of ranking putative crystal structures.
Collapse
Affiliation(s)
- Riccardo Montis
- Department of Chemical Engineering and Analytical ScienceThe University of ManchesterOxford RoadManchesterM13 9PLUK
| | - Roger J. Davey
- Department of Chemical Engineering and Analytical ScienceThe University of ManchesterOxford RoadManchesterM13 9PLUK
| | - Sarah E. Wright
- Department of Chemical Engineering and Analytical ScienceThe University of ManchesterOxford RoadManchesterM13 9PLUK
| | | | - Aurora J. Cruz‐Cabeza
- Department of Chemical Engineering and Analytical ScienceThe University of ManchesterOxford RoadManchesterM13 9PLUK
| |
Collapse
|
22
|
Reutzel-Edens SM, Bhardwaj RM. Crystal forms in pharmaceutical applications: olanzapine, a gift to crystal chemistry that keeps on giving. IUCRJ 2020; 7:955-964. [PMID: 33209310 PMCID: PMC7642794 DOI: 10.1107/s2052252520012683] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Figures] [Subscribe] [Scholar Register] [Received: 08/26/2020] [Accepted: 09/17/2020] [Indexed: 06/11/2023]
Abstract
This contribution reviews the efforts of many scientists around the world to discover and structurally characterize olanzapine crystal forms, clearing up inconsistencies in the scientific and patent literature and highlighting the challenges in identifying new forms amidst 60+ known polymorphs and solvates. Owing to its remarkable solid-state chemistry, olanzapine has emerged over the last three decades as a popular tool compound for developing new experimental and computational methods for enhanced molecular level understanding of solid-state structure, form diversity and crystallization outcomes. This article highlights the role of olanzapine in advancing the fundamental understanding of crystal forms, interactions within crystal structures, and growth units in molecular crystallization, as well as influencing the way in which drugs are developed today.
Collapse
Affiliation(s)
- Susan M. Reutzel-Edens
- Synthetic Molecule Design and Development, Eli Lilly and Company, Indianapolis, IN 46285, USA
| | - Rajni M. Bhardwaj
- Synthetic Molecule Design and Development, Eli Lilly and Company, Indianapolis, IN 46285, USA
| |
Collapse
|
23
|
Sun G, Jin Y, Li S, Yang Z, Shi B, Chang C, Abramov YA. Virtual Coformer Screening by Crystal Structure Predictions: Crucial Role of Crystallinity in Pharmaceutical Cocrystallization. J Phys Chem Lett 2020; 11:8832-8838. [PMID: 32969658 DOI: 10.1021/acs.jpclett.0c02371] [Citation(s) in RCA: 25] [Impact Index Per Article: 6.3] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 05/19/2023]
Abstract
One of the most popular strategies of the optimization of drug properties in the pharmaceutical industry appears to be a solid form changing into a cocrystalline form. A number of virtual screening approaches have been previously developed to allow a selection of the most promising cocrystal formers (coformers) for an experimental follow-up. A significant drawback of those methods is related to the lack of accounting for the crystallinity contribution to cocrystal formation. To address this issue, we propose in this study two virtual coformer screening approaches based on a modern cloud-computing crystal structure prediction (CSP) technology at a dispersion-corrected density functional theory (DFT-D) level. The CSP-based methods were for the first time validated on challenging cases of indomethacin and paracetamol cocrystallization, for which the previously developed approaches provided poor predictions. The calculations demonstrated a dramatic improvement of the virtual coformer screening performance relative to the other methods. It is demonstrated that the crystallinity contribution to the formation of paracetamol and indomethacin cocrystals is a dominant one and, therefore, should not be ignored in the virtual screening calculations. Our results encourage a broad utilization of the proposed CSP-based technology in the pharmaceutical industry as the only virtual coformer screening method that directly accounts for the crystallinity contribution.
Collapse
Affiliation(s)
- Guangxu Sun
- XtalPi Inc. (Shenzhen Jingtai Technology Co., Ltd.), Floor 4, No. 9, Hualian Industrial Zone, Dalang Street, Longhua District, Shenzhen 518100, China
| | - Yingdi Jin
- XtalPi Inc. (Shenzhen Jingtai Technology Co., Ltd.), Floor 4, No. 9, Hualian Industrial Zone, Dalang Street, Longhua District, Shenzhen 518100, China
| | - Sizhu Li
- XtalPi Inc. (Shenzhen Jingtai Technology Co., Ltd.), Floor 4, No. 9, Hualian Industrial Zone, Dalang Street, Longhua District, Shenzhen 518100, China
| | - Zhuocen Yang
- XtalPi Inc. (Shenzhen Jingtai Technology Co., Ltd.), Floor 4, No. 9, Hualian Industrial Zone, Dalang Street, Longhua District, Shenzhen 518100, China
| | - Baimei Shi
- XtalPi Inc. (Shenzhen Jingtai Technology Co., Ltd.), Floor 4, No. 9, Hualian Industrial Zone, Dalang Street, Longhua District, Shenzhen 518100, China
| | - Chao Chang
- XtalPi Inc. (Shenzhen Jingtai Technology Co., Ltd.), Floor 4, No. 9, Hualian Industrial Zone, Dalang Street, Longhua District, Shenzhen 518100, China
| | - 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
| |
Collapse
|
24
|
Montis R, Davey RJ, Wright SE, Woollam GR, Cruz‐Cabeza AJ. Transforming Computed Energy Landscapes into Experimental Realities: The Role of Structural Rugosity. Angew Chem Int Ed Engl 2020. [DOI: 10.1002/ange.202006939] [Citation(s) in RCA: 5] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/11/2022]
Affiliation(s)
- Riccardo Montis
- Department of Chemical Engineering and Analytical Science The University of Manchester Oxford Road Manchester M13 9PL UK
| | - Roger J. Davey
- Department of Chemical Engineering and Analytical Science The University of Manchester Oxford Road Manchester M13 9PL UK
| | - Sarah E. Wright
- Department of Chemical Engineering and Analytical Science The University of Manchester Oxford Road Manchester M13 9PL UK
| | | | - Aurora J. Cruz‐Cabeza
- Department of Chemical Engineering and Analytical Science The University of Manchester Oxford Road Manchester M13 9PL UK
| |
Collapse
|
25
|
Cui P, Svensson Grape E, Spackman PR, Wu Y, Clowes R, Day GM, Inge AK, Little MA, Cooper AI. An Expandable Hydrogen-Bonded Organic Framework Characterized by Three-Dimensional Electron Diffraction. J Am Chem Soc 2020; 142:12743-12750. [PMID: 32597187 PMCID: PMC7467715 DOI: 10.1021/jacs.0c04885] [Citation(s) in RCA: 45] [Impact Index Per Article: 11.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/27/2022]
Abstract
A molecular crystal of a 2-D hydrogen-bonded organic framework (HOF) undergoes an unusual structural transformation after solvent removal from the crystal pores during activation. The conformationally flexible host molecule, ABTPA, adapts its molecular conformation during activation to initiate a framework expansion. The microcrystalline activated phase was characterized by three-dimensional electron diffraction (3D ED), which revealed that ABTPA uses out-of-plane anthracene units as adaptive structural anchors. These units change orientation to generate an expanded, lower density framework material in the activated structure. The porous HOF, ABTPA-2, has robust dynamic porosity (SABET = 1183 m2 g-1) and exhibits negative area thermal expansion. We use crystal structure prediction (CSP) to understand the underlying energetics behind the structural transformation and discuss the challenges facing CSP for such flexible molecules.
Collapse
Affiliation(s)
- Peng Cui
- Department of Chemistry and Materials Innovation Factory, University of Liverpool, Liverpool L7 3NY, U.K
| | - Erik Svensson Grape
- Department of Materials and Environmental Chemistry, Stockholm University, Stockholm 106 91, Sweden
| | - Peter R Spackman
- Computational Systems Chemistry, School of Chemistry, University of Southampton, Southampton SO17 1BJ, U.K.,Leverhulme Research Centre for Functional Materials Design, University of Liverpool, Liverpool L7 3NY, U.K
| | - Yue Wu
- Department of Chemistry and Materials Innovation Factory, University of Liverpool, Liverpool L7 3NY, U.K
| | - Rob Clowes
- Department of Chemistry and Materials Innovation Factory, University of Liverpool, Liverpool L7 3NY, U.K
| | - Graeme M Day
- Computational Systems Chemistry, School of Chemistry, University of Southampton, Southampton SO17 1BJ, U.K.,Leverhulme Research Centre for Functional Materials Design, University of Liverpool, Liverpool L7 3NY, U.K
| | - A Ken Inge
- Department of Materials and Environmental Chemistry, Stockholm University, Stockholm 106 91, Sweden
| | - Marc A Little
- Department of Chemistry and Materials Innovation Factory, University of Liverpool, Liverpool L7 3NY, U.K
| | - Andrew I Cooper
- Department of Chemistry and Materials Innovation Factory, University of Liverpool, Liverpool L7 3NY, U.K.,Leverhulme Research Centre for Functional Materials Design, University of Liverpool, Liverpool L7 3NY, U.K
| |
Collapse
|
26
|
Shunnar AF, Dhokale B, Karothu DP, Bowskill DH, Sugden IJ, Hernandez HH, Naumov P, Mohamed S. Efficient Screening for Ternary Molecular Ionic Cocrystals Using a Complementary Mechanosynthesis and Computational Structure Prediction Approach. Chemistry 2020; 26:4752-4765. [PMID: 31793669 PMCID: PMC7187361 DOI: 10.1002/chem.201904672] [Citation(s) in RCA: 15] [Impact Index Per Article: 3.8] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/13/2019] [Indexed: 12/16/2022]
Abstract
The discovery of molecular ionic cocrystals (ICCs) of active pharmaceutical ingredients (APIs) widens the opportunities for optimizing the physicochemical properties of APIs whilst facilitating the delivery of multiple therapeutic agents. However, ICCs are often observed serendipitously in crystallization screens and the factors dictating their crystallization are poorly understood. We demonstrate here that mechanochemical ball milling is a versatile technique for the reproducible synthesis of ternary molecular ICCs in less than 30 min of grinding with or without solvent. Computational crystal structure prediction (CSP) calculations have been performed on ternary molecular ICCs for the first time and the observed crystal structures of all the ICCs were correctly predicted. Periodic dispersion-corrected DFT calculations revealed that all the ICCs are thermodynamically stable (mean stabilization energy=-2 kJ mol-1 ) relative to the crystallization of a physical mixture of the binary salt and acid. The results suggest that a combined mechanosynthesis and CSP approach could be used to target the synthesis of higher-order molecular ICCs with functional properties.
Collapse
Affiliation(s)
- Abeer F. Shunnar
- Department of ChemistryKhalifa University of Science and TechnologyP.O. Box 127788Abu DhabiUAE
| | - Bhausaheb Dhokale
- Department of ChemistryKhalifa University of Science and TechnologyP.O. Box 127788Abu DhabiUAE
| | | | - David H. Bowskill
- Molecular Systems Engineering GroupCentre for Process Systems EngineeringDepartment of Chemical EngineeringImperial College LondonLondonSW7 2AZUK
| | - Isaac J. Sugden
- Molecular Systems Engineering GroupCentre for Process Systems EngineeringDepartment of Chemical EngineeringImperial College LondonLondonSW7 2AZUK
| | - Hector H. Hernandez
- Department of Biomedical EngineeringCenter for Membrane and Advanced Water TechnologyKhalifa University of Science and TechnologyMasdar Campus P.O. Box 127788Abu DhabiUAE
| | - Panče Naumov
- New York University Abu DhabiP.O. Box 129188Abu DhabiUAE
| | - Sharmarke Mohamed
- Department of ChemistryKhalifa University of Science and TechnologyP.O. Box 127788Abu DhabiUAE
| |
Collapse
|
27
|
Liu Y, Gabriele B, Davey RJ, Cruz-Cabeza AJ. Concerning Elusive Crystal Forms: The Case of Paracetamol. J Am Chem Soc 2020; 142:6682-6689. [DOI: 10.1021/jacs.0c00321] [Citation(s) in RCA: 34] [Impact Index Per Article: 8.5] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Affiliation(s)
- Yumin Liu
- School of Chemical Engineering and Technology, Collaborative Innovation Center of Chemical Science and Chemical Engineering, Tianjin University, Tianjin 300072, P. R. China
- Department of Chemical Engineering and Analytical Science, University of Manchester, Manchester M13 9PL, United Kingdom
| | - Benjamin Gabriele
- Department of Chemical Engineering and Analytical Science, University of Manchester, Manchester M13 9PL, United Kingdom
| | - Roger J. Davey
- Department of Chemical Engineering and Analytical Science, University of Manchester, Manchester M13 9PL, United Kingdom
| | - Aurora J. Cruz-Cabeza
- Department of Chemical Engineering and Analytical Science, University of Manchester, Manchester M13 9PL, United Kingdom
| |
Collapse
|