Intermolecular Interactions of Edaravone in Aqueous Solutions of Ethaline and Glyceline Inferred from Experiments and Quantum Chemistry Computations

Edaravone, acting as a cerebral protective agent, is administered to treat acute brain infarction. Its poor solubility is addressed here by means of optimizing the composition of the aqueous choline chloride (ChCl)-based eutectic solvents prepared with ethylene glycol (EG) or glycerol (GL) in the three different designed solvents compositions. The slurry method was used for spectroscopic solubility determination in temperatures between 298.15 K and 313.15 K. Measurements confirmed that ethaline (ETA = ChCl:EG = 1:2) and glyceline (GLE = ChCl:GL = 1:2) are very effective solvents for edaravone. The solubility at 298.15 K in the optimal compositions was found to be equal x(E) = 0.158 (c(E) = 302.96 mg/mL) and x(E) = 0.105 (c(E) = 191.06 mg/mL) for glyceline and ethaline, respectively. In addition, it was documented that wetting of neat eutectic mixtures increases edaravone solubility which is a fortunate circumstance not only from the perspective of a solubility advantage but also addresses high hygroscopicity of eutectic mixtures. The aqueous mixture with 0.6 mole fraction of the optimal composition yielded solubility values at 298.15 K equal to x(E) = 0.193 (c(E) = 459.69 mg/mL) and x(E) = 0.145 (c(E) = 344.22 mg/mL) for glyceline and ethaline, respectively. Since GLE is a pharmaceutically acceptable solvent, it is possible to consider this as a potential new liquid form of this drug with a tunable dosage. In fact, the recommended amount of edaravone administered to patients can be easily achieved using the studied systems. The observed high solubility is interpreted in terms of intermolecular interactions computed using the Conductor-like Screening Model for Real Solvents (COSMO-RS) approach and corrected for accounting of electron correlation, zero-point vibrational energy and basis set superposition errors. Extensive conformational search allowed for identifying the most probable contacts, the thermodynamic and geometric features of which were collected and discussed. It was documented that edaravone can form stable dimers stabilized via stacking interactions between five-membered heterocyclic rings. In addition, edaravone can act as a hydrogen bond acceptor with all components of the studied systems with the highest affinities to ion pairs of ETA and GLE. Finally, the linear regression model was formulated, which can accurately estimate edaravone solubility utilizing molecular descriptors obtained from COSMO-RS computations. This enables the screening of new eutectic solvents for finding greener replacers of designed solvents. The theoretical analysis of tautomeric equilibria confirmed that keto-isomer edaravone is predominant in the bulk liquid phase of all considered deep eutectic solvents (DES).