Molecular Insights into the D1R Agonist and D2R/D3R Antagonist Effects of the Natural Product (−)-Stepholidine: Molecular Modeling and Dynamics Simulations

Abstract

(-)-Stepholidine (l-SPD), an active ingredient of the Chinese herb Stephania, is the first compound found to have a dual function as a dopamine receptor D1 agonist and D2 antagonist. The preliminary dynamical behaviors of D1R and D2R and their interaction modes with l-SPD were investigated in our previous study. Recently, the pharmacological effect of l-SPD on D3R was elucidated as an antagonist. This new discovery in combination with the explosion of structural biology in GPCR superfamily prompted us to perform a more comprehensive investigation on the special pharmacological profiles of l-SPD on dopamine receptors. In this study, the integration of homology modeling, automated molecular docking, and MD simulations was used to probe the agonistic and antagonistic mechanism of l-SPD on D1R, D2R, and D3R. Our analyses showed that hydrogen bonding of the hydroxyl group on the D ring of l-SPD with side chain of N6.55 which, in combination with hydrophobic stacking between I3.40, F6.44 and W6.48, is the key feature to mediate the agonist effect of l-SPD on D1R, whereas the absence of hydrophobic stacking between I3.40, F6.44, and W6.48 in D2R and D3R excludes receptor activation. Finally, the agonistic and antagonistic mechanisms of l-SPD and an activation model of D1R were proposed on the basis of these findings. The present study could guide future experimental works on these receptors and has the significance to the design of functionally selective drugs targeting dopamine receptors.

Tetrahydroprotoberberine alkaloids with dopamine and σ receptor affinity

“…40 As compared to prior work, we found similarities in the key residues that are important for affinity of stepholidine (and the C10 analogues), particularly Asp110, Ser192 and Phe345. Apart from the key protonated nitrogen-Asp110 salt bridge interaction, in both studies Ser192 interacts with an aromatic hydroxyl group via a hydrogen bond and Phe345 interacts with an aromatic ring via aromatic stacking (Fig 4a).…”

Computational investigation on the binding modes of Rimonabant analogs with CB1 and CB2

“…With the development of computer-aid drug discovery methods, the molecular modeling and dynamics simulation are becoming perfect complements in explaining experimental results and guiding further experiments. [16,17] The homology modeling is an important approach to get accurate 3D structures in the absence of experimental structures [18] ; the molecular docking is widely used in predicting the binding affinities of ligands and binding modes between ligands and receptors [19] ; the molecular dynamics simulations can provide detailed information about F I G U R E 1 Chemical structures of four selective inverse agonists/antagonists. Rimonabant, MJ08, and MJ15 are selective inverse agonists/antagonists of CB1, whereas AM630 is a selective inverse agonist/antagonist of CB2.…”

Improving prediction of helix–helix packing in membrane proteins using predicted contact numbers as restraints

“…HMPs are also very important therapeutic targets, about 50% of therapeutics on the market target HMPs . The availability of a three‐dimensional (3D) structure of a HMP not only improves our understanding of how the protein works at the atomic level but also facilitates the development of new therapeutics . Despite great progress in experimental techniques for determining HMP structures, only ∼2% structures in the protein databank are HMPs, highlighting the fact that HMP structure characterization is still a challenge.…”