In a step towards understanding the structure–property relationship among Synthetic Cathinones (SCs), a combined methodology based on Density Functional Theory (DFT), Administration, Distribution, Metabolism, Excretion, and Toxicity (ADMET) predictions, docking and molecular dynamics simulations have been applied to correlate physicochemical descriptors of various SCs to their biological activity. The results from DFT and molecular docking studies correlate well with each other explaining the biological activity trends of the studied SCs. Quantum mechanical descriptors viz. polarizability, electron affinity, ionization potential, chemical hardness, electronegativity, molecular electrostatic potential, and ion interaction studies unravel the distinguishingly reactive nature of Group D (pyrrolidine substituted) and Group E (methylenedioxy and pyrrolidine substituted) compounds. According to ADMET analysis, Group D and Group E molecules have a higher probability of permeating through the blood–brain barrier. Molecular docking results indicate that Phe76, Ala77, Asp79, Val152, Tyr156, Phe320, and Phe326 constitute the binding pocket residues of hDAT in which the most active ligands MDPV, MDPBP, and MDPPP are bound. Finally, to validate the derived quantum chemical descriptors and docking results, Molecular Dynamics (MD) simulations are performed with homology-modelled hDAT (human dopamine transporter). The MD simulation results revealed that the majority of SCs remain stable within the hDAT protein’s active sites via non-bonded interactions after 100 ns long simulations. The findings from DFT, ADMET analysis, molecular docking, and molecular dynamics simulation studies complement each other suggesting that pyrrolidine-substituted SCs (Group D and E), specifically, MPBP and PVN are proven potent SCs along with MDPV, validating various experimental observations. Communicated by Ramaswamy H. Sarma
Read full abstract