Long-Range Signal Transmission in Autocrine Relays

Intracellular signaling induced by peptide growth factors can stimulate secretion of these molecules into the extracellular medium. In autocrine and paracrine networks, this can establish a positive feedback loop between ligand binding and ligand release. When coupled to intercellular communication by autocrine ligands, this positive feedback can generate constant-speed traveling waves. To demonstrate that, we propose a mechanistic model of autocrine relay systems. The model is relevant to the physiology of epithelial layers and to a number of in vitro experimental formats. Using asymptotic and numerical tools, we find that traveling waves in autocrine relays exist and have a number of unusual properties, such as an optimal ligand binding strength necessary for the maximal speed of propagation. We compare our results to recent observations of autocrine and paracrine systems and discuss the steps toward experimental tests of our predictions.     

TLR and NKG2D Signaling Pathways Mediate CS-Induced Pulmonary Pathologies

Long-term exposure to cigarette smoke (CS) can have deleterious effects on lung epithelial cells including cell death and the initiation of inflammatory responses. CS-induced cell injury can elaborate cell surface signals and cellular byproducts that stimulate immune system surveillance. Our previous work has shown that the expression of ligands for the cytotoxic lymphocyte activating receptor NKG2D is enhanced in patients with COPD and that the induction of these ligands in a mouse model can replicate COPD pathologies. Here, we extend these findings to demonstrate a role for the NKG2D receptor in CS-induced pathophysiology and provide evidence linking nucleic acid-sensing endosomal toll-like receptor (TLR) signaling to COPD pathology through NKG2D activation. Specifically, we show that mice deficient in NKG2D exhibit attenuated pulmonary inflammation and airspace enlargement in a model of CS-induced emphysema. Additionally, we show that CS exposure induces the release of free nucleic acids in the bronchoalveolar lavage and that direct exposure of mouse lung epithelial cells to cigarette smoke extract similarly induces functional nucleic acids as assessed by TLR3, 7, and 9 reporter cell lines. We demonstrate that exposure of mouse lung epithelial cells to TLR ligands stimulates the surface expression of RAET1, a ligand for NKG2D, and that mice deficient in TLR3/7/9 receptor signaling do not exhibit CS-induced NK cell hyperresponsiveness and airspace enlargement. The findings indicate that CS-induced airway injury stimulates TLR signaling by endogenous nucleic acids leading to elevated NKG2D ligand expression. Activation of these pathways plays a major role in the altered NK cell function, pulmonary inflammation and remodeling related to long-term CS exposure.     

Complementary expression and repulsive signaling suggest that EphB2 and ephrin-B1 are possibly involved in epithelial boundary formation at the squamocolumnar junction in the rodent stomach

Eph receptors and ephrin ligands are cell–cell communication molecules with well-defined roles in cell adhesion, migration, and tissue boundary formation. However, their expression levels in the squamocolumnar epithelial junction region at the distal esophagus are completely unknown. We examined EphB2 and ephrin-B1 localization in the squamocolumnar epithelial junction region between the proximal and distal stomach of the rodents. Immunostaining showed complimentary expression patterns along the proximal-to-distal axis of the gastric epithelia across the junction: EphB2 expression was maximal around the epithelial junction and sharply decreased in the stratified squamous epithelium at a short distance from the junction, whereas ephrin-B1 was strongly expressed in the stratified squamous epithelium at a distance from the junction and sharply decreased toward the junction. These expression patterns suggest that EphB2/ephrin-B1 signaling occurs preferentially in the epithelia across the junction, where the receptor and ligand expression highly overlap. We also show that (1) EphB2 preferentially binds ephrin-B1, and (2) cell repulsion/lateral migration was induced in primary cultured gastric keratinocytes on ephrin-B1-Fc- and EphB2-Fc-coated surfaces. On the basis of these findings, we propose that EphB2 and ephrin-B1 are possibly involved in epithelial boundary formation at the squamocolumnar junction.     

Intra-membrane ligand diffusion and cell shape modulate juxtacrine patterning

A key problem in developmental biology is how pattern and planar polarity are transmitted in epithelial structures. Examples include Drosophila neuronal differentiation, ommatidia formation in the compound eye, and wing hair polarization. A key component for the generation of such patterns is direct cell-cell signalling by transmembrane ligands, called juxtacrine signalling. Previous models for this mode of communication have considered homogeneous distributions in the cell membrane, and the role of polarity has been largely ignored. In this paper we determine the role of inhomogeneous protein and receptor distributions in juxtacrine signalling. We explicitly include individual membrane segments, diffusive transport of proteins and receptors between these segments, and production terms with a combination of local and global responses to ligand binding. Our analysis shows that intra-membrane ligand transport is vital for the generation of long wavelength patterns. Moreover, with no ligand transport, there is no pattern formation for lateral induction, a process in which receptor activation up-regulates ligand production. Biased production of ligand also modulates patterning bifurcations and predicted wavelengths. In addition, biased ligand and receptor trafficking can lead to regular polarity across a lattice, in which each cell has the same orientation-directly analogous to patterns of hairs in the Drosophila wing. We confirm the trends in pattern wavelengths previously observed for patterns with cellular homogeneity-lateral inhibition tends to give short-range patterns, while lateral induction can give patterns with much longer wavelengths. Moreover, the original model can be recovered if intra-membrane bound receptor diffusion is included and rapid equilibriation between the sides is considered. Finally, we consider the role of irregular cell shapes and waves in such networks, including wave propagation past clones of non-signalling cells.     

Computational modeling of extracellular mechanotransduction

Mechanotransduction may occur through numerous mechanisms, including potentially through autocrine signaling in a dynamically changing extracellular space. We developed a computational model to analyze how alterations in the geometry of an epithelial lateral intercellular space (LIS) affect the concentrations of constitutively shed ligands inside and below the LIS. The model employs the finite element method to solve for the concentration of ligands based on the governing ligand diffusion-convection equations inside and outside of the LIS, and assumes idealized parallel plate geometry and an impermeable tight junction at the apical surface. Using the model, we examined the temporal relationship between geometric changes and ligand concentration, and the dependence of this relationship on system characteristics such as ligand diffusivity, shedding rate, and rate of deformation. Our results reveal how the kinetics of mechanical deformation can be translated into varying rates of ligand accumulation, a potentially important mechanism for cellular discrimination of varying rate-mechanical processes. Furthermore, our results demonstrate that rapid changes in LIS geometry can transiently increase ligand concentrations in underlying media or tissues, suggesting a mechanism for communication of mechanical state between epithelial and subepithelial cells. These results underscore both the plausibility and complexity of the proposed extracellular mechanotransduction mechanism.     

Competition in Notch Signaling with Cis Enriches Cell Fate Decisions

Notch signaling is involved in cell fate choices during the embryonic development of Metazoa. Commonly, Notch signaling arises from the binding of the Notch receptor to its ligands in adjacent cells driving cell-to-cell communication. Yet, cell-autonomous control of Notch signaling through both ligand-dependent and ligand-independent mechanisms is known to occur as well. Examples include Notch signaling arising in the absence of ligand binding, and cis-inhibition of Notch signaling by titration of the Notch receptor upon binding to its ligands within a single cell. Increasing experimental evidences support that the binding of the Notch receptor with its ligands within a cell (cis-interactions) can also trigger a cell-autonomous Notch signal (cis-signaling), whose potential effects on cell fate decisions and patterning remain poorly understood. To address this question, herein we mathematically and computationally investigate the cell states arising from the combination of cis-signaling with additional Notch signaling sources, which are either cell-autonomous or involve cell-to-cell communication. Our study shows that cis-signaling can switch from driving cis-activation to effectively perform cis-inhibition and identifies under which conditions this switch occurs. This switch relies on the competition between Notch signaling sources, which share the same receptor but differ in their signaling efficiency. We propose that the role of cis-interactions and their signaling on fine-grained patterning and cell fate decisions is dependent on whether they drive cis-inhibition or cis-activation, which could be controlled during development. Specifically, cis-inhibition and not cis-activation facilitates patterning and enriches it by modulating the ratio of cells in the high-ligand expression state, by enabling additional periodic patterns like stripes and by allowing localized patterning highly sensitive to the precursor state and cell-autonomous bistability. Our study exemplifies the complexity of regulations when multiple signaling sources share the same receptor and provides the tools for their characterization.     

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.     

Multiple TLRs activate EGFR via a signaling cascade to produce innate immune responses in airway epithelium

Toll-like receptors (TLRs) are critical for the recognition of inhaled pathogens that deposit on the airway epithelial surface. The epithelial response to pathogens includes signaling cascades that activate the EGF receptor (EGFR). We hypothesized that TLRs communicate with EGFR via epithelial signaling to produce certain innate immune responses. Airway epithelium expresses the highest levels of TLR2, TLR3, TLR5, and TLR6, and here we found that ligands for these TLRs increased IL-8 and VEGF production in normal human bronchial epithelial cells. These effects were prevented by treatment with a selective inhibitor of EGFR phosphorylation (AG-1478), a metalloprotease (MP) inhibitor, a reactive oxygen species (ROS) scavenger, and an NADPH oxidase inhibitor. In an airway epithelial cell line (NCI-H292), TNF-alpha-converting enzyme (TACE) small interfering RNA (siRNA) was used to confirm that TACE is the MP involved in TLR ligand-induced IL-8 and VEGF production. We show that transforming growth factor (TGF)-alpha is the EGFR ligand in this signaling cascade by using TGF-alpha neutralizing antibody and by showing that epithelial production of TGF-alpha occurs in response to TLR ligands. Dual oxidase 1 (Duox1) siRNA was used to confirm that Duox1 is the NADPH oxidase involved in TLR ligand-induced IL-8 and VEGF production. We conclude that multiple TLR ligands induce airway epithelial cell production of IL-8 and VEGF via a Duox1--> ROS--> TACE--> TGF-alpha--> EGFR phosphorylation pathway. These results show for the first time that multiple TLRs in airway epithelial cells produce innate immune responses by activating EGFR via an epithelial cell signaling cascade.     

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.     

A response calculus for immobilized T cell receptor ligands

To address the molecular mechanism of T cell receptor (TCR) signaling, we have formulated a model for T cell activation, termed the 2D-affinity model, in which the density of TCR on the T cell surface, the density of ligand on the presenting surface, and their corresponding two-dimensional affinity determine the level of T cell activation. When fitted to T cell responses against purified ligands immobilized on plastic surfaces, the 2D-affinity model adequately simulated changes in cellular activation as a result of varying ligand affinity and ligand density. These observations further demonstrated the importance of receptor cross-linking density in determining TCR signaling. Moreover, it was found that the functional two-dimensional affinity of TCR ligands was affected by the chemical composition of the ligand-presenting surface. This makes it possible that cell-bound TCR ligands, despite their low affinity in solution, are of optimal two-dimensional affinity thereby allowing effective TCR binding under physiological conditions, i.e. at low ligand densities in cellular interfaces.     

Juxtacrine Signaling Is Inherently Noisy

Juxtacrine signaling is an important class of signaling systems that plays a crucial role in various developmental processes ranging from coordination of differentiation between neighboring cells to guiding axon growth during neurogenesis. Such signaling systems rely on the interaction between receptors on one cell and trans-membrane ligands on the membrane of a neighboring cell. Like other signaling systems, the ability of signal-receiving cells to accurately determine the concentration of ligands, is affected by stochastic diffusion processes. However, it is not clear how restriction of ligand movement to the two-dimensional (2D) cell membrane in juxtacrine signaling affects the accuracy of ligand sensing. In this study, we use a statistical mechanics approach, to show that long integration times, from around one second to several hours, are required to reach high-sensing accuracy (better than 10%). Surprisingly, the accuracy of sensing cannot be significantly improved, neither by increasing the number of receptors above three to five receptors per contact area, nor by increasing the contact area between cells. We show that these results impose stringent constraints on the dynamics of processes relying on juxtacrine signaling systems, such as axon guidance mediated by Ephrins and developmental patterns mediated by the Notch pathway.     

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.     

A 3-D model of ligand transport in a deforming extracellular space

Cells communicate through shed or secreted ligands that traffic through the interstitium. Force-induced changes in interstitial geometry can initiate mechanotransduction responses through changes in local ligand concentrations. To gain insight into the temporal and spatial evolution of such mechanotransduction responses, we developed a 3-D computational model that couples geometric changes observed in the lateral intercellular space (LIS) of mechanically loaded airway epithelial cells to the diffusion-convection equations that govern ligand transport. By solving the 3-D fluid field under changing boundary geometries, and then coupling the fluid velocities to the ligand transport equations, we calculated the temporal changes in the 3-D ligand concentration field. Our results illustrate the steady-state heterogeneities in ligand distribution that arise from local variations in interstitial geometry, and demonstrate that highly localized changes in ligand concentration can be induced by mechanical loading, depending on both local deformations and ligand convection effects. The occurrence of inhomogeneities at steady state and in response to mechanical loading suggest that local variations in ligand concentration may have important effects on cell-to-cell variations in basal signaling state and localized mechanotransduction responses.     

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.     

Both LRP5 and LRP6 receptors are required to respond to physiological Wnt ligands in mammary epithelial cells and fibroblasts

A canonical Wnt signal maintains adult mammary ductal stem cell activity, and this signal requires the Wnt signaling reception, LRP5. However, previous data from our laboratory have shown that LRP5 and LRP6 are co-expressed in mammary basal cells and that LRP6 is active, leading us to question why LRP6 is insufficient to mediate canonical signaling in the absence of LRP5. Here, we show that at endogenous levels of LRP5 and LRP6 both receptors are required to signal in response to some Wnt ligands both in vitro (in mouse embryonic fibroblasts and mammary epithelial cells) and in vivo (in mammary outgrowths). This subgroup of canonical ligands includes Wnt1, Wnt9b, and Wnt10b; the latter two are expressed in mammary gland. In contrast, the ligand commonly used experimentally, Wnt3a, prefers LRP6 and requires just one receptor regardless of cellular context. When either LRP5 or LRP6 is overexpressed, signaling remains ligand-dependent, but the requirement for both receptors is abrogated (regardless of ligand type). We have documented an LRP5-6 heteromer using immiscible filtration assisted by surface tension (IFAST) immunoprecipitation. Together, our data imply that under physiological conditions some Wnt ligands require both receptors to be present to generate a canonical signal. We have designed a model to explain our results based on the resistance of LRP5-6 heteromers to a selective inhibitor of E1/2-binding Wnt-LRP6 interaction. These data have implications for stem cell biology and for the analysis of the oncogenicity of LRP receptors that are often overexpressed in breast tumors.     

Delta-like 1 expression promotes goblet cell differentiation in Notch-inactivated human colonic epithelial cells

Notch signaling has previously been implicated in the regulation of the cell fate of intestinal epithelial cells. However, the expression and function of Notch ligands in the human intestine remain largely unknown. In the present study, we showed that Notch ligands Delta-like 1 (Dll1) and Delta-like 4 (Dll4) are expressed in a goblet cell-specific manner in human colonic tissue. Additionally, we found that Dll1 and Dll4 expression was regulated in-parallel with Atoh1 and MUC2, which are both under the control of the Notch-Hes1 signaling pathway. Because knockdown of Dll1 expression completely abrogated the acquisition of the goblet cell phenotype in Notch-inactivated colonic epithelial cells, we postulate that Dll1 might function as a cis-acting regulatory element that induces undifferentiated cells to become goblet cells. Our results suggest a link between Dll1 expression and human goblet cell differentiation that might be mediated by a function that is distinct from its role as a Notch receptor ligand.