Mixed Analog/Digital Gonadotrope Biosynthetic Response to Gonadotropin-releasing Hormone

Mathematical models in GnRH research

Mathematical modelling is an indispensable tool in modern biosciences, enabling quantitative analysis and integration of biological data, transparent formulation of our understanding of complex biological systems, and efficient experimental design based on model predictions. This review article provides an overview of the impact that mathematical models had on GnRH research. Indeed, over the last 20 years mathematical modelling has been used to describe and explore the physiology of the GnRH neuron, the mechanisms underlying GnRH pulsatile secretion, and GnRH signalling to the pituitary. Importantly, these models have contributed to GnRH research via novel hypotheses and predictions regarding the bursting behaviour of the GnRH neuron, the role of kisspeptin neurons in the emergence of pulsatile GnRH dynamics, and the decoding of GnRH signals by biochemical signalling networks. We envisage that with the advent of novel experimental technologies, mathematical modelling will have an even greater role to play in our endeavour to understand the complex spatiotemporal dynamics underlying the reproductive neuroendocrine system.

Operating regimes in a single enzymatic cascade at ensemble-level

Single enzymatic cascade, ubiquitously found in cellular signaling networks, is a phosphorylation-dephosphorylation reaction cycle causing a transition between inactive and active states of a protein catalysed by kinase and phosphatase, respectively. Steady-state information processing ability of such a cycle (e.g., MAPK cascade) has been classified into four qualitatively different operating regimes, viz., hyperbolic (H), signal-transducing (ST), threshold-hyperbolic (TH) and ultrasensitive (U). These four regimes represent qualitatively different dose-response curves, that is, relationship between concentrations of input kinase (e.g., pMEK) and response activated protein (e.g., pERK). Regimes were identified using a deterministic model accounting for population-averaged behavior only. Operating regimes can be strongly influenced by the inherently present cell-to-cell variability in an ensemble of cells which is captured in the form of pMEK and pERK distributions using reporter-based single-cell experimentation. In this study, we show that such experimentally acquired snapshot pMEK and pERK distribution data of a single MAPK cascade can be directly used to infer the underlying operating regime even in the absence of a dose-response curve. This deduction is possible primarily due to the presence of a monotonic relationship between experimental observables R_IQR, ratio of the inter-quartile range of the pERK and pMEK distribution pairs and R_M, ratio of the medians of the distribution pair. We demonstrate this relationship by systematic analysis of a quasi-steady state approximated model superimposed with an input gamma distribution constrained by the stimulus strength specific pMEK distribution measured on Jurkat-T cells stimulated with PMA. As a first, we show that introduction of cell-to-cell variability only in the upstream kinase achieved by superimposition of an appropriate input pMEK distribution on the dose-response curve can predict bimodal response pERK distribution in ST regime. Implementation of the proposed method on the input-response distribution pair obtained in stimulated Jurkat-T cells revealed that while low-dosage PMA stimulation preserves the H regime observed in resting cells, high-dosage causes H to ST regime transition.

Single-cell stabilization method identifies gonadotrope transcriptional dynamics and pituitary cell type heterogeneity

Immediate-early response genes (IEGs) are rapidly and transiently induced following an extracellular signal. Elucidating the IEG response patterns in single cells (SCs) requires assaying large numbers of timed samples at high accuracy while minimizing handling effects. To achieve this, we developed and validated RNA stabilization Buffer for Examination of Single-cell Transcriptomes (RNA-Best), a versatile single-step cell and tissue preservation protocol that stabilizes RNA in intact SCs without perturbing transcription patterns. We characterize for the first time SC heterogeneity in IEG responses to pulsatile gonadotropin-releasing hormone (GnRH) stimuli in pituitary gonadotrope cells. Our study identifies a gene-specific hierarchical pattern of all-or-none transcript induction elicited by increasing concentrations of GnRH. This quantal pattern of gene activation raises the possibility that IEG activation, when accurately resolved at the SC level, may be mediated by gene bits that behave as pure binary switches.

Sharing of Phosphatases Promotes Response Plasticity in Phosphorylation Cascades

Sharing of positive or negative regulators between multiple targets is frequently observed in cellular signaling cascades. For instance, phosphatase sharing between multiple kinases is ubiquitous within the MAPK pathway. Here we investigate how such phosphatase sharing could shape robustness and evolvability of the phosphorylation cascade. Through modeling and evolutionary simulations, we demonstrate that 1) phosphatase sharing dramatically increases robustness of a bistable MAPK response, and 2) phosphatase-sharing cascades evolve faster than nonsharing cascades. This faster evolution is particularly pronounced when evolving from a monostable toward a bistable phenotype, whereas the transition speed of a population from a bistable to monostable response is not affected by phosphatase sharing. This property may enable the phosphatase-sharing design to adapt better in a changing environment. Analysis of the respective mutational landscapes reveal that phosphatase sharing reduces the number of limiting mutations required for transition from monostable to bistable responses, hence facilitating a faster transition to such response types. Taken together, using MAPK cascade as an example, our study offers a general theoretical framework to explore robustness and evolutionary plasticity of signal transduction cascades.

Gonadotropin-releasing hormone signaling: An information theoretic approach

Gonadotropin-releasing hormone (GnRH) is a peptide hormone that mediates central control of reproduction, acting via G-protein coupled receptors that are primarily G_q coupled and mediate GnRH effects on the synthesis and secretion of luteinizing hormone and follicle-stimulating hormone. A great deal is known about the GnRH receptor signaling network but GnRH is secreted in short pulses and much less is known about how gonadotropes decode this pulsatile signal. Similarly, single cell measures reveal considerable cell-cell heterogeneity in responses to GnRH but the impact of this variability on signaling is largely unknown. Ordinary differential equation-based mathematical models have been used to explore the decoding of pulse dynamics and information theory-derived statistical measures are increasingly used to address the influence of cell-cell variability on the amount of information transferred by signaling pathways. Here, we describe both approaches for GnRH signaling, with emphasis on novel insights gained from the information theoretic approach and on the fundamental question of why GnRH is secreted in pulses.

Regulatory Architecture of the LβT2 Gonadotrope Cell Underlying the Response to Gonadotropin-Releasing Hormone

The LβT2 mouse pituitary cell line has many characteristics of a mature gonadotrope and is a widely used model system for studying the developmental processes and the response to gonadotropin-releasing hormone (GnRH). The global epigenetic landscape, which contributes to cell-specific gene regulatory mechanisms, and the single-cell transcriptome response variation of LβT2 cells have not been previously investigated. Here, we integrate the transcriptome and genome-wide chromatin accessibility state of LβT2 cells during GnRH stimulation. In addition, we examine cell-to-cell variability in the transcriptional response to GnRH using Gel bead-in-Emulsion Drop-seq technology. Analysis of a bulk RNA-seq data set obtained 45 min after exposure to either GnRH or vehicle identified 112 transcripts that were regulated >4-fold by GnRH (FDR < 0.05). The top regulated transcripts constitute, as determined by Bayesian massive public data integration analysis, a human pituitary-relevant coordinated gene program. Chromatin accessibility [assay for transposase-accessible chromatin with high-throughput sequencing (ATAC-seq)] data sets generated from GnRH-treated LβT2 cells identified more than 58,000 open chromatin regions, some containing notches consistent with bound transcription factor footprints. The study of the most prominent open regions showed that 75% were in transcriptionally active promoters or introns, supporting their involvement in active transcription. Lhb, Cga, and Egr1 showed significantly open chromatin over their promoters. While Fshb was closed over its promoter, several discrete significantly open regions were found at -40 to -90 kb, which may represent novel upstream enhancers. Chromatin accessibility determined by ATAC-seq was associated with high levels of gene expression determined by RNA-seq. We obtained high-quality single-cell Gel bead-in-Emulsion Drop-seq transcriptome data, with an average of >4,000 expressed genes/cell, from 1,992 vehicle- and 1,889 GnRH-treated cells. While the individual cell expression patterns showed high cell-to-cell variation, representing both biological and measurement variation, the average expression patterns correlated well with bulk RNA-seq data. Computational assignment of each cell to its precise cell cycle phase showed that the response to GnRH was unaffected by cell cycle. To our knowledge, this study represents the first genome-wide epigenetic and single-cell transcriptomic characterization of this important gonadotrope model. The data have been deposited publicly and should provide a resource for hypothesis generation and further study.

Multifaceted Targeting of the Chromatin Mediates Gonadotropin-Releasing Hormone Effects on Gene Expression in the Gonadotrope

Gonadotropin-releasing hormone (GnRH) stimulates the expression of multiple genes in the pituitary gonadotropes, most notably to induce synthesis of the gonadotropins, luteinizing hormone (LH), and follicle-stimulating hormone (FSH), but also to ensure the appropriate functioning of these cells at the center of the mammalian reproductive endocrine axis. Aside from the activation of gene-specific transcription factors, GnRH stimulates through its membrane-bound receptor, alterations in the chromatin that facilitate transcription of its target genes. These include changes in the histone and DNA modifications, nucleosome positioning, and chromatin packaging at the regulatory regions of each gene. The requirements for each of these events vary according to the DNA sequence which determines the basal chromatin packaging at the regulatory regions. Despite considerable progress in this field in recent years, we are only beginning to understand some of the complexities involved in the role and regulation of this chromatin structure, including new modifications, extensive cross talk, histone variants, and the actions of distal enhancers and non-coding RNAs. This short review aims to integrate the latest findings on GnRH-induced alterations in the chromatin of its target genes, which indicate multiple and diverse actions. Understanding these processes is illuminating not only in the context of the activation of these hormones during the reproductive life span but may also reveal how aberrant epigenetic regulation of these genes leads to sub-fertility.

Exploring Dynamics and Noise in Gonadotropin-Releasing Hormone (GnRH) Signaling

Gonadotropin-releasing hormone (GnRH) acts via G-protein coupled receptors on pituitary gonadotropes. These are G_q-coupled receptors that mediate acute effects of GnRH on the exocytotic secretion of luteinizing hormone (LH) and follicle-stimulating hormone (FSH), as well as the chronic regulation of their synthesis. FSH and LH control steroidogenesis and gametogenesis in the gonads so GnRH mediates control of reproduction by the central nervous system. GnRH is secreted in short pulses and the effects of GnRH on its target cells are dependent on the dynamics of these pulses. Here we provide a brief overview of the signaling network activated by GnRH with emphasis on the use of high content imaging for their examination. We also describe computational approaches that we have used to simulate GnRH signaling in order to explore dynamics, noise, and information transfer in this system.

Mathematical modeling of gonadotropin-releasing hormone signaling

Gonadotropin-releasing hormone (GnRH) acts via G-protein coupled receptors on pituitary gonadotropes to control reproduction. These are G_q-coupled receptors that mediate acute effects of GnRH on the exocytotic secretion of luteinizing hormone (LH) and follicle-stimulating hormone (FSH), as well as the chronic regulation of their synthesis. GnRH is secreted in short pulses and GnRH effects on its target cells are dependent upon the dynamics of these pulses. Here we overview GnRH receptors and their signaling network, placing emphasis on pulsatile signaling, and how mechanistic mathematical models and an information theoretic approach have helped further this field.

Information Transfer via Gonadotropin-Releasing Hormone Receptors to ERK and NFAT: Sensing GnRH and Sensing Dynamics

Information theoretic approaches can be used to quantify information transfer via cell signaling networks. In this study, we do so for gonadotropin-releasing hormone (GnRH) activation of extracellular signal-regulated kinase (ERK) and nuclear factor of activated T cells (NFAT) in large numbers of individual fixed LβT2 and HeLa cells. Information transfer, measured by mutual information between GnRH and ERK or NFAT, was <1 bit (despite 3-bit system inputs). It was increased by sensing both ERK and NFAT, but the increase was <50%. In live cells, information transfer via GnRH receptors to NFAT was also <1 bit and was increased by consideration of response trajectory, but the increase was <10%. GnRH secretion is pulsatile, so we explored information gained by sensing a second pulse, developing a model of GnRH signaling to NFAT with variability introduced by allowing effectors to fluctuate. Simulations revealed that when cell–cell variability reflects rapidly fluctuating effector levels, additional information is gained by sensing two GnRH pulses, but where it is due to slowly fluctuating effectors, responses in one pulse are predictive of those in another, so little information is gained from sensing both. Wet laboratory experiments revealed that the latter scenario holds true for GnRH signaling; within the timescale of our experiments (1 to 2 hours), cell–cell variability in the NFAT pathway remains relatively constant, so trajectories are reproducible from pulse to pulse. Accordingly, joint sensing, sensing of response trajectories, and sensing of repeated pulses can all increase information transfer via GnRH receptors, but in each case the increase is small.

Reactive Oxygen Species Link Gonadotropin-Releasing Hormone Receptor Signaling Cascades in the Gonadotrope

Biological rhythms lie at the center of regulatory schemes that control many aspects of living systems. At the cellular level, meaningful responses to external stimuli depend on propagation and quenching of a signal to maintain vigilance for subsequent stimulation or changes that serve to shape and modulate the response. The hypothalamus-pituitary-gonad endocrine axis that controls reproductive development and function relies on control through rhythmic stimulation. Central to this axis is the pulsatile stimulation of the gonadotropes by hypothalamic neurons through episodic release of the neuropeptide gonadotropin-releasing hormone. Alterations in pulsatile stimulation of the gonadotropes result in differential synthesis and secretion of the gonadotropins LH and FSH and changes in the expression of their respective hormone subunit genes. The requirement to amplify signals arising from activation of the gonadotropin-releasing hormone (GnRH) receptor and to rapidly quench the resultant signal to preserve an adaptive response suggests the need for rapid activation and feedback control operating at the level of intracellular signaling. Emerging data suggest that reactive oxygen species (ROS) can fulfill this role in the GnRH receptor signaling through activation of MAP kinase signaling cascades, control of negative feedback, and participation in the secretory process. Results obtained in gonadotrope cell lines or other cell models indicate that ROS can participate in each of these regulatory cascades. We discuss the potential advantage of reactive oxygen signaling for modulating the gonadotrope response to GnRH stimulation and the potential mechanisms for this action. These observations suggest further targets of study for regulation in the gonadotrope.

Influence of gene copy number on self-regulated gene expression

Using an analytically solvable stochastic model, we study the properties of a simple genetic circuit consisting of multiple copies of a self-regulating gene. We analyse how the variation in gene copy number and the mutations changing the auto-regulation strength affect the steady-state distribution of protein concentration. We predict that one-reporter assay, an experimental method where the extrinsic noise level is inferred from the comparison of expression variance of a single and duplicated reporter gene, may give an incorrect estimation of the extrinsic noise contribution when applied to self-regulating genes. We also show that an imperfect duplication of an auto-activated gene, changing the regulation strength of one of the copies, may lead to a hybrid, binary+graded response of these genes to external signal. The analysis of relative changes in mean gene expression before and after duplication suggests that evolutionary accumulation of gene duplications may, at a given mean burst size, non-trivially depend on the inherent noisiness of a given gene, quantified by the inverse of the maximal mean frequency of bursts. Moreover, we find that the dependence of gene expression noise on gene copy number and auto-regulation strength may qualitatively differ, e.g. in monotonicity, depending on whether the noise is measured by Fano factor or coefficient of variation. Thus, experimentally-based hypotheses linking gene expression noise and evolutionary optimisation in the context of gene copy number variation may be ambiguous as they are dependent on the particular function chosen to quantify noise.

Molecular mechanisms of gonadotropin-releasing hormone signaling: integrating cyclic nucleotides into the network

Gonadotropin-releasing hormone (GnRH) is the primary regulator of mammalian reproductive function in both males and females. It acts via G-protein coupled receptors on gonadotropes to stimulate synthesis and secretion of the gonadotropin hormones luteinizing hormone and follicle-stimulating hormone. These receptors couple primarily via G-proteins of the Gq/ll family, driving activation of phospholipases C and mediating GnRH effects on gonadotropin synthesis and secretion. There is also good evidence that GnRH causes activation of other heterotrimeric G-proteins (Gs and Gi) with consequent effects on cyclic AMP production, as well as for effects on the soluble and particulate guanylyl cyclases that generate cGMP. Here we provide an overview of these pathways. We emphasize mechanisms underpinning pulsatile hormone signaling and the possible interplay of GnRH and autocrine or paracrine regulatory mechanisms in control of cyclic nucleotide signaling.

Signaling to extracellular signal-regulated kinase from ErbB1 kinase and protein kinase C: feedback, heterogeneity, and gating

Many extracellular signals act via the Raf/MEK/ERK cascade in which kinetics, cell-cell variability, and sensitivity of the ERK response can all influence cell fate. Here we used automated microscopy to explore the effects of ERK-mediated negative feedback on these attributes in cells expressing endogenous ERK or ERK2-GFP reporters. We studied acute rather than chronic stimulation with either epidermal growth factor (ErbB1 activation) or phorbol 12,13-dibutyrate (PKC activation). In unstimulated cells, ERK-mediated negative feedback reduced the population-average and cell-cell variability of the level of activated ppERK and increased its robustness to changes in ERK expression. In stimulated cells, negative feedback (evident between 5 min and 4 h) also reduced average levels and variability of phosphorylated ERK (ppERK) without altering the “gradedness” or sensitivity of the response. Binning cells according to total ERK expression revealed, strikingly, that maximal ppERK responses initially occur at submaximal ERK levels and that this non-monotonic relationship changes to an increasing, monotonic one within 15 min. These phenomena occur in HeLa cells and MCF7 breast cancer cells and in the presence and absence of ERK-mediated negative feedback. They were best modeled assuming distributive (rather than processive) activation. Thus, we have uncovered a novel, time-dependent change in the relationship between total ERK and ppERK levels that persists without negative feedback. This change makes acute response kinetics dependent on ERK level and provides a “gating” or control mechanism in which the interplay between stimulus duration and the distribution of ERK expression across cells could modulate the proportion of cells that respond to stimulation.

Interactions of the GnRH receptor with heterotrimeric G proteins

G protein-coupled receptors (GPCRs) are the largest class of integral membrane protein receptors in the human genome. We examined here the reports whether the GnRH receptor (GnRHR) interacts with a single or multiple types of G proteins. It seems that the GnRHR, as other GPCRs, alternates between various conformations and is stabilized by its ligands, other modulators and intracellular partners in selective conformations culminating in coupling with a single type or multiple G proteins in a cell- and context-specific manner.

Stat5 signaling specifies basal versus stress erythropoietic responses through distinct binary and graded dynamic modalities

Stat5 signaling in erythroblasts can assume either a binary, low-intensity form, essential for basal erythropoiesis, or a graded, high-intensity response, restricted to early erythroblasts and to erythropoietic stress.

Erythropoietin (Epo)-induced Stat5 phosphorylation (p-Stat5) is essential for both basal erythropoiesis and for its acceleration during hypoxic stress. A key challenge lies in understanding how Stat5 signaling elicits distinct functions during basal and stress erythropoiesis. Here we asked whether these distinct functions might be specified by the dynamic behavior of the Stat5 signal. We used flow cytometry to analyze Stat5 phosphorylation dynamics in primary erythropoietic tissue in vivo and in vitro, identifying two signaling modalities. In later (basophilic) erythroblasts, Epo stimulation triggers a low intensity but decisive, binary (digital) p-Stat5 signal. In early erythroblasts the binary signal is superseded by a high-intensity graded (analog) p-Stat5 response. We elucidated the biological functions of binary and graded Stat5 signaling using the EpoR-HM mice, which express a “knocked-in” EpoR mutant lacking cytoplasmic phosphotyrosines. Strikingly, EpoR-HM mice are restricted to the binary signaling mode, which rescues these mice from fatal perinatal anemia by promoting binary survival decisions in erythroblasts. However, the absence of the graded p-Stat5 response in the EpoR-HM mice prevents them from accelerating red cell production in response to stress, including a failure to upregulate the transferrin receptor, which we show is a novel stress target. We found that Stat5 protein levels decline with erythroblast differentiation, governing the transition from high-intensity graded signaling in early erythroblasts to low-intensity binary signaling in later erythroblasts. Thus, using exogenous Stat5, we converted later erythroblasts into high-intensity graded signal transducers capable of eliciting a downstream stress response. Unlike the Stat5 protein, EpoR expression in erythroblasts does not limit the Stat5 signaling response, a non-Michaelian paradigm with therapeutic implications in myeloproliferative disease. Our findings show how the binary and graded modalities combine to generate high-fidelity Stat5 signaling over the entire basal and stress Epo range. They suggest that dynamic behavior may encode information during STAT signal transduction.

Hormone signaling through the erythropoietin (Epo) pathway is required both for the continuous replacement of red blood cells (RBCs) that are lost through aging (a process known as "basal erythropoiesis") and to boost tissue oxygen when bleeding, in anemia or at high altitude ("stress erythropoiesis"). A key challenge lies in understanding how extracellular Epo concentration is translated into different intracellular signals that promote transcription of proteins that are specific to basal versus stress erythropoiesis. Binding of Epo to its receptor EpoR on the surface of an erythroblast (the precursors of RBCs) triggers the addition of phosphates to a target protein Stat5; the phosphorylated Stat5 becomes activated and induces transcription. We show that the dynamic properties of the Stat5 activation signal convey additional information that specifies either basal or stress responses. During basal conditions, the Stat5 signal is low and binary in nature—an on/off switch-like response. Stress, on the other hand, triggers a distinct Stat5 response consisting of a highintensity signal that increases in a graded fashion with rising Epo concentration. We found that a mouse bearing a truncated EpoR is restricted to the low-intensity binary Stat5 signal and correspondingly fails to initiate stress erythropoiesis. Ultimately, it is the Stat5 protein level in erythroblasts that determines their ability to generate the high-intensity graded Stat5 signal in response to high Epo. These findings have therapeutic potential: targeting Stat5's high-intensity graded signal may inhibit its aberrant function in blood cell cancers without affecting its important binary response in normal cells.

Mechanisms of GnRH-induced extracellular signal-regulated kinase nuclear localization

Gonadotropin-releasing hormone receptors (GnRHR) mediate activation and nuclear translocation of the extracellular signal regulated kinases 1 and 2 (ERK) by phosphorylation on the TEY motif. This is necessary for GnRH to initiate transcriptional programmes controlling fertility, but mechanisms that govern ERK targeting are unclear. Using automated microscopy to explore ERK regulation in single cells, we find that GnRHR activation induces marked redistribution of ERK to the nucleus and that this effect can be uncoupled from the level of TEY phosphorylation of ERK. Thus, 5 min stimulation with 100 nM GnRH increased phospho-ERK levels (from 89 ± 34 to 555 ± 45 arbitrary fluorescence units) and increased the nuclear:cytoplasmic (N:C) ERK ratio (from 1.36 ± 0.06 to 2.16 ± 0.05) in the whole cell population, but it also significantly increased N:C ERK in cells binned according to phospho-ERK levels. This phosphorylation unattributable component of the ERK translocation response occurs at a broad range of GnRHR expression levels, in the presence of tyrosine phosphatase and protein synthesis inhibitors, and in ERK mutants unable to undergo catalytic activation. It also occurred in mutants incapable of binding the DEF (docking site for ERK, F/Y-X-F/Y-P) domains found in many ERK binding partners. It was however, reduced by MEK or PKC inhibition and by mutations preventing TEY phosphorylation or that abrogate ERK binding to D (docking) domain partners. We therefore show that TEY phosphorylation of ERK is necessary, but not sufficient for the full nuclear localization response. We further show that this "phosphorylation unattributable" component of GnRH-mediated ERK nuclear translocation requires both PKC activity and association with partner proteins via the D-domain.

Decoding neurohormone pulse frequency by convergent signalling modules

GnRH (gonadotropin-releasing hormone) mediates control of reproduction. It is secreted in pulses and acts via intracellular effectors to activate gene expression. Submaximal GnRH pulse frequency can elicit maximal responses, yielding bell-shaped frequency-response curves characteristic of genuine frequency decoders. GnRH frequency decoding is therapeutically important (pulsatile GnRH can drive ovulation in assisted reproduction, whereas sustained activation can treat breast and prostate cancers), but the mechanisms are unknown. In the present paper, we review recent work in this area, placing emphasis on the regulation of transcription, and showing how mathematical modelling of GnRH effects on two effectors [ERK (extracellular-signal-regulated kinase) and NFAT (nuclear factor of activated T-cells)] reveals the potential for genuine frequency decoding as an emergent feature of the GnRH signalling network, rather than an intrinsic feature of a given protein or pathway within it.

Decoding GnRH neurohormone pulse frequency by convergent signalling modules

Gonadotropin-releasing hormone (GnRH) mediates control of reproduction. It is secreted in pulses and acts via intracellular effectors to activate gonadotrophin secretion and gene expression. Sub-maximal GnRH pulse frequency can elicit maximal responses, yielding bell-shaped frequency-response curves characteristic of genuine frequency decoders. GnRH frequency decoding is therapeutically important (pulsatile GnRH can drive ovulation in assisted reproduction whereas sustained activation can treat breast and prostate cancers), but the mechanisms are unknown. Here, we consider the possibility that it is due to convergence of distinct pulsatile signals at the transcriptome. We develop a model that mirrors wet-laboratory data for activation and nuclear translocation of GnRH effectors (extracellular signal regulated kinase and nuclear factors of activated T-cells) and incorporates transcription. The model predicts genuine frequency decoding when two transcription factors (TFs) converge at a cooperative gate, and shows how optimal pulse frequency could reflect TF activation kinetics and affinities. Importantly, this behaviour is revealed as an emergent feature of the network, rather than an intrinsic feature of a given protein or pathway, and since such network topology is extremely common, may well be widespread in biological systems.

Characterization of a MAPK scaffolding protein logic gate in gonadotropes

In the pituitary gonadotropes, both protein kinase C (PKC) and MAPK/ERK signaling cascades are activated by GnRH. Phosphoprotein-enriched in astrocytes 15 (PEA-15) is a cytosolic ERK scaffolding protein, which is expressed in LβT2 gonadotrope cells. Pharmacological inhibition of PKC and small interfering RNA-mediated silencing of Gαq/11 revealed that GnRH induces accumulation of phosphorylated PEA-15 in a PKC-dependent manner. To investigate the potential role of PEA-15 in GnRH signaling, we examined the regulation of ERK subcellular localization and the activation of ribosomal S6 kinase, a substrate of ERK. Results obtained by cellular fractionation/Western blot analysis and immunohistochemistry revealed that GnRH-induced accumulation of phosphorylated ERK in the nucleus was attenuated when PEA-15 expression was reduced. Conversely, in the absence of GnRH stimulation, PEA-15 anchors ERK in the cytosol. Our data suggest that GnRH-induced nuclear translocation of ERK requires its release from PEA-15, which occurs upon PEA-15 phosphorylation by PKC. Additional gene-silencing experiments in GnRH-stimulated cells demonstrated that ribosomal S6 kinase activation was dependent on both PEA-15 and PKC. Furthermore, small interfering RNA-mediated knockdown of PEA-15 caused a reduction in GnRH-stimulated expression of early response genes Egr2 and c-Jun, as well as gonadotropin FSHβ-subunit gene expression. PEA-15 knockdown increased LHβ and common α-glycoprotein subunit mRNAs, suggesting a possible role in differential regulation of gonadotropin subunit gene expression. We propose that PEA-15 represents a novel point of convergence of the PKC and MAPK/ERK pathways under GnRH stimulation. PKC, ERK, and PEA-15 form an AND logic gate that shapes the response of the gonadotrope cell to GnRH.

Encoding and decoding mechanisms of pulsatile hormone secretion

Ultradian pulsatile hormone secretion underlies the activity of most neuroendocrine systems, including the hypothalamic-pituitary adrenal (HPA) and gonadal (HPG) axes, and this pulsatile mode of signalling permits the encoding of information through both amplitude and frequency modulation. In the HPA axis, glucocorticoid pulse amplitude increases in anticipation of waking, and, in the HPG axis, changing gonadotrophin-releasing hormone pulse frequency is the primary means by which the body alters its reproductive status during development (i.e. puberty). The prevalence of hormone pulsatility raises two crucial questions: how are ultradian pulses encoded (or generated) by these systems, and how are these pulses decoded (or interpreted) at their target sites? We have looked at mechanisms within the HPA axis responsible for encoding the pulsatile mode of glucocorticoid signalling that we observe in vivo. We review evidence regarding the 'hypothalamic pulse generator' hypothesis, and describe an alternative model for pulse generation, which involves steroid feedback-dependent endogenous rhythmic activity throughout the HPA axis. We consider the decoding of hormone pulsatility by taking the HPG axis as a model system and focussing on molecular mechanisms of frequency decoding by pituitary gonadotrophs.

Modulation of gonadotropin-releasing hormone-induced extracellular signal-regulated kinase activation by dual-specificity protein phosphatase 1 in LbetaT2 gonadotropes

As the regulator of pituitary reproductive hormone synthesis, the hypothalamic neuropeptide GnRH is the central regulator of reproduction. A hallmark of GnRH action is the differential control of gene expression in pituitary gonadotropes through varied pulsatile stimulation. Among other signaling events, GnRH activation of the ERK family of MAPKs plays a significant role in the transcriptional regulation of the luteinizing hormone β-subunit gene and regulation of cap-dependent translation. We evaluated the ERK response to different GnRH pulse amplitudes in the gonadotrope cell line LβT2. We found that low-amplitude stimulation with GnRH invokes a rapid and transient ERK activation, whereas high-amplitude stimulation invokes a prolonged activation specifically in the cytoplasm fraction of LβT2 cells. Nuclear and cytoplasmic targets of ERK, Ets-like gene 1, and eukaryotic initiation factor 4E, respectively, are similarly activated. Feedback control of ERK activation occurs mainly through the dual-specificity protein phosphatases (DUSPs). DUSP1 is localized to the nucleus in LβT2 cells but DUSP4, another member implicated in GnRH feedback, exists in both the nucleus and cytoplasm. Manipulation of nuclear DUSP activity through overexpression or knockdown of Dusp1 modulates the ERK response to low and high GnRH pulse amplitudes and activation of the Lhb promoter. Dusp1 overexpression abolishes sustained ERK activation and inhibits Lhb promoter activity induced by high amplitude pulses. Conversely, Dusp1 knockdown enhances ERK activation by low-amplitude stimulation and increases stimulation of Lhb promoter activity. We conclude that DUSP1 feedback activity modulates ERK activation and the transcriptional response to GnRH.

Plato's cave algorithm: inferring functional signaling networks from early gene expression shadows

Improving the ability to reverse engineer biochemical networks is a major goal of systems biology. Lesions in signaling networks lead to alterations in gene expression, which in principle should allow network reconstruction. However, the information about the activity levels of signaling proteins conveyed in overall gene expression is limited by the complexity of gene expression dynamics and of regulatory network topology. Two observations provide the basis for overcoming this limitation: a. genes induced without de-novo protein synthesis (early genes) show a linear accumulation of product in the first hour after the change in the cell's state; b. The signaling components in the network largely function in the linear range of their stimulus-response curves. Therefore, unlike most genes or most time points, expression profiles of early genes at an early time point provide direct biochemical assays that represent the activity levels of upstream signaling components. Such expression data provide the basis for an efficient algorithm (Plato's Cave algorithm; PLACA) to reverse engineer functional signaling networks. Unlike conventional reverse engineering algorithms that use steady state values, PLACA uses stimulated early gene expression measurements associated with systematic perturbations of signaling components, without measuring the signaling components themselves. Besides the reverse engineered network, PLACA also identifies the genes detecting the functional interaction, thereby facilitating validation of the predicted functional network. Using simulated datasets, the algorithm is shown to be robust to experimental noise. Using experimental data obtained from gonadotropes, PLACA reverse engineered the interaction network of six perturbed signaling components. The network recapitulated many known interactions and identified novel functional interactions that were validated by further experiment. PLACA uses the results of experiments that are feasible for any signaling network to predict the functional topology of the network and to identify novel relationships.

Pulsatile and sustained gonadotropin-releasing hormone (GnRH) receptor signaling: does the ERK signaling pathway decode GnRH pulse frequency?

Gonadotropin-releasing hormone (GnRH) acts via G-protein-coupled receptors on gonadotrophs to stimulate synthesis and secretion of luteinizing hormone and follicle-stimulating hormone. It is secreted in pulses, and its effects depend on pulse frequency, but decoding mechanisms are unknown. Here we have used an extracellular signal regulated kinase-green fluorescent protein (ERK2-GFP) reporter to monitor GnRH signaling. GnRH caused dose-dependent ERK2-GFP translocation to the nucleus, providing a live-cell readout for activation. Pulsatile GnRH caused dose- and frequency-dependent ERK2-GFP translocation. These responses were rapid and transient, showed only digital tracking, and did not desensitize under any condition tested (dose, frequency, and receptor number varied). We also tested for the effects of cycloheximide (to prevent induction of nuclear-inducible MAPK phosphatases) and used GFP fusions containing ERK mutations (D319N, which prevents docking domain-dependent binding to MAPK phosphatases, and K52R, which prevents catalytic activity). These manipulations had little or no effect on the translocation responses, arguing against a role for MAPK phosphatases or ERK-mediated feedback in shaping ERK activation during pulsatile stimulation. GnRH also caused dose- and frequency-dependent activation of the alpha-gonadotropin subunit-, luteinizing hormone beta-, and follicle-stimulating hormone beta- luciferase reporters, and the latter response was inhibited by ERK1/2 knockdown. Moreover, GnRH caused frequency-dependent activation of an Egr1-luciferase reporter, but the response was proportional to cumulative pulse duration. Our data suggest that frequency decoding is not due to negative feedback shaping ERK signaling in this model.

Frequency-dependent regulation of follicle-stimulating hormone beta by pulsatile gonadotropin-releasing hormone is mediated by functional antagonism of bZIP transcription factors

Oscillatory synthesis and secretion of the gonadotropins, follicle-stimulating hormone (FSH) and luteinizing hormone (LH), under the control of pulsatile hypothalamic gonadotropin-releasing hormone (GnRH), is essential for normal reproductive development and fertility. The molecular mechanisms by which various patterns of pulsatile GnRH regulate gonadotrope responsiveness remain poorly understood. In contrast to the alpha and LH beta subunit genes, FSH beta subunit transcription is preferentially stimulated at low rather than high frequencies of pulsatile GnRH. In this study, mutation of a cyclic AMP response element (CRE) within the FSH beta promoter resulted in the loss of preferential GnRH stimulation at low pulse frequencies. We hypothesized that high GnRH pulse frequencies might stimulate a transcriptional repressor(s) to attenuate the action of CRE binding protein (CREB) and show that inducible cAMP early repressor (ICER) fulfills such a role. ICER was not detected under basal conditions, but pulsatile GnRH stimulated ICER to a greater extent at high than at low pulse frequencies. ICER binds to the FSH beta CRE site to reduce CREB occupation and abrogates both maximal GnRH stimulation and GnRH pulse frequency-dependent effects on FSH beta transcription. These data suggest that ICER production antagonizes the stimulatory action of CREB to attenuate FSH beta transcription at high GnRH pulse frequencies, thereby playing a critical role in regulating cyclic reproductive function.

The biology of gonadotroph regulation

PURPOSE OF REVIEW

To discuss recent progress in our understanding of pituitary gonadotroph development and gonadotropin gene regulation, with an emphasis on differential luteinizing hormone (LH) and follicle-stimulating hormone (FSH) secretion and subunit synthesis, and the implications this may have on female reproductive health.

RECENT FINDINGS

In the mature gonadotroph, there is an emerging concept that differential synthesis of gonadotropin beta-subunit genes, essential for cyclic reproductive function, is associated with modification of activation and/or stability of important regulatory proteins and transcription factors. Recent studies suggest that cellular events, which affect histone modification, play an essential role in both gonadotroph development and the ontogeny of gonadotropin subunit gene expression. Such dynamic events are under the orchestration of the hypothalamic neuropeptide gonadotropin-releasing hormone (GnRH), potentially through the ability of GnRH to activate several distinct signaling cascades within the gonadotroph.

SUMMARY

Greater insight into the cellular events that are key to gonadotroph physiology will contribute to our understanding of abnormal gonadotropin secretion in disorders such as hypothalamic amenorrhea and polycystic ovarian syndrome (PCOS), and provide a context for the design of novel therapeutic approaches.

Shared kinase fluctuations between two enzymatic reactions

Kinases serve crucial roles in many cellular signaling pathways that process and transfer information. When signaling kinases phosphorylate two targets, these can serve as branch points that distribute information among two pathways. Responses to stimuli transmitted by activated kinases show high levels of cell-to-cell variation that influence cellular function. We ask how fluctuations around a steady state, due to kinase fluctuations and intrinsic noise, are distributed between two reactions with substrates phosphorylated by a shared kinase. We develop the formalism to answer this question and, for a realistic set of biological constants, we illustrate various features of fluctuations and relaxation times to a steady state. We find that the steady-state response determines the size and range in enzyme concentration of phosphorylated substrate fluctuations, and that the choice of an operating point can have a large impact on how shared kinase noise is distributed among two available pathways.

Functional characterization of genetically labeled gonadotropes

Gonadotropes are crucial in the control of reproduction but difficult to isolate for functional analysis due to their scattered distribution in the anterior pituitary gland. We devised a binary genetic approach, and describe a new mouse model that allows visualization and manipulation of gonadotrope cells. Using gene targeting in embryonic stem cells, we generated mice in which Cre recombinase is coexpressed with the GnRH receptor, which is expressed in gonadotrope cells. We show that we can direct Cre-mediated recombination of a yellow fluorescent protein reporter allele specifically in gonadotropes within the anterior pituitary of these knock-in mice. More than 99% of gonadotropin-containing cells were labeled by yellow fluorescent protein fluorescence and readily identifiable in dissociated pituitary cell culture, allowing potentially unbiased sampling from the gonadotrope population. Using electrophysiology, calcium imaging, and the study of secretion on the single-cell level, the functional properties of gonadotropes isolated from male mice were analyzed. Our studies demonstrate a significant heterogeneity in the resting properties of gonadotropes and their responses to GnRH. About 50% of gonadotropes do not exhibit secretion of LH or FSH. Application of GnRH induced a broad range of both electrophysiological responses and increases in the intracellular calcium concentration. Our mouse model will also be able to direct expression of other Cre recombination-dependent reporter genes to gonadotropes and, therefore, represents a versatile new tool in the understanding of gonadotrope biology.

A proteomic comparison of immature and mature mouse gonadotrophs reveals novel differentially expressed nuclear proteins that regulate gonadotropin gene transcription and RNA splicing

The alphaT3-1 and LbetaT2 gonadotroph cell lines contain all the known factors required for expression of gonadotropin genes, yet only the LbetaT2 cells express the beta subunits. We hypothesized that comparison of their nuclear proteomes would reveal novel proteins and/or modifications that regulate expression of these genes. We identified nine proteins with different expression profiles in the two cell lines, of which several were chosen for further functional studies. Of those found at higher levels in alphaT3-1 nuclei, 1110005A23RIK was found associated with the Fshb gene promoter and repressed its expression. Transgelin 3 overexpression reduced transcript levels of Fshb, and its knockdown elevated Lhb and Cga transcript levels, indicating an ongoing repressive effect on these more highly expressed genes, possibly through altering levels of phosphorylated mitogen-activated protein kinase. Heterogeneous nuclear ribonucleoprotein A2/B1 repressed splicing of the Fshb primary transcript, which it binds in the first intron. Proteins at higher levels in LbetaT2 nuclei included prohibitin, the overexpression of which reduced promoter activity of all three gonadotropin subunits, and appeared to mediate the differential effect of GnRH on proliferation of the two cell lines; its knockdown also altered cell morphology. Two other splicing factors were also found at higher levels in LbetaT2 nuclei: the knockdown of PRPF19 or EIF4A3 decreased splicing of Lhb, or of both beta subunit transcripts, respectively. The levels of Eif4a3 mRNA were increased by activin, and both factors increased Fshb splicing. This study has revealed a number of novel factors that alter gonadotropin expression and gonadotroph function, and likely mediate or moderate effects of the regulatory hormones.

Noise propagation and scaling in regulation of gonadotrope biosynthesis

Reproductive physiology depends on the control of biosynthesis in the pituitary gonadotrope by hypothalamic gonadotropin-releasing hormone (GnRH). The responses to GnRH include activation of extracellular signal-regulated kinase (ERK) and induction of Egr1. Using population and single cell signaling assays, we investigated the signal and noise transmission through this signaling and gene circuit, analyzing data obtained from 43,775 individual cells in 40 experiments. After exposure to GnRH, phosphorylated ERK (pERK) is elevated in 50% of the cells at 1.7 (SD = 0.3) min. Studies of the cell-to-cell response showed that for both pERK and for Egr1 protein production the mean response (mu) and standard deviation (sigma) within individual cells were linearly related (sigma = kmu) and had similar values of k. To understand the basis for the scaling observed for noise propagation through this system, we determined the relationship between pERK and egr1 mRNA levels induced at varying concentration of GnRH. While both pERK and egr1 mRNA show a saturating sigmoidal relationship to the concentration of GnRH exposure, egr1 mRNA is linearly related to the levels of pERK. These results explain the basis for variation in cellular responses in an important mammalian signaling pathway leading to gene induction.

Competing docking interactions can bring about bistability in the MAPK cascade

Mitogen-activated protein kinases are crucial regulators of various cell fate decisions including proliferation, differentiation, and apoptosis. Depending on the cellular context, the Raf-Mek-Erk mitogen-activated protein kinase cascade responds to extracellular stimuli in an all-or-none manner, most likely due to bistable behavior. Here, we describe a previously unrecognized positive-feedback mechanism that emerges from experimentally observed sequestration effects in the core Raf-Mek-Erk cascade. Unphosphorylated/monophosphorylated Erk sequesters Mek into Raf-inaccessible complexes upon weak stimulation, and thereby inhibits cascade activation. Mek, once phosphorylated by Raf, triggers Erk phosphorylation, which in turn induces dissociation of Raf-inaccessible Mek-Erk heterodimers, and thus further amplifies Mek phosphorylation. We show that this positive circuit can bring about bistability for parameter values measured experimentally in living cells. Previous studies revealed that bistability can also arise from enzyme depletion effects in the Erk double (de)phosphorylation cycle. We demonstrate that the feedback mechanism proposed in this article synergizes with such enzyme depletion effects to bring about a much larger bistable range than either mechanism alone. Our results show that stable docking interactions and competition effects, which are common in protein kinase cascades, can result in sequestration-based feedback, and thus can have profound effects on the qualitative behavior of signaling pathways.

Pulse sensitivity of the luteinizing hormone beta promoter is determined by a negative feedback loop Involving early growth response-1 and Ngfi-A binding protein 1 and 2

The hypothalamic-pituitary-gonadal endocrine axis regulates reproduction through estrous phase-dependent release of the heterodimeric gonadotropic glycoprotein hormones, LH and FSH, from the gonadotropes of the anterior pituitary. Gonadotropin synthesis and release is dependent upon pulsatile stimulation by the hypothalamic neuropeptide GnRH. Alterations in pulse frequency and amplitude alter the relative levels of gonadotropin synthesis and release. The mechanism of interpretation of GnRH pulse frequency and amplitude by gonadotropes is not understood. We have examined gene expression in LbetaT2 gonadotropes under various pulse regimes in a cell perifusion system by microarray and identified 1127 genes activated by tonic or pulsatile GnRH. Distinct patterns of expression are associated with each pulse frequency, but the greatest changes occur at a 60-min or less interpulse interval. The immediate early gene mRNAs encoding early growth response (Egr)1 and Egr2, which activate the gonadotropin LH beta-subunit gene promoter, are stably induced at high pulse frequency. In contrast, mRNAs for the Egr corepressor genes Ngfi-A binding protein Nab1 and Nab2 are stably induced at low pulse frequency. We show that Ngfi-A binding protein members inhibit Egr-mediated frequency-dependent induction of the LH beta-subunit promoter. This pattern of expression suggests a model of pulse frequency detection that acts by suppressing activation by Egr family members at low frequency and allowing activation at sustained high-frequency pulses.

Subpopulations of somatotropes with differing intracellular calcium concentration responses to secretagogues

Multiple secretagogues stimulate the release of growth hormone (GH). The present studies examined the ability of chicken somatotropes to respond to GH secretagogues with increased intracellular calcium concentrations ([Ca2+]i). It was hypothesized that there are subsets of the somatotrope population with different responsiveness to the various secretagogues. Somatotropes were identified and distinguished from other adenohypophyseal cells, by their unique ability to respond to GH-releasing hormone with increased [Ca2+]i with immunocytochemistry used as a post-hoc confirmatory test. Large increases in [Ca2+]i (222 +/- 16 nM) were evoked by thyrotropin-releasing hormone in only 73% of the somatotropes. Similarly, [Ca2+]i was increased by perifusion with pituitary adenylate cyclase-activating peptide in 85% and by leptin but only in 51% of somatotropes. Ghrelin acutely increased [Ca2+]i in only 21% of somatotropes. Perfusion with gonadotropin-releasing hormone elevated [Ca2+]i, but in only 40% of somatotropes. The kinetics of calcium transients and the magnitude of the response differed from those observed in the presumptive gonadotropes. It is concluded that there are subsets of the somatotrope population in the anterior pituitary gland with differences in their ability to respond to various secretagogues.