Determination and Data Correlation of Solubility of Sofosbuvir Polymorphs in Ethyl Acetate + Toluene and Methyl tert -Butyl Ether Binary Solvents at the Temperature Range from 268.15 to 308.15 K Determination and Data Correlation of Solubility of Sofosbuvir Polymorphs in Ethyl Acetate + Toluene and Methyl tert -Butyl Ether Binary Solvents at the Temperature Range from 268.15 to 308.15 K

: A gravimetric method was used to experimentally determine the (solid + liquid) equilibrium of sofosbuvir of crystalline forms A and B in both ethyl acetate + toluene and methyl tert -butyl ether (MTBE) + toluene binary solvents systems at atmosphere pressure. Experiments were carried out at a temperature range of 268.15 − 308.15 K. Results show that the solubility of sofosbuvir increases with temperature, and the solubility of form B was higher than that of form A. The modiﬁed Apelblat model, the CNIBS / Redlich–Kister model, and the combined version of Jouyban–Acree model were employed to correlate the measured solubility data, respectively. Furthermore, an examination of the solid-state stability of the two polymorphs was conducted, ﬁnding that form A and form B exhibit good solid-state stability under high temperature, high humidity, and strong light exposure conditions. Abstract: A gravimetric method was used to experimentally determine the (solid + liquid) equilibrium of sofosbuvir of crystalline forms A and B in both ethyl acetate + toluene and methyl tert -butyl ether (MTBE) + toluene binary solvents systems at atmosphere pressure. Experiments were carried out at a temperature range of 268.15−308.15 K. Results show that the solubility of sofosbuvir increases with temperature, and the solubility of form B was higher than that of form A. The modified Apelblat model, the CNIBS/Redlich–Kister model, and the combined version of Jouyban–Acree model were employed to correlate the measured solubility data, respectively. Furthermore, an examination of the solid-state stability of the two polymorphs was conducted, finding that form A and form B exhibit good solid-state stability under high temperature, high humidity, and strong light exposure conditions.


Introduction
Sofosbuvir (C 22 H 29 FN 3 O 9 P, CAS No.: 1190307-88-0, presented in Figure 1) is a novel anti-HCV agent, offering a more satisfactory sustained virologic response rate than those that have ever been used in similar cases before [1,2]. It is indicated by our previous studies that ethyl acetate + toluene and MTBE + toluene binary systems can be used to crystallize sofosbuvir of form A or B [3]. However, there is no literature that has reported the solution thermodynamic data of sofosbuvir in these two systems. These are important for the optimization of the crystallization process, especially for the controlling of polymorphs, since the various solid forms generally present different physicochemical properties.

Introduction
Sofosbuvir (C22H29FN3O9P, CAS No.: 1190307-88-0, presented in Figure 1) is a novel anti-HCV agent, offering a more satisfactory sustained virologic response rate than those that have ever been used in similar cases before [1,2]. It is indicated by our previous studies that ethyl acetate + toluene and MTBE + toluene binary systems can be used to crystallize sofosbuvir of form A or B [3]. However, there is no literature that has reported the solution thermodynamic data of sofosbuvir in these two systems. These are important for the optimization of the crystallization process, especially for the controlling of polymorphs, since the various solid forms generally present different physicochemical properties.  In this study, a gravimetric method [4][5][6] was used to acquire the solubility data of sofosbuvir of forms A and B in both ethyl acetate + toluene binary solvents with ethyl acetate mole fraction Crystals 2020, 10

Modeling
Correlation of the experimental solubility data with different models is helpful for further understanding the thermodynamic properties in the measurement range. With this consideration, three models were used in this study. The modified Apelblat equation describes the dependence of solubility on temperature, the CNIBS/ Redlich-Kister model describes the dependence of solubility on solvent composition, while the Jouyban-Acree model describes the dependence of solubility on both parameters.

Modified Apelblat Equation
The Apelblat equation, as a well-known equation used to correlate solubility, was originally proposed by Apelblat and Manzurola; in it, the temperature dependence of the mole fraction solubility in different solvents is described as follows [7][8][9][10][11]: where x A is the mole fraction solubility and T represents the corresponding absolute temperature. A, B and C are semi-empirical parameters of the model.

CNIBS/Redlich-Kister Model
The governing equation of CNIBS/Redlich-Kister model [12,13] is defined as follows: where x B 0 and x C 0 are the mole fraction of each solvent in a binary solvent system, respectively; N is the number of "curve-fit" parameters and S i is the model constant, which is two for a binary solvent system in this case; thus, x C 0 can be represented as (1−x B 0 ). Substituting x C 0 , Equation (2) can be expressed as: Changing constants in Equation (3) with a constant term A i (A 0 , A 1 , A 2 , A 3 , A 4 ), the CNIBS/R-K model can be simplified as follows [14,15]:

Jouyban-Acree Model
The original function of Jouyban-Acree model, which was obtained by modifying the CNIBS/R-K model, is shown below [16]: in which J i is a constant of the function. In order to enlarge the applicability of solution behavior of non-ideal systems, the modified Apelblat equation can be used to substitute ln(x A ) i ; then Equation (6) is obtained, as in the CNIBS/R-K model: ln A final equation of the Jouyban-Acree model can be obtained by simplifying equation 6 as well, with a constant term B i (containing B 0 , B 1 , B 2 , B 3 , B 4 , B 5 , B 6 , B 7 , B 8 ) [17]:

Materials.
Sofosbuvir of forms A and B was prepared and identified using the methods published in previous studies [3]. Ethyl acetate, MTBE and toluene were of analytical grade and used without further purification. Detailed information on the above-mentioned materials is listed in Table 1.

Characterization of Sofosbuvir Polymorphs
Powder X-ray diffraction (PXRD) was used to identify the crystalline form of sofosbuvir, and the patterns were recorded using a Rigaku Ultima IV diffractometer with Cu Kα radiation (40 kV, 40 mA) scanned at 20 • /min over an angular range of 5-45 • of 2θ. The patterns are shown in Figure 2, and the main PXRD data are listed in Table S1, Supplementary Information.

Temperature
The two polymorphs of sofosbuvir were kept for 10 days at a temperature of 60 ± 2 °C and sampled on the 5th and 10th day, respectively. PXRD analysis shows no obvious changes, which means that the solid-state stability of the two polymorphs is good at high temperature conditions. The patterns are given in Figures S1 and S2, Supplementary Information.

Humidity
The two polymorphs of sofosbuvir were kept for 10 days at a humidity of RH 92.5 ± 5% and sampled on the 5th and 10th day, respectively. PXRD analysis shows no obvious changes, which means that the solid-state stability of the two polymorphs is good at high humidity conditions. The patterns are given in Figure S3 and S4, Supplementary Information.

Light
The two polymorphs of sofosbuvir were kept for 10 days under a light illumination of 4500 ± 500 lux and sampled on the 5th and 10th day, respectively. PXRD analysis shows no obvious changes, which means that the solid-state stability of the two polymorphs is good at strong light exposure conditions. The patterns are given in Figures S5 and S6, Supplementary Information.

Solubility Measurement
Sufficient amounts of sofosbuvir polymorphs were added to 25 ml solvent in a 50 ml jacket glass vessel with a magnetic stir bar and a thermometer, to form a saturated solution. A thermostatic bath was employed to control the temperature with an accuracy of ±0.1 K. The solution was agitated for 2 h to ensure that the system had reached the state of solid-liquid equilibrium. Then the stirring was stopped and the solution was kept stationary for 10 min. Subsequently, the upper clear solution (about 5 mL) was transferred into a small beaker using a membrane filter (0.45 μm). The beaker was dried in a vacuum oven at 333.15 K for about 48 h to make sure that the weight of the beakers became

Temperature
The two polymorphs of sofosbuvir were kept for 10 days at a temperature of 60 ± 2 • C and sampled on the 5th and 10th day, respectively. PXRD analysis shows no obvious changes, which means that the solid-state stability of the two polymorphs is good at high temperature conditions. The patterns are given in Figures S1 and S2, Supplementary Information.

Humidity
The two polymorphs of sofosbuvir were kept for 10 days at a humidity of RH 92.5 ± 5% and sampled on the 5th and 10th day, respectively. PXRD analysis shows no obvious changes, which means that the solid-state stability of the two polymorphs is good at high humidity conditions. The patterns are given in Figures S3 and S4, Supplementary Information.

Light
The two polymorphs of sofosbuvir were kept for 10 days under a light illumination of 4500 ± 500 lux and sampled on the 5th and 10th day, respectively. PXRD analysis shows no obvious changes, which means that the solid-state stability of the two polymorphs is good at strong light exposure conditions. The patterns are given in Figures S5 and S6, Supplementary Information.

Solubility Measurement
Sufficient amounts of sofosbuvir polymorphs were added to 25 ml solvent in a 50 ml jacket glass vessel with a magnetic stir bar and a thermometer, to form a saturated solution. A thermostatic bath was employed to control the temperature with an accuracy of ±0.1 K. The solution was agitated for 2 h to ensure that the system had reached the state of solid-liquid equilibrium. Then the stirring was stopped and the solution was kept stationary for 10 min. Subsequently, the upper clear solution (about 5 mL) was transferred into a small beaker using a membrane filter (0.45 µm). The beaker was dried in a vacuum oven at 333.15 K for about 48 h to make sure that the weight of the beakers became constant. Both of the polymorphs were proved by PXRD to be maintained over the temperature range. Measurements were conducted in triplicate at each temperature to minimize the relative deviation. All of the masses were measured using an analytical balance (Mettler Toledo XS105, Switzerland) with an accuracy of ±0.01 mg.
The solubility of the sofosbuvir polymorphs, described in mole fraction x in different systems, was calculated by the following equation [18]: in which m and m s represent the mass of sofosbuvir and the solvent, respectively; M and M s are the molecular mass of sofosbuvir and the solvent, respectively.

Solubility Data in Binary Solvents
Experimental data on solubility measurement are presented in Tables S2−S5, Supplementary Information, and graphically illustrated in Figure 3; Figure 4. It is obvious that the solubility of sofosbuvir increases with the temperature and the mole fraction of ethyl acetate or MTBE. Solubility of form B is higher than that of form A under all of the experimental conditions, which means that form A is more stable than form B at the temperature range of measurement. The relative thermodynamic stability is in line with our previous research [3]. constant. Both of the polymorphs were proved by PXRD to be maintained over the temperature range. Measurements were conducted in triplicate at each temperature to minimize the relative deviation. All of the masses were measured using an analytical balance (Mettler Toledo XS105, Switzerland) with an accuracy of ±0.01 mg. The solubility of the sofosbuvir polymorphs, described in mole fraction x in different systems, was calculated by the following equation [18]: = + ∑ ( ) (8) in which m and ms represent the mass of sofosbuvir and the solvent, respectively; M and Ms are the molecular mass of sofosbuvir and the solvent, respectively.

Solubility Data in Binary Solvents
Experimental data on solubility measurement are presented in Tables S2−S5, Supplementary Information, and graphically illustrated in Figure 3; Figure 4. It is obvious that the solubility of sofosbuvir increases with the temperature and the mole fraction of ethyl acetate or MTBE. Solubility of form B is higher than that of form A under all of the experimental conditions, which means that form A is more stable than form B at the temperature range of measurement. The relative thermodynamic stability is in line with our previous research [3].   constant. Both of the polymorphs were proved by PXRD to be maintained over the temperature range. Measurements were conducted in triplicate at each temperature to minimize the relative deviation. All of the masses were measured using an analytical balance (Mettler Toledo XS105, Switzerland) with an accuracy of ±0.01 mg. The solubility of the sofosbuvir polymorphs, described in mole fraction x in different systems, was calculated by the following equation [18]: = + ∑ ( ) (8) in which m and ms represent the mass of sofosbuvir and the solvent, respectively; M and Ms are the molecular mass of sofosbuvir and the solvent, respectively.

Solubility Data in Binary Solvents
Experimental data on solubility measurement are presented in Tables S2−S5, Supplementary Information, and graphically illustrated in Figure 3; Figure 4. It is obvious that the solubility of sofosbuvir increases with the temperature and the mole fraction of ethyl acetate or MTBE. Solubility of form B is higher than that of form A under all of the experimental conditions, which means that form A is more stable than form B at the temperature range of measurement. The relative thermodynamic stability is in line with our previous research [3].

Data Correlation
The Apelblat model, the CNIBS/Redlich-Kister model, and the Jouyban-Acree model were used to correlate the solubility data, based on the average relative deviation (ARD) defined by Equation (9), which was used to compare among three different models [19].
in which N is the number of the experimental measurement, x A,i is each measured solubility and x A,i cal is each calculated value. The 1stOpt program was applied to calculate the data using Equations (1), (4), and (7). Parameters of data correlation and values of ARD% for each individual modeling system are listed in Tables 2-7, with which the calculated solubility values were obtained and listed in Tables S2-S5, Supplementary Information. It was observed that experimental solubility data are satisfactory, fitting with the calculated solubility values. Values of ARD% calculated by the modified Apelblat model, the CNIBS/R-K model, and the Jouyban-Acree model are less than 2.41%, 0.88%, and 4.06%, respectively, indicating that the three models all correlate well in this system. After comparing calculated stabilities with corresponding experimental ones, the CNIBS/R-K model stood out to be more suitable with a higher accuracy than the other two models.

Conclusions
In this paper, the solubility data of sofosbuvir of forms A and B in ethyl acetate + toluene and MTBE + toluene binary solvent systems were measured at atmosphere pressure with T = 268.15−308.15 K by a gravimetric method. Three models were adopted to correlate the experimental solubility data and all of them manifested satisfactory consistency, especially for CNIBS/R-K, with a higher accuracy. It was found that the solubility of sofosbuvir increases with temperature and the mole fraction of ethyl acetate or MTBE. Solubility of form B was higher than that of form A under all of the experimental conditions, which means form A was more stable than form B in the measurement temperature range. A stability test was conducted as well, finding that the two polymorphs show good solid-state stability under high temperature, high humidity, and strong light exposure conditions. All of the experimental data may provide valuable guidance for the crystallization and purification process of sofosbuvir.