Templated Nucleation of Clotrimazole and Ketoprofen on Polymer Substrates.

The use of different template surfaces in crystallization experiments can directly influence the nucleation kinetics, crystal growth, and morphology of active pharmaceutical ingredients (APIs). Consequently, templated nucleation is an attractive approach to enhance crystal nucleation kinetics and preferentially nucleate desired crystal polymorphs for solid-form drug molecules, particularly large and flexible molecules that are difficult to crystallize. Herein, we investigate the effect of polymer templates on the crystal nucleation of clotrimazole and ketoprofen with both experiments and computational methods. Crystallization was carried out in toluene solvent for both APIs with a template library consisting of 12 different polymers. In complement to the experimental studies, we developed a computational workflow based on molecular dynamics (MD) and derived descriptors from the simulations to score and rank API-polymer interactions. The descriptors were used to measure the energy of interaction (EOI), hydrogen bonding, and rugosity (surface roughness) similarity between the APIs and polymer templates. We used a variety of machine learning models (14 in total) along with these descriptors to predict the crystallization outcome of the polymer templates. We found that simply rank-ordering the polymers by their API-polymer interaction energy descriptors yielded 92% accuracy in predicting the experimental outcome for clotrimazole and ketoprofen. The most accurate machine learning model for both APIs was found to be a random forest model. Using these models, we were able to predict the crystallization outcomes for all polymers. Additionally, we have performed a feature importance analysis using the trained models and found that the most predictive features are the energy descriptors. These results demonstrate that API-polymer interaction energies are correlated with heterogeneous crystallization outcomes.

Tale of Two Polymorphs: Investigating the Structural Differences and Dynamic Relationship between Nirmatrelvir Solid Forms (Paxlovid).

Two anhydrous polymorphs of the novel antiviral medicine nirmatrelvir were discovered during the development of Paxlovid, Pfizer's oral Covid-19 treatment. A comprehensive experimental and computational approach was necessary to distinguish the two closely related polymorphs, herein identified as Forms 1 and 4. This approach paired experimental methods, including powder X-ray diffraction and single-crystal X-ray diffraction, solid-state experimental methods, thermal analysis, solid-state nuclear magnetic resonance and Raman spectroscopy with computational investigations comprising crystal structure prediction, Gibbs free energy calculations, and molecular dynamics simulations of the polymorphic transition. Forms 1 and 4 were ultimately determined to be enantiotropically related polymorphs with Form 1 being the stable form above the transition temperature of ∼17 °C and designated as the nominated form for drug development. The work described in this paper shows the importance of using highly specialized orthogonal approaches to elucidate the subtle differences in structure and properties of similar solid-state forms. This synergistic approach allowed for unprecedented speed in bringing Paxlovid to patients in record time amidst the pandemic.

Improved Description of Intra- and Intermolecular Interactions through Dispersion-Corrected Second-Order Møller-Plesset Perturbation Theory.

ConspectusThe quantum chemical modeling of organic crystals and other molecular condensed-phase problems requires computationally affordable electronic structure methods which can simultaneously describe intramolecular conformational energies and intermolecular interactions accurately. To achieve this, we have developed a spin-component-scaled, dispersion-corrected second-order Møller-Plesset perturbation theory (SCS-MP2D) model. SCS-MP2D augments canonical MP2 with a dispersion correction which removes the uncoupled Hartree-Fock dispersion energy present in canonical MP2 and replaces it with a more reliable coupled Kohn-Sham treatment, all evaluated within the framework of Grimme's D3 dispersion model. The spin-component scaling is then used to improve the description of the residual (nondispersion) portion of the correlation energy.The SCS-MP2D model improves upon earlier corrected MP2 models in a few ways. Compared to the highly successful dispersion-corrected MP2C model, which is based solely on intermolecular perturbation theory, the SCS-MP2D dispersion correction improves the description of both inter- and intramolecular interactions. The dispersion correction can also be evaluated with trivial computational cost, and nuclear analytic gradients are computed readily to enable geometry optimizations. In contrast to earlier spin-component scaling MP2 models, the optimal spin-component scaling coefficients are only mildly sensitive to the choice of training data, and a single global parametrization of the model can describe both thermochemistry and noncovalent interactions.The resulting dispersion-corrected, spin-component-scaled MP2 (SCS-MP2D) model predicts conformational energies and intermolecular interactions with accuracy comparable to or better than those of many range-separated and double-hybrid density functionals, as is demonstrated on a variety of benchmark tests. Among the functionals considered here, only the revDSD-PBEP86-D3(BJ) functional gives consistently smaller errors in benchmark tests. The results presented also hint that further improvements of SCS-MP2D may be possible through a more robust fitting procedure for the seven empirical parameters.To demonstrate the performance of SCS-MP2D further, several applications to molecular crystal problems are presented. The three chosen examples all represent cases where density-driven delocalization error causes GGA or hybrid density functionals to artificially stabilize crystals exhibiting more extended π-conjugation. Our pragmatic strategy addresses the delocalization error by combining a periodic density functional theory (DFT) treatment of the infinite lattice with intramolecular/conformational energy corrections computed with SCS-MP2D. For the anticancer drug axitinib, applying the SCS-MP2D conformational energy correction produces crystal polymorph stabilities that are consistent with experiment, in contrast to earlier studies. For the crystal structure prediction of the ROY molecule, so named for its colorful red, orange, and yellow crystals, this approach leads to the first plausible crystal energy landscape, and it reveals that the lowest-energy polymorphs have already been found experimentally. Finally, in the context of photomechanical crystals, which transform light into mechanical work, these techniques are used to predict the structural transformations and extract design principles for maximizing the work performed.

Cocrystal Synthesis through Crystal Structure Prediction.

Crystal structure prediction (CSP) is an invaluable tool in the pharmaceutical industry because it allows to predict all the possible crystalline solid forms of small-molecule active pharmaceutical ingredients. We have used a CSP-based cocrystal prediction method to rank ten potential cocrystal coformers by the energy of the cocrystallization reaction with an antiviral drug candidate, MK-8876, and a triol process intermediate, 2-ethynylglyclerol. For MK-8876, the CSP-based cocrystal prediction was performed retrospectively and successfully predicted the maleic acid cocrystal as the most likely cocrystal to be observed. The triol is known to form two different cocrystals with 1,4-diazabicyclo[2.2.2]octane (DABCO), but a larger solid form landscape was desired. CSP-based cocrystal screening predicted the triol-DABCO cocrystal as rank one, while a triol-l-proline cocrystal was predicted as rank two. Computational finite-temperature corrections enabled determination of relative crystallization propensities of the triol-DABCO cocrystals with different stoichiometries and prediction of the triol-l-proline polymorphs in the free-energy landscape. The triol-l-proline cocrystal was obtained during subsequent targeted cocrystallization experiments and was found to exhibit an improved melting point and deliquescence behavior over the triol-free acid, which could be considered as an alternative solid form in the synthesis of islatravir.

Effect of Polymer Additives on the Crystal Habit of Metformin HCl.

The crystal habit can have a profound influence on the physical properties of crystalline materials, and thus controlling the crystal morphology is of great practical relevance across many industries. Herein, this work investigates the effect of polymer additives on the crystal habit of metformin HCl with both experiments and computational methods with the aim of developing a combined screening approach for crystal morphology engineering. Crystallization experiments of metformin HCl are conducted in methanol and in an isopropanol-water mixture (8:2 V/V). Polyethylene glycol, polyvinylpyrrolidone, Tween80, and hydroxypropyl methylcellulose polymer additives are used in low concentrations (1-2% w/w) in the experiments to study the effect they have on modifying the crystal habit. Additionally, this work has developed computational methods to characterize the morphology "landscape" and quantifies the overall effect of solvent and additives on the predicted crystal habits. Further analysis of the molecular dynamics simulations is used to rationalize the effect of additives on specific crystal faces. This work demonstrates that the effects of additives on the crystal habit are a result of their absorption and interactions with the slow growing {100} and {020} faces.

Absolute Configuration Determination of Chiral API Molecules by MicroED Analysis of Cocrystal Powders Formed Based on Cocrystal Propensity Prediction Calculations.

Establishing the absolute configuration of chiral active pharmaceutical ingredients (APIs) is of great importance. Single crystal X-ray diffraction (scXRD) has traditionally been the method of choice for such analysis, but scXRD requires the growth of large crystals, which can be challenging. Here, we present a method for determining absolute configuration that does not rely on the growth of large crystals. By examining microcrystals formed with chiral probes (small chiral compounds such as amino acids), absolute configuration can be unambiguously determined by microcrystal electron diffraction (MicroED). Our streamlined method employs three steps: (1) virtual screening to identify promising chiral probes, (2) experimental cocrystal screening and (3) structure determination by MicroED and absolute configuration assignment. We successfully applied this method to analyze two chiral API molecules currently on the market for which scXRD was not used to determine absolute configuration.

The interplay of intra- and intermolecular errors in modeling conformational polymorphs.

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.

How many more polymorphs of ROY remain undiscovered.

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.

Spin-component-scaled and dispersion-corrected second-order Møller-Plesset perturbation theory: a path toward chemical accuracy.

Second-order Møller-Plesset perturbation theory (MP2) provides a valuable alternative to density functional theory for modeling problems in organic and biological chemistry. However, MP2 suffers from known limitations in the description of van der Waals (London) dispersion interactions and reaction thermochemistry. Here, a spin-component-scaled, dispersion-corrected MP2 model (SCS-MP2D) is proposed that addresses these weaknesses. The dispersion correction, which is based on Grimme's D3 formalism, replaces the uncoupled Hartree-Fock dispersion inherent in MP2 with a more robust coupled Kohn-Sham treatment. The spin-component scaling of the residual MP2 correlation energy then reduces the remaining errors in the model. This two-part correction strategy solves the problem found in earlier spin-component-scaled MP2 models where completely different spin-scaling parameters were needed for describing reaction energies versus intermolecular interactions. Results on 18 benchmark data sets and two challenging potential energy curves demonstrate that SCS-MP2D considerably improves upon the accuracy of MP2 for intermolecular interactions, conformational energies, and reaction energies. Its accuracy and computational cost are competitive with state-of-the-art density functionals such as DSD-BLYP-D3(BJ), revDSD-PBEP86-D3(BJ), ωB97X-V, and ωB97M-V for systems with ∼100 atoms.

Predicting Density Functional Theory-Quality Nuclear Magnetic Resonance Chemical Shifts via Δ-Machine Learning.

First-principles prediction of nuclear magnetic resonance chemical shifts plays an increasingly important role in the interpretation of experimental spectra, but the required density functional theory (DFT) calculations can be computationally expensive. Promising machine learning models for predicting chemical shieldings in general organic molecules have been developed previously, though the accuracy of those models remains below that of DFT. The present study demonstrates how much higher accuracy chemical shieldings can be obtained via the Δ-machine learning approach, with the result that the errors introduced by the machine learning model are only one-half to one-third the errors expected for DFT chemical shifts relative to experiment. Specifically, an ensemble of neural networks is trained to correct PBE0/6-31G chemical shieldings up to the target level of PBE0/6-311+G(2d,p). It can predict 1H, 13C, 15N, and 17O chemical shieldings with root-mean-square errors of 0.11, 0.70, 1.69, and 2.47 ppm, respectively. At the same time, the Δ-machine learning approach is 1-2 orders of magnitude faster than the target large-basis calculations. It is also demonstrated that the machine learning model predicts experimental solution-phase NMR chemical shifts in drug molecules with only modestly worse accuracy than the target DFT model. Finally, the ability to estimate the uncertainty in the predicted shieldings based on variations within the ensemble of neural network models is also assessed.

Overcoming the difficulties of predicting conformational polymorph energetics in molecular crystals via correlated wavefunction methods.

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.

Accurate Noncovalent Interactions via Dispersion-Corrected Second-Order Møller-Plesset Perturbation Theory.

Noncovalent interactions govern many important areas of chemistry, ranging from biomolecules to molecular crystals. Here, an accurate and computationally inexpensive dispersion-corrected second-order Møller-Plesset perturbation theory model (MP2D) is presented. MP2D recasts the highly successful dispersion-corrected MP2C model in a framework based on Grimme's D3 dispersion correction, combining Grimme's D3 dispersion coefficients with new analogous uncoupled Hartree-Fock ones and five global empirical parameters. MP2D is faster than MP2C, and unlike MP2C, it is suitable for geometry optimizations and can describe both intra- and intermolecular noncovalent interactions with high accuracy. MP2D approaches the accuracy of higher-level ab initio wave function techniques and out-performs a widely used hybrid dispersion-corrected density functional on a range of intermolecular, intramolecular, and thermochemical benchmarks.