Thermodynamics of sublimation, crystal lattice energies, and crystal structures of racemates and enantiomers: (+)- and (+/-)-ibuprofen

Health is the greatest wealth: Quest for a diagnostic check for thermochemistry of pure drug compounds

The development of pharmaceutical formulations and the optimisation of drug synthesis are not possible without knowledge of thermodynamics. At the same time, the quantity and quality of the available data is not at a level that meets modern requirements. A convenient diagnostic approach is desirable to assess the quality of available experimental thermodynamic data of drugs. A comprehensive set of available data on phase transitions of profens family drugs was analysed using new complementary measurements and structure-property correlations. The consistent sets of solid-gas, liquid-gas and solid-liquid phase transitions were evaluated for twelve active pharmaceutical ingredients based on alkanoic acid derivatives and recommended for the calculations of the pharmaceutical processes. A "centerpiece approach" proposed in this work helped to perform the "health check" of the thermochemical data. The evaluated data on the sublimation enthalpies were used to derive the crystal lattice energies of the profens and to correlate the water solubilities with the sublimation vapour pressures and molecular parameters. A "paper-and-pen" approach proposed in this work can be extended to the diagnosis of "sick" or "healthy" thermodynamic data for drugs with a different structure than those studied in this work.

Quantum Chemistry and Pharmacy: Diagnostic Check of the Thermochemistry of Ibuprofen

The thermodynamic data on ibuprofen available in the literature shows that the disarray of experimental results is unacceptable for this very important drug. The data on ibuprofens available in the literature were collected, combined with our complementary experimental results and evaluated. The enthalpies of combustion and formation of the crystalline RS-(±)- and S-(+)-ibuprofens were measured using high-precision combustion calorimetry. The temperature dependence of the vapour pressure of S-(+)-ibuprofen was measured using the transpiration method and the enthalpy of vaporization was derived from this measurement. The enthalpies of fusion of both compounds were measured using DSC. The G4 calculations have been carried out to determine the enthalpy of formation in the gaseous state of the most stable conformer. Thermochemical properties of the compounds studied were evaluated and tested for consistency with the "centerpiece approach". A set of reliable and consistent values of thermodynamic properties of ibuprofens at 298.15 K is recommended for thermochemical calculations of the pharmaceutical processes. The diagnostic protocol was developed to distinguish between the "sick" or "healthy" thermodynamic data. This diagnostic is also applicable to other drugs with a different structure than ibuprofen.

Influence of Controlled Chirality on the Crystallization of Maleimide-Functionalized 3,4-Ethylenedioxythiophene (EDOT-MA) Monomers

Conjugated poly(alkoxythiophenes) such as poly(3,4-ethylenedioxythiophene) (PEDOT) have attracted considerable interest for use in a variety of applications such as biomedical devices, energy storage, and chemical sensing. Functionalized versions of the 3,4-ethylenedioxythiophene (EDOT) monomer make it possible to create polymers with properties tailored for specific applications. The maleimide functional group shows particular promise due to the wide variety of chemical modifications that it can undergo. Here, we examine the role that control of the chirality of the maleimide (MA) substituent has on the crystal structure and crystallization of the EDOT-MA monomer. We describe a method for the synthesis of a homochiral (S) variant of EDOT-MA and compare its crystallography, morphology, and thermal properties to that of the (R,S) EDOT-MA racemic compound. The conformation of the EDOT-MA molecule was substantially different, with the molecules adopting an "L" shape in the homochiral crystal, while in the racemic crystals, they were more colinear. The thermal stability of the homochiral crystals (Tm = 128.6 °C) was slightly higher than the racemic ones (Tm = 102.8 °C). We expect these results to be important in better understanding the solid-state assembly of the corresponding polymers prepared from these monomers.

Sublimation of Drugs from the Site of Application of Topical Products

The objective of the project was to investigate the plausibility of active pharmaceutical ingredients (APIs) to undergo sublimation from topical application following evaporation of solvent. Topical formulations with different APIs were subjected to a sublimation screening test. The APIs in the selected topical products were found to undergo sublimation to a different extent. The salicylic acid topical product was found to undergo a significant loss due to sublimation. The extent of sublimation of salicylic acid was significantly greater at skin temperature compared to room temperature. When the APIs were subjected to the sublimation screening test in their neat form at 32 ± 1 °C, the natural log of the rate of sublimation decreased linearly with the standard enthalpy of sublimation of compound (R² = 0.89). The formulation composition was found to have a significant impact on the extent of sublimation of the representative API, salicylic acid. The sublimation of APIs from the topical product was found to affect the mass balance studies in the case of the salicylic acid ointment. Furthermore, the results of the human studies agreed with the in vitro experimental results demonstrating the plausibility of loss of API due to sublimation from the site of application.

Structure and Glass Transition Temperature of Amorphous Dispersions of Model Pharmaceuticals with Nucleobases from Molecular Dynamics

Glass transition temperature (Tg) is an important material property, which predetermines the kinetic stability of amorphous solids. In the context of active pharmaceutical ingredients (API), there is motivation to maximize their Tg by forming amorphous mixtures with other chemicals, labeled excipients. Molecular dynamics simulations are a natural computational tool to investigate the relationships between structure, dynamics, and cohesion of amorphous materials with an all-atom resolution. This work presents a computational study, addressing primarily the predictions of the glass transition temperatures of four selected API (carbamazepine, racemic ibuprofen, indomethacin, and naproxen) with two nucleobases (adenine and cytosine). Since the classical non-polarizable simulations fail to reach the quantitative accuracy of the predicted Tg, analyses of internal dynamics, hydrogen bonding, and cohesive forces in bulk phases of pure API and their mixtures with the nucleobases are performed to interpret the predicted trends. This manuscript reveals the method for a systematic search of beneficial pairs of API and excipients (with maximum Tg when mixed). Monitoring of transport and cohesive properties of API-excipients systems via molecular simulation will enable the design of such API formulations more efficiently in the future.

Understanding dispersion interactions in molecular chemistry

This themed collection contains a selection of articles on the topic of Understanding Dispersion Interactions in Molecular Chemistry.

Adsorption of a poorly water-soluble drug onto porous calcium silicate by the sealed heating method

To improve the bioavailability of orally-administered drug, solubilization of poorly water-soluble drug is important. The solubility of a drug in its amorphous form is known to be higher than in its crystalline form. In this study, we attempted to adsorb a sublimable drug onto porous calcium silicate (Florite®, FLR) or non-porous calcium silicate (NPCS) by the sealed heating (SH) method and evaporated (EV) method. Ibuprofen (IBU) was used as the poorly water-soluble, sublimable drug. The physicochemical properties of samples were investigated using powder X-ray diffractometry (PXRD), differential scanning calorimetry (DSC), Fourier-transform infrared spectroscopy (FTIR), and dissolution test. From the PXRD and DSC results, IBU crystals were adsorbed and amorphized by the SH and EV methods with FLR. From the results of FTIR, the shift to a higher frequency by the carbonyl stretching vibration band of IBU suggests an intermolecular interaction between IBU and FLR. In SH with FLR, improved solubilization was observed. IBU adsorbed onto FLR showed a greater dissolution rate than the IBU crystals or NPCS. Thus, the petal-like structure of FLR may be an effective method to adsorb IBU onto FLR using the SH method.

Possible Physical Basis of Mirror Symmetry Effect in Racemic Mixtures of Enantiomers: From Wallach’s Rule, Nonlinear Effects, B–Z DNA Transition, and Similar Phenomena to Mirror Symmetry Effects of Chiral Objects

Effects associated with mirror symmetry may be underlying for a number of phenomena in chemistry and physics. Increase in the density and melting point of the 50%L/50%D collection of enantiomers of a different sign (Wallach’s rule) is probably based on a physical effect of the mirror image. The catalytic activity of metal complexes with racemic ligands differs from the corresponding complexes with enantiomers as well (nonlinear effect). A similar difference in the physical properties of enantiomers and racemate underlies L/D inversion points of linear helical macromolecules, helical nanocrystals of magnetite and boron nitride etc., B–Z DNA transition and phenomenon of mirror neurons may have a similar nature. Here we propose an explanation of the Wallach effect along with some similar chemical, physical, and biological phenomena related to mirror image.

Dissecting intermolecular interactions in the condensed phase of ibuprofen and related compounds: the specific role and quantification of hydrogen bonding and dispersion forces

Ibuprofen is a well-established non-steroidal anti-inflammatory drug, inhibiting the prostaglandin-endoperoxide synthase. One of the key features defining the ibuprofen structure is the doubly intermolecular O-HO[double bond, length as m-dash]C hydrogen bond in cyclic dimers as know from carboxylic acids and confirmed by X-ray analysis. Until now, there was neither information about the vaporization enthalpy of ibuprofen nor about how this thermal property is determined by the subtle balance between different types of intermolecular interaction. In this study we derive the vaporization enthalpy of ibuprofen from thermochemical experiments to be . We dissected the hydrogen bond energy, EHB = 45.0 kJ mol-1, exclusively from measured vaporization enthalpies of related aliphatic carboxylic acids, their homomorph methyl esters and alkyl acetates, respectively. This contribution from hydrogen bonding could be confirmed almost quantitatively from quantum chemical calculations of ibuprofen clusters, which also suggest dispersion interaction of similar order (Edisp = 47 kJ mol-1). Following the full analysis of the gas-vapor transition enthalpy, we studied the changing structural components from the solid to the liquid phase of ibuprofen by means of Attenuated Total Reflection Infrared (ATR-IR) spectroscopy. The cyclic dimers as observed in the X-ray patterns are essentially preserved in the liquid state just above the melting point. However, with increasing temperature the doubly hydrogen-bonded cyclic dimers are replaced by singly hydrogen-bonded linear dimers in the liquid ibuprofen. The transfer enthalpy from the temperature-dependent equilibria of both dimers as obtained from the IR intensity ratios of the vibrational bands quantifies for the first time the energy of the released, single hydrogen bond to be EHB = 21.0 kJ mol-1. Overall, we show that a combination of thermodynamics, infrared spectroscopy and quantum chemistry provides quantification and detailed understanding of structure and molecular interaction in ibuprofen and related compounds.

How much different are thermochemical properties of enantiomers and their racemates? Thermochemical properties of enantiopure and racemate of methyl- and butyl lactates

This work is a contribution to the molecular understanding of the thermodynamic properties of the chiral compounds. A comprehensive thermochemical study of the liquid enantiopure and racemate pairs of optically active alkyl lactates has been performed. Vapor pressures of DL-(±)-, L-(-)-methyl-, and DL-(±)-, L-(-)-n-butyl esters of lactic acid were measured by the transpiration method. The liquid phase standard molar enthalpies of formation of these esters were measured by using the high-precision combustion calorimetry. The standard molar enthalpies of vaporization of alkyl lactates at 298.15 K were derived from vapor pressure temperature dependencies. Thermochemical data of these compounds were collected, evaluated, and tested for internal and external consistency. The high-level G4 quantum-chemical method was used for mutual validation of the experimental and theoretical gas phase enthalpies of formation of alkyl lactates. A critical review of the available thermochemical data for the liquid and crystalline enantiopure and racemate pairs of optically active compounds has been performed. Useful general trends in energetics of sublimation, vaporization, and formation of optically active compounds have been revealed. This knowledge is required for evaluation of new and already available experimental data for the chiral compounds, and it can be helpful to assess volatility or feasibility of processes to separate enantiomers.

Impact of chirality on peculiar ibuprofen molecular dynamics: hydrogen bonding organization and syn vs. anti carboxylic group conformations

Impact of chirality (R and S enantiomers) on syn vs. anti carboxylic group conformations, hydrogen bond dimers and peculiar ibuprofen molecular dynamics.

Effects of crystal habit on the sticking propensity of ibuprofen-A case study

This study demonstrates the effect of active pharmaceutical ingredient (API) particle habit on the sticking propensity of ibuprofen. Four diverse crystal habits with similar physico chemical properties are reported and the sticking propensity was found to increase with shape regularity. The surface energy of the extreme habits were shown to be different where particles that were more regular in shape exhibited surface energies of 9mJ/m² higher than those that were needle-like in habit. Computational and experimental data reveals that the increase in surface energy of the regular shaped particles can be attributed to the increase in the specific (polar) component, which is due to greater presence of faces which contain the carboxylic acid functionality at the surface. The increase in the specific energy component is shown to correlate with the sticking propensity of ibuprofen. It is proposed that investigation of the chemical causality of sticking, for this API and others, using the techniques demonstrated in this paper will be of increasing importance.

Are the Sublimation Thermodynamics of Organic Molecules Predictable?

We compare a range of computational methods for the prediction of sublimation thermodynamics (enthalpy, entropy, and free energy of sublimation). These include a model from theoretical chemistry that utilizes crystal lattice energy minimization (with the DMACRYS program) and quantitative structure property relationship (QSPR) models generated by both machine learning (random forest and support vector machines) and regression (partial least squares) methods. Using these methods we investigate the predictability of the enthalpy, entropy and free energy of sublimation, with consideration of whether such a method may be able to improve solubility prediction schemes. Previous work has suggested that the major source of error in solubility prediction schemes involving a thermodynamic cycle via the solid state is in the modeling of the free energy change away from the solid state. Yet contrary to this conclusion other work has found that the inclusion of terms such as the enthalpy of sublimation in QSPR methods does not improve the predictions of solubility. We suggest the use of theoretical chemistry terms, detailed explicitly in the Methods section, as descriptors for the prediction of the enthalpy and free energy of sublimation. A data set of 158 molecules with experimental sublimation thermodynamics values and some CSD refcodes has been collected from the literature and is provided with their original source references.

Intermolecular and low-frequency intramolecular Raman scattering study of racemic ibuprofen

We report the low-temperature Raman scattering study of racemic ibuprofen. Detailed analysis of the racemic ibuprofen crystal symmetry, related to the vibrational properties of the system, has been presented. The first principle calculations of a single ibuprofen molecule dynamical properties are compered with experimental data. Nineteen, out of 26 modes expected for the spectral region below 200cm(-1), have been observed.

Overall stability for the ibuprofen racemate: experimental and topological results leading to the pressure-temperature phase relationships between its racemate and conglomerate

Enantiomer resolution is much sought after for pharmaceutical applications, because many optically active drug molecules have only one pharmaceutically active enantiomer. Although it is always possible to force separation, it will come at a cost. The present method, based on thermodynamics, provides a relatively easy approach to investigate whether separation can be thermodynamically spontaneous. A topological phase diagram of the binary enantiomer system at 0.5 mol-fraction is constructed as a function of temperature and pressure after analysis of pressure and heat related quantities. It is demonstrated that for ibuprofen, an optically active analgesic, the racemate is the only stable solid form; the phase relationship between the racemate and the conglomerate is analogous to dimorphism with overall monotropy in pure chemical compounds.

Molecular structure and vibrational spectra of ibuprofen using density function theory calculations

The molecular geometry and the theoretical harmonic frequencies and infrared intensities of ibuprofen were calculated for all the molecules using five different density functional methods (mPW1PW91, B3PW91, B3LYP, HCTH and LSDA) with five basic sets, including 6-311G, 6-311++G, 6-311+G (d, p), 6-311++G (d, p) and 6-311++G (2d, 2p). The purpose of this research was to compare the performance of different DFT methods at different basis sets in predicting geometry and vibration spectrum of ibuprofen. The optimized geometric band lengths and bond angles obtained by using mPW1PW91 at 6-311++G (d, p) and 6-311++G (2d, 2p) basic sets show the best agreement with the experimental data. Comparison of the observed fundamental vibrational frequencies of ibuprofen with calculated results indicates that the B3PW91/6-311++G (2d, 2p) level is superior to all the remaining levels for predicting all the vibration spectra on average for ibuprofen.

An examination of the thermodynamics of fusion, vaporization, and sublimation of ibuprofen and naproxen by correlation gas chromatography

The vaporization enthalpies of (S)-ibuprofen and (S)-naproxen measured by correlation gas chromatography at T = 298.15 K are reported and compared with literature values. Adjustment of the fusion enthalpies of (RS)- and (S)-ibuprofen and (S)-naproxen to T = 298.15 K and combined with the vaporization enthalpy of the (S)-enantiomer of both ibuprofen and naproxen also at T = 298.15 K resulted in the sublimation enthalpies of both (S)-enantiomers. On the assumption that the vaporization enthalpy of the racemic form of ibuprofen is within the experimental uncertainty of the chiral form, the sublimation enthalpy of racemic ibuprofen was also evaluated. The vaporization and sublimation enthalpies compare favorably to the most of the literature values for the racemic form of ibuprofen but differ from the value reported for chiral ibuprofen. The literature values of (S)-naproxen are somewhat smaller than the values measured in this work. The following vaporization enthalpies were measured for (S)-ibuprofen and (S)-naproxen, respectively: ΔH(vap) (298.15 K), 106.0 ± 5.5, 132.2 ± 5.0 kJ·mol(-1) . Sublimation enthalpies of 122.7 ± 5.6 and 155.2 ± 7.1 kJ·mol(-1) were calculated for the (S)-enantiomers of ibuprofen and naproxen and a value of 128.9 ± 5.8 kJ·mol(-1) was estimated for the racemic form of ibuprofen.

Temperature effect on the vibrational dynamics of cyclodextrin inclusion complexes: investigation by FTIR-ATR spectroscopy and numerical simulation

The vibrational dynamics of solid inclusion complexes of the nonsteroidal anti-inflammatory drug Ibuprofen (IBP) with beta-cyclodextrin (beta-CD) and methyl-beta-cyclodextrin (Me-beta-CD) has been investigated by using attenuated total reflection-Fourier transform infrared FTIR-ATR spectroscopy, in order to monitor the changes induced, as a consequence of complexation, on the vibrational spectrum of IBP, in the wavenumber range 600-4000 cm(-1). Quantum chemical calculations were performed on monomeric and dimeric structures of IBP, derived from symmetric hydrogen bonding of the two carboxylic groups, in order to unambiguously assign some characteristic IR bands in the IBP spectrum. The evolution in temperature from 250 to 340 K of the C horizontal lineO stretching vibration, described by a best-fit procedure, allowed us to extract the thermodynamic parameter DeltaH associated to the binding of IBP with betaCDs in the solid phase. By comparing these results, Me-beta-CD has been shown to be the most effective carrier for IBP.

Thermodynamic study of the solubility of ibuprofen in acetone and dichloromethane

Thermodynamic functions, Gibbs energy, enthalpy and entropy for the solution processes of ibuprofen (IBP) in acetone and dichloromethane (DCM) were calculated from solubility values obtained at temperatures ranging from 293.15 K to 313.15 K. The respective thermodynamic functions for mixing and solvation processes as well as the activity coefficients for the solute were calculated. IBP solubility was high and proved similar in both solvents but was greater in DCM than acetone. In addition, the thermodynamic quantities for the transfer process of this drug from cyclohexane to the organic solvents were also calculated in order to estimate the contributions of hydrogen-bonds or of other dipolar interactions. The results were discussed in terms of solute-solvent interactions.

The SAMPL2 blind prediction challenge: introduction and overview

The interactions between a molecule and the aqueous environment underpin any process that occurs in solution, from simple chemical reactions to protein-ligand binding to protein aggregation. Fundamental measures of the interaction between molecule and aqueous phase, such as the transfer energy between gas phase and water or the energetic difference between two tautomers of a molecule in solution, remain nontrivial to predict accurately using current computational methods. SAMPL2 represents the third annual blind prediction of transfer energies, and the first time tautomer ratios were included in the challenge. Over 60 sets of predictions were submitted, and each participant also attempted to estimate the error in their predictions, a task that proved difficult for most. The results of this blind assessment of the state of the field for transfer energy and tautomer ratio prediction both indicate where the field is performing well and point out flaws in current methods.

Initial salt screening procedures for manufacturing ibuprofen

The aim of this paper is to design initial salt screening procedures for manufacturing ibuprofen. Salt forms of a pharmaceutical acid racemic (R,S)-(+/-)-ibuprofen and their "developable" synthetic routes were ferreted out simultaneously through the screening of seven bases of sodium hydroxide, potassium hydroxide, L-arginine, L-histidine, L-lysine, diethanolamine, and tris(hydroxymethyl)aminomethane (THAM), and the match with the use of nine organic solvents of methanol, dimethyl sulfoxide, ethanol, N, N-dimethylformamide, acetonitrile, isopropyl alcohol, 1,4-dioxane, acetone, and tetrahydrofuran mainly in the presence of water in 20 mL scintillation vials. Racemic (R,S)-(+/-)-sodium ibuprofen dihydrate, a well-known ibuprofen salt and the newly discovered racemic (R,S)-(+/-)-THAM ibuprofen, appeared as white-squared powders with a molecular weight of 327.42 g/mol, a melting point of 160.17 degrees C, and the apparent solubility product, K'(sp), of 6.0 x 10(-4) M(2) at 25 degrees C were successfully synthesized by the initial salt screening methods. The new amine salt of ibuprofen was monoclinic and had a space group of P2(1)/c and lattice parameters of a = 17.578(8) degrees, b = 10.428(4) degrees, c = 9.991(4) A, alpha = 90.00 degrees , beta = 97.17(1) degrees, gamma = 90.00 degrees, and V = 1,817.05(244) A(3). The aspect ratio of the amine salt crystals of ibuprofen of approximately 1.0 implied that the crystals had a better flowability than the sodium salt counterparts. This amine salt of ibuprofen was more stable in moist or dried atmospheres and was more hydrophobic than the sodium salt of ibuprofen. Moreover, the slow dissolution of this amine salt of ibuprofen might have made it less bitter and more suitable as a sustained release drug than the sodium salt of ibuprofen. The future work is to search for the different polymorphs of this amine salt of ibuprofen and to extend the initial salt screening working logics to the formation of co-crystals.

Physicochemical properties of zwitterionic L- and DL-alanine crystals from their experimental and theoretical charge densities

The total experimental electron density distributions rho(r) of zwitterionic L- and DL-alanine crystals, as derived from extensive sets of X-ray diffracted intensities collected at 23 and 19 K, are compared to gain an insight into the different physical properties of the two related chiral compounds in the solid state and to explore the extent of the rho(r) transferability. Relevant parameters that characterize the two crystal forms are obtained, showing differences and similarities in terms of (i) geometric descriptors, (ii) topological indexes, (iii) molecular electrostatic potential Phi(r) distributions, (iv) atomic volumes and charges, (v) molecular electric moments, and (vi) electrostatic interaction energies. To assess the relative stability of the racemate with respect to the pure enantiomer, the crystal lattice energies, as obtained through DFT fully periodic calculations, are also discussed and compared with the experimental sublimation enthalpies after correction for the proton-transfer energies. In-crystal group charges, evaluated with the quantum theory of atoms in molecules, are found to be transferable between the racemic and the pure enantiomer, at variance with group volumes. Similarly, molecular first and third moments are not strictly transferable and indicate that for the zwitterionic alanine molecule the molecular charge distribution in the DL-crystal is more polarized in the c direction by about 10%. By contrast, quantitative agreement is observed for second and fourth moments. Significant differences arise from (1) the crystal packing of the dipole vectors, which are aligned in an antiparallel fashion in the L-crystal, to be compared with a parallel alignment in the racemate, due the polar space group Pna21 of the latter, (2) the strongly attractive electrostatic energy of a homochiral pair in the L-crystal, which is opposed to the corresponding heterochiral pair in the DL-crystal form. The difference between these Ees values amounts to 135-150 kJ mol(-1). Despite this, the two crystal forms are predicted as equally thermodynamically favored by the theoretical P-B3LYP estimates of the crystal lattice energies. Finally, the necessity of an upgrading of the dispersion and exchange-repulsion terms currently adopted within the experimental charge density approach to intermolecular interactions is recognized and discussed.

Predicting intrinsic aqueous solubility by a thermodynamic cycle

We report methods to predict the intrinsic aqueous solubility of crystalline organic molecules from two different thermodynamic cycles. We find that direct computation of solubility, via ab initio calculation of thermodynamic quantities at an affordable level of theory, cannot deliver the required accuracy. Therefore, we have turned to a mixture of direct computation and informatics, using the calculated thermodynamic properties, along with a few other key descriptors, in regression models. The prediction of log intrinsic solubility (referred to mol/L) by a three-variable linear regression equation gave r(2)=0.77 and RMSE=0.71 for an external test set comprising drug molecules. The model includes a calculated crystal lattice energy which provides a computational method to account for the interactions in the solid state. We suggest that it is not necessary to know the polymorphic form prior to prediction. Furthermore, the method developed here may be applicable to other solid-state systems such as salts or cocrystals.

In Silico Prediction of Drug Solubility. 3. Free Energy of Solvation in Pure Amorphous Matter

The solubility of drugs in water is investigated in a series of papers. In this work, we address the process of bringing a drug molecule from the vapor into a pure drug amorphous phase. This step enables us to actually calculate the solubility of amorphous drugs in water. In our general approach, we, on one hand, perform rigorous free energy simulations using a combination of the free energy perturbation and thermodynamic integration methods. On the other hand, we develop an approximate theory containing parameters that are easily accessible from conventional Monte Carlo simulations, thereby reducing the computation time significantly. In the theory for solvation, we assume that DeltaG* = DeltaGcav + ELJ + EC/2, where the free energy of cavity formation, DeltaGcav, in pure drug systems is obtained using a theory for hard-oblate spheroids, and ELJ and EC are the Lennard-Jones and Coulomb interaction energies between the chosen molecule and the others in the fluid. The theoretical predictions for the free energy of solvation in pure amorphous matter are in good agreement with free energy simulation data for 46 different drug molecules. These results together with our previous studies support our theoretical approach. By using our previous data for the free energy of hydration, we compute the total free energy change of bringing a molecule from the amorphous phase into water. We obtain good agreement between the theory and simulations. It should be noted that to obtain accurate results for the total process, high precision data are needed for the individual subprocesses. Finally, for eight different substances, we compare the experimental amorphous and crystalline solubility in water with the results obtained by the proposed theory with reasonable success.