Evaluating the Energetic Driving Force for Cocrystal Formation

A New Look at the Mechanism of Cocrystal Formation and Coformers Exchange in Processes Forced by Mechanical and/or Thermal Stimuli - ex situ and in situ Studies of Low-Melting Eutectic Mixtures

Three different devices: ball mill, hot stage melting, and magic angle spinning (MAS) NMR rotor were used for the preparation of ethenzamide (ET) cocrystals with glutaric acid (GLU), ethylmalonic acid (EMA) and maleic acid (MAL) as coformers. In each case, well-defined binary systems (ET:EMA, ET:GLU, ET:MAL) were obtained. The common features of the two solvent free methods of cocrystal formation (grinding, melting) are presented on the basis of arguments obtained by solid state NMR spectroscopy. Thermal analysis (Differential Scanning Calorimetry) proved that the eutectic phase arises over a wide range of molar ratios of components for each of the binary systems. NMR techniques, supported by theoretical calculations, allowed to provide details about the pathway of the reaction mechanism with atomic accuracy. It was found that the formation of ET cocrystals is a complex process that requires five steps. Each step has been recognized and described. Variable temperature 1D and 2D MAS NMR experiments allowed to track physicochemical processes taking place in a molten state. Moreover, it was found that in a multicomponent mixture consisting of all four components, ET, EMA, GLU, and MAL, ET in the molten phase behaves as a specific selector choosing only one partner to form binary cocrystals according to energy preferences. The process of exchange of coformers in binary systems during grinding, melting, and NMR measurements is described. The stabilization energies (Estab ) and molecular electrostatic potential (MEP) maps computed for the cocrystals under discussion and their individual components rationalize the selection rules and explain the relationships between individual species.

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.

Primary Structural Units "D+A-" Ion Pairs Dominating Near-Infrared Photothermal Conversion of Organic Ionic Cocrystals

The specific stacking mode of D/A blocks is often considered to largely determine the physicochemical properties of cocrystals. However, this rule may fail when encountering a large degree of (integer or near-integer) charge transfer situations. Herein, we explore the extensive correlations between the possible smallest structural units, stacking modes, and near-infrared photothermal conversion (NIR-PTC) properties of F4TCNQ-based cocrystals with typical features of integer-charge-transfer. Surprisingly, these cocrystals with distinct stacking modes display analogous D-A interactions, broad red-shift absorption, ultrafast (1-3 ps) relaxation dynamics of excited states, and excellent NIR-PTC properties. This supports that the resulting "D+A-" ion pairs from integer-charge-transfer may serve as the primary structural units beneath the secondary stacking modes to dominate the property of cocrystals. The stacking modes play an important but only secondary role. This work provides new insights into the structure-dynamics-property correlations and modular design of organic cocrystals for PTC and other applications.

Praziquantel Fifty Years on: A Comprehensive Overview of Its Solid State

This review discusses the entire progress made on the anthelmintic drug praziquantel, focusing on the solid state and, therefore, on anhydrous crystalline polymorphs, amorphous forms, and multicomponent systems (i.e., hydrates, solvates, and cocrystals). Despite having been extensively studied over the last 50 years, new polymorphs and the greater part of their cocrystals have only been identified in the past decade. Progress in crystal engineering science (e.g., the use of mechanochemistry as a solid form screening tool and more strategic structure-based methods), along with the development of analytical techniques, including Synchrotron X-ray analyses, spectroscopy, and microscopy, have furthered the identification of unknown crystal structures of the drug. Also, computational modeling has significantly contributed to the prediction and design of new cocrystals by considering structural conformations and interactions energy. Whilst the insights on praziquantel polymorphs discussed in the present review will give a significant contribution to controlling their formation during manufacturing and drug formulation, the detailed multicomponent forms will help in designing and implementing future praziquantel-based functional materials. The latter will hopefully overcome praziquantel's numerous drawbacks and exploit its potential in the field of neglected tropical diseases.

Coformer-Dependent Physical Stability in a Series of Naringenin-Based Coamorphous Materials with Caffeine, Theophylline, and Theobromine

PURPOSE

To investigate the production and physical stability of coamorphous materials (CAM) of naringenin (NAR) and coformers-caffeine, theophylline or theobromine (CAF/THY/THE, respectively). We independently assessed the impact of moisture and temperature on the physical stability of CAMs, and transformation products after destabilization were examined.

METHODS

Neat grinding, liquid assisted grinding and water slurry were selected to prepare multi-component materials with NAR and CAF, THY or THE. The physical stability of CAMs was investigated at 65°C/<10%RH, 21°C/85% RH and 21°C/<10% RH. Differential scanning calorimetry (DSC) and powder X-ray diffraction (PXRD) were employed to monitor for recrystallization during the stability studies. Glass forming ability of amorphous NAR was assessed to understand CAM formation and physical stability.

RESULTS

NAR:THY and NAR:THE CAMs showed physical stability for approximately nine months, under 21°C/<10% RH while NAR:CAF CAMs destabilized in 2.5 weeks. All CAMs recrystallized within a week at 65°C/<10%RH, and the physical stability at 21°C/85% RH was in the order of - NAR:THY > NAR:THE > NAR:CAF. NAR:THY produced 1:1 cocrystal under all storage conditions, while NAR:CAF destabilized to a 1:1 cocrystal at high RH but a physical mixture at high temperature. NAR:THE was found to recrystallize as a physical mixture in all conditions. NAR was found to be strong glass, with moderate kinetic fragility and good glass forming ability.

CONCLUSION

Five naringenin-based multi-component solids were generated in this study: 3 new CAMs, 1 new cocrystal, and 1 previously reported cocrystal. Destabilization of CAMs was found to be exposure specific and coformer dependent.

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.

Three-Dimensional Noncovalent Interaction Network within [NpO2Cl4]2- Coordination Compounds: Influence on Thermochemical and Vibrational Properties

Noncovalent interactions (NCIs) can influence the stability and chemical properties of pentavalent and hexavalent actinyl (AnO2+/2+) compounds. In this work, the impact of NCIs (actinyl-hydrogen and actinyl-cation interactions) on the enthalpy of formation (ΔHf) and vibrational features was evaluated using Np(VI) tetrachloro compounds as the model system. We calculated the ΔHf values of these solid-state compounds through density functional theory+ thermodynamics (DFT+ T) and validated the results against experimental ΔHf values obtained through isothermal acid calorimetry. Three structural descriptors were evaluated to develop predictors for ΔHf, finding a strong link between ΔHf and hydrogen bond energy (EHtotal) for neptunyl-hydrogen interactions and total electrostatic attraction energy (Eelectrostatictotal) for neptunyl-cation interactions. Finally, we used Raman spectroscopy together with bond order analysis to probe Np=O bond perturbation due to NCIs. Our results showed a strong correlation between the degree of NCIs by axial oxygen and red-shifting of Np=O symmetrical stretch (ν1) wavenumbers and quantitatively demonstrated that NCIs can weaken the Np=O bond. These properties were then compared to those of related U(VI) and Np(V) phases to evaluate the effects of subtle differences in the NCIs and overall properties. In general, the outcomes of our study demonstrated the role of NCIs in stabilizing actinyl solid materials, which consequently governs their thermochemical behaviors and vibrational signatures.

Synthesis, Characterization, and Density Functional Theory Investigation of the Solid-State [UO2Cl4(H2O)]2- Complex

A significant number of solid-state [UO2Cl4]2- coordination compounds have been synthesized and structurally characterized. Yet, despite their purposive relative abundance in aqueous solutions, characterization of aquachlorouranium(VI) complexes remain rare. In the current study, a solid-state uranyl aqua chloro complex ((C4H12N2)2[UO2Cl4(H2O)]Cl2) was synthesized using piperazinium as a charge-balancing ligand, and the structure was determined using single-crystal X-ray diffraction. Using periodic density functional theory, the electronic structure of the [UO2Cl4(H2O)]2- complex was compared to [UO2Cl4]2- to uncover the strengthening of the U═O bond in [UO2Cl4(H2O)]2-. Changes in the strength of the U═O bond were validated further with Raman and IR spectroscopy, where uranyl symmetrical (ν1) and asymmetrical (ν3) stretches were blue-shifted compared to the reference [UO2Cl4]2- complex. Furthermore, the formation energy of the solid-state (C4H12N2)2[UO2Cl4(H2O)]Cl2 complex was calculated to be -287.60 ± 1.75 kJ mol-1 using isothermal acid calorimetry. The demonstrated higher stability relative to the related [UO2Cl4]2- complex was related to the relative stoichiometry of the counterions.

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.

Guiding Principles for the Rational Design of Hybrid Materials: Use of DFT Methodology for Evaluating Non-Covalent Interactions in a Uranyl Tetrahalide Model System

Together with the synthesis and experimental characterization of 14 hybrid materials containing [UO2 X4 ]2- (X=Cl- and Br- ) and organic cations, we report on novel methods for determining correlation trends in their formation enthalpy (ΔHf ) and observed vibrational signatures. ΔHf values were analyzed through isothermal acid calorimetry and a Density Functional Theory+Thermodynamics (DFT+T) approach with results showing good agreement between theory and experiment. Three factors (packing efficiency, cation protonation enthalpy, and hydrogen bonding energy [E H , norm total ${{E}_{H,{\rm { norm}}}^{{\rm { total}}}}$ ]) were assessed as descriptors for trends in ΔHf . Results demonstrated a strong correlation betweenE H , norm total ${E_{{\rm{H}},{\rm{norm}}}^{{\rm{total}}} }$ and ΔHf , highlighting the importance of hydrogen bonding networks in determining the relative stability of solid-state hybrid materials. Lastly, we investigate how hydrogen bonding networks affect the vibrational characteristics of uranyl solid-state materials using experimental Raman and IR spectroscopy and theoretical bond orders and find that hydrogen bonding can red-shift U≡O stretching modes. Overall, the tightly integrated experimental and theoretical studies presented here bridge the trends in macroscopic thermodynamic energies and spectroscopic features with molecular-level details of the geometry and electronic structure. This modeling framework forms a basis for exploring 3D hydrogen bonding as a tunable design feature in the pursuit of supramolecular materials by rational design.

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.

New Lidocaine-Based Pharmaceutical Cocrystals: Preparation, Characterization, and Influence of the Racemic vs. Enantiopure Coformer on the Physico-Chemical Properties

This study describes the preparation, characterization, and influence of the enantiopure vs. racemic coformer on the physico-chemical properties of a pharmaceutical cocrystal. For that purpose, two new 1:1 cocrystals, namely lidocaine:dl-menthol and lidocaine:d-menthol, were prepared. The menthol racemate-based cocrystal was evaluated by means of X-ray diffraction, infrared spectroscopy, Raman, thermal analysis, and solubility experiments. The results were exhaustively compared with the first menthol-based pharmaceutical cocrystal, i.e., lidocaine:l-menthol, discovered in our group 12 years ago. Furthermore, the stable lidocaine/dl-menthol phase diagram has been screened, thoroughly evaluated, and compared to the enantiopure phase diagram. Thus, it has been proven that the racemic vs. enantiopure coformer leads to increased solubility and improved dissolution of lidocaine due to the low stable form induced by menthol molecular disorder in the lidocaine:dl-menthol cocrystal. To date, the 1:1 lidocaine:dl-menthol cocrystal is the third menthol-based pharmaceutical cocrystal, after the 1:1 lidocaine:l-menthol and the 1:2 lopinavir:l-menthol cocrystals reported in 2010 and 2022, respectively. Overall, this study shows promising potential for designing new materials with both improved characteristics and functional properties in the fields of pharmaceutical sciences and crystal engineering.

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.

Mechanistic Study on Transformation of Coamorphous Baicalein-Nicotinamide to Its Cocrystal Form

Recently, coamorphization and cocrystal technologies are of particular interest in the pharmaceutical industry due to their ability to improve the solubility/dissolution and bioavailability of poorly water-soluble drugs, while the coamorphous system often tends to convert into the stable crystalline form usually crystalline physical mixture of each component during formulation preparation or storage. In this paper, BCS II drug baicalein (BAI) along with nicotinamide (NIC) were prepared into a single homogeneous coamorphous system with a single transition temperature at 42.5 °C. Interestingly, instead of the physical mixture of crystalline BAI and NIC, coamorphous BAI-NIC would transform to its cocrystal form under stress of temperature and humidity. The transformation rate under isothermal condition was temperature-dependent, since the crystallinity of the cocrystal enhanced as the temperature increased. Further mechanic studies showed the activation energy for the transformation under non-isothermal condition was calculated to be 184.52 kJ/mol. Additionally, water vapor sorption tests with further solid characterizations indicated the transformation was faster under higher humidity condition due to the faster nucleation process of cocrystal BAI-NIC. This research not only discovered the mechanism of transformation from coamorphous BAI-NIC to cocrystal form, but also provided an unusual method for cocrystal preparation from its coamorphous form.

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.

Aryl Boronic Acids in Columnar Stacked Co-crystalline Materials: Key-Factors Governing the Assembly with Quinones

The assembly of aryl boronic acids B with quinones Q into columnar mixed stacked materials, as previously observed in the solid-state, has been here subjected to a detailed theoretical analysis focusing on the properties of the isolated synthons (HOMO-LUMO energies, electron affinity, ionization potential, reorganization hole/electron energies, partial Hirshfeld atomic charges and conformation stabilities) as well as those of the 1 : 1 adducts (Hirshfeld analysis, IRI surfaces, Hirshfeld atomic charges, hydrogen bond and slipped stacked π-π contributions). The overall picture obtained throught this study shows an intricate pattern of interconnected factors contributing to the formation and stability of the B_x Q_y adducts, and it unveils the importance of parameters such as HOMO-LUMO gap, polarization and charge transfer, in addition to the more evident hydrogen bond and slipped-stacked π-π interactions in the formation of 1 : 1 adducts. An explanation has been also given for the presence in some B_x Q_y adducts of the rare anti-anti conformation for the BO-H group with respect to the most studied and common anti-syn conformation.

Unexpected Factors Affecting the Kinetics of Guest Molecule Release from Investigation of Binary Chemical Systems Trapped in a Single Void of Mesoporous Silica Particles

In this work, we present results for loading of well-defined binary systems (cocrystal, solid solution) and untreated materials (physical mixtures) into the voids of MCM-41 mesoporous silica particles employing three different filling methods. The applied techniques belong to the group of "wet methods" (diffusion supported loading - DiSupLo) and "solvent-free methods" (mechanical ball-mill loading - MeLo, thermal solvent free - TSF). As probes for testing the guest1-guest2 interactions inside the MCM-41 pores we employed the benzoic acid (BA), perfluorobenzoic acid (PFBA), and 4-fluorobenzoic acid (4-FBA). The guests intermolecular contacts and phase changes were monitored employing magic angle spinning (MAS) NMR Spectroscopy techniques and powder X-ray diffraction (PXRD). Since mesoporous silica materials are commonly used in drug delivery system research, special attention has been paid to factors affecting guest release kinetics. It has been proven that not only the content and composition of binary systems, but also the loading technique have a strong impact on the rate of guests release. Innovative methods of visualizing differences in release kinetics are presented.

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.

An experimental and theoretical charge density study of theophylline and malonic acid cocrystallization

The pharmaceutical agent theophylline (THEO) is primarily used as a bronchodilator and is commercially available in both tablet and liquid dosage forms. THEO is highly hygroscopic, reducing its stability, overall shelf-life, and therefore usage as a drug. THEO and dicarboxylic acid cocrystals were designed by Trask et al. in an attempt to decrease the hygroscopic behaviour of THEO; cocrystallisation of THEO with malonic acid (MA) did not improve the hygroscopic stability of THEO in simulated atmospheric humidity testing. The current study employed high-resolution X-ray crystallography, and Density Functional Theory (DFT) calculations to examine the electron density distribution (EDD) changes between the cocrystal and its individual components. The EED changes identified the reasons why the THEO:MA cocrystal did not alter the hygroscopic profile of THEO. The cocrystal was equally porous, with atomic packing factors (APF) similar to those of THEO 0.73 vs. 0.71, respectively. The THEO:MA (1) cocrystal structure is held together by an array of interactions; a heterogeneous synthon between the imidazole and a carboxylic fragment stabilising the asymmetric unit, a pyrimidine-imidazole homosynthon, and an aromatic cycle stack between two THEO moieties have been identified, providing 9.7–12.9 kJ mol⁻¹ of stability. These factors did not change the overall relative stability of the cocrystal relative to its individual THEO and MA components, as shown by cocrystal (1) and THEO being equally stable, with calculated lattice energies within 2.5 kJ mol⁻¹ of one other. The hydrogen bond analysis and fragmented atomic charge analysis highlighted that the formation of (1) combined both the EDD of THEO and MA with no net chemical change, suggesting that the reverse reaction — (1) back to THEO and MA — is of equal potential, ultimately producing THEO hydrate formation, in agreement with the work of Trask et al. These results highlight that a review of the EDD change associated with a chemical reaction can aid in understanding cocrystal design. In addition, they indicate that cocrystal design requires further investigation before becoming a reliable process, with particular emphasis on identifying the appropriate balance of synthon engineering, weak interactions, and packing dynamics.

In this study, the 1:1 cocrystal of theophylline and malonic acid originally engineered by Trask undergoes charge density analysis to rationalise the chemical change process seen throughout crystallisation.

Surprising Chemistry of 6-Azidotetrazolo[5,1-a]phthalazine: What a Purported Natural Product Reveals about the Polymorphism of Explosives

6-Azidotetrazolo[5,1-a]phthalazine (ATPH) is a nitrogen-rich compound of surprisingly broad interest. It is purported to be a natural product, yet it is closely related to substances developed as explosives and is highly polymorphic despite having a nearly planar structure with little flexibility. Seven solid forms of ATPH have been characterized by single-crystal X-ray diffraction. The structures show diverse patterns of molecular organization, including both stacked sheets and herringbone packing. In all cases, N···N and C-H···N interactions play key roles in ensuring molecular cohesion. The high polymorphism of ATPH appears to arise in part from the ability of virtually every atom of nitrogen and hydrogen in the molecule to take part in close N···N and C-H···N contacts. As a result, adjacent molecules can adopt many different relative orientations that are energetically similar, thereby generating a polymorphic landscape with an unusually high density of potential structures. This landscape has been explored in detail by the computational prediction of crystal structures. Studying ATPH has provided insights into the field of energetic materials, where access to multiple polymorphs can be used to improve performance and clarify how it depends on molecular packing. In addition, our work with ATPH shows how valuable insights into molecular crystallization, often gleaned from statistical analyses of structural databases, can also come from in-depth empirical and theoretical studies of single compounds that show distinctive behavior.

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.

Quest: structure and properties of BTF–nitrobenzene cocrystals with different ratios of components

Using the methods of quantum chemistry and AAP, the structure of BTF cocrystals with nitrobenzene, 1,2-, 1,3-, 1,4-dinitrobenzene, 1,3,5-trinitrobenzene, and hexanitrobenzene with different ratios of components (1 : 1, 1 : 2, 1 : 3, 2 : 1, 3 : 1) is modeled.

The structure of methanol at 5.09 GPa: the fortuitous formation of a new high-pressure phase

The crystal structure of a new polymorph of methanol (the ε polymorph) has been determined at 5.09 GPa by single crystal X-ray diffraction.

Creating Organic Functional Materials beyond Chemical Bond Synthesis by Organic Cocrystal Engineering

Organic cocrystal engineering refers to two or more organic molecules stoichiometrically combined and held together by noncovalent intermolecular interactions, which differs from standard chemical synthesis involving covalent bond breakage and formation. Organic cocrystals have unique properties and offer a new strategy for creating enhanced organics. First, however, some key questions need to be addressed: How do diverse monomers affect the intermolecular interaction kinetics during cocrystallization? How do the intermolecular forces in cocrystals affect cocrystal functions? In this Perspective, the definition and advantages of organic cocrystal engineering, specifically in the construction of a reliable intermolecular interaction-stacking structure-performance relationship, are outlined. Additionally, recent developments in the field and the questions above are discussed. Finally, a brief conclusion and some hints on likely future developments are provided.

Using Desmotropes, Cocrystals, and Salts to Manipulate Reactivity in Mechanochemical Organic Reactions

Performing reactions in the solid state offers the largely unexplored possibility of influencing reactivity by manipulating the solid form of the starting reactants. In this work, we explore the use of various solid forms of barbituric acid and its effect on reaction paths and kinetics in a Knoevenagel condensation reaction with vanillin. Modifications of barbituric acid included the use of its desmotrope, a cocrystal, and a salt as the starting reactant. Comparing these reactions with the reaction starting from the commercial keto tautomer of barbituric acid, we find that the reaction kinetics could be accelerated or decelerated, together with a change in the reaction mechanism. Exploring solid forms of reactants can be used as general methodology for manipulating mechanochemical reactivity, further highlighting the benefits of conducting reactions in the solid state, because many of the modifications of solids become unavailable upon dissolution.

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.

Noncovalent Bonding Caught in Action: From Amorphous to Cocrystalline Molecular Thin Films

We demonstrate the solvent-free amorphous-to-cocrystalline transformations of nanoscale molecular films. Exposing amorphous films to vapors of a haloarene results in the formation of a cocrystalline coating. This transformation proceeds by gradual strengthening of halogen-bonding interactions as a result of the crystallization process. The gas-solid diffusion mechanism involves formation of an amorphous metastable phase prior to crystallization of the films. In situ optical microscopy shows mass transport during this process, which is confirmed by cross-section analysis of the final structures using focused ion beam milling combined with scanning electron microscopy. Nanomechanical measurements show that the rigidity of the amorphous films influences the crystallization process. This surface transformation results in molecular arrangements that are not readily obtained through other means. Cocrystals grown in solution crystallize in a monoclinic centrosymmetric space group, whereas the on-surface halogen-bonded assembly crystallizes into a noncentrosymmetric material with a bulk second-order nonlinear optical response.

First and Second Dissociation Enthalpies in Bi-Component Crystals Consisting of Maleic Acid and L-Phenylalanine

The energetics of the stepwise dissociation of a A:B₂ bi-component crystal, according to A:B₂(cr) → A:B(cr) + B(cr) and A:B(cr) → A(cr) + B(cr), was investigated using MA:Phe₂ and MA:Phe (MA = maleic acid; Phe = L-phenylalanine) as model systems. The enthalpy changes associated with these sequential processes and with the overall dissociation reaction A:B₂(cr) → A(cr) + 2B(cr) were determined by solution calorimetry. It was found that they are all positive, indicating that there is a lattice enthalpy gain when MA:Phe₂ is formed, either from the individual precursors or by adding Phe to MA:Phe. Single-crystal X-ray diffraction (SCXRD) analysis showed that MA:Phe₂ is best described as a protic salt containing a maleate anion (MA⁻) and two non-equivalent L-phenylalanine units, both linked to MA⁻ by NH···O hydrogen bonds (H-bond): one of these units is protonated (HPhe⁺) and the other zwitterionic (Phe^±). Only MA⁻ and HPhe⁺ molecules are present in the MA:Phe lattice. In this case, however, NH···O and OH···O H-bonds are formed between each MA⁻ unit and two HPhe⁺ molecules. Despite these structural differences, the enthalpy cost for the removal of the zwitterionic Phe^± unit from the MA:Phe₂ lattice to yield MA:Phe is only 0.9 ± 0.4 kJ mol⁻¹ higher than that for the dissociation of MA:Phe, which requires a proton transfer from HPhe⁺ to MA⁻ and the rearrangement of L-phenylalanine to the zwitterionic, Phe^±, form. Finally, a comparison of the dissociation energetics and structures of MA:Phe and of the previously reported glycine maleate (MA:Gly) analogue indicated that parameters, such as the packing coefficient, density, hydrogen bonds formed, or fusion temperature, are not necessarily good descriptors of dissociation enthalpy or lattice enthalpy trends when bi-component crystals with different molecular composition are being compared, even if the stoichiometry is the same.

Interplay of Halogen and Hydrogen Bonding through Co-Crystallization in Pharmacologically Active Dihydropyrimidines: Insights from Crystal Structure and Energy Framework

A solvent-assisted grinding method has been used to prepare co-crystals in substituted dihydropyrimidines (DHPM) that constitutes pharmacologically active compounds. These were characterized using FT-IR, PXRD, and single-crystal X-ray diffraction. In order to explore the possibility of formation of halogen (XB) and hydrogen bonding (HB) synthons in the solid state, co-crystallization attempts of differently substituted DHPM molecules, containing nitro, hydoxy, and chloro substituents, with different co-formers, such as 1,4-diiodo tetrafluorobenzene (1,4 DITFB) and 3-nitrobenzoic acid (3 NBA) were performed. The XB co-crystals (C2aXB, C2bXB, and C2cXB) prefer the formation of C-I⋅⋅⋅O/C-I⋅⋅⋅S XB synthon, whereas the HB co-crystal (C2dHB) is stabilized by N-H⋅⋅⋅O H-bond formation. Hirshfeld surface analysis revealed that the percentage contribution of intermolecular interactions for XB co-crystals prefer equal contribution of XB synthon along with HB synthon. Furthermore, the interaction energy was analyzed using energy frameworks, which suggests that their stability, a combination of electrostatics and dispersion, is enhanced through XB/HB in comparison to the parent DHPMs.

Theoretical research on cage-like furazan-based energetic compounds and its derivatives

In this manuscript, we reported the design and prediction of two furazan-based cage-like molecules and their derivatives using density function theory (DFT). The heat formation and detonation properties were calculated using Hess's law and Kamlet-Jacobs equations with the B3PW91 method. The molecular stability and geometry were analyzed using the M06-2X method, and molecular crystal structures were predicted based on Monte Carlo simulation, while chemical reactive sites were judged using the PBE0 method based on Fukui function. The theoretical calculation result proved that the designed molecules exhibit ideal symmetric cage-like geometry and show superior physicochemical and detonation properties. Compared with traditional energetic materials, the designed molecules display more positive solid heat formation and lower sensitivity. The designed molecules could be considered promising high energy density material candidates with potential synthesis and application value. Two designed molecules display superior detonation performance and ideal completely symmetric cage-like geometry, which were proved theoretically as a promising HEDM candidate. A series of derivatives also exhibited excellent crystal density and physicochemical properties, while with more stable structure.

Mechanosynthesis of Higher-Order Cocrystals: Tuning Order, Functionality and Size in Cocrystal Design*

The ability to rationally design and predictably construct crystalline solids has been the hallmark of crystal engineering research. To date, numerous examples of multicomponent crystals comprising organic molecules have been reported. However, the crystal engineering of cocrystals comprising both organic and inorganic chemical units is still poorly understood and mostly unexplored. Here, we report a new diverse set of higher-order cocrystals (HOCs) based on the structurally versatile-yet largely unexplored-phosph(V/V)azane heterosynthon building block. The novel ternary and quaternary cocrystals reported are held together by synergistic and orthogonal intermolecular interactions. Notably, the HOCs can be readily obtained either via sequential or one-pot mechanochemical methods. Computational modelling methods reveal that the HOCs are thermodynamically driven to form and that their mechanical properties strongly depend on the composition and intermolecular forces in the crystal, offering untapped potential for optimizing material properties.

Towards Computational Screening for New Energetic Molecules: Calculation of Heat of Formation and Determination of Bond Strengths by Local Mode Analysis

The reliable determination of gas-phase and solid-state heats of formation are important considerations in energetic materials research. Herein, the ability of PM7 to calculate the gas-phase heats of formation for CNHO-only and inorganic compounds has been critically evaluated, and for the former, comparisons drawn with isodesmic equations and atom equivalence methods. Routes to obtain solid-state heats of formation for a range of single-component molecular solids, salts, and co-crystals were also evaluated. Finally, local vibrational mode analysis has been used to calculate bond length/force constant curves for seven different chemical bonds occurring in CHNO-containing molecules, which allow for rapid identification of the weakest bond, opening up great potential to rationalise decomposition pathways. Both metrics are important tools in rationalising the design of new energetic materials through computational screening processes.

Detection and Stabilization of a Previously Unknown Two-Dimensional (Pseudo)polymorph using Lateral Nanoconfinement

We report on the detection and stabilization of a previously unknown two-dimensional (2D) pseudopolymorph of an alkoxy isophthalic acid using lateral nanoconfinement. The self-assembled molecular networks formed by the isophthalic acid derivative were studied at the interface between covalently modified graphite and an organic solvent. When self-assembled on graphite with moderate surface coverage of covalently bound aryl groups, a previously unknown metastable pseudopolymorph was detected. This pseudopolymorph, which was presumably "trapped" in between the surface bound aryl groups, underwent a time-dependent phase transition to the stable polymorph typically observed on pristine graphite. The stabilization of the pseudopolymorph was then achieved by using an alternative nanoconfinement strategy, where the domains of the pseudopolymorph could be formed and stabilized by restricting the self-assembly in nanometer-sized shallow compartments produced by STM-based nanolithography carried out on a graphite surface with a high density of covalently bound aryl groups. These experimental results are supported by molecular mechanics and molecular dynamics simulations, which not only provide important insight into the relative stabilities of the different structures, but also shed light onto the mechanism of the formation and stabilization of the pseudopolymorph under nanoscopic lateral confinement.

An industrial perspective on co-crystals: Screening, identification and development of the less utilised solid form in drug discovery and development

Active pharmaceutical ingredients are commonly marketed as a solid form due to ease of transport, storage and administration. In the design of a drug formulation, the selection of the solid form is incredibly important and is traditionally based on what polymorphs, hydrates or salts are available for that compound. Co-crystals, another potential solid form available, are currently not as readily considered as a viable solid form for the development process. Even though co-crystals are gaining an ever-increasing level of interest within the pharmaceutical community, their acceptance and application is still not as standard as other solid forms such as the ubiquitous pharmaceutical salt and stabilised amorphous formulations. Presented in this chapter is information that would allow for a co-crystal screen to be planned and conducted as well as scaled up using solution and mechanochemistry based methods commonly employed in both the literature and industry. Also presented are methods for identifying the formation of a co-crystal using a variety of analytical techniques as well as the importance of confirming the formation of co-crystals from a legal perspective and demonstrating the legal precedent by looking at co-crystalline products already on the market. The benefits of co-crystals have been well established, and presented in this chapter are a selection of examples which best exemplify their potential. The goal of this chapter is to increase the understanding of co-crystals and how they may be successfully exploited in early stage development.

Cocrystals of Praziquantel: Discovery by Network-Based Link Prediction

Cocrystallization has been promoted as an attractive early development tool as it can change the physicochemical properties of a target compound and possibly enable the purification of single enantiomers from racemic compounds. In general, the identification of adequate cocrystallization candidates (or coformers) is troublesome and hampers the exploration of the solid-state landscape. For this reason, several computational tools have been introduced over the last two decades. In this study, cocrystals of Praziquantel (PZQ), an anthelmintic drug used to treat schistosomiasis, are predicted with network-based link prediction and experimentally explored. Single crystals of 12 experimental cocrystal indications were grown and subjected to a structural analysis with single-crystal X-ray diffraction. This case study illustrates the power of the link-prediction approach and its ability to suggest a diverse set of new coformer candidates for a target compound when starting from only a limited number of known cocrystals.

Terahertz spectroscopic characterizations and DFT calculations of indomethacin cocrystals with nicotinamide and saccharin

Co-crystallization is an effective strategy to improve the drug properties such as solubility and stability. However, its thermodynamic backgrounds, especially lattice vibration, haven't been fully understood. In this work, indomethacin (IND) cocrystals formed with nicotinamide (NIC) and saccharin (SAC) are successfully characterized by using terahertz spectroscopy. DFT calculations at PBE-D3 level with and without constrained unit cell are performed to predict the absorption peaks at spectral range. The results suggest that the DFT calculations with constrained unit cell achieve a better agreement with experimental observations. Based on the optimized geometries and calculated phonons, the thermodynamic contributions from lattice vibrations to cocrystal formations are further evaluated. The findings reveal that the vibrational energy plays a comparable role with electronic energy, but has an opposite impact on these two cocrystal formations.

Microcrystal Electron Diffraction Elucidates Water-Specific Polymorphism-Induced Emission Enhancement of Bis-arylacylhydrazone

Aggregation-induced emission (AIE) phenomena have gained intense interest over the last decades because of its importance in solid-state emission. However, the elucidation of a working mechanism is difficult owing to the limited characterization methods on solid-state molecules, further complicated if dynamic structural changes occur. Here, a series of bis-arylacylhydrazones (BAHs) were synthesized, for which their AIE properties are only turned on by the reversible adsorption of water molecules. We used microcrystal electron diffraction (MicroED) to determine the molecular structures of two BAHs directly from bulk powders (without attempting to grow crystals) prepared in the absence or presence of water adsorption. This study reveals the unambiguous characterization of the dependence of crystal packing on the specific cocrystallization with hydrates. The structural analysis demonstrates that water molecules form strong hydrogen bonds with three neighboring BAH-1, resulting in the almost complete planarization and restriction of the intramolecular rotation of the molecule. MicroED plays an important role in providing a decisive clue for the reversible polymorphism changes induced by the adsorption of water molecules, regulating emissive properties.

Cocrystals Based on 4,4’-bipyridine: Influence of Crystal Packing on Melting Point

The reactions of piperonylic acid (HPip) and cinnamic acid (HCinn) with 4,4’-bipyridine (4,4’-bipy) have been assayed using the same synthetic methodology, yielding two binary cocrystals with different acid:4,4’-bipy molar ratios, (HPip)(4,4’-bipy) (1) and (HCinn)2(4,4’-bipy) (2). The melting point (m.p.) of these cocrystals have been measured and a remarkable difference (ΔT ≈ 78 °C) between them was observed. Moreover, the two cocrystals have been characterized by powder X-ray diffraction (PXRD), elemental analysis (EA), FTIR-ATR, 1H NMR spectroscopies, and single-crystal X-ray diffraction. The study of their structural packings via Hirshfeld surface analysis and energy frameworks revealed the important contribution of the π···π and C-H···π interactions to the formation of different structural packing motifs, this being the main reason for the difference of m.p. between them. Moreover, it has been observed that 1 and 2 presented the same packing motifs as the crystal structure of their corresponding carboxylic acids, but 1 and 2 showed lower m.p. than those of the carboxylic acids, which could be related to the lower strength of the acid-pyridine heterosynthons respect to the acid-acid homosynthons in the crystal structures.

Cocrystal Prediction Using Machine Learning Models and Descriptors

Cocrystals are of much interest in industrial application as well as academic research, and screening of suitable coformers for active pharmaceutical ingredients is the most crucial and challenging step in cocrystal development. Recently, machine learning techniques are attracting researchers in many fields including pharmaceutical research such as quantitative structure-activity/property relationship. In this paper, we develop machine learning models to predict cocrystal formation. We extract descriptor values from simplified molecular-input line-entry system (SMILES) of compounds and compare the machine learning models by experiments with our collected data of 1476 instances. As a result, we found that artificial neural network shows great potential as it has the best accuracy, sensitivity, and F1 score. We also found that the model achieved comparable performance with about half of the descriptors chosen by feature selection algorithms. We believe that this will contribute to faster and more accurate cocrystal development.

Mechanochemical reactivity inhibited, prohibited and reversed by liquid additives: examples from crystal-form screens

We demonstrate that liquid additives can exert inhibitive or prohibitive effects on the mechanochemical formation of multi-component molecular crystals, and report that certain additives unexpectedly prompt the dismantling of such solids into physical mixtures of their constituents. Computational methods were employed in an attempt to identify possible reasons for these previously unrecognised effects of liquid additives on mechanochemical transformations.

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.

One class classification as a practical approach for accelerating π-π co-crystal discovery

The implementation of machine learning models has brought major changes in the decision-making process for materials design. One matter of concern for the data-driven approaches is the lack of negative data from unsuccessful synthetic attempts, which might generate inherently imbalanced datasets. We propose the application of the one-class classification methodology as an effective tool for tackling these limitations on the materials design problems. This is a concept of learning based only on a well-defined class without counter examples. An extensive study on the different one-class classification algorithms is performed until the most appropriate workflow is identified for guiding the discovery of emerging materials belonging to a relatively small class, that being the weakly bound polyaromatic hydrocarbon co-crystals. The two-step approach presented in this study first trains the model using all the known molecular combinations that form this class of co-crystals extracted from the Cambridge Structural Database (1722 molecular combinations), followed by scoring possible yet unknown pairs from the ZINC15 database (21 736 possible molecular combinations). Focusing on the highest-ranking pairs predicted to have higher probability of forming co-crystals, materials discovery can be accelerated by reducing the vast molecular space and directing the synthetic efforts of chemists. Further on, using interpretability techniques a more detailed understanding of the molecular properties causing co-crystallization is sought after. The applicability of the current methodology is demonstrated with the discovery of two novel co-crystals, namely pyrene-6H-benzo[c]chromen-6-one (1) and pyrene-9,10-dicyanoanthracene (2).

Machine learning using one class classification on a database of existing co-crystals enables the identification of co-formers which are likely to form stable co-crystals, resulting in the synthesis of two co-crystals of polyaromatic hydrocarbons.