[Prediction of the Crystal Growth Mechanism of Aspirin Using Molecular Simulations]

Controlling the physicochemical properties of a drug formulation is important for proper drug efficacy, since in the gastrointestinal tract many drugs undergo dissolution, limiting their efficacy. Factors affecting a drug's physicochemical properties include its crystal habit. Therefore, we predicted the crystal habit by molecular simulation for the purpose of controlling crystal morphology. In this study, we used aspirin as a model compound. By performing simulations based on known crystal structure data, we trained the simulation algorithm to produce the cubic and plate-like morphologies of crystals actually obtained. By these methods, we showed that the crystal plane of the crystal form actually obtained coincides with the characteristic crystal plane obtained by simulation. Furthermore, to consider the influence of the crystallization solvent on crystal growth, we simulated adsorption of solvent molecules on characteristic crystal planes. The difference in adsorption energy of the solvent molecules prevents the aspirin molecules from attaching to the crystal plane. As a result, we concluded that the crystal habit was caused by the difference in growth rate of the crystal plane. By applying the methods developed in this research, the growth of crystal planes can be predicted by molecular simulation, making it possible to efficiently obtain crystal forms with optimal physical properties for drug development. We believe that further development of this approach will lead to dramatic decreases in the cost and duration of drug development.