Systematic thermodynamic analysis of apremilast polymorphs via solubility measurement with modeling: Mechanism evaluation through molecular simulation

The polymorphism of apremilast has been investigated. Two polymorphs have been identified and characterized by differential scanning calorimeter, fourier transform infrared spectroscopy, and powder X-ray diffractometer. Solubilities of apremilast forms B and E in three binary solvents of methanol-water, acetonitrile-water, and acetonitrile-methanol have been measured using the static method at a temperature ranging from 288.15 K to 328.15 K under standard atmospheric pressure. Subsequently, the solubility data have been analyzed using the Wilson, NRTL, and UNIQUAC thermodynamic models, respectively. Furthermore, the Gibbs energy of solution and the radial distribution function have been calculated using the molecular simulation method to evaluate the dissolution mechanism. The Gibbs energy of solution reveals that the rank of solute-solvent interaction correlated well with solubility order in binary solvent mixtures, and the radial distribution function indicates that weakening of solvent-solvent interaction led to an increase in solubility.