Efficient Screening for Ternary Molecular Ionic Cocrystals Using a Complementary Mechanosynthesis and Computational Structure Prediction Approach

Porous isoreticular non-metal organic frameworks

Metal-organic frameworks (MOFs) are useful synthetic materials that are built by the programmed assembly of metal nodes and organic linkers1. The success of MOFs results from the isoreticular principle2, which allows families of structurally analogous frameworks to be built in a predictable way. This relies on directional coordinate covalent bonding to define the framework geometry. However, isoreticular strategies do not translate to other common crystalline solids, such as organic salts3-5, in which the intermolecular ionic bonding is less directional. Here we show that chemical knowledge can be combined with computational crystal-structure prediction6 (CSP) to design porous organic ammonium halide salts that contain no metals. The nodes in these salt frameworks are tightly packed ionic clusters that direct the materials to crystallize in specific ways, as demonstrated by the presence of well-defined spikes of low-energy, low-density isoreticular structures on the predicted lattice energy landscapes7,8. These energy landscapes allow us to select combinations of cations and anions that will form thermodynamically stable, porous salt frameworks with channel sizes, functionalities and geometries that can be predicted a priori. Some of these porous salts adsorb molecular guests such as iodine in quantities that exceed those of most MOFs, and this could be useful for applications such as radio-iodine capture9-12. More generally, the synthesis of these salts is scalable, involving simple acid-base neutralization, and the strategy makes it possible to create a family of non-metal organic frameworks that combine high ionic charge density with permanent porosity.

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.

Perchlorate-Free Energetic Oxidizers Enabled by Ionic Cocrystallization

The search for a suitable replacement for the ubiquitous oxidizer ammonium perchlorate (AP) is a top priority to enable more sustainable solid rocket motors. The oxidizing salts ammonium nitrate (AN) and ammonium dinitramide (ADN) are regarded as potential green replacements for AP, but suffer from a plethora of handling and processing issues including poor stability and a needle-like crystal morphology which inhibits dense packing; these prevent their widespread use. In the present work, ionic cocrystallization is leveraged to produce the first cocrystals of these oxidizing salts with an energetic coformer and the first such cocrystals to maintain a positive oxygen balance. The azole-based energetic molecule 5,5'-dinitro-2H,2H'-3,3″-bi-1,2,4-triazole (DNBT) is successfully cocrystallized with AN to yield the cocrystal 2AN:DNBT. Differential scanning calorimetry data confirms that AN, which in its pure form suffers from a problematic solid-state phase transition, is stabilized in the cocrystal. Application of this cocrystallization strategy to ADN produces 2ADN:DNBT, which has the highest oxygen balance of any organic cocrystal.

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.

Overcoming the Hydrophobic Nature of Zinc Phenylacetate Through Co-Crystallization with Isonicotinamide

Zinc phenylacetate (Zn-PA), a substitute for sodium phenylacetate as an ammonia-scavenging drug is hydrophobic, which poses problems for drug dissolution and solubility. We were able to co-crystallize the zinc phenylacetate with isonicotinamide (INAM) and produce a novel crystalline compound (Zn-PA-INAM). The single crystal of this new crystal was obtained, and its structure is reported here for the first time. Zn-PA-INAM was characterized computationally by ab initio, Hirshfeld calculations, CLP-PIXEL lattice energy calculation, and BFDH morphology analysis, and experimentally by PXRD, Sc-XRD, FTIR, DSC, and TGA analyses. Structural and vibrational analyses showed a major modification in intermolecular interaction of Zn-PA-INAM compared to Zn-PA. The dispersion-based pi-stacking in Zn-PA is replaced by coulomb-polarization effect of hydrogen bonds. As a result, Zn-PA-INAM is hydrophilic, improving the wettability and powder dissolution of the target compound in an aqueous solution. Morphology analysis revealed, unlike Zn-PA, Zn-PA-INAM has polar groups exposed on its prominent crystalline faces, reducing the hydrophobicity of the crystal. The shift in average water droplet contact angle from 128.1° (Zn-PA) to 27.1° (Zn-PA-INAM) is strong evidence of a marked decrease in hydrophobicity of the target compound. Finally, HPLC was used to obtain the dissolution profile and solubility of Zn-PA-INAM compared to Zn-PA.

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.

Enhanced NSAIDs Solubility in Drug-Drug Formulations with Ciprofloxacin

Drug-drug salts are a kind of pharmaceutical multicomponent solid in which the two co-existing components are active pharmaceutical ingredients (APIs) in their ionized forms. This novel approach has attracted great interest in the pharmaceutical industry since it not only allows concomitant formulations but also has proved potential to improve the pharmacokinetics of the involved APIs. This is especially interesting for those APIs that have relevant dose-dependent secondary effects, such as non-steroidal anti-inflammatory drugs (NSAIDs). In this work, six multidrug salts involving six different NSAIDs and the antibiotic ciprofloxacin are reported. The novel solids were synthesized using mechanochemical methods and comprehensively characterized in the solid state. Moreover, solubility and stability studies, as well as bacterial inhibition assays, were performed. Our results suggest that our drug-drug formulations enhanced the solubility of NSAIDs without affecting the antibiotic efficacy.

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.

Formation of Salts and Molecular Ionic Cocrystals of Fluoroquinolones and α,ω-Dicarboxylic Acids

The cocrystallization of the fluoroquinolones ciprofloxacin (cip), norfloxacin (nor), and enrofloxacin (enro) with the α,ω-dicarboxylic acids glutaric acid (glu), adipic acid (adi), pimelic acid (pim), suberic acid (sub), azeliac acid (az), and sebacic acid (seb) resulted in 27 new molecular salts and ternary molecular ionic cocrystals of compositions A⁺B^-, A₂ ⁺B^2-, A₂ ⁺B^2-B, and A⁺B^-A. Depending on the solvent, different stoichiomorphs, solvates, or polymorphs were obtained. All salts and cocrystals contain the robust R₂NH₂ ^+...-OOC or R₃NH^+...-OOC synthon but have different supramolecular ring motifs. Moderate solubility enhancements over the parent fluoroquinolones were observed. Salts in the ratio of 1:1 and 2:1 were also prepared by ball-milling. The milled sample nor/az (1:1) was shown to gel the GRAS (generally recognized as safe) solvent propylene glycol, and enro/sub (1:1) was shown to gel both propylene glycol and water. Dynamic rheology measurements confirmed that nor/az and enro/sub behave like viscoelastic materials and supramolecular gels.

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.

Challenges and opportunities of pharmaceutical cocrystals: a focused review on non-steroidal anti-inflammatory drugs

The focus of the review is to discuss the relevant and essential aspects of pharmaceutical cocrystals in both academia and industry with an emphasis on non-steroidal anti-inflammatory drugs (NSAIDs). Although cocrystals have been prepared for a plethora of drugs, NSAID cocrystals are focused due to their humongous application in different fields of medication such as antipyretic, anti-inflammatory, analgesic, antiplatelet, antitumor, and anti-carcinogenic drugs. The highlights of the review are (a) background of cocrystals and other solid forms of an active pharmaceutical ingredient (API) based on the principles of crystal engineering, (b) why cocrystals are an excellent opportunity in the pharma industry, (c) common methods of preparation of cocrystals from the lab scale to bulk quantity, (d) some latest case studies of NSAIDs which have shown better physicochemical properties for example; mechanical properties (tabletability), hydration, solubility, bioavailability, and permeability, and (e) latest guidelines of the US FDA and EMA opening new opportunities and challenges.

Investigation of Poor Solubility of a Salt-Cocrystal Hydrate: A Case Study of the Common-Ion Effect in Betrixaban, an Anticoagulant Drug

Achieving the desired solubility and dissolution of active pharmaceutical ingredients (APIs) continues to be a big challenge in the pharmaceutical industry. In this regard, multicomponent solids of APIs such as salts and cocrystals have shown significant promise in resolving such solubility/dissolution issues. However, very little is known on how the APIs' solubility or dissolution is affected by the drug to coformer ratio in multicomponent solids. Betrixaban, is an anticoagulant drug approved in 2017 for the prevention of venous thromboembolism. During the alternate solid form development studies of the known betrixaban maleate, a rare multicomponent solid form, salt-cocrystal hydrate of betrixaban, was discovered and characterized thoroughly by spectroscopic, thermal, and X-ray crystallographic methods. Significantly, the new betrixaban maleate maleic acid hydrate (1:1:2:1) form has shown lower melting point (80 °C) as compared to its parent salt (197.5 °C). From such a large melting difference (117 °C) between the salt and salt-cocrystal hydrate of API, we anticipated substantially better solubility for the salt-cocrystal hydrate (low enthalpy). Furthermore, the predicted solubility also supported our anticipation. However, the powder dissolution tests at different pH conditions provided contrary results, that is, the salt-cocrystal hydrate showed 10 times lower solubility as compared to its salt. A detailed investigation, considering all the potential factors, revealed that "common-ion effect" could be a critical factor for the low solubility of the salt-cocrystal hydrate in which the API to coformer ratio is 1:3. To the best of our knowledge, this is the first case study on the solubility of pharmaceutical salt-cocrystal hydrates with an emphasis on "common-ion effect" or drug to coformer ratio.

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

One of the most popular strategies of the optimization of drug properties in the pharmaceutical industry appears to be a solid form changing into a cocrystalline form. A number of virtual screening approaches have been previously developed to allow a selection of the most promising cocrystal formers (coformers) for an experimental follow-up. A significant drawback of those methods is related to the lack of accounting for the crystallinity contribution to cocrystal formation. To address this issue, we propose in this study two virtual coformer screening approaches based on a modern cloud-computing crystal structure prediction (CSP) technology at a dispersion-corrected density functional theory (DFT-D) level. The CSP-based methods were for the first time validated on challenging cases of indomethacin and paracetamol cocrystallization, for which the previously developed approaches provided poor predictions. The calculations demonstrated a dramatic improvement of the virtual coformer screening performance relative to the other methods. It is demonstrated that the crystallinity contribution to the formation of paracetamol and indomethacin cocrystals is a dominant one and, therefore, should not be ignored in the virtual screening calculations. Our results encourage a broad utilization of the proposed CSP-based technology in the pharmaceutical industry as the only virtual coformer screening method that directly accounts for the crystallinity contribution.

Efficient Screening for Ternary Molecular Ionic Cocrystals Using a Complementary Mechanosynthesis and Computational Structure Prediction Approach

The discovery of molecular ionic cocrystals (ICCs) of active pharmaceutical ingredients (APIs) widens the opportunities for optimizing the physicochemical properties of APIs whilst facilitating the delivery of multiple therapeutic agents. However, ICCs are often observed serendipitously in crystallization screens and the factors dictating their crystallization are poorly understood. We demonstrate here that mechanochemical ball milling is a versatile technique for the reproducible synthesis of ternary molecular ICCs in less than 30 min of grinding with or without solvent. Computational crystal structure prediction (CSP) calculations have been performed on ternary molecular ICCs for the first time and the observed crystal structures of all the ICCs were correctly predicted. Periodic dispersion-corrected DFT calculations revealed that all the ICCs are thermodynamically stable (mean stabilization energy=-2 kJ mol-1 ) relative to the crystallization of a physical mixture of the binary salt and acid. The results suggest that a combined mechanosynthesis and CSP approach could be used to target the synthesis of higher-order molecular ICCs with functional properties.