Abstract
Effective rational drug discovery targeting a specific protein hinges on understanding their functional states and distinguishing it from homologs. However, for the G protein coupled receptors, both the activation-related conformational changes (ACCs) and the intrinsic divergence among receptors can be misled or obscured by ligand-induced conformational changes (LCCs). Here, we unraveled ACCs and intrinsic divergence from LCCs of the dopamine D3 and D2 receptors (D3R and D2R), by analyzing their experimentally determined structures and the molecular dynamics simulation results of the receptors bound with different ligands. In addition to the ACCs common to other aminergic receptors, we revealed unique ACCs for these two receptors including TM5e and TM6e shifting away from TM2e and TM3e, with a subtle rotation of TM5e. In identifying intrinsic divergence, we found pronounced outward tilting of TM6e in the D2R compared to the D3R in both experimental structures and simulations with ligands in different scaffolds. This tilting was drastically reduced in the simulations of the receptors bound with nonselective full agonist quinpirole, suggesting a misleading impact of LCCs. Further, in the quinpirole-bound simulations, TM1 showed a greater disparity between these receptors, indicating that LCCs may obscure intrinsic divergence. In addition, our analysis showed that the impact of the nonconserved TM1 propagated to conserved Trp(7.40) and Glu(2.65), both are ligand binding residues. We also found that the D2R exhibited heightened flexibility compared to the D3R in the extracellular portions of TMs 5, 6, and 7, potentially associated with its greater ligand binding site plasticity. Our results lay the groundwork for crafting ligands specifically targeting D2R or D3R with more precise pharmacological profiles.