key: cord-0991069-rhexh20f authors: Li, Xinbao; Cong, Yang; Li, Wentian; Yan, Pengyao; Zhao, Hongkun title: Thermodynamic modelling of solubility and preferential solvation for ribavirin (II) in co-solvent mixtures of (methanol, n-propanol, acetonitrile or 1,4-dioxane) + water date: 2017-07-24 journal: J Chem Thermodyn DOI: 10.1016/j.jct.2017.07.027 sha: 1bcfbd09c5338a3c2465f755505cdd93b1e6cbbb doc_id: 991069 cord_uid: rhexh20f The equilibrium solubility of ribavirin in solvent mixtures of {methanol (1) + water (2)}, {n-propanol (1) + water (2)}, {acetonitrile (1) + water (2)} and {1,4-dioxane (1) + water (2)} was determined experimentally by using isothermal dissolution equilibrium method within the temperature range from (278.15 to 318.15) K under atmospheric pressure (101.1 kPa). At the same temperature and mass fraction of methanol (n-propanol, acetonitrile or 1,4-dioxane), the mole fraction solubility of ribavirin is greater in (methanol + water) than in the other three solvent mixtures. The preferential solvation parameters were derived from their thermodynamic solution properties by means of the inverse Kirkwood–Buff integrals. The preferential solvation parameters for methanol, n-propanol, acetonitrile or 1,4-dioxane (δx(1,3)) were negative in the four solvent mixtures with a very wide compositions, which indicated that ribavirin was preferentially solvated by water. Temperature had little effect on the preferential solvation magnitudes. The higher solvation by water could be explained in terms of the higher acidic behaviour of water interacting with the Lewis basic groups of the ribavirin. Besides, the solubility of the drugs was mathematically represented by using the Jouyban-Acree model, van’t Hoff-Jouyban-Acree model and Apelblat-Jouyban-Acree model obtaining average relative deviations lower than 1.57% for correlative studies. It is noteworthy that the solubility data presented in this work contribute to expansion of the physicochemical information about the solubility of drugs in binary solvent mixtures and also allows the thermodynamic analysis of the respective dissolution and specific solvation process. The drugs solubility in solvent mixtures as a function of temperature and composition is evaluated essentially for the aims of raw material purification, design of liquid dosage forms, and understanding of the mechanisms regarding to the physical and chemical stability of pharmaceutical dissolutions [1, 2] . Therefore, the solubility of active ingredients is a significant physicochemical property to be considered in pharmaceutical design because it affects the drug efficacy, influencing several biopharmaceutical and pharmacokinetic properties [3, 4] . Moreover, temperaturedependence of the solubility allows evaluating the preferential solvation of the solute by the solvent components in mixtures, which gives a powerful tool in the understanding of molecular interactions in relation to the dissolution processes of drug [5] [6] [7] . Ribavirin (II) (CAS Registry No. 36791-04-5, chemical structure shown in Fig. S1 of Supporting material) is a well-known antiviral drug with a broad spectrum of action. Because of its activity against both RNA and DNA viruses [8] , it is recommended for the treatment of some infectious disorders including herpes [9] , respiratory syncytial viral disease, and acute respiratory infections caused by adenovirus [10] . At present, ribavirin is also part of a combination therapy for the treatment of hepatitis C [11] [12] [13] [14] [15] [16] and Lassa fever [14, 15] . Furthermore, it has been evaluated for severe acute respiratory syndrome (SARS) and as an anti-cancer drug with successful results for the treatment of leukemia and breast cancer [16] [17] [18] . Due to its good properties, ribavirin has been received great attention. Nevertheless, in spite of the usefulness of this drug, the physicochemical properties of ribavirin in aqueous and organic solutions have not been studied systematically in the previous works [19] [20] [21] [22] [23] [24] . Although several theoretical and semi-empirical models have employed to predict drug solubility in solvent mixtures, the http://dx.doi.org/10.1016/j.jct.2017.07.027 0021-9614/Ó 2017 Elsevier Ltd. availability of experimental data is still fundamental for the pharmaceutical scientists [1, 25] . It is well-known that solvent mixtures has been widely used in pharmacies long ago [1] [2] [3] [4] , only recently the mechanisms involving in the increase or decrease in drugs' solubility start to be approached from a deep thermodynamic point of view, including the analysis of the preferential solvation of solute by the components of mixtures [5] [6] [7] [26] [27] [28] [29] [30] . In terms of the considerations mentioned above, the main goal of this work is to report the equilibrium solubility of ribavirin in {methanol (1) + water (2)}, {n-propanol (1) + water (2)}, {acetonitrile (1) + water (2)}+ and {1,4-dioxane (1) + water (2)} mixtures at different temperatures in order to evaluate the respective thermodynamic quantities of the solution, as well as the preferential solvation of the drug by these organic solvents. This research expands the available solubility data about drug in neat organic solvents and solvent mixtures [1] and also allows the thermodynamic analysis of the respective dissolution and specific solvation process. It is noteworthy that methanol and 1,4-dioxane is not used to develop liquid medicines due to its high toxicity. But in some instances methanol is used in drug purification procedures [31] , as well as solvent in some drug microencapsulation techniques [32] . Moreover, methanol is widely used as mobile phase in high performance liquid chromatography [33] . In addition, 1,4-dioxane is widely used as a model co-solvent because it is miscible with water in all the range of compositions although it exhibits a low polarity as described by its dielectric constant [34] . Even more, the Jouyban-Acree model has been used to correlate the solubility of a lot of drugs in {methanol/1,4-dioxane (1) + water (2)} mixtures [35] . n-Propanol and acetonitrile are green solvents [36] and miscible with water in all compositions. In addition, although n-propanol is not widely used as co-solvent for design of liquid medicines, it has been used as solvent in the pharmaceutical industry for resins and cellulose esters [37, 38] . Acetonitrile is most widely utilized because of the high resolution attainable probably due to the low viscosity of water-acetonitrile mixtures [39] . Moreover, the wateracetonitrile mixtures may have a partially hydrophobic character similar to the methanol system since only the OH group in methanol is replaced by the CN group in acetonitrile. Ribavirin (II) was provided by Shanghai Sharing Technologies Co., Ltd., China with a mass fraction of 0.983. It was purified three times in methanol. The final composition of ribavirin (II) used for solubility measurement was 0.996 in mass fraction, which was analysed by using a high-performance liquid chromatography (HPLC, Agilent 1260). The four neat solvents (methanol, npropanol, acetonitrile and 1,4-dioxane) were provided by Sinopharm Chemical Reagent Co., Ltd., China, the mass fraction purities of which were all no less than 0.994 determined by gas chromatography {Smart (GC-2018)}. The water was twice distilled water prepared in our lab, which conductivity was less than 2 lSÁcm À1 . The detailed information of these chemicals was presented in Table S1 of Supporting material. The diagram of solubility determination apparatus was given in Fig. S2 of Supporting material. The experimental apparatus consisted of a 100 ml jacketed glass vessel with a circulating (water + isopropanol) system and a magnetic stirrer. The temperature of circulating (water + isopropanol) mixture was kept by a thermostatic bath (Model: QYHX-1030), which was provided by Shanghai Joyn Electronic CO., Ltd., China and had a standard uncertainty of 0.05 K. The actual temperature of solution was displayed by a mercury glass micro thermometer (standard uncertainty: 0.02 K) inserted in the inner chamber of the jacket glass vessel. A condenser was employed to avoid escaping of the solvent. An analytical balance (model: BSA224S) having a standard uncertainty of 0.0001 g and provided by Satorius Scientific Instrument (Beijing) was used to determine the mass of the solute, solvent, and saturated solution. The reliability of apparatus was verified in advance by determining the benzoic acid solubility in toluene [40, 41] . During the experiment, we use the analytical balance (model: BSA224S) to prepare the solvent mixtures. Each mixed solvent in the glass vessel was about 60 mL, and the relative standard uncertainty of which was estimated to be 0.0002. The mass fractions of organic solvent in the binary mixtures varied from 0 to 1. The glass vessel was covered with a stopper to prevent the solvent from escaping during the preparation procedure of solvent mixtures. During the experiment, the atmospheric pressure was about 101.1 kPa. In this work the solid-liquid equilibrium for ribavirin in binary solvent mixtures of (methanol + water), (n-propanol + water), (acetonitrile + water) and (1,4-dioxane + water) was determined by an isothermal dissolution equilibrium method [40] [41] [42] [43] [44] [45] , and the highperformance liquid phase chromatograph was employed to determine the solubility of ribavirin in equilibrium liquid phase. For each experiment, saturated solutions of ribavirin (II) were prepared in the jacketed glass vessel. An excessive amount of ribavirin (II) was introduced into the jacketed glass vessel filled with about 60 ml solvent mixtures. Continuous stirring was achieved by using a magnetic stirrer at a given temperature to mix the solution intensively. The temperature of the solutions was kept at a desired value by circulating (water + isopropanol) mixture from the thermostatic bath through the outer jacket. So as to obtain the equilibration time of the studied systems, about 0.5 ml liquid phase was withdrawn every one hour by using a 2 mL of preheated syringe equipped with a pore syringe filter (PTFE 0.2 lm), and then analysed by the high-performance liquid phase chromatograph (Agilent-1260). Once the composition of liquid phase didn't vary, the system was believed to reach equilibrium. In order to make sure that sampling was performed at equilibrium conditions, two types of experiments were carried out, one starting from a supersaturated solution, in which the solid phase precipitated to reach equilibrium and the other starting from a non-saturated solution, in which solid dissolved to reach equilibrium. The results revealed that it took about 13 h to reach equilibrium for all the studied systems. Once the mixture arrived at equilibrium, stirring was turned off to allow any undissolved solute to be precipitated. One hour later, the upper liquid phase was withdrawn using the 2 ml of preheated or precooled syringe, and transferred rapidly to a preweighed volumetric flask equipped with a rubber stopper. The volumetric flask filled with sample was weighed again with the analytical balance. Subsequently, the sample was diluted with methanol, and 1 ll of the solution was withdrawn for analysis by using the high-performance liquid chromatography.The equilibrium mole fraction solubility of ribavirin (II) ðx w;T Þ in the four binary solvent mixtures are obtained with Eq. (1), and the initial composition of binary solvent mixtures (w) is calculated with Eqs. (2) and (3). where m 1 is the mass of ribavirin; m 2 is the mass of methanol, npropanol, acetonitrile and 1,4-dioxane, and m 3 is the mass of water. M 1 , M 2 and M 3 are the corresponding molar mass. The concentration of ribavirin was analysed by the Agilent-1260 high-performance liquid chromatography (HPLC). The chromatographic column was a reverse phase column with a type of LP-C18 (250 mm  4.6 mm), which temperature was kept at 303 K. The wavelength of the UV detector was set to 210 nm, which was determined by continuous UV-Scanning. The mobile phase comprised two component, which were methanol and water with a flow rate of 0.8 mlÁmin À1 . The volume ratio of methanol to water was 25:75. Each analysis was carried out three times, and the average value of three measurements was regarded as the final value of the analysis. The relative standard uncertainty is estimated to be 4.5% in mole fraction solubility. The crystal of ribavirin was identified by X-ray powder diffraction (XPRD) carried out on a HaoYuan DX-2700B (HaoYuan, China) instrument. The samples were determined by Cu Ka radiation (k = 1.54184 nm), and the tube voltage and current were set at 40 kV and 30 mA, respectively. The data were collected at room temperature from 5°to 80°(2-Theta) at a scan speed of 6 degÁmin À1 under atmospheric pressure. Ribavirin has long been known to exist in two enantiotropically related polymorphic forms (I and II) [46] [47] [48] . Form II is the stable form at ambient temperature and melts at 439-441 K; while form I, the metastable polymorph, melts at (447-449) K [47, 48] . In order to demonstrate the existence of polymorph transformation or solvate formation of ribavirin during the solubility determination, the equilibrium solid phase is collected and analysed by XRD. The patterns of the raw material and the solids crystallized in neat solvents and solvent mixtures are plotted in Fig. S3 of Supporting material. It is confirmed by XRD pattern that all the XRD patterns of solid phase of ribavirin in equilibrium with its solution have the same characteristic peaks with the raw material, which corresponds to form II reported in literatures [46] [47] [48] . Therefore, it can be concluded that there was no polymorph transformation or solvate formation during the entire experiment. The determined mole fraction solubility of ribavirin in binary solvent mixtures of (methanol + water), (n-propanol + water), (acetonitrile + water) and (1,4-dioxane + water) are presented in Tables 1-4, respectively. In addition, the relationship between the mole fraction solubility and temperature and solvent composition are shown graphically in Figs. 1-4. It can be seen that, for the studied solvent mixtures, the ribavirin solubility is a function of temperature and solvent composition. The solubility of ribavirin increase with increasing temperature and mass fraction of water for the binary systems studied. The maximum solubility of ribavirin is observed in water. It can also be seen from Tables 1-4 that, at the same temperature and solvent composition, the solubility of ribavirin in (methanol + water) is greater than those in (n-propanol + water), (acetonitrile + water) and (1,4-dioxane + water). This behaviour may be due to the relative strong polarity of ribavirin molecule. At the same composition, the polarities of the four mixed solvents rank as (methanol + water) > (n-propanol + water) > (acetonitrile + water) > (1,4-dioxane + water) [49] . Therefore, the solubilities of ribavirin in (methanol + water) are larger than those in the other solvent mixtures. It is noteworthy that the mole-fraction solubility of ribavirin (II) in water at 298.15 K determined in this work is 0.01741, which is a little larger than 0.0103 (142 mgÁmLl À1 ) reported by Benet [19] , Goodarzi [20] , Dave [21] and Chen [22] , 0.0109 (149 mgÁmL À1 ) reported by Vasa [23] , and 0.0124 (169.99 mgÁmL À1 ) reported by Prasanthi [24] . The difference may be due to many factors, such as determination method, equilibration time, purity, sampling, analysis method and so on. Based on study made by Vasa and Wildfong [23] , ribavirin (II) reached a maximum solubility within 30 min. On the other hand, it can also be seen from the plots of the ribavirin (II) concentrations in distilled water as a function of time that, when the dissolution time is 180 min, the concentration of ribavirin (II) in water remains increased obviously. Therefore, it can be speculated that if the time is prolonged, the ribavirin (II) concentration will also increase. In other words, the concentration (149 mgÁml À1 ) at 30 min is not the equilibrium solubility of ribavirin (II) in water at studied temperature. The models presented from 1960 to 2007 for correlating the solubility of a solid in solvent mixtures have been reviewed in Ref. [25] . In the present work, three models are used to correlate the solubility of ribavirin in binary solvent mixtures of (methanol + water), (n-propanol + water), (acetonitrile + water) and (1,4dioxane + water) at different temperatures, which correspond to JouybanÀAcree model [25, 41, 50] , a combination of the Jouy-banÀAcree model with van't Hoff equation [41, 50] , and a combination of the JouybanÀAcree model with modified Apelblat equation [41, 50] . The Jouyban-Acree model is expressed as Eq. (4). This model can provide accurate mathematical description for the dependence of solute solubility on both temperature and solvent composition for binary and ternary mixed solvents [25, 41, 50] . where x w,T denotes the mole fraction solubility of solute in solvent mixtures at temperature T/K; w 1 and w 2 are the mass fraction of solvents 1 (methanol, n-propanol, acetonitrile or 1,4-dioxane) and 2 (water) in the absence of the solute (ribavirin), respectively; x 1,T and x 2,T are the solute solubility in mole fraction in neat solvent; and J i are the Jouyban-Acree model parameters. The van't Hoff equation is an ideal model, which is described as Combining Eq. (4) and Eq. (5), the van't Hoff-Jouyban-Acree model can be obtained [41, 50] and expressed as Eq. (6). A 1 , B 1 and A 2 , B 2 are equation parameters. The modified Apelblat equation is a semi-empirical model having three parameters. It may be employed to describe a nonlinear relationship between solubility ðln x T Þ in neat solvent and reciprocal of absolute temperature (1/T) and is described as where A, B, and C are equation parameters; and also x T is the mole fraction solubility of ribavirin in studied mixed solvents at absolute temperature T in Kelvin. By substituting Eq. (7) into Eq. (4), the modified Apelblat-Jouyban-Acree model can be obtained as Eq. (8) [41, 50] . The experimental solubility values of ribavirin in (methanol + water), (n-propanol + water), (acetonitrile + water) and (1,4dioxane + water) mixtures are correlated and calculated with Eqs. (4), (6) and (8) . The objective function is described as In order to show the error and evaluate the different models, the relative average deviation (RAD) and root-mean-square devia- 10) and (11) . where N is the number of experimental data points. x e w;T denotes the mole fraction solubility determined in this work; and x c w;T , the mole fraction solubility calculated with the corresponding solubility model. The parameters in Eqs. (4) to (8) were acquired by using Mathcad software based on the experimental solubility values. The obtained values of model parameters together with the RAD and the RMSD values are presented in Table S2 of the Supporting material. The solubility of ribavirin in the four binary mixtures of (methanol + water), (n-propanol + water), (acetonitrile + water) and (1,4-dioxane + water) are evaluated on the basis of the regressed parameters' values. The calculated solubilities using the JouybanÀAcree model are plotted in Figs. 1-4. It can be seen from Table S2 that for the studied binary solvent mixtures, the maximum value of relative average deviations (RAD) is 1.57%, which is obtained with the van't Hoff-Jouyban-Acree model for the system of (n-propanol + water). Besides, the root-meansquare deviations (RMSD) are no greater than 2.14  10 À4 . Comparison with the other models, the values of RAD and RMSD obtained with the Jouyban-Acree model is relatively small. On the whole, the three models can all be employed to correlate the solubility of ribavirin in the binary mixtures of (methanol + water), (npropanol + water), (acetonitrile + water) and (1,4-dioxane + water) at all initial composition ranges, and the Jouyban-Acree model provides better correlation results than the other two models. Thermodynamic property of a solute dissolved in mixed solvent can provide significant information for the dissolution process. The standard dissolution enthalpy ðDH o sol Þ for dissolution process of ribavirin in the studied solvent mixtures may be attained from the famous van't Hoff analysis [38, 41] . here R is the universal gas constant (8.314 JÁK À1 Ámol À1 ). T hm stands for the mean harmonic temperature which is calculated via Eq. (13). The values of DH o sol can be obtained from the slope of the curves of lnx w,T vs (1/T À 1/T hm ), which are shown graphically in Figs. S4-S7 of Supporting material. The standard dissolution enthalpy of ribavirin in the solvent mixtures of (methanol + water), (npropanol + water), (acetonitrile + water) and (1,4-dioxane + water) are computed from Eq. (12) and tabulated in Table S3 of Supporting material. It can be found that the values of standard molar enthalpy of dissolution are all positive, which demonstrates that the dissolution process of ribavirin in the four binary mixed solvents is endothermic, and the entropy is driving force for the dissolution process. Inverse Kirkwood-Buff integrals (IKBI) are widely used for evaluating the preferential solvation of non-electrolyte or nondissociated weak electrolyte compounds in solvent mixtures, describing the local solvent proportions around the solute with respect to the composition of the solvent mixtures [6, 7] . This treatment depends on the standard molar Gibbs energies of transfer of ribavirin (3) from neat water (2) to the solvent mixtures and also on the excess molar Gibbs energy of mixing for the solvent binary mixtures free of drug. Thus, the results are expressed in terms of the variation of preferential solvation parameter (dx 1,3 ) of the solute ribavirin (3) by the solvent molecules with the mixtures composition. The inverse Kirkwood-Buff integral equation is given by the following general expression: where, g i,3 is the pair correlation function for molecules of solvent i in the solvent mixtures around the solute ribavirin (3); r is the distance between the centers of molecules of ribavirin (3) and those of solvent (1) or water (2); and r cor is a correlation distance for which g i,3 (r > r cor ) % 1. Accordingly, for all distances r > r cor up to infinite, the value of the integral is essentially zero. The preferential solvation parameter of ribavirin (compound 3) by the solvent (compound 1) in solvent (1) + water (2) mixtures is defined as [5] [6] [7] [26] [27] [28] [29] [30] : where x L 1;3 is the local mole fraction of solvent (1) in the environment near to ribavirin (3) and x 1 is the bulk mole fraction composition of solvent (1) in the initial binary solvent. If dx 1,3 > 0 then the solute is preferentially solvated by solvent (1); on the contrary, if dx 1,3 is <0 the solute is preferentially solvated by water (2) . Values of dx 1,3 may be obtained from the inverse Kirkwood-Buff integrals for the individual solvent components analysed according to some thermodynamic quantities as shown in the following equations [5] [6] [7] [26] [27] [28] [29] [30] : with, In Eqs. (16) (19) , j T is the isothermal compressibility of the solvent (1) + water (2) mixtures (expressed in GPa À1 ); V 1 and V 2 are the partial molar volumes of the solvents in the mixtures (expressed in cm 3 Ámol À1 ); similarly, V 3 is the partial molar volume of ribavirin in these mixtures (expressed in cm 3 Ámol À1 ). The function D express by Eq. (20) is the derivative of the standard molar Gibbs energies of transfer of ribavirin from neat water (2) to solvent (1) + water (2) mixtures with respect to the solvent composition (expressed in kJÁmol À1 , as is RT). The function Q defined by Eq. (21) involves the second derivative of the excess molar Gibbs energy of mixing of the two solvents ðG Exc 1þ2 Þ with respect to the water proportion in the mixtures (also expressed in kJÁmol À1 ). V cor is the correlation volume and r 3 is the molecular radius of ribavirin calculated by using Eq. (22) with N Av as the Avogadro's number. Because of the dependence of j T on composition, this term is not known for all the systems investigated. In addition, due to the minor contribution of RTj T to the inverse Kirkwood-Buff integral, the dependence of j T upon composition can be calculated approximated as an additive property by using the mixture compositions and the reported values for neat solvents by [5, 6, [26] [27] [28] [29] [30] : D tr G 0 3;2!1þ2 ¼ RT ln The calculated values of D tr G 0 3;2!1þ2 are tabulated in Tables S4 and S5 Table S6 of Supplementary material. In similar manner, the partial molar volumes of both solvents in the mixtures are calculated from the reported density values of the solvent mixtures at the studied temperatures under study by Mikhail and Aliaj for methanol + water mixtures [53, 54] , by Pang for (n-propanol + water) mixtures [55] , by Saleh and Handa for {acetonitrile (1) + water (2)} mixtures [56, 57] , and by Ruidiaz for (1,4-dioxane + water) mixtures [58] using Eqs. (35) and (36) . In these equations V is the molar volume of the mixtures calculated as V =(x 1 ÁM 1 +x 2 ÁM 2 )/q. Here, M 1 is 32.04 gÁmol -1 for methanol, 60.06 gÁmol -1 for n-propanol, 41.05 gÁmol -1 for acetonitrile, 88.11 gÁmol -1 for 1,4-dioxane and M 2 is 18.02 gÁmol -1 for water. Because no partial molar volumes of ribavirin (3) in these mixtures are acquired in the literature, in this work this property is considered as similar to that for the pure compound as a good approximation [26] [27] [28] [29] [30] . In this way, the molar volume of ribavirin (3) is calculated by using Advanced Chemistry Development (ACD/ Labs) Software V 11.02 (sc1994-2016 ACD/Labs) as 117.0 cm 3 Ámol -1 . From this volume value, solute radius value (r 3 ) is calculated using Eq. (22) as 0.359 nm. The values of G 1,3 and G 2,3 can be obtained and shown in Tables S8-S11 of Supplementary material. Definitive correlation volume requires iteration because it depends on the local mole fractions around the solute. This iteration is done by replacing dx 1,3 and V cor in the Eqs. (15) , (16) and (19) to recalculate x L 1;3 until a non-variant value of V cor is acquired. The obtained values of correlation volume V cor and dx 1,3 are listed in Tables S12-S15 of Supplementary material for {methanol (1) + water (2) Maximum negative values are obtained with the composition x 1 = 0.75 with dx 1,3 = À2.314  10 À2 to À2.411  10 À2 for the methanol (1) + water, x 1 = 0.75 with dx 1,3 = À6.129  10 À2 to À6.488  10 À2 for the n-propanol (1) + water (2), x 1 = 0.80 with dx 1,3 = À61.20  10 À2 to À66.60  10 À2 for the {acetonitrile (1) + water (2)}, and x 1 = 0.60 to 0.65 with dx 1,3 = À78.98  10 À2 to À33.52  10 À2 for {1,4-dioxane (1) + water (2)} mixtures. In the other regions for the {methanol (1) + water (2)}, {n-propanol (1) + water (2)} and {1,4-dioxane (1) + water (2)} mixtures, the dx 1,3 values of ribavirin (3) are positive. However they are all less than 1.0  10 À2 . For the acetonitrile (1) + water (2) mixtures in water-rich composition, the dx 1,3 values are positive and slightly greater than 1.0  10 À2 (e.g. 1.47  10 À2 at 293.15 K). Nevertheless these values are much smaller than that in intermediate and acetonitrile-rich composition (e.g. -61.20  10 À2 at 293.15 K). Therefore they are just qualitative because of uncertainties propagation effects instead of preferential solvation [59] [60] [61] . So the conclusion can be made that ribavirin is preferentially solvated by water for the four solvent mixtures with a very wide compositions. Probably the structuring of water molecules around the ribavirin molecule contributes to lowering of the net dx 1,3 to negative values in the four solvent mixtures. It can also be seen from Fig. 5 that the temperature has little effect on the preferential solvation magnitudes of ribavirin in the four solvent mixtures. Ribavirin can act as a Lewis base because of the free electron pairs in nitrogen atoms of AN@ and AN and oxygen atoms of @O and AOA (Fig. S1 of Supporting materials), which interact with acidic hydrogen atoms of water. So it is conjecturable that ribavirin could be acting mainly as a Lewis base in front to water because the Kamlet-Taft hydrogen bond donor parameters are, a = 0.990 for methanol, 0.84 for n-pronanol, 0.19 for acetonitrile, 0.00 for 1,4-dioxane and 1.17 for water, respectively [51, 62] . The preferential solvation magnitudes of ribavirin by water is highest in acetonitrile mixtures, followed by 1,4-dioxane and n-propanol mixtures, and finally, by methanol mixtures in different solvent proportions, i.e. x 1 = 0.75, dx 1,3 = À2.367  10 À2 for {methanol (1) + water (2)}; x 1 = 0.75, dx 1,3 = À6.396  10 À2 for {nÀpropanol (1) + water (2)}; x 1 = 0.80, dx 1,3 = À62.50  10 À2 for {acetonitrile (1) + water (2)}; and x 1 = 0.60, dx 1,3 = À28.43  10 À2 for {1,4-dioxane (1) + water (2)}, respectively. Finally, it is noteworthy that in our previous studies [63] [64] [65] , the solute is preferentially solvated by co-solvent in co-solventrich mixtures. However, ribavirin (II) is preferentially solvated by water in this region. The reason is not clear at present. This case may be due to the chemical structure of solute. The equilibrium solubility of ribavirin in four binary mixed solvents of {methanol (1) + water (2)}, {n-propanol (1) + water (2)}, {acetonitrile (1) + water (2)} and {1,4-dioxane (1) + water (2)} at temperature range from (278.15 to 318.15) K were reported. At the same composition of methanol, n-propanol, acetonitrile or 1,4-dioxane and temperature, the mole fraction solubility of ribavirin was highest in {methanol (1) + water (2)} mixtures, and lowest in {1,4-dioxane (1) + water (2)} mixtures. By using the Jouyban-Acree model, van't Hoff-Jouyban-Acree model and Apelblat-Jouyban-Acree model, ribavirin solubility were well correlated obtaining ARD lower than 1.57% and RMSD lower than 2.14  10 À4 . Quantitative values for the local mole fraction of methanol (n-propanol, acetonitrile or 1,4-dioxane) and water around the ribavirin were calculated based on the IKBI method applied to the solubility values determined. Ribavirin was prefer- (2) to methanol (1) + water (2), n-propanol (1) + water (2), acetonitrile (1) + water (2) and 1,4-dioxane (1) + water (2) mixtures at several temperatures. (a), methanol (1) + water (2); (b), n-propanol (1) + water (2); (c), acetonitrile (1) + water (2); (d), 1,4-dioxane (1) + water (2). entially solvated by water for the four solvent mixtures over a very wide range of compositions. The preferential solvation magnitude of ribavirin by water was highest in acetonitrile mixtures, followed by 1,4-dioxane and n-propanol mixtures, and finally, by methanol mixtures. It is conjecturable that ribavirin could be acting mainly as a Lewis base in front to water. Handbook of Solubility Data for Pharmaceuticals Encyclopedia of Pharmaceutical Technology Absorption and drug development, solubility, Permeability and Charge State Pharmaceutics: The Science of Dosage Forms Design Further calculations on solubility of 3-amino-1-adamantanol in ethanol + water binary solvent mixtures at various temperatures Solvent Mixtures: Properties and Selective Solvation Preferential solvation in mixed solvents Design and synthesis of HCV agents with sequential triple inhibitory potentials Therapeutic efficacy of liposome-encapsulated ribavirin and muramyl tripeptide in experimental infection with influenza or herpes simplex virus Combination antiviral therapy for respiratory virus infections Early addition of ribavirin to interferon in chronic hepatitis C not responsive to interferon monotherapy A new effective method for the synthesis of 1,2,4-triazole-3-carboxamide and ribavirin derivatives An improved procedure for the preparation of ribavirin The arsenolysis reaction in the biotechnological synthesis of ribavirin. The in vitro and in vivo inhibition of influenza A virus replication with a combination of ribavirin and ozeltamivir Ribavirin for the treatment of chronic hepatitis C virus infection: a review of the proposed mechanisms of action Bioproduction of ribavirin by green microbial biotransformation Ribavirin as an anti-cancer therapy: acute myeloid leukemia and beyond Ribavirin treatment effects on breast cancers overexpressing eIF4E, a biomarker with prognostic specificity for luminal Btype breast cancer BDDCS Applied to over 900 Drugs Biowaiver monographs for immediate release solid oral dosage forms: Ribavirin Novel high/low solubility classification methods for new molecular entities Pharmaceutical cocrystals of ribavirin with reduced release rates Solid-state transformations of ribavirin as a result of high-shear mechanical processing Formulation and evaluation of immediate release film coated ribavirin tablets Review of the cosolvency models for predicting solubility of drugs in water-cosolvent mixtures Solubility and preferential solvation of the flavonoid naringenin in some aqueous/organic solvent mixtures Preferential solvation of sulfadiazine, sulfamerazine and sulfamethazine in ethanol + water solvent mixtures according to the IKBI method Solubility and preferential solvation of some n-alkylparabens in methanol + water mixtures at 298.15 K Solubility temperature and solvent dependence and preferential solvation of citrus flavonoid naringin in aqueous DMSO mixtures: an experimental and molecular dynamics simulation study Modelling the solubility and preferential solvation of gallic acid in cosolvent + water mixtures Study on enantiomeric separation of basic drugs by NACE in methanol-based medium using erythromycin lactobionate as a chiral selector Formation of nanoparticles of a hydrophilic drug using supercritical carbon dioxide and microencapsulation for sustained release HPLC for Pharmaceutical Scientists Physical Chemical Principles in the Pharmaceutical Sciences In silico prediction of drug solubility in water-dioxane mixtures using the Jouyban-Acree model Tools and techniques for solvent selection: green solvent selection guides Handbook of Aqueous Solubility Data Solubility and solution thermodynamics of some sulfonamides in 1-propanol + water mixtures Acetonitrile-protein interactions: amino acid solubility and preferential solvation Solubility measurement and thermodynamic modeling of 4-nitrophthalimide in twelve pure solvents at elevated temperatures ranging from (273.15 to 323.15) K Solubility determination and thermodynamic models for dehydroepiandrosterone acetate in mixed solvents of (ethyl acetate + methanol), (ethyl acetate + ethanol) and (ethyl acetate + isopropanol) Determination and modeling for solid-liquid phase equilibrium of ternary caprolactam + cyclohexanone oxime + methyl tert-butyl ether system Solubility determination and thermodynamic modeling of terephthaldialdehyde in ten organic solvents from T = (273.15 to 318.15) K and mixing properties of solutions Solubility of N-phenylanthranilic acid in nine organic solvents from T = (283.15 to 318.15) K: Determination and modeling Thermodynamic study of solubility for 2-amino-4-chloro-6-methoxypyrimidine in twelve organic solvents at temperatures from 273.15 K to 323 A new class of synthetic nucleoside analogues with broad-spectrum antiviral properties Design, synthesis, and broad spectrum antiviral activity of 1-betad-ribofuranosyl-1,2,4-triazole-3-carboxamide and related nucleosides An improved thermoanalytical approach to quantifying trace levels of polymorphic impurity in drug powders Handbook of Organic Solvent Properties Comment on ''Measurement and correlation of solubilities of (z)-2-(2-aminothiazol-4-yl)-2-methoxyiminoacetic acid in different pure solvents and binary mixtures of water + (ethanol, methanol, or glycol) The Properties of Solvents Preferential solvation of some structurally related sulfonamides in 1-propanol + water co-solvent mixtures Densities and viscosities of methanol-water mixtures Refractive indices, densities and excess molar volumes of binary systems methanol + water and ethanol + water at T = 293.15 K Densities and viscosities of aqueous solutions of 1-propanol and 2-propanol at temperatures from 293.15 K to 333.15 K Density, viscosity and thermodynamic activation of viscous flow of water1acetonitrile Thermodynamics of aqueous mixtures of nonelectrolytes. IV. Excess volumes of water-acetonitrile mixtures from 15 to 35°C Volumetric properties of the pharmaceutical model cosolvent system 1,4-dioxane + water at several temperatures Theory of preferential solvation of nonelectrolytes Preferential solvation in two-and in three-component systems Solubility and solvation in mixed solvent systems The solvatochromic comparison method. II. The alphascale of solvent hydrogen-bond donor (HBD) acidities Solubility modelling and preferential solvation of paclobutrazol in cosolvent mixtures of (ethanol, npropanol and 1,4-dioxane) + water Solubility and preferential solvation of econazole nitrate in binary solvent mixtures of methanol, ethanol and 1,4-dioxane in water Solubility, solution thermodynamics and preferential solvation for 4-bromopyrazole in mixed solvents of (methanol/ ethanol + water) from T = (283.15 to 318.15) K The authors would like to express their gratitude for the Priority Academic Program Development of Jiangsu Higher Education Institutions. Supplementary data associated with this article can be found, in the online version, at http://dx.doi.org/10.1016/j.jct.2017.07.027.