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 (_adjR² = 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.

Voltage-gated sodium channels (VGSCs) are heteromeric protein complexes that initiate action potentials in excitable cells. The voltage-gated sodium channel accessory subunit, Navβ1, allosterically modulates the α subunit pore structure upon binding. To date, the molecular determinants of the interface remain unknown. We made use of sequence, knowledge and structure-based methods to identify residues critical to the association of the α and β1 Nav1.4 subunits. The Navβ1 point mutant C43A disrupted the modulation of voltage dependence of activation and inactivation and delayed the peak current decay, the recovery from inactivation, and induced a use-dependent decay upon depolarisation at 1 Hz. The Navβ1 mutant R89A selectively delayed channel inactivation and recovery from inactivation and had no effect on voltage dependence or repetitive depolarisations. Navβ1 mutants Y32A and G33M selectively modified the half voltage of inactivation without altering the kinetics. Despite low sequence identity, highly conserved structural elements were identified. Our models were consistent with published data and may help relate pathologies associated with VGSCs to the Navβ1 subunit.