Agonist binding to G-protein-coupled receptors (GPCRs) triggers signal transduction cascades involving heterotrimeric G proteins as key players. A major obstacle for drug design is the limited knowledge of conformational changes upon agonist binding, the details of interaction with the different G proteins, and the transmission to movements within the G protein. Although a variety of different GPCR/G protein complex structures would be needed, the transient nature of this complex and the intrinsic instability against dissociation make this endeavor very challenging. We have previously evolved GPCR mutants that display higher stability and retain their interaction with G proteins. We aimed at finding all G-protein combinations that preferentially interact with neurotensin receptor 1 (NTR1) and our stabilized mutants. We first systematically analyzed by coimmunoprecipitation the capability of 120 different G-protein combinations consisting of αi1 or αsL and all possible βγ-dimers to form a heterotrimeric complex. This analysis revealed a surprisingly unrestricted ability of the G-protein subunits to form heterotrimeric complexes, including βγ-dimers previously thought to be nonexistent, except for combinations containing β5. A second screen on coupling preference of all G-protein heterotrimers to NTR1 wild type and a stabilized mutant indicated a preference for those Gαi1βγ combinations containing γ1 and γ11. Heterotrimeric G proteins, including combinations believed to be nonexistent, were purified, and complexes with the GPCR were prepared. Our results shed new light on the combinatorial diversity of G proteins and their coupling to GPCRs and open new approaches to improve the stability of GPCR/G-protein complexes.