1
|
A P V, O R S, T V V, G L P. Sublimation of pyridine derivatives: fundamental aspects and application for two-component crystal screening. Phys Chem Chem Phys 2024; 26:22558-22571. [PMID: 39150718 DOI: 10.1039/d4cp01442a] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 08/17/2024]
Abstract
The saturated vapour pressures of five heterocyclic compounds containing the pyridine fragment, namely, three isomers of aminopyridine (2-aminopyridine (2AmPy), 3-aminopyridine (3AmPy), and 4-aminopyridine (4AmPy)); 3-hydroxypyridine (3OHPy) and 2-(1H-imidazol-2-yl)pyridine (ImPy), were measured at appropriate temperature intervals using a transpiration (inert gas flow) method. The standard molar enthalpies, entropies, and Gibbs energies of sublimation for all the studied substances were determined. Among the compounds studied, the largest value of ΔH298sub was observed for ImPy. The influence of substitution and the effects of hydrogen bonds in the crystal lattices on sublimation parameters are discussed herein. The reliable dependences relating ΔG298sub to Tfus and ΔH298sub to ΔG298sub were plotted. A comparative analysis of several calculation schemes for the estimation of sublimation enthalpy and Gibbs free energy was carried out. Thermodynamic parameters obtained in this study were applied for the evaluation of cocrystallisation thermodynamic functions for two-component crystals (virtual screening) on the basis of the studied substituted pyridines.
Collapse
Affiliation(s)
- Voronin A P
- G.A. Krestov Institute of Solution Chemistry of the Russian Academy of Sciences, 1 Akademicheskaya St., Ivanovo, 153045, Russian Federation.
| | - Simonova O R
- G.A. Krestov Institute of Solution Chemistry of the Russian Academy of Sciences, 1 Akademicheskaya St., Ivanovo, 153045, Russian Federation.
| | - Volkova T V
- G.A. Krestov Institute of Solution Chemistry of the Russian Academy of Sciences, 1 Akademicheskaya St., Ivanovo, 153045, Russian Federation.
| | - Perlovich G L
- G.A. Krestov Institute of Solution Chemistry of the Russian Academy of Sciences, 1 Akademicheskaya St., Ivanovo, 153045, Russian Federation.
| |
Collapse
|
2
|
Ludík J, Kostková V, Kocian Š, Touš P, Štejfa V, Červinka C. First-Principles Models of Polymorphism of Pharmaceuticals: Maximizing the Accuracy-to-Cost Ratio. J Chem Theory Comput 2024; 20:2858-2870. [PMID: 38531828 PMCID: PMC11008097 DOI: 10.1021/acs.jctc.4c00099] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/24/2024] [Revised: 03/08/2024] [Accepted: 03/11/2024] [Indexed: 03/28/2024]
Abstract
Accuracy and sophistication of in silico models of structure, internal dynamics, and cohesion of molecular materials at finite temperatures increase over time. Applicability limits of ab initio polymorph ranking that would be feasible at reasonable costs currently represent crystals of moderately sized molecules (less than 20 nonhydrogen atoms) and simple unit cells (containing rather only one symmetry-irreducible molecule). Extending the applicability range of the underlying first-principles methods to larger systems with a real-life significance, and enabling to perform such computations in a high-throughput regime represent additional challenges to be tackled in computational chemistry. This work presents a novel composite method that combines the computational efficiency of density-functional tight-binding (DFTB) methods with the accuracy of density-functional theory (DFT). Being rooted in the quasi-harmonic approximation, it uses a cheap method to perform all of the costly scans of how static and dynamic characteristics of the crystal vary with respect to its volume. Such data are subsequently corrected to agree with a higher-level model, which must be evaluated only at a single volume of the crystal. It thus enables predictions of structural, cohesive, and thermodynamic properties of complex molecular materials, such as pharmaceuticals or organic semiconductors, at a fraction of the original computational cost. As the composite model retains the solid physical background, it suffers from a minimum accuracy deterioration compared to the full treatment with the costly approach. The novel methodology is demonstrated to provide consistent results for the structural and thermodynamic properties of real-life molecular crystals and their polymorph ranking.
Collapse
Affiliation(s)
- Jan Ludík
- Department of Physical Chemistry, University of Chemistry and Technology Prague, Technická 5, CZ-166 28 Prague 6, Czech Republic
| | - Veronika Kostková
- Department of Physical Chemistry, University of Chemistry and Technology Prague, Technická 5, CZ-166 28 Prague 6, Czech Republic
| | - Štefan Kocian
- Department of Physical Chemistry, University of Chemistry and Technology Prague, Technická 5, CZ-166 28 Prague 6, Czech Republic
| | - Petr Touš
- Department of Physical Chemistry, University of Chemistry and Technology Prague, Technická 5, CZ-166 28 Prague 6, Czech Republic
| | - Vojtěch Štejfa
- Department of Physical Chemistry, University of Chemistry and Technology Prague, Technická 5, CZ-166 28 Prague 6, Czech Republic
| | - Ctirad Červinka
- Department of Physical Chemistry, University of Chemistry and Technology Prague, Technická 5, CZ-166 28 Prague 6, Czech Republic
| |
Collapse
|
3
|
Hoja J, List A, Boese AD. Multimer Embedding Approach for Molecular Crystals up to Harmonic Vibrational Properties. J Chem Theory Comput 2024; 20:357-367. [PMID: 38109226 PMCID: PMC10782452 DOI: 10.1021/acs.jctc.3c01082] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/29/2023] [Revised: 11/16/2023] [Accepted: 11/20/2023] [Indexed: 12/20/2023]
Abstract
Accurate calculations of molecular crystals are crucial for drug design and crystal engineering. However, periodic high-level density functional calculations using hybrid functionals are often prohibitively expensive for the relevant systems. These expensive periodic calculations can be circumvented by the usage of embedding methods in which, for instance, the periodic calculation is only performed at a lower-cost level and then monomer energies and dimer interactions are replaced by those of the higher-level method. Herein, we extend such a multimer embedding approach to enable energy corrections for trimer interactions and the calculation of harmonic vibrational properties up to the dimer level. We evaluate this approach for the X23 benchmark set of molecular crystals by approximating a periodic hybrid density functional (PBE0+MBD) by embedding multimers into less expensive calculations using a generalized-gradient approximation functional (PBE+MBD). We show that trimer interactions are crucial for accurately approximating lattice energies within 1 kJ/mol and might also be needed for further improvement of lattice constants and hence cell volumes. Finally, the vibrational properties are already very well captured at the monomer and dimer level, making it possible to approximate vibrational free energies at room temperature within 1 kJ/mol.
Collapse
Affiliation(s)
- Johannes Hoja
- Department of Chemistry, University
of Graz, Heinrichstraße 28/IV, Graz 8010, Austria
| | - Alexander List
- Department of Chemistry, University
of Graz, Heinrichstraße 28/IV, Graz 8010, Austria
| | - A. Daniel Boese
- Department of Chemistry, University
of Graz, Heinrichstraße 28/IV, Graz 8010, Austria
| |
Collapse
|
4
|
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
|
5
|
Červinka C. Tuning the quasi-harmonic treatment of crystalline ionic liquids within the density functional theory. J Comput Chem 2021; 43:448-456. [PMID: 34958138 DOI: 10.1002/jcc.26804] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/12/2021] [Revised: 12/13/2021] [Accepted: 12/15/2021] [Indexed: 11/06/2022]
Abstract
Five ionic liquids are selected for benchmarking the performance of quasi-harmonic density functional theory (DFT) calculations of structural, phonon, and thermodynamic properties of their crystals. Data predicted by individual computational setups are sorted, establishing a distinct hierarchy among the first-principles approaches. PBE-D3 and B3LYP-D3 functionals are coupled with various plane wave and Gaussian-type orbital (GTO) basis sets. Propagation of the basis set superposition error and of the imperfections of both functionals into finite-temperature properties is discussed in detail. PBE-D3 together with a triple-zeta GTO basis set often yields the most accurate predictions of predicted molar volume and heat capacity with errors at 1% and 8%, respectively, representing the state-of-the-art for quasi-harmonic DFT calculations for crystalline ionic liquids. Fortuitous error cancellation between the basis-set superposition (overbinding) and PBE imperfection (overexpanding) strongly affects the overall accuracy, unlike the case of B3LYP/GTO calculations, impeding systematic convergence of the methodology towards higher accuracy.
Collapse
Affiliation(s)
- Ctirad Červinka
- Department of Physical Chemistry, University of Chemistry and Technology Prague, Prague, Czech Republic
| |
Collapse
|
6
|
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
|
7
|
Bowskill DH, Sugden IJ, Konstantinopoulos S, Adjiman CS, Pantelides CC. Crystal Structure Prediction Methods for Organic Molecules: State of the Art. Annu Rev Chem Biomol Eng 2021; 12:593-623. [PMID: 33770462 DOI: 10.1146/annurev-chembioeng-060718-030256] [Citation(s) in RCA: 20] [Impact Index Per Article: 6.7] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/09/2022]
Abstract
The prediction of the crystal structures that a given organic molecule is likely to form is an important theoretical problem of significant interest for the pharmaceutical and agrochemical industries, among others. As evidenced by a series of six blind tests organized over the past 2 decades, methodologies for crystal structure prediction (CSP) have witnessed substantial progress and have now reached a stage of development where they can begin to be applied to systems of practical significance. This article reviews the state of the art in general-purpose methodologies for CSP, placing them within a common framework that highlights both their similarities and their differences. The review discusses specific areas that constitute the main focus of current research efforts toward improving the reliability and widening applicability of these methodologies, and offers some perspectives for the evolution of this technology over the next decade.
Collapse
Affiliation(s)
- David H Bowskill
- Molecular Systems Engineering Group, Centre for Process Systems Engineering, Department of Chemical Engineering, and Institute for Molecular Science and Engineering, Imperial College London, South Kensington Campus, London SW7 2AZ, United Kingdom;
| | - Isaac J Sugden
- Molecular Systems Engineering Group, Centre for Process Systems Engineering, Department of Chemical Engineering, and Institute for Molecular Science and Engineering, Imperial College London, South Kensington Campus, London SW7 2AZ, United Kingdom;
| | - Stefanos Konstantinopoulos
- Molecular Systems Engineering Group, Centre for Process Systems Engineering, Department of Chemical Engineering, and Institute for Molecular Science and Engineering, Imperial College London, South Kensington Campus, London SW7 2AZ, United Kingdom;
| | - Claire S Adjiman
- Molecular Systems Engineering Group, Centre for Process Systems Engineering, Department of Chemical Engineering, and Institute for Molecular Science and Engineering, Imperial College London, South Kensington Campus, London SW7 2AZ, United Kingdom;
| | - Constantinos C Pantelides
- Molecular Systems Engineering Group, Centre for Process Systems Engineering, Department of Chemical Engineering, and Institute for Molecular Science and Engineering, Imperial College London, South Kensington Campus, London SW7 2AZ, United Kingdom;
| |
Collapse
|
8
|
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
|
9
|
Kamat K, Guo R, Reutzel-Edens SM, Price SL, Peters B. Diabat method for polymorph free energies: Extension to molecular crystals. J Chem Phys 2020; 153:244105. [PMID: 33380078 DOI: 10.1063/5.0024727] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
Lattice-switch Monte Carlo and the related diabat methods have emerged as efficient and accurate ways to compute free energy differences between polymorphs. In this work, we introduce a one-to-one mapping from the reference positions and displacements in one molecular crystal to the positions and displacements in another. Two features of the mapping facilitate lattice-switch Monte Carlo and related diabat methods for computing polymorph free energy differences. First, the mapping is unitary so that its Jacobian does not complicate the free energy calculations. Second, the mapping is easily implemented for molecular crystals of arbitrary complexity. We demonstrate the mapping by computing free energy differences between polymorphs of benzene and carbamazepine. Free energy calculations for thermodynamic cycles, each involving three independently computed polymorph free energy differences, all return to the starting free energy with a high degree of precision. The calculations thus provide a force field independent validation of the method and allow us to estimate the precision of the individual free energy differences.
Collapse
Affiliation(s)
- Kartik Kamat
- Department of Chemical Engineering, University of California-Santa Barbara, Santa Barbara, California 93106, USA
| | - Rui Guo
- Department of Chemistry, University College London, London WC1H 0AJ, United Kingdom
| | - Susan M Reutzel-Edens
- Small Molecule Design and Development, Eli Lilly and Company, Indianapolis, Indiana 46285, USA
| | - Sarah L Price
- Department of Chemistry, University College London, London WC1H 0AJ, United Kingdom
| | - Baron Peters
- Department of Chemical and Biomolecular Engineering, University of Illinois at Urbana-Champaign, Urbana, Illinois 61801, USA
| |
Collapse
|
10
|
Cook C, Beran GJO. Reduced-cost supercell approach for computing accurate phonon density of states in organic crystals. J Chem Phys 2020; 153:224105. [PMID: 33317313 DOI: 10.1063/5.0032649] [Citation(s) in RCA: 6] [Impact Index Per Article: 1.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/10/2023] Open
Abstract
Phonon contributions to organic crystal structures and thermochemical properties can be significant, but computing a well-converged phonon density of states with lattice dynamics and periodic density functional theory (DFT) is often computationally expensive due to the need for large supercells. Using semi-empirical methods like density functional tight binding (DFTB) instead of DFT can reduce the computational costs dramatically, albeit with noticeable reductions in accuracy. This work proposes approximating the phonon density of states via a relatively inexpensive DFTB supercell treatment of the phonon dispersion that is then corrected by shifting the individual phonon modes according to the difference between the DFT and DFTB phonon frequencies at the Γ-point. The acoustic modes are then computed at the DFT level from the elastic constants. In several small-molecule crystal test cases, this combined approach reproduces DFT thermochemistry with kJ/mol accuracy and 1-2 orders of magnitude less computational effort. Finally, this approach is applied to computing the free energy differences between the five crystal polymorphs of oxalyl dihydrazide.
Collapse
Affiliation(s)
- Cameron Cook
- 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
|
Abraham NS, Shirts MR. Statistical Mechanical Approximations to More Efficiently Determine Polymorph Free Energy Differences for Small Organic Molecules. J Chem Theory Comput 2020; 16:6503-6512. [PMID: 32877183 DOI: 10.1021/acs.jctc.0c00570] [Citation(s) in RCA: 8] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
Abstract
Methods to efficiently determine the relative stability of polymorphs of organic crystals are highly desired in crystal structure predictions (CSPs). Current methodologies include calculating the free energy of static lattice phonons, quasi-harmonic approximations (QHA), and computing the full thermodynamic cycle using replica exchange molecular dynamics (REMD). We found that 13 out of the 29 systems minimized from experimental crystal structures restructured to a lower energy minimum when heated and annealed using REMD, a phenomenon that QHA alone cannot capture. Here, we present a series of methods that are intermediate in accuracy and expense between QHA and computing the full thermodynamic cycle, which can save 42-80% of the computational cost and introduces, on this benchmark, a relatively small (0.16 ± 0.04 kcal/mol) error relative to the full thermodynamic cycle. In particular, a method that Boltzmann weights harmonic free energies from along the trajectory of REMD replica exchange appears to be an appropriate intermediate between QHA and the full thermodynamic cycle using MD when screening crystal polymorph stability.
Collapse
Affiliation(s)
- Nathan S Abraham
- Department of Chemical and Biological Engineering, University of Colorado Boulder, Boulder, Colorado 80309, United States
| | - Michael R Shirts
- Department of Chemical and Biological Engineering, University of Colorado Boulder, Boulder, Colorado 80309, United States
| |
Collapse
|
12
|
Červinka C, Štejfa V. Sublimation Properties of α,ω-Diamines Revisited from First-Principles Calculations. Chemphyschem 2020; 21:1184-1194. [PMID: 32243713 DOI: 10.1002/cphc.202000108] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.8] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/11/2020] [Revised: 04/02/2020] [Indexed: 11/06/2022]
Abstract
Sublimation enthalpies of alkane-α,ω-diamines exhibit an odd-even pattern within their homologous series. First-principles calculations coupled with the quasi-harmonic approximation for crystals and with the conformation mixing model for the ideal gas are used to explain this phenomenon from the theoretical point of view. Crystals of the odd and even alkane-α,ω-diamines distinctly differ in their packing motifs. However, first-principles calculations indicate that it is a delicate interplay of the cohesive forces, phonons, molecular vibrations and conformational equilibrium which governs the odd-even pattern of the sublimation enthalpies within the homologous series. High molecular flexibility of the alkane-α,ω-diamines predetermines higher sensitivity of the computational model to the quality of the optimized geometries and relative conformational energies. Performance of high-throughput computational methods, such as the density functional tight binding (DFTB, GFN2-xTB) and the explicitly correlated dispersion-corrected Møller-Plesset perturbative method (MP2C-F12), are benchmarked against the consistent state-of-the-art calculations of conformational energies and interaction energies, respectively.
Collapse
Affiliation(s)
- Ctirad Červinka
- Department of Physical Chemistry, University of Chemistry and Technology Prague, Technická 5, 166 28, Prague 6, Czech Republic
| | - Vojtěch Štejfa
- Department of Physical Chemistry, University of Chemistry and Technology Prague, Technická 5, 166 28, Prague 6, Czech Republic
| |
Collapse
|
13
|
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
|
14
|
DFT calculations of the structural, electronic, optical and vibrational properties of anhydrous orthorhombic L-threonine crystals. COMPUT THEOR CHEM 2019. [DOI: 10.1016/j.comptc.2019.112621] [Citation(s) in RCA: 8] [Impact Index Per Article: 1.6] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/11/2023]
|
15
|
Abstract
Since the introduction of the fragment molecular orbital method 20 years ago, fragment-based approaches have occupied a small but growing niche in quantum chemistry. These methods decompose a large molecular system into subsystems small enough to be amenable to electronic structure calculations, following which the subsystem information is reassembled in order to approximate an otherwise intractable supersystem calculation. Fragmentation sidesteps the steep rise (with respect to system size) in the cost of ab initio calculations, replacing it with a distributed cost across numerous computer processors. Such methods are attractive, in part, because they are easily parallelizable and therefore readily amenable to exascale computing. As such, there has been hope that distributed computing might offer the proverbial "free lunch" in quantum chemistry, with the entrée being high-level calculations on very large systems. While fragment-based quantum chemistry can count many success stories, there also exists a seedy underbelly of rarely acknowledged problems. As these methods begin to mature, it is time to have a serious conversation about what they can and cannot be expected to accomplish in the near future. Both successes and challenges are highlighted in this Perspective.
Collapse
Affiliation(s)
- John M Herbert
- Department of Chemistry and Biochemistry, The Ohio State University, Columbus, Ohio 43210, USA
| |
Collapse
|
16
|
Červinka C, Fulem M. Cohesive properties of the crystalline phases of twenty proteinogenic α-aminoacids from first-principles calculations. Phys Chem Chem Phys 2019; 21:18501-18515. [PMID: 31411212 DOI: 10.1039/c9cp03102b] [Citation(s) in RCA: 15] [Impact Index Per Article: 3.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022]
Abstract
Cohesive properties (lattice and cohesive energy of the crystal and corresponding sublimation enthalpy) of the complete set of twenty enantiopure anhydrous proteinogenic amino acids are investigated using first-principles calculations. In contrast to neutral amino acid molecules in the vapor phase, all amino acids form crystals in their zwitterionic form. Therefore, reliable ab initio calculations of the proton transfer energy are an indispensable step of such calculations. Simplifying procedures, designed to rationalize the computational cost of the quasi-harmonic approximation, which proves too demanding if performed fully at the given quantum level of theory, are presented and tested. For this purpose, atomic multipoles (up to the quadrupoles) for the amoeba force field are parametrized for all amino acid zwitterions. While the calculated lattice energies of the amino acids range from 235-458 kJ mol-1 in absolute value, the proton transfer energies typically amount to 100-220 kJ mol-1, which translates to sublimation enthalpies ranging from 117-202 kJ mol-1, appreciably exceeding the sublimation enthalpy values common for nonionic molecular crystals. Critically assessed experimental data on sublimation enthalpies are used as a benchmark for comparison of the data calculated in this work. Cohesive properties of most amino acids calculated in this work, combining the PBE-D3(BJ)/PAW and CCSD(T)-F12/aug-cc-pVDZ levels of theory used for predictions of the lattice energies and of the proton transfer energies, respectively, exhibit a reasonable agreement with the experiment. At the same time, this work contains the first published data on cohesive properties for several enantiopure amino acids.
Collapse
Affiliation(s)
- Ctirad Červinka
- Department of Physical Chemistry, University of Chemistry and Technology, Prague, Technická 5, CZ-166 28 Prague 6, Czech Republic.
| | | |
Collapse
|
17
|
McKinley JL, Beran GJO. Improving Predicted Nuclear Magnetic Resonance Chemical Shifts Using the Quasi-Harmonic Approximation. J Chem Theory Comput 2019; 15:5259-5274. [PMID: 31442040 DOI: 10.1021/acs.jctc.9b00481] [Citation(s) in RCA: 8] [Impact Index Per Article: 1.6] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/22/2022]
Abstract
Ab initio nuclear magnetic resonance chemical shift prediction plays an important role in the determination or validation of crystal structures. The ability to predict chemical shifts more accurately can translate to increased confidence in the resulting chemical shift or structural assignments. Standard electronic structure predictions for molecular crystal structures neglect thermal expansion, which can lead to an appreciable underestimation of the molar volumes. This study examines this volume error and its impact on 68 13C- and 28 15N-predicted chemical shifts taken from 20 molecular crystals. It assesses the ability to recover more realistic room-temperature crystal structures using the quasi-harmonic approximation and how refining the structures impacts the chemical shifts. Several pharmaceutical molecular crystals are also examined in more detail. On the whole, accounting for quasi-harmonic expansion changes the 13C and 15N chemical shifts by 0.5 and 1.0 ppm on average. This, in turn, reduces the root-mean-square errors relative to experiment by 0.3 ppm for 13C and 0.7 ppm for 15N. Although the statistical impacts are modest, changes in individual chemical shifts can reach multiple ppm. Accounting for thermal expansion in molecular crystal chemical shift prediction may not be needed routinely, but the systematic trend toward improved accuracy with the experiment could be useful in cases where discrimination between structural candidates is challenging, as in the pharmaceutical theophylline.
Collapse
Affiliation(s)
- Jessica L McKinley
- Department of Chemistry , University of California , Riverside , California 92521 , United States
| | - Gregory J O Beran
- Department of Chemistry , University of California , Riverside , California 92521 , United States
| |
Collapse
|
18
|
Červinka C, Klajmon M, Štejfa V. Cohesive Properties of Ionic Liquids Calculated from First Principles. J Chem Theory Comput 2019; 15:5563-5578. [DOI: 10.1021/acs.jctc.9b00625] [Citation(s) in RCA: 17] [Impact Index Per Article: 3.4] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Affiliation(s)
- Ctirad Červinka
- Department of Physical Chemistry, University of Chemistry and Technology Prague, Technická 5, CZ-166 28 Prague 6, Czech Republic
| | - Martin Klajmon
- Department of Physical Chemistry, University of Chemistry and Technology Prague, Technická 5, CZ-166 28 Prague 6, Czech Republic
| | - Vojtěch Štejfa
- Department of Physical Chemistry, University of Chemistry and Technology Prague, Technická 5, CZ-166 28 Prague 6, Czech Republic
| |
Collapse
|
19
|
Dolgonos GA, Hoja J, Boese AD. Revised values for the X23 benchmark set of molecular crystals. Phys Chem Chem Phys 2019; 21:24333-24344. [DOI: 10.1039/c9cp04488d] [Citation(s) in RCA: 15] [Impact Index Per Article: 3.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022]
Abstract
A revised reference value set for molecular crystals: X23b; new cell volumes and lattice energies including volumetric expansion due to zero-point energy and thermal effects.
Collapse
Affiliation(s)
| | - Johannes Hoja
- Institute of Chemistry
- University of Graz
- 8010 Graz
- Austria
| | | |
Collapse
|
20
|
Hoja J, Ko HY, Neumann MA, Car R, DiStasio RA, Tkatchenko A. Reliable and practical computational description of molecular crystal polymorphs. SCIENCE ADVANCES 2019; 5:eaau3338. [PMID: 30746448 PMCID: PMC6357866 DOI: 10.1126/sciadv.aau3338] [Citation(s) in RCA: 91] [Impact Index Per Article: 18.2] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Received: 05/30/2018] [Accepted: 11/28/2018] [Indexed: 05/12/2023]
Abstract
Reliable prediction of the polymorphic energy landscape of a molecular crystal would yield profound insight into drug development in terms of the existence and likelihood of late-appearing polymorphs. However, the computational prediction of molecular crystal polymorphs is highly challenging due to the high dimensionality of conformational and crystallographic space accompanied by the need for relative free energies to within 1 kJ/mol per molecule. In this study, we combine the most successful crystal structure sampling strategy with the most successful first-principles energy ranking strategy of the latest blind test of organic crystal structure prediction methods. Specifically, we present a hierarchical energy ranking approach intended for the refinement of relative stabilities in the final stage of a crystal structure prediction procedure. Such a combined approach provides excellent stability rankings for all studied systems and can be applied to molecular crystals of pharmaceutical importance.
Collapse
Affiliation(s)
- Johannes Hoja
- Physics and Materials Science Research Unit, University of Luxembourg, L-1511 Luxembourg, Luxembourg
| | - Hsin-Yu Ko
- Department of Chemistry, Princeton University, Princeton, NJ 08544, USA
| | - Marcus A. Neumann
- Avant-garde Materials Simulation Deutschland GmbH, Alte Str. 2, 79249 Merzhausen, Germany
| | - Roberto Car
- Department of Chemistry, Princeton University, Princeton, NJ 08544, USA
| | - Robert A. DiStasio
- Department of Chemistry and Chemical Biology, Cornell University, Ithaca, NY 14853, USA
| | - Alexandre Tkatchenko
- Physics and Materials Science Research Unit, University of Luxembourg, L-1511 Luxembourg, Luxembourg
- Corresponding author.
| |
Collapse
|
21
|
Beran GJO. Solid state photodimerization of 9-tert-butyl anthracene ester produces an exceptionally metastable polymorph according to first-principles calculations. CrystEngComm 2019. [DOI: 10.1039/c8ce01985a] [Citation(s) in RCA: 13] [Impact Index Per Article: 2.6] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/26/2022]
Abstract
Molecular crystal engineering seeks to tune the material properties by controlling the crystal packing.
Collapse
|
22
|
Červinka C, Beran GJO. Towards reliable ab initio sublimation pressures for organic molecular crystals - are we there yet? Phys Chem Chem Phys 2019; 21:14799-14810. [PMID: 31225538 DOI: 10.1039/c9cp01572h] [Citation(s) in RCA: 12] [Impact Index Per Article: 2.4] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022]
Abstract
Knowledge of molecular crystal sublimation equilibrium data is vital in many industrial processes, but this data can be difficult to measure experimentally for low-volatility species. Theoretical prediction of sublimation pressures could provide a useful supplement to experiment, but the exponential temperature dependence of sublimation (or any saturated vapor) pressure curve makes this challenging. An uncertainty of only a few percent in the sublimation enthalpy or entropy can propagate to an error in the sublimation pressure exceeding several orders of magnitude for a given temperature interval. Despite this fundamental difficulty, this paper performs some of the first ab initio predictions of sublimation pressure curves. Four simple molecular crystals (ethane, methanol, benzene, and imidazole) have been selected for a case study showing the currently achievable accuracy of quantum chemistry calculations. Fragment-based ab initio techniques and the quasi-harmonic approximation are used for calculations of cohesive and phonon properties of the crystals, while the vapor phase is treated by the ideal gas model. Ab initio sublimation pressure curves for model compounds are compared against their experimental counterparts. The computational uncertainties are estimated, weak points of the computational methodology are identified, and further improvements are proposed.
Collapse
Affiliation(s)
- Ctirad Červinka
- Department of Physical Chemistry, University of Chemistry and Technology Prague, Technická 5, CZ-166 28 Prague 6, Czech Republic.
| | - Gregory J O Beran
- Department of Chemistry, University of California, Riverside, California 92521, USA
| |
Collapse
|
23
|
Červinka C, Beran GJO. Ab initio prediction of the polymorph phase diagram for crystalline methanol. Chem Sci 2018; 9:4622-4629. [PMID: 29899955 PMCID: PMC5969506 DOI: 10.1039/c8sc01237g] [Citation(s) in RCA: 41] [Impact Index Per Article: 6.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/15/2018] [Accepted: 04/13/2018] [Indexed: 11/21/2022] Open
Abstract
Organic crystals frequently adopt multiple distinct polymorphs exhibiting different properties. The ability to predict not only what crystal forms might occur, but under what experimental thermodynamic conditions those polymorphs are stable would be immensely valuable to the pharmaceutical industry and others. Starting only from knowledge of the experimental crystal structures, this study successfully predicts the methanol crystal polymorph phase diagram from first-principles quantum chemistry, mapping out the thermodynamic regions of stability for three polymorphs over the range 0-400 K and 0-6 GPa. The agreement between the predicted and experimental phase diagrams corresponds to predicting the relative polymorph free energies to within ∼0.5 kJ mol-1 accuracy, which is achieved by employing fragment-based second-order Møller-Plesset perturbation theory and coupled cluster theory plus a quasi-harmonic treatment of the phonons.
Collapse
Affiliation(s)
- Ctirad Červinka
- Department of Physical Chemistry , University of Chemistry and Technology Prague , Technická 5 , CZ-16628 Prague 6 , Czech Republic .
| | - Gregory J O Beran
- Department of Chemistry , University of California , Riverside , California 92521 , USA .
| |
Collapse
|