# Experimental and Machine-Learning-Assisted Design of Pharmaceutically Acceptable Deep Eutectic Solvents for the Solubility Improvement of Non-Selective COX Inhibitors Ibuprofen and Ketoprofen

Piotr Cysewski, Tomasz Jeliński, Maciej Przybylek, Anna Mai, Julia Kulak

► To cite this version:

Piotr Cysewski, Tomasz Jeliński, Maciej Przybylek, Anna Mai, Julia Kulak. Experimental and Machine-Learning-Assisted Design of Pharmaceutically Acceptable Deep Eutectic Solvents for the Solubility Improvement of Non-Selective COX Inhibitors Ibuprofen and Ketoprofen. Molecules, 2024, 29 (10), pp.2296. 10.3390/molecules29102296. [hal-04575510](https://hal.science/hal-04575510)

HAL Id: hal-04575510

https://hal.science/hal-04575510

Submitted on 15 May 2024

HAL is a multi-disciplinary open access
archive for the deposit and dissemination of sci-
entific research documents, whether they are pub-
lished or not. The documents may come from
teaching and research institutions in France or
abroad, or from public or private research centers.

L'archive ouverte pluridisciplinaire HAL, est
destinée au dépôt et à la diffusion de documents
scientifiques de niveau recherche, publiés ou non,
émanant des établissements d'enseignement et de
recherche français ou étrangers, des laboratoires
publics ou privés.

<figure><img src="image_1.png" /></figure>

Distributed under a Creative Commons Attribution 4.0 International License
# Experimental and Machine-Learning-Assisted Design of Pharmaceutically Acceptable Deep Eutectic Solvents for the Solubility Improvement of Non-Selective COX Inhibitors Ibuprofen and Ketoprofen

Piotr Cysewski, Tomasz Jeliński, Maciej Przybylek, Anna Mai, Julia Kulak

<figure><img src="image_1.png" /></figure>

## To cite this version:

Piotr Cysewski, Tomasz Jeliński, Maciej Przybylek, Anna Mai, Julia Kulak. Experimental and Machine-Learning-Assisted Design of Pharmaceutically Acceptable Deep Eutectic Solvents for the Solubility Improvement of Non-Selective COX Inhibitors Ibuprofen and Ketoprofen. Molecules, 2024, 29 (10), pp.2296. 10.3390/molecules29102296. [hal-04575510](https://hal.science/hal-04575510v1)

HAL Id: hal-04575510

https://hal.science/hal-04575510v1

Submitted on 15 May 2024

HAL is a multi-disciplinary open access
archive for the deposit and dissemination of sci-
entific research documents, whether they are pub-
lished or not. The documents may come from
teaching and research institutions in France or
abroad, or from public or private research centers.

L'archive ouverte pluridisciplinaire HAL, est
destinée au dépôt et à la diffusion de documents
scientifiques de niveau recherche, publiés ou non,
émanant des établissements d'enseignement et de
recherche français ou étrangers, des laboratoires
publics ou privés.

<figure><img src="image_2.png" /></figure>
Article

# Experimental and Machine-Learning-Assisted Design of Pharmaceutically Acceptable Deep Eutectic Solvents for the Solubility Improvement of Non-Selective COX Inhibitors Ibuprofen and Ketoprofen

Piotr Cysewski* ID, Tomasz Jeliński ID, Maciej Przybyłek ID, Anna Mai and Julia Kułak

<figure><img src="image_2.png" /></figure>

check for
updates

**Citation:** Cysewski, P.; Jeliński, T.; Przybyłek, M.; Mai, A.; Kułak, J. Experimental and Machine-Learning-Assisted Design of Pharmaceutically Acceptable Deep Eutectic Solvents for the Solubility Improvement of Non-Selective COX Inhibitors Ibuprofen and Ketoprofen. *Molecules* **2024**, *29*, 2296. https://doi.org/10.3390/molecules29102296

Academic Editor: Reza Haghbakhsh

Received: 27 April 2024

Revised: 9 May 2024

Accepted: 12 May 2024

Published: 14 May 2024

<figure><img src="image_3.png" /></figure>

**Copyright:** © 2024 by the authors.
Licensee MDPI, Basel, Switzerland.
This article is an open access article
distributed under the terms and
conditions of the Creative Commons
Attribution (CC BY) license (https://
creativecommons.org/licenses/by/
4.0/).

Department of Physical Chemistry, Pharmacy Faculty, Collegium Medicum of Bydgoszcz, Nicolaus Copernicus University in Toruń, Kurpińskiego 5, 85-096 Bydgoszcz, Poland; tomasz.jelinski@cm.umk.pl (T.J.); m.przybyłek@cm.umk.pl (M.P.)

* Correspondence: piotr.cysewski@cm.umk.pl

**Abstract:** Deep eutectic solvents (DESs) are commonly used in pharmaceutical applications as excellent solubilizers of active substances. This study investigated the tuning of ibuprofen and ketoprofen solubility utilizing DESs containing choline chloride or betaine as hydrogen bond acceptors and various polyols (ethylene glycol, diethylene glycol, triethylene glycol, glycerol, 1,2-propanediol, 1,3-butanediol) as hydrogen bond donors. Experimental solubility data were collected for all DES systems. A machine learning model was developed using COSMO-RS molecular descriptors to predict solubility. All studied DESs exhibited a cosolvency effect, increasing drug solubility at modest concentrations of water. The model accurately predicted solubility for ibuprofen, ketoprofen, and related analogs (flurbiprofen, felbinac, phenylacetic acid, diphenylacetic acid). A machine learning approach utilizing COSMO-RS descriptors enables the rational design and solubility prediction of DES formulations for improved pharmaceutical applications.

**Keywords:** non-selective COX inhibitors; ibuprofen; ketoprofen; deep eutectic solvents; solubility; machine learning; COSMO-RS

## 1. Introduction

Ibuprofen (IP, dexibuprofen, CAS 51146-56-6) and ketoprofen (KP, dexketoprofen, CAS 22071-15-4) are nonsteroidal anti-inflammatory drugs (NSAIDs) [1,2] widely employed in clinical practice for their analgesic [3,4], antipyretic [5,6], and anti-inflammatory [7,8] properties. Chemically, both belong to the class of propionic acid derivatives. Ibuprofen, chemically known as (RS)-2-(4-(2-methylpropyl)phenyl)propanoic acid, and ketoprofen, 2-(3-benzoylphenyl)propionic acid, share structural similarities, differing primarily in their substituent groups. Their therapeutic efficacy stems from their ability to inhibit the cyclooxygenase (COX) enzymes, particularly COX-1 and COX-2, thereby impeding the synthesis of prostaglandins involved in inflammation, pain, and fever pathways [9–12]. The metabolism of ketoprofen and ibuprofen involves a series of enzymatic processes that primarily occur in the liver, leading to the formation of metabolites with varying pharmacological activities and elimination profiles. The two compounds undergo bio-transformation primarily via cytochrome P450 (CYP) enzymes [13–15], and the resulting metabolites include hydroxylated forms, such as 2-hydroxyibuprofen, 3-hydroxyibuprofen, 3-hydroxyketoprofen, and 4-hydroxyketoprofen. The two drugs find extensive application in the management of various conditions, including musculoskeletal disorders [16,17], rheumatic diseases [18,19], and prolonged conditions such as Alzheimer's [20], Parkinson's [21], and Machado-Joseph diseases [22]. However, their usage is not devoid of hazards. Gastrointestinal complications, such as ulceration and bleeding, constitute notable adverse effects associated with both drugs, albeit with varying incidences [23–25].
<figure><img src="image_3.png" /></figure>

Article

# Experimental and Machine-Learning-Assisted Design of Pharmaceutically Acceptable Deep Eutectic Solvents for the Solubility Improvement of Non-Selective COX Inhibitors Ibuprofen and Ketoprofen

Piotr Cysewski* ID, Tomasz Jeliński ID, Maciej Przybyłek ID, Anna Mai and Julia Kułak

Department of Physical Chemistry, Pharmacy Faculty, Collegium Medicum of Bydgoszcz, Nicolaus Copernicus University in Toruń, Kurpińskiego 5, 85-096 Bydgoszcz, Poland; tomasz.jelinski@cm.umk.pl (T.J.); m.przybyłek@cm.umk.pl (M.P.)

* Correspondence: piotr.cysewski@cm.umk.pl

**Abstract:** Deep eutectic solvents (DESs) are commonly used in pharmaceutical applications as excellent solubilizers of active substances. This study investigated the tuning of ibuprofen and ketoprofen solubility utilizing DESs containing choline chloride or betaine as hydrogen bond acceptors and various polyols (ethylene glycol, diethylene glycol, triethylene glycol, glycerol, 1,2-propanediol, 1,3-butanediol) as hydrogen bond donors. Experimental solubility data were collected for all DES systems. A machine learning model was developed using COSMO-RS molecular descriptors to predict solubility. All studied DESs exhibited a cosolvency effect, increasing drug solubility at modest concentrations of water. The model accurately predicted solubility for ibuprofen, ketoprofen, and related analogs (flurbiprofen, felbinac, phenylacetic acid, diphenylacetic acid). A machine learning approach utilizing COSMO-RS descriptors enables the rational design and solubility prediction of DES formulations for improved pharmaceutical applications.

**Keywords:** non-selective COX inhibitors; ibuprofen; ketoprofen; deep eutectic solvents; solubility; machine learning; COSMO-RS

<figure><img src="image_4.png" /></figure>

check for
updates

**Citation:** Cysewski, P.; Jeliński, T.; Przybyłek, M.; Mai, A.; Kułak, J. Experimental and Machine-Learning-Assisted Design of Pharmaceutically Acceptable Deep Eutectic Solvents for the Solubility Improvement of Non-Selective COX Inhibitors Ibuprofen and Ketoprofen. *Molecules* **2024**, *29*, 2296. https://doi.org/10.3390/molecules29102296

Academic Editor: Reza Haghbakhsh

Received: 27 April 2024

Revised: 9 May 2024

Accepted: 12 May 2024

Published: 14 May 2024

<figure><img src="image_5.png" /></figure>

**Copyright:** © 2024 by the authors. Licensee MDPI, Basel, Switzerland. This article is an open access article distributed under the terms and conditions of the Creative Commons Attribution (CC BY) license (https://creativecommons.org/licenses/by/4.0/).

## 1. Introduction

Ibuprofen (IP, dexibuprofen, CAS 51146-56-6) and ketoprofen (KP, dexketoprofen, CAS 22071-15-4) are nonsteroidal anti-inflammatory drugs (NSAIDs) [1,2] widely employed in clinical practice for their analgesic [3,4], antipyretic [5,6], and anti-inflammatory [7,8] properties. Chemically, both belong to the class of propionic acid derivatives. Ibuprofen, chemically known as (RS)-2-(4-(2-methylpropyl)phenyl)propanoic acid, and ketoprofen, 2-(3-benzoylphenyl)propionic acid, share structural similarities, differing primarily in their substituent groups. Their therapeutic efficacy stems from their ability to inhibit the cyclooxygenase (COX) enzymes, particularly COX-1 and COX-2, thereby impeding the synthesis of prostaglandins involved in inflammation, pain, and fever pathways [9–12]. The metabolism of ketoprofen and ibuprofen involves a series of enzymatic processes that primarily occur in the liver, leading to the formation of metabolites with varying pharmacological activities and elimination profiles. The two compounds undergo bio-transformation primarily via cytochrome P450 (CYP) enzymes [13–15], and the resulting metabolites include hydroxylated forms, such as 2-hydroxyibuprofen, 3-hydroxyibuprofen, 3-hydroxyketoprofen, and 4-hydroxyketoprofen. The two drugs find extensive application in the management of various conditions, including musculoskeletal disorders [16,17], rheumatic diseases [18,19], and prolonged conditions such as Alzheimer's [20], Parkinson's [21], and Machado-Joseph diseases [22]. However, their usage is not devoid of hazards. Gastrointestinal complications, such as ulceration and bleeding, constitute notable adverse effects associated with both drugs, albeit with varying incidences [23–25].
Additionally, the prolonged administration of high doses of these NSAIDs may precipitate
renal dysfunction [26] and cardiovascular events [27].

Solubility lies at the core of pharmaceutical science, exerting a profound influence on every stage of drug development and delivery [28–30]. Its importance resonates through the entire lifecycle of a medication, from initial formulation to therapeutic administration [31–33]. The solubility of a compound dictates its bioavailability, the rate and extent of absorption into the bloodstream, and consequently, its therapeutic efficacy [34,35]. Moreover, solubility profoundly impacts the stability of drug formulations, affecting shelf life and storage conditions [36,37]. Poor solubility poses significant challenges, necessitating higher doses to achieve therapeutic levels, which can increase the risk of adverse effects [38–40]. Therefore, pharmaceutical scientists dedicate considerable effort to understand [41,42], predict [43,44], and manipulate [45,46] the solubility of drug molecules. Techniques such as salt formation, cosolvent systems, cyclodextrin-based formulations, and nanoparticle delivery systems are employed to enhance solubility and thereby improve drug performance [47–50]. As it was said, solubility modeling and predictions play a pivotal role in the pharmaceutical industry by aiding in drug discovery and development processes [28,51–53]. Accurate predictions of solubility help researchers in selecting the most promising drug candidates, optimizing formulation strategies, and improving drug delivery systems. Among various predictive methods, artificial neural networks, with their ability to learn complex patterns from data, have emerged as powerful tools in solubility modeling and found widespread use in the pharmaceutical domain [54–57]. The integration of advanced computational methods like neural networks with traditional pharmaceutical research facilitates the design of more effective and safer medications.

The abovementioned need for safer pharmaceutical products, not only in terms of patient safety but also from an environmental perspective, led to a widespread application of the “green chemistry” framework [58,59]. In particular, the introduction of green solvents in the pharmaceutical industry could be regarded as a paradigm shift towards sustainability and environmental responsibility [60,61]. Green solvents, characterized by their low toxicity, renewable sourcing, and biodegradability, offer a compelling alternative to traditional organic solvents, which often pose significant environmental and health hazards [62,63]. Among the various techniques and chemicals used to meet the criteria of the concept of green chemistry, a notable group are deep eutectic solvents, abbreviated as DESs [64–68]. The majority of these systems are mixtures of two or more components, which include a hydrogen bond acceptor (HBA) and a hydrogen bond donor (HBD) that can be either solid or liquid. The melting point of the resulting mixture is lower than the melting points of its components, because of which they remain in liquid form even at room temperature [69,70]. DESs exhibit a unique combination of properties such as low toxicity, biodegradability, and importantly tunable physicochemical characteristics [71–73]. Despite the fact that the extent of their environmental friendliness has recently become a topic of debate [74–77], DESs showcase remarkable potential across various domains of pharmaceutical applications, including drug synthesis, formulation, and extraction processes [78–88]. Noteworthy, DESs have been shown to be effective solubilizers of ibuprofen [89,90] and ketoprofen [91]. As it can be seen, a substantial dataset has emerged, which serves as an intriguing starting point for modeling the solubilization properties of both examined active substances. Nevertheless, it is worth enriching this collection with new experimental results, especially the solubility values in DESs.

The aim of this paper is threefold. Firstly, the tuning of ibuprofen and ketoprofen solubility by designed DESs was undertaken. Then, machine learning was used for the formulation of the solubility model. Finally, predictions of selected analogues in DESs comprising choline chloride and betaine were provided.

**2. Results and Discussion**

The experimental part of this project involved the optimization of DES content and compositions at room temperature to find the most suitable solvent for ibuprofen and
Additionally, the prolonged administration of high doses of these NSAIDs may precipitate
renal dysfunction [26] and cardiovascular events [27].

Solubility lies at the core of pharmaceutical science, exerting a profound influence on every stage of drug development and delivery [28–30]. Its importance resonates through the entire lifecycle of a medication, from initial formulation to therapeutic administration [31–33]. The solubility of a compound dictates its bioavailability, the rate and extent of absorption into the bloodstream, and consequently, its therapeutic efficacy [34,35]. Moreover, solubility profoundly impacts the stability of drug formulations, affecting shelf life and storage conditions [36,37]. Poor solubility poses significant challenges, necessitating higher doses to achieve therapeutic levels, which can increase the risk of adverse effects [38–40]. Therefore, pharmaceutical scientists dedicate considerable effort to understand [41,42], predict [43,44], and manipulate [45,46] the solubility of drug molecules. Techniques such as salt formation, cosolvent systems, cyclodextrin-based formulations, and nanoparticle delivery systems are employed to enhance solubility and thereby improve drug performance [47–50]. As it was said, solubility modeling and predictions play a pivotal role in the pharmaceutical industry by aiding in drug discovery and development processes [28,51–53]. Accurate predictions of solubility help researchers in selecting the most promising drug candidates, optimizing formulation strategies, and improving drug delivery systems. Among various predictive methods, artificial neural networks, with their ability to learn complex patterns from data, have emerged as powerful tools in solubility modeling and found widespread use in the pharmaceutical domain [54–57]. The integration of advanced computational methods like neural networks with traditional pharmaceutical research facilitates the design of more effective and safer medications.

The abovementioned need for safer pharmaceutical products, not only in terms of patient safety but also from an environmental perspective, led to a widespread application of the “green chemistry” framework [58,59]. In particular, the introduction of green solvents in the pharmaceutical industry could be regarded as a paradigm shift towards sustainability and environmental responsibility [60,61]. Green solvents, characterized by their low toxicity, renewable sourcing, and biodegradability, offer a compelling alternative to traditional organic solvents, which often pose significant environmental and health hazards [62,63]. Among the various techniques and chemicals used to meet the criteria of the concept of green chemistry, a notable group are deep eutectic solvents, abbreviated as DESs [64–68]. The majority of these systems are mixtures of two or more components, which include a hydrogen bond acceptor (HBA) and a hydrogen bond donor (HBD) that can be either solid or liquid. The melting point of the resulting mixture is lower than the melting points of its components, because of which they remain in liquid form even at room temperature [69,70]. DESs exhibit a unique combination of properties such as low toxicity, biodegradability, and importantly tunable physicochemical characteristics [71–73]. Despite the fact that the extent of their environmental friendliness has recently become a topic of debate [74–77], DESs showcase remarkable potential across various domains of pharmaceutical applications, including drug synthesis, formulation, and extraction processes [78–88]. Noteworthy, DESs have been shown to be effective solubilizers of ibuprofen [89,90] and ketoprofen [91]. As it can be seen, a substantial dataset has emerged, which serves as an intriguing starting point for modeling the solubilization properties of both examined active substances. Nevertheless, it is worth enriching this collection with new experimental results, especially the solubility values in DESs.

The aim of this paper is threefold. Firstly, the tuning of ibuprofen and ketoprofen solubility by designed DESs was undertaken. Then, machine learning was used for the formulation of the solubility model. Finally, predictions of selected analogues in DESs comprising choline chloride and betaine were provided.

**2. Results and Discussion**

The experimental part of this project involved the optimization of DES content and compositions at room temperature to find the most suitable solvent for ibuprofen and
ketoprofen. Then, the most effective solvent systems were systematically studied at dif-
ferent temperatures and by dilution with water. This aspect of adding water to the DES
often leads to a surprising structure [92,93] and solubilizing properties [79,83]. Water is
considered an antisolvent because the aqueous solubility of various compounds is sev-
eral orders of magnitude lower compared to DESs, and dilution with water in excess
amounts significantly reduces the solubility of many solutes in such DES-water mixtures.
However, there are narrow ranges of water content that promote solubility compared to
pure DES. Nonetheless, due to the complicated and complex structure of aqueous DES
mixtures, water cannot be considered a cosolvent in such ranges. Using the example of
deep eutectic choline chloride/urea/water solvent mixtures [92], it was shown that the
DES nanostructure is formulated and preserved in a wide hydration range up to a remark-
ably high water content. This is attributed to the solvophobic sequestration of water in
nanostructured domains around cholinium cations. At higher water concentrations, such
separation becomes unfavorable, and the DES nanostructure is disrupted. That is why
DESs are often highly hygroscopic, and the control of the water content is so important. If
excess amounts of water are added, the water-water and DES-water interactions become
dominant, and only in such a regime can DES-water mixtures be considered as aqueous
solutions of DES components. The presence of these nanostructures might be regarded
as the structural origin of the observed increased solubility of the dissolved substances
compared to dry DES.

2.1. *Solubility of IP and KP in DESs*

In the initial phase of the experiments, a total of 36 eutectic systems were studied for both ibuprofen (IP) and ketoprofen (KP). One arrives at this number by combining two hydrogen bond acceptors (HBAs), namely, choline chloride (ChCl) and betaine (BI), with six hydrogen bond donors (HBDs), namely ethylene glycol (ETG), diethylene glycol (DEG), triethylene glycol (TEG), 1,3-butanediol (B3D), 1,2-propanediol (P2D), and glycerol (GLY), as well as three molar HBA:HBD proportions, namely 1:1, 1:2, and 1:4. This phase was aimed not only at determining the solubility of IP and KP in neat DES but also at selecting the best systems for next-stage studies involving aqueous mixtures. The composition screening was conducted at 25 °C.

In the case of ibuprofen, the following observations could be made. First of all, the
eutectic systems involving choline chloride were slightly more effective than those utilizing
betaine. Secondly, when comparing the effectiveness of HBD components, TEG turned
out to be the most effective in almost all cases, the only exception being the ChCl-GLY
among the 1:2 systems. Finally, when comparing molar ratios, the 1:2 ratio was the most
efficient, with the 1:4 ratio coming second and followed by the 1:1 ratio. For both choline
chloride and betaine in the 1:2 molar ratio, the following decreasing order of solubilizing
potential was noticed: TEG > GLY > DEG > ETG > B3D > P2D. This order varies only
slightly in other proportions, with GLY being the most effective in the case of 1:1 molar
ratio, ETG surpassing DEG in the 1:4 molar proportion, and P2D outperforming B3D on
two occasions. The actual solubility of IP at 25 °C in the ChCl-TEG (1:2) system was found
to be xIP = 0.353, while the ChCl-GLY (1:2) system amounted to a solubility of xIP = 0.349.
In the case of betaine, these values were xIP = 0.321 for BI-TEG (1:2) and xIP = 0.297 for
BI-GLY (1:2). Detailed solubility values of IP in the studied systems can be found in the
Supplementary Materials.

A quite similar picture arises when inspecting the results obtained for ketoprofen.
Again, the eutectics formed with choline chloride achieved a higher solubility of KP than
the ones with betaine. As in the case of IP, the 1:2 molar ratio was found to be the most
efficient, followed by the 1:4 and 1:1 molar ratios, both in the case of choline chloride and
betaine. There was, however, a difference in the effectiveness of the HBD compounds.
In the case of choline chloride, TEG was the best-performing hydrogen bond donor in
all molar ratios. For betaine, this position was taken by GLY, again in all molar ratios.
The decreasing solubility order of KP in eutectics with choline chloride is the following:
ketoprofen. Then, the most effective solvent systems were systematically studied at dif-
ferent temperatures and by dilution with water. This aspect of adding water to the DES
often leads to a surprising structure [92,93] and solubilizing properties [79,83]. Water is
considered an antisolvent because the aqueous solubility of various compounds is sev-
eral orders of magnitude lower compared to DESs, and dilution with water in excess
amounts significantly reduces the solubility of many solutes in such DES-water mixtures.
However, there are narrow ranges of water content that promote solubility compared to
pure DES. Nonetheless, due to the complicated and complex structure of aqueous DES
mixtures, water cannot be considered a cosolvent in such ranges. Using the example of
deep eutectic choline chloride/urea/water solvent mixtures [92], it was shown that the
DES nanostructure is formulated and preserved in a wide hydration range up to a remark-
ably high water content. This is attributed to the solvophobic sequestration of water in
nanostructured domains around cholinium cations. At higher water concentrations, such
separation becomes unfavorable, and the DES nanostructure is disrupted. That is why
DESs are often highly hygroscopic, and the control of the water content is so important. If
excess amounts of water are added, the water-water and DES-water interactions become
dominant, and only in such a regime can DES-water mixtures be considered as aqueous
solutions of DES components. The presence of these nanostructures might be regarded
as the structural origin of the observed increased solubility of the dissolved substances
compared to dry DES.

2.1. *Solubility of IP and KP in DESs*

In the initial phase of the experiments, a total of 36 eutectic systems were studied for both ibuprofen (IP) and ketoprofen (KP). One arrives at this number by combining two hydrogen bond acceptors (HBAs), namely, choline chloride (ChCl) and betaine (BI), with six hydrogen bond donors (HBDs), namely ethylene glycol (ETG), diethylene glycol (DEG), triethylene glycol (TEG), 1,3-butanediol (B3D), 1,2-propanediol (P2D), and glycerol (GLY), as well as three molar HBA:HBD proportions, namely 1:1, 1:2, and 1:4. This phase was aimed not only at determining the solubility of IP and KP in neat DES but also at selecting the best systems for next-stage studies involving aqueous mixtures. The composition screening was conducted at 25 °C.

In the case of ibuprofen, the following observations could be made. First of all, the
eutectic systems involving choline chloride were slightly more effective than those utilizing
betaine. Secondly, when comparing the effectiveness of HBD components, TEG turned
out to be the most effective in almost all cases, the only exception being the ChCl-GLY
among the 1:2 systems. Finally, when comparing molar ratios, the 1:2 ratio was the most
efficient, with the 1:4 ratio coming second and followed by the 1:1 ratio. For both choline
chloride and betaine in the 1:2 molar ratio, the following decreasing order of solubilizing
potential was noticed: TEG > GLY > DEG > ETG > B3D > P2D. This order varies only
slightly in other proportions, with GLY being the most effective in the case of 1:1 molar
ratio, ETG surpassing DEG in the 1:4 molar proportion, and P2D outperforming B3D on
two occasions. The actual solubility of IP at 25 °C in the ChCl-TEG (1:2) system was found
to be xIP = 0.353, while the ChCl-GLY (1:2) system amounted to a solubility of xIP = 0.349.
In the case of betaine, these values were xIP = 0.321 for BI-TEG (1:2) and xIP = 0.297 for
BI-GLY (1:2). Detailed solubility values of IP in the studied systems can be found in the
Supplementary Materials.

A quite similar picture arises when inspecting the results obtained for ketoprofen.
Again, the eutectics formed with choline chloride achieved a higher solubility of KP than
the ones with betaine. As in the case of IP, the 1:2 molar ratio was found to be the most
efficient, followed by the 1:4 and 1:1 molar ratios, both in the case of choline chloride and
betaine. There was, however, a difference in the effectiveness of the HBD compounds.
In the case of choline chloride, TEG was the best-performing hydrogen bond donor in
all molar ratios. For betaine, this position was taken by GLY, again in all molar ratios.
The decreasing solubility order of KP in eutectics with choline chloride is the following:
TEG > DEG > GLY > ETG > B3D > P2D, with glycerol replacing DEG in the 1:4 molar ratio. On the other hand, for the systems involving betaine, the order is GLY > TEG > DEG > ETG > B3D > P2D, with ETG replacing DEG in the 1:4 molar ratio. The solubilities of KP among the most effective systems for both choline chloride and betaine were x<sub>KP</sub> = 0.067 for ChCl-TEG (1:2), x<sub>KP</sub> = 0.058 for ChCl-DEG (1:2), x<sub>KP</sub> = 0.062 for BI-TEG (1:2), and x<sub>KP</sub> = 0.060 for BI-TEG (1:2). Detailed KP solubility values for the considered systems are presented in the Supplementary Materials.

Based on the performed initial study, the two best-performing eutectic systems were selected for further solubility measurements, one involving choline chloride (ChCl-TEG 1:2 for both IP and KP) and one with betaine (ChCl-TEG 1:2 for IP and ChCl-GLY 1:2 for KP). These systems were used in conjunction with water, forming ternary aqueous mixtures with varying amounts of DESs. The studies involved several temperatures ranging from 25 °C to 40 °C, and it was no surprise that elevated temperatures promoted the dissolution of the studied compounds. Again, it should be noted that the complex nature of aqueous DES systems means that they cannot be described in terms of simple cosolvency or antisolvency effects. Nonetheless, these apparent effects can still be regarded as a useful way of presenting and describing solubility data.

For both IP and KP, regardless of the applied DES, the addition of the eutectic to water increases the solubility of the active pharmaceutical ingredient up to a certain point at which the solubility is greater than in the neat DES. After achieving this infliction point, the solubility decreases to the level of the neat DES. Hence, it is useful to describe this behavior as an apparent cosolvency/antisolvency effect. The details of the solubility curves, however, vary among the studied systems and are presented in Figure 1. For both IP and KP, the systems involving choline chloride are characterized by the highest solubility obtained in the x\*<sub>DES</sub> = 0.9 composition. At 25 °C, the solubility of IP in the ChCl-TEG (1:2, x\*<sub>DES</sub> = 0.9) system was found to be x<sub>IP</sub> = 0.375, while the solubility of KP in the analogous system was x<sub>KP</sub> = 0.068. On the other hand, when inspecting the solubility of ibuprofen in betaine-containing aqueous DES mixtures, it was found that the optimal composition is x\*<sub>DES</sub> = 0.8. At this point, the IP solubility in the BI-TEG (1:2, x\*<sub>DES</sub> = 0.8) system is x<sub>IP</sub> = 0.380, which, interestingly, is greater than for the optimal system involving choline chloride. In the case of ketoprofen solubility, using aqueous mixtures of DESs with betaine results in the optimal composition of x\*<sub>DES</sub> = 0.9, although the difference in solubility with the x\*<sub>DES</sub> = 0.8 composition is less than 1%. For the BI-GLY (1:2, x\*<sub>DES</sub> = 0.9) system, the KP solubility equals x<sub>KP</sub> = 0.072 at 25 °C. Detailed solubility values can be obtained from the Supplementary Materials.

When comparing the obtained results with the literature data, it is evident that in the case of ibuprofen, the optimal aqueous DES mixtures outperform even the most efficient organic solvent, namely, chloroform, characterized by an IP solubility equal to x<sub>IP</sub> = 0.371 at 25 °C [94,95]. The difference is even greater when comparing the solubilities achieved in other most efficient organic solvents, such as octanol, dichloromethane, or acetone. Furthermore, chloroform can hardly be considered a green solvent, which makes the usage of the studied DES systems even more attractive. On the other hand, such high solubilization efficiency was not achieved in the case of ketoprofen. In this case, the optimal DES-water composition was outperformed even by methanol, which resulted in the solubility of KP equal to x<sub>IP</sub> = 0.158 at 25 °C [96]. The difference in the DESs' performance is visualized in Figure 2, which shows the logarithmic mole fraction solubility of the two studied compounds in various solvents. In general, ketoprofen is characterized by a lower solubility than ibuprofen in all systems, although in the case of water, the solubilities are comparably low. Applying DESs as solubilization media drastically increases the solubility compared to water in both cases. The differences between the two compounds are evident, however, when looking at the efficiency of DESs compared to methanol. While for ibuprofen they offer a substantial solubility advantage, their performance in the case of ketoprofen is not optimal.
TEG > DEG > GLY > ETG > B3D > P2D, with glycerol replacing DEG in the 1:4 molar ratio. On the other hand, for the systems involving betaine, the order is GLY > TEG > DEG > ETG > B3D > P2D, with ETG replacing DEG in the 1:4 molar ratio. The solubilities of KP among the most effective systems for both choline chloride and betaine were x<sub>KP</sub> = 0.067 for ChCl-TEG (1:2), x<sub>KP</sub> = 0.058 for ChCl-DEG (1:2), x<sub>KP</sub> = 0.062 for BI-TEG (1:2), and x<sub>KP</sub> = 0.060 for BI-TEG (1:2). Detailed KP solubility values for the considered systems are presented in the Supplementary Materials.

Based on the performed initial study, the two best-performing eutectic systems were selected for further solubility measurements, one involving choline chloride (ChCl-TEG 1:2 for both IP and KP) and one with betaine (ChCl-TEG 1:2 for IP and ChCl-GLY 1:2 for KP). These systems were used in conjunction with water, forming ternary aqueous mixtures with varying amounts of DESs. The studies involved several temperatures ranging from 25 °C to 40 °C, and it was no surprise that elevated temperatures promoted the dissolution of the studied compounds. Again, it should be noted that the complex nature of aqueous DES systems means that they cannot be described in terms of simple cosolvency or antisolvency effects. Nonetheless, these apparent effects can still be regarded as a useful way of presenting and describing solubility data.

For both IP and KP, regardless of the applied DES, the addition of the eutectic to water increases the solubility of the active pharmaceutical ingredient up to a certain point at which the solubility is greater than in the neat DES. After achieving this infliction point, the solubility decreases to the level of the neat DES. Hence, it is useful to describe this behavior as an apparent cosolvency/antisolvency effect. The details of the solubility curves, however, vary among the studied systems and are presented in Figure 1. For both IP and KP, the systems involving choline chloride are characterized by the highest solubility obtained in the x\*<sub>DES</sub> = 0.9 composition. At 25 °C, the solubility of IP in the ChCl-TEG (1:2, x\*<sub>DES</sub> = 0.9) system was found to be x<sub>IP</sub> = 0.375, while the solubility of KP in the analogous system was x<sub>KP</sub> = 0.068. On the other hand, when inspecting the solubility of ibuprofen in betaine-containing aqueous DES mixtures, it was found that the optimal composition is x\*<sub>DES</sub> = 0.8. At this point, the IP solubility in the BI-TEG (1:2, x\*<sub>DES</sub> = 0.8) system is x<sub>IP</sub> = 0.380, which, interestingly, is greater than for the optimal system involving choline chloride. In the case of ketoprofen solubility, using aqueous mixtures of DESs with betaine results in the optimal composition of x\*<sub>DES</sub> = 0.9, although the difference in solubility with the x\*<sub>DES</sub> = 0.8 composition is less than 1%. For the BI-GLY (1:2, x\*<sub>DES</sub> = 0.9) system, the KP solubility equals x<sub>KP</sub> = 0.072 at 25 °C. Detailed solubility values can be obtained from the Supplementary Materials.

When comparing the obtained results with the literature data, it is evident that in the case of ibuprofen, the optimal aqueous DES mixtures outperform even the most efficient organic solvent, namely, chloroform, characterized by an IP solubility equal to x<sub>IP</sub> = 0.371 at 25 °C [94,95]. The difference is even greater when comparing the solubilities achieved in other most efficient organic solvents, such as octanol, dichloromethane, or acetone. Furthermore, chloroform can hardly be considered a green solvent, which makes the usage of the studied DES systems even more attractive. On the other hand, such high solubilization efficiency was not achieved in the case of ketoprofen. In this case, the optimal DES-water composition was outperformed even by methanol, which resulted in the solubility of KP equal to x<sub>IP</sub> = 0.158 at 25 °C [96]. The difference in the DESs' performance is visualized in Figure 2, which shows the logarithmic mole fraction solubility of the two studied compounds in various solvents. In general, ketoprofen is characterized by a lower solubility than ibuprofen in all systems, although in the case of water, the solubilities are comparably low. Applying DESs as solubilization media drastically increases the solubility compared to water in both cases. The differences between the two compounds are evident, however, when looking at the efficiency of DESs compared to methanol. While for ibuprofen they offer a substantial solubility advantage, their performance in the case of ketoprofen is not optimal.
<figure><img src="image_4.png" /><figcaption>**Figure 1.** The solubility curves of ibuprofen (IP, **left panel**) and ketoprofen (KP, **right panel**) in aqueous DES mixtures at various temperatures, expressed as solvent-composition-related mole fractions. ChCl stands for eutectics with choline chloride, and BI stands for eutectics with betaine. $x^*_{DES}$ stands for mole fractions of solute-free DES in aqueous mixtures. For comparison, the room-temperature solubilities of IP in chloroform [94,95] and KP in methanol [96] are provided as reported in the literature.</figcaption></figure>

<figure><img src="image_5.png" /><figcaption>**Figure 2.** Comparison of IP and KP solubility in the most efficient DES and selected non-DES solvents at room temperature. DES1 and DES2 denote aqueous DES mixtures. For ibuprofen, DES1—ChCl + TEG 1:2 ($x^*_{DES} = 0.9$); DES2—BI + TEG 1:2 ($x^*_{DES} = 0.8$). For ketoprofen, DES1—ChCl + TEG 1:2 ($x^*_{DES} = 0.9$); DES2—BI + TEG 1:2 ($x^*_{DES} = 0.9$).</figcaption></figure>

## 2.2. Solubility Dataset

The experimentally determined solubility of both studied COX inhibitors proved the effectiveness of the DES design for optimizing the solubilizers' performance. It is interesting to provide a broader overview of the distribution of this physical property based on already published data. Hence, a careful literature search was undertaken to extend the solubility dataset. This was conducted by the inclusion of reported solubility in neat solvents and binary solvent mixtures. In the case of IP, such data are available for 18 neat solvents,
<figure><img src="image_6.png" /><figcaption>**Figure 1.** The solubility curves of ibuprofen (IP, **left panel**) and ketoprofen (KP, **right panel**) in aqueous DES mixtures at various temperatures, expressed as solvent-composition-related mole fractions. ChCl stands for eutectics with choline chloride, and BI stands for eutectics with betaine. $x^*_{DES}$ stands for mole fractions of solute-free DES in aqueous mixtures. For comparison, the room-temperature solubilities of IP in chloroform [94,95] and KP in methanol [96] are provided as reported in the literature.</figcaption></figure>

<figure><img src="image_7.png" /><figcaption>**Figure 2.** Comparison of IP and KP solubility in the most efficient DES and selected non-DES solvents at room temperature. DES1 and DES2 denote aqueous DES mixtures. For ibuprofen, DES1—ChCl + TEG 1:2 ($x^*_{DES} = 0.9$); DES2—BI + TEG 1:2 ($x^*_{DES} = 0.8$). For ketoprofen, DES1—ChCl + TEG 1:2 ($x^*_{DES} = 0.9$); DES2—BI + TEG 1:2 ($x^*_{DES} = 0.9$).</figcaption></figure>

## 2.2. Solubility Dataset

The experimentally determined solubility of both studied COX inhibitors proved the effectiveness of the DES design for optimizing the solubilizers' performance. It is interesting to provide a broader overview of the distribution of this physical property based on already published data. Hence, a careful literature search was undertaken to extend the solubility dataset. This was conducted by the inclusion of reported solubility in neat solvents and binary solvent mixtures. In the case of IP, such data are available for 18 neat solvents,
including methanol [94], ethanol [94,97–100], 1-propanol [97], 2-Propanol [94,97,101,102], 1-butanol [97], isobutanol [97], 1-pentanol [97], 3-methyl-1-butanol [97], 1-octanol [95,100], propylene glycol [98], ethyl acetate [94,97], isopropyl acetate [51], isopropyl myristate [51], acetone [94,103], 2-butanone [51], 4-methyl-2-pentanone [94], 3-methyl-2-butanone [51], acetonitrile [51], cyclohexane [95], toluene [94], heptane [51], dichloromethane [103], chloroform [94,95], and water [99,104]. Additionally, there are four reports documenting the solubility of IP in binary mixtures such as water–methanol [105], water–ethanol [99], water–propylene glycol [104], and ethanol–propylene glycol [98]. The solubility of S-enantiomer, (S)-(+)-ibuprofen, was measured [106] in methanol, 1-butylanol, isobutanol, 1-pentanol, and 1-octanol. The saturated solutions of ketoprofen were much less intensively studied, and there is only one comprehensive study [96] providing the solubility in the following neat solvents: methanol, isopropanol, n-butanol, acetonitrile, ethyl acetate, 1,4-dioxane, and toluene. The work by Soto et al. [96] is particularly interesting since the authors compared three different methods of experimental solubility determination. The concluded consensus between measured solubility data makes this report highly reliable. There are also other reports by Gantiva et al. [107,108] providing the solubility of KP in aqueous binary mixtures of ethanol and propylene glycol. However, the reliability of the provided values is debatable. Indeed, the reserve of these data was already emphasized by Filippa et al. [109]. Serious discrepancies were noticed and highlighted, suggesting that the aqueous solubility of ketoprofen at 298.15 K [107,108] (4.295 × 10⁻⁴ mol/L) is much lower than the value published by Yalkowsky et al. [110], which is 5.58 × 10⁻³ mol/L. The same reserve applies to the solubility of ketoprofen in cyclohexane, since, as reported in [111], the solubility of ketoprofen in cyclohexane is 5.50 × 10⁻⁴ mol/L at 298 K, while according to the study of Filippa et al. [109], the solubility of this compound is 4.16 × 10⁻³ mol/L at 300.15 K. Also, the solubility value determined for propylene glycol seems problematic, since according to the measurements of Filippa et al. [109], the solubility of ketoprofen in this solvent is 0.499 mol/L (300.15 K), while the value reported by Gantiva et al. [107,108], for similar temperature conditions (298.15 K), was 0.1832 mol/L. In order to resolve the discrepancies and provide a consistent set of solubility in the DES-water systems studied here, new measurements of both IP and KP in water were performed. The obtained data are shown in Figure 3, where the results of other studies are also provided.

As evidenced by the plots in Figure 3, there is a pretty solid consensus on IP solubility in water. The deviations between our measurements and those by Manrique et al. [104] are marginal and within standard deviations. To the contrary, serious discrepancies are observed for the solubility of KP in water. It happened that the results of our measurements agree with data provided by Jiménez et al. [105] for ambient conditions and are seriously higher compared to data provided by Gantiva et al. [107,108]. This raises the question whether the latter values correspond to the equilibrium state of the KP-water system. Although the congruency of solubility at 25 °C between our measurements and the available literature might serve as an argument, one can obtain additional clues based on the electron charge distributions of IP and KP molecules. As provided in Figure 4, the electron density distributions for both studied compounds are characterized mostly by non-polar regions (green), with the exception of the carboxylic group, which induces local positively (blue) and negatively (red) charged surfaces due to the presence of hydroxyl and keto groups, respectively. KP also has an additional negatively charged region due to the presence of the oxygen atom in the moiety. Hence, the observed low solubility of both COX inhibitors in water is due to the mostly apolar nature of these compounds. However, it seems to be rational to expect a higher solubility of KP in water compared to IP, as the oxygen atom can be actively involved in hydration. In addition, it is worth adding that both compounds are able to form stable dimers, which is rather expected bearing in mind that they belong to the group of carboxylic acids. Hence, as a result, self-association is highly probable, and stabilization is to be granted from the cyclic structure of the C₂ type, meaning that the two hydrogen bonds are involved in an eight-center ring creation. The occurrence of such strong associations is typical for all carboxylic acids. The involvement of carboxylic groups of two
including methanol [94], ethanol [94,97–100], 1-propanol [97], 2-Propanol [94,97,101,102], 1-butanol [97], isobutanol [97], 1-pentanol [97], 3-methyl-1-butanol [97], 1-octanol [95,100], propylene glycol [98], ethyl acetate [94,97], isopropyl acetate [51], isopropyl myristate [51], acetone [94,103], 2-butanone [51], 4-methyl-2-pentanone [94], 3-methyl-2-butanone [51], acetonitrile [51], cyclohexane [95], toluene [94], heptane [51], dichloromethane [103], chloroform [94,95], and water [99,104]. Additionally, there are four reports documenting the solubility of IP in binary mixtures such as water–methanol [105], water–ethanol [99], water–propylene glycol [104], and ethanol–propylene glycol [98]. The solubility of S-enantiomer, (S)-(+)-ibuprofen, was measured [106] in methanol, 1-butylanol, isobutanol, 1-pentanol, and 1-octanol. The saturated solutions of ketoprofen were much less intensively studied, and there is only one comprehensive study [96] providing the solubility in the following neat solvents: methanol, isopropanol, n-butanol, acetonitrile, ethyl acetate, 1,4-dioxane, and toluene. The work by Soto et al. [96] is particularly interesting since the authors compared three different methods of experimental solubility determination. The concluded consensus between measured solubility data makes this report highly reliable. There are also other reports by Gantiva et al. [107,108] providing the solubility of KP in aqueous binary mixtures of ethanol and propylene glycol. However, the reliability of the provided values is debatable. Indeed, the reserve of these data was already emphasized by Filippa et al. [109]. Serious discrepancies were noticed and highlighted, suggesting that the aqueous solubility of ketoprofen at 298.15 K [107,108] (4.295 × 10⁻⁴ mol/L) is much lower than the value published by Yalkowsky et al. [110], which is 5.58 × 10⁻³ mol/L. The same reserve applies to the solubility of ketoprofen in cyclohexane, since, as reported in [111], the solubility of ketoprofen in cyclohexane is 5.50 × 10⁻⁴ mol/L at 298 K, while according to the study of Filippa et al. [109], the solubility of this compound is 4.16 × 10⁻³ mol/L at 300.15 K. Also, the solubility value determined for propylene glycol seems problematic, since according to the measurements of Filippa et al. [109], the solubility of ketoprofen in this solvent is 0.499 mol/L (300.15 K), while the value reported by Gantiva et al. [107,108], for similar temperature conditions (298.15 K), was 0.1832 mol/L. In order to resolve the discrepancies and provide a consistent set of solubility in the DES–water systems studied here, new measurements of both IP and KP in water were performed. The obtained data are shown in Figure 3, where the results of other studies are also provided.

As evidenced by the plots in Figure 3, there is a pretty solid consensus on IP solubility in water. The deviations between our measurements and those by Manrique et al. [104] are marginal and within standard deviations. To the contrary, serious discrepancies are observed for the solubility of KP in water. It happened that the results of our measurements agree with data provided by Jiménez et al. [105] for ambient conditions and are seriously higher compared to data provided by Gantiva et al. [107,108]. This raises the question whether the latter values correspond to the equilibrium state of the KP–water system. Although the congruency of solubility at 25 °C between our measurements and the available literature might serve as an argument, one can obtain additional clues based on the electron charge distributions of IP and KP molecules. As provided in Figure 4, the electron density distributions for both studied compounds are characterized mostly by non-polar regions (green), with the exception of the carboxylic group, which induces local positively (blue) and negatively (red) charged surfaces due to the presence of hydroxyl and keto groups, respectively. KP also has an additional negatively charged region due to the presence of the oxygen atom in the moiety. Hence, the observed low solubility of both COX inhibitors in water is due to the mostly apolar nature of these compounds. However, it seems to be rational to expect a higher solubility of KP in water compared to IP, as the oxygen atom can be actively involved in hydration. In addition, it is worth adding that both compounds are able to form stable dimers, which is rather expected bearing in mind that they belong to the group of carboxylic acids. Hence, as a result, self-association is highly probable, and stabilization is to be granted from the cyclic structure of the C₂ type, meaning that the two hydrogen bonds are involved in an eight-center ring creation. The occurrence of such strong associations is typical for all carboxylic acids. The involvement of carboxylic groups of two
IP molecules is particularly important, since these are the only polar centers that might be involved in direct contact with polar solvent molecules. Indeed, as it is documented in Figure 4, a strong hydrogen bond is formed between the donating hydroxyl group of IP and the accepting oxygen atom in water. Since the rest of the molecule is non-polar, it is expected that the polarity of the IP dimer will be much smaller compared to the monomer. A similar pattern can also be noticed for ketoprofen. In either case, the dominant surface exposed for contacts with solvent molecules is mostly non-polar, and the involvement of a carboxylic group in the dimer formation further reduces the net polar region available for contacts with polar molecules.

<figure><img src="image_6.png" /><figcaption>**Figure 3.** The temperature dependence of IP and KP solubility in water. The newly measured data presented in this work were augmented with already published results, namely, [a] by Manrique et al. [104], [b] by Gantiva et al. [107,108], and [c] by Jiménez et al. [105].</figcaption></figure>

Also, in the case of dimers, ketoprofen seems to be more prone to interactions with polar solvents, since the oxygen atom attached to the carbon atom chain is not affected by dimerization. Hence, both qualitative analysis of electron charge distributions and congruency between our measurements and the ones provided by Jiménez et al. [105] suggest that the reports by Gantiva et al. [107,108] should be taken with reserve.

This conclusion drawn from the above discussion is important from the perspective of machine learning, since the quality of the trained model is crucially dependent on the consistency and representativeness of the data used for training purposes. That is why the questioned data were not included in the training dataset. However, to make the model as reliable and as general as possible, the dataset was extended by the inclusion of compounds belonging to the same group. This enabled the capture of a wider range of correlations between measurements and molecular features. There is, however, a limitation to the choice, since two criteria must be met, namely, not only reliable solubility data should be included, but fusion data must also be available. The latter requirement is related to the fact that the molecular descriptor of the first choice is the value of solubility computed using COSMOtherm (version 24.0.0). This seriously restricts the potential scope of compounds to be included, and, to our best knowledge, the solubility and fusion data were only recorded for two closely related compounds. The first one, flurbiprofen, was studied in 1-octanol [100,112], ethanol [100], and hexane [112]. For the second one, dexibuprofen, there are available solubility values reported by Wang et al. [106] for methanol, n-butyl
IP molecules is particularly important, since these are the only polar centers that might be involved in direct contact with polar solvent molecules. Indeed, as it is documented in Figure 4, a strong hydrogen bond is formed between the donating hydroxyl group of IP and the accepting oxygen atom in water. Since the rest of the molecule is non-polar, it is expected that the polarity of the IP dimer will be much smaller compared to the monomer. A similar pattern can also be noticed for ketoprofen. In either case, the dominant surface exposed for contacts with solvent molecules is mostly non-polar, and the involvement of a carboxylic group in the dimer formation further reduces the net polar region available for contacts with polar molecules.

<figure><img src="image_8.png" /><figcaption>**Figure 3.** The temperature dependence of IP and KP solubility in water. The newly measured data presented in this work were augmented with already published results, namely, [a] by Manrique et al. [104], [b] by Gantiva et al. [107,108], and [c] by Jiménez et al. [105].</figcaption></figure>

Also, in the case of dimers, ketoprofen seems to be more prone to interactions with polar solvents, since the oxygen atom attached to the carbon atom chain is not affected by dimerization. Hence, both qualitative analysis of electron charge distributions and congruency between our measurements and the ones provided by Jiménez et al. [105] suggest that the reports by Gantiva et al. [107,108] should be taken with reserve.

This conclusion drawn from the above discussion is important from the perspective of machine learning, since the quality of the trained model is crucially dependent on the consistency and representativeness of the data used for training purposes. That is why the questioned data were not included in the training dataset. However, to make the model as reliable and as general as possible, the dataset was extended by the inclusion of compounds belonging to the same group. This enabled the capture of a wider range of correlations between measurements and molecular features. There is, however, a limitation to the choice, since two criteria must be met, namely, not only reliable solubility data should be included, but fusion data must also be available. The latter requirement is related to the fact that the molecular descriptor of the first choice is the value of solubility computed using COSMOtherm (version 24.0.0). This seriously restricts the potential scope of compounds to be included, and, to our best knowledge, the solubility and fusion data were only recorded for two closely related compounds. The first one, flurbiprofen, was studied in 1-octanol [100,112], ethanol [100], and hexane [112]. For the second one, dexibuprofen, there are available solubility values reported by Wang et al. [106] for methanol, n-butyl
alcohol, isobutyl alcohol, n-amyl alcohol, and n-octyl alcohol. Hence, these data were included in the solubility training dataset.

<figure><img src="image_7.png" /><figcaption>**Figure 4.** The distributions of electron charge density of IP (**left**) and KP (**right**) molecules, as well as their dimers resulting from self-association.</figcaption></figure>

## 2.3. Machine Learning Model

The solubility model was constructed based on the best-performing regressors identified after comprehensive hyperparameter tuning and validation. During model selection, a multi-layer perceptron (MLP) regressor achieved significantly higher accuracy than other nonlinear algorithms for the COX inhibitors dataset. Consequently, the ensemble formulation step simply selected the best individual model, which was the MLP regressor.

The MLP model architecture was optimized in a flexible manner by allowing for the adjustment of both parameter values and the overall network structure. The optimal network consisted of an input layer, multiple hidden layers of variable sizes, and an output layer. Hyperparameter optimization was performed extensively using the Optuna Study framework. The best-performing model had the following hyperparameters: hidden layer sizes of 64, 17, 45, 54, 18, 14, 45, and 61, rectified linear unit (ReLU) activation, L2 regularization (alpha) of 0.045, an adaptive learning rate, a limited-memory Broyden-Fletcher-Goldfarb-Shanno (L-BFGS) solver, and 928,539 maximum iterations.

The multiple hidden layers enable complex nonlinear mapping between inputs and outputs. Larger initial layers (64–45 neurons) capture richer feature representations from descriptors, while smaller later layers (14–45 neurons) perform higher-level abstraction. The ReLU activation introduces nonlinearity while avoiding vanishing gradients. Adaptive learning optimizes the learning rate during training. An L2 regularization of 0.024 prevents overfitting by penalizing model complexity. The optimized MLP architecture and hyperparameters performed best for learning the solubility structure–property relationships from molecular fingerprints in a nonlinear fashion. As the most efficient optimizer, L-BFGS was selected as it can reach a local minimum gradient-based quasi-Newton style. The high maximum iteration count ensured convergence.

The accuracy of the formulated model is graphically illustrated in Figure 5, which documents the correlation between experimental and computed solubility expressed as decadal logarithms. The three dataset splits, namely, training (N = 406), test (N = 88), and
alcohol, isobutyl alcohol, n-amyl alcohol, and n-octyl alcohol. Hence, these data were included in the solubility training dataset.

<figure><img src="image_9.png" /><figcaption>**Figure 4.** The distributions of electron charge density of IP (**left**) and KP (**right**) molecules, as well as their dimers resulting from self-association.</figcaption></figure>

## 2.3. Machine Learning Model

The solubility model was constructed based on the best-performing regressors identified after comprehensive hyperparameter tuning and validation. During model selection, a multi-layer perceptron (MLP) regressor achieved significantly higher accuracy than other nonlinear algorithms for the COX inhibitors dataset. Consequently, the ensemble formulation step simply selected the best individual model, which was the MLP regressor.

The MLP model architecture was optimized in a flexible manner by allowing for the adjustment of both parameter values and the overall network structure. The optimal network consisted of an input layer, multiple hidden layers of variable sizes, and an output layer. Hyperparameter optimization was performed extensively using the Optuna Study framework. The best-performing model had the following hyperparameters: hidden layer sizes of 64, 17, 45, 54, 18, 14, 45, and 61, rectified linear unit (ReLU) activation, L2 regularization (alpha) of 0.045, an adaptive learning rate, a limited-memory Broyden-Fletcher-Goldfarb-Shanno (L-BFGS) solver, and 928,539 maximum iterations.

The multiple hidden layers enable complex nonlinear mapping between inputs and outputs. Larger initial layers (64–45 neurons) capture richer feature representations from descriptors, while smaller later layers (14–45 neurons) perform higher-level abstraction. The ReLU activation introduces nonlinearity while avoiding vanishing gradients. Adaptive learning optimizes the learning rate during training. An L2 regularization of 0.024 prevents overfitting by penalizing model complexity. The optimized MLP architecture and hyperparameters performed best for learning the solubility structure–property relationships from molecular fingerprints in a nonlinear fashion. As the most efficient optimizer, L-BFGS was selected as it can reach a local minimum gradient-based quasi-Newton style. The high maximum iteration count ensured convergence.

The accuracy of the formulated model is graphically illustrated in Figure 5, which documents the correlation between experimental and computed solubility expressed as decadal logarithms. The three dataset splits, namely, training (N = 406), test (N = 88), and
validation (N = 87), are plotted in black ink to distinguish the solubility value distributions determined by COSMO-RS for non-DES and DES systems. It is clearly visible that the back-computed values perfectly match the experimental data. This is not the case for solubility derived directly from COSMO-RS. Neither of the two divergent groups, non-DES and DES systems, can be characterized accurately enough by this approach alone. In the case of solubility, in many DESs, the analyzed solutes are predicted with deviations of several orders of magnitude in mole fraction. However, despite such inaccuracies, the utilization of COSMO-RS-computed solubility remains a valuable source of information for the nonlinear models. This directly justifies the adopted methodology herein. It is worth noting that excluding the computed solubility values from the pool of molecular descriptors led to models with much lower accuracy than the one presented in Figure 5.

<figure><img src="image_8.png" /><figcaption>**Figure 5.** The correlation between experimental and computed solubility values of COX inhibitors in DES and non-DES systems using the developed MLP model and COSMO-RS approach. The accuracy of the MLP model is as follows: RMSD(train) = 0.016, RMSD(test) = 0.118, and RMSD/validation) = 0.136.</figcaption></figure>

Based on the obtained model, a new prediction was made for selected analogues of IP and KP. As mentioned in the methodology section, this step is limited to those compounds for which solubility calculations are possible, i.e., for which fusion data are available. Four closely related compounds were therefore selected, namely, flurbiprofen, felbinac, phenylacetic acid, and diphenylacetic acid. For these compounds, the values of the molecular descriptors were calculated in the same way as for the training dataset, and the prediction was made using optimized values of the MLP architecture. As a result, the hypothetical solubilities were obtained, as shown in Figure 6. Based on these predictions, it is expected to be beneficial to perform actual solubility measurements for all screened analogues, as the solubility is likely to be improved in the choline chloride-TEG deep eutectic mixture. The highest solubility is expected for flurbiprofen in this DES.
validation (N = 87), are plotted in black ink to distinguish the solubility value distributions determined by COSMO-RS for non-DES and DES systems. It is clearly visible that the back-computed values perfectly match the experimental data. This is not the case for solubility derived directly from COSMO-RS. Neither of the two divergent groups, non-DES and DES systems, can be characterized accurately enough by this approach alone. In the case of solubility, in many DESs, the analyzed solutes are predicted with deviations of several orders of magnitude in mole fraction. However, despite such inaccuracies, the utilization of COSMO-RS-computed solubility remains a valuable source of information for the nonlinear models. This directly justifies the adopted methodology herein. It is worth noting that excluding the computed solubility values from the pool of molecular descriptors led to models with much lower accuracy than the one presented in Figure 5.

<figure><img src="image_10.png" /><figcaption>**Figure 5.** The correlation between experimental and computed solubility values of COX inhibitors in DES and non-DES systems using the developed MLP model and COSMO-RS approach. The accuracy of the MLP model is as follows: RMSD(train) = 0.016, RMSD(test) = 0.118, and RMSD/validation) = 0.136.</figcaption></figure>

Based on the obtained model, a new prediction was made for selected analogues of IP and KP. As mentioned in the methodology section, this step is limited to those compounds for which solubility calculations are possible, i.e., for which fusion data are available. Four closely related compounds were therefore selected, namely, flurbiprofen, felbinac, phenylacetic acid, and diphenylacetic acid. For these compounds, the values of the molecular descriptors were calculated in the same way as for the training dataset, and the prediction was made using optimized values of the MLP architecture. As a result, the hypothetical solubilities were obtained, as shown in Figure 6. Based on these predictions, it is expected to be beneficial to perform actual solubility measurements for all screened analogues, as the solubility is likely to be improved in the choline chloride-TEG deep eutectic mixture. The highest solubility is expected for flurbiprofen in this DES.
<figure><img src="image_9.png" /><figcaption>**Figure 6.** Results of predictions made by the solubility model for new closely related compounds at room temperature. DES stands for the deep eutectic involving choline chloride and TEG in a 1:2 molar ratio. Different compositions of the aqueous mixture are presented.</figcaption></figure>

# 3. Materials and Methods

## 3.1. Materials

Two chemical compounds, namely, racemic ibuprofen (IBU, CAS: 15687-27-1, MW = 206.28 g/mol) and racemic ketoprofen (KET, CAS: 22071-15-4, MW = 254.28 g/mol), were used as the main objects of this study. Both were supplied by Sigma Aldrich (Saint Louis, MO, USA) and had a purity of ≥98%. As the constituents of deep eutectic solvents, the following compounds were used: choline chloride (ChCl, CAS: 67-48-1), betaine (BI, CAS: 107-43-7), ethylene glycol (ETG, CAS: 107-21-1), diethylene glycol (DEG, CAS: 111-46-6), triethylene glycol (TEG, CAS: 112-27-6), glycerol (GLY, CAS: 56-81-5), 1,2-propanediol (P2D, CAS: 57-55-6), and 1,3-butanediol (B3D, CAS: 107-88-0). All the above chemicals were similarly obtained from Sigma Aldrich (St. Louis, MO, USA) and had a >99% purity. Also, methanol (CAS: 67-56-1) with a ≥99% purity supplied by Avantor Performance Materials (Gliwice, Poland) was used as a secondary solvent. Apart from choline chloride, which was dried before use, all other compounds were used as supplied without any additional procedures.

## 3.2. Preparation of the Calibration Curves

Stock solutions of ibuprofen and ketoprofen were prepared in order to prepare the two calibration curves in 100 mL volumetric flasks using methanol as a solvent. Subsequent dilutions of these initial solutions were made by transferring fixed amounts of the stock solution into 10 mL volumetric flasks and adding appropriate amounts of methanol. In this way, two series of solutions with decreasing concentrations were obtained. A total of eleven solutions were used, and their concentrations ranged from 0.0126 mg/mL to 0.0504 mg/mL in the case of ibuprofen and from 0.01008 mg/mL to 0.02016 mg/mL in the case of ketoprofen. These solutions were measured using an A360 spectrophotometer from AOE Instruments (Shanghai, China) in the wavelength range from 200 nm to 500 nm. The characteristic wavelengths corresponding to the absorbance maxima were found to be 223 nm and 254 nm for ibuprofen and ketoprofen, respectively. For each of the compounds, separate curves were prepared and averaged. The linear regression equations were obtained as $A = 27.268 \times C + 0.001$ for ibuprofen and as $A = 68.792 \times C + 0.002$ for ketoprofen. In order to validate the curves, the R² determination coefficients, limits of detection (LODs),
<figure><img src="image_11.png" /><figcaption>**Figure 6.** Results of predictions made by the solubility model for new closely related compounds at room temperature. DES stands for the deep eutectic involving choline chloride and TEG in a 1:2 molar ratio. Different compositions of the aqueous mixture are presented.</figcaption></figure>

**3. Materials and Methods**

**3.1. Materials**

Two chemical compounds, namely, racemic ibuprofen (IBU, CAS: 15687-27-1, MW = 206.28 g/mol) and racemic ketoprofen (KET, CAS: 22071-15-4, MW = 254.28 g/mol), were used as the main objects of this study. Both were supplied by Sigma Aldrich (Saint Louis, MO, USA) and had a purity of ≥98%. As the constituents of deep eutectic solvents, the following compounds were used: choline chloride (ChCl, CAS: 67-48-1), betaine (BI, CAS: 107-43-7), ethylene glycol (ETG, CAS: 107-21-1), diethylene glycol (DEG, CAS: 111-46-6), triethylene glycol (TEG, CAS: 112-27-6), glycerol (GLY, CAS: 56-81-5), 1,2-propanediol (P2D, CAS: 57-55-6), and 1,3-butanediol (B3D, CAS: 107-88-0). All the above chemicals were similarly obtained from Sigma Aldrich (St. Louis, MO, USA) and had a >99% purity. Also, methanol (CAS: 67-56-1) with a ≥99% purity supplied by Avantor Performance Materials (Gliwice, Poland) was used as a secondary solvent. Apart from choline chloride, which was dried before use, all other compounds were used as supplied without any additional procedures.

**3.2. Preparation of the Calibration Curves**

Stock solutions of ibuprofen and ketoprofen were prepared in order to prepare the two calibration curves in 100 mL volumetric flasks using methanol as a solvent. Subsequent dilutions of these initial solutions were made by transferring fixed amounts of the stock solution into 10 mL volumetric flasks and adding appropriate amounts of methanol. In this way, two series of solutions with decreasing concentrations were obtained. A total of eleven solutions were used, and their concentrations ranged from 0.0126 mg/mL to 0.0504 mg/mL in the case of ibuprofen and from 0.01008 mg/mL to 0.02016 mg/mL in the case of ketoprofen. These solutions were measured using an A360 spectrophotometer from AOE Instruments (Shanghai, China) in the wavelength range from 200 nm to 500 nm. The characteristic wavelengths corresponding to the absorbance maxima were found to be 223 nm and 254 nm for ibuprofen and ketoprofen, respectively. For each of the compounds, separate curves were prepared and averaged. The linear regression equations were obtained as $A = 27.268 \times C + 0.001$ for ibuprofen and as $A = 68.792 \times C + 0.002$ for ketoprofen. In order to validate the curves, the R² determination coefficients, limits of detection (LODs),
and limits of quantification (LOQs) were calculated. For the ibuprofen curve, the following parameters were obtained: R² = 0.9968, LOD = 0.00233 mg/mL, and LOQ = 0.00698 mg/mL. For the ketoprofen curve, these parameters were R² = 0.9979, LOD = 0.00052 mg/mL, and LOQ = 0.00156 mg/mL. The above values indicate the adequacy of the prepared curves for measuring the solubility of the two considered compounds.

### 3.3. Sample Preparation and Solubility Measurements

For the determination of the solubilities of ibuprofen and ketoprofen in the studied systems, the well-known and widely used shake-flask method was applied. This method has been used numerous times by many researchers, including our research team [113–116], and has proven its reliability through the comparison of obtained solubility values with the literature data. The choice of the eutectic components was based on our previous experiences, which resulted in the elevated solubility of many APIs [79,83,115].

The formulation of deep eutectic solvents involved using either choline chloride or betaine as a hydrogen bond acceptor (HBA) and one of six studied polyols, namely, TEG, DED, ETG, GLY, P2D, or B3D, as a hydrogen bond donor (HBD). The preparation of a particular DES system involved mixing the appropriate HBA with a selected HBD in glass vessels in appropriate molar ratios. Such mixtures were then placed on a heating plate and further mixed until they formed a homogenous solution. The created DESs were used both in their neat form and in aqueous mixtures with specifically calculated amounts of both water and eutectic.

The first step in the determination of the solubilities of ibuprofen and ketoprofen in the studied systems was the preparation of their saturated solutions. In order to achieve this, an excess amount of the compound was placed in a test tube, which was later filled with the pure DES or water-DES mixture. The solubilities were measured at 25 °C, 30 °C, 35 °C, and 40 °C. In order to achieve a stable temperature, the Orbital Shaker Incubator ES-20/60 from Biosan (Riga, Latvia) was used. The samples were incubated for 24 h with simultaneous mixing at 60 rev/min. In the next step, the samples were filtered using a syringe equipped with a 0.22 µm pore size PTFE filter. To prevent precipitation, the test tubes, syringes, pipette tips, and filters were initially heated at the same temperature as the measured sample. In order to remain in the linearity range of the curve, the filtered samples were diluted with methanol before spectrophotometric measurements. The spectra of the samples were recorded on the A360 spectrophotometer, with a 200 nm–500 nm wavelength range and a 1 nm resolution. Methanol was used for the calibration of the spectrophotometer. The absorbance at characteristic wavelengths, i.e., 223 nm for ibuprofen and 254 nm for ketoprofen, was recorded, and the concentration of the compounds in the samples was calculated based on the calibration curve obtained earlier. For the determination of the mole fractions of both compounds, the density of the samples was additionally determined. The procedure relied on weighing 1 mL of each solution in a 10 mL volumetric flask. For this task, a RADWAG (Radom, Poland) AS 110 R2.PLUS analytical balance was used with 0.1 mg precision. Three samples were prepared and measured for each system, and the final values are a result of averaging. In order to evaluate the quality of the obtained results, the standard deviation values were calculated. All experimental results are provided in Supplementary Tables S1 and S2.

### 3.4. Machine Learning Protocol

The training of nonlinear models was performed according to the scheme already applied to our previous projects [79,83,117,118]. Since details were already presented, here, only a brief synopsis is provided. Based on Python code, the comprehensive series of regressors were trained by means of hyperparameter tuning. The solubility prediction model was crafted utilizing in-house-developed Python (version 3.10, https://www.python.org/) code tailored for hyperparameter tuning across 36 regression models. These models encompass diverse algorithms, spanning linear models, boosting, ensembles, nearest neighbors, neural networks, and other regressor types. The hyperparameter space was meticulously explored
and limits of quantification (LOQs) were calculated. For the ibuprofen curve, the following parameters were obtained: R² = 0.9968, LOD = 0.00233 mg/mL, and LOQ = 0.00698 mg/mL. For the ketoprofen curve, these parameters were R² = 0.9979, LOD = 0.00052 mg/mL, and LOQ = 0.00156 mg/mL. The above values indicate the adequacy of the prepared curves for measuring the solubility of the two considered compounds.

### 3.3. Sample Preparation and Solubility Measurements

For the determination of the solubilities of ibuprofen and ketoprofen in the studied systems, the well-known and widely used shake-flask method was applied. This method has been used numerous times by many researchers, including our research team [113–116], and has proven its reliability through the comparison of obtained solubility values with the literature data. The choice of the eutectic components was based on our previous experiences, which resulted in the elevated solubility of many APIs [79,83,115].

The formulation of deep eutectic solvents involved using either choline chloride or betaine as a hydrogen bond acceptor (HBA) and one of six studied polyols, namely, TEG, DED, ETG, GLY, P2D, or B3D, as a hydrogen bond donor (HBD). The preparation of a particular DES system involved mixing the appropriate HBA with a selected HBD in glass vessels in appropriate molar ratios. Such mixtures were then placed on a heating plate and further mixed until they formed a homogenous solution. The created DESs were used both in their neat form and in aqueous mixtures with specifically calculated amounts of both water and eutectic.

The first step in the determination of the solubilities of ibuprofen and ketoprofen in the studied systems was the preparation of their saturated solutions. In order to achieve this, an excess amount of the compound was placed in a test tube, which was later filled with the pure DES or water-DES mixture. The solubilities were measured at 25 °C, 30 °C, 35 °C, and 40 °C. In order to achieve a stable temperature, the Orbital Shaker Incubator ES-20/60 from Biosan (Riga, Latvia) was used. The samples were incubated for 24 h with simultaneous mixing at 60 rev/min. In the next step, the samples were filtered using a syringe equipped with a 0.22 µm pore size PTFE filter. To prevent precipitation, the test tubes, syringes, pipette tips, and filters were initially heated at the same temperature as the measured sample. In order to remain in the linearity range of the curve, the filtered samples were diluted with methanol before spectrophotometric measurements. The spectra of the samples were recorded on the A360 spectrophotometer, with a 200 nm–500 nm wavelength range and a 1 nm resolution. Methanol was used for the calibration of the spectrophotometer. The absorbance at characteristic wavelengths, i.e., 223 nm for ibuprofen and 254 nm for ketoprofen, was recorded, and the concentration of the compounds in the samples was calculated based on the calibration curve obtained earlier. For the determination of the mole fractions of both compounds, the density of the samples was additionally determined. The procedure relied on weighing 1 mL of each solution in a 10 mL volumetric flask. For this task, a RADWAG (Radom, Poland) AS 110 R2.PLUS analytical balance was used with 0.1 mg precision. Three samples were prepared and measured for each system, and the final values are a result of averaging. In order to evaluate the quality of the obtained results, the standard deviation values were calculated. All experimental results are provided in Supplementary Tables S1 and S2.

### 3.4. Machine Learning Protocol

The training of nonlinear models was performed according to the scheme already applied to our previous projects [79,83,117,118]. Since details were already presented, here, only a brief synopsis is provided. Based on Python code, the comprehensive series of regressors were trained by means of hyperparameter tuning. The solubility prediction model was crafted utilizing in-house-developed Python (version 3.10, https://www.python.org/) code tailored for hyperparameter tuning across 36 regression models. These models encompass diverse algorithms, spanning linear models, boosting, ensembles, nearest neighbors, neural networks, and other regressor types. The hyperparameter space was meticulously explored
to identify optimal values using the Optuna study (version 3.2, https://optuna.org/), an open-source Python package for hyperparameter optimization [119]. Model refinement involved 5000 minimization trials employing the tree-structured Parzen estimator (TPE) as the search algorithm sampler. To assess each regression model's performance, a custom score function was defined, amalgamating multiple metrics to gauge accuracy and generalizability, as elucidated in a prior work [118]. The scoring function notably incorporates penalties derived from learning curve analysis (LCA) performed using the scikit-learn 1.2.2 library during parameter tuning. Due to the computational intensity of LCA, only two-point computations were executed, encompassing 50% and 100% of the total data. Subsequently, LCA evaluations of the final model entailed 20-point calculations within the 50–100% data range. The values integrated into the custom loss function correspond to the mean MAE values obtained at the largest training set size. Thus, the custom loss function integrates both accuracy and generalizability aspects, providing insights into the model's performance on unseen data. After determining the optimal parameters for each model, the highest-performing ones were selected for ensemble construction, wherein the weights of individual model contributions to the final solubility prediction were optimized. These weight values were derived through the minimization of the root mean square deviation between the predictions for the validation subset and the observed values. The efficacy of the final ensemble model was subsequently assessed using the test subset to ascertain its accuracy. All models which contribute less than 1% are not included in the ensemble.

### 3.5. Molecular Descriptors

The molecular information encoded in the appropriate series of descriptors is crucial from the perspective of an adequate representation of the relationship between measurements and features used for modeling. There are a variety of approaches intended to probe the molecular hyperspace from different perspectives. This project follows the main tenet that temperature- and concentration-dependent molecular characteristics are to be used for solubility modeling. This initial bias rules out many approaches for molecular descriptor determination, including SMILES-related methods, fingerprints, and also the direct utilization of many first-principle procedures. However, there is a very promising and well-recognized methodology combining both quantum chemistry computations and statistical mechanics to account for environmental influences on the thermodynamic properties of the systems under consideration. This approach, acronymed as COSMO-RS (Conductor-like Screening Model for Realistic Solvents) [120], is designed for bulk systems, in principle, of arbitrary compositions, concentrations, and temperature ranges. That is why this particular framework is selected as a default source of theoretically determined molecular characteristics, as it was already successfully documented in our previous projects [79,83,117,118,121–124]. In particular, the COSMO-RS theory implemented in the COSMOtherm package (version 24.0.0) offers a very convenient way to saturate system characteristics and is used as the source of molecular descriptors also in this project.

The first choice of molecular descriptors is the solubility computed with the aid of COSMOtherm. Unfortunately, this imposes limitations inherent to the approach, and its utilization requires providing the fusion data, such as heat of fusion, temperature of melting, and eventually the heat capacity change upon melting. These experimental parameters are crucial for including the lattice energy. Fortunately, these data are available for both solutes, and several estimates were reported of the melting temperature and heat of fusion. The following mean values were adopted in this work, based on the compilation of data from various sources, as available in the literature [125]: T<sub>m</sub> = 348.64 K and ΔH<sub>m</sub> = 26.69 kJ/mol for ibuprofen and dexibuprofen, T<sub>m</sub> = 367.99 K and ΔH<sub>m</sub> = 25.52 kJ/mol for ketoprofen, and T<sub>m</sub> = 387.34 K and ΔH<sub>m</sub> = 27.86 kJ/mol for flurbiprofen. The COSMOtherm program [126] was used for computing solubility using the newest BP_TZVPD_FINE_24.ctd parametrization. All structures of both solutes and solvent molecules were characterized by conformational diversity derived by comprehensive conformational analysis using the COSMOconf program (version 24.0.0) [127]. Solubilities were computed by fully solving
to identify optimal values using the Optuna study (version 3.2, https://optuna.org/), an open-source Python package for hyperparameter optimization [119]. Model refinement involved 5000 minimization trials employing the tree-structured Parzen estimator (TPE) as the search algorithm sampler. To assess each regression model's performance, a custom score function was defined, amalgamating multiple metrics to gauge accuracy and generalizability, as elucidated in a prior work [118]. The scoring function notably incorporates penalties derived from learning curve analysis (LCA) performed using the scikit-learn 1.2.2 library during parameter tuning. Due to the computational intensity of LCA, only two-point computations were executed, encompassing 50% and 100% of the total data. Subsequently, LCA evaluations of the final model entailed 20-point calculations within the 50–100% data range. The values integrated into the custom loss function correspond to the mean MAE values obtained at the largest training set size. Thus, the custom loss function integrates both accuracy and generalizability aspects, providing insights into the model's performance on unseen data. After determining the optimal parameters for each model, the highest-performing ones were selected for ensemble construction, wherein the weights of individual model contributions to the final solubility prediction were optimized. These weight values were derived through the minimization of the root mean square deviation between the predictions for the validation subset and the observed values. The efficacy of the final ensemble model was subsequently assessed using the test subset to ascertain its accuracy. All models which contribute less than 1% are not included in the ensemble.

### 3.5. Molecular Descriptors

The molecular information encoded in the appropriate series of descriptors is crucial from the perspective of an adequate representation of the relationship between measurements and features used for modeling. There are a variety of approaches intended to probe the molecular hyperspace from different perspectives. This project follows the main tenet that temperature- and concentration-dependent molecular characteristics are to be used for solubility modeling. This initial bias rules out many approaches for molecular descriptor determination, including SMILES-related methods, fingerprints, and also the direct utilization of many first-principle procedures. However, there is a very promising and well-recognized methodology combining both quantum chemistry computations and statistical mechanics to account for environmental influences on the thermodynamic properties of the systems under consideration. This approach, acronymed as COSMO-RS (Conductor-like Screening Model for Realistic Solvents) [120], is designed for bulk systems, in principle, of arbitrary compositions, concentrations, and temperature ranges. That is why this particular framework is selected as a default source of theoretically determined molecular characteristics, as it was already successfully documented in our previous projects [79,83,117,118,121–124]. In particular, the COSMO-RS theory implemented in the COSMOtherm package (version 24.0.0) offers a very convenient way to saturate system characteristics and is used as the source of molecular descriptors also in this project.

The first choice of molecular descriptors is the solubility computed with the aid of COSMOtherm. Unfortunately, this imposes limitations inherent to the approach, and its utilization requires providing the fusion data, such as heat of fusion, temperature of melting, and eventually the heat capacity change upon melting. These experimental parameters are crucial for including the lattice energy. Fortunately, these data are available for both solutes, and several estimates were reported of the melting temperature and heat of fusion. The following mean values were adopted in this work, based on the compilation of data from various sources, as available in the literature [125]: T<sub>m</sub> = 348.64 K and ΔH<sub>m</sub> = 26.69 kJ/mol for ibuprofen and dexibuprofen, T<sub>m</sub> = 367.99 K and ΔH<sub>m</sub> = 25.52 kJ/mol for ketoprofen, and T<sub>m</sub> = 387.34 K and ΔH<sub>m</sub> = 27.86 kJ/mol for flurbiprofen. The COSMOtherm program [126] was used for computing solubility using the newest BP_TZVPD_FINE_24.ctd parametrization. All structures of both solutes and solvent molecules were characterized by conformational diversity derived by comprehensive conformational analysis using the COSMOconf program (version 24.0.0) [127]. Solubilities were computed by fully solving
the Solid-Liquid Equilibria (SLE) problem rather than using an iterative procedure, since it fails for highly soluble cases.

The second type of molecular information was taken from σ-potential. These values were selected rather than σ-profiles since they are temperature-related. For each solute and solvent, the σ-potential profiles were computed for a pure single-component state at given temperature. The set of molecular descriptors used for machine learning was computed as the difference between the σ-potential of a pure solute and the σ-potential of a solvent at a given temperature. In the case of multicomponent solvents, the σ-potential was represented as the value weighted by the mole fraction of the solute-free mole fraction and then used for relative σ-potential profile determination. Typically, the COSMOtherm generates σ-potential profiles in the form of 61 points for σ between −0.03 and +0.03 e/Å² with 0.001 step. To reduce the number of descriptors to the most promising, the relationship between the experimental solubility and values of relative σ-potentials was inspected and quantified by R² values for every σ. This was conducted for two subsets, namely, only non-DES solvents, including neat ones and binary mixtures, and DES-only systems. The obtained results are provided in Figure 7. It documents that there is a modest correlation between experimental solubility and relative σ-potential profiles. However, the regions with the highest contribution are different for both subsets. For non-DES systems, the most significant is the region of σ ∈ [−0.01, +0.01] typically attributed to hydrophobicity (HYD) and, to a lesser extent, the fraction of the hydrogen bond donicity (HBD) defined by σ ∈ [+0.01, +0.03]. From the solubility perspective, the hydrogen acceptability (HBA) σ ∈ [−0.03, −0.01] is marginal in the case of non-DESs. Conversely, for DES systems, both HBD and HBA regions have the highest correlation with experimental solubility, and the contribution of the HYD region is marginal. This suggests that the total solubility dataset is very diverse and covers quite a large molecular structure representation. To include the diversity into the training protocol and, at the same time, reduce the number of data points of the relative σ-potentials, only such values were used which fulfilled the criterion of R² > 0.4 for either subset. This is marked by bold points in Figure 7. Hence, the values of computed σ-potentials within the range of σ between −0.017 and +0.021 with the exclusion of the region between 0.006 and 0.011 were used as molecular descriptors.

<figure><img src="image_10.png" /><figcaption>**Figure 7.** Correlation between relative solute–solvent σ-potentials and experimental solubility, expressed as a function of σ values. Two series represent correlations computed for subsets including only non-DES solvents (neat solvents and binary mixtures) or DES-only systems. Bold symbols define regions used as a source of molecular descriptors, where R² > 0.4 for either subset.</figcaption></figure>
the Solid-Liquid Equilibria (SLE) problem rather than using an iterative procedure, since it fails for highly soluble cases.

The second type of molecular information was taken from σ-potential. These values were selected rather than σ-profiles since they are temperature-related. For each solute and solvent, the σ-potential profiles were computed for a pure single-component state at given temperature. The set of molecular descriptors used for machine learning was computed as the difference between the σ-potential of a pure solute and the σ-potential of a solvent at a given temperature. In the case of multicomponent solvents, the σ-potential was represented as the value weighted by the mole fraction of the solute-free mole fraction and then used for relative σ-potential profile determination. Typically, the COSMOtherm generates σ-potential profiles in the form of 61 points for σ between −0.03 and +0.03 e/Å² with 0.001 step. To reduce the number of descriptors to the most promising, the relationship between the experimental solubility and values of relative σ-potentials was inspected and quantified by R² values for every σ. This was conducted for two subsets, namely, only non-DES solvents, including neat ones and binary mixtures, and DES-only systems. The obtained results are provided in Figure 7. It documents that there is a modest correlation between experimental solubility and relative σ-potential profiles. However, the regions with the highest contribution are different for both subsets. For non-DES systems, the most significant is the region of σ ∈ [−0.01, +0.01] typically attributed to hydrophobicity (HYD) and, to a lesser extent, the fraction of the hydrogen bond donicity (HBD) defined by σ ∈ [+0.01, +0.03]. From the solubility perspective, the hydrogen acceptability (HBA) σ ∈ [−0.03, −0.01] is marginal in the case of non-DESs. Conversely, for DES systems, both HBD and HBA regions have the highest correlation with experimental solubility, and the contribution of the HYD region is marginal. This suggests that the total solubility dataset is very diverse and covers quite a large molecular structure representation. To include the diversity into the training protocol and, at the same time, reduce the number of data points of the relative σ-potentials, only such values were used which fulfilled the criterion of R² > 0.4 for either subset. This is marked by bold points in Figure 7. Hence, the values of computed σ-potentials within the range of σ between −0.017 and +0.021 with the exclusion of the region between 0.006 and 0.011 were used as molecular descriptors.

<figure><img src="image_12.png" /><figcaption>**Figure 7.** Correlation between relative solute–solvent σ-potentials and experimental solubility, expressed as a function of σ values. Two series represent correlations computed for subsets including only non-DES solvents (neat solvents and binary mixtures) or DES-only systems. Bold symbols define regions used as a source of molecular descriptors, where R² > 0.4 for either subset.</figcaption></figure>
**4. Conclusions**

This project experimentally and theoretically screened solvent systems for enhancing the solubility of two COX inhibitors: ibuprofen and ketoprofen. Experimental DES formulations were optimized by varying solvent composition and content. Thirty-six DESs were explored by combining two hydrogen bond acceptors (choline chloride, betaine) with six donors (ethylene glycol, diethylene glycol, triethylene glycol, 1,3-butanediol, 1,2-propanediol, glycerol) at molar ratios of 1:1, 1:2, and 1:4. The most effective systems contained triethylene glycol at an HBA:HBD ratio of 1:2. Aqueous DES mixtures exhibited a pseudo-cosolvency effect, increasing solubility up to 0.2 mole fraction water content. This phenomenon is commonly observed due to water's role in stabilizing solute-rich nanodomains around cholinium cations.

For ibuprofen, the optimal aqueous DES outperformed the most effective organic solvent (chloroform), offering a more sustainable alternative. In contrast, ketoprofen solubility in DES-water was surpassed by methanol. Further DES optimization is warranted.

A machine learning model using COSMO-RS descriptors accurately predicted the solubility for test compounds. This approach enables rational DES formulation design and solubility screening to advance pharmaceutical applications, circumventing laborious experimental screening.

**Supplementary Materials:** The following supporting information can be downloaded at https://www.mdpi.com/article/10.3390/molecules29102296/s1, Table S1: Mole fraction solubility of ibuprofen (IP) in neat DES containing either choline chloride or betaine as an HBA and one of six studied polyols as an HBD in various molar proportions at 25 °C. Standard deviation values are given in parentheses. Table S2: Mole fraction solubility of ketoprofen (KP) in neat DES containing either choline chloride or betaine as an HBA and one of six studied polyols as an HBD in various molar proportions at 25 °C. Standard deviation values are given in parentheses. Table S3: Mole fraction solubility of ibuprofen (IP) in aqueous DES mixtures at various temperatures. In the first column, x\*DES denotes the mole fraction of DES in solute-free mixtures with water. Standard deviation values are given in parentheses. Table S4: Mole fraction solubility of ketoprofen (KP) in aqueous DES mixtures at various temperatures. In the first column, x\*DES denotes the mole fraction of DES in solute-free mixtures with water. Standard deviation values are given in parentheses.

**Author Contributions:** Conceptualization, P.C.; methodology, P.C., T.J., and M.P.; software, P.C.; validation, P.C., T.J., and M.P.; formal analysis, P.C., T.J., A.M., and J.K.; investigation, P.C., T.J., and M.P.; resources, P.C.; data curation, P.C., M.P., and T.J.; writing—original draft preparation, P.C., T.J., and M.P.; writing—review and editing, P.C., T.J., and M.P.; visualization, P.C., and T.J.; supervision, P.C.; project administration, P.C.; funding acquisition, P.C. All authors have read and agreed to the published version of the manuscript.

**Funding:** This research received no external funding.

**Institutional Review Board Statement:** Not applicable.

**Informed Consent Statement:** Not applicable.

**Data Availability Statement:** All data supporting the reported results are available on request from the corresponding author.

**Conflicts of Interest:** The authors declare no conflicts of interest.

**References**

1. Moses, V.S.; Bertone, A.L. Nonsteroidal anti-inflammatory drugs. *Vet. Clin. N. Am.-Equine Pract.* **2002**, *18*, 21–37. [CrossRef]

2. Vane, J.R.; Botting, R.M. Mechanism of Action of Nonsteroidal Anti-inflammatory Drugs. *Am. J. Med.* **1998**, *104*, 2S–8S. [CrossRef] [PubMed]

3. Ghanim, A.M.; Girgis, A.S.; Kariuki, B.M.; Samir, N.; Said, M.F.; Abdelnaser, A.; Nasr, S.; Bekheit, M.S.; Abdelhameed, M.F.; Almalki, A.J.; et al. Design and synthesis of ibuprofen-quinoline conjugates as potential anti-inflammatory and analgesic drug candidates. *Bioorg. Chem.* **2022**, *119*, 105557. [CrossRef] [PubMed]
**4. Conclusions**

This project experimentally and theoretically screened solvent systems for enhancing the solubility of two COX inhibitors: ibuprofen and ketoprofen. Experimental DES formulations were optimized by varying solvent composition and content. Thirty-six DESs were explored by combining two hydrogen bond acceptors (choline chloride, betaine) with six donors (ethylene glycol, diethylene glycol, triethylene glycol, 1,3-butanediol, 1,2-propanediol, glycerol) at molar ratios of 1:1, 1:2, and 1:4. The most effective systems contained triethylene glycol at an HBA:HBD ratio of 1:2. Aqueous DES mixtures exhibited a pseudo-cosolvency effect, increasing solubility up to 0.2 mole fraction water content. This phenomenon is commonly observed due to water's role in stabilizing solute-rich nanodomains around cholinium cations.

For ibuprofen, the optimal aqueous DES outperformed the most effective organic solvent (chloroform), offering a more sustainable alternative. In contrast, ketoprofen solubility in DES-water was surpassed by methanol. Further DES optimization is warranted.

A machine learning model using COSMO-RS descriptors accurately predicted the solubility for test compounds. This approach enables rational DES formulation design and solubility screening to advance pharmaceutical applications, circumventing laborious experimental screening.

**Supplementary Materials:** The following supporting information can be downloaded at https://www.mdpi.com/article/10.3390/molecules29102296/s1, Table S1: Mole fraction solubility of ibuprofen (IP) in neat DES containing either choline chloride or betaine as an HBA and one of six studied polyols as an HBD in various molar proportions at 25 °C. Standard deviation values are given in parentheses. Table S2: Mole fraction solubility of ketoprofen (KP) in neat DES containing either choline chloride or betaine as an HBA and one of six studied polyols as an HBD in various molar proportions at 25 °C. Standard deviation values are given in parentheses. Table S3: Mole fraction solubility of ibuprofen (IP) in aqueous DES mixtures at various temperatures. In the first column, x\*DES denotes the mole fraction of DES in solute-free mixtures with water. Standard deviation values are given in parentheses. Table S4: Mole fraction solubility of ketoprofen (KP) in aqueous DES mixtures at various temperatures. In the first column, x\*DES denotes the mole fraction of DES in solute-free mixtures with water. Standard deviation values are given in parentheses.

**Author Contributions:** Conceptualization, P.C.; methodology, P.C., T.J., and M.P.; software, P.C.; validation, P.C., T.J., and M.P.; formal analysis, P.C., T.J., A.M., and J.K.; investigation, P.C., T.J., and M.P.; resources, P.C.; data curation, P.C., M.P., and T.J.; writing—original draft preparation, P.C., T.J., and M.P.; writing—review and editing, P.C., T.J., and M.P.; visualization, P.C., and T.J.; supervision, P.C.; project administration, P.C.; funding acquisition, P.C. All authors have read and agreed to the published version of the manuscript.

**Funding:** This research received no external funding.

**Institutional Review Board Statement:** Not applicable.

**Informed Consent Statement:** Not applicable.

**Data Availability Statement:** All data supporting the reported results are available on request from the corresponding author.

**Conflicts of Interest:** The authors declare no conflicts of interest.

**References**

1. Moses, V.S.; Bertone, A.L. Nonsteroidal anti-inflammatory drugs. *Vet. Clin. N. Am.-Equine Pract.* **2002**, *18*, 21–37. [CrossRef]

2. Vane, J.R.; Botting, R.M. Mechanism of Action of Nonsteroidal Anti-inflammatory Drugs. *Am. J. Med.* **1998**, *104*, 2S–8S. [CrossRef] [PubMed]

3. Ghanim, A.M.; Girgis, A.S.; Kariuki, B.M.; Samir, N.; Said, M.F.; Abdelnaser, A.; Nasr, S.; Bekheit, M.S.; Abdelhameed, M.F.; Almalki, A.J.; et al. Design and synthesis of ibuprofen-quinoline conjugates as potential anti-inflammatory and analgesic drug candidates. *Bioorg. Chem.* **2022**, *119*, 105557. [CrossRef] [PubMed]
4. Maleškić Kapo, S.; Rakanović-Todić, M.; Burnazović-Ristić, L.; Kusturica, J.; Kulo Ćesić, A.; Ademović, E.; Loga-Zec, S.; Sarač-Hadžihalilović, A.; Aganović-Mušinović, I. Analgesic and anti-inflammatory effects of diclofenac and ketoprofen patches in two different rat models of acute inflammation. *J. King Saud Univ.-Sci.* **2023**, *35*, 102394. [CrossRef]

5. Wang, Y.; Han, Q.; Zhang, H.; Yan, Y. Evaluation of the binding interactions of p-acetylaminophenol, aspirin, ibuprofen and aminopyrine with norfloxacin from the view of antipyretic and anti-inflammatory. *J. Mol. Liq.* **2020**, *312*, 113397. [CrossRef]

6. Cabré, F.; Fernández, M.F.; Calvo, L.; Ferrer, X.; García, M.L.; Mauleón, D. Analgesic, Antiinflammatory, and Antipyretic Effects of S(+)-Ketoprofen In Vivo. *J. Clin. Pharmacol.* **1998**, *38*, 3S–10S. [CrossRef] [PubMed]

7. Panda, S.S.; Girgis, A.S.; Honkanadavar, H.H.; George, R.F.; Srour, A.M. Synthesis of new ibuprofen hybrid conjugates as potential anti-inflammatory and analgesic agents. *Future Med. Chem.* **2020**, *12*, 1369–1386. [CrossRef] [PubMed]

8. Espinosa-Cano, E.; Aguilar, M.R.; Portilla, Y.; Barber, D.F.; Román, J.S. Anti-Inflammatory Polymeric Nanoparticles Based on Ketoprofen and Dexamethasone. *Pharmaceutics* **2020**, *12*, 723. [CrossRef]

9. Næsdal, J.; Brown, K. NSAID-associated adverse effects and acid control aids to prevent them: A review of current treatment options. *Drug Saf.* **2006**, *29*, 119–132. [CrossRef]

10. Xia, W.; Zhang, S.; Li, Y.; Fan, J.; Liu, B.; Wang, L.; Peng, X. Ibuprofen-derived fluorescence inhibitor of COX-2 for breast cancer imaging, prevention and treatment. *Dye. Pigment* **2021**, *190*, 109326. [CrossRef]

11. Romero-Chávez, M.M.; Pineda-Urbina, K.; Pérez, D.J.; Obledo-Benicio, F.; Flores-Parra, A.; Gómez-Sandoval, Z.; Ramos-Organillo, Á. Organotin(IV) compounds derived from ibuprofen and cinnamic acids, an alternative into design of anti-inflammatory by the cyclooxygenases (COX-1 and COX-2) pathway. *J. Organomet. Chem.* **2018**, *862*, 58–70. [CrossRef]

12. Levoin, N.; Chrétien, F.; Lapicque, F.; Chapleur, Y. Synthesis and Biological Testing of Acyl-CoA–Ketoprofen Conjugates as Selective Irreversible Inhibitors of COX-2. *Bioorg. Med. Chem.* **2002**, *10*, 753–757. [CrossRef] [PubMed]

13. Knych, H.K.; McKemie, D.S.; Kass, P.H.; Stanley, S.D.; Blea, J. Ketoprofen in horses: Metabolism, pharmacokinetics, and effects on inflammatory biomarkers. *Drug Test. Anal.* **2024**, *16*, 289–302. [CrossRef] [PubMed]

14. Gaebler, A.J.; Schoretsanitis, G.; Ben Omar, N.; Haen, E.; Endres, K.; Hiemke, C.; Paulzen, M. Metamizole but not ibuprofen reduces the plasma concentration of sertraline: Implications for the concurrent treatment of pain and depression/anxiety disorders. *Br. J. Clin. Pharmacol.* **2021**, *87*, 1111–1119. [CrossRef] [PubMed]

15. Tiwari, S.; Mishra, M.; Salemi, M.R.; Phinney, B.S.; Newens, J.L.; Gomes, A.V. Gender-specific changes in energy metabolism and protein degradation as major pathways affected in livers of mice treated with ibuprofen. *Sci. Rep.* **2020**, *10*, 3386. [CrossRef]

16. Bettiol, A.; Marconi, E.; Vannacci, A.; Simonetti, M.; Magni, A.; Cricelli, C.; Lapi, F. Effectiveness of ibuprofen plus paracetamol combination on persistence of acute musculoskeletal disorders in primary care patients. *Int. J. Clin. Pharm.* **2021**, *43*, 1045–1054. [CrossRef] [PubMed]

17. Iolascon, G.; Giménez, S.; Mogyorósi, D. A review of aceclofenac: Analgesic and anti-inflammatory effects on musculoskeletal disorders. *J. Pain Res.* **2021**, *14*, 3651–3663. [CrossRef]

18. Yılmaz, M.; Gürses, D.; Tükenmez, G. The effectiveness and safety of ibuprofen and acetylsalicylic acid in acute rheumatic fever. *Pediatr. Int.* **2022**, *64*, e15133. [CrossRef]

19. Atzeni, F.; Masala, I.F.; Bagnasco, M.; Lanata, L.; Mantelli, F.; Sarzi-Puttini, P. Comparison of Efficacy of Ketoprofen and Ibuprofen in Treating Pain in Patients with Rheumatoid Arthritis: A Systematic Review and Meta-Analysis. *Pain Ther.* **2021**, *10*, 577–588. [CrossRef]

20. Wen, M.M.; Ismail, N.I.K.; Nasra, M.M.A.; El-Kamel, A.H. Repurposing ibuprofen-loaded microemulsion for the management of Alzheimer's disease: Evidence of potential intranasal brain targeting. *Drug Deliv.* **2021**, *28*, 1188–1203. [CrossRef]

21. Singh, A.; Tripathi, P.; Singh, S. Neuroinflammatory responses in Parkinson's disease: Relevance of Ibuprofen in therapeutics. *Inflammopharmacology* **2021**, *29*, 5–14. [CrossRef]

22. Mendonça, L.S.; Nóbrega, C.; Tavino, S.; Brinkhaus, M.; Matos, C.; Tomé, S.; Moreira, R.; Henriques, D.; Kaspar, B.K.; Pereira De Almeida, L. Ibuprofen enhances synaptic function and neural progenitors proliferation markers and improves neuropathology and motor coordination in Machado–Joseph disease models. *Hum. Mol. Genet.* **2019**, *28*, 3691–3703. [CrossRef] [PubMed]

23. Varrassi, G.; Pergolizzi, J.V.; Dowling, P.; Paladini, A. Ibuprofen Safety at the Golden Anniversary: Are all NSAIDs the Same? A Narrative Review. *Adv. Ther.* **2020**, *37*, 61–82. [CrossRef]

24. Zoubek, M.E.; Lucena, M.I.; Andrade, R.J.; Stephens, C. Systematic review: Ibuprofen-induced liver injury. *Aliment. Pharmacol. Ther.* **2020**, *51*, 603–611. [CrossRef]

25. Ibraimi, Q.; Driton Selmani, D.; Arjeta Shabani, D. The safety of ketoprofen usage in different age. *Int. J. Med. Healthc.* **2023**, *8*, 9–14.

26. Cumhur Cure, M.; Kucuk, A.; Cure, E. NSAIDs may increase the risk of thrombosis and acute renal failure in patients with COVID-19 infection. *Therapie* **2020**, *75*, 387. [CrossRef]

27. Schjerning, A.M.; McGettigan, P.; Gislason, G. Cardiovascular effects and safety of (non-aspirin) NSAIDs. *Nat. Rev. Cardiol.* **2020**, *17*, 574–584. [CrossRef]

28. Martínez, F.; Jouyban, A.; Acree, W.E. Pharmaceuticals solubility is still nowadays widely studied everywhere. *Pharm. Sci.* **2017**, *23*, 1–2. [CrossRef]

29. Savjani, K.T.; Gajjar, A.K.; Savjani, J.K. Drug Solubility: Importance and Enhancement Techniques. *ISRN Pharm.* **2012**, *2012*, 195727. [CrossRef]
4. Maleškić Kapo, S.; Rakanović-Todić, M.; Burnazović-Ristić, L.; Kusturica, J.; Kulo Ćesić, A.; Ademović, E.; Loga-Zec, S.; Sarač-Hadžihalilović, A.; Aganović-Mušinović, I. Analgesic and anti-inflammatory effects of diclofenac and ketoprofen patches in two different rat models of acute inflammation. *J. King Saud Univ.-Sci.* **2023**, *35*, 102394. [CrossRef]

5. Wang, Y.; Han, Q.; Zhang, H.; Yan, Y. Evaluation of the binding interactions of p-acetylaminophenol, aspirin, ibuprofen and aminopyrine with norfloxacin from the view of antipyretic and anti-inflammatory. *J. Mol. Liq.* **2020**, *312*, 113397. [CrossRef]

6. Cabré, F.; Fernández, M.F.; Calvo, L.; Ferrer, X.; García, M.L.; Mauleón, D. Analgesic, Antiinflammatory, and Antipyretic Effects of S(+)-Ketoprofen In Vivo. *J. Clin. Pharmacol.* **1998**, *38*, 3S–10S. [CrossRef] [PubMed]

7. Panda, S.S.; Girgis, A.S.; Honkanadavar, H.H.; George, R.F.; Srour, A.M. Synthesis of new ibuprofen hybrid conjugates as potential anti-inflammatory and analgesic agents. *Future Med. Chem.* **2020**, *12*, 1369–1386. [CrossRef] [PubMed]

8. Espinosa-Cano, E.; Aguilar, M.R.; Portilla, Y.; Barber, D.F.; Román, J.S. Anti-Inflammatory Polymeric Nanoparticles Based on Ketoprofen and Dexamethasone. *Pharmaceutics* **2020**, *12*, 723. [CrossRef]

9. Næsdal, J.; Brown, K. NSAID-associated adverse effects and acid control aids to prevent them: A review of current treatment options. *Drug Saf.* **2006**, *29*, 119–132. [CrossRef]

10. Xia, W.; Zhang, S.; Li, Y.; Fan, J.; Liu, B.; Wang, L.; Peng, X. Ibuprofen-derived fluorescence inhibitor of COX-2 for breast cancer imaging, prevention and treatment. *Dye. Pigment* **2021**, *190*, 109326. [CrossRef]

11. Romero-Chávez, M.M.; Pineda-Urbina, K.; Pérez, D.J.; Obledo-Benicio, F.; Flores-Parra, A.; Gómez-Sandoval, Z.; Ramos-Organillo, Á. Organotin(IV) compounds derived from ibuprofen and cinnamic acids, an alternative into design of anti-inflammatory by the cyclooxygenases (COX-1 and COX-2) pathway. *J. Organomet. Chem.* **2018**, *862*, 58–70. [CrossRef]

12. Levoin, N.; Chrétien, F.; Lapicque, F.; Chapleur, Y. Synthesis and Biological Testing of Acyl-CoA–Ketoprofen Conjugates as Selective Irreversible Inhibitors of COX-2. *Bioorg. Med. Chem.* **2002**, *10*, 753–757. [CrossRef] [PubMed]

13. Knych, H.K.; McKemie, D.S.; Kass, P.H.; Stanley, S.D.; Blea, J. Ketoprofen in horses: Metabolism, pharmacokinetics, and effects on inflammatory biomarkers. *Drug Test. Anal.* **2024**, *16*, 289–302. [CrossRef] [PubMed]

14. Gaebler, A.J.; Schoretsanitis, G.; Ben Omar, N.; Haen, E.; Endres, K.; Hiemke, C.; Paulzen, M. Metamizole but not ibuprofen reduces the plasma concentration of sertraline: Implications for the concurrent treatment of pain and depression/anxiety disorders. *Br. J. Clin. Pharmacol.* **2021**, *87*, 1111–1119. [CrossRef] [PubMed]

15. Tiwari, S.; Mishra, M.; Salemi, M.R.; Phinney, B.S.; Newens, J.L.; Gomes, A.V. Gender-specific changes in energy metabolism and protein degradation as major pathways affected in livers of mice treated with ibuprofen. *Sci. Rep.* **2020**, *10*, 3386. [CrossRef]

16. Bettiol, A.; Marconi, E.; Vannacci, A.; Simonetti, M.; Magni, A.; Cricelli, C.; Lapi, F. Effectiveness of ibuprofen plus paracetamol combination on persistence of acute musculoskeletal disorders in primary care patients. *Int. J. Clin. Pharm.* **2021**, *43*, 1045–1054. [CrossRef] [PubMed]

17. Iolascon, G.; Giménez, S.; Mogyorósi, D. A review of aceclofenac: Analgesic and anti-inflammatory effects on musculoskeletal disorders. *J. Pain Res.* **2021**, *14*, 3651–3663. [CrossRef]

18. Yılmaz, M.; Gürses, D.; Tükenmez, G. The effectiveness and safety of ibuprofen and acetylsalicylic acid in acute rheumatic fever. *Pediatr. Int.* **2022**, *64*, e15133. [CrossRef]

19. Atzeni, F.; Masala, I.F.; Bagnasco, M.; Lanata, L.; Mantelli, F.; Sarzi-Puttini, P. Comparison of Efficacy of Ketoprofen and Ibuprofen in Treating Pain in Patients with Rheumatoid Arthritis: A Systematic Review and Meta-Analysis. *Pain Ther.* **2021**, *10*, 577–588. [CrossRef]

20. Wen, M.M.; Ismail, N.I.K.; Nasra, M.M.A.; El-Kamel, A.H. Repurposing ibuprofen-loaded microemulsion for the management of Alzheimer's disease: Evidence of potential intranasal brain targeting. *Drug Deliv.* **2021**, *28*, 1188–1203. [CrossRef]

21. Singh, A.; Tripathi, P.; Singh, S. Neuroinflammatory responses in Parkinson's disease: Relevance of Ibuprofen in therapeutics. *Inflammopharmacology* **2021**, *29*, 5–14. [CrossRef]

22. Mendonça, L.S.; Nóbrega, C.; Tavino, S.; Brinkhaus, M.; Matos, C.; Tomé, S.; Moreira, R.; Henriques, D.; Kaspar, B.K.; Pereira De Almeida, L. Ibuprofen enhances synaptic function and neural progenitors proliferation markers and improves neuropathology and motor coordination in Machado–Joseph disease models. *Hum. Mol. Genet.* **2019**, *28*, 3691–3703. [CrossRef] [PubMed]

23. Varrassi, G.; Pergolizzi, J.V.; Dowling, P.; Paladini, A. Ibuprofen Safety at the Golden Anniversary: Are all NSAIDs the Same? A Narrative Review. *Adv. Ther.* **2020**, *37*, 61–82. [CrossRef]

24. Zoubek, M.E.; Lucena, M.I.; Andrade, R.J.; Stephens, C. Systematic review: Ibuprofen-induced liver injury. *Aliment. Pharmacol. Ther.* **2020**, *51*, 603–611. [CrossRef]

25. Ibraimi, Q.; Driton Selmani, D.; Arjeta Shabani, D. The safety of ketoprofen usage in different age. *Int. J. Med. Healthc.* **2023**, *8*, 9–14.

26. Cumhur Cure, M.; Kucuk, A.; Cure, E. NSAIDs may increase the risk of thrombosis and acute renal failure in patients with COVID-19 infection. *Therapie* **2020**, *75*, 387. [CrossRef]

27. Schjerning, A.M.; McGettigan, P.; Gislason, G. Cardiovascular effects and safety of (non-aspirin) NSAIDs. *Nat. Rev. Cardiol.* **2020**, *17*, 574–584. [CrossRef]

28. Martínez, F.; Jouyban, A.; Acree, W.E. Pharmaceuticals solubility is still nowadays widely studied everywhere. *Pharm. Sci.* **2017**, *23*, 1–2. [CrossRef]

29. Savjani, K.T.; Gajjar, A.K.; Savjani, J.K. Drug Solubility: Importance and Enhancement Techniques. *ISRN Pharm.* **2012**, *2012*, 195727. [CrossRef]
30. Coltescu, A.R.; Butnariu, M.; Sarac, I. The importance of solubility for new drug molecules. *Biomed. Pharmacol. J.* **2020**, *13*, 577–583.
[CrossRef]

31. Yang, Z.; Yang, Y.; Xia, M.; Dai, W.; Zhu, B.; Mei, X. Improving the dissolution behaviors and bioavailability of abiraterone acetate via multicomponent crystal forms. *Int. J. Pharm.* **2022**, *614*, 121460. [CrossRef] [PubMed]

32. Kalam, M.A.; Alshamsan, A.; Alkholief, M.; Alsarra, I.A.; Ali, R.; Haq, N.; Anwer, M.K.; Shakeel, F. Solubility Measurement and Various Solubility Parameters of Glipizide in Different Neat Solvents. *ACS Omega* **2020**, *5*, 1708–1716. [CrossRef] [PubMed]

33. Kim, H.-S.; Kim, C.-M.; Jo, A.-N.; Kim, J.-E. Studies on Preformulation and Formulation of JIN-001 Liquisolid Tablet with Enhanced Solubility. *Pharmaceuticals* **2022**, *15*, 412. [CrossRef] [PubMed]

34. Khadka, P.; Ro, J.; Kim, H.; Kim, I.; Kim, J.T.; Kim, H.; Cho, J.M.; Yun, G.; Lee, J. Pharmaceutical particle technologies: An approach to improve drug solubility, dissolution and bioavailability. *Asian J. Pharm. Sci.* **2014**, *9*, 304–316. [CrossRef]

35. Müller, C.E. Prodrug Approaches for Enhancing the Bioavailability of Drugs with Low Solubility. *Chem. Biodivers.* **2009**, *6*, 2071–2083. [CrossRef] [PubMed]

36. Das, T.; Mehta, C.H.; Nayak, U.Y. Multiple approaches for achieving drug solubility: An in silico perspective. *Drug Discov. Today* **2020**, *25*, 1206–1212. [CrossRef] [PubMed]

37. Tian, Y.; Shi, C.; Sun, Y.; Zhu, C.; Sun, C.C.; Mao, S. Designing micellar Nanocarriers with improved drug loading and stability based on solubility parameter. *Mol. Pharm.* **2015**, *12*, 816–825. [CrossRef] [PubMed]

38. Lipinski, C.A. Drug-like properties and the causes of poor solubility and poor permeability. *J. Pharmacol. Toxicol. Methods* **2000**, *44*, 235–249. [CrossRef]

39. Da Silva, F.L.O.; Marques, M.B.D.F.; Kato, K.C.; Carneiro, G. Nanonization techniques to overcome poor water-solubility with drugs. *Expert Opin. Drug Discov.* **2020**, *15*, 853–864. [CrossRef]

40. Das, B.; Baidya, A.T.K.; Mathew, A.T.; Yadav, A.K.; Kumar, R. Structural modification aimed for improving solubility of lead compounds in early phase drug discovery. *Bioorg. Med. Chem.* **2022**, *56*, 116614. [CrossRef]

41. Bergström, C.A.S.; Avdeef, A. Perspectives in solubility measurement and interpretation. *ADMET DMPK* **2019**, *7*, 88–105. [CrossRef]

42. Black, S.; Dang, L.; Liu, C.; Wei, H. On the measurement of solubility. *Org. Process Res. Dev.* **2013**, *17*, 486–492. [CrossRef]

43. Boobier, S.; Hose, D.R.J.; Blacker, A.J.; Nguyen, B.N. Machine learning with physicochemical relationships: Solubility prediction in organic solvents and water. *Nat. Commun.* **2020**, *11*, 5753. [CrossRef]

44. Lovrić, M.; Pavlović, K.; Žuvela, P.; Spataru, A.; Lučić, B.; Kern, R.; Wong, M.W. Machine learning in prediction of intrinsic aqueous solubility of drug-like compounds: Generalization, complexity, or predictive ability? *J. Chemom.* **2021**, *35*, e3349. [CrossRef]

45. Bhalani, D.V.; Nutan, B.; Kumar, A.; Singh Chandel, A.K. Bioavailability Enhancement Techniques for Poorly Aqueous Soluble Drugs and Therapeutics. *Biomedicines* **2022**, *10*, 2055. [CrossRef]

46. Manallack, D.T.; Yuriev, E.; Chalmers, D.K. The influence and manipulation of acid/base properties in drug discovery. *Drug Discov. Today Technol.* **2018**, *27*, 41–47. [CrossRef]

47. Merisko-Liversidge, E.; Liversidge, G.G. Nanosizing for oral and parenteral drug delivery: A perspective on formulating poorly-water soluble compounds using wet media milling technology. *Adv. Drug Deliv. Rev.* **2011**, *63*, 427–440. [CrossRef]

48. Brewster, M.E.; Loftsson, T. Cyclodextrins as pharmaceutical solubilizers. *Adv. Drug Deliv. Rev.* **2007**, *59*, 645–666. [CrossRef]

49. Korn, C.; Balbach, S. Compound selection for development—Is salt formation the ultimate answer? Experiences with an extended concept of the "100mg approach". *Eur. J. Pharm. Sci.* **2014**, *57*, 257–263. [CrossRef]

50. Seedher, N.; Kanojia, M. Co-solvent solubilization of some poorly-soluble antidiabetic drugs. *Pharm. Dev. Technol.* **2009**, *14*, 185–192. [CrossRef]

51. Hahnenkamp, I.; Graubner, G.; Gmehling, J. Measurement and prediction of solubilities of active pharmaceutical ingredients. *Int. J. Pharm.* **2010**, *388*, 73–81. [CrossRef] [PubMed]

52. Abraham, M.H.; Smith, R.E.; Luchtefeld, R.; Boorem, A.J.; Lou, R.; Acree, W.E. Prediction of solubility of drugs and other compounds in organic solvents. *J. Pharm. Sci.* **2010**, *99*, 1500–1515. [CrossRef]

53. Hewitt, M.; Cronin, M.T.D.; Enoch, S.J.; Madden, J.C.; Roberts, D.W.; Dearden, J.C. In silico prediction of aqueous solubility: The solubility challenge. *J. Chem. Inf. Model.* **2009**, *49*, 2572–2587. [CrossRef]

54. Tosca, E.M.; Bartolucci, R.; Magni, P. Application of artificial neural networks to predict the intrinsic solubility of drug-like molecules. *Pharmaceutics* **2021**, *13*, 1101. [CrossRef]

55. Deng, T.; Jia, G.z. Prediction of aqueous solubility of compounds based on neural network. *Mol. Phys.* **2020**, *118*, e1600754. [CrossRef]

56. Wesolowski, M.; Suchacz, B. Artificial Neural Networks: Theoretical Background and Pharmaceutical Applications: A Review. *J. AOAC Int.* **2012**, *95*, 652–668. [CrossRef]

57. Agatonovic-Kustrin, S.; Beresford, R. Basic concepts of artificial neural network (ANN) modeling and its application in pharmaceutical research. *J. Pharm. Biomed. Anal.* **2000**, *22*, 717–727. [CrossRef]

58. Lenoir, D.; Schramm, K.W.; Lalah, J.O. Green Chemistry: Some important forerunners and current issues. *Sustain. Chem. Pharm.* **2020**, *18*, 100313. [CrossRef]

59. Kopach, M.; Leahy, D.; Manley, J. The green chemistry approach to pharma manufacturing. *Innov. Pharm. Technol.* **2012**, *43*, 72–75.
30. Coltescu, A.R.; Butnariu, M.; Sarac, I. The importance of solubility for new drug molecules. *Biomed. Pharmacol. J.* **2020**, *13*, 577–583.
[CrossRef]

31. Yang, Z.; Yang, Y.; Xia, M.; Dai, W.; Zhu, B.; Mei, X. Improving the dissolution behaviors and bioavailability of abiraterone acetate via multicomponent crystal forms. *Int. J. Pharm.* **2022**, *614*, 121460. [CrossRef] [PubMed]

32. Kalam, M.A.; Alshamsan, A.; Alkholief, M.; Alsarra, I.A.; Ali, R.; Haq, N.; Anwer, M.K.; Shakeel, F. Solubility Measurement and Various Solubility Parameters of Glipizide in Different Neat Solvents. *ACS Omega* **2020**, *5*, 1708–1716. [CrossRef] [PubMed]

33. Kim, H.-S.; Kim, C.-M.; Jo, A.-N.; Kim, J.-E. Studies on Preformulation and Formulation of JIN-001 Liquisolid Tablet with Enhanced Solubility. *Pharmaceuticals* **2022**, *15*, 412. [CrossRef] [PubMed]

34. Khadka, P.; Ro, J.; Kim, H.; Kim, I.; Kim, J.T.; Kim, H.; Cho, J.M.; Yun, G.; Lee, J. Pharmaceutical particle technologies: An approach to improve drug solubility, dissolution and bioavailability. *Asian J. Pharm. Sci.* **2014**, *9*, 304–316. [CrossRef]

35. Müller, C.E. Prodrug Approaches for Enhancing the Bioavailability of Drugs with Low Solubility. *Chem. Biodivers.* **2009**, *6*, 2071–2083. [CrossRef] [PubMed]

36. Das, T.; Mehta, C.H.; Nayak, U.Y. Multiple approaches for achieving drug solubility: An in silico perspective. *Drug Discov. Today* **2020**, *25*, 1206–1212. [CrossRef] [PubMed]

37. Tian, Y.; Shi, C.; Sun, Y.; Zhu, C.; Sun, C.C.; Mao, S. Designing micellar Nanocarriers with improved drug loading and stability based on solubility parameter. *Mol. Pharm.* **2015**, *12*, 816–825. [CrossRef] [PubMed]

38. Lipinski, C.A. Drug-like properties and the causes of poor solubility and poor permeability. *J. Pharmacol. Toxicol. Methods* **2000**, *44*, 235–249. [CrossRef]

39. Da Silva, F.L.O.; Marques, M.B.D.F.; Kato, K.C.; Carneiro, G. Nanonization techniques to overcome poor water-solubility with drugs. *Expert Opin. Drug Discov.* **2020**, *15*, 853–864. [CrossRef]

40. Das, B.; Baidya, A.T.K.; Mathew, A.T.; Yadav, A.K.; Kumar, R. Structural modification aimed for improving solubility of lead compounds in early phase drug discovery. *Bioorg. Med. Chem.* **2022**, *56*, 116614. [CrossRef]

41. Bergström, C.A.S.; Avdeef, A. Perspectives in solubility measurement and interpretation. *ADMET DMPK* **2019**, *7*, 88–105. [CrossRef]

42. Black, S.; Dang, L.; Liu, C.; Wei, H. On the measurement of solubility. *Org. Process Res. Dev.* **2013**, *17*, 486–492. [CrossRef]

43. Boobier, S.; Hose, D.R.J.; Blacker, A.J.; Nguyen, B.N. Machine learning with physicochemical relationships: Solubility prediction in organic solvents and water. *Nat. Commun.* **2020**, *11*, 5753. [CrossRef]

44. Lovrić, M.; Pavlović, K.; Žuvela, P.; Spataru, A.; Lučić, B.; Kern, R.; Wong, M.W. Machine learning in prediction of intrinsic aqueous solubility of drug-like compounds: Generalization, complexity, or predictive ability? *J. Chemom.* **2021**, *35*, e3349. [CrossRef]

45. Bhalani, D.V.; Nutan, B.; Kumar, A.; Singh Chandel, A.K. Bioavailability Enhancement Techniques for Poorly Aqueous Soluble Drugs and Therapeutics. *Biomedicines* **2022**, *10*, 2055. [CrossRef]

46. Manallack, D.T.; Yuriev, E.; Chalmers, D.K. The influence and manipulation of acid/base properties in drug discovery. *Drug Discov. Today Technol.* **2018**, *27*, 41–47. [CrossRef]

47. Merisko-Liversidge, E.; Liversidge, G.G. Nanosizing for oral and parenteral drug delivery: A perspective on formulating poorly-water soluble compounds using wet media milling technology. *Adv. Drug Deliv. Rev.* **2011**, *63*, 427–440. [CrossRef]

48. Brewster, M.E.; Loftsson, T. Cyclodextrins as pharmaceutical solubilizers. *Adv. Drug Deliv. Rev.* **2007**, *59*, 645–666. [CrossRef]

49. Korn, C.; Balbach, S. Compound selection for development—Is salt formation the ultimate answer? Experiences with an extended concept of the "100mg approach". *Eur. J. Pharm. Sci.* **2014**, *57*, 257–263. [CrossRef]

50. Seedher, N.; Kanojia, M. Co-solvent solubilization of some poorly-soluble antidiabetic drugs. *Pharm. Dev. Technol.* **2009**, *14*, 185–192. [CrossRef]

51. Hahnenkamp, I.; Graubner, G.; Gmehling, J. Measurement and prediction of solubilities of active pharmaceutical ingredients. *Int. J. Pharm.* **2010**, *388*, 73–81. [CrossRef] [PubMed]

52. Abraham, M.H.; Smith, R.E.; Luchtefeld, R.; Boorem, A.J.; Lou, R.; Acree, W.E. Prediction of solubility of drugs and other compounds in organic solvents. *J. Pharm. Sci.* **2010**, *99*, 1500–1515. [CrossRef]

53. Hewitt, M.; Cronin, M.T.D.; Enoch, S.J.; Madden, J.C.; Roberts, D.W.; Dearden, J.C. In silico prediction of aqueous solubility: The solubility challenge. *J. Chem. Inf. Model.* **2009**, *49*, 2572–2587. [CrossRef]

54. Tosca, E.M.; Bartolucci, R.; Magni, P. Application of artificial neural networks to predict the intrinsic solubility of drug-like molecules. *Pharmaceutics* **2021**, *13*, 1101. [CrossRef]

55. Deng, T.; Jia, G.z. Prediction of aqueous solubility of compounds based on neural network. *Mol. Phys.* **2020**, *118*, e1600754. [CrossRef]

56. Wesolowski, M.; Suchacz, B. Artificial Neural Networks: Theoretical Background and Pharmaceutical Applications: A Review. *J. AOAC Int.* **2012**, *95*, 652–668. [CrossRef]

57. Agatonovic-Kustrin, S.; Beresford, R. Basic concepts of artificial neural network (ANN) modeling and its application in pharmaceutical research. *J. Pharm. Biomed. Anal.* **2000**, *22*, 717–727. [CrossRef]

58. Lenoir, D.; Schramm, K.W.; Lalah, J.O. Green Chemistry: Some important forerunners and current issues. *Sustain. Chem. Pharm.* **2020**, *18*, 100313. [CrossRef]

59. Kopach, M.; Leahy, D.; Manley, J. The green chemistry approach to pharma manufacturing. *Innov. Pharm. Technol.* **2012**, *43*, 72–75.
60. Becker, J.; Manske, C.; Randl, S. Green chemistry and sustainability metrics in the pharmaceutical manufacturing sector. *Curr. Opin. Green Sustain. Chem.* **2022**, *33*, 100562. [CrossRef]

61. Mishra, M.; Sharma, M.; Dubey, R.; Kumari, P.; Ranjan, V.; Pandey, J. Green synthesis interventions of pharmaceutical industries for sustainable development. *Curr. Res. Green Sustain. Chem.* **2021**, *4*, 100174. [CrossRef]

62. DeSimone, J.M. Practical approaches to green solvents. *Science* **2002**, *297*, 799–803. [CrossRef]

63. Häckl, K.; Kunz, W. Some aspects of green solvents. *Comptes Rendus Chim.* **2018**, *21*, 572–580. [CrossRef]

64. Santana-Mayor, Á.; Rodríguez-Ramos, R.; Herrera-Herrera, A.V.; Socas-Rodríguez, B.; Rodríguez-Delgado, M.Á. Deep eutectic solvents. The new generation of green solvents in analytical chemistry. *TrAC Trends Anal. Chem.* **2021**, *134*, 116108. [CrossRef]

65. Vanda, H.; Dai, Y.; Wilson, E.G.; Verpoorte, R.; Choi, Y.H. Green solvents from ionic liquids and deep eutectic solvents to natural deep eutectic solvents. *Comptes Rendus Chim.* **2018**, *21*, 628–638. [CrossRef]

66. Smith, E.L.; Abbott, A.P.; Ryder, K.S. Deep Eutectic Solvents (DESs) and Their Applications. *Chem. Rev.* **2014**, *114*, 11060–11082. [CrossRef]

67. Ruß, C.; König, B. Low melting mixtures in organic synthesis—An alternative to ionic liquids? *Green Chem.* **2012**, *14*, 2969–2982. [CrossRef]

68. Francisco, M.; Van Den Bruinhorst, A.; Kroon, M.C. Minireview. *Angew. Chemie Int. Ed.* **2013**, *11*, 3074–3085. [CrossRef]

69. Omar, K.A.; Sadeghi, R. Physicochemical properties of deep eutectic solvents: A review. *J. Mol. Liq.* **2022**, *360*, 119524. [CrossRef]

70. Paiva, A.; Craveiro, R.; Aroso, I.; Martins, M.; Reis, R.L.; Duarte, A.R.C. Natural Deep Eutectic Solvents–Solvents for the 21st Century. *ACS Sustain. Chem. Eng.* **2014**, *2*, 1063–1071. [CrossRef]

71. Espino, M.; de los Ángeles Fernández, M.; Gomez, F.J.V.; Silva, M.F. Natural designer solvents for greening analytical chemistry. *TrAC Trends Anal. Chem.* **2016**, *76*, 126–136. [CrossRef]

72. Xu, G.; Shi, M.; Zhang, P.; Tu, Z.; Hu, X.; Zhang, X.; Wu, Y. Tuning the composition of deep eutectic solvents consisting of tetrabutylammonium chloride and n-decanoic acid for adjustable separation of ethylene and ethane. *Sep. Purif. Technol.* **2022**, *298*, 121680. [CrossRef]

73. Cao, Y.; Tao, X.; Jiang, S.; Gao, N.; Sun, Z. Tuning thermodynamic properties of deep eutectic solvents for achieving highly efficient photothermal sensor. *J. Mol. Liq.* **2020**, *308*, 113163. [CrossRef]

74. Lomba, L.; Ribate, M.P.; Zaragoza, E.; Concha, J.; Garralaga, M.P.; Errazquin, D.; García, C.B.; Giner, B. Deep Eutectic Solvents: Are They Safe? *Appl. Sci.* **2021**, *11*, 10061. [CrossRef]

75. De Morais, P.; Gonçalves, F.; Coutinho, J.A.P.; Ventura, S.P.M. Ecotoxicity of Cholinium-Based Deep Eutectic Solvents. *ACS Sustain. Chem. Eng.* **2015**, *3*, 3398–3404. [CrossRef]

76. Macário, I.P.E.; Jesus, F.; Pereira, J.L.; Ventura, S.P.M.; Gonçalves, A.M.M.; Coutinho, J.A.P.; Gonçalves, F.J.M. Unraveling the ecotoxicity of deep eutectic solvents using the mixture toxicity theory. *Chemosphere* **2018**, *212*, 890–897. [CrossRef] [PubMed]

77. Nejrotti, S.; Antenucci, A.; Pontremoli, C.; Gontrani, L.; Barbero, N.; Carbone, M.; Bonomo, M. Critical Assessment of the Sustainability of Deep Eutectic Solvents: A Case Study on Six Choline Chloride-Based Mixtures. *ACS Omega* **2022**, *7*, 47449–47461. [CrossRef] [PubMed]

78. Pedro, S.N.; Freire, C.S.R.; Silvestre, A.J.D.; Freire, M.G. Deep Eutectic Solvents and Pharmaceuticals. *Encyclopedia 2021*, *1*, 942–963. [CrossRef]

79. Cysewski, P.; Jeliński, T.; Przybyłek, M. Experimental and Theoretical Insights into the Intermolecular Interactions in Saturated Systems of Dapsone in Conventional and Deep Eutectic Solvents. *Molecules* **2024**, *29*, 1743. [CrossRef] [PubMed]

80. Mustafa, N.R.; Spelbos, V.S.; Witkamp, G.J.; Verpoorte, R.; Choi, Y.H. Solubility and stability of some pharmaceuticals in natural deep eutectic solvents-based formulations. *Molecules* **2021**, *26*, 2645. [CrossRef]

81. Bazzo, G.C.; Pezzini, B.R.; Stulzer, H.K. Eutectic mixtures as an approach to enhance solubility, dissolution rate and oral bioavailability of poorly water-soluble drugs. *Int. J. Pharm.* **2020**, *588*, 119741. [CrossRef] [PubMed]

82. Kapre, S.; Palakurthi, S.S.; Jain, A.; Palakurthi, S. DES-igning the future of drug delivery: A journey from fundamentals to drug delivery applications. *J. Mol. Liq.* **2024**, *400*, 124517. [CrossRef]

83. Jeliński, T.; Przybyłek, M.; Mianowana, M.; Misiak, K.; Cysewski, P. Deep Eutectic Solvents as Agents for Improving the Solubility of Edaravone: Experimental and Theoretical Considerations. *Molecules* **2024**, *29*, 1261. [CrossRef] [PubMed]

84. Duarte, A.R.C.; Ferreira, A.S.D.; Barreiros, S.; Cabrita, E.; Reis, R.L.; Paiva, A. A comparison between pure active pharmaceutical ingredients and therapeutic deep eutectic solvents: Solubility and permeability studies. *Eur. J. Pharm. Biopharm.* **2017**, *114*, 296–304. [CrossRef]

85. Nguyen, C.-H.; Augis, L.; Fourmentin, S.; Barratt, G.; Legrand, F.-X. Deep Eutectic Solvents for Innovative Pharmaceutical Formulations. In *Deep Eutectic Solvents for Medicine, Gas Solubilization and Extraction of Natural Substances*; Springer: Berlin/Heidelberg, Germany, 2021; pp. 41–102.

86. Liu, Y.; Wu, Y.; Liu, J.; Wang, W.; Yang, Q.; Yang, G. Deep eutectic solvents: Recent advances in fabrication approaches and pharmaceutical applications. *Int. J. Pharm.* **2022**, *622*, 121811. [CrossRef] [PubMed]

87. Emami, S.; Shayanfar, A. Deep eutectic solvents for pharmaceutical formulation and drug delivery applications. *Pharm. Dev. Technol.* **2020**, *25*, 779–796. [CrossRef]

88. Pedro, S.N.; Freire, M.G.; Freire, C.S.R.; Silvestre, A.J.D. Deep eutectic solvents comprising active pharmaceutical ingredients in the development of drug delivery systems. *Expert Opin. Drug Deliv.* **2019**, *16*, 497–506. [CrossRef]
60. Becker, J.; Manske, C.; Randl, S. Green chemistry and sustainability metrics in the pharmaceutical manufacturing sector. *Curr. Opin. Green Sustain. Chem.* **2022**, *33*, 100562. [CrossRef]

61. Mishra, M.; Sharma, M.; Dubey, R.; Kumari, P.; Ranjan, V.; Pandey, J. Green synthesis interventions of pharmaceutical industries for sustainable development. *Curr. Res. Green Sustain. Chem.* **2021**, *4*, 100174. [CrossRef]

62. DeSimone, J.M. Practical approaches to green solvents. *Science* **2002**, *297*, 799–803. [CrossRef]

63. Häckl, K.; Kunz, W. Some aspects of green solvents. *Comptes Rendus Chim.* **2018**, *21*, 572–580. [CrossRef]

64. Santana-Mayor, Á.; Rodríguez-Ramos, R.; Herrera-Herrera, A.V.; Socas-Rodríguez, B.; Rodríguez-Delgado, M.Á. Deep eutectic solvents. The new generation of green solvents in analytical chemistry. *TrAC Trends Anal. Chem.* **2021**, *134*, 116108. [CrossRef]

65. Vanda, H.; Dai, Y.; Wilson, E.G.; Verpoorte, R.; Choi, Y.H. Green solvents from ionic liquids and deep eutectic solvents to natural deep eutectic solvents. *Comptes Rendus Chim.* **2018**, *21*, 628–638. [CrossRef]

66. Smith, E.L.; Abbott, A.P.; Ryder, K.S. Deep Eutectic Solvents (DESs) and Their Applications. *Chem. Rev.* **2014**, *114*, 11060–11082. [CrossRef]

67. Ruß, C.; König, B. Low melting mixtures in organic synthesis—An alternative to ionic liquids? *Green Chem.* **2012**, *14*, 2969–2982. [CrossRef]

68. Francisco, M.; Van Den Bruinhorst, A.; Kroon, M.C. Minireview. *Angew. Chemie Int. Ed.* **2013**, *11*, 3074–3085. [CrossRef]

69. Omar, K.A.; Sadeghi, R. Physicochemical properties of deep eutectic solvents: A review. *J. Mol. Liq.* **2022**, *360*, 119524. [CrossRef]

70. Paiva, A.; Craveiro, R.; Aroso, I.; Martins, M.; Reis, R.L.; Duarte, A.R.C. Natural Deep Eutectic Solvents–Solvents for the 21st Century. *ACS Sustain. Chem. Eng.* **2014**, *2*, 1063–1071. [CrossRef]

71. Espino, M.; de los Ángeles Fernández, M.; Gomez, F.J.V.; Silva, M.F. Natural designer solvents for greening analytical chemistry. *TrAC Trends Anal. Chem.* **2016**, *76*, 126–136. [CrossRef]

72. Xu, G.; Shi, M.; Zhang, P.; Tu, Z.; Hu, X.; Zhang, X.; Wu, Y. Tuning the composition of deep eutectic solvents consisting of tetrabutylammonium chloride and n-decanoic acid for adjustable separation of ethylene and ethane. *Sep. Purif. Technol.* **2022**, *298*, 121680. [CrossRef]

73. Cao, Y.; Tao, X.; Jiang, S.; Gao, N.; Sun, Z. Tuning thermodynamic properties of deep eutectic solvents for achieving highly efficient photothermal sensor. *J. Mol. Liq.* **2020**, *308*, 113163. [CrossRef]

74. Lomba, L.; Ribate, M.P.; Zaragoza, E.; Concha, J.; Garralaga, M.P.; Errazquin, D.; García, C.B.; Giner, B. Deep Eutectic Solvents: Are They Safe? *Appl. Sci.* **2021**, *11*, 10061. [CrossRef]

75. De Morais, P.; Gonçalves, F.; Coutinho, J.A.P.; Ventura, S.P.M. Ecotoxicity of Cholinium-Based Deep Eutectic Solvents. *ACS Sustain. Chem. Eng.* **2015**, *3*, 3398–3404. [CrossRef]

76. Macário, I.P.E.; Jesus, F.; Pereira, J.L.; Ventura, S.P.M.; Gonçalves, A.M.M.; Coutinho, J.A.P.; Gonçalves, F.J.M. Unraveling the ecotoxicity of deep eutectic solvents using the mixture toxicity theory. *Chemosphere* **2018**, *212*, 890–897. [CrossRef] [PubMed]

77. Nejrotti, S.; Antenucci, A.; Pontremoli, C.; Gontrani, L.; Barbero, N.; Carbone, M.; Bonomo, M. Critical Assessment of the Sustainability of Deep Eutectic Solvents: A Case Study on Six Choline Chloride-Based Mixtures. *ACS Omega* **2022**, *7*, 47449–47461. [CrossRef] [PubMed]

78. Pedro, S.N.; Freire, C.S.R.; Silvestre, A.J.D.; Freire, M.G. Deep Eutectic Solvents and Pharmaceuticals. *Encyclopedia 2021*, *1*, 942–963. [CrossRef]

79. Cysewski, P.; Jeliński, T.; Przybyłek, M. Experimental and Theoretical Insights into the Intermolecular Interactions in Saturated Systems of Dapsone in Conventional and Deep Eutectic Solvents. *Molecules* **2024**, *29*, 1743. [CrossRef] [PubMed]

80. Mustafa, N.R.; Spelbos, V.S.; Witkamp, G.J.; Verpoorte, R.; Choi, Y.H. Solubility and stability of some pharmaceuticals in natural deep eutectic solvents-based formulations. *Molecules* **2021**, *26*, 2645. [CrossRef]

81. Bazzo, G.C.; Pezzini, B.R.; Stulzer, H.K. Eutectic mixtures as an approach to enhance solubility, dissolution rate and oral bioavailability of poorly water-soluble drugs. *Int. J. Pharm.* **2020**, *588*, 119741. [CrossRef] [PubMed]

82. Kapre, S.; Palakurthi, S.S.; Jain, A.; Palakurthi, S. DES-igning the future of drug delivery: A journey from fundamentals to drug delivery applications. *J. Mol. Liq.* **2024**, *400*, 124517. [CrossRef]

83. Jeliński, T.; Przybyłek, M.; Mianowana, M.; Misiak, K.; Cysewski, P. Deep Eutectic Solvents as Agents for Improving the Solubility of Edaravone: Experimental and Theoretical Considerations. *Molecules* **2024**, *29*, 1261. [CrossRef] [PubMed]

84. Duarte, A.R.C.; Ferreira, A.S.D.; Barreiros, S.; Cabrita, E.; Reis, R.L.; Paiva, A. A comparison between pure active pharmaceutical ingredients and therapeutic deep eutectic solvents: Solubility and permeability studies. *Eur. J. Pharm. Biopharm.* **2017**, *114*, 296–304. [CrossRef]

85. Nguyen, C.-H.; Augis, L.; Fourmentin, S.; Barratt, G.; Legrand, F.-X. Deep Eutectic Solvents for Innovative Pharmaceutical Formulations. In *Deep Eutectic Solvents for Medicine, Gas Solubilization and Extraction of Natural Substances*; Springer: Berlin/Heidelberg, Germany, 2021; pp. 41–102.

86. Liu, Y.; Wu, Y.; Liu, J.; Wang, W.; Yang, Q.; Yang, G. Deep eutectic solvents: Recent advances in fabrication approaches and pharmaceutical applications. *Int. J. Pharm.* **2022**, *622*, 121811. [CrossRef] [PubMed]

87. Emami, S.; Shayanfar, A. Deep eutectic solvents for pharmaceutical formulation and drug delivery applications. *Pharm. Dev. Technol.* **2020**, *25*, 779–796. [CrossRef]

88. Pedro, S.N.; Freire, M.G.; Freire, C.S.R.; Silvestre, A.J.D. Deep eutectic solvents comprising active pharmaceutical ingredients in the development of drug delivery systems. *Expert Opin. Drug Deliv.* **2019**, *16*, 497–506. [CrossRef]
89. Lomba, L.; Garralaga, M.P.; Werner, Á.; Giner, B.; Baptista, P.M.; Sánchez-Romero, N. Ibuprofen solubility and cytotoxic study of deep eutectic solvents formed by xylitol, choline chloride and water. *J. Drug Deliv. Sci. Technol.* **2023**, *82*, 104327. [CrossRef]

90. Guinet, Y.; Paccou, L.; Hédoux, A. Analysis of xylitol—Citric acid system forming deep eutectic solvent: Application for dissolving poorly water-soluble drugs. A combination of calorimetric and Raman investigations. *J. Mol. Liq.* **2020**, *318*, 114317. [CrossRef]

91. Roda, A.; Paiva, A.; Duarte, A.R.C. Therapeutic Liquid Formulations Based on Low Transition Temperature Mixtures for the Incorporation of Anti-Inflammatory Drugs. *Pharmaceutics* **2021**, *13*, 1620. [CrossRef]

92. Hammond, O.S.; Bowron, D.T.; Edler, K.J.; Hammond, S.; Edler, K.J.; Bowron, D.T. The Effect of Water upon Deep Eutectic Solvent Nanostructure: An Unusual Transition from Ionic Mixture to Aqueous Solution. *Angew. Chemie Int. Ed.* **2017**, *56*, 9782–9785. [CrossRef] [PubMed]

93. Gabriele, F.; Chiarini, M.; Germani, R.; Tiecco, M.; Spreti, N. Effect of water addition on choline chloride/glycol deep eutectic solvents: Characterization of their structural and physicochemical properties. *J. Mol. Liq.* **2019**, *291*, 111301. [CrossRef]

94. Gracin, S.; Rasmuson, Å.C. Solubility of Phenylacetic Acid, p-Hydroxyphenylacetic Acid, p-Aminophenylacetic Acid, p-Hydroxybenzoic Acid, and Ibuprofen in Pure Solvents. *J. Chem. Eng. Data* **2002**, *47*, 1379–1383. [CrossRef]

95. Garzón, L.C.; Martínez, F. Temperature dependence of solubility for ibuprofen in some organic and aqueous solvents. *J. Solution Chem.* **2004**, *33*, 1379–1395. [CrossRef]

96. Soto, R.; Verma, V.; Sadeghi, M.; Rasmuson, Å.C. Ketoprofen Solubility in Pure Organic Solvents Using In Situ FTIR and UV–Vis and Analysis of Solution Thermodynamics. *Org. Process Res. Dev.* **2021**, *25*, 2403–2414. [CrossRef]

97. Wang, S.; Song, Z.; Wang, J.; Dong, Y.; Wu, M. Solubilities of Ibuprofen in Different Pure Solvents. *J. Chem. Eng. Data* **2010**, *55*, 5283–5285. [CrossRef]

98. Pacheco, D.P.; Manrique, Y.J.; Martínez, F. Thermodynamic study of the solubility of ibuprofen and naproxen in some ethanol+propylene glycol mixtures. *Fluid Phase Equilib.* **2007**, *262*, 23–31. [CrossRef]

99. Martínez, J.; Manrique, F. Solubility of Ibuprofen in Some Ethanol + Water Cosolvent Mixtures at Several Temperatures. *Lat. Am. J. Pharm.* **2007**, *26*, 344–354.

100. Domańska, U.; Pobudkowska, A.; Pelczarska, A.; Gierycz, P. p K a and Solubility of Drugs in Water, Ethanol, and 1-Octanol. *J. Phys. Chem. B* **2009**, *113*, 8941–8947. [CrossRef]

101. Filippa, M.A.; Gasull, E.I. Ibuprofen solubility in pure organic solvents and aqueous mixtures of cosolvents: Interactions and thermodynamic parameters relating to the solvation process. *Fluid Phase Equilib.* **2013**, *354*, 185–190. [CrossRef]

102. Stovall, D.M.; Givens, C.; Keown, S.; Hoover, K.R.; Rodriguez, E.; Acree, W.E.; Abraham, M.H. Solubility of crystalline nonelectrolyte solutes in organic solvents: Mathematical correlation of ibuprofen solubilities with the Abraham solvation parameter model. *Phys. Chem. Liq.* **2005**, *43*, 261–268. [CrossRef]

103. Aragón, D.M.; Rosas, J.E.; Martínez, F. Thermodynamic study of the solubility of ibuprofen in acetone and dichloromethane. *Brazilian J. Pharm. Sci.* **2010**, *46*, 227–235. [CrossRef]

104. Manrique, Y.J.; Pacheco, D.P.; Martínez, F. Thermodynamics of Mixing and Solvation of Ibuprofen and Naproxen in Propylene Glycol + Water Cosolvent Mixtures. *J. Solution Chem.* **2008**, *37*, 165–181. [CrossRef]

105. Jiménez, D.M.; Muñoz, M.M.; Rodríguez, C.J.; Cárdenas, Z.J.; Martínez, F. Solubility and preferential solvation of some non-steroidal anti-inflammatory drugs in methanol + water mixtures at 298.15 K. *Phys. Chem. Liq.* **2016**, *54*, 686–702. [CrossRef]

106. Wang, B.; Wang, L.; Zhang, J.; Li, W.; Wang, J. Measurement and correlation for solubility of dexibuprofen in different solvents from 263.15 to 293.15K. *Thermochim. Acta* **2012**, *540*, 91–97. [CrossRef]

107. Gantiva, M.; Yurquina, A.; Martínez, F. Solution Thermodynamics of Ketoprofen in Ethanol + Water Cosolvent Mixtures. *J. Chem. Eng. Data* **2010**, *55*, 113–118. [CrossRef]

108. Gantiva, M.; Martínez, F. Thermodynamic analysis of the solubility of ketoprofen in some propylene glycol+water cosolvent mixtures. *Fluid Phase Equilib.* **2010**, *293*, 242–250. [CrossRef]

109. Filippa, M.A.; Melo, G.M.; Gasull, E.I. Ketoprofen solubility in organic solvents and aqueous co-solvent systems: Interactions and thermodynamic parameters of solvation. *J. Pharm. Chem. Biol. Sci.* **2015**, *3*, 440–453.

110. Yalkowsky, S.H.; Valvani, S.C.; Roseman, T.J. Solubility and Partitioning VI: Octanol Solubility and Octanol-Water Partition Coefficients. *J. Pharm. Sci.* **1983**, *72*, 866–870. [CrossRef]

111. Jouyban, A. *Handbook of Solubility Data for Pharmaceuticals*; CRC Press: Boca Raton, FL, USA, 2009; ISBN 9781439804889.

112. Kurkov, S.V.; Perlovich, G.L. Thermodynamic studies of Fenbufen, Diflunisal, and Flurbiprofen: Sublimation, solution and solvation of biphenyl substituted drugs. *Int. J. Pharm.* **2008**, *357*, 100–107. [CrossRef]

113. Przybyłek, M.; Recki, Ł.; Mroczyńska, K.; Jeliński, T.; Cysewski, P. Experimental and theoretical solubility advantage screening of bi-component solid curcumin formulations. *J. Drug Deliv. Sci. Technol.* **2019**, *50*, 125–135. [CrossRef]

114. Jeliński, T.; Przybyłek, M.; Cysewski, P. Solubility advantage of sulfanilamide and sulfacetamide in natural deep eutectic systems: Experimental and theoretical investigations. *Drug Dev. Ind. Pharm.* **2019**, *45*, 1120–1129. [CrossRef] [PubMed]

115. Jeliński, T.; Stasiak, D.; Kosmalski, T.; Cysewski, P. Experimental and theoretical study on theobromine solubility enhancement in binary aqueous solutions and ternary designed solvents. *Pharmaceutics* **2021**, *13*, 1118. [CrossRef]

116. Cysewski, P.; Jeliński, T.; Przybyłek, M. Application of COSMO-RS-DARE as a Tool for Testing Consistency of Solubility Data: Case of Coumarin in Neat Alcohols. *Molecules* **2022**, *27*, 5274. [CrossRef]

117. Jeliński, T.; Przybyłek, M.; Różalski, R.; Cysewski, P. Solubility of dapsone in deep eutectic solvents: Experimental analysis, molecular insights and machine learning predictions. *Polym. Med.* **2024**, *online ahead of print*. [CrossRef]
89. Lomba, L.; Garralaga, M.P.; Werner, Á.; Giner, B.; Baptista, P.M.; Sánchez-Romero, N. Ibuprofen solubility and cytotoxic study of deep eutectic solvents formed by xylitol, choline chloride and water. *J. Drug Deliv. Sci. Technol.* **2023**, *82*, 104327. [CrossRef]

90. Guinet, Y.; Paccou, L.; Hédoux, A. Analysis of xylitol—Citric acid system forming deep eutectic solvent: Application for dissolving poorly water-soluble drugs. A combination of calorimetric and Raman investigations. *J. Mol. Liq.* **2020**, *318*, 114317. [CrossRef]

91. Roda, A.; Paiva, A.; Duarte, A.R.C. Therapeutic Liquid Formulations Based on Low Transition Temperature Mixtures for the Incorporation of Anti-Inflammatory Drugs. *Pharmaceutics* **2021**, *13*, 1620. [CrossRef]

92. Hammond, O.S.; Bowron, D.T.; Edler, K.J.; Hammond, S.; Edler, K.J.; Bowron, D.T. The Effect of Water upon Deep Eutectic Solvent Nanostructure: An Unusual Transition from Ionic Mixture to Aqueous Solution. *Angew. Chemie Int. Ed.* **2017**, *56*, 9782–9785. [CrossRef] [PubMed]

93. Gabriele, F.; Chiarini, M.; Germani, R.; Tiecco, M.; Spreti, N. Effect of water addition on choline chloride/glycol deep eutectic solvents: Characterization of their structural and physicochemical properties. *J. Mol. Liq.* **2019**, *291*, 111301. [CrossRef]

94. Gracin, S.; Rasmuson, Å.C. Solubility of Phenylacetic Acid, p-Hydroxyphenylacetic Acid, p-Aminophenylacetic Acid, p-Hydroxybenzoic Acid, and Ibuprofen in Pure Solvents. *J. Chem. Eng. Data* **2002**, *47*, 1379–1383. [CrossRef]

95. Garzón, L.C.; Martínez, F. Temperature dependence of solubility for ibuprofen in some organic and aqueous solvents. *J. Solution Chem.* **2004**, *33*, 1379–1395. [CrossRef]

96. Soto, R.; Verma, V.; Sadeghi, M.; Rasmuson, Å.C. Ketoprofen Solubility in Pure Organic Solvents Using In Situ FTIR and UV–Vis and Analysis of Solution Thermodynamics. *Org. Process Res. Dev.* **2021**, *25*, 2403–2414. [CrossRef]

97. Wang, S.; Song, Z.; Wang, J.; Dong, Y.; Wu, M. Solubilities of Ibuprofen in Different Pure Solvents. *J. Chem. Eng. Data* **2010**, *55*, 5283–5285. [CrossRef]

98. Pacheco, D.P.; Manrique, Y.J.; Martínez, F. Thermodynamic study of the solubility of ibuprofen and naproxen in some ethanol+propylene glycol mixtures. *Fluid Phase Equilib.* **2007**, *262*, 23–31. [CrossRef]

99. Martínez, J.; Manrique, F. Solubility of Ibuprofen in Some Ethanol + Water Cosolvent Mixtures at Several Temperatures. *Lat. Am. J. Pharm.* **2007**, *26*, 344–354.

100. Domańska, U.; Pobudkowska, A.; Pelczarska, A.; Gierycz, P. p K a and Solubility of Drugs in Water, Ethanol, and 1-Octanol. *J. Phys. Chem. B* **2009**, *113*, 8941–8947. [CrossRef]

101. Filippa, M.A.; Gasull, E.I. Ibuprofen solubility in pure organic solvents and aqueous mixtures of cosolvents: Interactions and thermodynamic parameters relating to the solvation process. *Fluid Phase Equilib.* **2013**, *354*, 185–190. [CrossRef]

102. Stovall, D.M.; Givens, C.; Keown, S.; Hoover, K.R.; Rodriguez, E.; Acree, W.E.; Abraham, M.H. Solubility of crystalline nonelectrolyte solutes in organic solvents: Mathematical correlation of ibuprofen solubilities with the Abraham solvation parameter model. *Phys. Chem. Liq.* **2005**, *43*, 261–268. [CrossRef]

103. Aragón, D.M.; Rosas, J.E.; Martínez, F. Thermodynamic study of the solubility of ibuprofen in acetone and dichloromethane. *Brazilian J. Pharm. Sci.* **2010**, *46*, 227–235. [CrossRef]

104. Manrique, Y.J.; Pacheco, D.P.; Martínez, F. Thermodynamics of Mixing and Solvation of Ibuprofen and Naproxen in Propylene Glycol + Water Cosolvent Mixtures. *J. Solution Chem.* **2008**, *37*, 165–181. [CrossRef]

105. Jiménez, D.M.; Muñoz, M.M.; Rodríguez, C.J.; Cárdenas, Z.J.; Martínez, F. Solubility and preferential solvation of some non-steroidal anti-inflammatory drugs in methanol + water mixtures at 298.15 K. *Phys. Chem. Liq.* **2016**, *54*, 686–702. [CrossRef]

106. Wang, B.; Wang, L.; Zhang, J.; Li, W.; Wang, J. Measurement and correlation for solubility of dexibuprofen in different solvents from 263.15 to 293.15K. *Thermochim. Acta* **2012**, *540*, 91–97. [CrossRef]

107. Gantiva, M.; Yurquina, A.; Martínez, F. Solution Thermodynamics of Ketoprofen in Ethanol + Water Cosolvent Mixtures. *J. Chem. Eng. Data* **2010**, *55*, 113–118. [CrossRef]

108. Gantiva, M.; Martínez, F. Thermodynamic analysis of the solubility of ketoprofen in some propylene glycol+water cosolvent mixtures. *Fluid Phase Equilib.* **2010**, *293*, 242–250. [CrossRef]

109. Filippa, M.A.; Melo, G.M.; Gasull, E.I. Ketoprofen solubility in organic solvents and aqueous co-solvent systems: Interactions and thermodynamic parameters of solvation. *J. Pharm. Chem. Biol. Sci.* **2015**, *3*, 440–453.

110. Yalkowsky, S.H.; Valvani, S.C.; Roseman, T.J. Solubility and Partitioning VI: Octanol Solubility and Octanol-Water Partition Coefficients. *J. Pharm. Sci.* **1983**, *72*, 866–870. [CrossRef]

111. Jouyban, A. *Handbook of Solubility Data for Pharmaceuticals*; CRC Press: Boca Raton, FL, USA, 2009; ISBN 9781439804889.

112. Kurkov, S.V.; Perlovich, G.L. Thermodynamic studies of Fenbufen, Diflunisal, and Flurbiprofen: Sublimation, solution and solvation of biphenyl substituted drugs. *Int. J. Pharm.* **2008**, *357*, 100–107. [CrossRef]

113. Przybyłek, M.; Recki, Ł.; Mroczyńska, K.; Jeliński, T.; Cysewski, P. Experimental and theoretical solubility advantage screening of bi-component solid curcumin formulations. *J. Drug Deliv. Sci. Technol.* **2019**, *50*, 125–135. [CrossRef]

114. Jeliński, T.; Przybyłek, M.; Cysewski, P. Solubility advantage of sulfanilamide and sulfacetamide in natural deep eutectic systems: Experimental and theoretical investigations. *Drug Dev. Ind. Pharm.* **2019**, *45*, 1120–1129. [CrossRef] [PubMed]

115. Jeliński, T.; Stasiak, D.; Kosmalski, T.; Cysewski, P. Experimental and theoretical study on theobromine solubility enhancement in binary aqueous solutions and ternary designed solvents. *Pharmaceutics* **2021**, *13*, 1118. [CrossRef]

116. Cysewski, P.; Jeliński, T.; Przybyłek, M. Application of COSMO-RS-DARE as a Tool for Testing Consistency of Solubility Data: Case of Coumarin in Neat Alcohols. *Molecules* **2022**, *27*, 5274. [CrossRef]

117. Jeliński, T.; Przybyłek, M.; Różalski, R.; Cysewski, P. Solubility of dapsone in deep eutectic solvents: Experimental analysis, molecular insights and machine learning predictions. *Polym. Med.* **2024**, online ahead of print. [CrossRef]
118. Cysewski, P.; Jeliński, T.; Przybyłek, M. Finding the Right Solvent: A Novel Screening Protocol for Identifying Environmentally Friendly and Cost-Effective Options for Benzenesulfonamide. *Molecules* **2023**, *28*, 5008. [CrossRef]

119. Akiba, T.; Sano, S.; Yanase, T.; Ohta, T.; Koyama, M. Optuna: A Next-generation Hyperparameter Optimization Framework. In Proceedings of the 25th ACM SIGKDD International Conference on Knowledge Discovery & Data Mining 2019, Anchorage, AK, USA, 4–8 August 2019; pp. 2623–2631. [CrossRef]

120. Klamt, A.; Schüürmann, G. COSMO: A new approach to dielectric screening in solvents with explicit expressions for the screening energy and its gradient. *J. Chem. Soc. Perkin Trans. 1993*, **2**, 799–805. [CrossRef]

121. Cysewski, P.; Przybyłek, M.; Jeliński, T. Predicting sulfanilamide solubility in the binary mixtures using a reference solvent approach. *Polym. Med.* 2024; *online ahead of print*. [CrossRef]

122. Cysewski, P.; Przybyłek, M.; Jeliński, T. Intermolecular Interactions as a Measure of Dapsone Solubility in Neat Solvents and Binary Solvent Mixtures. *Materials* **2023**, *16*, 6336. [CrossRef]

123. Przybyłek, M.; Jeliński, T.; Mianowana, M.; Misiak, K.; Cysewski, P. Exploring the Solubility Limits of Edaravone in Neat Solvents and Binary Mixtures: Experimental and Machine Learning Study. *Molecules* **2023**, *28*, 6877. [CrossRef]

124. Cysewski, P.; Jeliński, T.; Przybyłek, M. Intermolecular Interactions of Edaravone in Aqueous Solutions of Ethaline and Glyceline Inferred from Experiments and Quantum Chemistry Computations. *Molecules* **2023**, *28*, 629. [CrossRef]

125. Acree, W.; Chickos, J.S. Phase Transition Enthalpy Measurements of Organic and Organometallic Compounds and Ionic Liquids. Sublimation, Vaporization, and Fusion Enthalpies from 1880 to 2015. Part 2. C11–C192. *J. Phys. Chem. Ref. Data* **2017**, *46*, 013104. [CrossRef]

126. Dassault Systèmes. COSMOtherm; Version 24.0.0; Biovia: San Diego, CA, USA, 2022.

127. Dassault Systèmes. COSMOconf; Version 24.0.0; Dassault Systèmes, Biovia: San Diego, CA, USA, 2022.

**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.
118. Cysewski, P.; Jeliński, T.; Przybyłek, M. Finding the Right Solvent: A Novel Screening Protocol for Identifying Environmentally Friendly and Cost-Effective Options for Benzenesulfonamide. *Molecules* **2023**, *28*, 5008. [CrossRef]

119. Akiba, T.; Sano, S.; Yanase, T.; Ohta, T.; Koyama, M. Optuna: A Next-generation Hyperparameter Optimization Framework. In Proceedings of the 25th ACM SIGKDD International Conference on Knowledge Discovery & Data Mining 2019, Anchorage, AK, USA, 4–8 August 2019; pp. 2623–2631. [CrossRef]

120. Klamt, A.; Schüürmann, G. COSMO: A new approach to dielectric screening in solvents with explicit expressions for the screening energy and its gradient. *J. Chem. Soc. Perkin Trans. 1993*, **2**, 799–805. [CrossRef]

121. Cysewski, P.; Przybyłek, M.; Jeliński, T. Predicting sulfanilamide solubility in the binary mixtures using a reference solvent approach. *Polym. Med.* 2024; *online ahead of print*. [CrossRef]

122. Cysewski, P.; Przybyłek, M.; Jeliński, T. Intermolecular Interactions as a Measure of Dapsone Solubility in Neat Solvents and Binary Solvent Mixtures. *Materials* **2023**, *16*, 6336. [CrossRef]

123. Przybyłek, M.; Jeliński, T.; Mianowana, M.; Misiak, K.; Cysewski, P. Exploring the Solubility Limits of Edaravone in Neat Solvents and Binary Mixtures: Experimental and Machine Learning Study. *Molecules* **2023**, *28*, 6877. [CrossRef]

124. Cysewski, P.; Jeliński, T.; Przybyłek, M. Intermolecular Interactions of Edaravone in Aqueous Solutions of Ethaline and Glyceline Inferred from Experiments and Quantum Chemistry Computations. *Molecules* **2023**, *28*, 629. [CrossRef]

125. Acree, W.; Chickos, J.S. Phase Transition Enthalpy Measurements of Organic and Organometallic Compounds and Ionic Liquids. Sublimation, Vaporization, and Fusion Enthalpies from 1880 to 2015. Part 2. C11–C192. *J. Phys. Chem. Ref. Data* **2017**, *46*, 013104. [CrossRef]

126. Dassault Systèmes. COSMOtherm; Version 24.0.0; Biovia: San Diego, CA, USA, 2022.

127. Dassault Systèmes. COSMOconf; Version 24.0.0; Dassault Systèmes, Biovia: San Diego, CA, USA, 2022.

**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.