Efficient Screening of Coformers for Active Pharmaceutical Ingredient Cocrystallization

Sublimation of pyridine derivatives: fundamental aspects and application for two-component crystal screening

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.

Investigation of the Storage and Stability as Well as the Dissolution Rate of Novel Ilaprazole/Xylitol Cocrystal

Reflux esophagitis, a treatment for gastric ulcers known as Ilaprazole (Ila), is not stable during storage and handling at room temperature, requiring storage at 5 degrees Celsius. In this study, to address these issues with Ila, coformers rich in oxygen (O) and hydroxyl (OH) groups capable of forming hydrogen bonds with were selected. These coformers included Xylitol (Xyl), Meglumine (Meg), Nicotinic acid (Nic), L-Aspartic acid (Asp), and L-Glutamic acid (Glu). A 1:1 physical mixture of Ila and each coformer was prepared, and the potential for cocrystal formation was predicted using differential scanning calorimetry (DSC) screening. The results indicated the potential for cocrystal formation in the Ila/Xyl physical mixture. Subsequently, Ila and Xyl were mixed in ethyl acetate (EA) in a 1:1 ratio, and after 28 h of slurry, the formation of Ila/Xyl cocrystal was confirmed through solid-state CP/MAS 13C NMR spectrum analysis, showing intermolecular hydrogen bonding and conformational changes. Furthermore, the 1:1 ratio of Ila/Xyl cocrystal was confirmed through solution-state NMR (1H, 13C, and 2D) molecular structure analysis. To assess the stability of Ila/Xyl cocrystal at room temperature, it was stored and compared with Ila at 25 ± 2 °C and relative humidity (RH) of 65 ± 5% over three months. The results showed that the purity of Ila/Xyl cocrystal remained at 99.8% from the initial purity of 99.75% over the three months, while Ila was predicted to decrease from an initial 99.8% purity to 90% after three months. Additionally, at 25 ± 2 °C and RH 65 ± 5%, a specific impurity B in Ila/Xyl cocrystal was observed to be 0.03% over three months, whereas Ila was predicted to increase from an initial 0.032% to 2.28% after three months. To evaluate the dissolution rate of Ila/Xyl cocrystal, a formulation was prepared and compared with Ila at pH 10, with a dosage equivalent to 10 mg of Ila. The results showed that Ila/Xyl cocrystal reached 55% within 15 min and 100% within 45 min, while Ila was predicted to reach 32% at 15 min and 100% only after 60 min. However, overall, the Ila/Xyl cocrystal showed results equivalent to or exceeding the dissolution rate of Ila. Therefore, it is predicted that the Ila/Xyl cocrystal will maximize its effectiveness as a more convenient crystal structure for formulation development, allowing storage and preservation at room temperature without the need for the problematic 5 °C refrigeration during ambient conditions and storage, addressing the issues associated with Ila.

Frontiers of molecular crystal structure prediction for pharmaceuticals and functional organic materials

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.

Crystal engineering: from promise to delivery

Twenty years ago, I wrote a Chem. Commun. feature article entitled "Crystal Engineering: where from? Where to?": an update is in order. In this Highlight I argue that molecular crystal engineering, one of the areas of fast development of the field, has definitely reached the stage of "delivering the goods": new functional materials assembled via non-covalent interactions and/or improved properties of existing materials. As a proof of concept, the crystal engineering approach to tackle two contemporary emergencies, namely, urea fertilizer degradation and development of antimicrobial resistance by pathogens, is discussed and application-driven examples are provided.

In silico co-crystal design: Assessment of the latest advances

Pharmaceutical co-crystals represent a growing class of crystal forms in the context of pharmaceutical science. They are attractive to pharmaceutical scientists because they significantly expand the number of crystal forms that exist for an active pharmaceutical ingredient and can lead to improvements in physicochemical properties of clinical relevance. At the same time, machine learning is finding its way into all areas of drug discovery and delivers impressive results. In this review, we attempt to provide an overview of machine learning, deep learning and network-based recommendation approaches applied to pharmaceutical co-crystallization. We also present crystal structure prediction as an alternative to machine learning approaches.

Systematic Exploration of Structural Topologies in Hydrogen-Bonded Supramolecular Assemblies of Citric Acid with Different Heterocyclic Compounds

Hydrogen-bonded supramolecular assemblies of citric acid, CA, with some heterocyclic compounds (N-donor and N-oxide)-acridine (acr), phenazine (phenz), 1,10-phenanthroline (110phen), 1,7-phenanthroline (17phen), 4,7-phenanthroline (47phen), 1,4-diazabicyclo[2.2.2]octane (dabco), and 4,4'-bipyridyl-N,N'-dioxide (bpydo)-have been reported. Among these, only the N-donors phenz and N-oxide (bpydo) form neutral co-crystals, while the others form salts owing to the deprotonation of -COOH. Thus, depending on the nature of the aggregate (salt/co-crystal), recognition between the co-formers is established through O-H···N/N⁺-H···O/N⁺H···O^-heteromeric hydrogen bonding. Additionally, CA molecules establish homomeric interactions mediated by O-H···O hydrogen bonds. Moreover, CA forms a cyclic network with the co-formers or on its own, with a noteworthy feature of formation of host-guest networks in the assemblies with acr and phenz (solvated). In the assembly of acr, the CA molecules form a host network and captivate acr molecules as guest species, while in the case of phenz assembly, both the co-formers together encapsulate the solvent in the channels. However, the observed cyclic networks in the other structures form three-dimensional topologies in the form of ladders, a sandwich, lamellar layers, and interpenetrated networks. The structural features of the ensembles are evaluated unequivocally by the single-crystal X-ray diffraction method, while the homogeneity and phase purity are evaluated by using the powder X-ray diffraction method and differential scanning calorimetry. Further, conformational analysis of CA molecules reveals three types of conformations-T-shape (type I), syn-anti (type II), and syn (type III) as also observed in the literature for other CA co-crystals. In addition, the strength of the intermolecular interactions is quantified by performing Hirshfeld analysis.

Polymorphism of Carbamazepine Pharmaceutical Cocrystal: Structural Analysis and Solubility Performance

Polymorphism is a common phenomenon among single- and multicomponent molecular crystals that has a significant impact on the contemporary drug development process. A new polymorphic form of the drug carbamazepine (CBZ) cocrystal with methylparaben (MePRB) in a 1:1 molar ratio as well as the drug's channel-like cocrystal containing highly disordered coformer molecules have been obtained and characterized in this work using various analytical methods, including thermal analysis, Raman spectroscopy, and single-crystal and high-resolution synchrotron powder X-ray diffraction. Structural analysis of the solid forms revealed a close resemblance between novel form II and previously reported form I of the [CBZ + MePRB] (1:1) cocrystal in terms of hydrogen bond networks and overall packing arrangements. The channel-like cocrystal was found to belong to a distinct family of isostructural CBZ cocrystals with coformers of similar size and shape. Form I and form II of the 1:1 cocrystal appeared to be related by a monotropic relationship, with form II being proven to be the thermodynamically more stable phase. The dissolution performance of both polymorphs in aqueous media was significantly enhanced when compared with parent CBZ. However, considering the superior thermodynamic stability and consistent dissolution profile, the discovered form II of the [CBZ + MePRB] (1:1) cocrystal seems a more promising and reliable solid form for further pharmaceutical development.

Exploring the Supramolecular Interactions and Thermal Stability of Dapsone:Bipyridine Cocrystals by Combining Computational Chemistry with Experimentation

The application of computational screening methodologies based on H-bond propensity scores, molecular complementarity, molecular electrostatic potentials, and crystal structure prediction has guided the discovery of novel cocrystals of dapsone and bipyridine (DDS:BIPY). The experimental screen, which included mechanochemical and slurry experiments as well as the contact preparation, resulted in four cocrystals, including the previously known DDS:4,4'-BIPY (2:1, CC₄₄-B) cocrystal. To understand the factors governing the formation of the DDS:2,2'-BIPY polymorphs (1:1, CC₂₂-A and CC₂₂-B) and the two DDS:4,4'-BIPY cocrystal stoichiometries (1:1 and 2:1), different experimental conditions (such as the influence of solvent, grinding/stirring time, etc.) were tested and compared with the virtual screening results. The computationally generated (1:1) crystal energy landscapes had the experimental cocrystals as the lowest energy structures, although distinct cocrystal packings were observed for the similar coformers. H-bonding scores and molecular electrostatic potential maps correctly indicated cocrystallization of DDS and the BIPY isomers, with a higher likelihood for 4,4'-BIPY. The molecular conformation influenced the molecular complementarity results, predicting no cocrystallization for 2,2'-BIPY with DDS. The crystal structures of CC₂₂-A and CC₄₄-A were solved from powder X-ray diffraction data. All four cocrystals were fully characterized by a range of analytical techniques, including powder X-ray diffraction, infrared spectroscopy, hot-stage microscopy, thermogravimetric analysis, and differential scanning calorimetry. The two DDS:2,2'-BIPY polymorphs are enantiotropically related, with form B being the stable polymorph at room temperature (RT) and form A being the higher temperature form. Form B is metastable but kinetically stable at RT. The two DDS:4,4'-BIPY cocrystals are stable at room conditions; however, at higher temperatures, CC₄₄-A transforms to CC₄₄-B. The cocrystal formation enthalpy order, derived from the lattice energies, was calculated as follows: CC₄₄-B > CC₄₄-A > CC₂₂-A.

Virtual Screening, Structural Analysis, and Formation Thermodynamics of Carbamazepine Cocrystals

In this study, the existing set of carbamazepine (CBZ) cocrystals was extended through the successful combination of the drug with the positional isomers of acetamidobenzoic acid. The structural and energetic features of the CBZ cocrystals with 3- and 4-acetamidobenzoic acids were elucidated via single-crystal X-ray diffraction followed by QTAIMC analysis. The ability of three fundamentally different virtual screening methods to predict the correct cocrystallization outcome for CBZ was assessed based on the new experimental results obtained in this study and data available in the literature. It was found that the hydrogen bond propensity model performed the worst in distinguishing positive and negative results of CBZ cocrystallization experiments with 87 coformers, attaining an accuracy value lower than random guessing. The method that utilizes molecular electrostatic potential maps and the machine learning approach named CCGNet exhibited comparable results in terms of prediction metrics, albeit the latter resulted in superior specificity and overall accuracy while requiring no time-consuming DFT computations. In addition, formation thermodynamic parameters for the newly obtained CBZ cocrystals with 3- and 4-acetamidobenzoic acids were evaluated using temperature dependences of the cocrystallization Gibbs energy. The cocrystallization reactions between CBZ and the selected coformers were found to be enthalpy-driven, with entropy terms being statistically different from zero. The observed difference in dissolution behavior of the cocrystals in aqueous media was thought to be caused by variations in their thermodynamic stability.

Exploring the Cocrystal Landscape of Posaconazole by Combining High-Throughput Screening Experimentation with Computational Chemistry

The development of multicomponent crystal forms, such as cocrystals, represents a means to enhance the dissolution and absorption properties of poorly water-soluble drug compounds. However, the successful discovery of new pharmaceutical cocrystals remains a time- and resource-consuming process. This study proposes the use of a combined computational-experimental high-throughput approach as a tool to accelerate and improve the efficiency of cocrystal screening exemplified by posaconazole. First, we employed the COSMOquick software to preselect and rank cocrystal candidates (coformers). Second, high-throughput crystallization experiments (HTCS) were conducted on the selected coformers. The HTCS results were successfully reproduced by liquid-assisted grinding and reaction crystallization, ultimately leading to the synthesis of thirteen new posaconazole cocrystals (7 anhydrous, 5 hydrates, and 1 solvate). The posaconazole cocrystals were characterized by PXRD, ¹H NMR, Fourier transform-Raman, thermogravimetry-Fourier transform infrared spectroscopy, and differential scanning calorimetry. In addition, the prediction performance of COSMOquick was compared to that of two alternative knowledge-based methods: molecular complementarity (MC) and hydrogen bond propensity (HBP). Although HBP does not perform better than random guessing for this case study, both MC and COSMOquick show good discriminatory ability, suggesting their use as a potential virtual tool to improve cocrystal screening.