AbstractNanobody binding stabilizes the active state of G-protein-coupled receptor (GPCR) and modulates its affinity for bound ligands. However, the atomic level basis for this allosteric regulation remains elusive. Here, we investigate the conformational changes induced by the binding of a nanobody (Nb80) on the active-like β2 adrenergic receptor (β2AR) via enhanced sampling molecular dynamics simulations. Dimensionality reduction analysis shows that Nb80 stabilizes a highly active state of the β2AR with a ~14 Å outward movement of transmembrane helix 6 and close proximity of transmembrane (TM) helices 5 and 7. This is further supported by the residues located at hotspots located on TMs 5, 6 and 7, as shown by supervised machine learning methods. Besides, ligand-specific subtle differences in the conformations assumed by intercellular loop 2 and extracellular loop 2 are captured from the trajectories of various ligand-bound receptors in the presence of Nb80. Dynamic network analysis further reveals that Nb80 binding can enhance the communication between the binding sites of Nb80 and of the ligand. We identify unique allosteric signal transmission mechanisms between the Nb80-binding site and the extracellular domains in presence of full agonist and G-protein biased partial agonist, respectively. Altogether, our results provide insights into the effect of intracellular binding partners on the GPCR activation mechanism, which could be useful for structure-based drug discovery.TOCGraphical Table of Contents