Autocrine loops with positive feedback enable context-dependent cell signaling

We describe a mechanism for context-dependent cell signaling mediated by autocrine loops with positive feedback. We demonstrate that the composition of the extracellular medium can critically influence the intracellular signaling dynamics induced by extracellular stimuli. Specifically, in the epidermal growth factor receptor (EGFR) system, amplitude and duration of mitogen-activated protein kinase (MAPK) activation are modulated by the positive-feedback loop formed by the EGFR, the Ras-MAPK signaling pathway, and a ligand-releasing protease. The signaling response to a transient input is short-lived when most of the released ligand is lost to the cellular microenvironment by diffusion and/or interaction with an extracellular ligand-binding component. In contrast, the response is prolonged or persistent in a cell that is efficient in recapturing the endogenous ligand. To study functional capabilities of autocrine loops, we have developed a mathematical model that accounts for ligand release, transport, binding, and intracellular signaling. We find that context-dependent signaling arises as a result of dynamic interaction between the parts of an autocrine loop. Using the model, we can directly interpret experimental observations on context-dependent responses of autocrine cells to ionizing radiation. In human carcinoma cells, MAPK signaling patterns induced by a short pulse of ionizing radiation can be transient or sustained, depending on cell type and composition of the extracellular medium. On the basis of our model, we propose that autocrine loops in this, and potentially other, growth factor and cytokine systems may serve as modules for context-dependent cell signaling.     

Stochastic model of autocrine and paracrine signals in cell culture assays

Autocrine signaling systems are commonly studied under cell culture conditions. In a typical cell culture assay, a layer of liquid medium covers a random two-dimensional dispersion of cells, which secrete ligands. In a growing number of experiments, it is important to characterize the spatial range of autocrine and paracrine cell communication. Currently, the spatial distribution of diffusing signals can be analyzed only indirectly, from their effects on the intracellular signaling or physiological responses of autocrine cells. To directly characterize the spatial range of secreted ligands, we propose a stochastic model for autocrine cell cultures and analyze it using a combination of analytical and computational tools. The two main results derived within the framework of this model are 1), an expression for the fraction of autocrine trajectories, i.e., the probability for a ligand to be trapped by the same cell from which it has been secreted; and 2), an expression for the spatial distribution of trapping points of paracrine trajectories. We test these analytical results by stochastic simulations with efficient Brownian dynamics code and apply our model to analyze the spatial operation of autocrine epidermal growth factor receptor systems.     

Spatial range of autocrine signaling: modeling and computational analysis

Autocrine loops formed by growth factors and their receptors have been identified in a large number of developmental, physiological, and pathological contexts. In general, the spatially distributed and recursive nature of autocrine signaling systems makes their experimental analysis, and often even their detection, very difficult. Here, we combine Brownian motion theory, Monte Carlo simulations, and reaction-diffusion models to analyze the spatial operation of autocrine loops. Within this modeling framework, the ability of autocrine cells to recapture the endogenous ligand and the distances traveled by autocrine ligands are explicitly related to ligand diffusion coefficients, density of surface receptors, ligand secretion rate, and rate constants of ligand binding and endocytic internalization. Applying our models to study autocrine loops in the epidermal growth factor receptor system, we find that autocrine loops can be highly localized--even at the level of a single cell. We demonstrate how the variations in molecular and cellular parameters may "tune" the spatial range of autocrine signals over several orders of magnitude: from microns to millimeters. We argue that this versatile regulation of the spatial range of autocrine signaling enables autocrine cells to perceive a broad spectrum of environmental information.     

Suppression of α-catenin and adherens junctions enhances epithelial cell proliferation and motility via TACE-mediated TGF-α autocrine/paracrine signaling

Sustained cell migration is essential for wound healing and cancer metastasis. The epidermal growth factor receptor (EGFR) signaling cascade is known to drive cell migration and proliferation. While the signal transduction downstream of EGFR has been extensively investigated, our knowledge of the initiation and maintenance of EGFR signaling during cell migration remains limited. The metalloprotease TACE (tumor necrosis factor alpha converting enzyme) is responsible for producing active EGFR family ligands in the via ligand shedding. Sustained TACE activity may perpetuate EGFR signaling and reduce a cell’s reliance on exogenous growth factors. Using a cultured keratinocyte model system, we show that depletion of α-catenin perturbs adherens junctions, enhances cell proliferation and motility, and decreases dependence on exogenous growth factors. We show that the underlying mechanism for these observed phenotypical changes depends on enhanced autocrine/paracrine release of the EGFR ligand transforming growth factor alpha in a TACE-dependent manner. We demonstrate that proliferating keratinocyte epithelial cell clusters display waves of oscillatory extracellular signal–regulated kinase (ERK) activity, which can be eliminated by TACE knockout, suggesting that these waves of oscillatory ERK activity depend on autocrine/paracrine signals produced by TACE. These results provide new insights into the regulatory role of adherens junctions in initiating and maintaining autocrine/paracrine signaling with relevance to wound healing and cellular transformation.     

An integrated agent-mathematical model of the effect of intercellular signalling via the epidermal growth factor receptor on cell proliferation

We have previously developed Epitheliome, a software agent representation of the growth and repair characteristics of epithelial cell populations, where cell behaviour is governed by a number of simple rules. In this paper, we describe how this model has been extended to incorporate an example of a molecular 'mechanism' behind a rule-in this case, how signalling by both endogenous and exogenous ligands of the epidermal growth factor receptor (EGFR) can impact on the proliferation of cell agents. We have developed a mathematical model representing release of endogenous ligand by cells, three-dimensional diffusion of the secreted molecules through a volume of cell culture medium, ligand-receptor binding, and bound receptor internalization and trafficking. Information relating to quantities of molecular species associated with each cell agent is frequently exchanged between the agent and signalling models, and the ratio of bound to free receptors determines cell cycle progression and hence the proliferative behaviour of the cell agents. We have applied this integrated model to examine the effect of plating density on tissue growth via autocrine/paracrine signalling. This predicts that cell growth is dependent on the concentration of exogenous ligand, but where this is limited, then growth becomes dependent on cell density and the availability of endogenous ligand. We have further modified the calcium concentration of the medium to modulate the formation of intercellular bonds between cells and shown that the increased propensity for cells to form colonies in physiological calcium does not result in significantly different patterns of receptor occupancy. In conclusion, our approach demonstrates that by combining agent-based and mathematical modelling paradigms, it is possible to probe the complex feedback relationship between the behaviour of individual cells and their interaction with one another and their environment.     

Mathematical modeling of Erk activity waves in regenerating zebrafish scales

Erk signaling regulates cellular decisions in many biological contexts. Recently, we have reported a series of Erk activity traveling waves that coordinate regeneration of osteoblast tissue in zebrafish scales. These waves originate from a central source region, propagate as expanding rings, and impart cell growth, thus controlling tissue morphogenesis. Here, we present a minimal reaction-diffusion model for Erk activity waves. The model considers three components: Erk, a diffusible Erk activator, and an Erk inhibitor. Erk stimulates both its activator and inhibitor, forming a positive and negative feedback loop, respectively. Our model shows that this system can be excitable and propagate Erk activity waves. Waves originate from a pulsatile source that is modeled by adding a localized basal production of the activator, which turns the source region from an excitable to an oscillatory state. As Erk activity periodically rises in the source, it can trigger an excitable wave that travels across the entire tissue. Analysis of the model finds that positive feedback controls the properties of the traveling wavefront and that negative feedback controls the duration of Erk activity peak and the period of Erk activity waves. The geometrical properties of the waves facilitate constraints on the effective diffusivity of the activator, indicating that waves are an efficient mechanism to transfer growth factor signaling rapidly across a large tissue.     

Affinity regulates spatial range of EGF receptor autocrine ligand binding

Proper spatial localization of EGFR signaling activated by autocrine ligands represents a critical factor in embryonic development as well as tissue organization and function, and ligand/receptor binding affinity is among the molecular and cellular properties suggested to play a role in governing this localization. We employ a computational model to predict how receptor-binding affinity affects local capture of autocrine ligand vis-a-vis escape to distal regions, and provide experimental test by constructing cell lines expressing EGFR along with either wild-type EGF or a low-affinity mutant, EGF(L47M). The model predicts local capture of a lower affinity autocrine ligand to be less efficient when the ligand production rate is small relative to receptor appearance rate. Our experimental data confirm this prediction, demonstrating that cells can use ligand/receptor binding affinity to regulate ligand spatial distribution when autocrine ligand production is limiting for receptor signaling.     

Interrupting autocrine ligand-receptor binding: comparison between receptor blockers and ligand decoys.

Stimulation of cell behavioral functions by ligand/receptor binding can be accomplished in autocrine fashion, where cells secrete ligand capable of binding to receptors on their own surfaces. This proximal secretion of autocrine ligands near the surface receptors on the secreting cell suggests that control of these systems by inhibitors of receptor/ligand binding may be more difficult than for systems involving exogenous ligands. Hence, it is of interest to predict the conditions under which successful inhibition of cell receptor binding by the autocrine ligand can be expected. Previous theoretical work using a compartmentalized model for autocrine cells has elucidated the conditions under which addition of solution decoys for the autocrine ligand can interrupt cell receptor/ligand binding via competitive binding of the secreted molecules (Forsten, K. E., and D. A. Lauffenburger. 1992. Biophys. J. 61:1-12.) We now apply a similar modeling approach to examine the addition of solution blockers targeted against the cell receptor. Comparison of the two alternative inhibition strategies reveals that a significantly lower concentration of receptor blockers, compared to ligand decoys, will obtain a high degree of inhibition. The more direct interruption scheme characteristic of the receptor blockers may make them a preferred strategy when feasible.     

Time and length scales of autocrine signals in three dimensions

A model of autocrine signaling in cultures of suspended cells is developed on the basis of the effective medium approximation. The fraction of autocrine ligands, the mean and distribution of distances traveled by paracrine ligands before binding, as well as the mean and distribution of the ligand lifetime are derived. Interferon signaling by dendritic immune cells is considered as an illustration.     

Discrete models of autocrine cell communication in epithelial layers

Pattern formation in epithelial layers heavily relies on cell communication by secreted ligands. Whereas the experimentally observed signaling patterns can be visualized at single-cell resolution, a biophysical framework for their interpretation is currently lacking. To this end, we develop a family of discrete models of cell communication in epithelial layers. The models are based on the introduction of cell-to-cell coupling coefficients that characterize the spatial range of intercellular signaling by diffusing ligands. We derive the coupling coefficients as functions of geometric, cellular, and molecular parameters of the ligand transport problem. Using these coupling coefficients, we analyze a nonlinear model of positive feedback between ligand release and binding. In particular, we study criteria of existence of the patterns consisting of clusters of a few signaling cells, as well as the onset of signal propagation. We use our model to interpret recent experimental studies of the EGFR/Rhomboid/Spitz module in Drosophila development.     

A model for the behaviour of phosphorylase b. The generation of different binding sites via intermediate enzymatic states.

The model given in this paper can be applied to enzymatic systems which have more than two conformational states in equilibrium and which clearly exhibit heterogeneity in the binding of one ligand. The model we propose makes possible quantitative interpretation of our experimental results and of those of many other workers as well. In some cases calorimetric, dialysis and kinetic magnitudes, when plotted against ligand concentration, give multiregional or "stepwise" curves. We suggest that such a behaviour arises because total occupation of one class of binding sites completely moves the enzyme towards a different conformational state in which the affinity for the ligand is greatly increased by the formation of a new class of binding sites. Our calorimetric results for the interaction between some nucleotides and phosphorylase b closely conform to our model.     

Autocrine ligand binding to cell receptors. Mathematical analysis of competition by solution "decoys".

Autocrine ligands have been demonstrated to regulate cell proliferation, cell adhesion, and cell migration in a number of different systems and are believed to be one of the underlying causes of malignant cell transformation. Binding of these ligands to their cellular receptors can be compromised by diffusive transport of ligand away from the secreting cell. Exogenous addition of antibodies or solution receptors capable of competing with cellular receptors for these autocrine ligands has been proposed as a means of inhibiting autocrine-stimulated cell behavioral responses. Such "decoys" complicate cellular binding by offering alternative binding targets, which may also be capable of aiding or abating transport of the ligand away from the cell surface. We present a mathematical model incorporating autocrine ligand production and the presence of competing cellular and solution receptors. We elucidate effects of key system parameters including ligand diffusion rate, binding rate constants, cell density, and secretion rate on the ability of solution receptors to inhibit cellular receptor binding. Both plated and suspension cell systems are considered. An approximate analytical expression relating the key parameters to the critical concentration of solution "decoys" required for inhibition is derived and compared to the numerical calculations. We find that in order to achieve essentially complete inhibition of surface receptor binding, the concentration of decoys may need to be as much as four to eight orders of magnitude greater than the equilibrium disociation constant for ligand binding to surface receptors.     

A novel proteomic coculture model of prostate cancer cell growth

Chemotherapy and androgen ablation therapy are only temporarily effective against prostate cancer, and current studies are ongoing to test agents that target proteins responsible for autocrine and paracrine stimulated growth. Given limitations of current laboratory models to test the effect of these agents on cell growth and protein targets, we developed a coculture model that can distinguish paracrine stimulated growth and effects on proteins. We found that LNCaP prostate cancer cells and an immortalized rat prostate cell line transfected to overexpress the antiapoptotic resistance protein Bcl-2 were stimulated to grow (>2-fold increase, p < 0.01) through autocrine effects from additional cells in an upper chamber of our system. Using a proteomic approach with a two-dimensional differential in gel electrophoresis method to increase fidelity, four proteins were found to increase after autocrine induced growth stimulation. These proteins were all identified by mass spectrometry as enzymes in the glycolytic pathway, validating the ability of this system to detect both clonogenic growth and the effect on proteins. These data, therefore, demonstrate a novel coculture model for further study of agents that target proteins in pathways of paracrine or autocrine stimulated cell growth.     

ApoE-mediated cholesterol efflux from macrophages: separation of autocrine and paracrine effects

Macrophages in the vessel wall secrete high levels of apolipoprotein E (apoE). Cholesterol efflux from macrophages to apoE has been shown to decrease foam cell formation and prevent atherosclerosis. An apoE molecule can mediate cholesterol efflux from the macrophage that originally secreted it (autocrine effect) or from surrounding macrophages (paracrine effect). Traditional methodologies have not been able to separate these serial effects. The novel methodology presented here was developed to separate autocrine and paracrine effects by using a simple mathematical model to interpret the effects of dilution on apoE-mediated cholesterol efflux. Our results show that, at very dilute concentrations, the paracrine effect of apoE is not evident and the autocrine effect becomes the dominant mediator of efflux. However, at saturating concentrations, paracrine apoE causes 80–90% of the apoE-mediated cholesterol efflux, whereas autocrine apoE is responsible for the remaining 10–20%. These results suggest that the relative importance of autocrine and paracrine apoE depends on the size of the local distribution volume, a factor not considered in previous in vitro studies of apoE function. Furthermore, autocrine effects of apoE could be critical in the prevention of foam cell formation in vivo. This novel methodology may be applicable to other types of mixed autocrine/paracrine systems, such as signal transduction systems. autocrine/paracrine system; cholesterol acceptor; extracellular space; distribution volume     

Ligand accumulation in autocrine cell cultures

Cell-culture assays are routinely used to analyze autocrine signaling systems, but quantitative experiments are rarely possible. To enable the quantitative design and analysis of experiments with autocrine cells, we develop a biophysical theory of ligand accumulation in cell-culture assays. Our theory predicts the ligand concentration as a function of time and measurable parameters of autocrine cells and cell-culture experiments. The key step of our analysis is the derivation of the survival probability of a single ligand released from the surface of an autocrine cell. An expression for this probability is derived using the boundary homogenization approach and tested by stochastic simulations. We use this expression in the integral balance equations, from which we find the Laplace transform of the ligand concentration. We demonstrate how the theory works by analyzing the autocrine epidermal growth factor receptor system and discuss the extension of our methods to other experiments with cultured autocrine cells.     

Microfluidic perfusion for regulating diffusible signaling in stem cells

Background

Autocrine & paracrine signaling are widespread both in vivo and in vitro, and are particularly important in embryonic stem cell (ESC) pluripotency and lineage commitment. Although autocrine signaling via fibroblast growth factor-4 (FGF4) is known to be required in mouse ESC (mESC) neuroectodermal specification, the question of whether FGF4 autocrine signaling is sufficient, or whether other soluble ligands are also involved in fate specification, is unknown. The spatially confined and closed-loop nature of diffusible signaling makes its experimental control challenging; current experimental approaches typically require prior knowledge of the factor/receptor in order to modulate the loop. A new approach explored in this work is to leverage transport phenomena at cellular resolution to downregulate overall diffusible signaling through the physical removal of cell-secreted ligands.

Methodology/Principal Findings

We develop a multiplex microfluidic platform to continuously remove cell-secreted (autocrine\paracrine) factors to downregulate diffusible signaling. By comparing cell growth and differentiation in side-by-side chambers with or without added cell-secreted factors, we isolate the effects of diffusible signaling from artifacts such as shear, nutrient depletion, and microsystem effects, and find that cell-secreted growth factor(s) are required during neuroectodermal specification. Then we induce FGF4 signaling in minimal chemically defined medium (N2B27) and inhibit FGF signaling in fully supplemented differentiation medium with cell-secreted factors to determine that the non-FGF cell-secreted factors are required to promote growth of differentiating mESCs.

Conclusions/Significance

Our results demonstrate for the first time that flow can downregulate autocrine\paracrine signaling and examine sufficiency of extracellular factors. We show that autocrine\paracrine signaling drives neuroectodermal commitment of mESCs through both FGF4-dependent and -independent pathways. Overall, by uncovering autocrine\paracrine processes previously hidden in conventional culture systems, our results establish microfluidic perfusion as a technique to study and manipulate diffusible signaling in cell systems.     

Early embryonic vascular patterning by matrix-mediated paracrine signalling: a mathematical model study

During embryonic vasculogenesis, endothelial precursor cells of mesodermal origin known as angioblasts assemble into a characteristic network pattern. Although a considerable amount of markers and signals involved in this process have been identified, the mechanisms underlying the coalescence of angioblasts into this reticular pattern remain unclear. Various recent studies hypothesize that autocrine regulation of the chemoattractant vascular endothelial growth factor (VEGF) is responsible for the formation of vascular networks in vitro. However, the autocrine regulation hypothesis does not fit well with reported data on in vivo early vascular development. In this study, we propose a mathematical model based on the alternative assumption that endodermal VEGF signalling activity, having a paracrine effect on adjacent angioblasts, is mediated by its binding to the extracellular matrix (ECM). Detailed morphometric analysis of simulated networks and images obtained from in vivo quail embryos reveals the model mimics the vascular patterns with high accuracy. These results show that paracrine signalling can result in the formation of fine-grained cellular networks when mediated by angioblast-produced ECM. This lends additional support to the theory that patterning during early vascular development in the vertebrate embryo is regulated by paracrine signalling.     

TLR agonists regulate PDGF-B production and cell proliferation through TGF-beta/type I IFN crosstalk

Transforming growth factor-beta (TGF-beta) and type I interferon (IFN) autocrine/paracrine loops are recognized as key mediators of signaling cascades that control a variety of cellular functions. Here, we describe a novel mechanism by which Toll-like receptor (TLR) agonists utilize these two autocrine/paracrine loops to differentially regulate the induction of PDGF-B, a growth factor implicated in a number of diseases ranging from tumor metastasis to glomerulonephritis. We demonstrate that CpG-specific induction of PDGF-B requires activation of Smads through TGFbeta1 autocrine/paracrine signaling. In contrast, polyinosinic:polycytidylic acid strongly represses CpG's as well as its own intrinsic ability to induce PDGF-B mRNA through type I IFN-mediated induction of Smad7, a negative regulator of Smad3/4. Furthermore, we have shown that this crosstalk mechanism translates into similar regulation of mesangial cell proliferation. Thus, our results demonstrate the importance of crosstalk between TGF-beta and type I IFNs in determining the specificity of TLR-mediated gene induction.     

Glucocorticoids regulate the secretion of a 21kDa-IGF-binding protein by sheep adipose tissue explants

Using a solution phase assay we have demonstrated that sheep adipose tissue explants secrete insulin-like growth factor binding proteins (IGFBPs) when cultured in serum-free medium over a 24 h period. Further, we demonstrate that secretion of IGFBP(s) is inhibited (up to 50%) by incubation of the cultures in the presence of 10^–8M dexamethasone. This inhibitory effects is overcome when insulin (10 ng/ml) and ovine growth hormone (100 ng/ml) are incubated together (but not separately) with glucocorticoid. Further characterisation of this IGF binding activity by high performance size exclusion chromatography and Western ligand blot analysis indicated that under our culture conditions sheep adipose tissue explants secrete one predominant 21 kDa IGFBP and it is this BP which is hormonally regulated as described above. We discuss our results in the context of endocrine/paracrine/autocrine control of adipose tissue metabolism and differentiation.Abbreviations IGF insulin-like growth factor - IGF-BP insulin-like growth factor binding proteins - DX dexamethasone - GH ovine growth hormone - CM conditioned medium