In this article, density functional theory (DFT) calculations and 30 ns molecular dynamic (MD) simulations were performed to investigate the ability of α-, β- and γ-cyclodextrins (CDs) to form selective complexes with cathinone. DFT calculations in the gas phase, water, chloroform and methanol reveal that the solvents, reduce the stability of the complexes. Optimized structures confirm that α-CD cannot encapsulate cathinone, completely, while other CDs showed an opposite behavior. DFT calculations indicate that cathinone has the most stable complex with γ-CD in comparison to the α- and β-CDs. Natural bond orbital and quantum theory of atoms in molecules analyses reveal that the electrostatic interactions between cathinone and CDs are the driving force of the complex formation. MD simulations confirm that different solvents play an important role in the stability of the cathinone complexes and the obtained MD results are in good agreement with the DFT calculations.

Keywords