# Intermolecular Interactions as a Measure of Dapsone Solubility in Neat Solvents and Binary Solvent Mixtures

Piotr Cysewski, Maciej Przybylek, Tomasz Jeliński

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

## To cite this version:

Piotr Cysewski, Maciej Przybylek, Tomasz Jeliński. Intermolecular Interactions as a Measure of Dapsone Solubility in Neat Solvents and Binary Solvent Mixtures. Materials, 2023, 16 (18), pp.6336. 10.3390/ma16186336. [hal-04215566](https://hal.science/hal-04215566v1)

HAL Id: hal-04215566

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

Submitted on 22 Sep 2023

HAL is a multi-disciplinary open access archive for the deposit and dissemination of scientific research documents, whether they are published 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>

Distributed under a Creative Commons CC BY 4.0 - Attribution - International License
# Intermolecular Interactions as a Measure of Dapsone Solubility in Neat Solvents and Binary Solvent Mixtures

Piotr Cysewski * and Maciej Przybyłe k and Tomasz Jeliński

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

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

**Abstract:** Dapsone is an effective antibacterial drug used to treat a variety of conditions. However, the aqueous solubility of this drug is limited, as is its permeability. This study expands the available solubility data pool for dapsone by measuring its solubility in several pure organic solvents: N-methyl-2-pyrrolidone (CAS: 872-50-4), dimethyl sulfoxide (CAS: 67-68-5), 4-formylmorpholine (CAS: 4394-85-8), tetraethylene pentamine (CAS: 112-57-2), and diethylene glycol bis(3-aminopropyl) ether (CAS: 4246-51-9). Furthermore, the study proposes the use of intermolecular interactions as molecular descriptors to predict the solubility of dapsone in neat solvents and binary mixtures using machine learning models. An ensemble of regressors was used, including support vector machines, random forests, gradient boosting, and neural networks. Affinities of dapsone to solvent molecules were calculated using COSMO-RS and used as input for model training. Due to the polymorphic nature of dapsone, fusion data are not available, which prohibits the direct use of COSMO-RS for solubility calculations. Therefore, a consonance solvent approach was tested, which allows an indirect estimation of the fusion properties. Unfortunately, the resulting accuracy is unsatisfactory. In contrast, the developed regressors showed high predictive potential. This work documents that intermolecular interactions characterized by solute-solvent contacts can be considered valuable molecular descriptors for solubility modeling and that the wealth of encoded information is sufficient for solubility predictions for new systems, including those for which experimental measurements of thermodynamic properties are unavailable.

**Keywords:** dapsone; solubility; machine learning; intermolecular interactions; affinity; COSMO-RS; neat solvents; binary mixtures

## 1. Introduction

Dapsone (DAP, CAS number: 80-08-0, 4,4'-diaminodiphenylsulfone, C₁₂H₁₂N₂O₂S) is an effective antibacterial drug, being a part of the sulphone class, used to treat a variety of conditions [1]. It is used in both topical and systemic forms to treat a variety of conditions, including leprosy, malaria, dermatitis herpetiformis, acne, and diseases associated with AIDS [2–5]. Dapsone works by inhibiting the production of folic acid through the competitive inhibition of dihydropteroate synthetase, which ultimately disrupts the synthesis of nucleic acids essential for bacterial survival and reproduction, accounting for DAP's antibacterial activity [2,6]. Furthermore, the anti-inflammatory effect of dapsone can also be attributed to its ability to regulate the production of cytokines [7,8]. Dapsone also has antioxidant properties by limiting the generation of ROS and superoxide radicals. This is achieved by binding to NADPH oxidase [9]. The liver metabolizes dapsone through acetylation and hydroxylation, with primary elimination occurring via urine [10,11]. There are some negative effects associated with dapsone, which include, among others, some hematological problems and peripheral neuropathy, with the potential for hypersensitivity syndrome and a life-threatening drug reaction [12]. Due to its low water solubility and

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

check for
updates

**Citation:** Cysewski, P.; Przybyłe k, M.; Jeliński, T. Intermolecular Interactions as a Measure of Dapsone Solubility in Neat Solvents and Binary Solvent Mixtures. *Materials* **2023**, *16*, 6336. https://doi.org/10.3390/ma16186336

Academic Editor: Annalisa Paolone

Received: 7 September 2023

Revised: 19 September 2023

Accepted: 20 September 2023

Published: 21 September 2023

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

**Copyright:** © 2023 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/).
permeability, dapsone falls under class II of the Biopharmaceutics Classification System (BCS) [13]. Consequently, transdermal administration is favored over oral ingestion [14]. These limitations prompted various approaches to address dapsone's poor solubility and limited bioavailability [15,16]. Notably, an interesting way to enhance the dissolution behavior of dapsone is by using deep eutectic solvents (DESs) [17].

Solubility stands as a fundamental physicochemical property with significant im-
plications in both theoretical and practical aspects [18,19]. Beyond its clear influence on
bioavailability, solubility is an important technological factor [20]. This is understand-
able since, as estimated, most of the chemical compounds used in drug manufacturing
are solvents [21,22]. For this reason, theoretical solubility prediction methods are useful
tools that support the selection of an optimal solvent or solvent mixtures for technologi-
cal purposes, including in the pharmaceutical industry [23–25]. First principle methods
utilizing quantum chemistry computations augmented with statistical analysis, such as
conductor-like screening model for realistic solvents (COSMO-RS) [26–28], deserve special
attention. In general, this method has been found to be successful in modeling various phar-
maceutically relevant properties and characteristics, including solubility [29–32], partition
coefficients [33–35], acid–base properties [36–38], and co-crystallization abilities [27,39–43].
In recent years, the application of neural networks, including deep learning, to solubil-
ity modeling has increased [44–46]. As it was established in the previous works [47–50],
combining COSMO-RS with machine learning methods can achieve high prediction ac-
curacy. The purpose of this paper is threefold. Firstly, the available pool of dapsone
solubility–temperature profiles is carefully cured and augmented with new measurements
in neat solvents, such as N-methyl-2-pyrrolidone, dimethyl sulfoxide, 4-formylmorpholine,
tetraethylene-pentamine, and diethylene glycol bis(3-aminopropyl) ether, for extending the
diversity of the dissolution media. Secondly, the most probable solute–solvent contacts are
characterized on an advanced ab initio level. Finally, the methods of solubility screening
are provided, addressing the fundamental problem of the lack of the fusion thermodynamic
properties of dapsone.

## 2. Materials and Methods

2.1. *Materials*

Dapsone (DAP, CAS Number: 80-08-0) was obtained from Sigma Aldrich (Saint Louis, MO, USA) at a purity of over 99%. The solvents used in the study, namely N-methyl-2-pyrrolidone (NMP, CAS Number: 872-50-4), dimethyl sulfoxide (DMSO, CAS Number: 67-68-5), 4-formylmorpholine (4FM, CAS Number: 4394-85-8), tetraethylenepentamine (TEPA, CAS Number: 112-57-2), and diethylene glycol bis(3-aminopropyl) ether (B3APE, CAS Number: 4246-51-9). Moreover, methanol (CAS Number: 67-56-1) was used as an auxiliary solvent for dilution purposes. All solvents were similarly purchased from Sigma Aldrich, and their purity was no less than 99%.

## 2.2. *Solubility Measurements*

In this study, the shake-flak method of solubility determination was applied. The method has been used and validated in previous work by our group on various pharmaceuticals or drug-like compounds, including sulfa drugs (sulfanilamide [51], sulfamethizole [49]), amides (nicotinamide [52], benzamide, salicylamide, and ethenzamide [32]), acetanilide derivatives (phenacetin [53]), and nutraceuticals (coumarin [54]). Based on this protocol, the solubility determination of dapsone was preceded by the construction of a calibration curve. In the first step, a stock solution of DAP was prepared by dissolving 2.7 mg of dapsone in methanol in a 100 mL volumetric flask. Afterward, several dilutions were made in order to obtain solutions with decreasing solute concentrations in the range from 0.00135 mg/mL to 0.0135 mg/mL with a total of 14 points characterized by varying concentrations. The solutions prepared in this way were subjected to spectrophotometric measurements using an A360 spectrophotometer from AOE Instruments (Shanghai, China). The maximum absorbance of the samples corresponded to the wavelength of 295 nm. Three
separate curves were prepared, and the final one was the result of their averaging. The range of mean absorbance values was from 0.159 to 1.540. The linear regression equation of the final curve was $A = 113.8664 \cdot C + 0.0015$ (A—absorbance, C—concentration expressed in mg/mL), and the determination coefficient indicated a high degree of linearity with $R^2 = 0.999$. The calculated limits of detection and quantification were found to be $LOD = 3.71 \cdot 10^{-4}$ mg/mL and $LOQ = 1.11 \cdot 10^{-3}$ mg/mL, respectively.

To evaluate the solubility of dapsone in different solvents, its excess amounts were introduced into test tubes containing a specific solvent. The resulting saturated solutions were then placed in an Orbital Shaker Incubator ES-20/60 from Biosan (Riga, Latvia) and subjected to incubation at different temperatures over a 24 h period. Four temperature increments were used for this incubation, ranging from 298.15 K to 313.15 K at 5 K intervals. The temperature of the incubator remained within 0.1 degrees, and fluctuations within 0.5 degrees were observed throughout the 24 h cycle. Simultaneously, the samples were agitated at a rate of 60 revolutions per minute. Subsequently, the samples were filtrated using syringes equipped with PTFE filters featuring a pore size of 0.22 µm. To prevent any precipitation arising from temperature disparities between the solutions and the apparatus, all elements, including test tubes, pipette tips, syringes, and filters, underwent preheating. These items were placed within the same incubator as the samples, attaining the same temperature prior to manipulation. This step was particularly vital when dealing with higher temperatures. Following filtration, minor aliquots of the obtained filtrate were diluted within test tubes containing methanol and subjected to spectrophotometric measurement. The density of each solution was determined by weighing a 1 mL volume within 10 mL volumetric flasks, utilizing an Eppendorf Reference 2 pipette (Hamburg, Germany) with a systematic error of 6 µL. Additionally, the RADWAG AS 110 R2.PLUS analytical balance (Radom, Poland) with a precision of 0.1 mg was employed for this purpose. The process of solubility determination utilized the same spectrophotometer as was the case for the calibration curve preparation. Spectral data were captured within the wavelength span of 190 nm to 400 nm, maintaining a resolution of 1 nm. Throughout these procedures, methanol served both as the diluent for the samples and for the initial calibration of the spectrophotometer. The analytical wavelength was specifically set at 295 nm, and the absorbance at this wavelength was used to quantify the DAP concentration present in the samples, subsequently enabling the calculation of its mole fractions. Three distinct measurements were undertaken, with the resulting values being averaged for increased reliability.

### 2.3. Solubility Dataset Curation

Dapsone solubility values have been reported for twelve pure solvents [55] and four binary mixtures [56] as inferred from an extensive literature search. Unfortunately, there are two problems with this collection. First, there are some inconsistencies in the solubility values, which require careful consideration prior to nonlinear model formulation. Second, the solvent space is rather limited with insufficient diversity of solvent properties. In fact, DAP solubility has been collected mainly in polar protic solvents such as alcohols and polar aprotic solvents such as esters or acetone. The first aspect will be addressed by data curation and the latter by expanding the pool of solubility data with new measurements.

Although the diversities of reported DAP solubility in different sources are quite modest, they are still non-trivial and require standardization. This was carried out here by using the three-parameter van't Hoff equation, $\ln(x^{vHF}) = A + B \cdot T^{-1} + C \cdot T^{-2}$ [57–62]. This simple polynomial extension of the original model accounts for a non-constant enthalpy value [62] of the solubilization process. This simple model does not require any additional information other than solubility data and is very accurate in back-calculations, provided that the A, B, and C parameters are optimized. This was performed by using gradient optimization to minimize the values of RMSD (root means square deviation) using the solver implemented in MS Excel. The back-computed solubility values for the experimental temperatures were used for model building. The characteristics of the whole curated
solubility dataset are provided in the supporting materials (see Section S1 in File S1) together with the obtained consensus values used for model formulation (see column "log(x<sup>CONS</sup>)" in the "data" spreadsheet of the File S2). The values of solubility data in binary mixtures were used as reported in the literature [56] except for neat solvents, which were replaced with consensus data.

## 2.4. *Computations Protocol*

### 2.4.1. COSMO-RS Solubility Computations

The COSMO-RS approach [26–28] implemented in BIOVIA COSMOtherm 2021 (build: 21.0.0 [d1b290c105]) [63] was used for dapsone solubility computations. Since the iterative procedure occasionally fails in solubility computations, especially for cases with high solubility, the complete solution of solid–liquid equilibrium (SLE) was used by toggling the SLESOL option. Furthermore, since COSMOtherm can only treat liquids, a hypothetical subcooled liquid state is postulated in the case of solid solubility, and the thermodynamic contribution of an ordered solid transition to the random particle distribution in the subcooled liquid requires the provision of the values of the Gibbs free energy of fusion. Since neither the melting temperature, nor the heat of fusion, nor the change in heat capacity upon melting is experimentally undeterminable, the reference solvent approach was used. The main advantage of this practical approach is the ability to evaluate fusion thermodynamics from the provided experimental solubility in another solvent. Unfortunately, the calculated solubility values are strongly influenced by the choice of reference data. This is due to the fact that the errors of COMSO-RS are of similar magnitude for compounds with comparable structures. Therefore, screening is essential to optimize the number and type of reference solvents. This set of best-selected reference solvents is called conso-nance solvents and can minimize the overall error of solubility calculations [64]. For an adequate representation of the structure of dapsone and all solvents, the sets of relevant conformations were generated using COSMOconf 4.2 [65] and optimized with the aid of TURBOMOLE version TURBOMOLE V7-5-1 (V7-5-1 23 Dez 2020, Dassault Systèmes: Vélizy-Villacoublay, France,) [66] interfaced with BIOVIA TmoleX version 21.0.1 [67] as a default engine for geometries optimization. The obtained conformers used further for characteristics of bulk systems had their geometries fully optimized using BP functional and TZVP basis sets. All structures were generated both in the gas phase and including environmental effects via the COSMO-RS solvation model [28]. For solubility computations, the TZVPD-FINE level was used, which corresponds to single-point calculations with the TZVPD basis set and the same density functional based on previously generated geometries. The BP_TZVPD_FINE_21.ctd parameter set was used for all physicochemical property computations using COSMOtherm.

### 2.4.2. Affinity Characteristic of Solute–Solvent Systems

The conformational screening of solute–solute and solute–solvent bi-molecular systems was performed prior to affinity characteristics. The methodology is consistent with previously published work [50,68]; hence, only a brief description is given below. This step was aimed at finding geometries of the most probable clusters. For this purpose, the COSMOtherm program facilities were used by calculating the contact statistics based on the probability of interactions between molecule surface segments. Practically, it is performed by using the "CONTACT={1 2} sscProbability ssc_weak ssc Ang=15.0" command and automatic generation of contacts by alteration of the mutual orientation of the two contacting molecules with a 15° step rotation interval. Weak interactions are also included in the probability statistics as evidenced by the above prompt. Usually, this leads to a quite large number of potential structures whose geometries are far from optimal. Hence, the structure optimizations were performed using RI-DFT BP86 (B88-VWN-P86). In the final step, the number of pairs was reduced by comparing their energy values and RMSD values after cluster overlapping. Highly similar clusters and the ones exceeding the 2.5 kcal/mol threshold window of relative energy were discarded. The selected pairs of conformers
were prone to single-point energy computations using the def2-TZVPD basis set with the fine grid tetrahedron cavity and the inclusion of parameter sets with hydrogen bond interaction and the van der Waals dispersion term based on the “D3” method of Grimme et al. [69]. These final computations are performed to preserve consistency with monomer characteristics. The values of Gibbs free energies were determined using COSMOtherm with BP_TZVPD Fine_21.ctd parametrization. Cluster energies were additionally characterized by the inclusion of corrections accounting for zero-point vibrational energy (RI-DFT BP86 (B88-VWN-P86)/def2-TZVPD level) and electron correlation (RI-MP2/def2-QZVPP level). Moreover, the BSSE was estimated using the DFT-C approach [70], the formulation of which includes atom-atom many-body corrections and which is a parameterized geometry-based method. All energy corrections were calculated for each conformer of each cluster, averaged with weights corresponding to the population fraction estimated using Boltzmann probability. Hence, $\Delta E_g^{cor}$, $\Delta E_g^{ZPE}$, and $\Delta E_g^{BSSE}$ characterize the correction for a given pair including contributions coming from all conformations of a given cluster.

The solute-solvent affinity was represented by the values of the Gibbs free energy corresponding reactions, $\Delta G_r$, of pair formation X + Y = XY, where X and Y represent either solute or solvent molecules. In the case of X = Y, dapsone dimers are formed, and in the case of X ≠ Y, solute-solvent heteromolecular binary contacts are considered. It is worth mentioning that computations of $\Delta G_r$ might result either in a concentration-independent measure defined based on activity values, $\Delta G_r(a)$, or concentration-dependent items defined using a mole fraction, $\Delta G_r(x)$. The two are interrelated via the activity coefficient product. For the purpose of machine learning, the latter was used, but for overall affinity characteristics, the former is most suitable. In addition, the values of the Gibbs free energies of a solution of dapsone, dapsone dimers, and dapsone-solvent pairs in a given solvent were extracted from the COSMOtherm output files. These values were used as additional molecular descriptors besides the values of Gibbs free energies of cluster formation. The enthalpic and entropic contributions to the affinities were also included in the pool of molecular descriptors. The whole descriptors dataset is provided in supporting materials (see columns L-T in the “data” spreadsheet of the File S2).

### 2.4.3. Machine Learning Protocol

The solubility prediction model was formulated using the in-house Python (ver.3.10, https://www.python.org/) code developed for hyperparameter tuning of 36 regression models. They use a variety of algorithms including linear models, boosting, ensembles, nearest neighbors, neural networks, and also some other types of regressors. The hyperparameter space was explored to find their optimal values using the Optuna study (ver.3.2, https://optuna.org/), a freely available Python package for hyperparameter optimization [71]. Model tuning was performed through 5000 minimization trials using the tree-structured Parzen estimator (TPE) as the search algorithm sampler. This computationally efficient model-based optimization algorithm uses a probability density function to model the relationship between hyperparameters and performance metrics. To evaluate the performance of each regression model, a new custom score function was developed that combines multiple metrics to account for both the accuracy and generalizability of the model as defined in the previous work [47], where the mathematical details are provided. The most important aspect of the applied scoring function is the inclusion of the penalties obtained from the learning curve analysis (LCA) of the scikit-learn 1.2.2 library during parameter tuning. Since LCA can be computationally expensive, only two-point computations were performed here by including 50% and 100% of the total data. The LCA evaluations of the final model were performed using 20-point calculations in the range of 50%–100%. The values included in the custom loss correspond to the average MAE values obtained at the largest training set size. Thus, such a custom loss function combines the two types of components and provides information about the model’s accuracy and ability to generalize to new, unseen data. The final performance of all modes was evaluated using the loss values characterizing the test and validation subsets. The ensemble model (EM)
was defined by including the subset of regression modes with the lowest values of both criteria, and the final predictions were averaged over selected models.

### 3. Results and Discussion

#### 3.1. New Data of Dapsone Solubility

In order to extend the available solubility data for dapsone, a series of measurements were performed with five different organic solvents, namely N-methyl-2-pyrrolidone (NMP), dimethylsulfoxide (DMSO), 4-formylmorpholine (4FM), tetraethylene pentamine (TEPA), and diethylene glycol bis(3-aminopropyl) ether (B3APE), at a temperature range of 298.15 K to 313.15 K. The results of these measurements are summarized in Table 1.

<table><caption>**Table 1.** The solubility of dapsone in neat organic solvents at various temperatures expressed as molar concentration s<sub>DAP</sub> and mole fraction (x<sub>DAP</sub>, 10<sup>4</sup>). Standard deviation values are given in parentheses.</caption><thead><tr><th rowspan="2">T (K)</th><th colspan="5">s<sub>DAP</sub> (mol/dm<sup>3</sup>)</th><th colspan="4">x<sub>DAP</sub>·10<sup>4</sup></th></tr><tr><th>4FM</th><th>DMSO</th><th>TEPA</th><th>NMP</th><th>B3APE</th><th>4FM</th><th>DMSO</th><th>TEPA</th><th>NMP</th><th>B3APE</th></tr></thead><tbody><tr><td>298.15</td><td>0.042<br>(±0.000)</td><td>0.256<br>(±0.009)</td><td>0.056<br>(±0.001)</td><td>0.108<br>(±0.003)</td><td>0.068<br>(±0.001)</td><td>43.11<br>(±0.25)</td><td>187.57<br>(±7.02)</td><td>105.69<br>(±2.23)</td><td>105.13<br>(±3.20)</td><td>147.51<br>(±1.43)</td></tr><tr><td>303.15</td><td>0.074<br>(±0.004)</td><td>0.344<br>(±0.002)</td><td>0.079<br>(±0.002)</td><td>0.243<br>(±0.005)</td><td>0.109<br>(±0.011)</td><td>73.90<br>(±4.18)</td><td>251.58<br>(±1.54)</td><td>149.90<br>(±4.00)</td><td>236.59<br>(±5.16)</td><td>231.14<br>(±22.69)</td></tr><tr><td>308.15</td><td>0.126<br>(±0.001)</td><td>0.462<br>(±0.012)</td><td>0.104<br>(±0.000)</td><td>0.525<br>(±0.009)</td><td>0.142<br>(±0.002)</td><td>125.57<br>(±0.89)</td><td>342.88<br>(±9.24)</td><td>196.54<br>(±0.47)</td><td>524.74<br>(±10.19)</td><td>297.60<br>(±4.48)</td></tr><tr><td>313.15</td><td>0.191<br>(±0.001)</td><td>0.598<br>(±0.011)</td><td>0.141<br>(±0.002)</td><td>1.040<br>(±0.118)</td><td>0.191<br>(±0.003)</td><td>191.48<br>(±0.55)</td><td>444.27<br>(±9.20)</td><td>262.29<br>(±4.27)</td><td>1113.63<br>(±17.53)</td><td>397.51<br>(±6.06)</td></tr></tbody></table>

Some interesting observations can be made when analyzing the presented results. At a temperature of 298.15 K, the highest solubility of dapsone was found in the case of DMSO with x<sub>DAP</sub> = 187.57·10<sup>-4</sup>, and a decreasing solubility trend of DMSO > B3APE > TEPA > NMP > 4FM is present. With increasing temperature, the DAP solubility also increases, although not in an even manner for all studied solvents. The highest solubility increase from 298.15 K to 313.15 K was observed for NMP (of 10.6 times in mole fractions) with the lowest in the case of DMSO (of 2.4 times in mole fractions). This resulted in a change in the overall solubility trend with NMP yielding the highest solubility of dapsone at 313.15 K, amounting to x<sub>DAP</sub> = 1113.63·10<sup>-4</sup>.

#### 3.2. Solubility Prediction Using Consonance Solvents

The COSMO-RS framework is often used for theoretical solubility evaluations. However, since the approach is defined for fluids, its application to saturated solid-liquid systems is only possible if information on the melting properties is available. Unfortunately, in many cases, these data are not available due to a lack of reported measurements or, even worse, due to the impossibility of making such measurements. Many solids cannot be studied by using the classical DSC approach due to thermal instabilities or complex phase transitions [72–74]. This is exactly the case for dapsone [55,75]. Solid DAP exhibits quite complex equilibria as it can exist in five anhydrous forms [76]. It happens so that polymorph III is predominantly present in commercial products as the most stable in ambient conditions. Unfortunately, this polymorph undergoes a transition to form II below the melting point [55,75]. Therefore, it is not possible to accurately measure the melting temperature of this form, which exists in the equilibrated saturated solutions under conditions of solubility determination. Consequently, the experimental determination of the enthalpy of fusion and the heat capacity is also prohibited due to this circumstance. This inherent limitation can be overcome by estimating T<sub>m</sub> and H<sub>fus</sub> using, for example, the group additive approach, as documented by Li et al. [55]. However, this introduces errors in the solubility computations, the significance of which is difficult to estimate. Therefore,
for the purpose of validating the ability of COSMO-RS to predict the solubility of DAP, a different strategy was adopted using a method called the consonance solvent approach [64]. This approach is based on the identification of reference solvents for which the solubility is known and using these data to indirectly calculate the thermodynamic properties of the melt. Therefore, for the purpose of this study, all possible combinations of reference solvents were considered to estimate the solubility of dapsone. The results are presented in the form of a heat map, as visualized in Figure 1. The heat map collects the values of the MAPE determined for all possible combinations of solvent–reference solvents. Interestingly, this allows for the grouping of solvents with a similar accuracy of solubility calculations and then for finding the best set of matching solvents. According to chemical intuition, COSMO-RS-derived solubilities are expected to be of similar accuracy for closely related chemical structures of solvents. Indeed, all proton-donating solvents were included in the first subset. Therefore, using alcohol as a reference solvent provided the best estimate of solubility for alcohols. In this case, isobutanol (iBuOH) seems to be the best choice. For DAP solubility estimation in proton-accepting solvents, it is logically better not to use alcohol but rather the solvent belonging to the same subgroup. Here, ethylene propionate (EtPr) seems to be the best choice for this class of solvents. The last subset identified in Figure 1 includes other polar non-protic solvents, and TEPA was identified as the best choice. Thus, the procedure led to the determination of the minimal set of optimal selection of three reference solvents for minimizing the error of the COSMO-RS derived solubility. A closer inspection of the lower left rectangle of the heat map clearly shows a similarity between NMP and 4FM. The remaining three solvents classified in this region, namely DMSO, TEPA, and B3APE, should represent a distinct subset since the values of percentage error have the opposite sign. To keep the number of congruent solvents to a minimum, only three classes are accepted. However, acetone is in a class of its own and can hardly be considered a reasonable reference solvent, and vice versa—no other solvent can serve as a reference for acetone. Therefore, acetone was excluded from the pool of reference solvents, or to put it another way, it can be expected that the class of other ketones might emerge from this analysis. In the end, four reference solvents were defined that resembled the structural diversity of the solvents in the dapsone solubility dataset. However, despite the efforts made to optimize the solubility calculations, the accuracy of the values obtained is still hardly acceptable, which can be directly deduced from the MAPE values provided. Considering that this statistical measure is calculated on the basis of the decadal logarithm of the mole fraction, the overall accuracy of solubility calculations using the COSMO-RS approach is at best qualitatively related to the measured values. This was the direct impulse for the development of nonlinear models and the formulation of a quantitative tool. To further illustrate this need, the values computed using the optimized set of consonant solvents are confronted with the predictions made by the ensemble of regressor models in the next section.

### 3.3. Solute–Solvent Intermolecular Interactions

The selection of the molecular descriptors is a crucial step in the machine learning protocol. There are many possible parameters, which are used by different authors [77–81], for machine learning. However, the actual selection is limited by the fact that the temperature dependence has to be explicitly included for an adequate description of the pool of experimental dapsone solubility. Fortunately, there are many valuable molecular descriptors that could be derived using the COSMO-RS approach, such as calculated values of chemical potential, activities, or different energetic contributions. It has already been documented that they can serve as quite reasonable quantities for training purposes [47]. Nevertheless, an alternative set of molecular properties was tested here. Solute–solute and solute–solvent affinities were determined as described in the ‘Materials and Methods’ section. Thus, the most likely homo- and heteromolecular pairs were screened, for which the values of the Gibbs free energies of the synthesis reaction from monomers were determined. The first step was to analyze the self-association of two DAP molecules. Although dapsone can be
considered a molecule with potential proton acceptor and proton donor centers, the self-association of two DAP molecules adopts a stacking conformation in all solvents analyzed rather than stabilization by hydrogen bonding. Conformational analysis identified two types of stacking complexes in which either a parallel or antiparallel orientation of two DAP molecules occurs, as documented in Figure 2.

<figure><img src="image_5.png" /><figcaption>**Figure 1.** The collection of MAPE (mean absolute percentage error) values for all combinations of solvent-reference solvent solubility calculations. The left column denotes the solvent in which the solubility is calculated, and the first row enumerates the reference solvents. Thus, the diagonal represents solubility calculations in a given solvent with itself as a reference, and, therefore, the MAPE is equal to zero. In red color, the names of the consonant solvents are marked, which minimize the MAPE for the subset of solvents enclosed by a rectangle.</figcaption></figure>

<figure><img src="image_6.png" /><figcaption>**Figure 2.** The self-association properties of DAP. The electron charge density profiles of the two most probable conformations in all studied solutions were supplemented with energetic contributions of self-association.</figcaption></figure>
Despite the fact that in the antiparallel orientation, the hydrogen bond is formed between the sulfanyl group of one molecule and the amino side group of the other, this structure is less stable by about 1.1 kcal/mol. This suggests a marginal contribution to the overall thermodynamic properties of DAP. The necessary distortion of the amino group is responsible for such a destabilizing effect. Moreover, the shape of the HB is not expected to be strong, as suggested by the low value of the O···H–N angle.

It is worth mentioning that the stacking conformation is stable in all saturated systems considered and is the most favorable intermolecular complex among all solute–solvent contacts considered here. A fairly linear trend between solubility and DAP stacking affinity is observed, suggesting that this could be a promising descriptor for solubility prediction via machine learning. It is also worth mentioning that the dispersion forces are very strong in the case of DAP stacking, as documented by the values of the electron correlation contribution to the total energy. The contributions from ZPE and BSSE are of opposite signs but only slightly reduce the DAP self-affinity. It is also worth noting that there is a modest correlation of the $\Delta G_r$ of DAP self-association estimated in different solvents with solubility ($R^2 = 0.7$), which is a good prognosis for using these values as molecular descriptors in the nonlinear model training. As expected, the self-affinity of DAP is highest in the case of water, which has the lowest solubility. In contrast, the highest solubility observed in DMSO is associated with the weakest self-affinity of dapsone. Since the stacking of dapsone is similar to the interactions in the crystal, it is reasonable to expect that the promotion of self-association will lead to clustering, which will make it difficult to disperse in the bulk solvents, eventually leading to precipitation. Therefore, the coincidence of solute and solvent in this hydrophobic region seems to be the most important solubility factor. To further explore this suggestion, the heteromolecular complexes potentially present in the solvents analyzed were studied.

Dapsone has weak hydrogen-donating and -accepting properties, allowing interactions with a variety of solvent molecules. The former interactions are provided by the peripheral amine groups, while the latter are provided by the two oxygen atoms on the sulfonyl group. In addition, as mentioned in the context of self-association, DAP has wide delocalized electron clouds that allow for non-specific interactions. This is granted by the pocket-like region formed by two phenyl groups, which form a nonpolar interior attractive for all molecular fractions of a nonpolar nature. Therefore, it is expected that dapsone will be stabilized by different intermolecular interactions with solvents depending on their nature. In fact, upon closer inspection of the structure of the most probable pairs, it is possible to distinguish three classes of systems. In Figure 3, the examples of pairs are presented in which DAP acts as a proton acceptor via the sulfanyl group. These affinities vary from −7.7 kcal/mol in the case of water up to −13.5 kcal/mol for nPeOH (n-pentanol). In all cases, hydrogen bonds are formed that could be classified as strong based on their geometric parameters. The presented structures show that the region of the DAP molecule active in self-association is not affected by all these solute–solvent interactions, suggesting that higher-order trimers or clusters are quite likely. Unfortunately, the size of these contacts prohibits the advanced calculations used in this project. However, this speculation might suggest that the solubility in these proton-donating solvents is expected to be at most modest due to possible dapsone self-association driven by non-specific interactions.
<figure><img src="image_7.png" /><figcaption>**Figure 3.** The characteristics of solute–solvent contacts of DAP in proton-donating solvents like water and alcohols. Prime denotes atoms on HBD counter partner in hydrogen bonding.</figcaption></figure>
In contrast, the second class of solvents collected in Figure 4 includes contacts in which DAP acts as a proton donor. This class includes acetone, DGE (diethylene glycol bis(3-aminopropyl) ether), and TETA (tetraethylene pentamine). It is worth noting that the molecules of the latter two solvents, due to their chain size, are also attracted to the apolar pocket of dapsone. Therefore, it is expected that there is a coincidence between self-association and solvents, which could be used as a rationale for the very high solubility of DAP in these solvents. The molecular mechanism is quite clear. By blocking the nonpolar region, DAP is less prone to self-associate and form larger clusters, promoting dispersion rather than precipitation. Even a small molecule of acetone makes DAP dimerization difficult by forming very strong hydrogen bonds and partially occupying the nonpolar pocket.

<figure><img src="image_8.png" /><figcaption>**Figure 4.** The solute–solvent contact characteristics of DAP in polar proton-accepting solvents such as acetone and esters.</figcaption></figure>

The third class represents solvents that interact with DAP via non-HB interactions. These compounds are typical polar aprotic solvents such as those used in this study, namely
NMP, 4FM, and DMSO, as well as others taken from the literature, namely MeAc (methyl acetate), EtAc (ethyl acetate), iPrAc (isopropyl acetate), BuAc (butyl acetate), and EtPr (ethyl propionate). Interestingly, the ΔGr values calculated for the complexes formed by DAP and these solvents (Figure 5) generally do not differ significantly from the previously discussed classes.

<figure><img src="image_9.png" /><figcaption>**Figure 5.** The characteristics of solute–solvent contacts of DAP in polar non-proton-containing solvents in which dapsone interacts mainly via non-HB interactions.</figcaption></figure>
3.4. Ensemble Model for Solubility Prediction

Although the consonance solvent approach allows solubility calculations for dapsone despite the lack of fusion thermodynamic properties, the accuracy of this approach is far from acceptable. Therefore, machine learning techniques were used to develop a more reliable solubility prediction model. The methodology used in this study is similar to our previous work [47], although a new set of molecular descriptors was used. As discussed above, the solubility of dapsone can be rationalized in terms of intermolecular interactions in the saturated solutions. Consequently, the machine learning protocol included affinities as well as entropic and enthalpic contributions for regressor training. All molecular descriptors, predicted solubilities, and model details are provided in the supporting materials (see the File S2).

In particular, the focus of the regressor model tuning effort went beyond minimizing the deviation between computed and experimental data. The primary goal was to improve the predictive power of the developed model. To accomplish this, a custom loss function was employed that included a penalty estimated via learning curve analysis. While this approach may slightly compromise the fitting accuracy, it significantly improves the predictive potential of the model. From the 36 tunable models included in our Python code, a final set of nine regressor models was selected for ensemble construction. The inclusion criterion was derived from the learning curve analysis using five-fold cross-validation. Figure 6a illustrates the relationship between the area under the curve (AUC) for the validation and training subsets, clearly showing a distinct cluster of nine regressors. The order in which the models are listed reflects the increasing value of the AUC derived from the root mean square (RMS) plot obtained by systematically increasing the data portion from 50% to 100% using learning curve analysis. The ensemble of regressors developed in this study includes a diverse set of machine learning algorithms, each with unique characteristics. The aggregate consists of NuSVR, SVR, MLPRegressor, CatBoostRegressor, RandomForestRegressor, BaggingRegressor, HistGradientBoostingRegressor, KNeighborsRegressor, and ExtraTreeRegressor models. These models offer various learning strategies and approaches. Support vector machines (SVMs), such as NuSVR and SVR, excel at finding hyperplanes that separate data or predict continuous values, providing flexibility and adaptability for solving different regression problems. Random forests and extra trees employ ensemble learning techniques based on decision trees, while gradient boosting methods (e.g., Hist-GradientBoostingRegressor or CatBoostRegressor) sequentially build ensembles of weak learners. K-nearest neighbors considers data point proximity for predictions. The CatBoostRegressor is a powerful machine learning algorithm that is well suited for solubility modeling due to its ability to effectively detect complex relationships between molecular descriptors and solubility. The MLPRegressor takes advantage of the artificial neural networks, which provide a flexible infrastructure for capturing nonlinear relationships by allowing the specification of the number of hidden layers and their respective sizes. This flexibility enables the model to effectively represent complex patterns and dependencies present in the solubility data. To optimize the performance of all of these models, an extensive tuning process was performed using an Optuna study. This involved optimizing the hyperparameters of the models to identify the best-performing configurations. As a result, the ensemble exhibits diversity in terms of the hyperparameter settings chosen across the different models. By integrating these different models into an ensemble, one can leverage their individual strengths while mitigating their weaknesses. The use of an ensemble approach improves predictive performance by aggregating the predictions of multiple models, allowing a wider range of patterns to be captured and reducing the impact of individual model biases. In the case of the dapsone solubility model, the ensemble is constructed simply by averaging the predictions provided by each regressor without optimizing their weight. It is important to note that when predicting solubility expressed in terms of the decadal logarithm of the mole fraction, a formal constraint is imposed on the solution. Therefore, only negative predicted values are considered for averaging. Fortunately, the ensemble components consistently met this requirement. The developed
ensemble of regressors provides a robust framework that leverages the strengths of machine
learning algorithms, resulting in highly accurate predictions, as evidenced by the results
shown in Figure 6b. The performance of the aggregated regressors exceeds the accuracy
of solutions obtained using the consonance solvent approach. To illustrate the individual
performance of the regressors, the compilation is offered in the supporting materials (see
Section S2 in File S1). For illustrative purposes, the results of the best-performing model
are documented in Figure 7.

<figure><img src="image_10.png" /><figcaption>**Figure 6.** Results of machine learning analysis of dapsone solubility in pure solvents and binary mixtures. (a) AUC (area under the curve) derived from the LCA using 5-fold cross-validation. (b) Correlation between calculated and measured solubility data. The black circles indicate the predictions of the ensemble, while the open circles document the results of consonance solvation using four reference solvents (MeAc, iBuOH, TETA, and Ace).</figcaption></figure>

The power of even a single regressor to accurately predict dapsone solubility is evident. This is demonstrated not only by the near-perfect agreement between calculated and experimental values but also by the smooth trend observed in the learning curve analysis (LCA) plot. As the percentage of data included in the LCA increases, there is a systematic and gradual improvement in the resulting R2 and MEA, indicating the robustness and stability of the model predictions.
<figure><img src="image_11.png" /><figcaption>**Figure 7.** Graphical representation of the performance of the NuSVR regressor, which was found to be best suited for dapsone solubility calculations based on molecular descriptors derived from solute–solvent intermolecular interactions. The top panel shows a correlation between calculated and measured solubility data and the applicability domain plot. Bottom panel illustrates the results of learning curve analysis (LCA) and area under the curve (AUC) determination.</figcaption></figure>

4. Conclusions

The dissolution of dapsone, as an important aspect from both practical and theoretical points of view, can be rationalized based on its ability to form intermolecular complexes in saturated solutions. The structure of DAP suggests that this solute can interact with solvent molecules in a variety of ways. Solute–solvent contacts can be stabilized by proton acceptor and proton donor centers present in the DAP molecule, suggesting the promotion of dissolution in solvents, the structures of which provide a counterpart for hydrogen bonding. Moreover, broad delocalized electron clouds of aromatic rings are favorable regions for interactions with nonpolar solvents or those containing nonpolar fragments. It is therefore quite surprising that DAP is not too soluble in most organic solvents and their mixtures. The clue to this apparent paradox may be provided by the study of the self-association of two DAP molecules. The results presented suggest that the homomolecular contacts adopt a stacking conformation in all pure solvents analyzed and that this type of cluster is the most stable among all other types of binary complexes. The strong tendency to self-associate is due to nonpolar interactions similar to the close contacts observed in the solid state. There is a very high contribution from electron correlation, which further enhances the stability of DAP dimers. Thus, the high self-affinity of dapsone
is responsible for its aggregation, which eventually leads to precipitation regardless of the solvent. However, solvents that can disrupt dapsone stacking by at least partially blocking the apolar regions act as better dispersants and allow higher concentrations in saturated conditions. It is important to note that the elucidation of this mechanism may have been due to the fact that the solvent space was significantly expanded by providing new solubility data in such solvents. Without such new experimental data, the mechanism underlying the dissolution of dapsone would not be clear, as there are no accurate models to predict solubility values in solvents in which DAP has not been measured. This has been documented by COSMO-RS calculations, whose predictions have at most qualitative accuracy. Furthermore, other existing models are not helpful in this respect. For example, recently, an interesting approach using the 2D structural data (SMILES code) of the solvent and solute for the prediction of the solubility at molar concentrations has appeared [45,82]. In general, the concept of avoiding quantum chemical computation is very attractive due to its efficiency. Unfortunately, there are unsatisfactory high deviations (>100%) when comparing the online prediction results with new experimental data presented in Table 1.

Therefore, the ensemble of nonlinear models was defined based on nine regressors. The effectiveness of the Optuna study allows for hyperparameter tuning and tailoring of the ensemble of models to solve a particular problem. It is interesting to note that the values of intermolecular parameters such as contact affinities and solvability carried sufficient information to formulate a very accurate model of dapsone solubility in neat solvents and their binary mixtures. Despite the very high accuracy of the provided model, it has some disadvantages related to the set of molecular descriptors, since the characteristics of intermolecular interactions are computationally expensive. However, the value of the obtained information about the structure, affinities, and molecular descriptors is worth the effort. In principle, the developed model allows screening in any type of solvent if only the contact characteristics are provided.

**Supplementary Materials:** The following are available online at https://www.mdpi.com/article/10.3390/ma16186336/s1, (a) in File S1, dapsone solubility data curation and normalization and regressor model performance are provided; (b) in the File S2, all experimental and computed solubility data, molecular descriptors, and optimized parameters of nine regressors are collected.

**Author Contributions:** Conceptualization, P.C.; methodology, P.C.; validation, P.C., M.P. and T.J.; formal analysis, P.C., M.P. and T.J.; investigation, P.C., M.P. and T.J.; resources, P.C., M.P. and T.J.; data curation, P.C.; writing—original draft preparation, P.C., M.P. and T.J.; writing—review and editing, P.C., M.P. and T.J.; visualization, P.C.; supervision, P.C.; project administration, 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 conflict of interest.

## References

1. Madanipour, M.R.; Fatehi-zardalou, M.; Rahimi, N.; Hemmati, S.; Alaeddini, M.; Etemad-Moghadam, S.; Shayan, M.; Dabiri, S.; Dehpour, A.R. The anti-inflammatory effect of dapsone on ovalbumin-induced allergic rhinitis in balb/c mice. *Life Sci.* **2022**, *297*, 120449. [CrossRef] [PubMed]

2. Zhu, Y.I.; Stiller, M.J. Dapsone and sulfones in dermatology: Overview and update. *J. Am. Acad. Dermatol.* **2001**, *45*, 420–434. [CrossRef] [PubMed]

3. May, S.M.; Motosue, M.S.; Park, M.A. Dapsone is often tolerated in HIV-infected patients with history of sulfonamide antibiotic intolerance. *J. Allergy Clin. Immunol. Pract.* **2017**, *5*, 831–833. [CrossRef] [PubMed]

4. Wozel, G.; Blasum, C. Dapsone in dermatology and beyond. *Arch. Dermatol. Res.* **2014**, *306*, 103. [CrossRef]
5. Moreno, E.; Calvo, A.; Schwartz, J.; Navarro-Blasco, I.; González-Peñas, E.; Sanmartín, C.; Irache, J.; Espuelas, S. Evaluation of Skin Permeation and Retention of Topical Dapsone in Murine Cutaneous Leishmaniasis Lesions. *Pharmaceutics* **2019**, *11*, 607. [CrossRef]

6. Swain, S.S.; Paidesetty, S.K.; Dehury, B.; Sahoo, J.; Vedithi, S.C.; Mahapatra, N.; Hussain, T.; Padhy, R.N. Molecular docking and simulation study for synthesis of alternative dapsone derivative as a newer antileprosy drug in multidrug therapy. *J. Cell. Biochem.* **2018**, *119*, 9838–9852. [CrossRef]

7. Roman, C.; Dima, B.; Muyshont, L.; Schurmans, T.; Gilliaux, O. Indications and efficiency of dapsone in IgA vasculitis (Henoch-Schonlein purpura): Case series and a review of the literature. *Eur. J. Pediatr.* **2019**, *178*, 1275–1281. [CrossRef]

8. Ghaoui, N.; Hanna, E.; Abbas, O.; Kibbi, A.G.; Kurban, M. Update on the use of dapsone in dermatology. *Int. J. Dermatol.* **2020**, *59*, 787–795. [CrossRef]

9. Ríos, C.; Orozco-Suarez, S.; Salgado-Ceballos, H.; Mendez-Armenta, M.; Nava-Ruiz, C.; Santander, I.; Barón-Flores, V.; Caram-Salas, N.; Diaz-Ruiz, A. Anti-Apoptotic Effects of Dapsone After Spinal Cord Injury in Rats. *Neurochem. Res.* **2015**, *40*, 1243–1251. [CrossRef]

10. Tingle, M.; Mahmud, R.; Maggs, J.; Pirmohamed, M.; Park, B. Comparison of the metabolism and toxicity of dapsone in rat, mouse and man. *J. Pharmacol. Exp. Ther.* **1997**, *283*, 817–823.

11. Mitra, A.K.; Thummel, K.E.; Kalhorn, T.F.; Kharasch, E.D.; Unadkat, J.D.; Slattery, J.T. Metabolism of dapsone to its hydroxylamine by CYP2E1 in vitro and in vivo. *Clin. Pharmacol. Ther.* **1995**, *58*, 556–566. [CrossRef] [PubMed]

12. Molinelli, E.; Paolinelli, M.; Campanati, A.; Brisigotti, V.; Offidani, A. Metabolic, pharmacokinetic, and toxicological issues surrounding dapsone. *Expert Opin. Drug Metab. Toxicol.* **2019**, *15*, 367–379. [CrossRef] [PubMed]

13. Jouyban, A.; Rahimpour, E.; Karimzadeh, Z.; Zhao, H. Simulation of dapsone solubility data in mono- and mixed-solvents at various temperatures. *J. Mol. Liq.* **2022**, *345*, 118223. [CrossRef]

14. Schneider-Rauber, G.; Argenta, D.F.; Caon, T. Emerging Technologies to Target Drug Delivery to the Skin—The Role of Crystals and Carrier-Based Systems in the Case Study of Dapsone. *Pharm. Res.* **2020**, *37*, 240. [CrossRef]

15. Wu, Y.; Hao, X.; Li, J.; Guan, A.; Zhou, Z.; Guo, F. New insight into improving the solubility of poorly soluble drugs by preventing the formation of their hydrogen-bonds: A case of dapsone salts with camphorsulfonic and 5-sulfosalicylic acid. *CrystEngComm* **2021**, *23*, 6191–6198. [CrossRef]

16. Paredes da Rocha, N.; de Souza, A.; Nishitani Yukuyama, M.; Lopes Barreto, T.; Macedo, L.D.O.; Löbenberg, R.; Lima Barros de Araújo, G.; Ishida, K.; Araci Bou-Chacra, N. Highly water-soluble dapsone nanocrystals: Towards innovative preparations for an undermined drug. *Int. J. Pharm.* **2023**, *630*, 122428. [CrossRef]

17. Trombino, S.; Siciliano, C.; Procopio, D.; Curcio, F.; Laganà, A.S.; Di Gioia, M.L.; Cassano, R. Deep Eutectic Solvents for Improving the Solubilization and Delivery of Dapsone. *Pharmaceutics* **2022**, *14*, 333. [CrossRef]

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

19. Savjani, K.T.; Gajjar, A.K.; Savjani, J.K. Drug Solubility: Importance and Enhancement Techniques. *ISRN Pharm.* **2012**, *2012*, 195727. [CrossRef]

20. Yaseen, G.; Ahmad, M.; Zafar, M.; Akram, A.; Sultana, S.; Kilic, O.; Sonmez, G.D. Current status of solvents used in the pharmaceutical industry. In *Green Sustainable Process for Chemical and Environmental Engineering and Science*; Elsevier: Amsterdam, The Netherlands, 2021; pp. 195–219.

21. Parmentier, M.; Gabriel, C.M.; Guo, P.; Isley, N.A.; Zhou, J.; Gallou, F. Switching from organic solvents to water at an industrial scale. *Curr. Opin. Green Sustain. Chem.* **2017**, *7*, 13–17. [CrossRef]

22. Constable, D.J.C.; Jimenez-Gonzalez, C.; Henderson, R.K. Perspective on Solvent Use in the Pharmaceutical Industry. *Org. Process Res. Dev.* **2007**, *11*, 133–137. [CrossRef]

23. Lovette, M.A. Solubility Model to Guide Solvent Selection in Synthetic Process Development. *Cryst. Growth Des.* **2022**, *22*, 4404–4420. [CrossRef]

24. Modarresi, H.; Conte, E.; Abildskov, J.; Gani, R.; Crafts, P. Model-Based Calculation of Solid Solubility for Solvent Selection—A Review. *Ind. Eng. Chem. Res.* **2008**, *47*, 5234–5242. [CrossRef]

25. Moodley, K.; Rarey, J.; Ramjugernath, D. Model evaluation for the prediction of solubility of active pharmaceutical ingredients (APIs) to guide solid–liquid separator design. *Asian J. Pharm. Sci.* **2018**, *13*, 265–278. [CrossRef] [PubMed]

26. Klamt, A.; Eckert, F. Fast Solvent Screening via Quantum Chemistry: COSMO-RS Approach. *AIChE J.* **2002**, *48*, 369–385.

27. Klamt, A. Solvent-screening and co-crystal screening for drug development with COSMO-RS. *J. Cheminform.* **2012**, *4*, O14. [CrossRef]

28. 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**, 799–805. [CrossRef]

29. Palmelund, H.; Andersson, M.P.; Asgreen, C.J.; Boyd, B.J.; Rantanen, J.; Löbmann, K. Tailor-made solvents for pharmaceutical use? Experimental and computational approach for determining solubility in deep eutectic solvents (DES). *Int. J. Pharm. X* **2019**, *1*, 100034. [CrossRef]

30. Klajmon, M. Purely Predicting the Pharmaceutical Solubility: What to Expect from PC-SAFT and COSMO-RS. *Mol. Pharm.* **2022**, *19*, 4212–4232. [CrossRef]
31. Klamt, A.; Eckert, F.; Hornig, M.; Beck, M.E.; Bürger, T. Prediction of aqueous solubility of drugs and pesticides with COSMO-RS. *J. Comput. Chem.* **2002**, *23*, 275–281. [CrossRef]

32. Przybyłek, M.; Miernicka, A.; Nowak, M.; Cysewski, P. New Screening Protocol for Effective Green Solvents Selection of Benzamide, Salicylamide and Ethenzamide. *Molecules* **2022**, *27*, 3323. [CrossRef]

33. Loschen, C.; Klamt, A. Prediction of Solubilities and Partition Coefficients in Polymers Using COSMO-RS. *Ind. Eng. Chem. Res.* **2014**, *53*, 11478–11487. [CrossRef]

34. Buggert, M.; Cadenà, C.; Mokrushina, L.; Smirnova, I.; Maginn, E.J.; Arlt, W. COSMO-RS Calculations of Partition Coefficients: Different Tools for Conformation Search. *Chem. Eng. Technol.* **2009**, *32*, 977–986. [CrossRef]

35. Roy, D.; Patel, C. Revisiting the Use of Quantum Chemical Calculations in LogOctanol-water Prediction. *Molecules* **2023**, *28*, 801. [CrossRef]

36. Eckert, F.; Klamt, A. Accurate prediction of basicity in aqueous solution with COSMO-RS. *J. Comput. Chem.* **2006**, *27*, 11–19. [CrossRef] [PubMed]

37. Panić, M.; Radović, M.; Cvjetko Bubalo, M.; Radošević, K.; Rogošić, M.; Coutinho, J.A.P.; Radojčić, Redovniković, I.; Jurinjak Tušek, A. Prediction of pH Value of Aqueous Acidic and Basic Deep Eutectic Solvent Using COSMO-RS σ Profiles' Molecular Descriptors. *Molecules* **2022**, *27*, 4489. [CrossRef]

38. Andersson, M.P.; Jensen, J.H.; Stipp, S.L.S. Predicting pK a for proteins using COSMO-RS. *PeerJ* **2013**, *1*, e198. [CrossRef]

39. Guidetti, M.; Hilfiker, R.; Kuentz, M.; Bauer-Brandl, A.; Blatter, F. Exploring the Cocrystal Landscape of Posaconazole by Combining High-Throughput Screening Experimentation with Computational Chemistry. *Cryst. Growth Des.* **2023**, *23*, 842–852. [CrossRef]

40. Deng, Y.; Liu, S.; Jiang, Y.; Martins, I.C.B.; Rades, T. Recent Advances in Co-Former Screening and Formation Prediction of Multicomponent Solid Forms of Low Molecular Weight Drugs. *Pharmaceutics* **2023**, *15*, 2174. [CrossRef]

41. Przybyłek, M.; Ziółkowska, D.; Mroczyńska, K.; Cysewski, P. Applicability of Phenolic Acids as Effective Enhancers of Cocrystal Solubility of Methylxanthines. *Cryst. Growth Des.* **2017**, *17*, 2186–2193. [CrossRef]

42. Li, C.; Wu, D.; Li, J.; Ji, X.; Qi, L.; Sun, Q.; Wang, A.; Xie, C.; Gong, J.; Chen, W. Multicomponent crystals of clotrimazole: A combined theoretical and experimental study. *CrystEngComm* **2021**, *23*, 6977–6993. [CrossRef]

43. Li, J.; Wu, D.; Xiao, Y.; Li, C.; Ji, X.; Sun, Q.; Chang, D.; Zhou, L.; Jing, D.; Gong, J.; et al. Salts of 2-hydroxybenzylamine with improvements on solubility and stability: Virtual and experimental screening. *Eur. J. Pharm. Sci.* **2022**, *169*, 106091. [CrossRef]

44. Lee, S.; Lee, M.; Gyak, K.-W.; Kim, S.D.; Kim, M.-J.; Min, K. Novel Solubility Prediction Models: Molecular Fingerprints and Physicochemical Features vs Graph Convolutional Neural Networks. *ACS Omega* **2022**, *7*, 12268–12277. [CrossRef] [PubMed]

45. Panapitiya, G.; Girard, M.; Hollas, A.; Sepulveda, J.; Murugesan, V.; Wang, W.; Saldanha, E. Evaluation of Deep Learning Architectures for Aqueous Solubility Prediction. *ACS Omega* **2022**, *7*, 15695–15710. [CrossRef] [PubMed]

46. Vermeire, F.H.; Chung, Y.; Green, W.H. Predicting Solubility Limits of Organic Solutes for a Wide Range of Solvents and Temperatures. *J. Am. Chem. Soc.* **2022**, *144*, 10785–10797. [CrossRef]

47. 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]

48. Cysewski, P.; Jeliński, T.; Przybyłek, M.; Nowak, W.; Olczak, M. Solubility Characteristics of Acetaminophen and Phenacetin in Binary Mixtures of Aqueous Organic Solvents: Experimental and Deep Machine Learning Screening of Green Dissolution Media. *Pharmaceutics* **2022**, *14*, 2828. [CrossRef]

49. Cysewski, P.; Przybyłek, M.; Rozalski, R. Experimental and Theoretical Screening for Green Solvents Improving Sulfamethizole Solubility. *Materials* **2021**, *14*, 5915. [CrossRef]

50. Cysewski, P.; Jeliński, T.; Cymerman, P.; Przybyłek, M. Solvent Screening for Solubility Enhancement of Theophylline in Neat, Binary and Ternary NADES Solvents: New Measurements and Ensemble Machine Learning. *Int. J. Mol. Sci.* **2021**, *22*, 7347. [CrossRef]

51. Jeliński, T.; Bugalska, N.; Koszucka, K.; Przybyłek, M.; Cysewski, P. Solubility of sulfanilamide in binary solvents containing water: Measurements and prediction using Buchowski-Ksiazczak solubility model. *J. Mol. Liq.* **2020**, *319*, 114342. [CrossRef]

52. Cysewski, P.; Przybyłek, M.; Kowalska, A.; Tymorek, N. Thermodynamics and intermolecular interactions of nicotinamide in neat and binary solutions: Experimental measurements and COSMO-RS concentration dependent reactions investigations. *Int. J. Mol. Sci.* **2021**, *22*, 7365. [CrossRef] [PubMed]

53. Przybyłek, M.; Kowalska, A.; Tymorek, N.; Dziaman, T.; Cysewski, P. Thermodynamic Characteristics of Phenacetin in Solid State and Saturated Solutions in Several Neat and Binary Solvents. *Molecules* **2021**, *26*, 4078. [CrossRef] [PubMed]

54. 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]

55. Li, W.; Ma, Y.; Yang, Y.; Xu, S.; Shi, P.; Wu, S. Solubility measurement, correlation and mixing thermodynamics properties of dapsone in twelve mono solvents. *J. Mol. Liq.* **2019**, *280*, 175–181. [CrossRef]

56. Li, H.; Xie, Y.; Xue, Y.; Zhu, P.; Zhao, H. Comprehensive insight into solubility, dissolution properties and solvation behaviour of dapsone in co-solvent solutions. *J. Mol. Liq.* **2021**, *341*, 117403. [CrossRef]

57. Shi, Y.; Wang, S.; Wang, J.; Liu, T.; Qu, Y. Solubility Determination and Thermodynamic Modeling of Amitriptyline Hydrochloride in 13 Pure Solvents at Temperatures of 283.15–323.15 K. *J. Chem. Eng. Data* **2021**, *66*, 1877–1889. [CrossRef]
58. Wang, S.; Cheng, X.; Liu, B.; Du, Y.; Wang, J. Temperature dependent solubility of sodium cyclamate in selected pure solvents and binary methanol+water mixed solvents. *Fluid Phase Equilib.* **2015**, *390*, 1–6. [CrossRef]

59. Liang, A.; Wang, S.; Qu, Y. Determination and Correlation of Solubility of Phenylbutazone in Monosolvents and Binary Solvent Mixtures. *J. Chem. Eng. Data* **2017**, *62*, 864–871. [CrossRef]

60. Yang, X.; Wang, J.; Ma, M.; Qu, Y. Determination and Modeling of Artesunate Solubility in 13 Pure Solvents at 283.15–323.15 K. *J. Chem. Eng. Data* **2022**, *67*, 3734–3747. [CrossRef]

61. Yang, X.; Wang, S.; Wang, J. Measurement and Correlation of Solubility of Loratadine in Different Pure Solvents and Binary Mixtures. *J. Chem. Eng. Data* **2017**, *62*, 391–397. [CrossRef]

62. Galaon, T.; David, V. Deviation from van't Hoff dependence in RP-LC induced by tautomeric interconversion observed for four compounds. *J. Sep. Sci.* **2011**, *34*, 1423–1428. [CrossRef] [PubMed]

63. Dassault Systèmes. *Biovia COSMOtherm*, version 22.0.0; Dassault Systèmes: Vélizy-Villacoublay, France, 2022.

64. Cysewski, P. Application of the consonance solvent concept for accurate prediction of buckminster solubility in 180 net solvents using COSMO-RS approach. *Symmetry* **2019**, *11*, 828. [CrossRef]

65. Dassault Systèmes. *Biovia COSMOconf*, version 20.0.0; Dassault Systèmes: Vélizy-Villacoublay, France, 2020.

66. TURBOMOLE GmbH. *TURBOMOLE*, version 7.5.1; Dassault Systèmes: Vélizy-Villacoublay, France, 2020.

67. Dassault Systèmes. *Biovia TmoleX*, version 21.0.1; Dassault Systèmes: Vélizy-Villacoublay, France, 2020.

68. Jeliński, T.; Kubsik, M.; Cysewski, P. Application of the Solute-Solvent Intermolecular Interactions as Indicator of Caffeine Solubility in Aqueous Binary Aprotic and Proton Acceptor Solvents: Measurements and Quantum Chemistry Computations. *Materials* **2022**, *15*, 2472. [CrossRef] [PubMed]

69. Grimme, S.; Antony, J.; Ehrlich, S.; Krieg, H. A consistent and accurate ab initio parametrization of density functional dispersion correction (DFT-D) for the 94 elements H-Pu. *J. Chem. Phys.* **2010**, *132*, 154104. [CrossRef]

70. Witte, J.; Neaton, J.B.; Head-Gordon, M. Effective empirical corrections for basis set superposition error in the def2-SVPD basis: gCP and DFT-C. *J. Chem. Phys.* **2017**, *146*, 234105. [CrossRef]

71. 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, Anchorage, AK, USA, 4–8 August 2019; pp. 2623–2631.

72. Ueda, H.; Osaki, H.; Miyano, T. Baloxavir Marboxil Shows Anomalous Conversion of Crystal Forms from Stable to Metastable through Formation of Specific Solvate Form. *J. Pharm. Sci.* **2023**, *112*, 158–165. [CrossRef]

73. Do, H.T.; Chua, Y.Z.; Kumar, A.; Pabsch, D.; Hallermann, M.; Zaitsau, D.; Schick, C.; Held, C. Melting properties of amino acids and their solubility in water. *RSC Adv.* **2020**, *10*, 44205–44215. [CrossRef]

74. Tsioptsias, C.; Tsivintzelis, I. On the Thermodynamic Thermal Properties of Quercetin and Similar Pharmaceuticals. *Molecules* **2022**, *27*, 6630. [CrossRef]

75. Braun, D.E.; Krüger, H.; Kahlenberg, V.; Griesser, U.J. Molecular level understanding of the reversible phase transformation between forms III and II of dapsone. *Cryst. Growth Des.* **2017**, *17*, 5054–5060. [CrossRef]

76. Braun, D.E.; Vickers, M.; Griesser, U.J. Dapsone Form V: A Late Appearing Thermodynamic Polymorph of a Pharmaceutical. *Mol. Pharm.* **2019**, *16*, 3221–3236. [CrossRef]

77. Keyvanpour, M.R.; Shirzad, M.B. An Analysis of QSAR Research Based on Machine Learning Concepts. *Curr. Drug Discov. Technol.* **2021**, *18*, 17–30. [CrossRef] [PubMed]

78. Wang, J.; Xu, P.; Ji, X.; Li, M.; Lu, W. Feature Selection in Machine Learning for Perovskite Materials Design and Discovery. *Materials* **2023**, *16*, 3134. [CrossRef] [PubMed]

79. How, W.B.; Wang, B.; Chu, W.; Tkatchenko, A.; Prezhdo, O.V. Significance of the Chemical Environment of an Element in Nonadiabatic Molecular Dynamics: Feature Selection and Dimensionality Reduction with Machine Learning. *J. Phys. Chem. Lett.* **2021**, *12*, 12026–12032. [CrossRef] [PubMed]

80. Zhang, K.; Zhang, H. Machine Learning Modeling of Environmentally Relevant Chemical Reactions for Organic Compounds. *ACS EST Water* **2022**. [CrossRef]

81. Dhal, P.; Azad, C. A comprehensive survey on feature selection in the various fields of machine learning. *Appl. Intell.* **2022**, *52*, 4543–4581. [CrossRef]

82. RMG. Available online: https://rmg.mit.edu/database/solvation/searchSolubility/ (accessed on 4 September 2023).

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