Francesco De Pascali
email: francesco.depascali@jefferson.edu

Mohammed Akli Ayoub
email: mayoub@uaeu.ac.ae

Riccardo Benevelli
email: riccardo.benevelli95@gmail.com

Silvia Sposini
email: silvia.sposini@u-bordeaux.fr

Jordan Lehoux
email: jordan.lehoux.chem@gmail.com

Nathalie Gallay

Pauline Raynaud
email: pauline.raynaud@inrae.fr

Frédéric Flavie Landomiel
email: flavie.landomiel@gmail.com

Christophe Gauthier
email: christophe.gauthier@inrae.fr

Francesco De Pasca

Silvia Sposi

Frédéric Jean-Alphonse

Lucie P Pellissier
email: lucie.pellissier@inrae.fr

Pascale Crépieux
email: pascale.crepieux@inrae.fr

Anne Poup
email: anne.poupon@inrae.fr

Asuka Inou

Nicolas Joubert
email: nicolas.joubert@univ-tours.fr

Marie-Claude Viaud-Massuard
email: marieclaude.viaudmassuard@mcsaf.fr

Livio Casarini
email: livio.casarini@unimore.it

Manuela Simoni
email: manuela.simoni@unimore.it

Aylin C Hanyaloglu
email: a.hanyaloglu@imperial.ac.uk

Selva G Nataraja

Henry N Yu

Stephen S Palmer
email: stephen.palmer2@bcm.edu

Romain Yvinec
email: romain.yvinec@inrae.fr

Eric Reiter
email: eric.reiter@inrae.fr

De Pascali

F ; Ayoub

M A ; Benevelli

R ; Sposini

S ; Lehoux

Pharmacological characterization of low molecular weight biased agonists at the follicle stimulating hormone receptor

Keywords: FSHR, biased signaling, allosteric ligands, system bias, operational model

à la diffusion de documents scientifiques de niveau recherche, publiés ou non, émanant des établissements d'enseignement et de recherche français ou étrangers, des laboratoires publics ou privés.

Introduction

Follicle-stimulating hormone receptor (FSHR), together with its cognate endogenous ligand follicle-stimulating hormone (FSH), plays a key role in regulating gametogenesis in both female and male reproduction [START_REF] Gromoll | Functional and Clinical Consequences of Mutations in the FSH Receptor[END_REF][START_REF] Kaprara | The Hypothalamus-Pituitary-Gonad Axis: Tales of Mice and Men[END_REF]. FSHR is a G protein-coupled receptor (GPCR) made of a large glycosylated extracellular domain (ECD) containing multiple leucine-rich repeats connected to seven spanning transmembrane (TMD) α-helices joined together by three extracellular loops and three intracellular loops and terminated by a C-terminal intracellular tail. The horse-shoe-shaped ECD is prototypical to the glycoprotein hormone receptor subfamily where the orthosteric binding site is located [START_REF] Fan | Structure of Human Follicle-Stimulating Hormone in Complex with Its Receptor[END_REF]. Once activated, FSHR induces a plethora of intracellular mechanisms [START_REF] Gloaguen | Mapping the Follicle-Stimulating Hormone-Induced Signaling Networks[END_REF][START_REF] Landomiel | Biased Signaling and Allosteric Modulation at the FSHR[END_REF], amongst which are the Gα s /PKA/cAMP [START_REF] Dattatreyamurty | Physical and Functional Association of Follitropin Receptors with Cholera Toxin-Sensitive Guanine Nucleotide-Binding Protein[END_REF][START_REF] Salvador | Follicle-Stimulating Hormone Stimulates Protein Kinase A-Mediated Histone H3 Phosphorylation and Acetylation Leading to Select Gene Activation in Ovarian Granulosa Cells[END_REF] and extracellular signal-regulated kinase (ERK) pathways [START_REF] Cottom | Follicle-Stimulating Hormone Activates Extracellular Signal-Regulated Kinase but Not Extracellular Signal-Regulated Kinase Kinase through a 100-KDa Phosphotyrosine Phosphatase[END_REF]. Like most GPCRs, β-arrestin recruitment and binding to FSHR lead to desensitization and internalization [START_REF] De Pascali | β-Arrestins and Biased Signaling in Gonadotropin Receptors[END_REF] while also assembling signaling modules [START_REF] Kara | A Phosphorylation Cluster of Five Serine and Threonine Residues in the C-Terminus of the Follicle-Stimulating Hormone Receptor Is Important for Desensitization but Not for Beta-Arrestin-Mediated ERK Activation[END_REF][START_REF] Tranchant | Preferential β-Arrestin Signalling at Low Receptor Density Revealed by Functional Characterization of the Human FSH Receptor A189 V Mutation[END_REF][START_REF] Tréfier | G Protein-Dependent Signaling Triggers a β-Arrestin-Scaffolded P70S6K/ RpS6 Module That Controls 5'TOP MRNA Translation[END_REF]. In addition to Gα s , FSHR displays coupling promiscuity [START_REF] Okashah | Variable G Protein Determinants of GPCR Coupling Selectivity[END_REF] with Gα q/11 [START_REF] Quintana | Follitropin (FSH) and a Phorbol Ester Stimulate the Phosphorylation of the FSH Receptor in Intact Cells[END_REF] and Gα i [START_REF] Crépieux | The ERK-Dependent Signalling Is Stage-Specifically Modulated by FSH, during Primary Sertoli Cell Maturation[END_REF]. The FSHR/FSH system is implicated in several reproduction-related pathologies both in women and men. These include premature ovarian insufficiency (POI), polycystic ovarian syndrome (PCOS), and male idiopathic infertility. Another severe condition is ovarian hyperstimulation syndrome (OHSS). Indeed, FSHR is the target of assisted reproductive techniques (ART), in which ovulation is induced through the injection of FSH preparations in a procedure called controlled ovarian stimulation. This procedure may lead to OHSS. Despite the fact that all aforementioned pathologies are well known and clinically described, the active role of FSHR, as well as the molecular mechanisms causing these disorders, are not entirely known. FSH induces multiple signaling pathways upon FSHR activation, making it difficult to understand which specific signaling module may be implicated in a specific pathology. In recent years, a growing body of evidence is highlighting the importance of dissecting GPCRs signaling pathways for two main reasons. [START_REF] Gromoll | Functional and Clinical Consequences of Mutations in the FSH Receptor[END_REF] In many cases, a particular GPCR-associated pathology is driven by a specific signaling pathway, and (2) current drugs are balanced agonists that result in both beneficial and unwanted side effects. We believe that having molecules capable of biasing signaling at the FSHR constitutes an important pharmacological tool to understand the specific FSHR signaling pathway contributing to pathologies (such as POI, PCOS) as well as to develop new generation drugs devoid of the adverse side effects of current clinical treatments (OHSS) [START_REF] Nataraja | Allosteric Regulation of the Follicle-Stimulating Hormone Receptor[END_REF].

Accumulating structural, biophysical, and pharmacological evidence has transformed our understanding of GPCR activation and pharmacological modulation. The fact that multiple active and inactive receptor conformations co-exist is supported by overwhelming structural and biophysical evidence [START_REF] Kobilka | Structural Insights into Adrenergic Receptor Function and Pharmacology[END_REF][START_REF] Nygaard | The Dynamic Process of β(2)-Adrenergic Receptor Activation[END_REF][START_REF] Wacker | Structural Features for Functional Selectivity at Serotonin Receptors[END_REF]. Consequently, efficacy is now considered multi-dimensional and explicitly incorporates the notion that receptors can engage distinct subsets of their full signaling repertoire [START_REF] Galandrin | The Evasive Nature of Drug Efficacy: Implications for Drug Discovery[END_REF]. This implies that different subsets of conformations can be stabilized by different ligands with distinct transduction mechanisms. This is the concept of pharmacological bias [START_REF] Kenakin | Ligand-Selective Receptor Conformations Revisited: The Promise and the Problem[END_REF][START_REF] Violin | Beta-Arrestin-Biased Ligands at Seven-Transmembrane Receptors[END_REF][START_REF] Reiter | Molecular Mechanism of β-Arrestin-Biased Agonism at Seven-Transmembrane Receptors[END_REF]. This class of ligands can modulate signaling in the presence (positive/negative allosteric modulators) or in the absence (allosteric agonists) of the endogenous ligand. Importantly, synthetic allosteric modulators of GPCRs can also lead to pharmacological bias [START_REF] Changeux | Allosteric Modulation as a Unifying Mechanism for Receptor Function and Regulation[END_REF][START_REF] Wootten | Emerging Paradigms in GPCR Allostery: Implications for Drug Discovery[END_REF]. It is also established that ligand bias must be set apart from system bias or observational bias [START_REF] Smith | Biased Signalling: From Simple Switches to Allosteric Microprocessors[END_REF].

To date, different classes of low molecular weight (LMW) molecules capable of modulating FSHR have been identified [START_REF] Nataraja | Discovery and Development of Small Molecule Allosteric Modulators of Glycoprotein Hormone Receptors[END_REF][START_REF] Ulloa-Aguirre | FSH Receptor Signaling: Complexity of Interactions and Signal Diversity[END_REF], including negative allosteric modulators (NAMs) [START_REF] Dias | A Negative Allosteric Modulator Demonstrates Biased Antagonism of the Follicle Stimulating Hormone Receptor[END_REF][START_REF] Dias | Inhibition of Follicle-Stimulating Hormone-Induced Preovulatory Follicles in Rats Treated with a Nonsteroidal Negative Allosteric Modulator of Follicle-Stimulating Hormone Receptor1[END_REF] and agonist/positive allosteric modulators (PAMs) [START_REF] Van Koppen | Signaling of an Allosteric, Nanomolar Potent, Low Molecular Weight Agonist for the Follicle-Stimulating Hormone Receptor[END_REF][START_REF] Sriraman | Investigation of a Thiazolidinone Derivative as an Allosteric Modulator of Follicle Stimulating Hormone Receptor: Evidence for Its Ability to Support Follicular Development and Ovulation[END_REF][START_REF] Yu | Discovery of Substituted Benzamides as Follicle Stimulating Hormone Receptor Allosteric Modulators[END_REF][START_REF] Nataraja | Discovery and Preclinical Development of Orally Active Small Molecules That Exhibit Highly Selective Follicle Stimulating Hormone Receptor Agonism[END_REF]. Small molecules have also been reported to behave as competitive antagonists [START_REF] Danesi | Clinical and Experimental Evidence of Inhibition of Testosterone Production by Suramin[END_REF][START_REF] Yanofsky | Allosteric Activation of the Follicle-Stimulating Hormone (FSH) Receptor by Selective, Nonpeptide Agonists[END_REF][START_REF] Anderson | Small Molecule Follicle-Stimulating Hormone Receptor Agonists and Antagonists[END_REF]. Of particular interest to this study, thiazolidinone derivatives were the first chemical molecules identified by combinatorial screening to specifically activate FSHR with no cross-reactivity with the other glycoprotein hormone receptors [START_REF] Wrobel | 5-Alkylated Thiazolidinones as Follicle-Stimulating Hormone (FSH) Receptor Agonists[END_REF]. Thiazolinidone derivatives did not compete with radiolabeled FSH for binding to FSHR, suggesting an allosteric mode of action. Moreover, slight changes in their side chains switched FSHR-coupling from Gα s to Gα i [START_REF] Arey | Differing Pharmacological Activities of Thiazolidinone Analogs at the FSH Receptor[END_REF]. However, bias factors were not determined. Compound 5 (renamed T1 thereafter) was able to induce preovulatory follicular development in immature rats but was not further pursued due to unfavorable pharmacokinetic parameters [START_REF] Sriraman | Investigation of a Thiazolidinone Derivative as an Allosteric Modulator of Follicle Stimulating Hormone Receptor: Evidence for Its Ability to Support Follicular Development and Ovulation[END_REF]. A separate class of ligands was identified through structure-activity relationship studies on substituted benzamides. Candidates capable of enhancing FSHR-mediated responses in the presence of EC 20 of FSH, with no cross-reactivity with luteinizing hormone choriogonadotropin receptor (LHCGR) or thyroid stimulating hormone receptor (TSHR), were identified [START_REF] Yu | Discovery of Substituted Benzamides as Follicle Stimulating Hormone Receptor Allosteric Modulators[END_REF]. This class of ligands has not been studied yet for their ability to activate FSHR in the absence of FSH. Interestingly, we recently reported that two chemically distinct LMW agonists belonging to these two chemical classes (i.e., benzamide, termed B3, and thiazolidinone, termed T1) were capable of pharmacologically programming endosomal signaling and post-endocytic trafficking of the FSHR [START_REF] Sposini | Pharmacological Programming of Endosomal Signaling Activated by Small Molecule Ligands of the Follicle Stimulating Hormone Receptor[END_REF]. In the present study, we systematically explored the ability of four LMW allosteric ligands belonging to these two chemical classes (i.e., B1, B2, B3, and T1, Figure 1) to induce biased signaling at the FSHR in the absence of FSH.

Results

Diversity in the Recruitment of G Proteins, β-Arrestin 2, and in cAMP Production

We first assessed the specificity of the effects mediated by B1, B2, B3, and T1. Mock transfected HEK293 cells treated with either LMW ligand or FSH showed no cAMP response, while forskolin led to a significant increase in cAMP (Supplementary Figure S1A). In HEK293 cells expressing FSHR, cAMP response induced by either FSH or LMW ligands could be antagonized using a previously validated antagonist (T2, Figure 1; Supplementary Figure S1B-F) [START_REF] Arey | Differing Pharmacological Activities of Thiazolidinone Analogs at the FSH Receptor[END_REF]. We next analyzed the recruitment of Gαs, Gαq, Gαi, and β-arrestin 2 to FSHR using Bioluminescence Resonance Energy Transfer (BRET). This technique allows monitoring protein-protein interaction by virtue of changes occurring in the energy transfer between a donor (Renilla luciferase) and an acceptor (fluorescent protein) when both are close enough (~10-100 angstroms) to generate a signal. Additionally, we assessed cAMP production using an intramolecular BRET sensor (both the donor and the acceptor are fused to a cAMP binding domain) that undergoes a conformational change, hence BRET change, upon cAMP binding. HEK293 cells expressing FSHR were challenged with increasing concentrations of FSH or LMW ligands over a 30 min period. Kinetics were generated for each concentration (Supplementary Figure S2), and the area under the curves was plotted as a concentration/activity sigmoid (Figure 2). We normalized each concentration/activity curve with respect to the FSH maximal response, as described in the materials and methods. As shown in Figure 2, LMW ligands displayed distinctive pharmacological profiles. Overall, small molecules were significantly less potent compared to FSH in all assays, as shown by the EC 50 comparison (Supplementary Table S1). Interestingly, B3 and T1 displayed significantly higher E max than FSH for the Gα s and Gα i subunit recruitment (Figure 2A,C), whereas they were both full agonists for Gαq recruitment (Figure 2B). In contrast, B1 and B2 behaved as full agonists for Gαs but as partial agonists for Gα q recruitment (Figure 2A,B). For Gα i recruitment, B1 was neutral with a slightly inverse agonism at the very high concentration, whereas B2 was a partial agonist (Figure 2C). Measurement of β-arrestin 2 recruitment revealed striking differences: T1 displayed full agonism, B1 and B2 behaved as weak partial agonists, whereas B3 recruited β-arrestin 2 with an almost two-fold higher E max (231.8 ± 4.9) compared to FSH, qualifying as a super-agonist for β-arrestin 2 recruitment (Figure 2D). For cAMP production, no major differences in E max were found compared to FSH, except for B1 displaying partial agonism (Figure 2E).

Differences in Kinetics That Affect Potencies and Efficacies

Comparing the kinetics across the different LMW ligands and concentrations, we observed markedly different kinetics compared to FSH (Supplementary Figure S2). Figure 3 shows time courses of G proteins and β-arrestin 2 recruitment to FSHR with each ligand used at its EC 80 as sub-maximal concentration and normalized to FSH response at 20 min (Figure 3). We extrapolated T 1/2 by using the general rise and fall exponential curve fitting equation as described in the materials and methods (Supplementary Tables S2 andS3). B3 led to significantly faster kinetics than FSH in recruiting Gα s , Gα q , Gα i , and β-arrestin 2. T1-induced kinetics were comparable to FSH for Gα s and Gα i recruitments but were slower for Gα q and β-arrestin recruitments as well as cAMP production. B1 displayed a weaker and more transient profile with all the read-outs considered. B2 displayed significantly slower but sustained activation kinetics with all read-outs compared with FSH, except for Gα s . It has been reported by others that kinetics can impact apparent biases in signaling [START_REF] Klein Herenbrink | The Role of Kinetic Context in Apparent Biased Agonism at GPCRs[END_REF]. In this context, we compared the AUC (area under the curve) for each dose with endpoints at 2-, 7-and 20-min. Signals from each LMW ligand were normalized over the FSH signal taken at the same time point. Consistent with Klein Herenbrinck et al., we found that the different kinetics induced by T1 could indeed affect potency and efficacy as a function of the endpoint considered. Figure 4 demonstrates that T1-induced E max and potencies were significantly different from FSH for both Gα s recruitment and cAMP production as a function of the stimulation time.

Diversity of FSHR-Mediated ERK Phosphorylation and cAMP-Responsive Element-Dependent Transcription

We next assessed the effect of each LMW ligand on two integrated cellular responses: (i) extracellular signal-regulated kinase (ERK1/2) phosphorylation, assessed by Western blot, and (ii) cAMP-responsive element (CRE) transcription using a reporter gene (Figure 5). ERK1/2 phosphorylation levels were evaluated after stimulating FSHR-transfected HEK293 cells for 5 min with EC 80 of each compound (Supplementary Figures S3 andS4) and were normalized over total ERK. The 5 min time corresponded to the peak of phosphorylated ERK activation in our cellular settings, as previously described [START_REF] Piketty | Follicle-Stimulating Hormone (FSH) Activates Extracellular Signal-Regulated Kinase Phosphorylation Independently of Beta-Arrestin-and Dynamin-Mediated FSH Receptor Internalization[END_REF][START_REF] Sposini | Integration of GPCR Signaling and Sorting from Very Early Endosomes via Opposing APPL1 Mechanisms[END_REF]. As shown in Figure 5, FSH and LMW ligands induced ERK phosphorylation compared to basal. Interestingly, B3 induced a higher level of ERK phosphorylation than the other LMW ligands and FSH. The Gαs/cAMP/PKA pathway is canonically activated upon FSH binding to the receptor. In particular, the second messenger cAMP binds to PKA, which, in turn, activates the cAMP response element-binding protein (CREB) in the nucleus. CREB is a transcription factor that binds to the cAMP response element (CRE) of the promoters of its target genes and leads to their activation. In this context, we were interested in analyzing the impact of LMW ligands on the FSHR-mediated Gαs/cAMP/PKA-dependent gene expression. To do this, CRE-dependent transcription was assessed after 6 h of stimulation using luciferase as a reporter gene. The data revealed striking differences compared to Gα s recruitment and cAMP production. B3 and T1 were partial agonists compared to FSH, whereas B1 and B2 displayed very poor potencies, inducing weak luciferase expression only at high concentrations (Figure 5C; Supplementary Table S4).

Profiling in HEK293/∆Gα s and HEK293/∆ARRB Cells

To further dissect the transduction mechanisms triggered by the different ligands, but also to assess the potential impact of system biases, we carried out the same profiling in HEK293 cells depleted in two key molecular effectors of FSHR, namely either Gα s or β-arrestins. The levels of FSHR expression at the plasma membrane were quantified across the three cellular models by flow cytometry and were found to be similar between wild-type HEK293 and HEK293/∆Gα s , but slightly lower for HEK293/∆ARRB (Supplementary Figure S5).

Gα q , Gα i , and β-arrestin 2 recruitment, and CRE-dependent transcription were measured in HEK293/∆Gα s cells stimulated with increasing doses of FSH or LMW ligands (Figure 6, Supplementary Figure S6). FSH-induced signals were arbitrarily set at 100%. The lack of Gα s affected the relative efficacies of some ligands compared to FSH. In particular, B3 and T1 were less efficacious in recruiting Gα i compared to the unmodified cells, while being more efficacious for Gα q recruitment (Supplementary Tables S1 andS5). Interestingly, both B3 and T1 were also more efficacious than FSH in recruiting β-arrestin 2. B2 displayed a higher efficacy for Gα q and β-arrestin 2 recruitments but was ineffective for Gα i recruitment. Interestingly, compound B1 did not display any notable efficacy regardless of the read-out in these cells. As expected, Gα s depletion strongly impaired CRE-dependent transcription, yet a clear residual activation remained and was analyzed. Compound B3 displayed weak efficacy and potency, whereas T1 induced CRE-dependent transcription with higher efficacy compared to wild-type cells. Interestingly, B1 and B2 did not trigger any measurable responses. Kinetics revealed that Gα i recruitment was severely affected for all ligands and FSH compared to the wild-type cells (Supplementary Figure S6).

We next investigated the effects of our ligands in a system where both β-arrestin 1 and 2 were genetically depleted (i.e., HEK293/∆ARRB). Again, FSHR-transfected cells were challenged with increasing doses of FSH or LMW ligands. G proteins, cAMP production, and CRE-luciferase expression were analyzed (Figure 7; Supplementary Table S5). B1 reached higher efficacy in Gα s recruitment, while B3 and T1 showed less efficacy in recruiting Gα i in HEK293/∆ARRB compared to the wild-type counterpart. However, the general overall behavior of these three LMW ligands on upstream signaling was maintained in HEK293/∆ARRB compared to wild-type HEK293. Interestingly, the lack of β-arrestin significantly affected B3-induced reporter gene expression as shown by the decreased efficacy compared to wild-type HEK293 (71.0 ± 2.7 vs. 53.4 ± 7.1, respectively), consistent with this compound being the most efficacious at inducing β-arrestin 2 recruitment in wild-type cells. The kinetics indicated that β-arrestins depletion of β-arrestins did not impair Gα i recruitment the way Gα s depletion did (Supplementary Figure S7).

Analysis of Bias Factors and Multi-Dimensional Comparison of Pharmacological Efficacies

To visualize the pharmacological activity induced by the LMW ligands, we first graphically represented the values for the transduction coefficients (Supplementary Tables S6-S8) and the maximal efficacies (Supplementary Tables S1A andS5A,C) in radial graphs (Figure 8). Such representation visually correlates an intrinsic pharmacological parameter (i.e., the transduction coefficient) to the maximal efficacy of each ligand for the read-outs considered. As expected, FSH-induced transduction coefficients and maximal efficacies were balanced in the three cellular models considered. On the contrary, LMW ligands showed imbalance for different read-outs, highlighting a high pharmacological diversity. It is interesting to note that no B1-induced transduction coefficients in the HEK293/∆Gα s cells were statistically identifiable except for the read-out Gαq; therefore, the radial plot was not determined (nd).

We then calculated bias factors for each compound compared to the endogenous ligand FSH and for each pair of read-outs across the three cellular models, using the operational model of agonism, as described in the materials and methods (Supplementary Tables S9-S11). In Figure 9, bias values were graphically represented as arrows pointing to the direction of the bias (i.e., one read-out in each pair considered). In several cases, as described in the materials and methods, we were not able to accurately determine the bias value, reflecting a concentration-activity curve that was very weak or flat, for at least one of the two read-outs. In such cases, the fitting of the sigmoid is not robust and can lead to unreliable bias values. Since there was a clear bias towards one read-out, where we could not perform statistical analysis, we decided to represent such cases with empty arrows (observed bias) as opposed to full arrows (statistically significant bias). As shown in Figure 9, each compound displayed a unique profile, indicating pharmacological diversity. We also noticed that bias profiles varied as a function of the cell model considered, confirming the importance of the cellular context when assessing biases. Noticeably, all ligands were biased towards Gα s compared to cAMP in both wild-type and HEK293/∆ARRB cells. There was also a general trend in favor of a Gα s -oriented bias for all ligands in HEK293/∆ARRB cells as well. It should be noted that B1 appeared as the most balanced ligand, mostly due to its low efficacy and potency in many read-outs that generated no bias values. This phenomenon is particularly visible in wild-type HEK293 and HEK293/∆Gα s cells in which B1 did not display any activity except a significant bias in favor of Gα s over cAMP in wild-type HEK293. Consistent with its strong dependence on Gα s , B1 was also largely biased in favor of Gα s recruitment in HEK293/∆ARRB. Interestingly, B3 was significantly biased for β-arrestin 2 compared to cAMP and CRE-luciferase in wild-type, while T1 was significantly biased for arrestin 2, only compared to CRE-luciferase. Moreover, in wild-type cells, ligands B3 and T1 were significantly biased towards Gα i compared to cAMP, while B3, B2, and T1 were significantly biased towards cAMP compared to Gα q . These two biases remained consistent in HEK293/∆ARRB cells.

To better discriminate pharmacological efficacies of FSH and LMW ligands across the multiple dimensions in the three cellular models explored, we performed principal component analyses (PCA) using the calculated bias values. Since FSH was the reference ligand, all the values entered for FSH in the PCAs were set to 0, meaning no bias. In all three PCAs shown in Figure 10, the two principal components, PC1+PC2, accounted for 91.44 to 99.14% of the variability across the ligands, hence highlighting differences in ligand behavior. Furthermore, we performed hierarchical clustering based on the PCA results to group ligands in clusters (Figure 10). Interestingly, ligands did not group in the same clusters (dotted lines) in the different cellular models, suggesting that the gene knockouts and the associated compensatory rewiring [START_REF] Luttrell | Manifold Roles of β-Arrestins in GPCR Signaling Elucidated with SiRNA and CRISPR/Cas9[END_REF] led to system biases. In wild-type HEK293, B3, FSH, and T1 were clustered together; B1, T1, and FSH were in the same cluster in HEK293/∆Gα s , whereas two separate clusters, composed by B3 and B2 (left cluster) and by FSH and T1 (right cluster), were found in HEK293/∆ARRB. The apparent superposition of B1 and FSH in HEK293/∆Gα s was clearly misleading and is not at all indicative of pharmacological equivalence of these two ligands in this cellular system. Instead, it reflects the fact that FSH, being the reference, had no bias in any of the cellular systems (hence all values = 0), yet it clearly triggered robust responses for all read-outs measured in HEK293/∆Gαs. On the contrary, B1 displayed no activity at all in HEK293/∆Gαs, hence, no bias (all values = 0). This is an illustration of the type of caveats that can arise when running such systematic bias analyses. , and (E) cAMP production were measured. HEK293 cells transfected with the FSHR and the appropriate BRET sensors (Rluc stands for renilla luciferase, YFP for yellow fluorescent protein, and EPAC for exchange protein directly activated by cAMP) were challenged with increasing doses of FSH (10 -12.5 to 10 -6.5 M) or LMW ligands (10 -10 to 10 -4 M). Concentration/activity curves were measured for the different read-outs. Signals were monitored for 30 min and normalized, considering the FSH maximal response as 100%. The area under the concentration/activity curves generated by each compound was plotted and fitted using nonlinear regression. Data are represented as mean ± SEM from at least six independent experiments.

Differences in Kinetics that Affect Potencies and Efficacies

Comparing the kinetics across the different LMW ligands and concentrations, we observed markedly different kinetics compared to FSH (Supplementary Figure S2). Figure 3 shows time courses of G proteins and β-arrestin 2 recruitment to FSHR with each ligand used at its EC80 as sub-maximal concentration and normalized to FSH response at 20 min

Diversity of FSHR-Mediated ERK Phosphorylation and cAMP-Responsive Element-Dependent Transcription

We next assessed the effect of each LMW ligand on two integrated cellular responses: (i) extracellular signal-regulated kinase (ERK1/2) phosphorylation, assessed by Western blot, and (ii) cAMP-responsive element (CRE) transcription using a reporter gene (Figure 5). ERK1/2 phosphorylation levels were evaluated after stimulating FSHR-transfected HEK293 cells for 5 min with EC80 of each compound (Supplementary Figures S3 andS4) and were normalized over total ERK. The 5 min time corresponded to the peak of phosphorylated ERK activation in our cellular settings, as previously described [START_REF] Piketty | Follicle-Stimulating Hormone (FSH) Activates Extracellular Signal-Regulated Kinase Phosphorylation Independently of Beta-Arrestin-and Dynamin-Mediated FSH Receptor Internalization[END_REF][START_REF] Sposini | Integration of GPCR Signaling and Sorting from Very Early Endosomes via Opposing APPL1 Mechanisms[END_REF]. As shown in Figure 5, FSH and LMW ligands induced ERK phosphorylation compared to basal. Interestingly, B3 induced a higher level of ERK phosphorylation than the other LMW ligands and FSH. The Gαs/cAMP/PKA pathway is canonically activated upon FSH x FOR PEER REVIEW 9 of 24 inducing weak luciferase expression only at high concentrations (Figure 5C; Supplementary Table S4). transcription, yet a clear residual activation remained and was analyzed. Compound B3 displayed weak efficacy and potency, whereas T1 induced CRE-dependent transcription with higher efficacy compared to wild-type cells. Interestingly, B1 and B2 did not trigger any measurable responses. Kinetics revealed that Gαi recruitment was severely affected for all ligands and FSH compared to the wild-type cells (Supplementary Figure S6). superposition of B1 and FSH in HEK293/ΔGαs was clearly misleading and is not at all indicative of pharmacological equivalence of these two ligands in this cellular system. Instead, it reflects the fact that FSH, being the reference, had no bias in any of the cellular systems (hence all values = 0), yet it clearly triggered robust responses for all read-outs measured in HEK293/ΔGαs. On the contrary, B1 displayed no activity at all in HEK293/ΔGαs, hence, no bias (all values = 0). This is an illustration of the type of caveats that can arise when running such systematic bias analyses.

Discussion

FSHR has long been a target for pharmaceutical companies since FSH is administered to patients in ovarian stimulation protocols [START_REF] Riccetti | Genetics of Gonadotropins and Their Receptors as Markers of Ovarian Reserve and Response in Controlled Ovarian Stimulation[END_REF]. Over the years, several attempts were made to identify suitable orally-active LMW substitutes of FSH [START_REF] Nataraja | Discovery and Development of Small Molecule Allosteric Modulators of Glycoprotein Hormone Receptors[END_REF][START_REF] Nataraja | Discovery and Preclinical Development of Orally Active Small Molecules That Exhibit Highly Selective Follicle Stimulating Hormone Receptor Agonism[END_REF][START_REF] Yanofsky | Allosteric Activation of the Follicle-Stimulating Hormone (FSH) Receptor by Selective, Nonpeptide Agonists[END_REF][START_REF] Gerrits | Oral Follicle-Stimulating Hormone Agonist Tested in Healthy Young Women of Reproductive Age Failed to Demonstrate Effect on Follicular Development but Affected Thyroid Function[END_REF]. However,

Discussion

FSHR has long been a target for pharmaceutical companies since FSH is administered to patients in ovarian stimulation protocols [START_REF] Riccetti | Genetics of Gonadotropins and Their Receptors as Markers of Ovarian Reserve and Response in Controlled Ovarian Stimulation[END_REF]. Over the years, several attempts were made to identify suitable orally-active LMW substitutes of FSH [START_REF] Nataraja | Discovery and Development of Small Molecule Allosteric Modulators of Glycoprotein Hormone Receptors[END_REF][START_REF] Nataraja | Discovery and Preclinical Development of Orally Active Small Molecules That Exhibit Highly Selective Follicle Stimulating Hormone Receptor Agonism[END_REF][START_REF] Yanofsky | Allosteric Activation of the Follicle-Stimulating Hormone (FSH) Receptor by Selective, Nonpeptide Agonists[END_REF][START_REF] Gerrits | Oral Follicle-Stimulating Hormone Agonist Tested in Healthy Young Women of Reproductive Age Failed to Demonstrate Effect on Follicular Development but Affected Thyroid Function[END_REF]. However, no LMW ligand-exhibiting undisputable and well-characterized biased agonisms at the FSHR have been reported so far. In the current study, we took advantage of previously described thiazolidinone and benzamide derivatives and investigated their ability to trigger biased signaling compared to FSH [START_REF] Yu | Discovery of Substituted Benzamides as Follicle Stimulating Hormone Receptor Allosteric Modulators[END_REF][START_REF] Wrobel | 5-Alkylated Thiazolidinones as Follicle-Stimulating Hormone (FSH) Receptor Agonists[END_REF]. Of note, our panel included three benzamide derivatives that presented a common core structure. We capitalized on highly-sensitive BRET sensors to carry out multiplexed measurements encompassing the main known coupling mechanisms of the FSHR: Gα s , Gα q , Gα i , and β-arrestin 2 recruitments [START_REF] Gloaguen | Mapping the Follicle-Stimulating Hormone-Induced Signaling Networks[END_REF][START_REF] Wan | Mini G Protein Probes for Active G Protein-Coupled Receptors (GPCRs) in Live Cells[END_REF]. We also assessed cAMP production, ERK phosphorylation, and CRE-dependent transcription-three downstream cellular responses known to be induced by FSH. Concentration-activity experiments were conducted for all read-outs but ERK phosphorylation to achieve quantitative comparisons (efficacy and potency) between FSH and the different ligands. Then, we applied operational model-derived analysis to the whole dataset to calculate the transduction coefficients and bias factors of all the possible pairs of read-outs.

Our data confirmed that FSH and the LMW ligands were able to recruit Gα q and Gα i , in addition to the canonical Gα s protein; thereby, supporting the existence of G protein promiscuity at the FSHR [START_REF] Quintana | Follitropin (FSH) and a Phorbol Ester Stimulate the Phosphorylation of the FSH Receptor in Intact Cells[END_REF]. Importantly, our results reveal that the different LMW ligands lead to distinct biases at the FSHR. Indeed, B3 behaved as a super-agonist for Gα s and β-arrestin 2 recruitment compared to FSH. Moreover, it was statistically biased towards β-arrestin 2 when compared to cAMP. Consistent with this finding, we observed that B3 induced more ERK phosphorylation than FSH in wild-type HEK293 cells and that its CRE-dependent transcription was more affected than the one elicited by FSH or other LMW ligands in HEK293/∆ARRB cells. T1 also showed β-arrestin 2 bias compared to cAMP and CRE-dependent transcription. Similar to B3, T1 also demonstrated twofold higher efficacy for Gα s protein recruitments while being slightly more efficacious at inducing cAMP production compared to FSH. However, its β-arrestin recruitment efficacy was significantly reduced compared to B3. In contrast, B1 displayed full agonism for Gα s and partial agonism for cAMP production, respectively, in wild-type HEK293 cells; whereas, it showed little to no pharmacological activity in all other assays, including β-arrestin recruitment. In line with this Gα s -biased profile, B1 failed in inducing any significant response in HEK293/∆Gα s cells. Compound B2, on the other hand, displayed full agonism for Gα s recruitment and cAMP production but was only partial agonist for all other read-outs. These data demonstrated an unprecedented pharmacological diversity at the FSHR. These ligands represent novel molecular tools that can serve to further dissect FSHR signaling.

We noticed that LMW ligands displayed less bias towards cAMP production compared with Gα and β-arrestin recruitment. The exact mechanisms that determine the orientation of the bias remain largely unknown, even though more studies are attempting to link signaling behaviors to receptor structural conformations. In the case of cAMP production vs. Gα and β-arrestin recruitment, the stoichiometry of the elements at play is vastly different: for Gα or β-arrestin recruitment to the receptor, the stoichiometry is 1:1, meaning that one activated receptor recruits one Gα protein/β-arrestin. In contrast, cAMP production is strongly amplified: activation of one receptor leads to plenty of cAMP molecules produced, making cAMP response less sensitive to pharmacological modulation. We previously reported that only a few percent of occupied FSHRs still generate maximum cAMP responses, whereas, Gα or β-arrestin recruitment are perfectly correlated with receptor occupancy [START_REF] Ayoub | Assessing Gonadotropin Receptor Function by Resonance Energy Transfer-Based Assays[END_REF]. This discrepancy may explain why LMW compounds showed less bias towards cAMP production.

Importantly, distinct pharmacological profiles were obtained with three molecules that are relatively similar. Indeed, B1, B2, and B3 share about 90% of their chemical structure [START_REF] Yu | Discovery of Substituted Benzamides as Follicle Stimulating Hormone Receptor Allosteric Modulators[END_REF], yet their pharmacological profiles were widely different. The same is true with T1 and T2, which also share 90% of their chemical structure. However, T1 was a biased agonist, while T2 was an antagonist [START_REF] Arey | Differing Pharmacological Activities of Thiazolidinone Analogs at the FSH Receptor[END_REF]. This observation is in line with the constantly growing number of examples of slight chemical modifications leading to profound pharmacological biases at other GPCRs.

Another interesting finding of our study is that, despite some of the LMW ligands (i.e., B3 and T1) being full-or even super-agonists for G protein recruitments, cAMP production and β-arrestin 2 recruitment compared to FSH was less efficient for CRE-dependent transcription. From this observation, we hypothesized that the different kinetic profiles elicited at the FSHR by the different ligands could account for the differences measured at the transcriptional level as well as for ligand-specific receptor assembly at the membrane. Signal compartmentalization could also play a role as it was recently shown that FSHR and LHCGR produce a large proportion of their cAMPs from endosomes [START_REF] Sposini | Pharmacological Programming of Endosomal Signaling Activated by Small Molecule Ligands of the Follicle Stimulating Hormone Receptor[END_REF][START_REF] Sposini | Integration of GPCR Signaling and Sorting from Very Early Endosomes via Opposing APPL1 Mechanisms[END_REF]. It is therefore conceivable that LMW ligands lead to differential activations of the plasma-membrane vs. endosomal cAMP production. Indeed, increased T1-induced cAMP signaling is endosomal and, at the same time, is insensitive to the negative regulation of APPL1. Interestingly, the importance of kinetics in pharmacological efficacy at the dopamine D 1 receptor has been recently explored [START_REF] Klein Herenbrink | The Role of Kinetic Context in Apparent Biased Agonism at GPCRs[END_REF]. In particular, the D 1 receptor did not display the same kinetics of activation with the different ligands analyzed. Similarly, we found that concentration/activity curves made from different endpoints could lead to different relative efficacies and potencies compared to FSH. Looking at the full kinetics, we found that LMW ligands displayed different profiles that were slow, fast, or transient compared to FSH. Our study further highlights the importance of considering kinetics when pharmacologically profiling ligands.

Given the biased properties of the different LMW ligands in our panel, one can wonder what effects they would have on reproduction. For example, the selective activation of a specific FSHR-mediated signaling pathway by a biased compound may enhance the ovulation output in ART while reducing the risk for OHSS. In a similar way, such biased properties may contribute to dissect the molecular mechanism leading to reproduction-related disorders. Interestingly, T1 was already tested in immature rats.

Despite good potency and efficacy in vitro, serum T1 was undetectable in vivo. Effects on ovulation were achieved only when T1 was continuously perfused with a pump. Moreover, some genotoxic effects were reported [START_REF] Sriraman | Investigation of a Thiazolidinone Derivative as an Allosteric Modulator of Follicle Stimulating Hormone Receptor: Evidence for Its Ability to Support Follicular Development and Ovulation[END_REF]. On the contrary, B1, B2, and B3 were not yet evaluated in vivo. Other benzamide derivatives from the same series were tested in immature rats. They displayed acceptable pharmacokinetic parameters but were not tested for their effects on ovulation [START_REF] Yu | Discovery of Substituted Benzamides as Follicle Stimulating Hormone Receptor Allosteric Modulators[END_REF]. In this context, evaluating B1, B2, and B3 in vivo could help decipher the respective contributions of G proteins and β-arrestin on reproduction, thereby assessing the potential of FSHR-targeted biased drugs for clinical applications in the field of reproductive medicine.

System bias has been identified as a potential confusion factor for pharmacological profiling of ligand bias [START_REF] Smith | Biased Signalling: From Simple Switches to Allosteric Microprocessors[END_REF]. System bias is often evoked when comparing different cell types and/or physiological contexts. Here we reasoned that the advent of CRISPR Cas9 knock-out models of HEK293 could be an interesting opportunity to look at system bias from a different perspective. Indeed, it was recently reported that such gene depletion forced an important rewiring of signaling pathways in order for the knock-out cells to adapt and compensate for the missing factor [START_REF] Luttrell | Manifold Roles of β-Arrestins in GPCR Signaling Elucidated with SiRNA and CRISPR/Cas9[END_REF]. We compared, side-by-side, three models of HEK293 cells wild-type or depleted for Gα s and both β-arrestin isoforms. We clearly found that our ligands did not display the same pharmacological profiles compared to FSH in the different models, further confirming the importance of system bias. We demonstrated that even a discrete alteration affecting one transduction partner could lead to pathway rewiring, completely changing the pharmacological profiles, even though the receptor was expressed at similar or slightly lower levels and the assays were identical.

In addition, we observed that Gα i recruitment was profoundly reduced for all LMW ligands and for FSH in HEK293/∆Gα s compared to wild-type and ∆ARRB HEK293 cells. Even though the residual effects were sufficient to determine potencies and efficacies for Gα i recruitment, the drastic reduction of the overall coupling to Gα i in the absence of Gα s suggests that a possible Gα s /Gα i switching occurs at the FSHR. A switching from Gα s to Gα i coupling, subsequent to phosphorylation of β2-adrenergic by PKA, was previously reported [START_REF] Daaka | Switching of the Coupling of the Beta2-Adrenergic Receptor to Different G Proteins by Protein Kinase A[END_REF][START_REF] Zamah | Protein Kinase A-Mediated Phosphorylation of the Beta 2-Adrenergic Receptor Regulates Its Coupling to Gs and Gi. Demonstration in a Reconstituted System[END_REF]. This is yet another illustration of the importance and confounding nature of system bias.

In addition to the importance of kinetics and system biases, our study pointed to some limitations in the way pharmacological biases are currently statistically measured. Stateof-the-art methods rely on the Black and Leff operational model of agonism to quantify bias [START_REF] Rajagopal | Quantifying Ligand Bias at Seven-Transmembrane Receptors[END_REF][START_REF] Kenakin | A Simple Method for Quantifying Functional Selectivity and Agonist Bias[END_REF]. The main limitation of this approach is that it is based on fitting a sigmoid to concentration/activity curves. It works robustly when actual responses can be measured. However, as expected in the case of biased ligands, some ligands do not exhibit clear efficacy in some assays. In such a case, the operational model still fits a sigmoid, but the inflection points are not reliably identified, and this generates aberrant biased values. A modification of the operational model was recently proposed to deal with such cases of extreme biases [START_REF] Stahl | A Novel Method for Analyzing Extremely Biased Agonism at G Protein-Coupled Receptors[END_REF]. However, this approach only works for orthosteric ligands, as it relies on competition; therefore, it is not applicable to allosteric ligands. Further studies will be necessary to address this important limitation in the future and allow reliable pharmacological profiling in all situations.

Overall, we report for the first time a panel of LMW agonists presenting varied biased profiles at the FSHR. This finding will help decipher the relative contributions to the reproduction-related physiopathology of the different signaling pathways triggered by the FSHR. In that sense, our study potentially represents a significant step towards the development of improved therapeutics in reproductive medicine. Such molecules could be used to gain a better understanding of FSHR signaling, a pre-requisite to the discovery of novel FSHR-targeting drugs. In structure-activity relationship studies, these compounds may also be used to determine the chemical elements responsible for a particular drug effect.

Materials and Methods

Ligands and Materials

The FSHR low molecular weight ligands (B1, B2, B3, T1, T2) were provided by To-copheRx (Burlington, VT, USA). Classification and chemical structures are displayed in Figure 1. Ligands were diluted in sterile dimethyl sulfoxide (DMSO, Sigma-Aldrich, Saint-Louis, MO, USA) for conservation at -20 • C then diluted in PBS (Eurobio, Les Ulis, France) at appropriate concentrations for stimulations. Recombinant FSH was kindly donated by Merck (Darmstadt, Germany) and diluted in PBS. Ninety-six-well white microplates were acquired from Greiner Bio-one (Courtaboeuf, France). Coelenterazine H was purchased from Interchim (Montluçon, France).

Cell Culture and Transfection

HEK293/∆Gαs were generated by introducing null mutations into the GNAS gene, and the GNAL genes and HEK293/∆ARRB were mutated in the ARRB1 and the ARRB2 genes, using the CRISPR/Cas9 genome editing technology as described elsewhere [START_REF] O'hayre | Genetic Evidence That β-Arrestins Are Dispensable for the Initiation of B2-Adrenergic Receptor Signaling to ERK[END_REF][START_REF] Stallaert | Purinergic Receptor Transactivation by the B2-Adrenergic Receptor Increases Intracellular Ca2+ in Nonexcitable Cells[END_REF][START_REF] Grundmann | Lack of Beta-Arrestin Signaling in the Absence of Active G Proteins[END_REF]. Cells were cultured in a complete DMEM medium, supplemented with 10% (v/v) fetal bovine serum, 4.5 g/L glucose, 100 U/mL penicillin, 0.1 mg/mL streptomycin, and 1mM glutamine (all from Thermo Fisher Scientific, Waltham, MA, USA). Transient transfections were performed in suspension in a 96-well plate using Metafectene PRO (Biontex, München, Germany) following the manufacturer's protocol.

BRET Sensors

To evaluate Gα s , Gα q , Gα i recruitment, cells were transiently transfected with previously described FSHR C-terminally fused with donor Rluc8 (FSHR-Rluc8) and acceptor NES-Venus-mGs, or NES-Venus-mGsq or NES-Venus-mGsi (kindly provided by Pr. N.A. Lambert, Augusta University, Augusta, GA, USA) [START_REF] Wan | Mini G Protein Probes for Active G Protein-Coupled Receptors (GPCRs) in Live Cells[END_REF][START_REF] Ayoub | Assessing Gonadotropin Receptor Function by Resonance Energy Transfer-Based Assays[END_REF]. In order to measure real-time cAMP response in living cells, HEK293 cells were transiently transfected with two plasmids coding for FSHR [START_REF] Tranchant | Preferential β-Arrestin Signalling at Low Receptor Density Revealed by Functional Characterization of the Human FSH Receptor A189 V Mutation[END_REF] and the BRET-based cAMP sensor CAMYEL (kindly provided by L.I. Jiang, University of Texas, Austin, TX, USA) [START_REF] Jiang | Use of a CAMP BRET Sensor to Characterize a Novel Regulation of CAMP by the Sphingosine 1-Phosphate/G13 Pathway[END_REF]. For the assessment of β-arrestin 2 recruitment, HEK293 cells were transiently co-transfected with plasmids coding for FSHR-Rluc8 and for β-arrestin 2 N-terminally fused to the BRET acceptor yPET (kindly provided by Dr. M.G. Scott, Cochin Institute, Paris, France) [START_REF] Ayoub | Assessing Gonadotropin Receptor Function by Resonance Energy Transfer-Based Assays[END_REF].

Cell Stimulation and BRET Measurement

Forty-eight hours after transfection, BRET measurements were immediately performed upon addition of increasing concentrations of FSH (10 -12.5 to 10 -6.5 M) or of LMW ligands (10 -10 to 10 -4 M) and 5 µM of coelenterazine H diluted in PBS. Signals were recorded for 30 min in a Mithras LB 943 plate reader (Berthold Technologies GmbH & Co., Wildbad, Germany).

Western Blot Analysis

Cells were stimulated with 80% effective concentrations (EC 80 ) of FSH or LMW ligands for 5 min, and then the supernatants were removed. Cells were then lysed in 2X Laemmli buffer (250 mM tris-HCl (pH 6.8), 2% (w/v) SDS, 10% (v/v) glycerol, 0.01% (w/v) bromophenol blue, and 5% (v/v) β-mercaptoethanol) and equal amounts of protein were analysed by Western blot. The membranes were incubated overnight at 4 • C with anti-phospho-ERK 1/2 (P-p44/42 MAPK; Rabbit mAB; Ref:4370S; Lot 17 or Lot 23; 1:3000 or 1:1000; Cell Signalling Technologies, Danvers, MA, USA) and were reprobed with primary polyclonal antibody against total ERK2 (ERK2(C-14); sc-154; Lot #K2012; rabbit polyclonal IgG; 1:10,000; Santa Cruz Biotechnology, Dallas, TX, USA) or anti-p42/44 ERK (Cell Signalling Technology, catalog number: 9102, 1:10,000, Ref: 9107, Lot 28) followed by IRDye conjugates (IRDye800CW, 926-68071, Lot C70901-13/ IRDye680RD, 926-32211, Lot 70926-05; goat anti-rabbit; 1:300; Li-Cor Biosciences, Lincoln, NE, USA) or Horseradish Peroxidase HRP conjugated mouse anti-rabbit (Santa Cruz, sc-2357, Lot L1218, 1:10,000) as secondary antibodies. Protein bands were detected with a LI-Cor Odyssey CLX far-infrared scanner (Li-Cor Biosciences, Lincoln, NE, USA), or membranes were subsequently exposed to autoradiography films (GE Healthcare Amersham TM Hyperfilm TM ECL, 10534205) in a dark room. Bands were quantified by densitometry with Image Studio Ver. 2.1 (Li-Cor Biosciences, Lincoln, NE, USA).

Flow Cytometry

Cells (wild type HEK293, HEK293/∆Gα s, and HEK293∆ARRB) were transiently transfected with either plasmid encoding for an empty pcDNA3.1 vector (mock) or N-terminallytagged FLAG-FSHR. Cells were detached, washed, and re-suspended in working buffer (PBS w/o Ca 2+ and Mg 2+ , 1% BSA, 2 mM EDTA) 48 h after transfection. Cells were then incubated with PE-conjugated anti-FLAG antibody (anti-DYKDDDDK-PE, Order n.: 130-101-577, Lot: 5190225429, 1:100, Milteny Biotech, Bergisch Gladbach, Germany) for 1 h at 4 • C. Cells were washed twice and re-suspended in working buffer before analysis with MACSQuant Analyzer 10 Flow cytometer (Milteny Biotech, Bergisch Gladbach, Germany). Data were analyzed and plotted with FlowJo Ver. 9.6 (FlowJo, Ashland, OR, USA).

Cre-Dependent Reporter Assay

HEK293 cells were transiently transfected with FSHR and the pSOM-Luc plasmid expressing the firefly luciferase reporter gene under the control of the cAMP Responsive Element of the somatostatin promoter region [START_REF] Troispoux | Involvement of G Protein-Coupled Receptor Kinases and Arrestins in Desensitization to Follicle-Stimulating Hormone Action[END_REF]. After 48-h transfection, cells were stimulated for 6 h with increasing doses of FSH (10 -12.5 to 10 -6.5 M) or LMW ligands (10 -10 to 10 -4 M). Luciferase activity was measured by adding Bright-Glo Luciferase assay substrate (Promega, Madison, WI, USA), and the emitted light was measured in a Mithras LB 943 plate reader. Values were expressed in relative luciferase activity units (RLU).

Statistical Analysis of Concentration/Activity Curves

The results are shown as mean ± SEM from at least six independent experiments. Data were analyzed and plotted using GraphPad Prism Ver.9.0 (San Diego, CA, USA). For normalizations, the values of all replicates were divided by the mean of FSH-induced maximal responses and multiplied by 100 for any given read-out. Statistical significance was assessed by an unpaired t-test with Welch's correction, based on the number of datasets in each given experiment. p-values were considered significant when <0.05.

Statistical Analysis of Kinetic Curves

The results are shown as mean ± SEM from at least six independent experiments. T 1/2 were obtained after curve fitting by using the general rise and fall exponential equation in GraphPad Prism:

E = C/(K 1 -K 2 ) × (e -K 2 t -e -K 1 t )
where T 1/2 = log(2)/K 1 , as previously described by Hoare et al., 2018 [START_REF] Hoare | Kinetic Operational Models of Agonism for G-Protein-Coupled Receptors[END_REF].

Bias Calculation and Statistics

Transduction coefficients (R = τ/Ka) and bias factor (BF) were obtained after statistical fitting (Supplementary Figures S8-S10) of the operational model for each dose-response curve using normalized data from Figures 2, 6 and7. Starting from the procedure described in van der Westhuizen et al. [START_REF] Van Der Westhuizen | Quantification of Ligand Bias for Clinically Relevant B2-Adrenergic Receptor Ligands: Implications for Drug Taxonomy[END_REF], performed in GraphPad Prism, we adapted the methodology using the D2D-software on Matlab Vers. 2016b (MathWorks, Natick, MA, USA) [START_REF] Raue | Data2Dynamics: A Modeling Environment Tailored to Parameter Estimation in Dynamical Systems[END_REF] to improve the statistical significance and avoid possible misleading procedures in data fitting due to parameter non-identifiability. For a given response and a given ligand, the corresponding concentration/activity curve was fitted using the operational model given by Equation ( 1):

E = Basal + (E max -Basal) × (τ[L])ˆn/((τ[L])ˆn + (Ka + [L])ˆn), (1)
where [L] denotes the concentration of ligand and E is the quantification of its effect (such as BRET ratio, for instance). The basal parameter is the baseline of the response; E max is the maximal possible response of the system; τ is the efficacy; Ka is the functional equilibrium dissociation constant of the agonist; n is the Hill slope of the transducer function that links occupancy to response. Equation ( 1) was poorly parameterized to yield proper identification of the transduction coefficient, hence using standard algebraic manipulations, Equation ( 1) was transformed into Equation ( 2):

E = Basal + (E max -Basal)/(1+(([L]/10ˆLogKa+1)/(10ˆLogR × [L]))ˆn), (2)
Factor n was set to 1 for all concentration/activity curves in order to improve parameter identifiability.

For a given pair of responses, E1 and E2, we systematically computed the bias between each LMW ligand and FSH. Equation ( 2) was used with Ka set to 1 for both responses when referring to FSH:

E1_FSH = Basal1 + (E max 1 -Basal1)/(1 + ([L] + 1)/(10ˆLogR1 × [L])), (3)
E2_FSH = Basal2 + (E max 2 -Basal2)/(1 + ([L] + 1)/(10ˆLogR2 × [L])), (4)
Then, we used Equation ( 2) with the same Basal and E max value, but with modified functional equilibrium dissociation constant and transduction coefficient and applied this to each LMW compound individually. Considering LMW as a generic compound belonging to the panel, equations are as follows:

E1_LMW = Basal1 + (E max 1 -Basal1)/(1 + ([L]/10ˆLogKa1_LMW + 1)/(10ˆ(LogR1 + ∆LogR1_LMW)) × [L])), (5)
E2_LMW = Basal2 + (E max 2 -Basal2)/(1 + ([L]/10ˆLogKa2_LMW + 1)/(10ˆ(LogR2 + ∆LogR1_LMW + ∆∆LogBF_LMW)) × [L])), (6)
Parameter ∆LogR1_LMW is the log ratio of the transduction coefficient between FSH and LMW for response E1, that is:

∆LogR1_LMW = LogR1_LMW -LogR1 (7)
where ∆∆LogBF_LMW is the log bias between FSH and LMW and between response E1 and E2, that is:

∆∆LogBF_LMW = ∆LogR2_LMW-∆LogR1_LMW, (8)
∆LogR2_LMW = LogR2_LMW -LogR2 (9)
Performing the fit with Equations ( 5) and ( 6) allows a better estimation of the confidence interval of the log bias ∆∆LogBF_LMW as it directly enters the equation as a free parameter.

Then, we collectively fit the Equations ( 3) and ( 4) for FSH and Equations ( 5) and ( 6) for B1, B2, B3, and T1, yielding maximum likelihood estimates for all parameter values (Basal1, Emax1, R1, Basal2, Emax2, R2, Ka1_B1, Ka1_B2, Ka1_B3, Ka1_T1, Ka1_T2, Ka2_B1, Ka2_B2, Ka2_B3, Ka2_T1, Ka2_T2, ∆LogR1_B1, ∆LogR1_B2, ∆LogR1_B3, ∆LogR1_T1, ∆LogR1_T2, ∆∆LogBF_B1, ∆∆LogBF_B2, ∆∆LogBF_B3, ∆∆LogBF_T1, ∆∆LogBF_T2) in a single global fit (1000 multi-start deterministic optimization with default setting performed in D2D). For each response, we choose an additive Gaussian noise with unknown variance as the error model.

Since multi-start deterministic optimization requires a priori bounds for each parameter, we used the interval (-5, +5) for each log bias parameter value.

The 95% parameter confidence intervals for log bias values (∆∆LogBF_B1, ∆∆LogBF_B2, ∆∆LogBF_B3, ∆∆LogBF_T1, ∆∆LogBF_T2) were determined using the profile likelihood method [START_REF] Raue | Comparison of Approaches for Parameter Identifiability Analysis of Biological Systems[END_REF]. This approach is considered more robust than the asymptotic approach and considers finite sample size as well as parameter non-identifiability due to possible parameter correlations.

In order to test the statistical significance of a non-zero value for the log bias, we used the following likelihood ratio test for each compound: we defined the null model associated with LMW ligands as the very same model associated with Equations ( 3) and ( 4) for FSH and Equations ( 5) and ( 6) for each LMW allosteric ligand, except that ∆∆LogBF_LMW was set to 0. The maximum likelihood was found for the null model using the same optimization procedure as above (1000 multi-start deterministic optimization with default setting in D2D). Then, we used the likelihood ratio test to compare the negative log ratio of the likelihood of the null model over the likelihood of the full model with a chi-square distribution with one degree of freedom. Thus, we obtained the associated p-value for the significance of a non-zero log bias value vs. a zero-log bias value, taking into account possible parameter compensations when the log bias is fixed to zero. We performed this procedure for each compound individually. Importantly, these procedures (global data fitting, profile likelihood estimate, and likelihood ratio test) were performed for every pair of responses in each cell model.

Dealing with Parameter Non-Identifiability in Bias Calculations

Interpretation of systematic bias calculations requires caution. In particular, a number of concentration/activity curves have a flat profile for which precise estimation of the transduction coefficient is impossible, yielding very large or infinite confidence intervals. Several possibilities were considered based on visual inspection of concentration/activity curves and corresponding data fitting of the operational model (Equations ( 3) and ( 4) for FSH and Equations ( 5) and ( 6) for LMW ligands).

Principal Component Analysis (PCA) and Hierarchical Clustering

To analyze the high dimensional results obtained with the bias calculations, we performed a standard Principal Component Analysis (PCA) using the R package "FactoMineR" (Ver. 3.6.2, R Foundation for Statistical Computing, Vienna, Austria). Hierarchical clustering was performed based on the PCA results [START_REF] Kassambara | Practical Guide to Cluster Analysis in R: Unsupervised Machine Learning[END_REF] using the HCPC function in the R package "FactoMineR" as well as the fviz_dend function in the R package "factoextra".

Supplementary Materials:

The following are available online at https://www.mdpi.com/article/10 .3390/ijms22189850/s1.

Figure 1 .

 1 Figure 1. Chemical structures of LMW ligands used in the study. B1, B2, B3 are Benzamide derivatives; T1 is a Thiazolidinone derivative, and T2 is an antagonist obtained by chemically modifying T1.

Figure 1 .

 1 Figure 1. Chemical structures of LMW ligands used in the study. B1, B2, B3 are Benzamide derivatives; T1 is a Thiazolidinone derivative, and T2 is an antagonist obtained by chemically modifying T1.

Figure 2 .

 2 Figure 2. Profiling of LMW ligands vs. FSH in HEK293 cells. (A) Gαs recruitment (light blue), (B) Gαq recruitment (green), (C) Gαi recruitment (purple), (D) β-arrestin 2 recruitment (grey), and (E) cAMP production were measured. HEK293 cells transfected with the FSHR and the appropriate BRET sensors (Rluc stands for renilla luciferase, YFP for yellow fluorescent protein, and EPAC for exchange protein directly activated by cAMP) were challenged with increasing doses of FSH (10 -12.5 to 10 -6.5 M) or LMW ligands (10 -10 to 10 -4 M). Concentration/activity curves were measured for the different read-outs. Signals were monitored for 30 min and normalized, considering the FSH maximal response as 100%. The area under the concentration/activity curves generated by each compound was plotted and fitted using nonlinear regression. Data are represented as mean ± SEM from at least six independent experiments.

Figure 2 .

 2 Figure 2. Profiling of LMW ligands vs. FSH in HEK293 cells. (A) Gαs recruitment (light blue), (B) Gαq recruitment (green), (C) Gαi recruitment (purple), (D) β-arrestin 2 recruitment (grey), and (E) cAMP production were measured. HEK293 cells transfected with the FSHR and the appropriate BRET sensors (Rluc stands for renilla luciferase, YFP for yellow fluorescent protein, and EPAC for exchange protein directly activated by cAMP) were challenged with increasing doses of FSH (10 -12.5to 10 -6.5 M) or LMW ligands (10 -10 to 10 -4 M). Concentration/activity curves were measured for the different read-outs. Signals were monitored for 30 min and normalized, considering the FSH maximal response as 100%. The area under the concentration/activity curves generated by each compound was plotted and fitted using non-linear regression. Data are represented as mean ± SEM from at least six independent experiments.

Figure 3 .

 3 Figure 3. Kinetics curves induced by EC80 of FSH and LMW ligands. (A) Gαs recruitment, (B) Gαq recruitment, (C) Gαi recruitment, (D) β-arrestin 2 recruitment, and (E) cAMP production were measured. HEK293 cells transfected with FSHR and appropriate BRET sensors were stimulated with EC80 of FSH or LMW ligands. Signals were monitored for 20 to 40 min, depending on the read-out. Kinetic curves corresponding to each read-out were normalized considering the signal elicited by FSH at 20 min as 100%. Kinetics curves were plotted as a function of time. Data are represented as mean ± SEM from at least six independent experiments.

Figure 3 .

 3 Figure 3. Kinetics curves induced by EC80 of FSH and LMW ligands. (A) Gαs recruitment, (B) Gαq recruitment, (C) Gαi recruitment, (D) β-arrestin 2 recruitment, and (E) cAMP production were measured. HEK293 cells transfected with FSHR and appropriate BRET sensors were stimulated with EC 80 of FSH or LMW ligands. Signals were monitored for 20 to 40 min, depending on the read-out. Kinetic curves corresponding to each read-out were normalized considering the signal elicited by FSH at 20 min as 100%. Kinetics curves were plotted as a function of time. Data are represented as mean ± SEM from at least six independent experiments.

  tration/activity curves of T1 vs. FSH at different time points. Gαs recruitment (A) and cAMP production (B) at 2, 7, and 20 min of stimulation. Efficacies and potencies were compared (C). The area under the tivity curves generated by T1 at the different endpoints considered (2, 7, and 20 min) were calculated. FSH e at each endpoint was set at 100%, and T1-induced values were normalized accordingly. Normalized values fitted using non-linear regression. Emax values were represented as the % of FSH Emax ± SEM. Potency esented as the positive logarithm of the ligand EC50 concentration ± SEM for each time considered. Statistical assessed by unpaired t-test with Welch's correction. n = 6, * p < 0.05; *** p < 0.001.

Figure 4 .

 4 Figure 4. Concentration/activity curves of T1 vs. FSH at different time points. Gαs recruitment (A) and cAMP production (B) were measured at 2, 7, and 20 min of stimulation. Efficacies and potencies were compared (C). The area under the concentration/activity curves generated by T1 at the different endpoints considered (2, 7, and 20 min) were calculated. FSH maximal response at each endpoint was set at 100%, and T1-induced values were normalized accordingly. Normalized values were plotted and fitted using non-linear regression. Emax values were represented as the % of FSH Emax ± SEM. Potency values were represented as the positive logarithm of the ligand EC 50 concentration ± SEM for each time considered. Statistical significance was assessed by unpaired t-test with Welch's correction. n = 6, * p < 0.05; *** p < 0.001.

Figure 5 .

 5 Figure 5. Activity of FSH vs. LMW ligands on downstream signaling. (A) HEK293 cells expressing the FSHR were stimulated with EC80 concentration of FSH or LMW ligands for 5 min. ERK phosphorylation was analyzed by western blotting. (B) Images were quantified and represented as normalized phosphorylated ERK/ERK1/2 values. (C) HEK293 co-expressing the FSHR and CREluciferase reporter gene were stimulated with increasing concentrations of FSH (10 -12.5 to 10 -6.5 M) or LMW ligands (10 -10 to 10 -4 M). Luciferase activities were recorded 6 h after stimulation, and signals were normalized over FSH maximal response. Values were plotted and fitted by non-linear regression. Data are represented as mean ± SEM from at least six independent experiments.

Figure 5 .

 5 Figure 5. Activity of FSH vs. LMW ligands on downstream signaling. (A) HEK293 cells expressing the FSHR were stimulated with EC 80 concentration of FSH or LMW ligands for 5 min. ERK phosphorylation was analyzed by western blotting. (B) Images were quantified and represented as normalized phosphorylated ERK/ERK1/2 values. (C) HEK293 co-expressing the FSHR and CRE-luciferase reporter gene were stimulated with increasing concentrations of FSH (10 -12.5 to 10 -6.5 M) or LMW ligands (10 -10 to 10 -4 M). Luciferase activities were recorded 6 h after stimulation, and signals were normalized over FSH maximal response. Values were plotted and fitted by non-linear regression. Data are represented as mean ± SEM from at least six independent experiments.

Figure 6 .

 6 Figure 6. Profiling of LMW ligands vs. FSH in HEK293/ΔGαs cells. (A) Gαq recruitment, (B) Gαi recruitment, (C) β-arrestin 2 recruitment, and (D) CRE reporter gene induction were measured.Transfected HEK293/ΔGαs cells were challenged with increasing doses of FSH (10 -12.5 to 10 -6.5 M) or LMW ligands (10 -10 to 10 -4 M). Signals were monitored for 30 min. Concentration/activity curves were calculated by considering the area under the curves. Values were normalized considering FSH maximal response as 100%, plotted, and fitted using non-linear regression. Data are represented as mean ± SEM from at least six independent experiments.

Figure 6 .

 6 Figure 6. Profiling of LMW ligands vs. FSH in HEK293/∆Gαs cells. (A) Gαq recruitment, (B) Gαi recruitment, (C) β-arrestin 2 recruitment, and (D) CRE reporter gene induction were measured. Transfected HEK293/∆Gαs cells were challenged with increasing doses of FSH (10 -12.5 to 10 -6.5 M) or LMW ligands (10 -10 to 10 -4 M). Signals were monitored for 30 min. Concentration/activity curves were calculated by considering the area under the curves. Values were normalized considering FSH maximal response as 100%, plotted, and fitted using non-linear regression. Data are represented as mean ± SEM from at least six independent experiments.

Figure 7 .

 7 Figure 7. Profiling of LMW ligands vs. FSH in HEK293/ΔARRB cells. (A) Gαs recruitment, (B) Gαq recruitment, (C) Gαi recruitment, (D) cAMP production, and (E) CRE reporter gene induction were measured. Transfected HEK293/ΔARRB cells were challenged with increasing doses of FSH (10 -12.5 to 10 -6.5 M) or LMW ligands (10 -10 to 10 -4 M). Signals were monitored for 30 min.

Figure 7 .

 7 Figure 7. Profiling of LMW ligands vs. FSH in HEK293/∆ARRB cells. (A) Gαs recruitment, (B) Gαq recruitment, (C) Gαi recruitment, (D) cAMP production, and (E) CRE reporter gene induction were measured. Transfected HEK293/∆ARRB cells were challenged with increasing doses of FSH (10 -12.5 to 10 -6.5 M) or LMW ligands (10 -10 to 10 -4 M). Signals were monitored for 30 min. Concentration/activity curves were calculated by considering the area under the curves. Values were normalized considering FSH maximal response as 100%, plotted, and fitted using non-linear regression. Data are represented as mean ± SEM from at least six independent experiments.

Figure 8 .

 8 Figure 8. Radial graph representation of FSH and LMW ligand-induced transduction coefficients (TC) and maximal efficacies (Emax).Transduction coefficients and maximal efficacies were obtained after curve-fitting of FSH and LMW ligand-induced concentration/activity curves and are represented as radial graphs. Radius for TC is in logarithmic scale while radius for Emax is normalized, considering FSH maximal response as 100%, and is shown in linear scale. B1-induced transduction coefficients in HEK293/ΔGαs were not identifiable, and the relative radial graph was labeled as "not determined" (nd).

Figure 8 . 24 Figure 9 .

 8249 Figure 8. Radial graph representation of FSH and LMW ligand-induced transduction coefficients (TC)and maximal efficacies (Emax). Transduction coefficients and maximal efficacies were obtained after curve-fitting of FSH and LMW ligand-induced concentration/activity curves and are represented as radial graphs. Radius for TC is in logarithmic scale while radius for Emax is normalized, considering FSH maximal response as 100%, and is shown in linear scale. B1-induced transduction coefficients in HEK293/∆Gαs were not identifiable, and the relative radial graph was labeled as "not determined" (nd). Int. J. Mol. Sci. 2021, 22, x FOR PEER REVIEW 14 of 24

Figure 9 .

 9 Figure 9. Graphic representation of the biased pharmacological behavior induced by LMW ligands. The biased pharmacological behavior of each LMW (represented in a color code, B1 = green, B2 = orange, B3 = blue and T1 = red) ligand is represented in the shape of an arrow pointing to the direction of bias for each of the pair of read-outs considered in the three cellular contexts. Statistically significant biases (full arrow) as well as observed biases (empty arrow) are shown.

Figure 9 .

 9 Figure 9. Graphic representation of the biased pharmacological behavior induced by LMW ligands. The biased pharmacological behavior of each LMW (represented in a color code, B1 = green, B2 = orange, B3 = blue and T1 = red) ligand is represented in the shape of an arrow pointing to the direction of bias for each of the pair of read-outs considered in the three cellular contexts. Statistically significant biases (full arrow) as well as observed biases (empty arrow) are shown.

Figure 10 .

 10 Figure 10. Hierarchical clustering and principal component analysis of LMW ligand-induced bias factors. Continuous lines represent the cut-off set up of a number of clusters equal to three. Clusters of ligands are highlighted with dotted lines. Note that in the PCA graph, corresponding to HEK293/ΔGαs, FSH and B1 are perfectly overlaid.

Figure 10 .

 10 Figure 10. Hierarchical clustering and principal component analysis of LMW ligand-induced bias factors. Continuous lines represent the cut-off set up of a number of clusters equal to three. Clusters of ligands are highlighted with dotted lines. Note that in the PCA graph, corresponding to HEK293/∆Gαs, FSH and B1 are perfectly overlaid.

Author Contributions:

  Participation in research design: F.D.P., M.A.A., F.J.-A., L.P.P., P.C., A.P., L.C., M.S., A.C.H., R.Y. and E.R.; conducted experiments: F.D.P., M.A.A., R.B., S.S., N.G., F.L., C.G. and P.R.; development of new reagents or analytical tools: J.L., A.I., N.J., M.-C.V.-M., L.P.P., S.G.N., H.N.Y. and S.S.P.; funding and supervision of the study: P.C., R.Y. and E.R.; writing of the manuscript: F.D.P., M.A.A., R.Y. and E.R. All authors reviewed and edited the manuscript. All authors have read and agreed to the published version of the manuscript. Funding: This publication was funded with support from the French National Research Agency under the program "Investissements d'avenir" Grant Agreement LabEx MabImprove: ANR-10-LABX-53; "ARD2020 Biomédicaments" and "APR-IR MODUPHAC" grants from Région Centre Val de Loire. F.D.P. was the recipient of a doctoral fellowship from INRAE and Région Centre Val de Loire. R.B. was the recipient of an Erasmus plus fellowship from the University of Modena and Reggio Emilia and the European Union. M.A.A., F.J.-A. and M.S. were LE STUDIUM RESEARCH FELLOWS, Loire Valley Institute for Advanced Studies, Orléans & Tours, France, receiving funding from the European Union's Horizon 2020 research and innovation programme under the Marie Skłodowska-Curie grant agreement No 665790. A.I. was funded by the PRIME 18gm5910013 and the LEAP 18gm0010004 from the Japan Agency for Medical Research and Development (AMED) and KAKENHI 17K08264 by the Japan Society for the Promotion of Science (JSPS). A.C.H. was supported by a Biotechnology and Biological Sciences Research Council (BBSRC) project grant (BB/S001565/1) and a Society for Endocrinology Early Career Grant to S.S.

Data Availability Statement: Not applicable.

Conflicts of Interest: S.G.N., H.N.Y. and S.S.P. were founders and employees of TocopheRx Inc.