Abstract
Development of effective bivalent ligands has become the focus of intensive research toward modulation of G protein-coupled receptor (GPCR) oligomers, particularly in the field of GPCR pharmacology. Experimental studies have shown that they increased binding affinity and signaling potency compared to their monovalent counterparts, yet underlying molecular mechanism remains elusive. To address this, we performed accelerated molecular dynamics simulations on bivalent-ligand bound Adenosine 2A receptor (A(2A)R) dimer in the context of a modeled tetramer, which consists of A(2A)R and dopamine 2 receptor (D(2)R) homodimers and their cognate G proteins. Our results demonstrate that bivalent ligand impacted interactions between pharmacophore groups and ligand binding residues, thus modulating allosteric communication network and water channel formed within the receptor. Moreover, it also strengthens contacts between receptor and G protein, by modulating the volume of ligand binding pocket and intracellular domain of the receptor. Importantly, we showed that impact evoked by the bivalent ligand on A(2A)R dimer was also transmitted to apo D(2)R, which is part of the neighboring D(2)R dimer. To the best of our knowledge, this is the first study that provides a mechanistic insight into the impact of a bivalent ligand on dynamics of a GPCR oligomer. Consequently, this will pave the way for development of effective ligands for modulation of GPCR oligomers and hence treatment of crucial diseases such as Parkinson's disease and cancer.