Allosteric pathways govern Gα protein coupling selectivity at promiscuous GPCRs

G protein-coupled receptors (GPCRs) regulate diverse physiological responses by engaging distinct heterotrimeric G proteins, yet the basis of G selectivity in promiscuous receptors remains unclear. Although ligand bias holds therapeutic promise, selectivity has been assumed to reside mainly in the ligand-binding site (LBS) or G protein interface (GPI). Here, we combine whole-receptor mutagenesis, functional G protein assays, molecular dynamics simulations, interpretable machine learning method, and Bayesian network modeling to identify residue networks governing Gq/11 and G12/13 coupling and to define the molecular basis of G protein preference and promiscuity at two vasopressor GPCRs, the angiotensin II type 1 and prostaglandin F2 receptors. While residues within the LBS and GPI domains contribute to coupling efficiency and subtype discrimination, we find that long-range allosteric communication across the receptor, including from structurally unresolved domains, is the principal determinant of G protein preference and promiscuity. These allosteric pathways integrate multiple receptor domains, confer signaling robustness to mutation, and hierarchically govern coupling preferences. Our findings suggest that G protein selectivity is an allosterically encoded property of GPCRs and provide a conceptual framework for designing ligands and receptors with tailored G protein-biased signaling.