Insights into pharmaceutical nanocrystal dissolution: A molecular dynamics simulation study on aspirin

The presented molecular dynamics simulations are the first simulations to reveal dynamic dissolution of a pharmaceutical crystal in its experimentally determined shape. Continuous dissolution at constant undersaturation of the surrounding medium is ensured by introducing a plane of sticky dummy atoms into the water slab. These atoms have a strong interaction potential with dissolved aspirin molecules, but interactions with water are excluded from the calculations. Thus, the number of aspirin molecules diffusing freely in solution is kept at a low value and continuous dissolution of the aspirin crystal is monitored. Further insight into face-specific dissolution is drawn. The dissolution mechanism of receding edges is found for the (001) plane. These findings are in good agreement with experimental results. While the proposed dissolution mechanism for the (100) plane is terrace sinking on a rough surface, no pronounced dissolution of the perfectly flat face is seen in the present work. Molecular simulations of pharmaceuticals in their experimentally obtained structure therefore have shown to be especially suited for the investigation of dissolving faces, where the edges have a pronounced effect. In contrast to previous studies a propagation of the dissolution front into the crystal face is reported, and the crystal bulk is stable over the whole simulation time of 150 ns.