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

Abstract

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⁻¹. The positional isomerism appears to have little effect on the enthalpy of the primary supramolecular synthon (65–70 kJ·mol⁻¹), 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.

1. Introduction

The design of pharmaceutical multicomponent solids, primarily relies on the utilization of supramolecular synthons, as reliable building blocks that enable the construction of molecular assemblies with desired structural features and properties in a predictable manner [1,2,3]. Even though various types of intermolecular interactions concurrently exist in the crystal structure and contribute to its thermodynamic stability [4,5,6,7], intermolecular hydrogen bonds (H-bonds) play a more prominent role than others, due to their strength and directionality [8], allowing crystal engineering principals to be used [9,10,11,12]. However, the possibility to properly regulate or anticipate the formation of multiple supramolecular synthons in multicomponent crystals, consisting of molecules with a diverse spectrum of functional groups, remains a challenge, due to a lack of understanding of the nature and strength of intermolecular interactions responsible for competition between different synthons in organic solids [13,14,15]. Therefore, it is of great importance to provide a quantitative evaluation of the binding strength and nature of synthons, to rationalize their stabilizing and structure-directing roles in the formation of multicomponent crystals [16,17,18,19,20,21].

Riluzole (2-amino-6-trifluoromethoxy-benzothiazol) is one of the few approved options to treat amyotrophic lateral sclerosis [22]. Treatment with riluzole (RLZ) leads to slower degradation of neurons and extended survival without tracheostomy [23,24]. The mechanism of action of riluzole is complex, involving glutamatergic neurotransmission enhancement [25], inhibition of fast and persistent sodium currents [26,27], neuroprotective activity [28], as well as neuroplasticity stimulation via stimulation of the neurotrophic factors [29,30]. RLZ was shown to exhibit neuroprotective, anticonvulsant, and antidepressant activity in multiple in vivo studies [31]. In addition, RLZ can be used for the treatment of severe mood, anxiety and impulsive disorders, including obsessive-compulsive disorder, unipolar and bipolar depression, and generalized anxiety disorder [32,33]. Recent studies have demonstrated the potential of repurposing RLZ as an antiproliferative agent for many types of cancer, including skin, breast, pancreas, colon, liver, bone, brain, lung, and nasopharynx cancers ([34] and references therein). Because of its ability to cross the blood–brain barrier, RLZ is one of the few promising possibilities for glioma therapy [35,36] and brain metastasis in late-stage malignant melanoma [37]. The low systemic bioavailability of RLZ (60% [38,39]) is limited by its poor aqueous solubility (350 mg/L at 37 °C and pH 7 [40]), which is typical for a BCS Class II drug [41]. The potential use of RLZ in psychiatry and oncology stimulates the search for novel forms of this drug, with improved dissolution performance and reduced side effects. During its 20 years of use, many approaches for optimizing the drug delivery of RLZ have been proposed, including salt/cocrystal formation [40,42,43], encapsulation in nanoparticles [44], amorphous form stabilization [45], the use of oral suspensions [46], biodegradable sublingual foils [47], and microemulsions for intranasal delivery [48].

From the crystal engineering point of view, RLZ is known to form multicomponent crystals, with amino and carboxylic acids, aromatic amides, etc. [40,42,43]. Dihydroxybenzoic acids are promising coformers/counterions of RLZ due to their safety and biological activity. 3,4-Dihydroxybenzoic acid (syringic acid) exhibits antibacterial, antifungal, and antiprotozoal activity, 2,5-dihydroxybenzoic acid (gentisic acid) displays antioxidant properties [49], while 2,6-dihydroxybenzoic acid and other derivatives of salicylic acid display anti-inflammatory activity [50] (Scheme 1).

Based on a survey of multicomponent crystals of RLZ, deposited in the Cambridge Structural Database (CSD) [51] (Supplementary Materials, Section S1), we concluded that in multicomponent crystals with salicylic acid and its hydroxyl derivatives, RLZ exists in the protonated state (as a cation), whereas acids exist in the anionic form. Depending on the position of the OH group in hydroxyl derivatives of salicylic acid, either intra- or intermolecular H-bonds can form. Because RLZ can only generate three H-bonds in the protonated state, the emergence of a new proton donor, with an excess of unshared electron pairs, may cause a dramatic shift in the H-bond network in a multicomponent crystal [52]. In addition, a change in the number of intra- and intermolecular H-bonds can result in a significant change in the lattice energy value [53]. Thus, the position of the second OH group in the salicylic acid derivative is likely to play an essential structure-directing role in riluzole salts. In other words, hydroxyl derivatives of salicylic acid are envisaged as counterions for the design of pharmaceutical salts of riluzole.

This work aims to study the effect of positional isomerism in hydroxyl derivatives of salicylic acid on the H-bond network and supramolecular organization in RLZ salts. Collaborative research, using X-ray diffraction and periodic (solid-state) density functional theory (DFT), was conducted to highlight the critical significance of the mutual arrangement of OH groups in the packing of these multicomponent crystals.

2. Materials and Methods

2.1. Experimental

Riluzole of 98% purity was obtained from TCI Chemicals. All the benzoic acid derivatives were purchased from Sigma-Aldrich and were used as received. The purity of materials was checked by differential scanning calorimetry (DSC). The solvents were of analytical or chromatographic grade.

2.1.1. Mechanochemical Experiments

The grinding experiments were performed using a Fritsch planetary micro mill, model Pulverisette 7, in 12 mL agate grinding jars with ten 5 mm agate balls at a rate of 500 rpm for 60 min. In a typical experiment, 50–60 mg of a physical mixture, with a 1:1 molar ratio, was placed into a grinding jar, and 50 μL of solvent (water, methanol, ethanol, acetonitrile, or ethyl acetate) was added with a micropipette. The resulting powder samples were analyzed by powder X-ray diffraction (PXRD), DSC, TGA, and FTIR (see Supplementary Materials, Section S2).

2.1.2. Single Crystal XRD

Single crystal diffraction data were collected using a Bruker D8 Venture diffractometer (Bruker AXS, Karlsruhe, Germany) with graphite-monochromated Mo-Kα radiation (λ = 0.71073 Å). It was found that crystals [RLZ+2,4DHBA] degrade below 170 K and thus reflection intensities were collected at 175 K. Absorption corrections, based on measurements of equivalent reflections, were applied [54]. The structures were solved by direct methods and refined by full matrix least-squares on F² with anisotropic thermal parameters for all non-hydrogen atoms [55]. In the structure [RLZ+2,6DHBA], the trifluoromethoxy group is disordered over two positions, with occupancy ratio 0.785(3)/0.225(3), and all fluorine atoms were refined with restrained C–F and F···F distances (SADI instruction). All hydrogen atoms were found from difference Fourier synthesis and refined isotropically. Experimental details are listed in Table S2. The crystallographic data were deposited with the Cambridge Crystallographic Data Centre as supplementary publications under the CCDC numbers 2238373, 2238374. This information can be obtained free of charge from the Cambridge Crystallographic Data Centre via www.ccdc.cam.ac.uk/data_request/cif (accessed on 21 February 2023).

2.2. Computational Details

2.2.1. Periodic DFT Calculations

B3LYP [56,57] and PBE-D3 [58,59] are the most popular functionals in periodic (solid-state) DFT calculations of molecular crystals [60,61]. B3LYP and PBE-D3 provide a grounded trade-off between the accuracy and the rate of calculations of experimentally observed properties of molecular crystals [62,63,64,65,66,67,68]. Computations were conducted using the CRYSTAL17 package [69]. The 6-31G** basis set was used. The space groups and the unit cell parameters of the crystal obtained from the X-ray diffraction experiment were fixed, and the structural relaxations were limited to the positional parameters of the atoms [70,71]. Tolerance on energy controlling the self-consistent field convergence for geometry optimizations and frequency computations was set to 10⁻¹⁰ and 10⁻¹¹ Hartree, respectively. The shrinking factor of the reciprocal space net was set to three. The structures optimized at the PBE-D3/6-31G** level were found to correspond to the potential energy surface’s minimum point. Therefore, PBE-D3 was used in estimating the lattice energy using Equation (3).

2.2.2. Enthalpy of the Intermolecular H-Bonds

The enthalpy ΔH_HB of the intermolecular H-bond in the considered crystals was estimated using the Rozenberg approach [72]:

−ΔH_HB [kJ·mol⁻¹] = 0.134·R(H∙∙∙O)^−3.05,(1)

where the R(H∙∙∙O) is the H∙∙∙O distance (nm). Equation (1) describes −ΔH_HB varying from 10 to 80 kJ·mol⁻¹ [72]. The R(H∙∙∙O) values were obtained as a result of geometry optimization.

2.2.3. Estimation of the Crystal Lattice Energy

The energy of intermolecular noncovalent interaction E_int in the considered crystal was evaluated according to ref. [73] as follows:

E_int [kJ·mol⁻¹] = 1124·G_b [atomic units](2)

where G_b is the positively defined local electronic kinetic energy density at the (3,−1) bond critical point, corresponding to the considered intermolecular noncovalent interaction [74]. The methodology of the calculation is presented elsewhere [75]. Periodic wave functions computed at the PBE-D3/6-31G** level were used in Topond14 [76]. Equation (2) overestimates energies for the short (strong) intermolecular H-bonds [77]. Therefore, the use of Equation (1) is preferable for systems with short H-bonds, both in crystals [77] and in solutions [78].

On the other hand, Equation (2) enables one to estimate the lattice energy of a two-component organic crystal using the following approach [75]:

E_latt = ∑_i∑_j<i(E_int,j,i),(3)

where indices j and i denote the atoms belonging to different molecules. Equation (3) is BSSE-free. Details of the E_latt computations in the two-component crystals using Equation (3) are given elsewhere [75]. In the present study, we were interested in the relative E_latt values estimated for riluzole salts. Hence, overestimating the values of E_int given by Equation (2) does not play a significant role.

2.2.4. Molecular Electrostatic Potential

The molecular electrostatic potential surfaces (MEPs) [79] of the RLZ and selected (di)hydroxybenzoic acids were generated at the B3LYP/def2-tzvp level of theory, with the latest version of the D3 dispersion correction, using Gaussian09 [80]. Calculations were carried out using the optimized molecular geometries. The local maxima and minima sites on the molecular electrostatic potential surfaces were extracted using the Multiwfn software [81].

3. Results

In the present work, two crystals of riluzole, namely, RLZ+2,6-dihydroxybenzoic acid [RLZ+2,6DHBA] (1:1) and RLZ+2,4-dihydroxybenzoic acid [RLZ+2,4DHBA] (1:1), were obtained and structurally characterized (see Table S2, Figures S1, S2, S4, S6 and S7). The details of the screening experiments and phase characterization using powder XRD, DSC, thermogravimetric analysis (TGA), and Fourier-transform infrared spectroscopy (FTIR) are given as Supplementary Materials (Sections S2 and S3). The formation of the novel phase was also observed in the RLZ+2,3-dihydroxybenzoic acid system (Figure S3), which was also proved to be a salt, based on its FTIR spectrum (Figure S4) and DSC/TGA profile (Figure S5). However, we were unable to grow single crystals of sufficient size for structural characterization. For comparison, crystalline RLZ+2,5-dihydroxybenzoic acid [RLZ+2,5DHBA] (1:1) and RLZ+salicylic acid [RLZ+2HBA] (1:1), obtained earlier [43], were also considered in the present study.

All the crystals under consideration contain a robust amino-N-heterocycle carboxylate supramolecular heterosynthon, of $R_2^2\left(8\right)$ topology. This synthon is structurally identical to 2-aminopyridinium carboxylate [82], which is known to play a structure-directing role in various multicomponent crystals [83,84,85,86,87,88]. The mutual orientation of molecules in the synthon is determined by the side H-bonds, and varies from coplanar in [RLZ+2,5DHBA] and [RLZ+2,6DHBA] to the twisted one in [RLZ+2HBA] and [RLZ+2,4DHBA]. The addition of a second hydroxyl group in the salicylic acid changes the H-bonding pattern considerably, in a non-systematic manner, see Figure 1 and Figure 2 and Table 1. The second hydrogen of the amino group of RLZ, denoted as N_a-H_c, forms in these salts an H-bond with the lone electron pair of the OH group of the neighboring acid molecule, rather than with the oxygen of the CO₂− group, c.f. Scheme 4 in Ref. [73], Figure 1 and Figure 2. The supramolecular organization in crystalline [RLZ+2,6DHBA] (1:1) differs from that in [RLZ+2,4DHBA] (1:1) and [RLZ+2,5DHBA] (1:1). Because the second OH group is absent in the sixth position, intermolecular charge-assisted hydrogen bonds are formed in the latter crystals (Table 1). This isomerism leads to the appearance of a two-dimensional network of H-bonds, c.f. Figure 1 and Figure 2. The CO₂− group of the counterion in all RLZ salts, forms four H-bonds. However, the resulting H-bond network is determined by the number of intramolecular H-bonds. [RLZ+2,6DHBA] (1:1) is the only crystal with two such H-bonds. As a consequence, 2,6-dihydroxybenzoic acid can form only three “unique” intermolecular H-bonds, while other OH derivatives of salicylic acid form four intermolecular H-bonds.

As can be seen from the crystal structure data, different positions of the hydroxyl group in the counterion lead to variations in the H-bond topology and mutual orientation of H-bonded units, without significant changes in the interaction energy (Table 1). A structural similarity analysis, performed with the XPac program [89,90], revealed that the salts of RLZ with salicylic and isomeric dihydroxybenzoic acids have no common packing motifs that could be associated with any kind of non-covalent interaction. Such a variety of crystal structures of structurally related compounds, based on a single “basic” supramolecular heterosynthon, highlights the gap between the topology of H-bonds and the actual packing architecture. The differences in the packing of heterodimers observed in the [RLZ+2,4DHBA] salt, and its solvate with methanol [42], are explained by the participation of the solvent molecule in the H-bond network. No similarity was found between [RLZ+2,6DHBA] and the reported [RLZ+2HBA] crystal [43], even taking into account the fact that the number and topology of H-bonds in these crystals are very similar. In the latter structure, the RLZ-acid dimer is twisted from planarity, which greatly affects the packing. The non-planar structure of the dimer allows the keto oxygen atom to participate in a side N–H···O H-bond with the amide group of RLZ. For [RLZ+2,6DHBA], this site is inaccessible for intermolecular contacts due to the strong O–H···O intramolecular bond.

The studied crystals share the same T-shaped mutual orientation of heterodimers in the hydrogen-bonded chain, despite having varying symmetry and packing. In contrast to the two-centered $R_2^2\left(8\right)$ motif, that facilitates the coplanar orientation of molecules within the heterodimer, the side N–H···O interaction is single-centered, so the heterodimers can rotate around the bond. An indicative parameter of the T-packing is the dihedral angle between the two nearest heterodimers, which varies from 67° in [RLZ+2,6DHBA] to 87° in [RLZ+2HBA]. The difference in packing angles reflects the difference in the mutual positions of donors and acceptors of side H-bonds within the heterodimer and the influence of weak interactions on the packing architecture.

The variety of H-bonds formed in the salts under study can be rationalized based on the comparison of the minima on the MEP surface of RLZ and its counterions (Figure 3). The relative strength of the acceptor sites at the hydroxyl oxygen atom in the ortho-position (O2), is larger than that of hydroxyl oxygens in other positions (O4 and O5), which makes it a preferred target for H-bonding with the N_a–H_c donor group remaining after the main heterosynthon is formed.

According to the results of the MEP analysis, the molecules have several donor and acceptor sites, that may be rated by strength in accordance with their electrostatic potential values. The isolated RLZ molecule has three acceptors and two H-bond donors, with the heterocyclic nitrogen being considerably stronger than the ether oxygen and the amine nitrogen. 2HBA has one strong donor and two acceptor sites at the carbonyl oxygen and hydroxyl oxygen atoms, while dihydroxybenzoic acids have three strong acceptor sites. Due to the formation of the second intramolecular hydrogen bond, the number of available donors and acceptors in the 2,6DHBA molecule is lower by one than in 2,4DHBA and 2,5DHBA. The competition between acceptor sites with close MEP values at the carbonyl oxygen atom of the carboxylic group and the oxygen atom of the side hydroxyl groups, is observed in the (di)hydroxybenzoic acids. In 2,4DHBA and 2,5DHBA, the hydroxyl oxygen in the ortho-position is stronger than the oxygen atoms of the free hydroxyl groups.

Due to steric factors, the preferred interaction in the crystal is the carboxyl–amino-N-heterocycle heterosynthon [82,91] with the two strongest donor/acceptor sites in the RLZ molecule, e.g., N_a–H and N_p–H (for atom labeling see Figure 1 and Figure 2). Compared to (di)hydroxybenzoic acids, RLZ has only one good acceptor site (the N atom of the thiazole ring), which becomes a donor after protonation. After the main amino-N-heterocycle carboxylate heterosynthon is formed, RLZ can only participate in the H-bonding, since its remaining acceptors (the lone pairs at the –OCF₃ oxygen and the –NH₂ nitrogen) have lower absolute electrostatic potential values than the acceptor sites of the hydroxybenzoic acids. The formation of H-bonds occurs in accordance with Etter’s third rule, as the strongest donors bind to the strongest available acceptors [92]. The remaining strong donor in the RLZ molecule (N_a–H_c) is then attached to the stronger of two competing H-bond acceptors: the O2 atom in the salts with 2,4DHBA and 2,5DHBA, and the O6 atom in the salt with 2,6DHBA.

The enthalpies, ΔH_HB, of the intermolecular H-bonds in the considered crystals, evaluated using Equation (1), are given in Table 1. As follows from the presented data, the total enthalpy of the amino-N-heterocycle carboxylate supramolecular heterosynthon (65–70 kJ·mol⁻¹) does not change significantly depending on the number of OH groups in the substituted benzoic acid or the number of intramolecular hydrogen bonds. The obtained values of the total enthalpy are in good agreement with the available literature data [42,52,93]. On the other hand, the [RLZ+2,4DHBA] (1:1) and [RLZ+2,5DHBA] (1:1) salts feature an additional relatively short O–H···O⁻ bond, with an enthalpy above 30 kJ·mol⁻¹ (Table 1). As a result, the total enthalpy of intermolecular H-bonds in [RLZ+2,6DHBA] (92 kJ·mol⁻¹) is lower compared to that of [RLZ+2,4DHBA] (125 kJ·mol⁻¹) and [RLZ+2,5DHBA] (130 kJ·mol⁻¹), because of the lower number of intermolecular H-bonds. It follows from Table 1, that both approaches reproduce the experimental values of distances between heavy atoms involved in intermolecular H-bonds quite well.

The total energy of H-bonds is about 45% of the crystal lattice energy of pharmaceutical multicomponent crystals [63,65,94]. This fact suggests that the energy of the [RLZ+2,6DHBA] crystal lattice is significantly lower than the corresponding values for [RLZ+2,4DHBA] and [RLZ+2,5DHBA] (Table 2). The point is, that the lattice energy of [RLZ+2,6DHBA] is lower than that of [RLZ+2HBA].

4. Discussion

The formation of salts in the RLZ-2,6DHBA and RLZ-2,4DHBA systems supports our hypothesis on the influence of the pK_a of components on governing the ionization state. As seen from Table S1, the ΔpK_a values for RLZ-2,6DHBA and RLZ-2,4DHBA pairs are positive, which agrees with the trend observed for the known RLZ crystals. In brief, hydroxybenzoic acids with –OH groups in positions other than 2 and 6 have pK_a values larger than protonated RLZ, and form cocrystals, whereas ortho-hydroxybenzoic acids result in positive ΔpK_a with RLZ, and form salts.

Although the solubility of a single-component organic compound can be estimated using the Yalkowsky equation [95], or other QSAR models [96,97], with a fair degree of precision, forecasting the solubility of pharmaceutical salts containing different types of counterions is not a trivial task. In fact, only a handful of studies have looked at this issue in depth, examining the effect of structural and physicochemical features of organic salts on their solubility in a systematic manner [98,99,100,101]. The generally accepted view on this topic, implies that the change in the free energy of salt dissolution can be represented as the sum of contributions from the values of the free energy of hydration of the cation and anion, as well as the free energy of the crystal [101,102]:

Δ_solG(salt) = Δ_hydG(cation) + Δ_hydG(anion) − Δ_lattG(salt),(4)

In the case of salts of riluzole with hydroxyl derivatives of salicylic acid, the first term on the right side of Equation (4) is a constant, and the second term can also be considered a constant in the first approximation. Thus, the differences in the solubility of the considered salts are mainly due to the third contribution on the right side of Equation (4). From the point of view of the “enthalpy-determined drug solubility processes” [103], Δ_lattG(salt) can be replaced by Δ_lattH(salt). The latter value is proportional to the crystal lattice energy of the salt. Given that [RLZ+2,6DHBA] has the lowest crystal lattice energy among all the salts under consideration, it might be reasonable to assume that [RLZ+2,6DHBA] is likely to have a better solubility compared to other salts. This conclusion is consistent with the literature data. For example, the solubility of the ethionamide salt with 2,6DHBA, was found to be higher than that with 2,3-dihydroxybenzoic and 2,4-dinitrobenzoic acids [104]. To summarize, 2,6DHBA appears to be a promising compound to consider during cocrystal or salt screening trials, due to its ability to modify the crystalline environment of the target drug without compromising total lattice energy gain and, thus, solubility.

5. Conclusions

In this work, specific features of hydrogen bond topology, induced by the introduction of a second OH group into the counterion molecule, were revealed through a detailed structural analysis of multicomponent crystals of riluzole with 2-hydroxybenzoic acid and 2,4-, 2,5-, and 2,6-dihydroxybenzoic acids. In the solid forms containing 2-hydroxybenzoic acid and 2,6-dihydroxybenzoic acid, the hydrogen bonding network is sustained solely by the N–H groups of riluzole, resulting in the formation of 1D chains. The competition between the acceptor centers of the counterion causes a change in the topology of hydrogen bonds and differences in the packing arrangements of salts with 2-hydroxybenzoic and 2,6-dihydroxybenzoic acids. The presence of an additional H-bond donor at position 4 or 5, leads to the formation of charge-assisted H-bonds and the combination of isolated hydrogen-bonded chains into 2D sheets.

According to periodic DFT calculations, the number and position of the OH-groups in the counterion molecule appear to have little effect on the enthalpy of the primary supramolecular synthon (65–70 kJ·mol⁻¹). A novel intermolecular charge-assisted H-bond in the salts containing 2,4- and 2,5-dihydroxybenzoic acids, stabilizes the corresponding crystal structure and increases the lattice energy by about 30 kJ·mol⁻¹. When compared to riluzole salts with 2,4- and 2,5-dihydroxybenzoic acids, the salt with 2,6-dihydroxybenzoic acid lacks an intermolecular charged-assisted O–H···O⁻ H-bond, resulting in a lower lattice energy and expected higher solubility. In addition, 2,6-dihydroxybenzoic acid is preferred over salicylic acid, because the lattice energy of riluzole salt with 2,6-dihydroxybenzoic acid is lower than that of riluzole salt with 2-hydroxybenzoic acid. We do hope that this prediction encourages experimentalists to conduct further investigations of riluzole multicomponent crystals with dihydroxybenzoic acids, especially 2,6-dihydroxybenzoic acid.

Author Contributions

Conceptualization, A.O.S. and M.V.V.; experimental methodology, A.O.S. and A.P.V.; theoretical methodology, M.V.V.; investigation, A.P.V., A.V.C. and M.V.V.; single-crystal XRD experiment, A.V.C.; writing—original draft preparation, A.P.V. and M.V.V.; writing—review and editing, A.O.S. and M.V.V.; visualization, A.P.V.; supervision, M.V.V. and A.O.S. All authors have read and agreed to the published version of the manuscript.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

Samples of salts obtained in this study, PXRD, DSC, and FTIR data and I/O files are available from the respective author upon reasonable request. Crystal structures of [RLZ+2,4DHBA] and [RLZ+2,6DHBA] have been deposited with the Cambridge Crystallographic Data Centre, under the CCDC numbers 2238373-2238374.

Acknowledgments

We thank “the Upper Volga Region Centre of Physicochemical Research” for technical assistance with PXRD, TGA, and FTIR experiments. The single-crystal X-ray diffraction studies were performed at the Centre of Shared Equipment of IGIC RAS. M.V.V. thanks S.V. Artobolevskii, O.A. Alatortsev, and D.E. Makhrov for pointing his attention to this topic.

Conflicts of Interest

The authors declare no conflict of interest.

Footnotes

Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Figures, Scheme and Tables

Scheme 1. Structural formulas of riluzole and (di)hydroxybenzoic acids considered in this study.

Figure 1. The H-bonded chain in crystalline [RLZ+2,6DHBA] (1:1). The asymmetric unit contains two inequivalent RLZ+2,6DHBA heterodimers. H-bonds are shown as dotted lines. Numbers indicate the enthalpies of intermolecular H-bonds estimated using Equation (1) at the PBE-D3 level. Hc denotes the H atom of the amino group of the RLZ molecule in the cis-position with respect to the S atom (Scheme 1).

Figure 2. Fragment of the hydrogen bond network in crystalline [RLZ+2,4DHBA] (1:1). H-bonds are shown as dotted lines. Numbers indicate the enthalpies of intermolecular H-bonds estimated using Equation (1) at the PBE-D3 level.

Figure 3. MEP mapped over total electron density, ρ = 0.002, isosurface for relaxed molecules of RLZ and the selected counterions. The isosurfaces are color-coded in the range between −0.3 a.u. (blue) and +0.3 a.u. (red). The numbers display the selected electrostatic potential values of minima (cyan points) and maxima (yellow points) of selected electrostatic potentials values, in kJ mol−1.

Table 1

The difference between the calculated (B3LYP/6-31G** and PBE-D3/6-31G**) and experimental values of the distances between heavy atoms ΔD in the considered crystals. The enthalpies ΔH_HB (kJ·mol⁻¹) of these H-bonds, evaluated using the Rozenberg approach, see Equation (1). The PBE-D3 values are given in italics.

	[RLZ+2HBA]		[RLZ+2,6DHBA] (b)		[RLZ+2,5DHBA]		[RLZ+2,4DHBA]
Fragment (a)	ΔD, Å	−ΔHHB	ΔD, Å	−ΔHHB	ΔD, Å	−ΔHHB	ΔD, Å	−ΔHHB
Np–H∙∙∙Op	0.0460.003	32.136.3	−0.026−0.024	35.836.6	0.001−0.038	34.438.7	−0.001−0.016	35.438.0
Na–H∙∙∙Oa	0.009−0.017	29.330.3	−0.008−0.028	28.530.4	0.004−0.005	30.131.1	0.014−0.002	25.626.9
Na–Hc∙∙∙O’p	0.002−0.017	27.428.8	-	-	-	-	-	-
Na–Hc∙∙∙O’6	-	-	0.0820.030	22.524.6	-	-	-	-
Na–Hc∙∙∙O’2	-	-	-	-	−0.022−0.051	26.928.8	−0.019−0.042	22.724.4
O’5–H∙∙∙Op	-	-	-	-	−0.021−0.042	29.731.9	-	-
O’4–H∙∙∙Op	-	-	-	-	-	-	−0.019−0.051	32.635.6

^(a) See Figure 1 and Figure 2 for atom numbering; H_c denotes the H atom of the amino group of the RLZ molecule in the cis-position with respect to the S atom (Scheme 1). ^(b) The values for [RLZ+2,6DHBA] are averaged over two heterodimers in the asymmetric unit.

Table 2

Total energy of intermolecular H-bonds in the crystals of RLZ salts obtained using Equation (2), and lattice energy estimated using Equation (3). Geometric parameters were calculated at the PBE-D3 level. The units are kJ·mol⁻¹.

Crystal	[RLZ+2HBA]	[RLZ+2,4DHBA]	[RLZ+2,5DHBA]	[RLZ+2,6DHBA]
ΣEHB	95.4	124.9	130.7	91.6
Elatt	229.0	281.9	266.0	209.4

Supplementary Materials

The following supporting information can be downloaded at: https://www.mdpi.com/article/10.3390/pharmaceutics15030878/s1, Section S1: Analysis of the literature data on crystal structures of riluzole multicomponent crystals with benzoic acid derivatives; Table S1: ΔpKa and hydrogen bond parameters in known crystals of riluzole with benzoic acid derivatives; Section S2: Details of PXRD, DSC, TGA, and FTIR experiments; Section S3: Phase identification and crystal structure analysis; Figures S1–S3: PXRD patterns of RLZ mixtures with dihydroxybenzoic acids after liquid-assisted grinding; Table S2: Crystallographic data for [RLZ+2,4DHBA] and [RLZ+2,6DHBA]; Figure S4: FTIR spectra of RLZ and its salts; Figures S5–S7: DSC and TGA traces of RLZ salts. References [105,106] are cited in the supplementary materials.

References

1. Gunawardana, C.A.; Aakeröy, C.B. Co-crystal synthesis: Fact, fancy, and great expectations. Chem. Commun.; 2018; 54, pp. 14047-14060. [DOI: https://dx.doi.org/10.1039/C8CC08135B] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/30488919]

2. Kavanagh, O.N.; Croker, D.M.; Walker, G.M.; Zaworotko, M.J. Pharmaceutical cocrystals: From serendipity to design to application. Drug Discov. Today; 2019; 24, pp. 796-804. [DOI: https://dx.doi.org/10.1016/j.drudis.2018.11.023] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/30521935]

3. Bolla, G.; Sarma, B.; Nangia, A.K. Crystal Engineering of Pharmaceutical Cocrystals in the Discovery and Development of Improved Drugs. Chem. Rev.; 2022; 122, pp. 11514-11603. [DOI: https://dx.doi.org/10.1021/acs.chemrev.1c00987] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/35642550]

4. Shattock, T.R.; Arora, K.K.; Vishweshwar, P.; Zaworotko, M.J. Hierarchy of Supramolecular Synthons: Persistent Carboxylic Acid·Pyridine Hydrogen Bonds in Cocrystals That also Contain a Hydroxyl Moiety. Cryst. Growth Des.; 2008; 8, pp. 4533-4545. [DOI: https://dx.doi.org/10.1021/cg800565a]

5. Mir, N.A.; Dubey, R.; Desiraju, G.R. Strategy and Methodology in the Synthesis of Multicomponent Molecular Solids: The Quest for Higher Cocrystals. Acc. Chem. Res.; 2019; 52, pp. 2210-2220. [DOI: https://dx.doi.org/10.1021/acs.accounts.9b00211]

6. Topić, F.; Rissanen, K. Systematic Construction of Ternary Cocrystals by Orthogonal and Robust Hydrogen and Halogen Bonds. J. Am. Chem. Soc.; 2016; 138, pp. 6610-6616. [DOI: https://dx.doi.org/10.1021/jacs.6b02854]

7. Aakeröy, C.B.; Panikkattu, S.; Chopade, P.D.; Desper, J. Competing hydrogen-bond and halogen-bond donors in crystal engineering. CrystEngComm; 2013; 15, pp. 3125-3136. [DOI: https://dx.doi.org/10.1039/C2CE26747K]

8. Desiraju, G.R. Hydrogen Bridges in Crystal Engineering:  Interactions without Borders. Acc. Chem. Res.; 2002; 35, pp. 565-573. [DOI: https://dx.doi.org/10.1021/ar010054t]

9. Etter, M.C. Hydrogen bonds as design elements in organic chemistry. J. Phys. Chem.; 1991; 95, pp. 4601-4610. [DOI: https://dx.doi.org/10.1021/j100165a007]

10. Aakeröy, C.B.; Seddon, K.R. The hydrogen bond and crystal engineering. Chem. Soc. Rev.; 1993; 22, pp. 397-407. [DOI: https://dx.doi.org/10.1039/CS9932200397]

11. Moulton, B.; Zaworotko, M.J. From Molecules to Crystal Engineering:  Supramolecular Isomerism and Polymorphism in Network Solids. Chem. Rev.; 2001; 101, pp. 1629-1658. [DOI: https://dx.doi.org/10.1021/cr9900432] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/11709994]

12. Karas, L.J.; Wu, C.-H.; Das, R.; Wu, J.I.C. Hydrogen bond design principles. WIREs Comput. Mol. Sci.; 2020; 10, e1477. [DOI: https://dx.doi.org/10.1002/wcms.1477] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/33936251]

13. Desiraju, G.R. Crystal Engineering: A Holistic View. Angew. Chem. Int. Ed.; 2007; 46, pp. 8342-8356. [DOI: https://dx.doi.org/10.1002/anie.200700534] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/17902079]

14. Lo Presti, L. On the significance of weak hydrogen bonds in crystal packing: A large databank comparison of polymorphic structures. CrystEngComm; 2018; 20, pp. 5976-5989. [DOI: https://dx.doi.org/10.1039/C8CE00674A]

15. Vologzhanina, A.V. Intermolecular Interactions in Functional Crystalline Materials: From Data to Knowledge. Crystals; 2019; 9, 478. [DOI: https://dx.doi.org/10.3390/cryst9090478]

16. Moragues-Bartolome, A.M.; Jones, W.; Cruz-Cabeza, A.J. Synthon preferences in cocrystals of cis-carboxamides:carboxylic acids. CrystEngComm; 2012; 14, pp. 2552-2559. [DOI: https://dx.doi.org/10.1039/c2ce06241k]

17. Aakeröy, C.B.; Baldrighi, M.; Desper, J.; Metrangolo, P.; Resnati, G. Supramolecular Hierarchy among Halogen-Bond Donors. Chem. A Eur. J.; 2013; 19, pp. 16240-16247. [DOI: https://dx.doi.org/10.1002/chem.201302162]

18. Otero-de-la-Roza, A.; Johnson, E.R. A benchmark for non-covalent interactions in solids. J. Chem. Phys.; 2012; 137, 054103. [DOI: https://dx.doi.org/10.1063/1.4738961]

19. Vener, M.V.; Levina, E.O.; Koloskov, O.A.; Rykounov, A.A.; Voronin, A.P.; Tsirelson, V.G. Evaluation of the Lattice Energy of the Two-Component Molecular Crystals Using Solid-State Density Functional Theory. Cryst. Growth Des.; 2014; 14, pp. 4997-5003. [DOI: https://dx.doi.org/10.1021/cg5005243]

20. Taylor, R. Which intermolecular interactions have a significant influence on crystal packing?. CrystEngComm; 2014; 16, pp. 6852-6865. [DOI: https://dx.doi.org/10.1039/C4CE00452C]

21. Miller, D.K.; Chernyshov, I.Y.; Torubaev, Y.V.; Rosokha, S.V. From weak to strong interactions: Structural and electron topology analysis of the continuum from the supramolecular chalcogen bonding to covalent bonds. Phys. Chem. Chem. Phys.; 2022; 24, pp. 8251-8259. [DOI: https://dx.doi.org/10.1039/D1CP05441D] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/35320823]

22. Johnson, S.A.; Fang, T.; De Marchi, F.; Neel, D.; Van Weehaeghe, D.; Berry, J.D.; Paganoni, S. Pharmacotherapy for Amyotrophic Lateral Sclerosis: A Review of Approved and Upcoming Agents. Drugs; 2022; 82, pp. 1367-1388. [DOI: https://dx.doi.org/10.1007/s40265-022-01769-1]

23. Bensimon, G.; Lacomblez, L.; Meininger, V. A Controlled Trial of Riluzole in Amyotrophic Lateral Sclerosis. N. Engl. J. Med.; 1994; 330, pp. 585-591. [DOI: https://dx.doi.org/10.1056/NEJM199403033300901] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/8302340]

24. Lacomblez, L.; Bensimon, G.; Leigh, P.N.; Guillet, P.; Meininger, V. Dose-ranging study of riluzole in amyotrophic lateral sclerosis. Amyotrophic Lateral Sclerosis/Riluzole Study Group II. Lancet; 1996; 347, pp. 1425-1431. [DOI: https://dx.doi.org/10.1016/S0140-6736(96)91680-3] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/8676624]

25. Brennan, B.P.; Hudson, J.I.; Jensen, J.E.; McCarthy, J.; Roberts, J.L.; Prescot, A.P.; Cohen, B.M.; Pope, H.G.; Renshaw, P.F.; Öngür, D. Rapid Enhancement of Glutamatergic Neurotransmission in Bipolar Depression Following Treatment with Riluzole. Neuropsychopharmacology; 2010; 35, pp. 834-846. [DOI: https://dx.doi.org/10.1038/npp.2009.191] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/19956089]

26. Urbani, A.; Belluzzi, O. Riluzole inhibits the persistent sodium current in mammalian CNS neurons. Eur. J. Neurosci.; 2000; 12, pp. 3567-3574. [DOI: https://dx.doi.org/10.1046/j.1460-9568.2000.00242.x]

27. Benoit, E.; Escande, D. Riluzole specifically blocks inactivated Na channels in myelinated nerve fibre. Pflügers Arch.; 1991; 419, pp. 603-609. [DOI: https://dx.doi.org/10.1007/BF00370302]

28. Centonze, D.; Calabresi, P.; Pisani, A.; Marinelli, S.; Marfia, G.A.; Bernardi, G. Electrophysiology of the neuroprotective agent riluzole on striatal spiny neurons. Neuropharmacology; 1998; 37, pp. 1063-1070. [DOI: https://dx.doi.org/10.1016/S0028-3908(98)00081-1]

29. Mizuta, I.; Ohta, M.; Ohta, K.; Nishimura, M.; Mizuta, E.; Kuno, S. Riluzole stimulates nerve growth factor, brain-derived neurotrophic factor and glial cell line-derived neurotrophic factor synthesis in cultured mouse astrocytes. Neurosci. Lett.; 2001; 310, pp. 117-120. [DOI: https://dx.doi.org/10.1016/S0304-3940(01)02098-5]

30. Van Den Bosch, L.; Van Damme, P.; Bogaert, E.; Robberecht, W. The role of excitotoxicity in the pathogenesis of amyotrophic lateral sclerosis. Biochim. Biophys. Acta (BBA) Mol. Basis Dis.; 2006; 1762, pp. 1068-1082. [DOI: https://dx.doi.org/10.1016/j.bbadis.2006.05.002]

31. Cheah, B.C.; Vucic, S.; Krishnan, A.V.; Kiernan, M.C. Riluzole, Neuroprotection and Amyotrophic Lateral Sclerosis. Curr. Med. Chem.; 2010; 17, pp. 1942-1959. [DOI: https://dx.doi.org/10.2174/092986710791163939] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/20377511]

32. Zarate, C.A., Jr.; Manji, H.K. Riluzole in psychiatry: A systematic review of the literature. Expert Opin. Drug Metab. Toxicol.; 2008; 4, pp. 1223-1234. [DOI: https://dx.doi.org/10.1517/17425255.4.9.1223] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/18721116]

33. Pittenger, C.; Coric, V.; Banasr, M.; Bloch, M.; Krystal, J.H.; Sanacora, G. Riluzole in the Treatment of Mood and Anxiety Disorders. CNS Drugs; 2008; 22, pp. 761-786. [DOI: https://dx.doi.org/10.2165/00023210-200822090-00004] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/18698875]

34. Blyufer, A.; Lhamo, S.; Tam, C.; Tariq, I.; Thavornwatanayong, T.; Mahajan, S.S. Riluzole: A neuroprotective drug with potential as a novel anti-cancer agent (Review). Int. J. Oncol.; 2021; 59, 95. [DOI: https://dx.doi.org/10.3892/ijo.2021.5275] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/34713302]

35. Khan, A.J.; LaCava, S.; Mehta, M.; Schiff, D.; Thandoni, A.; Jhawar, S.; Danish, S.; Haffty, B.G.; Chen, S. The glutamate release inhibitor riluzole increases DNA damage and enhances cytotoxicity in human glioma cells, in vitro and in vivo. Oncotarget; 2019; 10, pp. 2824-2834. [DOI: https://dx.doi.org/10.18632/oncotarget.26854] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/31073373]

36. Benavides-Serrato, A.; Saunders, J.T.; Holmes, B.; Nishimura, R.N.; Lichtenstein, A.; Gera, J. Repurposing Potential of Riluzole as an ITAF Inhibitor in mTOR Therapy Resistant Glioblastoma. Int. J. Mol. Sci.; 2020; 21, 344. [DOI: https://dx.doi.org/10.3390/ijms21010344]

37. Yu, L.J.; Wall, B.A.; Chen, S. The current management of brain metastasis in melanoma: A focus on riluzole. Expert Rev. Neurother.; 2015; 15, pp. 779-792. [DOI: https://dx.doi.org/10.1586/14737175.2015.1055321]

38. Liboux, A.L.; Lefebvre, P.; Roux, Y.L.; Truffinet, P.; Aubeneau, M.; Kirkesseli, S.; Montay, G. Single- and Multiple-Dose Pharmacokinetics of Riluzole in White Subjects. J. Clin. Pharmacol.; 1997; 37, pp. 820-827. [DOI: https://dx.doi.org/10.1002/j.1552-4604.1997.tb05630.x]

39. Bryson, H.M.; Fulton, B.; Benfield, P. Riluzole. Drugs; 1996; 52, pp. 549-563. [DOI: https://dx.doi.org/10.2165/00003495-199652040-00010]

40. Mondal, P.K.; Rao, V.; Mittapalli, S.; Chopra, D. Exploring Solid State Diversity and Solution Characteristics in a Fluorine-Containing Drug Riluzole. Cryst. Growth Des.; 2017; 17, pp. 1938-1946. [DOI: https://dx.doi.org/10.1021/acs.cgd.6b01894]

41. Kawabata, Y.; Wada, K.; Nakatani, M.; Yamada, S.; Onoue, S. Formulation design for poorly water-soluble drugs based on biopharmaceutics classification system: Basic approaches and practical applications. Int. J. Pharm.; 2011; 420, pp. 1-10. [DOI: https://dx.doi.org/10.1016/j.ijpharm.2011.08.032] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/21884771]

42. Yadav, B.; Balasubramanian, S.; Chavan, R.B.; Thipparaboina, R.; Naidu, V.G.M.; Shastri, N.R. Hepatoprotective Cocrystals and Salts of Riluzole: Prediction, Synthesis, Solid State Characterization, and Evaluation. Cryst. Growth Des.; 2018; 18, pp. 1047-1061. [DOI: https://dx.doi.org/10.1021/acs.cgd.7b01514]

43. Thomas, S.P.; Kumar, V.; Alhameedi, K.; Guru Row, T.N. Non-Classical Synthons: Supramolecular Recognition by S⋅⋅⋅O Chalcogen Bonding in Molecular Complexes of Riluzole. Chem. A Eur. J.; 2019; 25, pp. 3591-3597. [DOI: https://dx.doi.org/10.1002/chem.201805131] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/30576020]

44. Verma, S.K.; Arora, I.; Javed, K.; Akhtar, M.; Samim, M. Enhancement in the Neuroprotective Power of Riluzole Against Cerebral Ischemia Using a Brain Targeted Drug Delivery Vehicle. ACS Appl. Mater. Interfaces; 2016; 8, pp. 19716-19723. [DOI: https://dx.doi.org/10.1021/acsami.6b01776] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/27378322]

45. Patel, R.J.; Patel, A.A.; Patel, H.P. Stabilized amorphous state of riluzole by immersion-rotavapor method with synthesized mesoporous SBA-15 carrier to augment in-vitro dissolution. J. Drug Deliv. Sci. Technol.; 2021; 61, 102270. [DOI: https://dx.doi.org/10.1016/j.jddst.2020.102270]

46. Povedano Panades, M.; Couratier, P.; Sidle, K.; Sorarù, G.; Tsivgoulis, G.; Ludolph, A.C. Administration of Riluzole Oral Suspension During the Different Stages of Amyotrophic Lateral Sclerosis. Front. Neurol.; 2021; 12, 633854. [DOI: https://dx.doi.org/10.3389/fneur.2021.633854]

47. Wymer, J.; Apple, S.; Harrison, A.; Hill, B.A. Pharmacokinetics, Bioavailability, and Swallowing Safety With Riluzole Oral Film. Clin. Pharmacol. Drug Dev.; 2023; 12, pp. 57-64. [DOI: https://dx.doi.org/10.1002/cpdd.1168]

48. Parikh, R.H.; Patel, R.J. Nanoemulsions for Intranasal Delivery of Riluzole to Improve Brain Bioavailability: Formulation Development and Pharmacokinetic Studies. Curr. Drug Deliv.; 2016; 13, pp. 1130-1143. [DOI: https://dx.doi.org/10.2174/1567201813666151202195729]

49. Sroka, Z.; Cisowski, W. Hydrogen peroxide scavenging, antioxidant and anti-radical activity of some phenolic acids. Food Chem. Toxicol.; 2003; 41, pp. 753-758. [DOI: https://dx.doi.org/10.1016/S0278-6915(02)00329-0]

50. Reid, J.; Watson, R.D.; Cochran, J.B.; Sproull, D.H.; Clayton, B.E.; Prunty, F.T. Sodium gamma-resorcylate in rheumatic fever. Br. Med. J.; 1951; 2, pp. 321-326. [DOI: https://dx.doi.org/10.1136/bmj.2.4727.321]

51. Taylor, R.; Wood, P.A. A Million Crystal Structures: The Whole Is Greater than the Sum of Its Parts. Chem. Rev.; 2019; 119, pp. 9427-9477. [DOI: https://dx.doi.org/10.1021/acs.chemrev.9b00155]

52. Surov, A.O.; Vasilev, N.A.; Churakov, A.V.; Parashchuk, O.D.; Artobolevskii, S.V.; Alatortsev, O.A.; Makhrov, D.E.; Vener, M.V. Two Faces of Water in the Formation and Stabilization of Multicomponent Crystals of Zwitterionic Drug-Like Compounds. Symmetry; 2021; 13, 425. [DOI: https://dx.doi.org/10.3390/sym13030425]

53. Karamertzanis, P.G.; Day, G.M.; Welch, G.W.A.; Kendrick, J.; Leusen, F.J.J.; Neumann, M.A.; Price, S.L. Modeling the interplay of inter- and intramolecular hydrogen bonding in conformational polymorphs. J. Chem. Phys.; 2008; 128, 244708. [DOI: https://dx.doi.org/10.1063/1.2937446] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/18601366]

54. Sheldrick, G. SADABS, Program for Scaling and Correction of Area Detector Data; University of Göttingen: Göttingen, Germany, 1997.

55. Sheldrick, G. A short history of SHELX. Acta Crystallogr. Sect. A Found. Crystallogr.; 2008; 64, pp. 112-122. [DOI: https://dx.doi.org/10.1107/S0108767307043930] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/18156677]

56. Becke, A.D. Density-functional thermochemistry. III. The role of exact exchange. J. Chem. Phys.; 1993; 98, pp. 5648-5652. [DOI: https://dx.doi.org/10.1063/1.464913]

57. Lee, C.; Yang, W.; Parr, R.G. Development of the Colle-Salvetti correlation-energy formula into a functional of the electron density. Phys. Rev. B; 1988; 37, pp. 785-789. [DOI: https://dx.doi.org/10.1103/PhysRevB.37.785] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/9944570]

58. Perdew, J.P.; Burke, K.; Ernzerhof, M. Generalized Gradient Approximation Made Simple. Phys. Rev. Lett.; 1996; 77, pp. 3865-3868. [DOI: https://dx.doi.org/10.1103/PhysRevLett.77.3865]

59. Grimme, S.; Ehrlich, S.; Goerigk, L. Effect of the damping function in dispersion corrected density functional theory. J. Comput. Chem.; 2011; 32, pp. 1456-1465. [DOI: https://dx.doi.org/10.1002/jcc.21759]

60. Deringer, V.L.; George, J.; Dronskowski, R.; Englert, U. Plane-Wave Density Functional Theory Meets Molecular Crystals: Thermal Ellipsoids and Intermolecular Interactions. Acc. Chem. Res.; 2017; 50, pp. 1231-1239. [DOI: https://dx.doi.org/10.1021/acs.accounts.7b00067]

61. Cutini, M.; Civalleri, B.; Corno, M.; Orlando, R.; Brandenburg, J.G.; Maschio, L.; Ugliengo, P. Assessment of Different Quantum Mechanical Methods for the Prediction of Structure and Cohesive Energy of Molecular Crystals. J. Chem. Theory Comput.; 2016; 12, pp. 3340-3352. [DOI: https://dx.doi.org/10.1021/acs.jctc.6b00304] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/27304925]

62. Mohaček-Grošev, V.; Grdadolnik, J.; Stare, J.; Hadži, D. Identification of hydrogen bond modes in polarized Raman spectra of single crystals of α-oxalic acid dihydrate. J. Raman Spectrosc.; 2009; 40, pp. 1605-1614. [DOI: https://dx.doi.org/10.1002/jrs.2308]

63. Surov, A.O.; Manin, A.N.; Voronin, A.P.; Churakov, A.V.; Perlovich, G.L.; Vener, M.V. Weak Interactions Cause Packing Polymorphism in Pharmaceutical Two-Component Crystals. The Case Study of the Salicylamide Cocrystal. Cryst. Growth Des.; 2017; 17, pp. 1425-1437. [DOI: https://dx.doi.org/10.1021/acs.cgd.7b00019]

64. Zhang, F.; Kambara, O.; Tominaga, K.; Nishizawa, J.-I.; Sasaki, T.; Wang, H.-W.; Hayashi, M. Analysis of vibrational spectra of solid-state adenine and adenosine in the terahertz region. RSC Adv.; 2014; 4, pp. 269-278. [DOI: https://dx.doi.org/10.1039/C3RA44285C]

65. Surov, A.O.; Voronin, A.P.; Vener, M.V.; Churakov, A.V.; Perlovich, G.L. Specific features of supramolecular organisation and hydrogen bonding in proline cocrystals: A case study of fenamates and diclofenac. CrystEngComm; 2018; 20, pp. 6970-6981. [DOI: https://dx.doi.org/10.1039/C8CE01458B]

66. Masunov, A.E.; Torres, K.; Dyakov, A.A.; Yushina, I.D.; Bartashevich, E.V. First-Principles Crystal Engineering of Nonlinear Optical Materials. II. Effect of Halogen Bonds on the Structure and Properties of Triiodobenzenes. J. Phys. Chem. C; 2018; 122, pp. 22622-22631. [DOI: https://dx.doi.org/10.1021/acs.jpcc.8b04932]

67. Vener, M.V.; Chernyshov, I.Y.; Rykounov, A.A.; Filarowski, A. Structural and spectroscopic features of proton hydrates in the crystalline state. Solid-state DFT study on HCl and triflic acid hydrates. Mol. Phys.; 2018; 116, pp. 251-262. [DOI: https://dx.doi.org/10.1080/00268976.2017.1380860]

68. Vener, M.V.; Kharlanov, O.G.; Sosorev, A.Y. High-Mobility Naphthalene Diimide Derivatives Revealed by Raman-Based In Silico Screening. Int. J. Mol. Sci.; 2022; 23, 13305. [DOI: https://dx.doi.org/10.3390/ijms232113305]

69. Dovesi, R.; Erba, A.; Orlando, R.; Zicovich-Wilson, C.M.; Civalleri, B.; Maschio, L.; Rérat, M.; Casassa, S.; Baima, J.; Salustro, S. et al. Quantum-mechanical condensed matter simulations with CRYSTAL. Wiley Interdiscip. Rev. Comput. Mol. Sci.; 2018; 8, e1360. [DOI: https://dx.doi.org/10.1002/wcms.1360]

70. Surov, A.O.; Voronin, A.P.; Drozd, K.V.; Gruzdev, M.S.; Perlovich, G.L.; Prashanth, J.; Balasubramanian, S. Polymorphic forms of antiandrogenic drug nilutamide: Structural and thermodynamic aspects. Phys. Chem. Chem. Phys.; 2021; 23, pp. 9695-9708. [DOI: https://dx.doi.org/10.1039/D1CP00793A]

71. Rogers, F.J.M.; Radhanpura, K.; Horvat, J.; Farrant, D. On the use of a volume constraint to account for thermal expansion effects on the low-frequency vibrations of molecular crystals. Phys. Chem. Chem. Phys.; 2022; 24, pp. 10408-10419. [DOI: https://dx.doi.org/10.1039/D1CP05718A]

72. Rozenberg, M.; Loewenschuss, A.; Marcus, Y. An empirical correlation between stretching vibration redshift and hydrogen bond length. Phys. Chem. Chem. Phys.; 2000; 2, pp. 2699-2702. [DOI: https://dx.doi.org/10.1039/b002216k]

73. Mata, I.; Alkorta, I.; Espinosa, E.; Molins, E. Relationships between interaction energy, intermolecular distance and electron density properties in hydrogen bonded complexes under external electric fields. Chem. Phys. Lett.; 2011; 507, pp. 185-189. [DOI: https://dx.doi.org/10.1016/j.cplett.2011.03.055]

74. Bader, R.F.W. A quantum theory of molecular structure and its applications. Chem. Rev.; 1991; 91, pp. 893-928. [DOI: https://dx.doi.org/10.1021/cr00005a013]

75. Melikova, S.M.; Voronin, A.P.; Panek, J.; Frolov, N.E.; Shishkina, A.V.; Rykounov, A.A.; Tretyakov, P.Y.; Vener, M.V. Interplay of π-stacking and inter-stacking interactions in two-component crystals of neutral closed-shell aromatic compounds: Periodic DFT study. RSC Adv.; 2020; 10, pp. 27899-27910. [DOI: https://dx.doi.org/10.1039/D0RA04799F]

76. Gatti, C.; Casassa, S.M. TOPOND14 User’s Manual; CNR-ISTM Milano: Milano, Italy, 2013.

77. Medvedev, A.G.; Churakov, A.V.; Navasardyan, M.A.; Prikhodchenko, P.V.; Lev, O.; Vener, M.V. Fast Quantum Approach for Evaluating the Energy of Non-Covalent Interactions in Molecular Crystals: The Case Study of Intermolecular H-Bonds in Crystalline Peroxosolvates. Molecules; 2022; 27, 4082. [DOI: https://dx.doi.org/10.3390/molecules27134082]

78. Vener, M.V.; Makhrov, D.E.; Voronin, A.P.; Shalafan, D.R. Molecular Dynamics Simulation of Association Processes in Aqueous Solutions of Maleate Salts of Drug-like Compounds: The Role of Counterion. Int. J. Mol. Sci.; 2022; 23, 6302. [DOI: https://dx.doi.org/10.3390/ijms23116302]

79. Murray, J.S.; Politzer, P. Hydrogen Bonding: A Coulombic σ-Hole Interaction. J. Indian Inst. Sci.; 2020; 100, pp. 21-30. [DOI: https://dx.doi.org/10.1007/s41745-019-00139-3]

80. Frisch, M.J.; Trucks, G.W.; Schlegel, H.B.; Scuseria, G.E.; Robb, M.A.; Cheeseman, J.R.; Scalmani, G.; Barone, V.; Mennucci, B.; Petersson, G.A. et al. Gaussian 09, Revision B.01; Gaussian Inc.: Wallingford, CT, USA, 2009.

81. Lu, T.; Chen, F. Multiwfn: A multifunctional wavefunction analyzer. J. Comput. Chem.; 2012; 33, pp. 580-592. [DOI: https://dx.doi.org/10.1002/jcc.22885]

82. Bis, J.A.; Zaworotko, M.J. The 2-Aminopyridinium-carboxylate Supramolecular Heterosynthon:  A Robust Motif for Generation of Multiple-Component Crystals. Cryst. Growth Des.; 2005; 5, pp. 1169-1179. [DOI: https://dx.doi.org/10.1021/cg049622c]

83. Ebenezer, S.; Muthiah, P.T. Design of Co-crystals/Salts of Aminopyrimidines and Carboxylic Acids through Recurrently Occurring Synthons. Cryst. Growth Des.; 2012; 12, pp. 3766-3785. [DOI: https://dx.doi.org/10.1021/cg3005954]

84. da Silva, C.C.; Cirqueira, M.d.L.; Martins, F.T. Lamivudine salts with 1,2-dicarboxylic acids: A new and a rare synthon with double pairing motif fine-tuning their solubility. CrystEngComm; 2013; 15, pp. 6311-6317. [DOI: https://dx.doi.org/10.1039/c3ce40844b]

85. Voronin, A.P.; Surov, A.O.; Churakov, A.V.; Parashchuk, O.D.; Rykounov, A.A.; Vener, M.V. 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. Molecules; 2020; 25, 2386. [DOI: https://dx.doi.org/10.3390/molecules25102386] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/32455564]

86. Mirzaei, M.; Sadeghi, F.; Molčanov, K.; Zarȩba, J.K.; Gomila, R.M.; Frontera, A. Recurrent Supramolecular Motifs in a Series of Acid–Base Adducts Based on Pyridine-2,5-Dicarboxylic Acid N-Oxide and Organic Bases: Inter- and Intramolecular Hydrogen Bonding. Cryst. Growth Des.; 2020; 20, pp. 1738-1751. [DOI: https://dx.doi.org/10.1021/acs.cgd.9b01475]

87. Vener, M.V.; Churakov, A.V.; Voronin, A.P.; Parashchuk, O.D.; Artobolevskii, S.V.; Alatortsev, O.A.; Makhrov, D.E.; Medvedev, A.G.; Filarowski, A. Comparison of Proton Acceptor and Proton Donor Properties of H₂O and H₂O₂ in Organic Crystals of Drug-like Compounds: Peroxosolvates vs. Crystallohydrates. Molecules; 2022; 27, 717. [DOI: https://dx.doi.org/10.3390/molecules27030717]

88. Garg, U.; Azim, Y.; Kar, A.; Pradeep, C.P. Cocrystals/salt of 1-naphthaleneacetic acid and utilizing Hirshfeld surface calculations for acid–aminopyrimidine synthons. CrystEngComm; 2020; 22, pp. 2978-2989. [DOI: https://dx.doi.org/10.1039/D0CE00106F]

89. Gelbrich, T.; Hursthouse, M.B. A versatile procedure for the identification, description and quantification of structural similarity in molecular crystals. CrystEngComm; 2005; 7, pp. 324-336. [DOI: https://dx.doi.org/10.1039/b502484f]

90. Gelbrich, T.; Hursthouse, M.B. Systematic investigation of the relationships between 25 crystal structures containing the carbamazepine molecule or a close analogue: A case study of the XPac method. CrystEngComm; 2006; 8, pp. 448-460. [DOI: https://dx.doi.org/10.1039/b600753h]

91. Shan, N.; Bond, A.D.; Jones, W. Supramolecular synthons in the co-crystal structures of 2-aminopyrimidine with diols and carboxylic acids. Tetrahedron Lett.; 2002; 43, pp. 3101-3104. [DOI: https://dx.doi.org/10.1016/S0040-4039(02)00511-7]

92. Etter, M.C. Encoding and decoding hydrogen-bond patterns of organic compounds. Acc. Chem. Res.; 1990; 23, pp. 120-126. [DOI: https://dx.doi.org/10.1021/ar00172a005]

93. Tao, Q.; Hao, Q.-Q.; Voronin, A.P.; Dai, X.-L.; Huang, Y.; Perlovich, G.L.; Lu, T.-B.; Chen, J.-M. Polymorphic Forms of a Molecular Salt of Phenazopyridine with 3,5-Dihydroxybenzoic Acid: Crystal Structures, Theoretical Calculations, Thermodynamic Stability, and Solubility Aspects. Cryst. Growth Des.; 2019; 19, pp. 5636-5647. [DOI: https://dx.doi.org/10.1021/acs.cgd.9b00618]

94. Voronin, A.P.; Perlovich, G.L.; Vener, M.V. Effects of the crystal structure and thermodynamic stability on solubility of bioactive compounds: DFT study of isoniazid cocrystals. Comput. Theor. Chem.; 2016; 1092, pp. 1-11. [DOI: https://dx.doi.org/10.1016/j.comptc.2016.07.022]

95. Jain, N.; Yalkowsky, S.H. Estimation of the aqueous solubility I: Application to organic nonelectrolytes. J. Pharm. Sci.; 2001; 90, pp. 234-252. [DOI: https://dx.doi.org/10.1002/1520-6017(200102)90:2<234::AID-JPS14>3.0.CO;2-V] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/11169540]

96. Palmer, D.S.; O’Boyle, N.M.; Glen, R.C.; Mitchell, J.B.O. Random Forest Models To Predict Aqueous Solubility. J. Chem. Inf. Model.; 2007; 47, pp. 150-158. [DOI: https://dx.doi.org/10.1021/ci060164k] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/17238260]

97. Wang, J.; Hou, T. Recent Advances on Aqueous Solubility Prediction. Comb. Chem. High Throughput Screen.; 2011; 14, pp. 328-338. [DOI: https://dx.doi.org/10.2174/138620711795508331]

98. Sacchi, P.; Loconte, L.; Macetti, G.; Rizzato, S.; Lo Presti, L. Correlations of Crystal Structure and Solubility in Organic Salts: The Case of the Antiplasmodial Drug Piperaquine. Cryst. Growth Des.; 2019; 19, pp. 1399-1410. [DOI: https://dx.doi.org/10.1021/acs.cgd.8b01794]

99. de Moraes, L.S.; Edwards, D.; Florence, A.J.; Johnston, A.; Johnston, B.F.; Morrison, C.A.; Kennedy, A.R. Aqueous Solubility of Organic Salts. Investigating Trends in a Systematic Series of 51 Crystalline Salt Forms of Methylephedrine. Cryst. Growth Des.; 2017; 17, pp. 3277-3286. [DOI: https://dx.doi.org/10.1021/acs.cgd.7b00255]

100. Arlin, J.-B.; Florence, A.J.; Johnston, A.; Kennedy, A.R.; Miller, G.J.; Patterson, K. Systematic Data Set for Structure−Property Investigations: Solubility and Solid-State Structure of Alkaline Earth Metal Salts of Benzoates. Cryst. Growth Des.; 2011; 11, pp. 1318-1327. [DOI: https://dx.doi.org/10.1021/cg101547r]

101. Yang, Q.; Ren, T.; Yang, S.; Li, X.; Chi, Y.; Yang, Y.; Gu, J.; Hu, C. Synthesis and Pharmacokinetic Study of Three Gemfibrozil Salts: An Exploration of the Structure–Property Relationship. Cryst. Growth Des.; 2016; 16, pp. 6060-6068. [DOI: https://dx.doi.org/10.1021/acs.cgd.6b01100]

102. Serajuddin, A.T.M. Salt formation to improve drug solubility. Adv. Drug Deliv. Rev.; 2007; 59, pp. 603-616. [DOI: https://dx.doi.org/10.1016/j.addr.2007.05.010]

103. Perlovich, G.L. Thermodynamic Approaches to the Challenges of Solubility in Drug Discovery and Development. Mol. Pharm.; 2014; 11, pp. 1-11. [DOI: https://dx.doi.org/10.1021/mp400460r]

104. Nechipadappu, S.K.; Trivedi, D.R. Pharmaceutical salts of ethionamide with GRAS counter ion donors to enhance the solubility. Eur. J. Pharm. Sci.; 2017; 96, pp. 578-589. [DOI: https://dx.doi.org/10.1016/j.ejps.2016.10.035] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/27815138]

105. Cruz-Cabeza, A.J. Acid–base crystalline complexes and the pKa rule. CrystEngComm; 2012; 14, pp. 6362-6365. [DOI: https://dx.doi.org/10.1039/c2ce26055g]

106. Vecchio, S.; Brunetti, B. Thermochemical study of 2,4-, 2,6- and 3,4-dihydroxybenzoic acids in the liquid phase using a TG apparatus. Thermochim. Acta; 2011; 515, pp. 84-90. [DOI: https://dx.doi.org/10.1016/j.tca.2011.01.001]