Virtual Coformer Screening by Crystal Structure Predictions: Crucial Role of Crystallinity in Pharmaceutical Cocrystallization

Computational predictions of cocrystal formation: A benchmark study of 28 assemblies comparing five methods from high-throughput to advanced models

Cocrystals are assemblies of more than one type of molecule stabilized through noncovalent interactions. They are promising materials for improved drug formulation in which the stability, solubility, or biocompatibility of the active pharmaceutical ingredient (API) is improved by including a coformer. In this work, a range of density functional theory (DFT) and density functional tight binding (DFTB) models are systematically compared for their ability to predict the lattice enthalpy of a broad range of existing pharmaceutically relevant cocrystals. These range from cocrystals containing model compounds 4,4'-bipyridine and oxalic acid to those with the well benchmarked APIs of aspirin and paracetamol, all tested with a large set of alternative coformers. For simple cocrystals, there is a general consensus in lattice enthalpy calculated by the different DFT models. For the cocrystals with API coformers the cocrystals, enthalpy predictions depend strongly on the DFT model. The significantly lighter DFTB models predict unrealistic values of lattice enthalpy even for simple cocrystals.

Driving Forces in the Formation of Paracetamol Cocrystals and Solvate with Naphthalene, Quinoline and Acridine

Paracetamol is an important analgesic and antipyretic drug showing poor tabletability. Among the various approaches used to improve this property, understanding the forces that govern the crystal packing is revealed to be crucial. We prepared three stable compounds: (par)2∙(nap) (1), (par)∙(quin) (2), and (par)∙(acr) (3) (nap-naphthalene, quin-quinoline, acr-acridine) being cocrystals or solvate. The structural studies showed that all the reported compounds are composed of alternately arranged layers of paracetamol and coformer. Several supramolecular motifs in the paracetamol layer were identified: R44(22) in (1); R64(20) and R22(8) in (2); and R22(8), R42(12), and R44(26) rings in (3). The stability of the crystal network was studied by interactions analysis performed by Hirshfeld surface and fingerprint approaches and the energy between the closest units in the crystal network was calculated. It showed that the strongest interactions were found between blocks connected by N-H⋯O=C and O-H⋯O/N hydrogen bonds due to an important coulombic factor. The dispersive energy becomes important for tail-to-tail (and head-to-tail) arranged paracetamol units, and it prevails in the case of stacking interactions between coformer molecules. The importance of dispersive forces increases with the size of the aromatic system of the coformer. XAS studies confirmed the successful preparation of compounds and provided some details about electron structure.

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.

Cocrystal Prediction of Nifedipine Based on the Graph Neural Network and Molecular Electrostatic Potential Surface

Nifedipine (NIF) is a dihydropyridine calcium channel blocker primarily used to treat conditions such as hypertension and angina. However, its low solubility and low bioavailability limit its effectiveness in clinical practice. Here, we developed a cocrystal prediction model based on Graph Neural Networks (CocrystalGNN) for the screening of cocrystals with NIF. And scoring 50 coformers using CocrystalGNN. To validate the reliability of the model, we used another prediction method, Molecular Electrostatic Potential Surface (MEPS), to verify the prediction results. Subsequently, we performed a second validation using experiments. The results indicate that our model achieved high performance. Ultimately, cocrystals of NIF were successfully obtained and all cocrystals exhibited better solubility and dissolution characteristics compared to the parent drug. This study lays a solid foundation for combining virtual prediction with experimental screening to discover novel water-insoluble drug cocrystals.

Thermodynamic Stability Is a Poor Indicator of Cocrystallization in Models of Organic Molecules

Cocrystallizing a given molecule with another can be useful for adjusting the physical properties of molecules in the solid state. However, most combinations of molecules do not readily cocrystallize but form either one-component crystals or amorphous solids. Computational methods of crystal structure prediction can, in principle, identify the thermodynamically stable cocrystal and thus predict if molecules will cocrystallize or not. However, the pronounced polymorphism and tendency of many organic molecules to form disordered solids suggest that kinetic factors can play an important role in cocrystallization. The question remains: if a binary system of molecules has a thermodynamically stable cocrystal, will it indeed cocrystallize? To address this question, we simulate the crystallization of more than 2600 distinct pairs of chiral model molecules of similar size in 2D and calculate accurate crystal energy landscapes for all of them. Our analysis shows that thermodynamic criteria alone are unreliable in the prediction of cocrystallization. While the vast majority of cocrystals that form in our simulations are thermodynamically favorable, most coformer systems that have a thermodynamically stable cocrystal do not cocrystallize. We furthermore show that cocrystallization rates increase 3-fold when coformers are used that do not form well-ordered single-component crystals. Our results suggest that kinetic factors of cocrystallization are decisive in many cases.

Computational and Experimental Screening Approaches to Aripiprazole Salt Crystallization

INTRODUCTION

The screening of multicomponent crystal system (MCC) is a key method for improving physicochemical properties of active pharmaceutical ingredients (APIs). The challenges associated with experimental salt screening include a large number of potential counterions and solvent systems and tendency to undergo disproportionation to produce free form during crystallization. These challenges may be mitigated by a combination of experimental and computational approaches to salt screening. The goal of this study is to evaluate performance of the counterion screening methods and propose and validate novel approaches to virtual solvent screening for MCC crystallization.

METHODS

The actual performance of the ΔpKa > 3 rule for counterion selection was validated using multiple screenings reports. Novel computational models for virtual solvent screening to avoid MCC incongruent crystallization were proposed. Using the ΔpKa rule, 10 acid counterions were selected for experimental aripiprazole (APZ) salt screening using 10 organic solvents. The experimental results were used to validate the proposed novel virtual solvent screen models.

RESULTS

Experimental APZ salt screening resulted in a total of eight MCCs which included glucuronate, mesylate, oxalate, tartrate, salicylate and mandelate. The new model to virtually screen solvents provided a general agreement with APZ experimental findings in terms of selecting the optimal solvent for MCC crystallization.

CONCLUSION

The rational selection of counterions and organic solvents for MCC crystallization was presented using combined novel computational model as well as experimental studies. The current virtual solvent screen model was successfully implemented and validated which can be easily applied to newly discovered APIs.

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.

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.

Recent Advances in Co-Former Screening and Formation Prediction of Multicomponent Solid Forms of Low Molecular Weight Drugs

Multicomponent solid forms of low molecular weight drugs, such as co-crystals, salts, and co-amorphous systems, are a result of the combination of an active pharmaceutical ingredient (API) with a pharmaceutically acceptable co-former. These solid forms can enhance the physicochemical and pharmacokinetic properties of APIs, making them increasingly interesting and important in recent decades. Nevertheless, predicting the formation of API multicomponent solid forms in the early stages of formulation development can be challenging, as it often requires significant time and resources. To address this, empirical and computational methods have been developed to help screen for potential co-formers more efficiently and accurately, thus reducing the number of laboratory experiments needed. This review provides a comprehensive overview of current screening and prediction methods for the formation of API multicomponent solid forms, covering both crystalline states (co-crystals and salts) and amorphous forms (co-amorphous). Furthermore, it discusses recent advances and emerging trends in prediction methods, with a particular focus on artificial intelligence.

Pharmaceutical Cocrystal Discovery via 3D-SMINBR: A New Network Recommendation Tool Augmented by 3D Molecular Conformations

Cocrystals have significant potential in various fields such as chemistry, material, and medicine. For instance, pharmaceutical cocrystals have the ability to address issues associated with physicochemical and biopharmaceutical properties. However, it can be challenging to find proper coformers to form cocrystals with drugs of interest. Herein, a new in silico tool called 3D substructure-molecular-interaction network-based recommendation (3D-SMINBR) has been developed to address this problem. This tool first integrated 3D molecular conformations with a weighted network-based recommendation model to prioritize potential coformers for target drugs. In cross-validation, the performance of 3D-SMINBR surpassed the 2D substructure-based predictive model SMINBR in our previous study. Additionally, the generalization capability of 3D-SMINBR was confirmed by testing on unseen cocrystal data. The practicality of this tool was further demonstrated by case studies on cocrystal screening of armillarisin A (Arm) and isoimperatorin (iIM). The obtained Arm-piperazine and iIM-salicylamide cocrystals present improved solubility and dissolution rate compared to their parent drugs. Overall, 3D-SMINBR augmented by 3D molecular conformations would be a useful network-based tool for cocrystal discovery. A free web server for 3D-SMINBR can be freely accessed at http://lmmd.ecust.edu.cn/netcorecsys/.

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.

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.

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.

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.

Efficient Screening of Coformers for Active Pharmaceutical Ingredient Cocrystallization

Controlling the physical properties of solid forms for active pharmaceutical ingredients (APIs) through cocrystallization is an important part of drug product development. However, it is difficult to know a priori which coformers will form cocrystals with a given API, and the current state-of-the-art for cocrystal discovery involves an expensive, time-consuming, and, at the early stages of pharmaceutical development, API material-limited experimental screen. We propose a systematic, high-throughput computational approach primarily aimed at identifying API/coformer pairs that are unlikely to lead to experimentally observable cocrystals and can therefore be eliminated with only a brief experimental check, from any experimental investigation. On the basis of a well-established crystal structure prediction (CSP) methodology, the proposed approach derives its efficiency by not requiring any expensive quantum mechanical calculations beyond those already performed for the CSP investigation of the neat API itself. The approach and assumptions are tested through a computational investigation on 30 potential 1:1 multicomponent systems (cocrystals and solvate) involving 3 active pharmaceutical ingredients and 9 coformers and one solvent. This is complemented with a detailed experimental investigation of all 30 pairs, which led to the discovery of five new cocrystals (three API-coformer combinations, a polymorphic cocrystal example, and one with different stoichiometries) and a cis-aconitic acid polymorph. The computational approach indicates that, for some APIs, a significant proportion of all potential API/coformer pairs could be investigated with only a brief experimental check, thereby saving considerable experimental effort.

Emerging Landscape of Computational Modeling in Pharmaceutical Development

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.

Insight of the various in silico screening techniques developed for assortment of cocrystal formers and their thermodynamic characterization

Most of the widely used drugs have problems associated with their oral bioavailability either due to their poor aqueous solubility or due to their poor permeability. Co-crystallization is an efficient and economically feasible approach that offers a great opportunity for improvement in physicochemical properties such as solubility, stability, and bioavailability of such type of therapeutic agent. Selection of the best co-former plays a major role in co-crystallization. Various approaches have been developed for the selection of suitable co-formers with API. In recent years in silico screening, a computational tool paying more attention for screening of co-formers has been developed. Numerous approaches can be used for in silico screening such as the Autodocking tool, COSMORS, COSMOTHERM, etc. Autodocking can predict several numbers of co-former effectively screened in silico method to identify a suitable co-former with an API. Prediction of solubility and dissolution is also important for the development of co-crystal. In this review, we discuss in silico screening of coformer and thermodynamic approaches to determine the dissolution and solubility of co-crystal specially with reference to the drugs belonging to BCS class II group.

Novel Cocrystals of Vonoprazan: Machine Learning-Assisted Discovery

Vonoprazan (VPZ) is the first-in-class potassium-competitive acid blocker (P-CAB), and has many advantages over proton pump inhibitors (PPIs). It is administered as a fumarate salt for the treatment of acid-related diseases, including reflux esophagitis, gastric ulcer, and duodenal ulcer, and for eradication of Helicobacter pylori. To discover novel cocrystals of VPZ, we adopted an artificial neural network (ANN)-based machine learning model as a virtual screening tool that can guide selection of the most promising coformers for VPZ cocrystals. Experimental screening by liquid-assisted grinding (LAG) confirmed that 8 of 19 coformers selected by the ANN model were likely to create new solid forms with VPZ. Structurally similar benzenediols and benzenetriols, i.e., catechol (CAT), resorcinol (RES), hydroquinone (HYQ), and pyrogallol (GAL), were used as coformers to obtain phase pure cocrystals with VPZ by reaction crystallization. We successfully prepared and characterized three novel cocrystals: VPZ–RES, VPZ–CAT, and VPZ–GAL. VPZ–RES had the highest solubility among the novel cocrystals studied here, and was even more soluble than the commercially available fumarate salt of VPZ in solution at pH 6.8. In addition, novel VPZ cocrystals had superior stability in aqueous media than VPZ fumarates, demonstrating their potential for improved pharmaceutical performance.

The extensive solid-form landscape of sulfathiazole: hydrogen-bond topology and node shape

Patterns of hydrogen bonds are described in a set of 101 crystal structures containing sulfathiazole. Topological analysis of the hydrogen-bond nets is augmented by comparison of the shapes of the nodes extracted from each net.

The extensive solid-form landscape of sulfathiazole: geometrical similarity and interaction energies

Sulfathiazole shows one of the most extensive solid-form landscapes known to date for an active pharmaceutical ingredient. A standardised structure set of 5 polymorphs, 59 co-crystals, 29 salts, and 3 other structures is established.

A study to discover novel pharmaceutical cocrystals of pelubiprofen with a machine learning approach compared

Pelubiprofen (PF), a biopharmaceutical classification system (BCS) class II non-steroidal anti-inflammatory drug, has been on the market only in its crystalline form. To discover the first cocrystal form(s) of the...

Molecular interactions of indomethacin and amino acids: Computational approach

Molecular interactions of indomethacin (IND) and amino acids (AA) were investigated in this work by employing the computational approaches. To this aim, the models of IND-AA were stabilized by performing density functional theory (DFT) calculations yielding the most favorable configurations regarding the energy values. Next, the approach of quantum theory of atoms in molecules (QTAIM) was used to recognize the roles of interactions and their significance in the bimolecular models. The results of interaction energies indicate that tryptophan (TRP) and phenylalanine (PHE) could be considered for participating in strong interactions with the IND substance. The results of QTAIM indicated that not only the electronegative atomic centers, but also homo-atomic centers could play significant roles in formations of IND-AA bimolecular models.

Coupling complementary strategy to flexible graph neural network for quick discovery of coformer in diverse co-crystal materials

Cocrystal engineering have been widely applied in pharmaceutical, chemistry and material fields. However, how to effectively choose coformer has been a challenging task on experiments. Here we develop a graph neural network (GNN) based deep learning framework to quickly predict formation of the cocrystal. In order to capture main driving force to crystallization from 6819 positive and 1052 negative samples reported by experiments, a feasible GNN framework is explored to integrate important prior knowledge into end-to-end learning on the molecular graph. The model is strongly validated against seven competitive models and three challenging independent test sets involving pharmaceutical cocrystals, π-π cocrystals and energetic cocrystals, exhibiting superior performance with accuracy higher than 96%, confirming its robustness and generalization. Furthermore, one new energetic cocrystal predicted is successfully synthesized, showcasing high potential of the model in practice. All the data and source codes are available at https://github.com/Saoge123/ccgnet for aiding cocrystal community.

Pharmaceutical cocrystals: A review of preparations, physicochemical properties and applications

Pharmaceutical cocrystals are multicomponent systems in which at least one component is an active pharmaceutical ingredient and the others are pharmaceutically acceptable ingredients. Cocrystallization of a drug substance with a coformer is a promising and emerging approach to improve the performance of pharmaceuticals, such as solubility, dissolution profile, pharmacokinetics and stability. This review article presents a comprehensive overview of pharmaceutical cocrystals, including preparation methods, physicochemical properties, and applications. Furthermore, some examples of drug cocrystals are highlighted to illustrate the effect of crystal structures on the various aspects of active pharmaceutical ingredients, such as physical stability, chemical stability, mechanical properties, optical properties, bioavailability, sustained release and therapeutic effect. This review will provide guidance for more efficient design and manufacture of pharmaceutical cocrystals with desired physicochemical properties and applications.

Selecting a stable solid form of remdesivir using microcrystal electron diffraction and crystal structure prediction

Therapeutic options in response to the coronavirus disease 2019 (COVID-19) outbreak are urgently needed. In this communication, we demonstrate how to support selection of a stable solid form of an antiviral drug remdesivir in quick time using the microcrystal electron diffraction (MicroED) technique and a cloud-based and artificial intelligence implemented crystal structure prediction platform. We present the MicroED structures of remdesivir forms II and IV and conclude that form II is more stable than form IV at ambient temperature in agreement with experimental observations. The combined experimental and theoretical study can serve as a template for formulation scientists in the pharmaceutical industry.

Cocrystal engineering of pharmaceutical solids: therapeutic potential and challenges

This highlight presents an overview of pharmaceutical cocrystal production and its potential in reviving problematic properties of drugs in different dosage forms. The challenges and future outlook of its translational development are discussed.

Virtual coformer screening by a combined machine learning and physics-based approach

Cocrystals as a solid form technology for improving physicochemical properties have gained increasing popularity in the pharmaceutical, nutraceutical, and agrochemical industries.