Open Access Article
This Open Access Article is licensed under a
Creative Commons Attribution 3.0 Unported Licence

Solvation and aggregation of heteroaromatic drugs in the Reline deep eutectic solvent – a combined molecular dynamics simulation and DFT study

Narges Karimi a, Maryam Heydari Dokoohaki a, Amin Reza Zolghadr *a and Axel Klein *ab
aDepartment of Chemistry, School of Science, Shiraz University, Shiraz 71946-84795, Iran. E-mail: nargec_karimi@yahoo.com; maryamhdi5844@gmail.com; arzolghadr@shirazu.ac.ir
bUniversity of Cologne, Faculty for Mathematics and Natural Sciences, Department of Chemistry and Biochemistry, Institute for Inorganic Chemistry and Materials Chemistry, Köln, Germany. E-mail: axel.klein@uni-koeln.de

Received 5th April 2025 , Accepted 30th June 2025

First published on 30th June 2025


Abstract

Molecular dynamics (MD) simulation and density functional theory (DFT) analyses were carried out on the solubilization of the heteroaromatic drugs allopurinol, losartan, and omeprazole in reline, which is a deep eutectic solvent (DES) composed of choline chloride/urea in 1[thin space (1/6-em)]:[thin space (1/6-em)]2 molar ratio. The aim was to gain a deep understanding of the molecular properties of the different components in the solutions and the dynamic interactions between them. Especially, the impact of hydrogen bonding and π-stacking interactions between rings planes in drug⋯drug interactions leading to unwanted aggregations was carefully evaluated using combined angular/radial distribution functions and also natural bond orbital theory. Addition of the drugs to the DES led to decreased conductivity and diffusion coefficients (ranging from 1.52 to 1.07 × 10−12 m2 s−1). The effects of the small allopurinol on the DES structure are more pronounced compared with the larger losartan and omeprazole molecules. This is in line with the formation of aggregates, which are markedly smaller for allopurinol (∼2.6 molecules in average) than for omeprazole (∼4.3) and losartan (∼5.5). DFT-calculated dimers of the drug molecules are stabilized by π-stacking with energies ranging from −10 (allopurinol) to −32 kcal mol−1 (losartan) and show interplanar distances ranging from 0.36 to 0.47 nm. Our observations might pave the way for designing optimized DES for drug delivery.


1. Introduction

The term drug delivery systems (DDS) summarizes chemical vehicles by which drugs are transported and delivered to target cells, tissues, or organs.1–3 The common goal of DDS is to promote the application of therapeutic drugs by overcoming obstacles such as lack of selectivity and limited solubility. Limited solubility leads to low bioavailability and poor biodistribution. On a molecular level, low solubility leads to drug aggregation and thus reduced activity.1 Thus, solubility is one of the most important parameters for achieving the desired drug concentration in the target locus.3

Aromatic and heteroaromatic moieties are frequently found in drugs and an estimate says that oral drugs have an average of 1.6 aromatic rings in their structures.4 Detailed analyses have emphasized that rising the number of aromatic rings in their molecular structure of drugs lowers their overall developability/applicability, by decreasing their solubility in water and raising their lipophilicity. However, an important factor of aromatic and heteroaromatic structures in drug discovery is the observed enhancement of compound potency by increasing the ligand⋯receptor binding energy through reducing the entropy term.3

The biopharmaceutical classification system (BCS)5,6 sorts drugs along their solubility and their permeability (Fig. 1) and for drug development both aspects must be considered. The use of co-solvents to increase the solubility of drugs with inherently low water solubility has been previously employed including ionic liquids (IL)7 and deep eutectic solvents (DES).8–11 Various use of DESs for therapeutic applications has been reported in the last two decades.10–13 including the solubilization of heteroaromatic drugs such as the nonsteroidal anti-inflammatory drugs celecoxib or piroxicam by the use of choline-based DESs.11–14 Very recently, the mechanism of solubility improvement of some drugs belonging to BCS class II was studied in natural DESs based on solid organic acids, proline, choline chloride, and sugars, which revealed the formation of drug-DES supramolecular assemblies through hydrogen bonding (H bonding).15 Therapeutic DESs (TDES) and volatile DESs (VODES) have recently been explored as novel approaches for enhancing the solubility of active pharmaceutical ingredients (APIs) and controlling APIs crystal growth in a scalable way at room temperature, with the additional possibility of control over polymorphism.2,16


image file: d5cp01312g-f1.tif
Fig. 1 Biopharmaceutical Classification System of drugs (BCS).

Molecular dynamics (MD) simulations have evolved as a powerful tool to describe the complex behavior of the components in such DES.12,17–26 An early MD study was devoted to the structure and dynamics of a 1[thin space (1/6-em)]:[thin space (1/6-em)]2 eutectic mixture of choline chloride (ChCl) and ethylene glycol.23 A first-principles MD study on reline and its mixture with water was followed in 201825 and recently refined.19 Very similar to these studies, an MD study on dimethylsulfoxide (DMSO) as an aprotic solvent and reline as DES showed that DMSO has only minor effects on the intermolecular interactions between the reline components in comparison to water.24 Furthermore, while water molecules preferentially interacted with chloride ions, the main interactions of DMSO were shown to be with the urea molecules.24 In a series of ab initio molecular dynamics (AIMD) studies, the solvation shell structures of ammonia, aromatic and aliphatic compounds, and also CO2 and SO2 in DES were reported recently.20–22

The dissolution of heteroaromatic drugs in DES was the topic of a theoretical study on the dissolution of lidocaine (Scheme 1) in two selected DES.26 Very recently, the solubility of celecoxib (Scheme 1) in two ChCl-based DES containing ethylene glycol and 1,2-propanediol hydrogen bond (H bond) donors (HBD) was investigated using MD simulations showing that attractive interactions between the drug and the DES were prevalent over water⋯celecoxib interactions.12


image file: d5cp01312g-s1.tif
Scheme 1 Schematic structures of choline chloride (ChCl), lidocaine, celecoxib, losartan, omeprazole, and allopurinol.

We selected losartan, omeprazole, and allopurinol (Scheme 1), all containing heteroaromatic groups, from the list of the 100 most important drugs.27 Losartan is used as a medicine to treat high blood pressure and prevent heart damage and belongs to class III drugs in the BCS classification (Fig. 1).28 Omeprazole is a stomach acid inhibitor and is found in Class II of the BCS classification.6,15,29 Omeprazole has low stability in acidic solutions, poor water solubility, and undergoes first pass metabolism markedly reducing its concentration before reaching the systemic circulation. All of these factors result in low bioavailability.29,30 Allopurinol is effective in treating gout by reducing the amount of uric acid in the blood and belongs to class IV of BCS classification.31 Due to its low permeability and limited solubility, allopurinol has limited therapeutic effects when used through oral or rectal administration.32 A previous experimental study has shown that the aqueous solubility of allopurinol in various ChCl-based deep eutectic solvents (DESs) combined with urea, glycerol, citric acid, and glycol increases as the weight fraction of DESs increases.33 In contrast to a good number of experimental studies,8–10,13–15,17,33,34 insight on the molecular level mechanism of drug-DES interactions and the role of DESs in the solubility of various drugs through molecular dynamics (MD) studies remained scarce.12,18,19,26

We herein report on an MD study on the solvation properties of losartan, omeprazole, and allopurinol in the reline DES (ChCl/urea 1[thin space (1/6-em)]:[thin space (1/6-em)]2 mixture), where we observed some degree of self-association through π-stacking and further intermolecular interactions of the heteroaromatic structures.34–37 However, the aggregates are transient and limited in size, rather than persistent or large. Density functional theory (DFT) computational methods were also used in this study to gain insight into the hydrophobic π-stacking interactions. To our knowledge, this is the first study to provide a molecular-level picture of how solvation and aggregation mechanisms vary across different drug scaffolds in a DES medium. By integrating MD simulations with DFT, our work offers new insight into the subtle balance between solvation, self-aggregation, and drug⋯solvent interactions in DES, which is crucial for guiding rational drug formulation and delivery system design.

2. Computational methods

2.1. MD simulations

The structures of drugs and DES components were optimized at the B3LYP/6-31+G* level using Gaussian 09,38 and partial charges were computed through the electrostatic potential surface with the CHELPG procedure.39 DFT-calculated molecular structures of the DES components and the drug molecules are shown with their numbering in Fig. 2.
image file: d5cp01312g-f2.tif
Fig. 2 DFT-optimized structures of the DES components (a and b) and the heteroaromatic drugs (c–e) with atom labels.

The MD topology for each component was generated using the automated topology builder (ATB). MD modeling was performed with the GROMACS 5.1.2 program40 using the GROMOS96 53A6 force field, which was chosen for its reliability in capturing H bonding and solvation behavior in polar and H bond-rich environments,41–45 employed periodic boundary conditions in xyz directions with a short-range cut-off of 1.2 nm. The PME (Particle Mesh Ewald) procedure was utilized to model the electrostatic interactions of atom pairs further apart than 1.5 nm.46

A ChCl/urea 1[thin space (1/6-em)]:[thin space (1/6-em)]2 molar ratio was modeled using 1000 ion pairs of ChCl and 2000 urea molecules, randomly distributed in a low-density box. This system size was chosen to provide good statistical sampling while maintaining the correct eutectic ratio. The system was energy-minimized and equilibrated to allow the DES structure to form. After equilibration, the simulation box was enlarged to accommodate drug molecules. 12 drug molecules were then randomly introduced into the DES simulation box to construct a representative model system as a compromise between capturing statistically relevant intermolecular interactions and maintaining computational treatability. A second equilibration was performed to allow the drugs to interact with the solvent before starting the production simulations. As previous studies have shown, only a small number of molecules must be treated within the simulation box to achieve relevant data.41–45,47–51

To equilibrate each system five steps were executed: (1) Primary energy minimization was carried out to eliminate close contacts. (2) An equilibration NVT (constant number of particles, volume, temperature) run was executed beyond 5 ns in 1 × 10−3 ps time steps, while the system was linked to a velocity-rescaling thermostat using 0.1 ps coupling constants.52 (3) Consequent configurations were utilized to launch the NPT (constant number of particles, pressure, temperature) runs for 5 ns using the Parrinello–Rahman barostat.53 (4) After the primary equipoise, the systems were captured as the initial configurations for the production simulations of 10 ns, with 2 × 10−3 ps time steps at 450 K, ensuring that the mobility of molecules is high enough for achieving a reliable equilibration. (5) Simulations were carried out for 20 ns at 298 K. Considering the simulation annealing steps in which the temperature of each system was decreased at intervals of 50 K from 400 K to 298 K, the simulation was performed for 500 ps under constant NPT conditions at each temperature. The overlap of RDF graphs considering various atom⋯atom interactions at different simulation times proves that the systems are well equilibrated (see Fig. S1–S3, ESI).42,48–51

To address potential biases from temperature scaling in the NVT/NPT ensembles, we performed additional simulations in the NVE (constant N, V, and E) ensemble for 20 ns. The NVE ensemble allows for energy conservation without external thermostats, providing a more direct evaluation of the system's dynamical properties.

Radial distribution function (RDF) analysis was carried out by means of the gmx rdf tool in Gromacs40 to provide insight into the probability distribution of each DES or drug molecule around the other molecule.

Combined radial/angular distribution function (CDF) plots are a potent analysis for the study of H bonding and CDF analyses were performed by combining radial and angular distribution functions by trajectory analyzer and visualizer (TRAVIS) package.54 With this method, both inter- and intramolecular interactions have been previously studied.55

The formation of H bonds between the DES components as well as in DES⋯drug and drug⋯drug interactions was calculated using the gmx hbond criterium through the GROMACS package.40 The lifetime of H bonds was evaluated through time autocorrelation functions (TCF), using the TRAVIS software.54 The size of aggregations of the drug molecules was estimated using gmx clustsize through the GROMACS package.40 A histogram of aggregate (cluster) sizes was obtained from the simulation trajectory files with a threshold of lower than 0.35 nm indicative for aggregates.

To examine the dynamics of the drugs in the DES, the center of mass (COM) mean-square displacement (MSD) dataset was obtained during the MD trajectories using gmx msd through the GROMACS program.40 Also, the diffusion coefficients Di were computed. In general, the average distance squared when an atom or molecule moves away from beginning point throughout the time interval τ (eqn (1)).

 
image file: d5cp01312g-t1.tif(1)
In which the position of the COM of ith atom or molecule at time t is represented by [r with combining right harpoon above (vector)]ci(t) and the average of ensemble throughout the time origins is identified by the angular brackets. The diffusion coefficients (Di), for each trajectory, were calculated from the slope of the linear domain of the MSD curves. Di was obtained from the Einstein relation (eqn (2)):
 
image file: d5cp01312g-t2.tif(2)
Cationic and anionic transference numbers (ti) for the choline cations and chloride anions were calculated using eqn (3):
 
image file: d5cp01312g-t3.tif(3)
wherein (+) and (−) are related to the cations and anions, respectively.

The molar electrical conductivity of the solutions was determined using the Nernst–Einstein (eqn (4)).

 
image file: d5cp01312g-t4.tif(4)
with zi = ionic charges, F = Faraday constant (96[thin space (1/6-em)]485 C mol−1), R = gas constant (8.314 J mol−1 K−1).

Viscosity calculations were performed for each system using the Einstein–Stokes equation (eqn (5)).

 
image file: d5cp01312g-t5.tif(5)
where kB is the Boltzmann constant and c is a constant. In higher viscous liquids, the value of c is equal to 4.

2.2. DFT calculations

The non-covalent drug⋯drug interactions, were studied through a quantum mechanical DFT approach. Monomeric and dimeric models of the studied drugs were calculated by using the Gaussian 09 software package38 and utilizing the ωB97XB functional with the 6−311++G(d,p) basis set. This functional has an empirical dispersion correction term and accurately predicts π–π stacking and van der Waals interactions.56–58Eqn (6) was used to obtain the interaction energies for dimers.59
 
Interaction energy = EDimer − 2EMonomer(6)
wherein, EDimer is the single point energy of the dimer and EMonomer is the single point energy of the monomeric drug molecule. Natural bond orbital (NBO) analysis60,61 was used to study the π-stacking interactions in more detail. Based on this analysis, second-order perturbation stabilization energies as electron donor–acceptor bond energies E(2) were calculated to gauge the intensity of orbital overlap between electron donors and electron acceptors using eqn (7).
 
image file: d5cp01312g-t6.tif(7)
wherein F(i,j) is the off-diagonal NBO Fock matrix element, qi is the donor orbital occupancy, and εi and εj are diagonal elements. The larger the value of the bond energy is, the stronger is the interaction.

3. Results and discussion

3.1. MD simulation of structures in reline and the DES/drug systems

To probe for the general solubility of the drugs in the reline DES, we first studied the final simulated snapshots (Fig. S4 in the ESI). They showed that ChCl and urea are thoroughly mixed in the absence of drug molecules and homogeneously distributed in the simulation box. For the drug containing systems, the drug molecules are partially aggregated in the DES media (Fig. 3). For the larger molecules, losartan and omeprazole (Fig. 3b and c), pronounced aggregations was found, while the smaller allopurinol showed a far lower degree of aggregation (Fig. 3a).
image file: d5cp01312g-f3.tif
Fig. 3 Snapshots of the equilibrated systems: (a) DES/allopurinol, (b) DES/losartan and (c) DES/omeprazole. Urea (red), choline (orange), and Cl (green). C, H, N, S, and O atoms of drug molecules are depicted as pink, white, dark blue, yellow and red spheres, respectively.

In general, the simulated snapshots as pictures of individual configurations or even averaged structures do not comprise any numerical details of full simulation length. To explore the structure of the solutions in detail and to examine how the drugs are dissolved in DES on a molecular level, we used radial distribution function (RDF), spatial distribution function (SDF), and combined radial/angular distribution function (CDF) analyses.

3.1.1. Pure DES. When computing the atom⋯atom RDFs between urea, Cl, and choline, we chose the H(Urea), N(Urea), O(Urea), O(Chol), H(Chol) and N(Chol), as reference sites (Fig. 4). The strongest correlation in the RDF plot is that of H(Chol)⋯Cl at 0.20 nm distance followed in intensity by H(Urea)⋯Cl at around 0.22 nm with a pronounced second component at around 0.35 nm. Only a fraction of intensity and a very sharp peak is found for the very short H(Chol)⋯O(Urea) interaction at 0.18 nm. Markedly broader peaks and thus higher dynamics was found for the H(Urea)⋯O(Chol) interactions at 0.21 nm and the H(Chol)⋯N(Urea) interactions at 0.35 nm. This is in line with the OH group of the choline cation maintaining H bonds of moderate electrostatic (0.15 to 0.22 nm) to weak electrostatic (>0.22 nm) character with the Cl ions and urea molecules.12,18,19,62,63
image file: d5cp01312g-f4.tif
Fig. 4 Atom⋯atom RDFs between DES components in the pure system.

These interactions align closely with findings from prior investigations into pure reline. In these studies, the essential role of reline in achieving the eutectic properties is attributed to the specific H bonding interactions between chloride anions and the hydroxyl H atom of the choline cations, as well as with the H atoms of urea.12,18,19,62,63

3.1.2. DES/allopurinol. When adding allopurinol to the DES, the H(Chol)⋯Cl and H(Chol)⋯O(Urea) interactions are still dominant, as in the pure DES system (Fig. 5a), but the H(Chol)⋯Cl interaction has massively gained in intensity compared with the H(Chol)⋯O(Urea) interaction. In addition, we observe a short H(Chol)⋯O1 interaction at 0.16 nm with an intensity in the range of the shortest H(Chol)⋯O(Urea) RDF. Medium intensity H(Urea)⋯O1 interactions are observed at a distance of 0.20 nm (Fig. 5b). Comparison with pure DES shows that the H(Urea)⋯Cl RDF exhibits only minor variations in peak height and position, indicating that allopurinol does not significantly disrupt the urea⋯chloride H bonding network.
image file: d5cp01312g-f5.tif
Fig. 5 Atom⋯atom RDFs for the DES/allopurinol system (a–d).

The H3 atoms of allopurinol show a very high RDF with the Cl anions slightly below 0.2 nm, while further, less intense interactions of H3 and H2 with O(Urea) and O(Chol) atoms cluster range between 0.18 (H3⋯O(Urea)) and 0.25 nm (H2⋯O(Urea) (Fig. 5c). While these drug⋯DES interactions are short-ranged, the drug⋯drug interactions O1, N1, and N4 around H2 and H3 are slightly longer but probably dominant as can be seen from their far higher RDFs (Fig. 5d and Fig. S2, ESI). The strongest peaks are found for H2⋯N3 and H3⋯N2 at distances above 0.35 nm. A third, less intense RDF is observed for H2⋯O1 at 0.40 nm. Some intense RDFs are found for C1⋯C1, C1⋯C5 and C2⋯C3 at 0.47 nm (Fig. S5, ESI), which lies in the typical range for intermolecular π-stacking interactions of either π–π or CH–π type.62,63 These RDFs exhibit lower intensity when compared with the prevalent drug⋯drug H bond interactions.

3.1.3. DES/losartan. In the presence of losartan, the RDFs between the DES components are very similar to the DES/allopurinol system in that H(Chol)⋯Cl and H(Chol)⋯O(Urea) are dominant (Fig. 6b). The strongest interaction of losartan (Fig. 6a) appears at a very short distance 0.17 nm and represents the H(Urea)⋯O1 (OH group of losartan) interaction. The H(Chol)⋯O1 interaction has similar length but is far weaker (Fig. 6b). Further important solvent-drug interactions are those of Cl, O(Chol), and O(Urea) with H1 and H6 (NH group of losartan) (Fig. S6, ESI). Strong interactions of the H1 atom with Cl2 and the triazine N2 atom (Fig. S6, ESI), as well as the very intense RDFs for H6⋯N5 (tetrazole) at ∼0.5 nm (Fig. 6c) can be interpreted as intermolecular H bonds. Finally, the C atoms of the losartan biphenyl moieties show interactions with the biphenyl rings and also with C20 of the tetrazole ring (Fig. 6d). The strongest is found for C11 at 0.47 nm in line with intermolecular π-interactions. As in the DES/allopurinol system, the DES⋯drug H bonding interactions are markedly stronger than the drug⋯drug H bonding interactions. Also similar in both systems, are the relatively strong π-interactions. They seem more pronounced for losartan than for allopurinol.
image file: d5cp01312g-f6.tif
Fig. 6 Atom⋯atom RDF plots for the DES/losartan system (a–d).
3.1.4. DES/omeprazole. The presence of the sulfinyl (SO) group in omeprazole is very prominent, with an intense RDF peak for the H(Chol)⋯O2 interaction, appearing at a shorter distance than the H(Chol)⋯Cl (0.16 vs. 0.20 nm) RDF peak, which is found also in all other systems (Fig. 7a).
image file: d5cp01312g-f7.tif
Fig. 7 Atom⋯atom RDF plots of the DES/omeprazole system (a–d).

Additionally, there are peaks of low probability for H(Chol)⋯O1 and H(Chol)⋯O3 at 0.18 nm. The interaction of omeprazole with urea is less important for the solvation, but the H(Urea)⋯O2, ⋯O3, and ⋯N3 interactions are all found around 0.18 nm and represent an extended H bond network (Fig. 7b). The acidic N–H1 maintains similar weak, but short interactions with Cl, N(Urea), and O(Urea) (Fig. S7, ESI). As for allopurinol and losartan, the drug⋯drug H bonding network in the DES/omeprazol system seems to be dominant over the DES⋯drug H bonds as can be interpreted from the short H1⋯O2 and H1⋯N2 interactions that comes along with a longer but more intense H1⋯N3 interaction (Fig. 7c). π-stacking interactions can be inferred from intense the C3⋯C6 RDF peak and further smaller peaks in the range 0.35 to 0.5 nm (Fig. 7d). The highest RDF approaches a level of intensity similar to that of the most potent drug⋯drug H bond.

The sharp rise in the first peak intensity in the H(chol)⋯Cl RDF plot in the presence of drug molecules (Fig. 4, 5a, 6b, and 7a) can be attributed to the reorganization of the solvent structure caused by drug⋯solvent interactions. The drugs seem to induce stronger H bonding or van der Waals interactions with the choline cation and chloride anions, leading to a denser local environment. This results in a higher first peak in the RDF, suggesting more frequent or stronger interactions between H(chol) and Cl. We found that the first peak height in the H(chol)⋯Cl RDF was lowest for losartan, slightly larger for allopurinol, and markedly larger for omeprazol. Thus, the size of the drug molecules is not decisive.

The intense and strong peaks in the RDFs for the H bonding between the H atoms of drug molecules and chloride align with previous studies, which revealed that H bonding between chloride and various H atoms in β-cyclodextrin, MeOH, EtOH and water occurs at 0.2 nm.62–66

More details of the H bonding were revealed by CDF analyses and in the following, an H bond is defined when the D–H⋯A angle is between 160 and 180° and the A⋯H distance is lower than 0.3 nm. Pronounced H bonds are found in the DES/allopurinol system between allopurinol H2, H3, and O1 and the Cl, O(Urea), H(Urea), and H(Chol) DES components (Fig. 8 and Fig. S8, ESI). The H3⋯Cl distance is found around 0.19 nm with a linear N3–H3⋯Cl (180°) H bond (Fig. 8a). The strongest drug⋯urea interaction is found for H3⋯O(Urea) with an N3–H3⋯O(Urea) angle of 178° (Fig. 8b). A little bit less pronounced drug⋯urea and drug⋯choline H bond interactions are found for H(Urea)⋯O1 (Fig. 8c) and H(Chol)⋯O1 (Fig. 8d) with distances less than 0.2 nm and angles of 175°. This is in complete agreement with the conclusions drawn from the RDF plots (see above).


image file: d5cp01312g-f8.tif
Fig. 8 CDF plots for DES/allopurinol system (a–d). Red relates to low probability and dark blue relates to high normalized probability.

The same is true for the DES/losartan system (Fig. S9, ESI). Here, the most pronounced interactions occur between the DES H atoms H(Chol)/H(Urea) and the N1/O1/N4 atoms of losartan. By far the highest probabilities were found for the O1 atom with H(Chol) and N(Urea) with the O(Chol)–H(Chol)⋯O1 angle at 178° and the N(Urea)–H(Urea)⋯O1 angle at 180°. In turn, for the acidic H1 and H6 atoms of losartan, the strongest H bonds are found with the chloride ion of the DES and O of the choline cation (H1⋯Cl and H6⋯Cl at 0.18 nm with angles ranging from 165 to 180°, and 0.17 nm with an 180° angle for H6⋯O(Chol)). Overall, these H bonds are similar in character to those of allopurinol. However, there is a tendency of higher probability for the corresponding N–H⋯Cl interaction and an additional very strong N–H⋯O(Chol) H bond for the losartan systems, which is unmatched in the DES/allopurinol system.

CDF plots of the DES/omeprazole system show strong O(Chol)–H(Chol)⋯O2 interactions with angles close to 180° and sub-0.2 nm distances, confirming the important role of the S[double bond, length as m-dash]O group for H bonding (Fig. S10, ESI). Less intense interactions were detected for N1–H1⋯O(Urea) and O(Chol)–H(Chol)⋯O3, both showing 178° angles and distances below 0.2 nm.

Previous studies on solvation of rather hydrophilic (ammonia)20 and hydrophobic (methane)67 solutes showed, that both are mainly surrounded by the Chol hydroxyl groups, with NH3 H bonds having specific geometric criteria (≤0.25 nm and H(Chol)-O(Chol)-N(NH3) angle of ≤30°). Most OH groups in Chol show a tangential orientation at 60° to 80° and are about 0.4 nm from the hydrophobic solute. In arginine-based DESs, anesthetic drug solubility studies revealed strong HBA (hydrogen bond acceptor) and HBD (hydrogen bond donor) interactions at 0.19 nm and 160°.68

Insight into the three-dimensional space came from SDF analysis (Fig. 9 and Fig. S11, S12, ESI). In the DES/allopurinol system, the spatial density of the choline cation is predominantly located near the H3 and N3 atoms of allopurinol. The chloride anions are exclusively grouped close to the allopurinol H3 atom. Urea approaches the drug molecule clearly at its O1 atom. Allopurinol molecules are grouped in a clamp-like fashion around the drug molecule with prevailing presence in well-defined regions above and below the double-ring plane. So, for allopurinol, all intermolecular interactions (H bond and π-stacking) are rather localized at specific parts of the molecule.


image file: d5cp01312g-f9.tif
Fig. 9 The SDF iso-surfaces plots identifying the distribution of the first nearest neighbor molecules of the DES components and drug: choline (blue), Cl (green), urea (red) and drug (yellow) for the DES/allopurinol system.

The SDF analysis of DES/losartan (Fig. S11, ESI) is less simple. Choline, Cl, and urea show several components of spatial density around the losartan molecule. For the chloride ion the major component is located very close to the OH group of losartan, in line with the RDF (Fig. S6, ESI) and CDF (Fig. S9a, ESI) plots. Urea SDF densities are vastly delocalized, which is in line with the plethora of H bonds of the carbonyl H-accepting and the NH2 H-donating units with various OH, NH and even CH groups of the drug, that can be assumed from the RDF and CDF plots. The SDF density of the drug found in close proximity to the alkyl chain and in the space far from the losartan molecule, can be interpreted as dimerization, leading to special closeness of the alkyl chains.

The SDF pictures of the DES/omeprazole system (Fig. S12, ESI) show several density components for the choline cation, the chloride, and urea. Interestingly, the urea density clamps the drug molecule between the S[double bond, length as m-dash]O and the MeOPh groups, while the density of omeprazole is essentially located in the binding pocket of the concavely arranged omeprazole molecule, which also points to marked dimerization.

These observations match an earlier study on NH3, showing urea and chloride near the NH3 H atoms and Chol near the N atom.20 Chloride ions in alcohol–reline mixtures primarily interact with the alcohol hydroxyl group H atoms.65 This is fully in line with the isosurface distribution of chloride in the current study. Arginine-based DESs around anesthetic drugs show localized HBA and HBD molecules with minimal neighboring drug distribution.68 In ChCl:lactic acid DES, Chol cations surround the lidocaine phenyl, NH, and CO groups, (Scheme 1) with less prominent chloride ions, and no significant drug self-clustering.26 Lidocaine, due to its amphiphilic nature, could potentially replace one of the DES components and form gel-like structures under certain conditions, as suggested by a previous study.69 While this was not the focus of the present study, the possibility of gel formation in DESs incorporating lidocaine represents an interesting future research direction. This could have implications for improving drug release properties and stability in DES-based drug delivery systems.

The formation of ChCl⋯urea, ChCl⋯drug, and urea⋯drug H bonds was calculated using the gmx hbond criterium through the GROMACS package (Fig. S13 and Table S1, ESI).40 In the presence of drug molecules, the number of H bonds between solvent components is only slightly reduced from 812 (per time frame) to 808 for omeprazole, to 795 for allopurinol and to 773 for losartan, while essentially urea⋯drug (27 to 38) and to a smaller extent ChCl⋯drug (2 to 16) H bonds are formed (Table S1, ESI). Lifetimes of these H bonds as determined through time autocorrelation functions (Fig. S14, ESI) range from 3.8 to 5.5 ms and slightly decrease along the series: choline > urea > chloride (Table S2, ESI). Omeprazole forms significantly more H bonds with both urea and ChCl, compared to allopurinol and losartan. Especially the number of H bonds with ChCl is far higher (16 compared with 2). This is probably due to the unique and very polar sulfinyl SO group, adding to the acidic NH function. The latter is also found in allopurinol and losartan.

To further study possible π-stacking interactions, the CDFs between different aromatic and heteroaromatic rings are explored by monitoring both the distance (d) and the angle (α) (Fig. 10 and Fig. S15, S16, ESI). As in previous studies,70–73 we defined π–π stacking along two criteria: (i) d < 0.5 nm and (ii) 0° < α < 40° or 150° < α < 180°.


image file: d5cp01312g-f10.tif
Fig. 10 CDFs (angular radial distribution functions) of omeprazole ring planes in six combinations. Here, colors from bright red to dark blue indicate the probability distribution of stacking angles. CoR and RN are the center of ring and ring normal, respectively, for the labeling of the rings, see Fig. 2.

For allopurinol, the highest CDFs lie in the range of 0.75 to 1 nm distance with 40° and ∼140° angles (Fig. S15, ESI). This is in line with the intense atom⋯atom RDF, observed for C2⋯C5 at 0.75 nm (Fig. S5, ESI). Only a very weak RN1 (RNi = ring normal of ring i)/CoR1⋯CoR1 (CoRi = center of ring i) interaction is observed at around <0.5 nm with angles around 50 and 130°. For losartan, very similarly to allopurinol, most of the interactions lie markedly above 0.5 nm, only the RN1 (ring normal, ring 1)/CoR1⋯CoR1 (center of ring 1), CoR1⋯CoR3 (center of ring 3), and CoR1⋯CoR4 (center of ring 4) interactions lies at around 0.47 nm with 0° < α < 40° or 150° < α < 180° (Fig. S16, ESI). The interaction of the RN2/CoR2⋯CoR4 occurs at 0.47 nm with 0° < α < 30° or 160° < α < 180° in line with the intense RDF C11⋯C20 of the tetrazole ring (R4) and biphenyl (R2) at 0.47 nm (Fig. 6d). For the biphenyl moiety, a strong interaction of the RN2/CoR2⋯CoR2 is observed. For omeprazole, well defined maxima slightly above 0.4 nm with angles of around 50 and 140° for the RN2/CoR2⋯CoR2 interaction fulfill roughly our π–π stacking criteria (Fig. 10). This is fully in line with the intense C6⋯C6 RDF peak (Fig. 7d). The same is true for the intense RN1/CoR1⋯CoR1 interaction observed at distances <45 nm and angles varying from 10 to 170°, which fully match the shorter RDF peak of C9⋯C11 at around 0.42 nm. Thus, it seems that omeprazole and losartan exhibit pronounced π–π stacking interactions in the DES solution, while allopurinol shows only weak π–π stacking.

The average size of aggregations of the drug molecules was estimated using gmx clustsize through the GROMACS program (Fig. S17, ESI).40 Production MD simulations were extended to 100 ns to ensure sufficient sampling of interactions, particularly those involving transient aggregation events. The average drug cluster sizes were 2.585 for allopurinol (136.112 g mol−1), 5.511 for losartan (422.91 g mol−1), and 4.255 for omeprazole (345.42 g mol−1). It should be noted that these values reflect aggregation behavior under the specific simulation conditions used in this study (12 drug molecules per box), and would be expected to vary with changes in drug concentration.

Overall, the aggregation numbers remain small throughout the simulations, suggesting that although transient self-association occurs, the extent and stability of the aggregation is limited. This highlights the efficiency of reline DES in solvating these drugs and preventing the formation of large or insoluble aggregates. The smaller aggregate size of allopurinol suggests that it remains more dispersed in solution, which can facilitate mass transfer, whereas larger aggregates of losartan and omeprazole may reduce overall mobility within the DES.

These findings are consistent with previous studies showing that the solvation of doxorubicin varies across solvents, affecting its aggregation size,42 and that tetracycline forms larger clusters in salt solution compared to EtOH due to hydrophobic interactions.74 Additionally, 1-ethyl-3-methyl imidazolium acetate IL effectively disrupts Dopamine aggregation, resulting in smaller aggregates and improved dissolution.75

Experimental studies have consistently demonstrated that DESs can support the formation of nanostructures through self-assembly of both surfactants and drug molecules confirmed by Tyndall scattering and dynamic light scattering (DLS). For instance, the self-aggregation behavior of various surfactants in reline-based DES has revealed a broad range of aggregate sizes (1–200 nm), with some micellar systems comprising as many as ∼60–75 surfactant molecules.76–80 In the context of pharmaceuticals, experimental studies have shown that depending on the nature of the drug and the DES system, drug molecules may either remain molecularly dissolved or form colloidal aggregates.80 Further, it was reported that several poorly water-soluble drugs form colloidal dispersions (10–200 nm) in thymol-based DES media.81 In contrast, other drugs did not exhibit aggregation and remained truly dissolved.81 Our MD simulations, though limited to small cluster sizes considering a model system, offer insight into the molecular-level interactions that may initiate such behavior, and capture the initial stages of such aggregation phenomena.

3.2. Dynamic properties of DES and DES/drug systems

The diffusion coefficients of all components obtained from MSD diagrams listed in Table 1, show increasing dynamics for the DES components in the order Cl < urea < choline in both the pure DES and the DES/drug systems and generally lower dynamics for the drug molecules compare with those of the DES components. The highest diffusion coefficient of 0.48 × 10−12 m2 s−1 was found for allopurinol. Moreover, the diffusion coefficients of the DES components decrease in the presence of drug molecules, following the trend: DES > DES/allopurinol > DES/losartan > DES/omeprazole. These trends suggest that an increasing size and degree of aggregation of the drug molecules decreases the overall dynamics of the DES components, resulting in slower diffusion. This could also affect the mass transfer rates within the system and potentially the solubility and bioavailability of the drugs. The diffusion coefficients represent the drug mobility in DES, highlighting the influence of molecular size and interactions on solubilization and transport.9,10,12,15,17
Table 1 Diffusion coefficients (× 10−12 m2 s−1), electrical conductivity Λ (10−6 S m2 mol−1), and transference numbers (t) calculated from NPT simulationsa
System D chol+ D Cl D urea D drug t chol+ t Cl Λ
a Diffusion coefficients obtained from the MSD-time curves. Molar electrical conductivity Λ determined using the Nernst–Einstein equation (eqn (4)). t = transference (or transport) numbers outlined according to the anionic and cationic diffusion coefficient values image file: d5cp01312g-t7.tif.
DES 1.52 ± 0.15 0.76 ± 0.02 0.93 ± 0.11 0.67 0.33 8.57 ± 0.23
DES/allopurinol 1.14 ± 0.05 0.61 ± 0.06 0.80 ± 0.09 0.48 ± 0.06 0.65 0.35 6.58 ± 0.14
DES/losartan 1.12 ± 0.13 0.58 ± 0.04 0.76 ± 0.12 0.28 ± 0.03 0.66 0.34 6.39 ± 0.09
DES/omeprazole 1.07 ± 0.08 0.55 ± 0.07 0.74 ± 0.04 0.21 ± 0.01 0.66 0.34 6.09 ± 0.12


Previously reported experimental values for the diffusion coefficients of the studied drugs across different media are as follows: allopurinol in neutral aqueous medium (930.7 × 10−12 m2 s−1),82 omeprazole in MeOH (188.0 × 10−12 m2 s−1)83 and an aqueous buffer supporting electrolytes (231.0 × 10−12 m2 s−1),84 and losartan in aqueous solution (260.0 × 10−12 m2 s−1), DMSO (160.0 × 10−12 m2 s−1) and in a micellar environment (43.0 × 10−12 m2 s−1).85 The lower diffusion coefficients observed in our simulations are expected, given the significantly higher viscosity and complex H bonding network of the reline DES compared to typical organic solvents. Additionally, strong drug⋯solvent and drug⋯drug interactions in DES, such as H bonding and π-stacking, can further hinder molecular mobility. Nevertheless, the observed trends remain qualitatively consistent and useful for comparing relative mobility among the drugs in the DES environment.

To assess the reliability of our results, we compared the dynamical properties obtained from NPT (Table 1) and NVE (Table S3, ESI) simulations. The NVE simulations further confirmed the accuracy of the observed dynamic properties, demonstrating that temperature fluctuations do not significantly alter the diffusion coefficients or other transport parameters. This reinforces the reliability of our results and the validity of the NPT ensemble choice in capturing the key behaviors of the studied systems. The differences in the diffusion coefficients and other dynamic properties were minimal, suggesting that our original findings are robust and unaffected by the choice of ensemble.

As shown in Table 1, the values of cationic transport numbers (ti) for the choline cations are larger than 0.5 both in the presence or absence of the drug molecules in the DES. This means that the cation plays a major role in carrying electric current in these systems.

The calculated molar electrical conductivity of the solutions showed that addition of drug molecules to the DES decreases the mobility of ionic species, which negatively impacts the ionic conductivity of the mixtures. The dissolution of drug molecules lowered the electrical conductivity of the solutions along the series: pure DES > DES/allopurinol > DES/losartan > DES/omeprazole. Allopurinol has a relatively minor effect on the ionic mobility of the DES components due to its smaller size. Losartan and omeprazole show larger aggregation, which possibly leads to a stronger disruption of the ion movement in the solution, and reducing the free ions available for conductivity.86 The decrease in conductivity upon drug addition, further induces alterations in the DES structure, suggesting potential effects on the mass transfer of the solutes.9

Reline is known for its high viscosity, which can significantly impact the drug solubility and diffusion dynamics. Hence, viscosity calculations were performed for each system using the Einstein–Stokes equation (eqn (5)). The calculated viscosity for the pure DES of 1037 mPa s at 298.15 K, lies within the range of experimental values (750–1750 mPa s), which validates the computational model used in this study.87–90

Upon the addition of the drugs, the calculated viscosity increased to 1277 mPa s for allopurinol, 1336 mPa s for losartan, and 1395 mPa s for omeprazole. The larger the aggregates, the more resistance there is to the molecular mobility, leading to higher viscosity values. The dissolution of drug molecules also likely affects the solvent structure of the DES, leading to a more ordered or structured arrangement of the DES components around the drug molecules. This structuring effect can also increase viscosity, as the solvent components become less mobile. The increased viscosity upon the addition of the drugs, particularly the larger drugs, implies a more “viscous” environment that could slow down the mass transfer process.91–93

3.3. DFT study on the intermolecular π-stacking interactions

The stacking interaction for heterocyclic dimers with different conformations was modelled through DFT-geometry optimization of dimers in the gas phase (Fig. 11).
image file: d5cp01312g-f11.tif
Fig. 11 DFT-calculated geometry-optimized (ωB97XD/6-311++G(d,p)) dimer structures for (a) allopurinol, (b) losartan, and (c) omeprazole and their interaction energies in kcal mol−1. Atoms are represented as gray (C), red (O), dark blue (N), white (H), green (Cl), and yellow (S) balls.

For the bicyclic allopurinol structure, which comprises a pyrazole ring fused to a hydroxy-substituted pyrimidine ring, we observed an intermolecular rotation (∼120°) in the stacked dimer with an energy of about −10 kcal mol−1 compared with the monomer. The interplanar distance was calculated to be around 0.36 nm. The shortest CoR1⋯CoR2 interaction is found at around 0.45 nm in the DES, based on the RDF and CDF plots (Fig. S5 and S15, ESI). In the dimer model for losartan, the molecules were stacked through the biphenyl rings resulting in a markedly higher interaction energy of about −32 kcal mol−1 compared with that of allopurinol. The interplanar distances were calculated to lie between 0.36 and 0.47 nm, in good agreement with RDF and CDF interacting distances of 0.47 nm in the DES (Fig. 6 and Fig. S13, ESI). For the omeprazole dimer, both the substituted pyridine and the substituted benzimidazole rings are stacking with each other, the energy gain being about −26 kcal mol−1. The interplanar distances ranging between 0.34 and 0.42 nm, closely match with the primary π–π stacking interactions depicted in Fig. 7d and 10. It is important to note, that the dimer interaction energies were computed in the gas phase, and thus do not account for solvent effects present in the reline DES environment. As a result, the reported values likely represent upper-bound estimates of interaction strength. These gas-phase calculations provide qualitative insight into the relative tendencies of the drugs to self-associate, which were further supported by MD simulations performed in the presence of solvent.

The formation of dimers is thus favorable in all three systems with an increasing energy gain along the series: allopurinol < omeprazole < losartan. Among the three drugs studied, losartan also showed the highest average aggregation number in MD simulations, which can be attributed to its larger molecular size, extended aromatic system, and hydrophobic character. The biphenyl system in losartan provides a large, conjugated, and rigid π-surface that favors π–π stacking, a key driving force for self-aggregation in non-aqueous environments like DES. This structural feature, combined with the presence of additional polar and H bonding groups, promotes stable dimer and small aggregate formation. Notably, an experimental study94 has also reported the formation of dimeric degradation products of losartan under acidic or neutral conditions, in which nucleophilic substitution reactions occur between functional groups of losartan molecules. While these degradation pathways differ from our observed non-covalent aggregates, they suggest that losartan has a predisposition for molecular self-association, especially under certain solvent conditions.

A previous DFT-AIM study found stacking interaction energies between drug fragments such as indole, benzofuran, and benzothiophene and the natural nucleobases guanine, adenine, thymine, and cytosine falling below 12 kcal mol−1, particularly those involving one or two rings.95 N-based organic fused heterocycle stacked dimeric systems like allopurinol have been reported to exhibit effective π–π stacking interactions in the gas phase, with interaction energies typically below 10 kcal mol−1.59 Calculations on imidazole stacking interactions, have shown inter-ring distances ranging from 0.38 to 0.43 nm,96 which aligns with the π–π stacking distances between the benzimidazole rings in omeprazole calculated by us.

To assess the physico-chemical properties of the dimer models, we compared the HOMO–LUMO energies (Table S4, ESI). The HOMO–LUMO gap, representing the overall reactivity, decreased in the order: allopurinol < omeprazole < losartan. This is in line with the highest energy-gain through dimerization of losartan within this series. The value of global hardness for the losartan dimer (1.487 eV) compared with allopurinol (1.893 eV) and omeprazole (1.598 eV) agrees also with the idea that the negative dimerization energies go along with reduced reactivity. The losartan dimer structure showed the highest softness within the series and less resistance against changing in its electronic configuration.

The most important donor–acceptor interactions within the drug dimers obtained from NBO analysis are reported in Tables 2 to 4 as well as in Tables S5, and S6, ESI. π–π stacking interactions are found for all the dimers, whereas CH–π stacking interactions are observed for losartan and omeprazole (Tables S5 and S6, ESI). As expected, the overall strength of the π–interactions were considerably larger for losartan and omeprazole compared with allopurinol, in good agreement with the RDF and CDF results (Fig. 5, 6, 7, 8, 10 and Fig. S5 to S10, ESI) from the MD study.

Table 2 Electron donor and acceptor orbitals and second-order interaction energy E(2) for the allopurinol dimera
Dimer structure

image file: d5cp01312g-u1.tif

π–π stacking interactions
Donor Acceptor E (2) (kcal mol−1)
a BD(2) and BD*(2) are bonding and anti-bonding natural orbital of π and π* type, respectively. a and b refer to the two molecules as shown.
aBD(2) N1⋯C1 bBD*(2) C5⋯N4 0.36
bBD(2) C5⋯N4 aBD*(2) C2⋯O1 0.11
bBD(2) C3⋯C4 aBD*(2) C3⋯C4 0.14
bBD(2) N1⋯C1 aBD*(2) C5⋯N4 0.36
aBD(2) C5⋯N4 bBD*(2) C2⋯O1 0.11


Table 3 Electron donor and acceptor orbitals and second-order interaction energy E(2) for the losartan dimera
Dimer structure

image file: d5cp01312g-u2.tif

π–π stacking interactions
Donor Acceptor E (2) (kcal mol−1)
a BD(2) and BD*(2) are bonding and anti-bonding natural orbital of π and π* type, respectively. a and b refer to the two molecules as shown.
aBD(2) C12⋯C13 bBD*(2) C11⋯C12 2.29
aBD(2) C12⋯C13 bBD*(2) C9⋯C10 3.94
aBD(2) C12⋯C13 bBD*(2) C17⋯C18 1.00
aBD(2) C12⋯C13 bBD*(2) N4⋯N5 2.66
aBD(2) C12⋯C13 bBD*(2) N3⋯C20 1.50
aBD(2) C8⋯C9 bBD*(2) C11⋯C12 0.48
aBD(2) C10⋯C11 bBD*(2) C9⋯C10 0.13
bBD(2) C14⋯C19 aBD*(2) C12⋯C13 0.80
bBD(2) C14⋯C19 aBD*(2) C17⋯C18 0.89
aBD(2) C15⋯C16 bBD*(2) C17⋯C18 0.41
bBD(2) C14⋯C19 aBD*(2) C15⋯C16 1.06


Table 4 Electron donor and acceptor orbitals and second-order interaction energy E(2) for the omeprazole dimera
Dimer structure

image file: d5cp01312g-u3.tif

π–π stacking interactions
Donor Acceptor E (2) (kcal mol−1)
a BD(2) and BD*(2) are bonding and anti-bonding natural orbital of π and π* type, respectively. a and b refer to the two molecules as shown.
aBD(2) C11⋯C12 bBD*(2) C9⋯C10 0.19
aBD(2) C6⋯N3 bBD*(2) C5⋯C7 13.77
aBD(2) C6⋯N3 bBD*(2) C3⋯C4 43.62
aBD(2) C6⋯N3 bBD*(2) C2⋯C8 9.05
aBD(2) C5⋯C7 bBD*(2) C3⋯C4 4.42
aBD(2) C5⋯C7 bBD*(2) C2⋯C8 0.71
aBD(2) C3⋯C4 bBD*(2) C5⋯C7 2.37
aBD(2) C3⋯C4 bBD*(2) C2⋯C8 1.99
aBD(2) C3⋯C4 bBD*(2) C3⋯C4 28.61
aBD(2) C2⋯C8 bBD*(2) C5⋯C7 0.24
aBD(2) C2⋯C8 bBD*(2) C3⋯C4 1.19
aBD(2) C2⋯C8 bBD*(2) C2⋯C8 0.28


4. Conclusions

The solvation properties of the heteroaromatic drugs allopurinol, losartan, and omeprazole in the deep eutectic solvent (DES) reline, composed of choline chloride/urea in 1[thin space (1/6-em)]:[thin space (1/6-em)]2 molar ratio were studied combining molecular dynamics (MD) simulations and density functional theory (DFT) methods. Radial distribution function (RDF), combined radial/angular distribution function (CDF), and spatial distribution function (SDF) analyses showed multiple DES⋯DES, DES⋯drug, and drug⋯drug interactions. In all cases, the drug⋯drug interactions were slightly superior to DES⋯drug interactions, in line with aggregation observed already in the MD simulation snapshots. The average size of the aggregates was found to range between 2.6 and 5.5 molecules. In a real solution, aggregates of 2 to 5 can still be considered as “dissolved”. Intermolecular interactions, which are of H bonding and π-stacking character were studied in detail through RDF, CDF, and SDF analyses combined with calculation of the diffusion coefficients of the DES components and drugs, as well as through DFT-geometry optimizations for dimeric models of the three drugs. The dynamics for the DES components in each system is higher than that of the drugs. The calculated diffusion coefficients reach 0.48 × 10−12 m2 s−1 for allopurinol, 0.28 × 10−12 m2 s−1 for losartan, and 0.21 × 10−12 m2 s−1 for omeprazole, which match reasonably well experimental values. The π–stacking is especially important in losartan and omeprazole, in line with their lower coefficients. DFT-calculated dimers of the drug molecules are stabilized by π-stacking with energies ranging from −10 (allopurinol) to −32 kcal mol−1 (losartan) and show interplanar distances ranging from 0.36 to 0.47 nm, which lies in the typical range for heteroaromatic π–π-stacking. Overall, all three drugs show good solubility in the DES reline and we could show that reline maintains a plethora of drug⋯DES interactions improving the solubility of drugs such as losartan and omeprazole containing at the same time polar heteroaromatic units such as tetrazoles, purines, pyridines, and (benz)imidazoles, non-polar diphenyl groups, and also non-polar alkyl chains (losartan), methyl and methoxy substituents (omeprazole).

Conflicts of interest

The authors declare no competing interest.

Data availability

The data supporting this article have been included as part of the ESI.

Acknowledgements

The authors are indebted to the research council of the Shiraz University for the financial supports. We also thank “Iran Science Elites Federation” for partial support of this work. A. K. acknowledges the German Academic Exchange Service (DAAD) KD_0001052598-2 for a short-time guest lectureship at the Shiraz University.

References

  1. C. Li, J. Wang, Y. Wang, H. Gao, G. Wei, Y. Huang, H. Yu, Y. Gan, Y. Wang and L. Mei, Recent Progress in Drug Delivery, Acta Pharm. Sin. B, 2019, 9, 1145–1162 CrossRef.
  2. L. S. Taylor, D. E. Braun and J. W. Steed, Crystals and Crystallization in Drug Delivery Design, Mol. Pharmaceutics, 2021, 18, 751–753 CrossRef.
  3. C. A. S. Bergström, W. N. Charman and C. J. H. Porter, Computational Prediction of Formulation Strategies for Beyond-Rule-of-5 Compounds, Adv. Drug Delivery Rev., 2016, 101, 6–21 CrossRef.
  4. P. Ertl, Magic Rings: Navigation in the Ring Chemical Space Guided by the Bioactive Rings, J. Chem. Inf. Model., 2022, 62, 2164–2170 CrossRef CAS PubMed.
  5. J. M. Butler, The Application of Biopharmaceutics Classification Systems to Modified-Release Formulations, in Oral Drug Delivery for Modified Release Formulations, ed. E. S. Kostewicz, M. Vertzoni, H. A. E. Benson and M. S. Roberts, Wiley, 2022, vol. 1, pp. 103–117 Search PubMed.
  6. V. Gigante, G. M. Pauletti, S. Kopp, M. Xu, I. Gonzalez-Alvarez, V. Merino, M. P. McIntosh, A. Wessels, B.-J. Lee, K. R. Rezende, G. K. E. Scriba, G. P. S. Jadaun and M. Bermejo, Global Testing of a Consensus Solubility Assessment to Enhance Robustness of the WHO Biopharmaceutical Classification System, ADMET & DMPK, 2021, 9, 23–39 Search PubMed.
  7. C. Liu, B. Chen, W. Shi, W. Huang and H. Qian, Ionic Liquids for Enhanced Drug Delivery: Recent Progress and Prevailing Challenges, Mol. Pharmaceutics, 2022, 19, 1033–1046 CrossRef CAS PubMed.
  8. F. Oyoun, A. Toncheva, L. Castillo Henríquez, R. Grougnet, F. Laoutid, N. Mignet, K. Alhareth and Y. Corvis, Deep Eutectic Solvents: An Eco-friendly Design for Drug Engineering, ChemSusChem, 2023, 16, e202300669 Search PubMed.
  9. M. M. Abdelquader, S. Li, G. P. Andrews and D. S. Jones, Therapeutic deep eutectic solvents: A comprehensive review of their thermodynamics, microstructure and drug delivery applications, Eur. J. Pharm. Biopharm., 2023, 186, 85–104 Search PubMed.
  10. A. Sharma, Y. R. Park, A. Garg and B.-S. Lee, Deep Eutectic Solvents Enhancing Drug Solubility and Its Delivery, J. Med. Chem., 2024, 67, 14807–14819 CrossRef CAS PubMed.
  11. F. Santos and A. R. C. Duarte, Therapeutic Deep Eutectic Systems for the Enhancement of Drug Bioavailability, Deep Eutectic Solvents for Medicine, Gas Solubilization and Extraction of Natural Substances, ed. S. Fourmentin, M. C. Gomes, E. Lichtfouse, Springer International Publishing, Cham, 2021, vol. 56, pp. 103–129 Search PubMed.
  12. T. Amini, A. Haghtalab and J. Y. Seyf, Celecoxib solubility in two choline chloride-based deep eutectic solvents-Experimental study, thermodynamic modeling, and molecular dynamic simulation, Fluid Phase Equilib., 2023, 568, 113761 CrossRef CAS.
  13. M. Mokhtarpour, H. Shekaari, M. T. Zafarani-Moattar and S. Golgoun, Solubility and solvation behavior of some drugs in choline based deep eutectic solvents at different temperatures, J. Mol. Liq., 2020, 297, 111799 CrossRef CAS.
  14. H. Shekaari, M. T. Zafarani-Moattar, A. Shayanfar and M. Mokhtarpour, Effect of Choline Chloride/Ethylene Glycol or Glycerol as Deep Eutectic Solvents on the Solubility and Thermodynamic Properties of Acetaminophen, J. Mol. Liq., 2018, 249, 1222–1235 CrossRef CAS.
  15. B. Albertini, S. Bertoni, S. Sangiorgi, G. Nucci, N. Passerini and E. Mezzina, NaDES as a Green Technological Approach for the Solubility Improvement of BCS Class II APIs: An Insight into the Molecular Interactions, Int. J. Pharm., 2023, 634, 122696 CrossRef CAS PubMed.
  16. J. Potticary, C. Hall, V. Hamilton, J. F. McCabe and S. R. Hall, Crystallization from Volatile Deep Eutectic Solvents, Cryst. Growth Des., 2020, 20, 2877–2884 CrossRef CAS.
  17. Y. Vora, O. Sethi, S. N. Bariya, S. L. Gawali, S. S. Soni, T. S. Kang, P. A. Hassan and K. Kuperkar, Self-assembly of cationic surfactant in choline chloride-based deep eutectic solvents: structural solvation and dynamics, Phys. Chem. Chem. Phys., 2025, 27, 1171–1186 RSC.
  18. L. Yi, J. Wang, J. Liu, H. Luo, X. Wu and W.-Y. Li, Extraction mechanism of phenolic compounds by a choline chloride/glycerol solvent: DFT and molecular dynamics studies, Phys. Chem. Chem. Phys., 2024, 26, 29140–29149 RSC.
  19. H. Monteiro, A. Paiva, A. R. C. Duarte and N. Galamba, On the not so anomalous water-induced structural transformations of choline chloride–urea (reline) deep eutectic system, Phys. Chem. Chem. Phys., 2023, 25, 439–454 RSC.
  20. A. Malik and H. K. Kashyap, Solvation Shell Structures of Ammonia in Reline and Ethaline Deep Eutectic Solvents, J. Phys. Chem. B, 2023, 127, 2499–2510 CrossRef CAS PubMed.
  21. N. Kumar, P. K. Naik and T. Banerjee, Molecular Dynamic Insights into the Distinct Solvation Structures of Aromatic and Aliphatic Compounds in Monoethanolamine-Based Deep Eutectic Solvents, J. Phys. Chem. B, 2022, 126, 4925–4938 CrossRef CAS PubMed.
  22. A. Malik, H. S. Dhattarwal and H. K. Kashyap, Distinct Solvation Structures of CO2 and SO2 in Reline and Ethaline Deep Eutectic Solvents Revealed by AIMD Simulations, J. Phys. Chem. B, 2021, 125, 1852–1860 CrossRef CAS.
  23. Y. Zhang, D. Poe, L. Heroux, H. Squire, B. W. Doherty, Z. Long, M. Dadmun, B. Gurkan, M. E. Tuckerman and E. J. Maginn, Liquid Structure and Transport Properties of the Deep Eutectic Solvent Ethaline, J. Phys. Chem. B, 2020, 124, 5251–5264 CrossRef CAS.
  24. D. Shah, U. Mansurov and F. S. Mjalli, Intermolecular Interactions and Solvation Effects of Dimethylsulfoxide on Type III Deep Eutectic Solvents, Phys. Chem. Chem. Phys., 2019, 21, 17200–17208 RSC.
  25. E. O. Fetisov, D. B. Harwood, I.-F. W. Kuo, S. E. E. Warrag, M. C. Kroon, C. J. Peters and J. I. Siepmann, First-Principles Molecular Dynamics Study of a Deep Eutectic Solvent: Choline Chloride/Urea and its Mixture with Water, J. Phys. Chem. B, 2018, 122, 1245–1254 CrossRef CAS.
  26. A. Gutiérrez, M. Atilhan and S. Aparicio, A Theoretical Study on Lidocaine Solubility in Deep Eutectic Solvents, Phys. Chem. Chem. Phys., 2018, 20, 27464–27473 RSC.
  27. A. Hitchins, D. Lonsdale, D. Burrage and E. Baker, The Top 100 Drugs: Clinical Pharmacology and Practical Prescribing, 3rd edn, Elsevier, Amsterdam, 2022 Search PubMed.
  28. B. Pitt, R. Segal, F. A. Martinez, G. Meurers, A. J. Cowley, I. Thomas, P. C. Deedwania, D. E. Ney, D. B. Snavely and P. I. Chang, Randomised Trial of Losartan versus Captopril in Patients over 65 with Heart Failure Evaluation of Losartan in the Elderly Study, ELITE, Lancet, 1997, 349, 747–752 CrossRef CAS PubMed.
  29. M. G. Russo, M. I. Sancho, L. M. Silva, H. A. Baldoni, T. Venancio, J. Ellena and G. E. Narda, Looking for the Interactions Between Omeprazole and Amoxicillin in a Disordered Phase. An Experimental and Theoretical Study, Spectrochim. Acta, Part A, 2016, 156, 70–77 CrossRef CAS.
  30. S. S. Al-Nimry, K. A. Alkhamis and B. M. Altaani, Solid Self-Nanoemulsifying Drug Delivery System Filled in Enteric Coated Hard Gelatin Capsules for Enhancing Solubility and Stability of Omeprazole Hydrochloride, Pharm. Dev. Technol., 2020, 25, 588–600 CrossRef CAS PubMed.
  31. X.-L. Dai, J. Yao, C. Wu, J.-H. Deng, Y.-H. Mo, T.-B. Lu and J.-M. Chen, Solubility and Permeability Improvement of Allopurinol by Cocrystallization, Cryst. Growth Des., 2020, 20, 5160–5168 CrossRef CAS.
  32. S.-L. Chang, W. G. Kramer, S. Feldman, R. Ballentine and L. S. Frankel, Bioavailability of Allopurinol Oral and Rectal Dosage Forms, Am. J. Hosp. Pharm., 1981, 38, 365–368 CAS.
  33. W. Li, N. Wang, Z. Liu, Z. Fei, M. Zheng, J. Zhao and H. Zhao, Solubility Modelling and Thermodynamic Properties of Allopurinol in Aqueous Solutions of Four Deep Eutectic Solvents, J. Chem. Thermodyn., 2019, 132, 363–372 CrossRef CAS.
  34. S. K. Mann, T. N. Pham, L. L. McQueen, J. R. Lewandowski and S. P. Brown, Revealing Intermolecular Hydrogen Bonding Structure and Dynamics in a Deep Eutectic Pharmaceutical by Magic-Angle Spinning NMR Spectroscopy, Mol. Pharmaceutics, 2020, 17, 622–631 CrossRef CAS PubMed.
  35. P. D. Leeson, A. P. Bento, A. Gaulton, A. Hersey, E. J. Manners, C. J. Radoux and A. R. Leach, Target-Based Evaluation of “Drug-Like” Properties and Ligand Efficiencies, J. Med. Chem., 2021, 64, 7210–7230 CrossRef CAS PubMed.
  36. R. G. Huber, M. A. Margreiter, J. E. Fuchs, S. von Grafenstein, C. S. Tautermann, K. R. Liedl and T. Fox, Heteroaromatic π-Stacking Energy Landscapes, J. Chem. Inf. Model., 2014, 54, 1371–1379 CrossRef CAS PubMed.
  37. J. Graton, J.-Y. Le Questel, P. Maxwell and P. Popelier, Hydrogen-Bond Accepting Properties of New Heteroaromatic Ring, Chemical Motifs: A Theoretical Study, J. Chem. Inf. Model., 2016, 56, 322–334 CrossRef CAS.
  38. M. J. Frisch, G. W. Trucks, H. B. Schlegel, G. E. Scuseria, M. A. Robb, J. R. Cheeseman, G. Scalmani, V. Barone, B. Mennucci, G. A. Petersson, H. Nakatsuji, X. Li, M. Caricato, A. V. Marenich, J. Bloino, B. G. Janesko, R. Gomperts, H. P. Hratchian, J. V. Ortiz, A. F. Izmaylov, J. L. Sonnenberg, D. Williams-Young, F. Ding, F. Lipparini, F. Egidi, J. Goings, B. Peng, A. Petrone, T. Henderson, D. Ranasinghe, V. G. Zakrzewski, J. Gao, N. Rega, G. Zheng, W. Liang, M. Hada, M. Ehara, K. Toyota, R. Fukuda, J. Hasegawa, M. Ishida, T. Nakajima, Y. Honda, O. Kitao, H. Nakai, T. Vreven, K. Throssell, J. A. Montgomery, Jr., J. E. Peralta, F. Ogliaro, M. J. Bearpark, J. J. Heyd, E. N. Brothers, K. N. Kudin, V. N. Staroverov, T. A. Keith, R. Kobayashi, J. Normand, K. Raghavachari, A. P. Rendell, J. C. Burant, S. S. Iyengar, J. Tomasi, M. Cossi, J. M. Millam, M. Klene, C. Adamo, R. Cammi, J. W. Ochterski, R. L. Martin, K. Morokuma, O. Farkas, J. B. Foresman and D. J. Fox, Gaussian 09, Gaussian Inc., Wallingford, CT, USA, 2009 Search PubMed.
  39. C. M. Breneman and K. B. Wiberg, Determining Atom-Centered Monopoles from Molecular Electrostatic Potentials. The Need for High Sampling Density in Formamide Conformational Analysis, J. Comput. Chem., 1990, 11, 361–373 CrossRef CAS.
  40. S. Pronk, S. Páll, R. Schulz, P. Larsson, P. Bjelkmar, R. Apostolov, M. R. Shirts, J. C. Smith, P. M. Kasson, D. van der Spoel, B. Hess and E. Lindahl, GROMACS 4.5: a High-Throughput and Highly Parallel Open Source Molecular Simulation Toolkit, Bioinformatics, 2013, 29, 845–854 CrossRef CAS.
  41. H. Karimi, M. Heydari Dokoohaki, A. R. Zolghadr and M. H. Ghatee, The Interactions of an Aβ Protofibril with a Cholesterol-Enriched Membrane and Involvement of Neuroprotective Carbazolium-Based Substances, Phys. Chem. Chem. Phys., 2019, 21, 11066–11078 RSC.
  42. S. Shirazi-Fard, A. R. Zolghadr and A. Klein, How Does Aggregation of Doxorubicin Molecules Affect Its Solvation and Membrane Penetration?, New J. Chem., 2023, 47, 22063–22077 RSC.
  43. A. R. Zolghadr and M. Heydari Dokoohaki, How Does the P7C3-Series of Neuroprotective Small Molecules Prevent Membrane Disruption?, J. Chem. Inf. Model., 2017, 57, 2009–2019 CrossRef CAS.
  44. Z. Sheybani, M. Heydari Dokoohaki, M. Negahdaripour, M. Dehdashti, H. Zolghadr, M. Moghadami, S. M. Masoompour and A. R. Zolghadr, The Interactions of Folate with the Enzyme Furin: A Computational Study, RSC Adv., 2021, 11, 23815–23824 RSC.
  45. R. Ahmadi, B. Hemmateenejad, A. Safavi, Z. Shojaeifard, A. Shahsavar, A. Mohajeri, M. Heydari Dokoohaki and A. R. Zolghadr, Deep Eutectic–Water Binary Solvent Associations Investigated by Vibrational Spectroscopy and Chemometrics, Phys. Chem. Chem. Phys., 2018, 20, 18463–18473 RSC.
  46. T. Darden, D. York and L. Pedersen, Particle Mesh Ewald: An N logN Method for Ewald Sums in Large Systems, J. Chem. Phys., 1993, 98, 10089–10092 CrossRef CAS.
  47. A. R. Zolghadr and M. Heydari Dokoohaki, Self-Assembly of Neuroprotective Carbazolium Based Small Molecules at Octane/Water Interface: A Simulation Investigation, Chem. Phys., 2016, 480, 1–11 CrossRef CAS.
  48. M. Heydari Dokoohaki, A. R. Zolghadr and A. Klein, Impact of the Chemical Structure on the Distribution of Neuroprotective N-alkyl-9H-Carbazoles at Octanol/Water Interfaces, New J. Chem., 2020, 44, 1211–1220 RSC.
  49. A. R. Zolghadr and S. Boroomand, Spontaneous Assembly of HSP90 Inhibitors at Water/Octanol Interface: A Molecular Dynamics Simulation Study, Chem. Phys. Lett., 2017, 669, 130–136 CrossRef CAS.
  50. N. Azari, A. R. Zolghadr and M. Heydari Dokoohaki, The Use of 1-Ethyl-3-Methylimidazolium Iodide Ionic Liquid in Dye Sensitized Solar Cells: A Joint Experimental and Computational Perspective, J. Mol. Liq., 2022, 364, 119982 CrossRef.
  51. M. Heydari Dokoohaki, F. Mohammadpour and A. R. Zolghadr, Dye-Sensitized Solar Cells Based on Deep Eutectic Solvent Electrolytes: Insights from Experiment and Simulation, J. Phys. Chem. C, 2021, 125, 15155–15165 CrossRef CAS.
  52. M. Parrinello and A. Rahman, Polymorphic Transitions in Single Crystals: A New Molecular Dynamics Method, J. Appl. Phys., 1981, 52, 7182–7190 CrossRef CAS.
  53. G. Bussi, D. Donadio and M. Parrinello, Canonical Sampling Through Velocity Rescaling, J. Chem. Phys., 2007, 126, 014101 CrossRef PubMed.
  54. M. Brehm and B. Kirchner, TRAVIS-a Free Analyzer and Visualizer for Monte Carlo and Molecular Dynamics Trajectories, J. Chem. Inf., 2012, 4, 1–2 Search PubMed.
  55. R. Macchieraldo, L. Esser, R. Elfgen, P. Voepel, S. Zahn, B. M. Smarsly and B. Kirchner, Hydrophilic Ionic Liquid Mixtures of Weakly and Strongly Coordinating Anions with and without Water, ACS Omega, 2018, 3, 8567–8582 CrossRef CAS PubMed.
  56. J.-D. Chai and M. Head-Gordon, Long-Range Corrected Hybrid Density Functionals with Damped Atom-Atom Dispersion Corrections, Phys. Chem. Chem. Phys., 2008, 10, 6615–6620 RSC.
  57. L. M. da Costa, L. M. da Costa, S. R. Stoyanov, S. Gusarov, X. Tan, M. R. Gray, J. M. Stryker, R. R. Tykwinski, J. W. de, M. Carneiro, P. R. Seidl and A. Kovalenko, Density Functional Theory Investigation of the Contributions of π–π Stacking and Hydrogen-Bonding Interactions to the Aggregation of Model Asphaltene Compounds, Energy Fuels, 2012, 26, 2727–2735 CrossRef CAS.
  58. S. Grimme, Accurate Description of van der Waals Complexes by Density Functional Theory Including Empirical Corrections, J. Comput. Chem., 2004, 25, 1463–1473 CrossRef CAS PubMed.
  59. I. Ali, B. Ahmed, P. Barukial, M. J. Bezbaruah, M. Upadhyaya, A. Talukdar, K. Das and B. Bezbaruah, A Comparative DFT Study on the π–π Stacking Interaction of Some N-Based Organic Fused Heterocycles with Eclipsed and Staggered Conformations, J. Indian Chem. Soc., 2022, 99, 100698 CrossRef CAS.
  60. P. Bandyopadhyay, A. Karmakar, J. Deb, U. Sarkar and M. M. Seikh, Non-Covalent Interactions between Epinephrine and Nitroaromatic Compounds: A DFT Study, Spectrochim. Acta, Part A, 2020, 228, 117827 Search PubMed.
  61. A. E. Reed, L. A. Curtiss and F. Weinhold, Intermolecular Interactions from a Natural Bond Orbital, Donor-Acceptor Viewpoint, Chem. Rev., 1988, 88, 899–926 CrossRef CAS.
  62. P. Kumari, Shobhna, S. Kaur and H. K. Kashyap, Influence of Hydration on the Structure of Reline Deep Eutectic Solvent: A Molecular Dynamics Study, ACS Omega, 2018, 3, 15246–15255 CrossRef CAS PubMed.
  63. S. Kaur, S. Sharma and H. K. Kashyap, Bulk and Interfacial Structures of Reline Deep Eutectic Solvent: A Molecular Dynamics Study, J. Chem. Phys., 2017, 147, 194507 Search PubMed.
  64. A. Triolo, F. Lo Celso and O. Russina, Structural Features of β-cyclodextrin Solvation in the Deep Eutectic Solvent, Reline, J. Phys. Chem. B, 2020, 124, 2652–2660 CrossRef CAS PubMed.
  65. K. Kumar, A. Bharti and S. Mogurampelly, Insights on Choline Chloride–Based Deep Eutectic Solvent Reline + Primary Alcohol Mixtures: A Molecular Dynamics Simulation Study, J. Mol. Model., 2022, 28, 1–14 CrossRef.
  66. A. T. Celebi, T. J. Vlugt and O. A. Moultos, Structural, Thermodynamic, and Transport Properties of Aqueous Reline and Ethaline Solutions from Molecular Dynamics Simulations, J. Phys. Chem. B, 2019, 123, 11014–11025 CrossRef CAS PubMed.
  67. P. Nanavare, A. R. Choudhury, S. Sarkar, A. Maity and R. Chakrabarti, Structure and Orientation of Water and Choline Chloride Molecules around a Methane Hydrophobe: A Computer Simulation Study, Chem. Phys. Chem., 2022, 23, e202200446 CrossRef CAS PubMed.
  68. A. Gutiérrez, M. Atilhan and S. Aparicio, Theoretical Study on Deep Eutectic Solvents as Vehicles for the Delivery of Anesthetics, J. Phys. Chem. B, 2020, 124, 1794–1805 Search PubMed.
  69. J. Wu and T. Yin, Amphiphilic Deep Eutectic Solvent Based on Lidocaine and Lauric acid: Formation of Microemulsion and Gel, Langmuir, 2022, 38, 1170–1177 CrossRef CAS PubMed.
  70. C. R. Martinez and B. L. Iverson, Rethinking the Term ‘‘Pi-Stacking’, Chem. Sci., 2012, 3, 2191–2201 RSC.
  71. H. Weber and B. Kirchner, Complex Structural and Dynamical Interplay of Cyano-Based Ionic Liquids, J. Phys. Chem. B, 2016, 120, 2471–2483 Search PubMed.
  72. Y. Huang, Y. Ji, M. Zhang and D. Ouyang, How Imidazolium-Based Ionic Liquids Solubilize the Poorly Soluble Ibuprofen? A theoretical study, AIChE J., 2020, 66, e16940 CrossRef CAS.
  73. R. Thakuria, N. K. Nath and B. K. Saha, The Nature and Applications of π–π Interactions: A Perspective, Cryst. Growth Des., 2019, 19, 523–528 CrossRef CAS.
  74. M. Nabavi, M. R. Housaindokht, M. R. Bozorgmehr and A. Sadeghi, Self-Association Process of Tetracycline Antibiotic in Different Aqueous Solutions: A Joint Experimental Study and Molecular Dynamics Simulation, J. Iran. Chem. Soc., 2020, 17, 2997–3007 CrossRef CAS.
  75. A. Singh, T. G. Mason, Z. Lu, A. J. Hill, S. J. Pas, B. M. Teo, B. D. Freeman and E. I. Izgorodina, Structural Elucidation of Polydopamine Facilitated by Ionic Liquid Solvation, Phys. Chem. Chem. Phys., 2023, 25, 14700–14710 RSC.
  76. M. Pal, R. Rai, A. Yadav, R. Khanna, G. A. Baker and S. Pandey, Self-Aggregation of Sodium Dodecyl Sulfate within (Choline Chloride + Urea) Deep Eutectic Solvent, Langmuir, 2014, 30, 13191–13198 CrossRef CAS PubMed.
  77. A. Kumar, S. Singh and T. S. Kang, Deep Eutectic Solvents as Modulator on the Micellization Behaviour of Cationic Surfactants and Potential Application in Human Serum Albumin Aggregation, J. Mol. Liq., 2021, 344, 117864 CrossRef.
  78. R. Hakkinen, O. Alshammari, V. Timmermann, C. D’Agostino and A. Abbott, Nanoscale Clustering of Alcoholic Solutes in Deep Eutectic Solvents Studied by Nuclear Magnetic Resonance and Dynamic Light Scattering, ACS Sustain. Chem. Eng., 2019, 7, 15086–15092 CrossRef CAS.
  79. G. Komal, G. Singh, A. Singh and T. S. Kang, Aggregation Behavior of Sodium Dioctyl Sulfosuccinate in Deep Eutectic Solvents and Their Mixtures with Water: An Account of Solvent's Polarity, Cohesiveness, and Solvent Structure, ACS Omega, 2018, 3, 13387–13398 CrossRef CAS PubMed.
  80. F. Mazzotta, M. Curcio, L. Maiuolo and R. Cassano, Deep Eutectic Systems as Novel Vehicles for Assisting Drug Transdermal Delivery, Pharmaceutics, 2022, 14, 2265 CrossRef PubMed.
  81. J. Zhu, Y. Wei, J. Zhang, S. Qian, Y. Gao and W. Heng, Are All Poorly Soluble Drugs Dissolved in Deep Eutectic Solvents True Solutions?, J. Colloid Interface Sci., 2023, 645, 813–822 CrossRef CAS.
  82. V. Jiménez, J. B. Alderete, E. J. Delgado, J. Belmar and J. Gavín, On the Complexation of Allopurinol with β-cyclodextrin, Struct. Chem., 2006, 17, 217–223 CrossRef.
  83. F. Belal, N. El-Enany and M. Rizk, Anodic Polarographic Determination of Lansoprazole and Omeprazole in Pure Form and in Pharmaceutical Dosage Forms, J. Food Drug Anal., 2004, 12, 12 Search PubMed.
  84. S. M. A. Jorge, A. D. R. Pontinha and A. M. Oliveira-Brett, Electrochemical Redox Behavior of Omeprazole Using a Glassy Carbon Electrode, Electroanalysis, 2010, 22, 625–631 CrossRef CAS.
  85. M. Zervou, Z. Cournia, C. Potamitis, G. Patargias, S. Durdagi, S. G. Grdadolnik and T. Mavromoustakos, Insights into the Molecular Basis of Action of the AT1 Antagonist Losartan Using a Combined NMR Spectroscopy and Computational Approach, Biochim. Biophys. Acta, Biomembr., 2014, 1838, 1031–1046 CrossRef CAS PubMed.
  86. B. Kirchner, F. Malberg, D. S. Firaha and O. Hollóczki, Ion Pairing in Ionic Liquids, J. Phys.: Condens. Matter, 2015, 27, 463002 CrossRef PubMed.
  87. V. Agieienko and R. Buchner, Densities, Viscosities, and Electrical Conductivities of Pure Anhydrous Reline and Its Mixtures with Water in the Temperature Range 293.15 to 338.15 K, J. Chem. Eng. Data, 2019, 64, 4763–4774 CrossRef CAS.
  88. D. Lapeña, F. Bergua, L. Lomba, B. Giner and C. Lafuente, A Comprehensive Study of the Thermophysical Properties of Reline and Hydrated Reline, J. Mol. Liq., 2020, 303, 112679 CrossRef.
  89. H. Shekaari, M. T. Zafarani-Moattar and B. Mohammadi, Thermophysical Characterization of Aqueous Deep Eutectic Solvent Choline Chloride/Urea Solutions in Full Ranges of Concentration at T = 293.15–323.15 K, J. Mol. Liq., 2017, 243, 451–461 CrossRef CAS.
  90. D. Dhingra and S. Pandey, Effect of Lithium Chloride on the Density and Dynamic Viscosity of Choline Chloride/Urea Deep Eutectic Solvent in the Temperature Range 303.15–358.15 K, J. Chem. Thermodyn., 2019, 130, 166–172 CrossRef CAS.
  91. Z. Coin, S. Dangwal, M. H. Irwin, T. Knight, R. Bhave, G. Rother, R. Sacci, Md Arifuzzaman, V. Bocharova, I. Ivanov, S. S. Sholl and S. Z. Islam, Effect of Viscosity of a Deep Eutectic Solvent on CO2 Capture Performance in an Energy-Efficient Membrane Contactor-Based Process, ACS Omega, 2025, 10, 3407–3417 CrossRef CAS PubMed.
  92. C. Ferreira and M. Sarraguça, A Comprehensive Review on Deep Eutectic Solvents and Its Use to Extract Bioactive Compounds of Pharmaceutical Interest, Pharmaceuticals, 2024, 17, 124 CrossRef CAS.
  93. P. Makoś-Chełstowska, E. Słupek, S. Fourmentin and J. Gębicki, Supramolecular Deep Eutectic Solvents in Extraction Processes: A Review, Environ. Chem. Lett., 2025, 23, 41–65 CrossRef.
  94. Z. Zhao, Q. Wang, E. W. Tsai, X. Z. Qin and D. Ip, Identification of Losartan Degradates in Stressed Tablets by LC-MS and LC-MS/MS, J. Pharm. Biomed. Anal., 1999, 20, 129–136 CrossRef CAS PubMed.
  95. H. A. Toupkanloo and Z. Rahmani, An In-Depth Study on Noncovalent Stacking Interactions between DNA Bases and Aromatic Drug Fragments Using DFT Method and AIM Analysis: Conformers, Binding Energies, and Charge Transfer, Appl. Biol. Chem., 2018, 61, 209–226 CrossRef CAS.
  96. S. Ray and A. Das, Studies on the π–π Stacking Features of Imidazole Units Present in a Series of 5-Amino-1-Alkylimidazole-4-Carboxamides, J. Mol. Struct., 2015, 1089, 146–152 CrossRef CAS.

Footnote

Electronic supplementary information (ESI) available. See DOI: https://doi.org/10.1039/d5cp01312g

This journal is © the Owner Societies 2025
Click here to see how this site uses Cookies. View our privacy policy here.