Structural context of GnRHR background mutations.

A homology model of M. musculus GnRHR (mGnRHR) and a crystal structure of H. sapiens GnRHR (hGnRHR) provide structural context for two mutants which disrupt folding. A) The homology model of mGnRHR (green) structurally aligns with the crystal structure of hGnRHR (gray). B) A side view of hGnRHR shows that T277, which is analogous to V276 in mGnRHR, is located in transmembrane domain 6 and is exposed to the lipid bilayer. C) A top view of hGnRHR shows that W107 in extracellular loop 1 forms a hydrogen bond network with two other residues.

Plasma Membrane Expression of mGnRHR variants.

Surface immunostaining and flow cytometry were used to compare the plasma membrane expression of select mGnRHR variants expressed in HEK293T cells. A histogram from one representative biological replicate depicts the distribution of plasma membrane expression intensities among cells expressing WT (black), V276T (red), or W107A (blue) mGnRHR.

Plasma membrane expression of mGnRHR cellular libraries.

mGnRHR mutant libraries (red), as well as their respective controls (WT or V276T/W107A single mutants, black) were expressed in mammalian cells, and the amount of mGnRHR expressed at the plasma membrane was measured during fluorescence-activated cell sorting. Histograms show the distributions of plasma membrane expression for mGnRHR mutants in the (A) WT library, (B) V276T library, and (C) W107A library.

Relative plasma membrane expression of mGnRHR mutants.

A) A histogram shows the distribution of plasma membrane expression measurements relative to their respective controls (WT or V276T/W107A single mutants) for the WT (gray), V276T (red) and W107A (blue) mutational libraries. A dashed line at the normalized plasma membrane expression value of 1.0, representing no mutational effect relative to the control, is displayed for reference. B & C) The average relative intensities across each domain are compared between the V276T (B, red) and W107A (C, blue) double mutant libraries and the single mutant library (gray). The size of each dot represents the number of variants measured, and significant p-values are also displayed. D – F) Average relative intensities for the (D) WT library, (E) V276T library, and (F) W107A library were projected onto a homology model of mGnRHR. Residues 276 and 107 are displayed in green for reference. Individual variant scores can be found in Table S4.

Epistasis in the mGnRHR double mutant libraries.

A histogram depicts the distribution of epistasis scores (ε) for interactions the subset of 251 high quality secondary mutations form with V276T (red) and W107A (blue).

General trends in mGnRHR epistasis.

Trends associated with the observed pairwise epistasis within mGnRHR are identified using unsupervised learning. A) Uniform manifold approximation projection (UMAP) was used to differentiate variants based on differences in their relative expression in the V276T, W107A, or WT background. Variants are projected onto an arbitrary two-dimensional coordinate based on the results and are colored according to whether they were assigned to cluster 1 (green), cluster 2 (purple), cluster 3 (orange), or were designated as outliers (gray) by HDBSCAN. The percentage of the mutations that fall within TMDs or loops are shown for reference. B) A box and whisker plot depicts the statistical distributions of relative PME values among variants within each cluster in the context of V276T (red), W107A (blue), or WT (Gray) mGnRHR. Select clusters of variants that exhibit statistically different expression profiles according to a Mann-Whitney U-test are indicated (*, p < 0.001). A value of 1 corresponds to mutations that have no effect on the PME of mGnRHR in the indicated genetic background. C) A box and whisker plot depicts the statistical distribution of epistasis scores associated with the interactions between the mutations within each cluster and either V276T (red) or W107A (blue). P-values for select Mann-Whitney U-tests comparing the interactions of these mutations with V276T and W107A are indicated. A value of 0 indicates that the effects of the two mutations are additive. D) A box and whisker plot depicts the statistical distribution of Rosetta ΔΔG values among mutations within each cluster. P-values for select Mann-Whitney U-tests comparing the ΔΔG values across clusters are shown for reference. For panels B-D, the edges of the boxes correspond to the 75th and 25th percentile values while the whiskers reflect the values of the 90th and 10th percentile. The central hash and square within the box represent the average and median values, respectively. These analyses were carried out on a subset of 243 variants with high-quality expression measurements with calculable Rosetta ΔΔG values.

Context-dependent impacts of stabilizing and destabilizing mutations.

Line plots depict the relationship between the equilibrium fraction of folded protein (ffold) and the free energy of unfolding (ΔGunf). A) A moderately stable protein (ΔGunf = +2 kcal/mol) exhibits a greater change in ffold after acquiring one or more destabilizing mutations (total ΔΔGunf = -2 kcal/mol, red) than after acquiring one or more mutations that stabilize the protein to the same extent (total ΔΔGunf = +2 kcal/mol, blue). B) A moderately unstable protein (ΔGunf = -2 kcal/mol) exhibits a greater change in ffold after acquiring one or more stabilizing mutations (total ΔΔGunf = +2 kcal/mol, blue) than after acquiring one or more mutations that destabilize the protein to the same extent (total ΔΔGunf = -2 kcal/mol, red).