Induced fit, ensemble binding space docking and Monte Carlo simulations of MDMA ‘ecstasy’ and 3D pharmacophore design of MDMA derivatives on the human serotonin transporter (hSERT)

The popular recreational drug MDMA (3,4-methylenedioxy-methamphetamine) has a documented potential as a psychopharmacological clinical and research tool. This is due to its unique ability to promote reprocessing of traumatic memories, empathetic and pro-social states. Although it is established that MDMA exerts its behavioural effects via the serotonin transporter (SERT), the ligand-protein molecular interplay remains elusive. In order to shed light on the binding of MDMA and its primary congeneric entactogens (MDA, MBDB and MDAI), we first combined induced fit with Monte Carlo simulations. The computed interaction energies of the models correlated well with experimental activities ((adj)R(2) = 0.78). Then we carried out ‘ensemble binding space docking’ on trajectories generated by interpolation of experimentally derived structures of the hSERT from the outward-open, and the occluded, to the inward-open states. This approach revealed low-energy alternative binding modes, suggesting high occupancy of the central site, yet considerable MDMA mobility within it, favouring the paroxetine-like orientation. Finally, we designed a pharmacophore that may be used to recognise hSERT-mediated serotonin releasers and uptake inhibitors of diverse chemical structure, identifying their active conformations and interacting residues. We conclude that the conserved amine-Asp98 ionic and edge-to-face π-π interactions are crucial to the mode of action of MDMA on the hSERT, underscoring the contributions of Tyr95 and gating residues Phe341, Tyr176 and Phe335. Amenable to experimental testing, our modelling may aid the rational design of novel entactogenic compounds and contribute to the understanding of an action mechanism, common and typical of psychotropic agents.