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.

Supramolecular Organization in Salts of Riluzole with Dihydroxybenzoic Acids-The Key Role of the Mutual Arrangement of OH Groups.

Intermolecular interactions, in particular hydrogen bonds, play a key role in crystal engineering. The ability to form hydrogen bonds of various types and strengths causes competition between supramolecular synthons in pharmaceutical multicomponent crystals. In this work, we investigate the influence of positional isomerism on the packing arrangements and the network of hydrogen bonds in multicomponent crystals of the drug riluzole with hydroxyl derivatives of salicylic acid. The supramolecular organization of the riluzole salt containing 2,6-dihydroxybenzoic acid differs from that of the solid forms with 2,4- and 2,5-dihydroxybenzoic acids. Because the second OH group is not at position 6 in the latter crystals, intermolecular charge-assisted hydrogen bonds are formed. According to periodic DFT calculations, the enthalpy of these H-bonds exceeds 30 kJ·mol-1. The positional isomerism appears to have little effect on the enthalpy of the primary supramolecular synthon (65-70 kJ·mol-1), but it does result in the formation of a two-dimensional network of hydrogen bonds and an increase in the overall lattice energy. According to the results of the present study, 2,6-dihydroxybenzoic acid can be treated as a promising counterion for the design of pharmaceutical multicomponent crystals.

Influence of Co-amorphization on the Physical Stability and Dissolution Performance of an Anthelmintic Drug Flubendazole.

In this work, the co-amorphization approach was applied to flubendazole (FluBZ), resulting in the formation of two novel solid forms of FluBZ with l-phenylalanine (Phe) and l-tryptophan (Trp). A variety of physicochemical techniques have been used to describe new systems, including powder X-ray diffraction, thermal methods, infrared spectroscopy, and scanning electron microscopy. Co-amorphization has been shown to suppress crystallization tendency and considerably increase the shelf-life storage of amorphous flubendazole solid across a wide range of relative humidities. The dissolution behavior of the amorphous forms in biorelevant media at pH = 1.6, pH = 6.5, and 37 °C has been studied in terms of Cmax (maximum FluBZ concentration), Tmax (time to attain peak drug concentration), and AUC (concentration area under the curve during dissolution). At pH = 6.5, a continuous supersaturation and the highest AUC value of all examined systems were observed for the FluBZ-Phe (1:1) system. The phase solubility diagrams revealed that the reason for the better dissolution performance of FluBZ-Phe (1:1) at pH = 6.5 is a complexation between the components in a solution. This work highlights the applicability of co-amorphous systems in improving the physical stability and dissolution performance of drug compounds with poor biopharmaceutical characteristics.

Molecular Dynamics Simulation of Association Processes in Aqueous Solutions of Maleate Salts of Drug-like Compounds: The Role of Counterion.

The study of the formation of microstructures during the interaction of a protonated drug-like compound (API) with a maleic acid monoanion sheds light on the assembly processes in an aqueous solution at the molecular level. Molecular dynamics (MD) simulations coupled with density functional theory (DFT) calculations made it possible to find initial hydrogen bonding motifs during the assembly process, leading to the formation of heterodimers and trimers. The process of trimer formation [protonated API-maleic acid monoanion-protonated API] proceeds through the formation of three intermolecular H-bonds by the CO2- group of the maleic acid monoanion in both systems. The total enthalpy/energy of these H-bonds is more than 70 kJ/mol. Thus, the maleic acid monoanion plays a key role in the processes of association in aqueous solution, and the interaction of the maleic acid monoanion with API is more preferable than the interaction of API molecules with each other. DFT computations in the discrete continuum approximation reveal the spectral features of heterodimers and trimers, and the ATR-IR spectra confirmed these findings. MD simulations followed by DFT calculations made it possible to describe the initial stages of the formation of pharmaceutical cocrystals in an aqueous solution.

Comparison of Proton Acceptor and Proton Donor Properties of H2O and H2O2 in Organic Crystals of Drug-like Compounds: Peroxosolvates vs. Crystallohydrates.

Two new peroxosolvates of drug-like compounds were synthesized and studied by a combination of X-ray crystallographic, Raman spectroscopic methods, and periodic DFT computations. The enthalpies of H-bonds formed by hydrogen peroxide (H2O2) as a donor and an acceptor of protons were compared with the enthalpies of analogous H-bonds formed by water (H2O) in isomorphic (isostructural) hydrates. The enthalpies of H-bonds formed by H2O2 as a proton donor turned out to be higher than the values of the corresponding H-bonds formed by H2O. In the case of H2O2 as a proton acceptor in H-bonds, the ratio appeared reversed. The neutral O∙∙∙H-O/O∙∙∙H-N bonds formed by the lone electron pair of the oxygen atom of water were the strongest H-bonds in the considered crystals. In the paper, it was found out that the low-frequency Raman spectra of isomorphous crystalline hydrate and peroxosolvate of N-(5-Nitro-2-furfurylidene)-1-aminohydantoin are similar. As for the isostructural hydrate and peroxosolvate of the salt of protonated 2-amino-nicotinic acid and maleic acid monoanion, the Raman spectra are different.

A Combined Experimental and Theoretical Study of Nitrofuran Antibiotics: Crystal Structures, DFT Computations, Sublimation and Solution Thermodynamics.

Single crystal of furazolidone (FZL) has been successfully obtained, and its crystal structure has been determined. Common and distinctive features of furazolidone and nitrofurantoin (NFT) crystal packing have been discussed. Combined use of QTAIMC and Hirshfeld surface analysis allowed characterizing the non-covalent interactions in both crystals. Thermophysical characteristics and decomposition of NFT and FZL have been studied by differential scanning calorimetry (DSC), thermogravimetric analysis (TG) and mass-spectrometry. The saturated vapor pressures of the compounds have been measured using the transpiration method, and the standard thermodynamic functions of sublimation were calculated. It was revealed that the sublimation enthalpy and Gibbs energy of NFT are both higher than those for FZL, but a gain in the crystal lattice energy of NFT is leveled by an entropy increase. The solubility processes of the studied compounds in buffer solutions with pH 2.0, 7.4 and in 1-octanol was investigated at four temperatures from 298.15 to 313.15 K by the saturation shake-flask method. The thermodynamic functions of the dissolution and solvation processes of the studied compounds have been calculated based on the experimental data. Due to the fact that NFT is unstable in buffer solutions and undergoes a solution-mediated transformation from an anhydrate form to monohydrate in the solid state, the thermophysical characteristics and dissolution thermodynamics of the monohydrate were also investigated. It was demonstrated that a combination of experimental and theoretical methods allows performing an in-depth study of the relationships between the molecular and crystal structure and pharmaceutically relevant properties of nitrofuran antibiotics.

Polymorphic forms of antiandrogenic drug nilutamide: structural and thermodynamic aspects.

Attempts to obtain new cocrystals of nonsteroidal antiandrogenic drug nilutamide produced alternative polymorphic forms of the compound (Form II and Form III) and their crystal structures were elucidated by single-crystal X-ray diffraction. Apart from the cocrystallization technique, lyophilization was found to be an effective strategy for achieving polymorph control of nilutamide, which was difficult to obtain by other methods. The physicochemical properties and relative stability of the commercial Form I and newly obtained Form II were comprehensively investigated by a variety of analytical methods (thermal analysis, solution calorimetry, solubility, and sublimation), whereas for Form III, only a handful of experimental parameters were obtained due to the elusive nature of the polymorph. Form I and Form II were found to be monotropically related, with Form I being confirmed as the thermodynamically most stable solid phase. In addition, the performance of different DFT-D and semi-empirical schemes for lattice energy calculation and polymorph energy ranking was compared and analysed. Lattice energy calculations using periodic DFT at B3LYP-D3/6-31(F+)G(d,p) and PBEh-3c/def2-mSVP levels of theory were found to provide the most accurate lattice energy values for Form I against experimental data, while PIXEL and PBEh-3c/def2-mSVP were the only methods that predicted the correct order of stability of Forms I and II.

Novel cocrystals of itraconazole: Insights from phase diagrams, formation thermodynamics and solubility.

In this work, the cocrystallization approach was applied to itraconazole (ITR), a very slightly soluble triazole antifungal drug, which led to the formation of two new solid forms of ITR with 4-aminobenzoic acid (4AmBA) and 4-hydroxybenzamide (4OHBZA). A thermodynamic analysis of the solid-liquid binary phase diagrams for the (ITR + 4AmBA) and (ITR + 4OHBZA) systems provided conclusive evidence of the cocrystal stoichiometry: 1:1 for the cocrystal with 4-aminobenzoic acid, and 1:2 for the cocrystal with 4-hydroxybenzamide. Powder X-Ray diffraction analysis confirmed the formation of two different polymorphic forms of the [ITR + 4OHBZA] (1:2) cocrystal obtained either through solution or melt crystallization. Cocrystal formation and polymorphic transition processes were investigated in detail by the DSC and HSM methods. The thermodynamic functions of cocrystal formation were estimated from the solubility of the cocrystals and the corresponding solubility of the pure compounds at different temperatures. The combination of ITR and 4OHBZA was found to be more favorable than the reaction between ITR and 4AmBA in terms of both Gibbs energy and enthalpy. The pH-solubility behavior of the cocrystals was investigated at different pH values using eutectic concentrations of the components and the cocrystal solubility advantage was estimated. It was found that the cocrystallization of itraconazole with 4OHBZA and 4AmBA can potentially increase the drug solubility at pH1.2 and 37 °C by 225 and 64 times, respectively. The cocrystal dissolution behavior in biorelevant media was analyzed in terms of Cmax, σmax parameters (the maximum ITR concentration and supersaturation), and AUC (the concentration area under the curve during the dissolution - supersaturation - precipitation process). The cocrystals had similar σmax values during the dissolution and sustained supersaturation for up to 6 h, which gave them an advantage in the AUC values (13-37 times higher) over the drug. The differences in the dissolution profiles of the cocrystals were rationalized in terms of their dissolution rate values.

Combined X-ray Crystallographic, IR/Raman Spectroscopic, and Periodic DFT Investigations of New Multicomponent Crystalline Forms of Anthelmintic Drugs: A Case Study of Carbendazim Maleate.

Synthesis of multicomponent solid forms is an important method of modifying and fine-tuning the most critical physicochemical properties of drug compounds. The design of new multicomponent pharmaceutical materials requires reliable information about the supramolecular arrangement of molecules and detailed description of the intermolecular interactions in the crystal structure. It implies the use of a combination of different experimental and theoretical investigation methods. Organic salts present new challenges for those who develop theoretical approaches describing the structure, spectral properties, and lattice energy Elatt. These crystals consist of closed-shell organic ions interacting through relatively strong hydrogen bonds, which leads to Elatt > 200 kJ/mol. Some technical problems that a user of periodic (solid-state) density functional theory (DFT) programs encounters when calculating the properties of these crystals still remain unsolved, for example, the influence of cell parameter optimization on the Elatt value, wave numbers, relative intensity of Raman-active vibrations in the low-frequency region, etc. In this work, various properties of a new two-component carbendazim maleate crystal were experimentally investigated, and the applicability of different DFT functionals and empirical Grimme corrections to the description of the obtained structural and spectroscopic properties was tested. Based on this, practical recommendations were developed for further theoretical studies of multicomponent organic pharmaceutical crystals.

Interplay of π-stacking and inter-stacking interactions in two-component crystals of neutral closed-shell aromatic compounds: periodic DFT study.

This paper bridges the gap between high-level ab initio computations of gas-phase models of 1 : 1 arene-arene complexes and calculations of the two-component (binary) organic crystals using atom-atom potentials. The studied crystals consist of electron-rich and electron-deficient compounds, which form infinite stacks (columns) of heterodimers. The sublimation enthalpy of crystals has been evaluated by DFT periodic calculations, while intermolecular interactions have been characterized by Bader analysis of the periodic electronic density. The consideration of aromatic compounds without a dipole moment makes it possible to reveal the contribution of quadrupole-quadrupole interactions to the π-stacking energy. These interactions are significant for heterodimers formed by arenes with more than 2 rings, with absolute values of the traceless quadrupole moment (Q zz) larger than 10 D Å. The further aggregation of neighboring stacks is due to the C-H⋯F interactions in arene/perfluoroarene crystals. In crystals consisting of arene and an electron-deficient compound such as pyromellitic dianhydride, aggregation occurs due to the C-H⋯O interactions. The C-H⋯F and C-H⋯O inter-stacking interactions make the main contribution to the sublimation enthalpy, which exceeds 150 kJ mol-1 for the two-component crystals formed by arenes with more than 2 rings.

Solving the enigma of weak fluorine contacts in the solid state: a periodic DFT study of fluorinated organic crystals.

The nature and strength of weak interactions with organic fluorine in the solid state are revealed by periodic density functional theory (periodic DFT) calculations coupled with experimental data on the structure and sublimation thermodynamics of crystalline organofluorine compounds. To minimize other intermolecular interactions, several sets of crystals of perfluorinated and partially fluorinated organic molecules are considered. This allows us to establish the theoretical levels providing an adequate description of the metric and electron-density parameters of the C-F⋯F-C interactions and the sublimation enthalpy of crystalline perfluorinated compounds. A detailed comparison of the C-F⋯F-C and C-H⋯F-C interactions is performed using the relaxed molecular geometry in the studied crystals. The change in the crystalline packing of aromatic compounds during their partial fluorination points to the structure-directing role of C-H⋯F-C interactions due to the dominant electrostatic contribution to these contacts. C-H⋯F-C and C-H⋯O interactions are found to be identical in nature and comparable in energy. The factors that determine the contribution of these interactions to the crystal packing are revealed. The reliability of the results is confirmed by considering the superposition of the electrostatic potential and electron density gradient fields in the area of the investigated intermolecular interactions.

Pharmaceutical salts of emoxypine with dicarboxylic acids.

New salt forms of the antioxidant drug emoxypine (EMX, 2-ethyl-6-methylpyridin-3-ol) with pharmaceutically acceptable maleic (Mlt), malonic (Mln) and adipic (Adp) acids were obtained {emoxypinium maleate, C8H12NO+·C4H3O4-, [EMX+Mlt], emoxypinium malonate, C8H12NO+·C3H3O4-, [EMX+Mln], and emoxypinium adipate, C8H12NO+·C6H9O4-, [EMX+Adp]} and their crystal structures determined. The molecular packing in the three EMX salts was studied by means of solid-state density functional theory (DFT), followed by QTAIMC (quantum theory of atoms in molecules and crystals) analysis. It was found that the major contribution to the packing energy comes from pyridine-carboxylate and hydroxy-carboxylate heterosynthons forming infinite one-dimensional ribbons, with [EMX+Adp] additionally stabilized by hydrogen-bonded C(9) chains of Adp- ions. The melting processes of the [EMX+Mlt] (1:1), [EMX+Mln] (1:1) and [EMX+Adp] (1:1) salts were studied and the fusion enthalpy was found to increase with the increase of the calculated lattice energy. The dissolution process of the EMX salts in buffer (pH 7.4) was also studied. It was found that the formation of binary crystals of EMX with dicarboxylic acids increases the EMX solubility by more than 30 times compared to its pure form.

Influence of Secondary Interactions on the Structure, Sublimation Thermodynamics, and Solubility of Salicylate:4-Hydroxybenzamide Cocrystals. Combined Experimental and Theoretical Study.

Cocrystal screening of 4-hydroxybenzamide with a number of salicylates (salicylic acid, SA; 4-aminosalicylic acid, PASA; acetylsalicylic acid, ASA; and salicylsalicylic acid, SSA) was conducted to confirm the formation of two cocrystals, [SA+4-OHBZA] (1:1) and [PASA+4-OHBZA] (1:1). Their structures were determined using single-crystal X-ray diffraction, and the hydrogen-bond network topology was studied. Thermodynamic characteristics of salicylic acid cocrystal sublimation were obtained experimentally. It was proved that PASA cocrystallization with 4-OHBZA makes the drug more stable and prevents the irreversible process of decarboxylation of PASA resulting in formation of toxic 3-aminophenol. The pattern of non-covalent interactions in the cocrystals is described quantitatively using solid-state density functional theory followed by Bader analysis of the periodic electron density. It has been found that the total energy of secondary interactions between synthon atoms and the side hydroxyl group of the acid molecule in [SA+4-OHBZA] (1:1) and [PASA+4-OHBZA] (1:1) cocrystals is comparable to the energy of the primary acid-amide heterosynthon. The theoretical value of the sublimation enthalpy of [SA+4-OHBZA], 231 kJ/mol, agrees fairly well with the experimental one, 272 kJ/mol. The dissolution experiments with [SA+4-OHBZA] have proved that the relatively large cocrystal stability in relation to the stability of its components has a negative effect on the dissolution rate and equilibrium solubility. The [PASA+4-OHBZA] (1:1) cocrystal showed an enhancement of apparent solubility compared to that of the corresponding pure active pharmaceutical ingredient, while their intrinsic dissolution rates are comparable.

Pharmaceutical salts of ciprofloxacin with dicarboxylic acids.

New salts of antibiotic drug ciprofloxacin (CIP) with pharmaceutically acceptable maleic (Mlt), fumaric (Fum) and adipic (Adp) acids were obtained and their crystal structures were determined. The crystal lattices of the fumarate and adipate salts were found to accommodate the water molecules, while the maleate salt was anhydrous. The dehydration and melting processes were analyzed by means of differential scanning calorimetry and thermogravimetric analysis. Solubility and intrinsic dissolution rates of the salts were measured in pharmaceutically relevant buffer solutions with pH 1.2 and pH 6.8. Under acidic conditions, the salts were found to be less soluble than the parent form of drug, while the [CIP+Fum+H2O] and [CIP+Mlt] solids showed enhanced dissolution rate when compared to a commercially available ciprofloxacin hydrochloride hydrate. In the pH 6.8 solution, all the salts demonstrated solubility improvement and faster dissolution rate with respect to pure CIP.

Cocrystal screening of hydroxybenzamides with benzoic acid derivatives: a comparative study of thermal and solution-based methods.

The main problem occurring at the early stages of cocrystal search is the choice of an effective screening technique. Among the most popular techniques of obtaining cocrystals are crystallization from solution, crystallization from melt and solvent-drop grinding. This paper represents a comparative analysis of the following screening techniques: DSC cocrystal screening method, thermal microscopy and saturation temperature method. The efficiency of different techniques of cocrystal screening was checked in 18 systems. Benzamide and benzoic acid derivatives were chosen as model systems due to their ability to form acid-amide supramolecular heterosynthon. The screening has confirmed the formation of 6 new cocrystals. The screening by the saturation temperature method has the highest screen-out rate but the smallest range of application. DSC screening has a satisfactory accuracy and allows screening over a short time. Thermal microscopy is most efficient as an additional technique used to interpret ambiguous DSC screening results. The study also included an analysis of the influence of solvent type and component solubility on cocrystal formation.

Pharmaceutical cocrystals of diflunisal and diclofenac with theophylline.

Pharmaceutical cocrystals of nonsteroidal anti-inflammatory drugs diflunisal (DIF) and diclofenac (DIC) with theophylline (THP) were obtained, and their crystal structures were determined. In both of the crystal structures, molecules form a hydrogen bonded supramolecular unit consisting of a centrosymmetric dimer of THP and two molecules of active pharmaceutical ingredient (API). Crystal lattice energy calculations showed that the packing energy gain of the [DIC + THP] cocrystal is derived mainly from the dispersion energy, which dominates the structures of the cocrystals. The enthalpies of cocrystal formation were estimated by solution calorimetry, and their thermal stability was studied by differential scanning calorimetry. The cocrystals showed an enhancement of apparent solubility compared to the corresponding pure APIs, while the intrinsic dissolution rates are comparable. Both cocrystals demonstrated physical stability upon storing at different relative humidity.

Salicylamide cocrystals: screening, crystal structure, sublimation thermodynamics, dissolution, and solid-state DFT calculations.

A new cocrystal of 2-hydroxybenzamide (A) with 4-acetamidobenzoic acid (B) has been obtained by the DSC screening method. Thermophysical analysis of the aggregate [A:B] has been conducted and a fusion diagram has been plotted. Cocrystal formation from melts was studied by using thermomicroscopy. A cocrystal single-crystal was grown and its crystal structure was determined. The pattern of noncovalent interactions has been quantified using the solid-state DFT computations coupled with the Bader analysis of the periodic electron density. The sublimation processes of A-B cocrystal have been studied and its thermodynamic functions have been calculated. The classical method of substance transfer by inert gas-carrier was chosen to investigate sublimation processes experimentally. The lattice energy is found to be 143 ± 4 kJ/mol. It is lower than the sum of the corresponding values of the cocrystal pure components. The theoretical value of the lattice energy, 156 kJ/mol, is in reasonable agreement with the experimental one. A ternary phase diagram of solubility (A-B-ethanol) has been plotted and the areas with solutions for growing thermodynamically stable cocrystals have been determined.