Experimental measurement, thermodynamic analysis, and mathematical modeling for budesonide solubility in 1-propanol + water mixtures at T = (293.2 to 313.2) K

Budesonide (BDS) a steroid-based anti-inflammatory drug widely prescribed for various diseases, has a low aqueous solubility. In this study, we investigated cosolvency approach to study the thermodynamic specifications related to the solubility of BDS at the temperature range of 293.2–313.2 K in (1-propanol + water) mixtures applying the shaking flask method. The predictive power of different mathematical models for experimental data in the cosolvency systems was evaluated. For this purpose, the linear and nonlinear mathematical equations such as van’t Hoff model (as a linear model), Buchowski-Ksiazczak equation (as a non-linear), CNIBS/R–K and MRS models (as a linear model for solvent composition at an isothermal condition), modified Wilson model (as a non-linear model for isothermal condition), the Jouyban-Acree model (as a model that considers temperature and solvent composition), and Jouyban-Acree-van’t Hoff model (as a model with no further input data) were studied. Also, the Williams-Amidon excess Gibbs energy model was investigated. In addition, the related apparent thermodynamics of the BDS dissolution process in the desired temperature such as Gibbs free energy, enthalpy, and entropy, were computed by the corresponding equations. Moreover, based on the inverse Kirkwood-Buff integrals, it is demonstrated that BDS is preferentially solvated by water in water-rich mixtures. The accuracy of the fitness was evaluated with mean relative deviations (MRDs%) for back-calculated molar BDS solubility data. The result showed that the maximum solubility of BDS was obtained at 0.7 mass fraction of 1-propanol at all temperatures. Thermodynamic studies demonstrated that BDS dissolution procedures were obtained as endothermic and entropy-driven in almost all cases. The overall MRDs% values for the back-computed BDS solubility in the aqueous mixture of 1-propanol based on van’t Hoff model, Buchowski-Ksiazczak equation, CNIBS/R–K model, modified Wilson model, Jouyban-Acree model, Jouyban-Acree-van’t Hoff model, MRS model, and Williams-Amidon excess Gibbs energy model were found 1.93%, 1.80%, 11.68%, 33.32%, 12.30%, 9.24%, 10.70%, and 6.57%, respectively.

Graphical Abstract