Solubility Study of Acetylsalicylic Acid in Ethanol + Water Mixtures: Measurement, Mathematical Modeling, and Stability Discussion

Solubility determination of poorly water-soluble drugs is pivotal for formulation scientists when they want to develop a liquid formulation. Performing such a test with different ratios of cosolvents with water is time-consuming and costly. The scarcity of solubility data for poorly water-soluble drugs increases the importance of developing correlation and prediction equations for these mixtures. Therefore, the aim of the current research is to determine the solubility of acetylsalicylic acid in binary mixtures of ethanol+water at 25 and 37°C. Acetylsalicylic acid is non-stable in aqueous solutions and readily hydrolyze to salicylic acid. So, the solubility of acetylsalicylic acid is measured in ethanolic mixtures by HPLC to follow the concentration of produced salicylic acid as well. Moreover, the solubility of acetylsalicylic acid is modeled using different cosolvency equations. The measured solubility data were also predicted using PC-SAFT EOS model. DSC results ruled out any changes in the polymorphic form of acetylsalicylic acid after the solubility test, whereas XRPD results showed some changes in crystallinity of the precipitated acetylsalicylic acid after the solubility test. Fitting the solubility data to the different cosolvency models showed that the mean relative deviation percentage for the Jouyban-Acree model was less than 10.0% showing that this equation is able to obtain accurate solubility data for acetylsalicylic acid in mixtures of ethanol and water. Also, the predicted data with an average mean relative deviation percentage (MRD%) of less than 29.65% show the capability of the PC-SAFT model for predicting solubility data. A brief comparison of the solubilities of structurally related solutes to acetylsalicylic acid was also provided.


INTRODUCTION
It has been reported that non-steroidal anti-inflammatory drugs are among the most commonly used drugs (1). Based on the biopharmaceutics classification system (BCS) proposed by Amidon et al. (2), these drugs belong to class II, characterized by low solubility and high permeability (3). Therefore, dissolution/solubility plays a key role in better absorption and the fast dissolution of BCS class II drugs in the gastrointestinal tract (GIT) following oral administration (4). Because of high membrane permeability, the degree of absorption of these drugs could have been as high as 100% if they did not suffer poor solubility/dissolution in the GIT. However, this may be difficult because the drugs display poor aqueous solubility (5). Poor aqueous solubility could be a result of high lipophilicity and strong intermolecular interactions which make the solubilization of the solid energetically costly (6).
The solubility of a drug in the GIT is one of the main parameters in formulating a drug into a dosage form as it can control bioavailability. So, the poor solubility of a drug in the GIT can present a key issue during drug formulation development. Furthermore, the solubility of drugs in water cannot be ignored in the formulation development of liquid, parenteral formulations, and soft gelatine capsules. In the case of poorly water-soluble drugs, the solubility should be enhanced to enable the formulator to develop liquid dosage forms. This is the main reason that researchers are interested in exploring the solubility profile of a specific drug in a specific solvent as it permits the scientists to select the ideal solvent medium for drug formulation. In addition, the solubility knowledge can help overcome problems arising in the re-crystallization of drugs by changing solvent polarity and/or temperature of the solution. In addition, cosolvency has wider application throughout the chemical industry. It has been used for soil remediation to increase the solubility of contaminants in water for easier removal of toxic contaminations such as pesticides (7).
Aqueous solubility of acetylsalicylic acid at 278.15-345.15 K was reported by Apelblat and Manzurola (8). The enhancement effects of a number of surfactants on the aqueous solubility of acetylsalicylic acid at 37°C (9) and also the solubility of acetylsalicylic acid in a number of monosolvents at various temperatures were reported (10,11). Effects of four hydrotropes on the solubilization of acetylsalicylic acid at different temperatures were also investigated (12).
The calculation of drug solubility in binary mixtures of solvents is paramount important as such the knowledge could provide useful information to researchers to find the binary solvents which are capable of dissolving more drugs. Generally, in the formulation development of liquid dosage forms, the optimum concentration of the cosolvent in water is normally determined by trial and error experimentations. But this approach not only is time-consuming but also is a costly process. To overcome these two drawbacks, cosolvency data can be modeled for prediction purposes. Such cosolvency models often predict a given drug's solubility in various fractions of a single cosolvent-water mixture based on the known solubility of the drug in each of the two neat solvents and solubility in several cosolvent-water fractions (13). Examples of these cosolvency models are the extended Hildebrand solubility approach (EHS) (14), the log-linear model of Yalkowsky (15), the excess free energy approach (16), the Jouyban-Acree model (17) as an extended version of the combined nearly ideal binary solvent/Redlich-Kister equation (CNIBS/R-K) (18), the general single model (GSM) (19), and the mixture response surface methods (MR-S) (20). Although the aforementioned cosolvency models are very useful, some of them lack generalization to other drugs' solubility in the same cosolvent-water mixtures, or extrapolation to other cosolvents, when the solubility of the drug in the neat solvents and the solubilizing power of the cosolvents are not known.
There have been efforts to predict cosolvency of various solvents by using physicochemical properties of the solvents (21,22), which allow prediction of a given drug's solubility in various cosolvent-water mixtures. Moreover, the solubilization profiles of a set of cosolvents (mixed with water) towards a group of chemically unrelated drugs have been shown to differ according to the hydrophobicity of the solutes (23) or Abraham parameters of the solvents and drug (24). By employing the quantitative structure-property relationship (QSPR), it is possible to estimate the solubility of drugs from the properties of the molecular structures of the drugs. Comparison of QSPR models for drug solubility in different solvents can identify solute features determining solubility in various solvents or solvent mixtures. One such investigation by Ghafourian and Bozorgi (25) has shown that the impact of drug hydrophobicity (log P) on solubility is reduced at higher volume fractions of PEG in the PEG/water binary mixtures. In addition to the success of predicting solubility data with the mentioned models, researchers have recently paid special attention to the perturbed-chain statistical associating fluid theory (PC-SAFT) model. The main advantage of the PC-SAFT model is that it uses only data from quantum chemical calculations, thus enabling predictions when there are no experimental data available. Until now, the drug solubility in the different solutions such as molecular solvents, CO 2 , ionic liquids, and deep eutectic solvents has been modeled using PC-SAFT (26)(27)(28)(29)(30). However, such studies are few and mostly considered the simplified models or cubic equation of states (EOSs) (25). In this study, using the PC-SAFT EOS, the solubility of acetylsalicylic acid in the presence of binary mixtures of ethanol + water at 25 and 37°C was estimated. Acetylsalicylic acid is non-stable in aqueous solutions and readily hydrolyzes to salicylic acid (31)(32)(33)(34). The hydrolysis rate could be decreased by the addition of polyethylene glycol (PEG) 6000, povidone, or sorbitol (35 (33). It has also been shown that the degradation rate of acetylsalicylic acid in ethanol + water mixtures is increased by increasing ethanol fraction and temperature (31). Concerning this information, the determination of acetylsalicylic acid solubility in ethanol + water mixtures is a challenging topic and needs further consideration. The aim of the current study was to use acetylsalicylic acid as a non-stable model drug and waterethanol as safe solvents to use in pharmaceutical formulation to calculate the solubility of the drug with a minimal error.

MATERIALS AND METHODS
Acetylsalicylic acid (ASA or aspirin) powder with a mass fraction purity of >99% was obtained from Sigma (Sigma-Aldrich, Gillingham, UK). Absolute ethanol was purchased from Fisher Scientific (Loughborough, UK).

Solubility Studies
Eleven combinations of binary solvent mixtures were prepared using suitable volume fractions of water-ethanol increasing consecutively from 0.0 to 1.0 ethanol: water, the neat solvents. A shake-flask technique with spectrophotometry was used to determine drug solubility in these solvents. Excess amounts of drug were added to the mixtures to allow saturation concentration to be reached. The eleven samples were then incubated at 37°C and 25°C in a thermostatic water bath shaker (Cambridge, Crafton) at atmospheric pressure while constantly shaking at a speed of 200 rpm, for a minimum of 72 h (the preliminary results showed that 72 h was sufficient to reach equilibrium condition).
The suspensions were allowed to settle for 1 h, and then the supernatants were filtered to remove the excess solid using a syringe-driven filter unit (pore size 0.20 μm). A volume of either 0.2 mL or 1 mL (depending on the concentration of drug in the saturated sample) of the filtrate from each sample was immediately diluted quantitatively using an appropriate amount of the same ethanol: water solvent mixture. Dilutions in the range of 10-10,000 times were made depending on the concentration of the drug in the filtered solutions. A Knauer HPLC instrument (Berlin, Germany) composed of a K-1001 HPLC pump, a BioTech. degasser, a sample loop (20 μL), and a K-2600 ultraviolet detector was used for the determination of acetylsalicylic acid and salicylic acid. The chromatographic data processing was performed by employing the Chromgate software (version 3.1). The separation was performed by a stationary phase of C 18 XBridge analytical column (5 μm × 250 mm × 4.6 mm) from Waters Co. (Ireland), and the mobile phase consisted of phosphoric acid:acetonitrile:water (2:400:600 v/v). The freshly prepared mobile phase was filtered using a vacuum filter system equipped with a 0.45-μm membrane filter (Millipore Corp., Billerica, Massachusetts) and degassed by ultrasonic for 15 min. Chromatography was run at 25°C by pumping the mobile phase at a flow rate of 1.5 mL/min. The UV detector recorded the column effluent at 254 nm (36). The calibration curve between the peak area and the concentration in the range of 10-1000 mg.L −1 for acetylsalicylic acid is Y = 3336.8 C ASA + 181,813 and in the range of 0.5-20 mg.L −1 for salicylic acid is Y = 39,158 C SA + 14,652. The solubility data are obtained from the interpolation of these plots.

Density Measurement of the Solutions
In order to measure the true density of each saturated solution, 1 mL of each filtrate was weighed accurately using an accurate balance with a precision of 0.1 mg (Sartorius, Ireland). The true density values are needed to convert molarity and mole fraction. Three repeats of the procedure were completed for all seven drugs at 37°C and 25°C, and the results averaged. It should be noted that the density of the saturated solution of acetylsalicylic acid is overestimated since there is a considerable amount of salicylic acid which increases the density value.

Differential Scanning Calorimetry
A differential scanning calorimeter (DSC7, Mettler Toledo, Switzerland) was employed to investigate the thermal behavior (enthalpy and melting point) of acetylsalicylic acid before and after the solubility test. This information allows us to identify whether a different polymorphic form is produced during the solubility test. The samples studied through the DSC machine were the pure drug sample and the samples obtained after equilibration with 0, 0.5, and 1 ethanol fractions at both 25 and 37°C temperatures. The acetylsalicylic acid particles left in the solubility test were collected and dried. The dried samples were placed in DSC pans and heated between 25 and 300°C at a scanning rate of 10°C/min under nitrogen gas (50 mL/min). After obtaining the DSC traces for each sample, the melting points and enthalpies of fusion were calculated by the software provided.

Fourier Transform Infrared Spectroscopy
In order to explore any changes in the structure of the solid extracted from the solubility test as a result of possible hydrolysis of acetylsalicylic acid and precipitation of any byproducts, FT-IR was employed (Perkin Elmer's Spectrum One, Shelton, CT, USA). Briefly, methanol was used to clean the instrument to remove any residual matter left on the apparatus, followed by placing a few milligrams of each of the separated solid particles after the solubility test. The sample was pressed with a pressure of 100 bar followed by scanning the sample three times over a range of 4000 cm −1 to 500 cm −1 to obtain spectra with appropriate resolution.

X-ray Powder Diffraction
The XRPD patterns were obtained for all samples including original acetylsalicylic acid using a D2 Phaser diffractometer (Bruker AXS GmbH, Karlsruhe, Germany). All samples produced were scanned in Bragg-Brentano geometry, over a scattering (Bragg, 2θ) angle range from 5 to 50°, in 0.02°steps at 1.5°min −1 (37). Microsoft Excel was used to analyze and plot the collected XRPD patterns. The crystallinity of the samples was also determined to elucidate the effect of the type of solvent (water and ethanol) on the crystallinity of the recovered acetylsalicylic acid samples. The area under the curve for the "distinctive crystalline peaks" at 7.8 and 15.6 2θ angles was measured for each XRPD diffractogram and used in the determination of crystallinity (%) using the equation 1 (38).
Mathematical Modeling of the Solubility Data The solubility values determined for acetylsalicylic acid in ethanol + water mixtures are correlated and backcalculated utilizing the mathematical cosolvency models such as the Yalkowsky, Jouyban-Acree, modified Wilson, and PC-SAFT models, and details for each studied model are discussed below.

Yalkowsky Model
The Yalkowsky model was employed to express the natural logarithm of solubility in a mixture of solvent + cosolvent (39).
where C 1 and C 2 are solubility data in mono-solvents 1 and 2 in molar fraction unit, x m is the solubility of the drug in the solvent mixture, and f 1 and f 2 are volume fractions of solvents 1 and 2 in the absence of the drug. After modification of Eq. 2 (i.e., substitution of f 2 with (1−f 1 ) and subsequent rearrangements), Eq. 3 can be obtained as (40): the σ is the model constant. A linear relationship has been shown between the logarithm of octanol-to-water partition coefficient of the solute (log P) and σ as below (41): where M and N are the cosolvent constants. After replacing Eq. 4 in Eq. 3, a predictive mathematical equation is attained (41).
By employing M and N values obtained from the literature for used cosolvent (ethanol) and log P of a drug, the solubility of the drug in the solvent mixture can be computed only using solubility data in water.

Jouyban-Acree Model
The Jouyban-Acree model as a simple linear cosolvency model used for binary mixtures of solvents at various temperatures can be presented by Eq. 6(17): J i is the model parameter calculated using linear and other model parameters have the same meanings as those of the above model.

The Modified Wilson Model
In addition to linear models employed for fitting and prediction of solubility values, the non-linear model of modified Wilson is also utilized for modeling the solubility data in the solvent mixtures at isothermal conditions. The equation is as (42): λ 12 and λ 21 are the model constants computing using nonlinear analysis.

PC-SAFT Model
The perturbed chain SAFT equation of state (EOS) or PC-SAFT was first proposed and developed by Gross and Sadowski in 2001 (43) as an alternative to the original version of SAFT derived by Chapman et al. (44). The residual molar Helmholtz energy of the PC-SAFT (a res ) obtained by the Helmholtz energy contributions from the reference system hard chain (a hc ), dispersion force (a disp ), and hydrogen bonding (a assoc ) is obtained as follows: In PC-SAFT, pure components can be described using five pure-component parameters: (i)m, number of segments per chain; (ii)σ, diameter of each segment in Angstrom (Å); (iii)ε, energy parameter for each segment in Joules (J); (iv) κ AiBi , effective volume of the association (Å 3 ); (v)ε AiBi , energy parameter of the association (bar.l/mol); ). The parameters of each component are reported in Table I (43,45). The interaction parameters for binary systems (ethanol + water), (acetylsalicylic acid + water) and (acetylsalicylic acid + ethanol), k ij , for the purely predictive model were set to zero.
The fugacity coefficient for component k (ϕ k ) and compressibility factor (z) using the PC-SAFT EOS are computed as follows: where ρ is the molar density. Using PC-SAFT, the activity coefficients are calculated from the fugacity coefficients through Eq. 11: where ϕ i and 0 i are the fugacity coefficients of component i in the mixture and that of the pure compound, respectively. In solid-liquid equilibria, the solid solubility in the liquid phase is calculated according to the following expression (46): where x i and γ i represent the solubility and activity coefficient of compound i. In this study, the activity coefficient of compound i (γ i ) was determined via Eq. 11. Since the activity coefficient depends on solubility in mole fraction (x i ), solubility must be determined from the iterations with Eq. 12. In the equation above, ΔH m and T m represent fusion enthalpy and melting point temperature, respectively, and their values are presented in Table II. All explained models are correlated to the measured solubility values of acetylsalicylic acid, and the mean relative deviation (MRD%) (Eq. 13) is used to obtain the model's accuracy.

MRD%
N is the number of data points in each set.

Thermal Analysis
The melting point and enthalpy of fusion were determined for pure acetylsalicylic acid as well as for the samples that had undergone solubility experiments following treatment with 0%, 50%, and 100% ethanol. Since the solid phase remaining at the end of the solubility experiments could have altered to the different polymorphic forms as a result of the equilibrium state during the solubility experiments, the enthalpy and melting point of the solid particles obtained after each solubility test were compared to the original drug crystals. Acetylsalicylic acid has two polymorphic forms known as polymorph I and polymorph II. Polymorph I which is the stable form of acetylsalicylic acid melts above 137°C (47), whereas polymorph II melts between 128 and 130°C (48). The DSC traces showed that the concentration of ethanol + water had no effect on the type of polymorphic form of acetylsalicylic acid since the filtered acetylsalicylic acid after the solubility experiments showed a melting peak above 137°C. The results of Table II ruled out the presence of polymorph II in the acetylsalicylic acid samples as the melting points of all samples are above 137°C hence portraying the absence of changes in the polymorphic form of acetylsalicylic acid before and after the solubility test.
Although some authors suggested that the onset melting temperature of the solute is preferred over the peak melting temperature (49), others such as El-Badryet al. (50) employed the peak temperature to represent the melting point of the drugs. Likewise, this study utilized the peak melting  other. This could be due to changes in the particle size or the crystallinity of the original acetylsalicylic acid particles after the solubility test where some of the dissolved acetylsalicylic acid particles can undergo recrystallization leading to different particle sizes and different crystallinity percentages. It has been shown that, generally, smaller particles should melt at a lower temperature than when larger particles are used for the given material (51). This is true due to the high surface area of small particles that would be in contact with the heat during DSC hence lowering the melting point. In order to explore any changes in the structure of acetylsalicylic acid crystals as a result of hydrolysis (degradation), FT-IR was carried out on the extracted acetylsalicylic acid after the solubility test. FT-IR spectra of all samples showed that all the diagnostic peaks for acetylsalicylic acid can be also detected in the samples collected after the solubility test at different temperatures ( Figure 1). Furthermore, the FT-IR spectra of the collected samples after the solubility test did not show any peak relevant to the presence of salicylic acid. This indicates that salicylic acid is not reached above the saturated solubility to precipitate. Table II shows that the enthalpy of fusion is different compared to the enthalpy of untreated acetylsalicylic acid. For example, the enthalpy of fusion of unprocessed acetylsalicylic acid was 165.6 ± 6.2 J/g, whereas this value was 130.7±2.8 J/g when 100% water was used in the solubility test. Bustamante et al. (52) suggested that the small differences in enthalpy of fusion values observed after contact with the solvents are not sufficiently important and thermodynamic activity can be assumed constant in the solubility model which could be true when 50:50 water:ethanol or 100% ethanol was used in the solubility test as their enthalpy is closer to the enthalpy of untreated acetylsalicylic acid. Table II shows that the enthalpy difference between unprocessed acetylsalicylic acid and when 100% water was used in the solubility test is significant. It has been reported that the significant difference in the enthalpy could be due to changes in the crystallinity of acetylsalicylic acid after the solubility test. This has been shown when acetylsalicylic acid has been crystallized from various organic solvents (53). The enthalpy data in Table II shows that the acetylsalicylic acid crystals   a f 1 is volume fraction of ethanol in the ethanol and water mixtures in the absence of acetylsalicylic acid b Molar concentration produced from hydrolysis of acetylsalicylic acid. Data in bracket are related to saturated concentration for salicylic acid in ethanol + water mixtures and taken from a reference (56) filtered after the solubility test in 100% water have less crystallinity compared to untreated acetylsalicylic acid.
Since the crystallinity of acetylsalicylic acid could change the solubility, XRPD was used to investigate any changes in the crystallinity of the collected acetylsalicylic acid samples compared to the original acetylsalicylic acid after the solubility test (Figure 2). The XRPD of the samples shows that the only difference between the samples is the intensity of peaks at 7.8 and 15.6 2θ. The change in the intensity of the peaks could be due to changes in the crystallinity of the samples collected from different solvents at different temperatures. The relative crystallinity of the samples was calculated (Table III), and the results showed that the crystallinity of the samples collected at 25°C has not significantly changed (ranging from 32.3 to 35.5%), but this was not the case for the samples collected at 37°C. The crystallinity of the acetylsalicylic acid samples collected from 100% water was significantly lower than that of other samples which could increase the solubility of acetylsalicylic acid during the solubility test. This is in agreement with melting point and enthalpy data (Table II) where the sample with the lowest crystallinity showed the lowest melting point and enthalpy (Table II; Figures 3 and 4). Since drugs with low crystallinity could have a higher solubility, the current study suggests that it is worth measuring the crystallinity of samples after the  solubility test for better interpretation of the data in the solubility test. In conclusion, the finding of this study regarding DSC and XRPD proved that if a drug is prone to any changes in its crystallinity during the solubility test, then the obtained solubility data could not be an exact reflection of the solubility data for the original crystals which was discussed further in the solubility section of the manuscript.

Solubility Studies
Binary aqueous-cosolvent systems are usually employed in practice to enhance the solubility of poorly water-soluble drugs in order to formulate them as a liquid dosage form for patients. Organic cosolvents particularly ethanol are among the most powerful pharmaceutical solubilizing agents which are usually used in the formulation of elixirs. The prediction of solubility trends in aqueous mixtures of ethanol is of most attention and it facilitates finding all cosolvent systems (54). The molar fraction solubility of acetylsalicylic acid in the mixed solutions with various ethanol fractions at 25 and 37°C along with standard deviations of repeated measurements are reported in Table IV. This table shows that when the volume fraction of ethanol enhances, the solubility of the drug increases and reaches a peak at the ratio of 90:10 ethanol:water at both temperatures (25 and 37°C), after which the solubility decreases. Moreover, the molar fraction solubility of acetylsalicylic acid increases when the temperature of the solution increases from 25 to 37°C at the same composition of the mixed solutions. Molar solubility value in neat water in this study at 25°C (0.0246 mol.L −1 ) and 37°C (0.0308 mol.L −1 ) is in good agreement with that reported in the literature for 25°C (0.0255 mol.L −1 ) and 37°C (0.0357 mol.L −1 ) (55) and the small observed difference can be related to the person to person error and the employed methodology. Furthermore, Table IV shows also salicylic acid concentration in the acetylsalicylic acid saturated mixtures at the investigated temperatures. The saturated concentrations of salicylic acid in the ethanol + water (56) mixtures are also reported in this table. As can be seen, salicylic acid in these mixed solutions is under saturated concentration. However, the presence of a considerable concentration of salicylic acid in the saturated solutions of acetylsalicylic acid in ethanol + water mixtures may pass some further changes in its intrinsic solubility values. Figure 5 compares the solubility profiles of acetylsalicylic acid (used in the current study) and other related compounds in these mixtures (i.e., salicylic acid (56), 5-amino salicylic acid (57), 3-amino salicylic acid (58), and 3,5-dinitrosalicylic acid (59)) at 25°C. As can be seen, the solubility difference between 5-amino salicylic acid and 3-amino salicylic acid is relatively small and between ethanol mass fractions of 0.6 and 1.0, and the solubility of 3-amino salicylic acid exceeded 5amino salicylic acid. The solubility of 5-amino salicylic acid and 3-amino salicylic acid is weaker than that of acetylsalicylic acid, 3,5-dinitrosalicylic acid, and salicylic acid in the ethanol and water mixtures and shows a maximum solubility in lower ethanol mass fractions, whereas acetylsalicylic acid represents a maximum solubility at f 1 = 0.9 (or mass fraction of 0.877) of ethanol and 3,5-dinitrosalicylic acid and salicylic acid show it at neat ethanol. For salicylic acid with log P of 2.26, and 3,5dinitrosalicylic acid with log P of 1.71, these behaviors can be easily explained. However, for acetylsalicylic acid (log P of 1.18) and 5-amino salicylic acid with similar log P (=1.2), another parameter can affect the solubility such as steric hindrance. Generally, the solubility profile of a solute in a mixed solution is influenced by the factors of polarity, steric hindrance, hydrogen bonding, intermolecular interactions, and van der Waals' forces between solute-solvent, solutesolute, and solvent-solvent, etc. and it is too difficult to elucidate the solubility behavior of a solute in terms of a single reason. All log P values were taken from DrugBank database (https://go.drugbank.com). Although the density is used for converting molar solubility data to mole fraction data, in the current study considering the instability of solute, salicylic acid also exists in the saturated solutions. So, the measured density data can not exactly be related to acetylsalicylic acid saturated solutions. However, in this section and for comparing its solubility with other salicylate derivatives, the authors used these density data for converting the molar solubility to mole fraction solubility for expressing in the same unit with other systems. It should be noted that the measured density data are considered as approximate values for acetylsalicylic acid saturated solutions considering the assumption that salicylic acid does not have a significant effect on the density value. a The measured density data are considered as approximate values for acetylsalicylic acid saturated solutions and used for converting the molar data to mole fraction data with considering this assumption that as salicylic acid is not reached above the saturated solubility in these solutions, it does not have any effect on the density value b f 1 is volume fraction of ethanol in the ethanol and water mixtures in the absence of acetylsalicylic acid can be seen that introducing the solubility parameters of Abraham is able to enhance the prediction capability of the model from 33.9 to 16.8%.

CONCLUSION
In the current research, the solubility profile of acetylsalicylic acid in the binary mixtures of ethanol + water at 25 and 37°C are measured by a simple shake-flask method and the experimental data are correlated and back-calculated by Yalkowsky, Jouyban-Acree, and the modified Wilson models. The results proved that the solubility of acetylsalicylic acid at various temperatures in binary mixtures of ethanol + water can be correlated well using the Jouyban-Acree and the modified Wilson models. This is important for the pharmaceutical industry as it can save time and reduce the cost to perform the solubility test. TheJouyban-Acree-based general cosolvency models are also extended to predict solubility data of acetylsalicylic acid by considering the fact that none of the data points in the mixed solvents of this study is used in the training process of models. During the last decades, the accuracy of the cosolvency models was improved by a factor of 3-4 as it is confirmed in this work employing solubility data of acetylsalicylic acid where the prediction error of the model presented in 1998 (i.e., Eq. 4 with the prediction error of 44.5%) was reduced to 16.8% for the model presented in 2007 (i.e., Eq. 15). It is obvious that 16.8% prediction error could be well tolerated in many practical process design computations in the pharmaceutical industry. Since PC-SAFT EOS shows an acceptable function in all ranges of systems, its applicability was explored for pharmaceuticals. The results of the PC-SAFT EOS for predicting the solubility data of acetylsalicylic acid in the binary mixtures of ethanol + water at 25 and 37°C are in agreement with experimental data, suggesting the good performance of the PC-SAFT EOS for pharmaceutical purposes. For most pharmaceutical scientists, the complicated computations associated with PC-SAFT EOS are a limiting factor for its practical applications.