Correlated receptor transport processes buffer single-cell heterogeneity

The ability to sense the environment is heterogeneously distributed in cell populations

Channel capacity of signaling networks quantifies their fidelity in sensing extracellular inputs. Low estimates of channel capacities for several mammalian signaling networks suggest that cells can barely detect the presence/absence of environmental signals. However, given the extensive heterogeneity and temporal stability of cell state variables, we hypothesize that the sensing ability itself may depend on the state of the cells. In this work, we present an information-theoretic framework to quantify the distribution of sensing abilities from single-cell data. Using data on two mammalian pathways, we show that sensing abilities are widely distributed in the population and most cells achieve better resolution of inputs compared to an 'average cell'. We verify these predictions using live-cell imaging data on the IGFR/FoxO pathway. Importantly, we identify cell state variables that correlate with cells' sensing abilities. This information-theoretic framework will significantly improve our understanding of how cells sense in their environment.

The ability to sense the environment is heterogeneously distributed in cell populations

Channel capacity of signaling networks quantifies their fidelity in sensing extracellular inputs. Low estimates of channel capacities for several mammalian signaling networks suggest that cells can barely detect the presence/absence of environmental signals. However, given the extensive heterogeneity and temporal stability of cell state variables, we hypothesize that the sensing ability itself may depend on the state of the cells. In this work, we present an information theoretic framework to quantify the distribution of sensing abilities from single cell data. Using data on two mammalian pathways, we show that sensing abilities are widely distributed in the population and most cells achieve better resolution of inputs compared to an " average cell ". We verify these predictions using live cell imaging data on the IGFR/FoxO pathway. Importantly, we identify cell state variables that correlate with cells' sensing abilities. This information theoretic framework will significantly improve our understanding of how cells sense in their environment.

Modeling Cellular Signaling Variability Based on Single-Cell Data: The TGFβ-SMAD Signaling Pathway

Nongenetic heterogeneity is key to cellular decisions, as even genetically identical cells respond in very different ways to the same external stimulus, e.g., during cell differentiation or therapeutic treatment of disease. Strong heterogeneity is typically already observed at the level of signaling pathways that are the first sensors of external inputs and transmit information to the nucleus where decisions are made. Since heterogeneity arises from random fluctuations of cellular components, mathematical models are required to fully describe the phenomenon and to understand the dynamics of heterogeneous cell populations. Here, we review the experimental and theoretical literature on cellular signaling heterogeneity, with special focus on the TGFβ/SMAD signaling pathway.

Data-based stochastic modeling reveals sources of activity bursts in single-cell TGF-β signaling

Cells sense their surrounding by employing intracellular signaling pathways that transmit hormonal signals from the cell membrane to the nucleus. TGF-β/SMAD signaling encodes various cell fates, controls tissue homeostasis and is deregulated in diseases such as cancer. The pathway shows strong heterogeneity at the single-cell level, but quantitative insights into mechanisms underlying fluctuations at various time scales are still missing, partly due to inefficiency in the calibration of stochastic models that mechanistically describe signaling processes. In this work we analyze single-cell TGF-β/SMAD signaling and show that it exhibits temporal stochastic bursts which are dose-dependent and whose number and magnitude correlate with cell migration. We propose a stochastic modeling approach to mechanistically describe these pathway fluctuations with high computational efficiency. Employing high-order numerical integration and fitting to burst statistics we enable efficient quantitative parameter estimation and discriminate models that assume noise in different reactions at the receptor level. This modeling approach suggests that stochasticity in the internalization of TGF-β receptors into endosomes plays a key role in the observed temporal bursting. Further, the model predicts the single-cell dynamics of TGF-β/SMAD signaling in untested conditions, e.g., successfully reflects memory effects of signaling noise and cellular sensitivity towards repeated stimulation. Taken together, our computational framework based on burst analysis, noise modeling and path computation scheme is a suitable tool for the data-based modeling of complex signaling pathways, capable of identifying the source of temporal noise.

Second-generation non-hematopoietic erythropoietin-derived peptide for neuroprotection

Erythropoietin (EPO) is a well-known erythropoietic cytokine having a tissue-protective effect in various tissues against hypoxic stress, including the brain. Thus, its recombinants may function as neuroprotective compounds. However, despite considerable neuroprotective effects, the EPO-based therapeutic approach has side effects, including hyper-erythropoietic and tumorigenic effects. Therefore, some modified forms and derivatives of EPO have been proposed to minimize the side effects. In this study, we generated divergently modified new peptide analogs derived from helix C of EPO, with several amino acid replacements that interact with erythropoietin receptors (EPORs). This modification resulted in unique binding potency to EPOR. Unlike recombinant EPO, among the peptides, ML1-h3 exhibited a potent neuroprotective effect against oxidative stress without additional induction of cell-proliferation, owing to a differential activating mode of EPOR signaling. Furthermore, it inhibited neuronal death and brain injury under hypoxic stress in vitro and in an in vivo ischemic brain injury model. Therefore, the divergent modification of EPO-derivatives for affinity to EPOR could provide a basis for a more advanced and optimal neuroprotective strategy.

•

Short peptides derived from helix C of EPO have a neuroprotective effect.

•

Divergent modification of EPO-derived peptides has a differential affinity to EPOR.

•

ML1 and its analogs have differential cell protective and proliferative effects.

•

ML1-h3 protects neurons by suppressing in vitro oxidative stress.

•

ML1-h3 mitigates brain injury in the in vivo mouse ischemic model without hematopoietic effect.

Identification of key genes in the tumor microenvironment of lung adenocarcinoma

The tumor microenvironment plays an important role in tumor development and progression, but the role of immune and stromal cells in this environment has not been sufficiently studied. In this study, we aimed to identify key genes associated with the microenvironment of lung adenocarcinoma (LUAD). Raw data for stromal and immune cells in malignant tumors were downloaded from The Cancer Genome Atlas (TCGA). These expression data were used to identify the differentially expressed genes (DEGs) in tissue samples of LUAD with high and low immune scores. A protein-protein interaction (PPI) network based on genes with significant differential expression was constructed. Additionally, Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway analyses were performed to functionally annotate putative hub genes. These genes were assessed via Kaplan Meier analysis to determine their correlation with overall survival. In total, we identified 216 DEGs which were correlated with immune and stromal scores, including 30 hub genes which were identified based on the PPI network. Further analysis suggested that the expression levels of 10 of these genes were significantly correlated with overall survival of LUAD patients. These key hub genes included CCR2, CCR5, CD53, CYBB, HCK, IRF8, LCP2, PLEK, PTPRC, and TLR7. Moreover, the expression level of CCR2 was found to have strong prognostic value for LUAD patients. Additionally, high expression of CYBB was also correlated with better survival of patients with LUAD. The results of this study open several new avenues to explore in the treatment of LUAD.

Maximum Entropy Framework for Predictive Inference of Cell Population Heterogeneity and Responses in Signaling Networks

Predictive models of signaling networks are essential for understanding cell population heterogeneity and designing rational interventions in disease. However, using computational models to predict heterogeneity of signaling dynamics is often challenging because of the extensive variability of biochemical parameters across cell populations. Here, we describe a maximum entropy-based framework for inference of heterogeneity in dynamics of signaling networks (MERIDIAN). MERIDIAN estimates the joint probability distribution over signaling network parameters that is consistent with experimentally measured cell-to-cell variability of biochemical species. We apply the developed approach to investigate the response heterogeneity in the EGFR/Akt signaling network. Our analysis demonstrates that a significant fraction of cells exhibits high phosphorylated Akt (pAkt) levels hours after EGF stimulation. Our findings also suggest that cells with high EGFR levels predominantly contribute to the subpopulation of cells with high pAkt activity. We also discuss how MERIDIAN can be extended to accommodate various experimental measurements.