No abstract

No abstract

As a first step in the computational prediction of drug solubility the free energy of hydration, DeltaG*(vw) in TIP4P water has been computed for a data set of 48 drug molecules using the free energy of perturbation method and the optimized potential for liquid simulations all-atom force field. The simulations were performed in two steps, where first the Coulomb and then the Lennard-Jones interactions between the solute and the water molecules were scaled down from full to zero strength to provide physical understanding and simpler predictive models. The results have been interpreted using a theory assuming DeltaG*(vw) = A(MS)gamma + E(LJ) + E(C)/2 where A(MS) is the molecular surface area, gamma is the water-vapor surface tension, and E(LJ) and E(C) are the solute-water Lennard-Jones and Coulomb interaction energies, respectively. It was found that by a proper definition of the molecular surface area our results as well as several results from the literature were found to be in quantitative agreement using the macroscopic surface tension of TIP4P water. This is in contrast to the surface tension for water around a spherical cavity that previously has been shown to be dependent on the size of the cavity up to a radius of approximately 1 nm. The step of scaling down the electrostatic interaction can be represented by linear response theory.

In the present paper, we have studied particle dissolution and crystal growth of the poorly water soluble drug felodipine, using fluorescence as a probe for the amount of crystalline material. Dissolution kinetics is essentially diffusion-controlled, while the rate of crystal growth is significantly slower compared to the diffusion-controlled limit. The deviation from diffusion control was characterized by the effective length, lambda, related to the kinetics of a surface integration process. Amorphous nanoparticles may be highly unstable in the presence of small amounts of crystalline particles. This is due to the fact that the molecular solubility from the amorphous nanoparticles often is at least an order of magnitude higher than the corresponding crystalline solubility. In a mixed system where crystalline nanoparticles have been added to an amorphous nanosuspension, the bulk will have a monomer concentration intermediate between the amorphous and crystalline solubilities, and is thus supersaturated with respect to the crystalline particles while being undersaturated with respect to the amorphous particles. As a consequence, the amorphous particles spontaneously dissolve, while crystalline particles grow, in a combined process which is similar to Ostwald ripening. By knowing the parameters describing dissolution and crystal growth, respectively, it was possible to simulate the outcome of controlled seeding experiments, where a small amount of crystalline nanoparticles was added to a dispersion of amorphous nanoparticles. A good agreement between model calculations and experiments was obtained including how the crystal growth rate varied with the amounts of added crystalline seeds.

The solubility of drugs in water is investigated in a series of papers and in the current work. The free energy of solvation, DeltaG*(vl), of a drug molecule in its pure drug melt at 673.15 K (400 degrees C) has been obtained for 46 drug molecules using the free energy perturbation method. The simulations were performed in two steps where first the Coulomb and then the Lennard-Jones interactions were scaled down from full to no interaction. The results have been interpreted using a theory assuming that DeltaG*(vl) = DeltaG(cav) + E(LJ) + E(C)/2 where the free energy of cavity formation, DeltaG(cav), in these pure drug systems was obtained using hard body theories, and E(LJ) and E(C) are the Lennard-Jones and Coulomb interaction energies, respectively, of one molecule with the other ones. Since the main parameter in hard body theories is the volume fraction, an equation of state approach was used to estimate the molecular volume. Promising results were obtained using a theory for hard oblates, in which the oblate axial ratio was calculated from the molecular surface area and volume obtained from simulations. The Coulomb term, E(C)/2, is half of the Coulomb energy in accord with linear response, which showed good agreement with our simulation results. In comparison with our previous results on free energy of hydration, the Coulomb interactions in pure drug systems are weaker, and the van der Waals interactions play a more important role.