Increasing the Equilibrium Solubility of Meloxicam in Aqueous Media by Using Dimethyl Sulfoxide as a Cosolvent: Correlation, Dissolution Thermodynamics and Preferential Solvation

: Meloxicam is widely prescribed as an analgesic and anti-inﬂammatory drug in human therapeutics. Owing the very low aqueous solubility of meloxicam, this property has been studied in dimethyl sulfoxide (DMSO)-aqueous solvent systems at several temperatures from 273.15 to 313.15 K to expand the solubility database about analgesic drugs in mixed solvents. The ﬂask shake method followed by ultraviolet-visible (UV-vis) spectrophotometry analysis were used for meloxicam solubility determinations. A number of cosolvency models, including the Jouyban–Acree model, were challenged for solubility correlation/prediction of this drug in these mixtures. The van’t Hoff and Gibbs equations were employed to calculate the apparent standard thermodynamic quantities relative to dissolution and mixing processes. The inverse Kirkwood–Buff integral method was employed for calculating the preferential solvation parameters of meloxicam by DMSO in the mixtures. Meloxicam solubility increases with increasing temperature and maximum solubilities are observed in neat DMSO at all temperatures studied. Dissolution processes were endothermic in all cases and entropy-driven in the composition interval of 0.40 ≤ x 1 ≤ 1.00. A nonlinear enthalpy–entropy relationship was observed in the plot of enthalpy vs. Gibbs energy for drug transfer processes. Meloxicam is preferentially solvated by water in water-rich mixtures but preferentially solvated by DMSO in the composition interval of 0.21 < x 1 < 1.00.


Introduction
Meloxicam (molecular structure shown in Figure 1, IUPAC name: 4-hydroxy-2-methyl-N-(5-methyl-2-thiazolyl)-2H-1,2-benzothiazine-3-carboxamide-1,1-dioxide, molar mass 351.40 g·mol −1 , CAS number: 71125-38-7, PubChem CID: 54677470) is a non-steroidal anti-inflammatory drug used commonly for pain and inflammatory treatments [1][2][3][4][5]. From a physicochemical point of view, meloxicam exhibits very low aqueous solubility, which influences negatively in vivo dissolution rates, affecting its biological performance. Otherwise, because of the very low aqueous solubility of this drug, all the duties relative to research and development of homogeneous liquid dosage forms, such as peroral or Otherwise, because of the very low aqueous solubility of this drug, all th to research and development of homogeneous liquid dosage forms, su injectable products, based on this drug, are very long and hard at an ind ceutical level. In order to overcome the drawbacks mentioned above, som have been intended to increase the aqueous equilibrium solubility of m investigations were mainly based in the evaluation of the solubilizing eff mon pharmaceutical cosolvents, as has recently been summarized in a p nication of our research group [6]. More recently, some other aqueous mi different cosolvents, including choline-based deep eutectic solvents, hav ied and reported [7][8][9][10][11]. It is noteworthy that very good meloxicam solu has been reported, reaching more than 1000-fold in some cases, as follow bitol ® [9], 1399 with N-methylformamide [6], and 11,233 with N,N-dim [6]. Dimethyl sulfoxide (molar mass: 78.13 g·mol −1 , CAS number: 67-68-679) is a polar, aprotic solvent, miscible with water and with a wide rang vents in all possible compositions, which makes it adequate for dissolvi both polar and nonpolar nature [1,12]. It exhibits less toxicity than N-m and N-methyl-2-pyrrolidone and has been reported about its good pow the skin and other membranes without damaging them. This is the rea high potential to increase the penetration of less soluble active ingredi DMSO exhibits by itself some analgesic and anti-inflammatory propert there are some references in the literature that suggest the potential us possible pharmacotherapeutic agent in the management of pain and othe Moreover, owing its high solubilizing power for different active pharma ents, it has been used as an excipient in topical and parenteral medication veterinary use [15]. On the other hand, DMSO has also been considered in medicinal chemistry for studying intermolecular effects of different d way, aqueous mixtures of DMSO have been studied for solubilizing seve alike compounds, and other organic chemicals, including bergenin [17 dinitramide [18], naringin [19], sinapic acid [20], p-nitrobenzamide [21], micoflavin [23], phenformin [24], baricitinib [25], a pyridazinone derivat otinamide [27], among others.
As it is well-described in the specialized chemical and pharmaceuti the physicochemical data about the equilibrium solubility of drugs or pounds in aqueous cosolvent mixtures, as well, as the deep understandi tive dissolution mechanisms, are very important for both theoretical and of view in pharmaceutical and chemical sciences. This is because the mea and analyzed solubility values expand the respective solubility database useful for theoretical and practical purposes in both the pharmaceutical Dimethyl sulfoxide (molar mass: 78.13 g·mol −1 , CAS number: 67-68-5, PubChem CID 679) is a polar, aprotic solvent, miscible with water and with a wide range of organic solvents in all possible compositions, which makes it adequate for dissolving substances of both polar and nonpolar nature [1,12]. It exhibits less toxicity than N-methylformamide and N-methyl-2-pyrrolidone and has been reported about its good power of penetrating the skin and other membranes without damaging them. This is the reason why it has a high potential to increase the penetration of less soluble active ingredients. Otherwise, DMSO exhibits by itself some analgesic and anti-inflammatory properties [13]. Indeed, there are some references in the literature that suggest the potential use of DMSO as a possible pharmacotherapeutic agent in the management of pain and other conditions [14]. Moreover, owing its high solubilizing power for different active pharmaceutical ingredients, it has been used as an excipient in topical and parenteral medications for human and veterinary use [15]. On the other hand, DMSO has also been considered as solvent model in medicinal chemistry for studying intermolecular effects of different drugs [16]. In this way, aqueous mixtures of DMSO have been studied for solubilizing several drugs, drug-alike compounds, and other organic chemicals, including bergenin [17], N-guanylurea dinitramide [18], naringin [19], sinapic acid [20], p-nitrobenzamide [21], d-histidine [22], micoflavin [23], phenformin [24], baricitinib [25], a pyridazinone derivative [26], and nicotinamide [27], among others.
As it is well-described in the specialized chemical and pharmaceutical literature, all the physicochemical data about the equilibrium solubility of drugs or drug-alike compounds in aqueous cosolvent mixtures, as well, as the deep understanding of the respective dissolution mechanisms, are very important for both theoretical and practical points of view in pharmaceutical and chemical sciences. This is because the measured, reported, and analyzed solubility values expand the respective solubility databases, which is very useful for theoretical and practical purposes in both the pharmaceutical and chemical industries, as indicated above [28][29][30][31].
Therefore, the main aims of this research were as follows: (i) to determine and analyze the effects of both mixtures' composition and temperature on the solubility of meloxicam in {DMSO (1) + water (2)} mixtures; (ii) to correlate equilibrium solubility data with several well-known thermodynamic models; (iii) to calculate the apparent standard dissolution and mixing thermodynamic parameters; and (iv) to study the preferential solvation parameters of meloxicam in binary mixtures conformed by DMSO and water. Therefore, this research is a continuation of some other similar ones reported earlier in the literature [6][7][8][9]11,32] about the meloxicam equilibrium solubility in other aqueous cosolvent systems of pharmaceutical interest for design of dosage forms and its quality control analysis, but using another commonly used cosolvent because of its high solubilizing power.

Preparation of Solvent Mixtures
All the {DMSO (1) + water (2)} binary solvent mixtures were prepared gravimetrically by using an Ohaus Pioneer TM PA214 analytical balance (USA, sensitivity ± 0.1 mg), in quantities of 50.00 g. The mole fractions of DMSO of the nine mixtures prepared, varied by 0.10 in steps from x 1 = 0.10 to x 1 = 0.90.

Solubility Determinations
Equilibrium meloxicam solubilities were determined by using the shake-flask method [33], followed by UV-spectrophotometric analysis, as follows: an excess amount of meloxicam was added to 50.00 g of each binary solvent mixture or neat DMSO in dark glass pharmaceutical flasks. The stoppered flasks were putted in an ultrasonic bath (Elma ® E60H Elmasonic, Fremont, CA, USA) during 15 min and were later transferred to thermostatic mechanical shakers (Julabo SW23, Seelbach, Germany) or re-circulating thermostatic baths (Neslab RTE 10 Digital One Thermo Electron Company, Waltham, MA, USA) and kept at 313.15 K for at least four days to ensure that the drug saturation had been achieved. After that, the supernatant solutions were isothermally filtered (Millipore Corp. Swinnex ® -13, Burlington, MA, USA) to remove undissolved solid particles before sampling. Meloxicam concentrations were determined after appropriate gravimetric dilution with a 0.10 mol·dm −3 NaOH solution, by measuring the UV light absorbance at the maximum absorbance wavelength, λ max = 361 nm (UV/VIS BioMate 3 Thermo Electron Company spectrophotometer, USA), followed by interpolation from a previously validated UV spectrophotometric gravimetric calibration curve prepared in NaOH 0.10 mol·dm −3 . The respective equation was: Absorbance = 0.0073 + 52.508·C, with C expressed as µg·g −1 . Later, the thermostatic baths temperature was decreased from 313.15 K to 308.15 K allowing the meloxicam excess precipitation during two days, following with the same procedures mentioned above to determine the new meloxicam concentrations at saturation. All these procedures were performed successively until solid-liquid equilibrium was achieved at 293.15 K. All the solubility experiments were performed at least three times and the respective results were averaged. The density of the saturated solutions was measured by using a digital density meter (DMA 45 Anton Paar, Graz, Austria) connected to a re-circulating thermostatic bath (Neslab RTE 10 Digital One Thermo Electron Company, USA) in order to transform the obtained gravimetric solubility values into volumetric concentration scales. The density meter was calibrated at every temperature by using air and water as standards as indicated in the respective instruction manuals [34].

X-ray Diffraction (XRD) Analysis
To determine the crystal nature of the solid meloxicam samples both before and after the saturation in neat water, in the mixture of x 1 = 0.50, and in neat DMSO, the respective X-ray powder diffraction analyses were performed by using a PANalytical Xpert Pro X-ray diffractometer. The equipment is provided with CuKα radiation λ = 1.5418 Å. Generator setting: 40 kV and 40 mA and Bragg-Brentano geometry. Data were collected at 2θ from 5 • to 70 • and angle variation of 0.02 • with detector data acquisition time of 9.46 min operating under room temperature.

Fourier Transform Infrared (FTIR) Analysis
In additional to XRD analyses, in order to confirm the nature of the solid meloxicam samples, both before and after the saturation in neat water, in the mixture of x 1 = 0.50, and in neat DMSO, FTIR analyses were also performed. The meloxicam solid samples were ground with quantities from 10 to 100 times its bulk of pure potassium bromide and the resulting mixtures were pressed into discs by using a special mold and a manual hydraulic press (Specac ® , Fort Washington, PA, USA). The respective spectra were obtained in an FTIR spectrophotometer (IRAffinity-1, Shimadzu, Kyoto, Japan).

Experimental Mole Fraction and Molarity Solubility
Tables 1 and 2 summarize the experimental equilibrium solubilities of meloxicam in all {DMSO (1) + water (2)} solvent systems at 293.15 ≤ T/K ≤ 313.15, as expressed in mole fraction and molarity (mol·dm −3 ), respectively. The studied temperature interval includes what is commonly known as "room temperature" for products storage, as well as, the normal human body temperature. It is worth mentioning that all solubility values in neat water were taken from reference [32]. The average relative uncertainty obtained in the reported solubility values was 2.4%, which is in agreement with that commonly observed in this kind of experiments. If the mole fraction scale is considered, at T = 298.15 K, Table 1 shows that the meloxicam solubility increased 6956 times from x 3 = 1.137 × 10 −6 in neat water to x 3 = 7.909 × 10 −3 in neat DMSO, where maximum solubility is obtained. A deep comparison about the meloxicam solubility in neat water has been reported and discussed earlier in one of our previous communications [6]. Regarding the meloxicam mole fraction equilibrium solubility in neat DMSO Sathesh-Babu et al. reported a value of x 3 = 5.496 × 10 −3 at T = 298.15 K [35], which is in good agreement with the one obtained in this research (i.e., x 3 = 7.909 × 10 −2 , Table 1). Moreover, when considering the molarity scale, a value of C = 1.516 × 10 −2 mol·dm −3 was reported by Castro et al. at T = 298.15 K [36], which differs in almost one order of magnitude regarding our value (i.e., C = 0.1089 mol·dm −3 , Table 2). The observed differences in solubility values could be attributed to several reasons such as different polymorphic states, different saturation times, or different analytical procedures, among others, as described earlier [31]. Up to the best of our knowledge, no solubility values of meloxicam in aqueous mixtures of DMSO have been reported and no more comparisons are possible. Figure 2 depicts the meloxicam solubility profiles as function of the Hildebrand solubility parameters (δ 1+2 ) of {DMSO (1) + water (2)} mixtures at T = 298.15 K. As widely described, δ 1+2 is a very important polarity descriptor of cosolvent mixtures [28][29][30][31]. This descriptor was calculated considering the Hildebrand solubility parameter of both pure solvents (δ 1 = 26.6 MPa 1/2 for DMSO and δ 2 = 47.8 MPa 1/2 for water [37,38]), and the volume fraction (f i ) of each solvent, as [29,39]:    (3). Mean uncertainty in x 1 , u(x 1 ) = 0.0005. b Standard uncertainty in p is u(p) = 3.0 kPa. Average relative uncertainty in C, u r (C) = 0.024. Standard uncertainty in T is u(T) = 0.10 K. c Data taken from Delgado et al. [32].

Solid Phases' Analyses
X-ray diffraction (XRD) spectra for meloxicam of the original untreated sam after its saturation in neat water, neat DMSO, and the aqueous-DMSO mixture of are shown in Figure 3. Because of the high similarity among all obtained spectra, be concluded that changes of the crystalline form of meloxicam are not observed dissolution and saturation in these four solvent systems. Moreover, Table 3 sum the position, 2q spacing, peak height (in counts), and the relative intensity of peak iting values higher than 10% for the original untreated meloxicam sample. These are in good coincidence with those reported by Wu et al. [43]. Moreover, all the o XRD spectra of this research are very similar to that reported earlier for the poly of meloxicam [32,[44][45][46][47]. Finally, the FTIR spectra of all solid meloxicam samples in Figure 4 are also coincident with those reported in the literature, which allows cate that the three bottom-solid phases, obtained after drug saturation, have the s ture as the original untreated sample [47][48][49]. Therefore, meloxicam did not suffe polymorphic transitions or solvates formation after saturation in these dissolution iments.

Solid Phases' Analyses
X-ray diffraction (XRD) spectra for meloxicam of the original untreated sample and after its saturation in neat water, neat DMSO, and the aqueous-DMSO mixture of x 1 = 0.50 are shown in Figure 3. Because of the high similarity among all obtained spectra, it could be concluded that changes of the crystalline form of meloxicam are not observed after its dissolution and saturation in these four solvent systems. Moreover, Table 3 summarizes the position, 2q spacing, peak height (in counts), and the relative intensity of peaks exhibiting values higher than 10% for the original untreated meloxicam sample. These values are in good coincidence with those reported by Wu et al. [43]. Moreover, all the obtained XRD spectra of this research are very similar to that reported earlier for the polymorph I of meloxicam [32,[44][45][46][47]. Finally, the FTIR spectra of all solid meloxicam samples shown in Figure 4 are also coincident with those reported in the literature, which allows to indicate that the three bottom-solid phases, obtained after drug saturation, have the same nature as the original untreated sample [47][48][49]. Therefore, meloxicam did not suffer crystal polymorphic transitions or solvates formation after saturation in these dissolution experiments. Liquids 2022, 2, FOR PEER REVIEW 7          Table 4 summarizes the asymmetrical activity coefficients (γ 3 ) of meloxicam in all the {DMSO (1) + water (2)} solvent systems. These values were calculated as the quotient x id 3 /x 3 from the ideal (x id 3 ) and experimental solubilities of meloxicam summarized in Table 1. It is noteworthy that ideal solubilities were taken from the literature [32]. As observed, the γ 3 values vary from 2708 in neat water (where the lower drug solubilities are observed) to 0.389 in neat DMSO at T = 298.15 K, where the maximum drug solubility is observed at this temperature. At all temperatures, meloxicam exhibits γ 3 values higher than unity in neat water and the mixtures of x 1 ≤ 0.70, but lower than unity in DMSO-rich mixtures and neat DMSO. On the other hand, in neat water and the mixtures of x 1 = 0.10 and 0.20, the γ 3 values increase with an increase in temperature. This implies some distancing from the ideal dissolution behavior with increasing temperature. In mixtures of 0.30 ≤ x 1 ≤ 0.50, the γ 3 values are almost independent on temperature. On the contrary, in mixtures of x 1 ≥ 0.60 and neat DMSO, the γ 3 values decrease with the increase in temperature.  (3). Mean uncertainty in x 1 , u(x 1 ) = 0.0005. b Standard uncertainty in p is u(p) = 3.0 kPa. Average relative uncertainty in γ 3 is u r (γ 3 ) = 0.027. Standard uncertainty in T is u(T) = 0.10 K. c Data taken from Delgado et al. [32].

Activity Coefficients in Mixed Solvents
Moreover, Equation (2) allows a rough estimate of the magnitudes of solute-solvent intermolecular interactions from γ 3 values [50].
where subscript 1 stands for the solvent system that here corresponds to neat solvents or aqueous DMSO binary mixtures, e ss, e 33 , and e s3 represent the magnitudes of solventsolvent, solute-solute, and solvent-solute interaction energies, respectively. However, it is important to keep in mind that in the case of ternary systems, such as DMSO-watermeloxicam studied here, some water-cosolvent interactions are present, which could also play an important role in drug solubilities and dissolution rates. V 3 denotes the molar volume of the super-cooled liquid meloxicam, whereas, ϕ s denotes the volume fraction of every solvent system. For low x 3 solubility values, V 3 ϕ s 2 /RT may be considered as constant despite of the solvent system. Hence, γ 3 values would depend mainly on e ss, e 33 , and e s3 [50]. As is well-known, e ss and e 33 are unfavorable for drug solubility and dissolution, whereas e s3 favors the respective drug solubility and dissolution rate increasing.
The contribution of e 33 could be considered as constant in the different solvent systems studied. Thus, from a qualitative viewpoint, based on the energetic quantities described in Equation (2), the following analysis could be established: because e ss is highest in neat water (δ 2 = 47.8 MPa 1/2 ) and lowest in neat DMSO (δ 1 = 26.6 MPa 1/2 ) [37,38], neat water and water-rich mixtures (exhibiting γ 3 values higher than 2300) would imply high e ss and low e s3 values, whereas, in DMSO-rich mixtures (exhibiting γ 3 values near the unity), the e ss values are relatively low, and therefore, the e s3 values would be high. In this way, a higher solvation of meloxicam by DMSO in DMSO-rich mixtures is expected.

Solubility Modeling
Among the available cosolvency models presented for calculation of drug solubilities in mixed solvents at isothermal condition or at various temperatures [51,52], the Yalkowsky model is the simplest one [53] and requires only two experimental solubility determinations (in neat solvents) to predict the solubility at other solvent compositions. The Yalkowsky model is commonly represented as: where N is the number of experimental data points. As mentioned above, Equation (3) is capable of estimating drug solubility in cosolvent mixtures at individual T using only the drug solubility data in the mono-solvents at this T. However, it could be extended to obtain: to be applied at various temperatures (x 3(1+2),T ) using a single equation. In Equation (5), A and B terms are the model constants [54]. The trained model for solubility of meloxicam in {DMSO (1) + water (2)} mixtures is: which resulted in the MPD of 44.4%. The main limitation of the Yalkowsky model is that it does not consider any more interaction term after mixing the solutions and considers the mixing behavior as an ideal one. The Jouyban-Acree model includes as many as required interaction terms (J i terms) to describe the non-ideality of the mixing and produced the most accurate results in correlating drug solubility data in binary solvent mixtures at various temperatures. The model is presented as [51]: Liquids 2022, 2

170
where J i terms are the respective model constants that are computed by using a nonintercept least square analysis [31]. Accordingly, the generated solubility values of meloxicam in {DMSO (1) + water (2)} were fitted to Equation (7) and the obtained trained model was: The F value of Equation (8) was 791 and the adjusted correlation coefficient (R) was 0.979, whereas the correlation and the model constants were significant with p < 0.0005. Equation (8) is valid for calculating the solubility of meloxicam in different {DMSO (1) + water (2)} mixtures at various temperatures or mixtures-composition of interest, by employing the solubility data of meloxicam in neat DMSO and neat water at each temperature. The obtained MPD for the back-calculated solubility data of meloxicam when using Equation (8) was 9.6%.
Although Equation (8) provided an accurate correlation for the solubility of meloxicam in these {DMSO (1) + water (2)} mixtures, it requires the experimental solubility data in mono-solvents (i.e., x 3(1),T and x 3(2),T ) at any temperature of interest to calculate the solubility of meloxicam in the required binary solvent mixtures. However, one may combine the trained version of Equation (5) with Equation (7) to provide a full predictive model, obtaining the following: Equation (9) allows to calculate the solubilities of meloxicam in all these binary mixtures at various temperatures with an MPD of 9.9% and it does not require any experimental input data. The F and R values for Equation (9) were 787 and 0.978, respectively. For practical applications of Equation (9), one may train the model using the minimum number of seven experimental solubility points, and then, predict the rest of required data in any aqueous-DMSO mixture composition and temperature of interest, as has been exemplified in the literature [55]. When the model was trained with the solubility data in DMSO and water at T = (293.15 and 313.15) K (i.e., the lowest and highest temperatures) and in x 1 = 0.3, 0.5 and 0.7 at T = 298.15 K (total 7 data points), the rest of the data points were predicted with the MPD of 18.7% (N = 48).
In a previous work [56], generally trained version of the Jouyban-Acree-Abraham and Jouyban-Acree-Hansen models were presented for predicting the solubility of meloxicam in various binary solvent mixtures. These models are: and ln x m,T = w 1 ln x 1,T + w 2 ln x 2,T The c, e, s, a, b, and v are Abraham solvent's coefficients and δ d1 , δ p1 , δ h1 , and δ d2 , δ p2 , and δ h2 are the Hansen parameters for solvents 1 and 2, respectively [56]. Equations (10) and (11) predicted the solubility of meloxicam in (DMSO + water) mixtures with the MPDs of 74.3 and 35.6%. Although the prediction errors are relatively large, these equations only require the solubility data in the neat solvent.

Apparent Thermodynamic Functions of Dissolution
All the apparent standard thermodynamic quantities relative to meloxicam dissolution processes were calculated at the mean harmonic temperature, T hm = 303.0 K, which was calculated by using Equation (12) [57].
where n is the number of temperatures under study (i.e., 5 in this case). Hence, all the apparent standard enthalpy changes relative to dissolution processes (∆ soln H • ) were obtained by means of the modified van't Hoff equation, as shown in Equation (13) [58]: The apparent standard Gibbs energy changes relative to all the meloxicam dissolution processes (∆ soln G • ) were calculated by means of [58][59][60]: The intercepts used in Equation (14) were those obtained as the result of the linear regressions of ln x 3 as function of (1/T − 1/T hm ). Therefore, Figure 5 depicts the meloxicam solubility linear van't Hoff behavior in all the {DMSO (1) + water (2)} mixtures, as well as in both neat solvents. It is noteworthy that linear regressions with r 2 > 0.993 were observed in all the solvent systems [61][62][63].
Finally, the apparent standard changes in entropy, for all the studied meloxicam dissolution processes (∆ soln S • ) were calculated based on the respective ∆ soln H • and ∆ soln G • values by using [59,60]: Table 5 summarizes all the apparent standard thermodynamic quantities for the dissolution processes of meloxicam in all the {DMSO (1) + water (2)} mixtures at T hm = 303.0 K, including those corresponding to dissolution processes in neat water and DMSO. It is notable that all apparent standard dissolution thermodynamic quantities for neat water were taken from previous research results [32]. As expected, all the apparent standard Gibbs energies and apparent enthalpies of dissolution of meloxicam are positive in every case in these DMSO-aqueous systems. Otherwise, the apparent standard entropies of dissolution were negative in neat water, as well as in the mixtures of x 1 ≤ 0.30 but positive from the mixture of x 1 = 0.40 to neat DMSO. Thus, the global dissolution processes of meloxicam are always endothermic in nature and entropy-driven for those occurring in the composition interval of 0.40 ≤ x 1 ≤ 1.00; whereas, in neat water and the mixtures of x 1 ≤ 0.30, neither entropy or enthalpy-driving are observed because ∆ soln H • > 0 and ∆ soln S • < 0). All ∆ soln G • values decrease continuously from neat water (where highest ∆ soln G • value is obtained) to reach the lowest value in neat DMSO. Otherwise, ∆ soln H • decreases from neat water to reach the lowest value in the mixture of x 1 = 0.10, and later, it increases with the DMSO proportion to reach a quasi-plateau in the mixtures of  Finally, the apparent standard changes in entropy, for all the studied meloxicam d solution processes (∆solnS°) were calculated based on the respective ∆solnH° and ∆solnG° v ues by using [59,60]:     Moreover, to calculate the magnitude contributions by enthalpy (ζ H ) and entropy (ζ TS ) toward the dissolution processes, the following equations were used [64]: As observed in Table 5, the higher contribution to the positive apparent standard molar Gibbs energies of meloxicam dissolution is given by the positive enthalpy. This demonstrates that in almost all the mixtures, the main contributor to this positive standard molar Gibbs energy of solution of meloxicam (reflected in the low meloxicam solubility) is the enthalpy except for neat water and the mixture of x 1 = 0.10, where ζ H = 0.224 and 0.177, respectively, and thus, entropy is the dominant function in these two cases.

Apparent Thermodynamic Quantities of Mixing
The overall dissolution processes of meloxicam in all {DMSO (1) + water (2)} solvent systems may be represented with the following hypothetical stages: Solute (Solid state) at T hm → Solute (Solid state) at T fus → Solute (Liquid state) at T fus → Solute (Liquid state) at T hm → Solute (Solution state) at T hm Here, the hypothetical stages are considered as follows: (i) the heating and fusion of meloxicam at T fus = 536.7 K, (ii) the cooling of the liquid fused meloxicam to the considered temperature (i.e., T hm = 303.0 K), and (iii) the subsequent mixing of both the hypothetical super-cooled liquid meloxicam and the liquid aqueous-DMSO solvent system at T hm = 303.0 K [65]. This treatment allowed us to calculate every individual thermodynamic contribution toward the overall meloxicam dissolution processes by means of: where ∆ fus H T hm and ∆ fus S T hm represent the thermodynamic quantities relative to meloxicam melting and its cooling at T hm = 303.0 K, which, in turn, are calculated by means of [66]: ∆ fus S T hm = ∆ fus S T fus − ∆C p ln T fus T hm (21) where ∆C p denotes the difference of heat capacities of liquid and solid states at the temperature of melting. Owing the difficulties in ∆C p experimental determinations, the entropy of fusion (∆S f ) is used instead [66]. Table 6  In the mixtures 0.30 ≤ x 1 ≤ 0.50, neither enthalpy nor entropydriving is observed for mixing. Finally, in the interval of 0.60 ≤ x 1 ≤ 1.00, entropy-driven is observed. Furthermore, in order to evaluate the relative contributions by enthalpy (ζ H ) and entropy (ζ TS ) to the mixing processes in all these solvent systems, two equations analogous to Equations (16) and (17) were also employed. As observed, in water-rich and DMSO-rich mixtures, the main contributor to Gibbs energies of mixing is the entropy, but in the mixtures' composition interval of 0.30 ≤ x 1 ≤ 0.70, it is the enthalpy. As described earlier in the literature, the net variation of ∆ mix H • values regarding the aqueous-cosolvent mixtures' composition depends on the contribution of different kinds of intermolecular interactions. Hence, the cavity formation in the solvent system, required for the solute accommodation, is endothermic because some quantity of energy must be supplied against the respective cohesive forces of the solvent. This contribution diminishes the drug solubility as mentioned above. Oppositely, the solvent-solute interactions, resulting mainly from van der Waals and Lewis acid-base interactions, such as hydrogen bonding, are clearly exothermic in nature. This effect increases the drug solubility and dissolution rate as also indicated before. Even more, the structuring of water molecules around the phenylene ring and the methyl group of meloxicam structure (Figure 1) would be contributing to diminish the net ∆ mix H • quantity to small or even negative values in water-rich mixtures [67]. This event is clearly observed with meloxicam in aqueous-DMSO mixtures as indicated in Table 5 for solvent systems from neat water to the mixture of x 1 = 0.50. Table 6. Apparent thermodynamic functions relative to mixing processes of meloxicam (3) in {DMSO (1) + water (2)} mixtures at T hm = 303.0 K and p = 96 kPa a,b .

Enthalpy-Entropy Compensation Analysis
Classical extra-thermodynamic studies, in particular those based on the enthalpyentropy compensation analyses, provide a powerful physicochemical tool to find or identify similar mechanisms responsible for some physical and chemical processes involving organic compounds [68,69]. Some well-known literature reports demonstrated the presence of nonlinear-enthalpy-entropy compensation effects in the dissolution processes of many drugs and drug-alike compounds in different aqueous cosolvent binary systems. These extra-thermodynamic studies have usually been performed by different research groups to identify the main mechanisms involved in the cosolvent action for solubility increasing or decreasing, depending on the mixtures' composition [70][71][72]. As shown in Figure 6, meloxicam exhibits a nonlinear ∆ soln H • vs. ∆ soln G • trend with negative slopes from neat water to the mixture of x 1 = 0.10 and from the mixture of x 1 = 0.60 to neat DMSO, whereas, in the interval of 0.10 ≤ x 1 ≤ 0.30 and from x 1 = 0.50 to x 1 = 0.60, positive slopes are observed; finally, in the interval of 0.30 ≤ x 1 ≤ 0.50, a plateau is observed. In the first cases, the driving mechanism for transferring meloxicam from the most polar solvent systems to less polar solvent systems is the entropy. For the intervals exhibiting positive slopes, the drug transfer is driven by the enthalpy. Otherwise, in the interval of 0.30 ≤ x 1 ≤ 0.50, the function driving the transfer is not clear. Nevertheless, it is not easy to identify the molecular effects involved, owing the complexity of aqueous-DMSO mixtures as indicated earlier [41,42].
observed; finally, in the interval of 0.30 ≤ x1 ≤ 0.50, a plateau is observed. In the first cases, the driving mechanism for transferring meloxicam from the most polar solvent systems to less polar solvent systems is the entropy. For the intervals exhibiting positive slopes, the drug transfer is driven by the enthalpy. Otherwise, in the interval of 0.30 ≤ x1 ≤ 0.50, the function driving the transfer is not clear. Nevertheless, it is not easy to identify the molecular effects involved, owing the complexity of aqueous-DMSO mixtures as indicated earlier [41,42].

Preferential Solvation Analysis
The preferential solvation parameter of meloxicam (component 3) by DMSO (component 1) in the {DMSO (1) + water (2)} mixtures at saturation is defined as: where x L 1,3 is the local mole fraction of DMSO in the molecular environment around meloxicam and x 1 is the bulk mole fraction of DMSO in the initial aqueous-DMSO mixture in the absence of meloxicam. If δx 1,3 values were positive, meloxicam would be preferentially solvated by DMSO, but if they were negative, meloxicam would be preferentially solvated by water. Thus, the respective δx 1,3 values were obtained by means of the inverse Kirkwood-Buff integrals (IKBI) for the solvent components based on [73][74][75]: with, Here, κ T represents the isothermal compressibility of every {DMSO (1) + water (2)} mixture. V 1 and V 2 denote the partial molar volumes of DMSO and water in the aqueous-DMSO mixtures. V 3 denotes the partial molar volume of meloxicam. The function D corresponds to the first derivative of the variation of standard molar Gibbs energies of transfer of meloxicam from neat water to {DMSO (1) + water (2)} mixtures regarding the DMSO-proportion in the mixtures free of solute, as shown in Equation (27). The function Q involves the second derivative of the variation of excess molar Gibbs energy of mixing of DMSO and water (G Exc 1+2 ) regarding the water-proportion in the aqueous-DMSO mixtures, as shown in Equation (28). V cor is the correlation volume and r 3 is the hydrodynamic molecular radius of meloxicam, which is commonly calculated by means of Equation (29), where N Av is the number of Avogadro.
As exposed in the literature, the definitive V cor values require iteration because they depend on the local mole fractions of DMSO and water around the meloxicam molecules in the equilibrated solutions. Hence, these iterations are performed by substituting δx 1,3 and V cor values in Equations (22), (23), and (26) in order to recalculate the x L 1,3 value until almost invariant values of V cor are obtained. Figure 7 depicts the apparent Gibbs energies of transfer of meloxicam from neat water to {DMSO (1) + water (2) ∆ tr G o 3,2→1+2 values were correlated according to the regular fourth degree polynomial presented as Equation (31), with adjusted r 2 = 0.998, typical error = 0.317, and F value = 1484. In this way, the D values shown in Table 7 were calculated from the first derivative of Equation (31)    In this way, the D values shown in Table 7 were calculated from the first derivative of Equation (31) by considering the variation of aqueous-DMSO mixtures composition in incremental x 1 = 0.05 steps through all the mixtures' composition interval. Otherwise, the required Q, RTκ T , V 1 , and V 2 values corresponding to {DMSO (1) + water (2)} mixtures were taken from the literature [76]. Because V 3 values are not available for meloxicam in {DMSO (1) + water (2)} mixtures, these values were considered as the one calculated based on Fedors' method, i.e., 183.3 cm 3 ·mol −1 [6]. G 1,3 and G 2,3 shown in Table 7 are negative in all cases, indicating the affinity of meloxicam by DMSO and water. The approximated hydrodynamic radius of meloxicam (r 3 ) was calculated as 0.417 nm by means of Equation (29). In turn, the preferential solvation parameters of meloxicam by DMSO molecules are also summarized in Table 7. According to Figure 8, initially adding of DMSO to water makes the δx 1,3 values of meloxicam negative in the interval from neat water to the mixture of x 1 = 0.21. The maximum negative value of this parameter is obtained in the mixture of x 1 = 0.10, with δx 1,3 = −1.36 × 10 −2 , which is higher than 1.00 × 10 −2 if the absolute value is considered; therefore, it could be a consequence of real preferential solvation effects by water on meloxicam, rather than a consequence of the uncertainty propagation in the respective IKBI calculations [77,78]. The cosolvent action of DMSO for increasing the meloxicam solubility in these water-rich mixtures could be associated to the breaking of the ordered structure exhibited by water molecules, such as "icebergs", around the non-polar moieties of meloxicam, which, in turn, would be increasing the meloxicam solubility and solvation.
In mixtures of 0.21 < x 1 < 1.00, the δx 1,3 values are positive indicating preferential solvation of meloxicam by DMSO. Maximum δx 1,3 value was obtained in the mixture of x 1 = 0.45 (δx 1,3 = 3.34 × 10 −2 ). This maximum positive δx 1,3 value is also higher than |1.00 × 10 −2 | being a consequence of real preferential solvation effects by DMSO [77,78]. From a mechanistic viewpoint, in the mixtures' composition region of 0.19 < x 1 < 1.00, it is adequately conjecturable that meloxicam could be acting as a Lewis acid in front of the DMSO molecules owing the unshared electrons of the sulfoxide oxygen atom of this cosolvent. Notably, this cosolvent is more basic than water, as remarkable by the magnitude of their Kamlet-Taft hydrogen bond acceptor parameters, namely β = 0.76 for DMSO and 0.47 for water [38]. Moreover, Figure 8 allows the comparison of preferential solvation of meloxicam by DMSO and N,N-dimethylformamide in their respective aqueous mixtures [6]. As observed, the cosolvent regions of preferential solvation are similar, as well as the magnitudes of preferential solvation by both cosolvents. Nevertheless, preferential hydration of meloxicam is higher with N,N-dimethylformamide, which could be a consequence of its lower polarity (δ 1 = 24.1 MPa 1/2 ) regarding DMSO (δ 1 = 26.6 MPa 1/2 ) [37,38]. These interesting behaviors could be a consequence of the high water-association effects around the non-polar groups of this drug, which could be favored by the more hydrophobic moieties present in the cosolvents, as they exhibit less polar nature. In turn, the hydrophobic groups present in the cosolvents could also be acting as water-association promotors depending on their respective molecular sizes [11]. Finally, from all the physicochemical analyses reported, it is noteworthy to mention that this investigation expands the equilibrium solubility database about non-steroidal anti-inflammatory drugs in commonly used aqueous cosolvent mixtures [79].
Liquids 2022, 2, FOR PEER REVIEW Because 3 V values are not available for meloxicam in {DMSO (1) + water (2)} m tures, these values were considered as the one calculated based on Fedors' method, 183.3 cm 3 ·mol −1 [6]. G1,3 and G2,3 shown in Table 7 are negative in all cases, indicating affinity of meloxicam by DMSO and water. The approximated hydrodynamic radiu meloxicam (r3) was calculated as 0.417 nm by means of Equation (29). In turn, the pre ential solvation parameters of meloxicam by DMSO molecules are also summarized Table 7. According to Figure 8, initially adding of DMSO to water makes the δx1,3 val of meloxicam negative in the interval from neat water to the mixture of x1 = 0.21. maximum negative value of this parameter is obtained in the mixture of x1 = 0.10, w δx1,3 = −1.36 × 10 −2 , which is higher than 1.00 × 10 −2 if the absolute value is conside therefore, it could be a consequence of real preferential solvation effects by water meloxicam, rather than a consequence of the uncertainty propagation in the respec IKBI calculations [77,78]. The cosolvent action of DMSO for increasing the meloxicam ubility in these water-rich mixtures could be associated to the breaking of the orde structure exhibited by water molecules, such as "icebergs", around the non-polar moie of meloxicam, which, in turn, would be increasing the meloxicam solubility and solvat In mixtures of 0.21 < x1 < 1.00, the δx1,3 values are positive indicating preferential s ation of meloxicam by DMSO. Maximum δx1,3 value was obtained in the mixture of 0.45 (δx1,3 = 3.34 × 10 −2 ). This maximum positive δx1,3 value is also higher than |1.00 × 1 being a consequence of real preferential solvation effects by DMSO [77,78]. From a me anistic viewpoint, in the mixtures' composition region of 0.19 < x1 < 1.00, it is adequa conjecturable that meloxicam could be acting as a Lewis acid in front of the DMSO m cules owing the unshared electrons of the sulfoxide oxygen atom of this cosolvent. N bly, this cosolvent is more basic than water, as remarkable by the magnitude of their K let-Taft hydrogen bond acceptor parameters, namely β = 0.76 for DMSO and 0.47 for w [38]. Moreover, Figure 8 allows the comparison of preferential solvation of meloxicam DMSO and N,N-dimethylformamide in their respective aqueous mixtures [6]. As served, the cosolvent regions of preferential solvation are similar, as well as the mag tudes of preferential solvation by both cosolvents. Nevertheless, preferential hydratio meloxicam is higher with N,N-dimethylformamide, which could be a consequence o lower polarity (δ1 = 24.1 MPa 1/2 ) regarding DMSO (δ1 = 26.6 MPa 1/2 ) [37,38]. These inter ing behaviors could be a consequence of the high water-association effects around non-polar groups of this drug, which could be favored by the more hydrophobic moie present in the cosolvents, as they exhibit less polar nature. In turn, the hydropho

Conclusions
Equilibrium molar and mole fraction solubilities of analgesic drug meloxicam in different {DMSO (1) + water (2)} mixtures at five temperatures from 293.15 to 313.15 K were determined by using the shake flask method followed by UV-vis drug quantification, reported and analyzed. Meloxicam mole fraction solubility in these mixtures was adequately correlated with some well-known correlation models obtaining mean percentage deviations (MPDs) of 9.6 to 9.9%. In addition, a number of predictive models which were already trained using published datasets and by employing the minimum number of measured experimental data from this project were produced with MPDs of 35.6 to 74.3%. Apparent standard thermodynamic quantities of dissolution and mixing processes were calculated observing endothermal dissolution processes in all cases and favored in DMSO-rich mixtures. Nonlinear enthalpy-entropy compensation was observed indicating different mechanisms for the cosolvent action. IKBI treatment demonstrated preferential hydration of meloxicam in water-rich mixtures but preferential solvation by DMSO in mixtures of 0.21 < x 1 < 1.00.