Dopamine D₁ receptor and serotonin 5-HT_1A receptor agonist effects of the natural product (–)-stepholidine: molecular modelling and dynamics simulations

Abstract

In this study, by homology modelling and molecular dynamics (MD) simulation, models of l-stepholidine (l-SPD) activating the 5-HT_1A and D₁ receptors were constructed. In 100-ns MD simulations, the D₁ and 5-HT_1A receptors were activated by the partial agonist l-SPD, conforming with the global toggle switch activation model and the sequential activation model. The residues Y^7.53 and Y^5.58 swing significantly between different transmembrane (TM) domains after activation. Similarities between D₁ and 5-HT_1A included (1) the outward motion of TM-5; (2) the ionic lock was independent of the tilt of TM-6 and (3) there was an apparent bending of TM-6, and the ring of l-SPD formed strong π–π interactions with residue W^6.48. Differences between the two included the following: (1) in 5-HT_1A, l-SPD formed a hydrogen bond with Ala172^5.46 of TM-5, and the intracellular end of TM-5 moved outward slowly; that hydrogen bond did not form with the D₁ receptor; (2) l-SPD formed stronger interactions with D^3.32 and W^6.48 in the D₁ receptor than in the 5-HT_1A receptor and (3) the hydrogen bonding network was somewhat different in SPD-5-HT_1A and SPD-D₁ receptors. We propose the interaction between l-SPD and D^3.32 or/and W^6.48 is the original driving force during the whole activation process.

Keywords:

• molecular dynamics
• l-SPD (l-stepholidine)
• G-protein-coupled receptor (GPCR)
• D₁ receptor
• 5-HT_1A receptor