In this work, the polymorphism of the active pharmaceutical ingredient carbamazepine (CBZ) was investigated by using molecular dynamics simulations with an enhanced sampling scheme. A single molecule of CBZ attaching to flat surfaces of different polymorphs was used as a model for secondary nucleation in water. A novel approach was developed to compute the free energy profile characterizing the adsorption of molecules with orientation aligned with the crystal structure of the surface. The distribution of states that showed alignment was used to rescale the adsorption free energy to include only the contribution that is consistent with crystal growth. The resulting free energy surfaces showed favorable thermodynamics for the most stable form, Form III and the second most stable form, Form I. The primary crystallization product, a dihydrate, was found to be less favorable, implying a nonclassical crystallization pathway. We suggest that a major contribution determining the energetics is the hydrophobicity of the surface. This thermodynamic ranking provides valuable information about the molecular pathways of polymorph growth and will further contribute to the understanding of the crystallization process of CBZ, which is imperative since polymorph formation can alter the physical properties of a drug significantly.
Read full abstract