Systems biology. Accurate information transmission through dynamic biochemical signaling networks

Information transmission in a cell monolayer: A numerical study

Motivated by the spatiotemporal waves of MAPK/ERK activity, crucial for long-range communication in regenerating tissues, we investigated stochastic homoclinic fronts propagating through channels formed by directly interacting cells. We evaluated the efficiency of long-range communication in these channels by examining the rate of information transmission. Our study identified the stochastic phenomena that reduce this rate: front propagation failure, new front spawning, and variability in the front velocity. We found that a trade-off between the frequencies of propagation failures and new front spawning determines the optimal channel width (which geometrically determines the front length). The optimal frequency of initiating new waves is determined by a trade-off between the input information rate (higher with more frequent initiation) and the fidelity of information transmission (lower with more frequent initiation). Our analysis provides insight into the relative timescales of intra- and intercellular processes necessary for successful wave propagation.

Macrophage memory emerges from coordinated transcription factor and chromatin dynamics

Cells of the immune system operate in dynamic microenvironments where the timing, concentration, and order of signaling molecules constantly change. Despite this complexity, immune cells manage to communicate accurately and control inflammation and infection. It is unclear how these dynamic signals are encoded and decoded and if individual cells retain the memory of past exposure to inflammatory molecules. Here, we use live-cell analysis, ATAC sequencing, and an in vivo model of sepsis to show that sequential inflammatory signals induce memory in individual macrophages through reprogramming the nuclear factor κB (NF-κB) network and the chromatin accessibility landscape. We use transcriptomic profiling and deep learning to show that transcription factor and chromatin dynamics coordinate fine-tuned macrophage responses to new inflammatory signals. This work demonstrates how macrophages retain the memory of previous signals despite single-cell variability and elucidates the mechanisms of signal-induced memory in dynamic inflammatory conditions like sepsis.

Oncogenic PIK3CA corrupts growth factor signaling specificity

Technical limitations have prevented understanding of how growth factor signals are encoded in distinct activity patterns of the phosphoinositide 3-kinase (PI3K)/AKT pathway, and how this is altered by oncogenic pathway mutations. We introduce a kinetic, single-cell framework for precise calculations of PI3K-specific information transfer for different growth factors. This features live-cell imaging of PI3K/AKT activity reporters and multiplexed CyTOF measurements of PI3K/AKT and RAS/ERK signaling markers over time. Using this framework, we found that the PIK3CAH1047R oncogene was not a simple, constitutive activator of the pathway as often presented. Dose-dependent expression of PIK3CAH1047R in human cervical cancer and induced pluripotent stem cells corrupted the fidelity of growth factor-induced information transfer, with preferential amplification of epidermal growth factor receptor (EGFR) signaling responses compared to insulin-like growth factor 1 (IGF1) and insulin receptor signaling. PIK3CAH1047R did not only shift these responses to a higher mean but also enhanced signaling heterogeneity. We conclude that oncogenic PIK3CAH1047R corrupts information transfer in a growth factor-dependent manner and suggest new opportunities for tuning of receptor-specific PI3K pathway outputs for therapeutic benefit.

Emergence of temporal noise hierarchy in co-regulated genes of multi-output feed-forward loop

Natural variations in gene expression, called noise, are fundamental to biological systems. The expression noise can be beneficial or detrimental to cellular functions. While the impact of noise on individual genes is well-established, our understanding of how noise behaves when multiple genes are co-expressed by shared regulatory elements within transcription networks remains elusive. This lack of understanding extends to how the architecture and regulatory features of these networks influence noise. To address this gap, we study the multi-output feed-forward loop motif. The motif is prevalent in bacteria and yeast and influences co-expression of multiple genes by shared transcription factors (TFs). Focusing on a two-output variant of the motif, the present study explores the interplay between its architecture, co-expression (symmetric and asymmetric) patterns of the two genes, and the associated noise dynamics. We employ a stochastic modeling approach to investigate how the binding affinities of the TFs influence symmetric and asymmetric expression patterns and the resulting noise dynamics in the co-expressed genes. This knowledge could guide the development of strategies for manipulating gene expression patterns through targeted modulation of TF binding affinities.

Dynamics of positional information in the vertebrate neural tube

In developing embryos, cells acquire distinct identities depending on their position in a tissue. Secreted signalling molecules, known as morphogens, act as long-range cues to provide the spatial information that controls these cell fate decisions. In several tissues, both the level and the duration of morphogen signalling appear to be important for determining cell fates. This is the case in the forming vertebrate nervous system where antiparallel morphogen gradients pattern the dorsal-ventral axis by partitioning the tissue into sharply delineated domains of molecularly distinct neural progenitors. How information in the gradients is decoded to generate precisely positioned boundaries of gene expression remains an open question. Here, we adopt tools from information theory to quantify the positional information in the neural tube and investigate how temporal changes in signalling could influence positional precision. The results reveal that the use of signalling dynamics, as well as the signalling level, substantially increases the precision possible for the estimation of position from morphogen gradients. This analysis links the dynamics of opposing morphogen gradients with precise pattern formation and provides an explanation for why time is used to impart positional information.

Quantification of cell cycle re-entry during dedifferentiation of primary adipocytes in vitro

Dedifferentiated adipose tissue (DFAT) has been proposed as a promising source of patient-specific multipotent progenitor cells (MPPs). During induced dedifferentiation, adipocytes exhibit profound gene expression and cell morphology changes. However, dedifferentiation of post-mitotic cells is expected to enable proliferation, which is critical if enough MPPs are to be obtained. Here, lineage tracing was employed to quantify cell proliferation in mouse adipocytes subjected to a dedifferentiation-inducing protocol commonly used to obtain DFAT cells. No evidence of cell proliferation in adipocyte-derived cells was observed, in contrast to the robust proliferation of non-adipocyte cells present in adipose tissue. We conclude that proliferative MPPs derived using the ceiling culture method most likely arise from non-adipocyte cells in adipose tissue.

Single-cell stimulus-response gene expression trajectories reveal the stimulus specificities of dynamic responses by single macrophages

Macrophages induce the expression of hundreds of genes in response to immune threats. However, current technology limits our ability to capture single-cell inducible gene expression dynamics. Here, we generated high-resolution time series single-cell RNA sequencing (scRNA-seq) data from mouse macrophages responding to six stimuli, and imputed ensembles of real-time single-cell gene expression trajectories (scGETs). We found that dynamic information contained in scGETs substantially contributes to macrophage stimulus-response specificity (SRS). Dynamic information also identified correlations between immune response genes, indicating biological coordination. Furthermore, we showed that the microenvironmental context of polarizing cytokines profoundly affects scGETs, such that measuring response dynamics offered clearer discrimination of the polarization state of individual macrophage cells than single time-point measurements. Our findings highlight the important contribution of dynamic information contained in the single-cell expression responses of immune genes in characterizing the SRS and functional states of macrophages.

Distinct negative-sense RNA viruses induce a common set of transcripts encoding proteins forming an extensive network

The large group of negative-strand RNA viruses (NSVs) comprises many important pathogens. To identify conserved patterns in host responses, we systematically compared changes in the cellular RNA levels after infection of human hepatoma cells with nine different NSVs of different virulence degrees. RNA sequencing experiments indicated that the amount of viral RNA in host cells correlates with the number of differentially expressed host cell transcripts. Time-resolved differential gene expression analysis revealed a common set of 178 RNAs that are regulated by all NSVs analyzed. A newly developed open access web application allows downloads and visualizations of all gene expression comparisons for individual viruses over time or between several viruses. Most of the genes included in the core set of commonly differentially expressed genes (DEGs) encode proteins that serve as membrane receptors, signaling proteins and regulators of transcription. They mainly function in signal transduction and control immunity, metabolism, and cell survival. One hundred sixty-five of the DEGs encode host proteins from which 47 have already been linked to the regulation of viral infections in previous studies and 89 proteins form a complex interaction network that may function as a core hub to control NSV infections.IMPORTANCEThe infection of cells with negative-strand RNA viruses leads to the differential expression of many host cell RNAs. The differential spectrum of virus-regulated RNAs reflects a large variety of events including anti-viral responses, cell remodeling, and cell damage. Here, these virus-specific differences and similarities in the regulated RNAs were measured in a highly standardized model. A newly developed app allows interested scientists a wide range of comparisons and visualizations.

Pathway to independence: perspectives on the future

In this Perspective, our 2024 Pathway to Independence Fellows provide their thoughts on the future of their field. Covering topics as diverse as plant development, tissue engineering and adaptation to climate change, and using a wide range of experimental organisms, these talented postdocs showcase some of the major open questions and key challenges across the spectrum of developmental biology research.

Dynamical and combinatorial coding by MAPK p38 and NFκB in the inflammatory response of macrophages

Macrophages sense pathogens and orchestrate specific immune responses. Stimulus specificity is thought to be achieved through combinatorial and dynamical coding by signaling pathways. While NFκB dynamics are known to encode stimulus information, dynamical coding in other signaling pathways and their combinatorial coordination remain unclear. Here, we established live-cell microscopy to investigate how NFκB and p38 dynamics interface in stimulated macrophages. Information theory and machine learning revealed that p38 dynamics distinguish cytokine TNF from pathogen-associated molecular patterns and high doses from low, but contributed little to information-rich NFκB dynamics when both pathways are considered. This suggests that immune response genes benefit from decoding immune signaling dynamics or combinatorics, but not both. We found that the heterogeneity of the two pathways is surprisingly uncorrelated. Mathematical modeling revealed potential sources of uncorrelated heterogeneity in the branched pathway network topology and predicted it to drive gene expression variability. Indeed, genes dependent on both p38 and NFκB showed high scRNAseq variability and bimodality. These results identify combinatorial signaling as a mechanism to restrict NFκB-AND-p38-responsive inflammatory cytokine expression to few cells.

An intermediate Rb-E2F activity state safeguards proliferation commitment

Tissue repair, immune defence and cancer progression rely on a vital cellular decision between quiescence and proliferation1,2. Mammalian cells proliferate by triggering a positive feedback mechanism3,4. The transcription factor E2F activates cyclin-dependent kinase 2 (CDK2), which in turn phosphorylates and inactivates the E2F inhibitor protein retinoblastoma (Rb). This action further increases E2F activity to express genes needed for proliferation. Given that positive feedback can inadvertently amplify small signals, understanding how cells keep this positive feedback in check remains a puzzle. Here we measured E2F and CDK2 signal changes in single cells and found that the positive feedback mechanism engages only late in G1 phase. Cells spend variable and often extended times in a reversible state of intermediate E2F activity before committing to proliferate. This intermediate E2F activity is proportional to the amount of phosphorylation of a conserved T373 residue in Rb that is mediated by CDK2 or CDK4/CDK6. Such T373-phosphorylated Rb remains bound on chromatin but dissociates from it once Rb is hyperphosphorylated at many sites, which fully activates E2F. The preferential initial phosphorylation of T373 can be explained by its relatively slower rate of dephosphorylation. Together, our study identifies a primed state of intermediate E2F activation whereby cells sense external and internal signals and decide whether to reverse and exit to quiescence or trigger the positive feedback mechanism that initiates cell proliferation.

Reliable ligand discrimination in stochastic multistep kinetic proofreading: First passage time vs. product counting strategies

Cellular signaling, crucial for biological processes like immune response and homeostasis, relies on specificity and fidelity in signal transduction to accurately respond to stimuli amidst biological noise. Kinetic proofreading (KPR) is a key mechanism enhancing signaling specificity through time-delayed steps, although its effectiveness is debated due to intrinsic noise potentially reducing signal fidelity. In this study, we reformulate the theory of kinetic proofreading (KPR) by convolving multiple intermediate states into a single state and then define an overall "processing" time required to traverse these states. This simplification allows us to succinctly describe kinetic proofreading in terms of a single waiting time parameter, facilitating a more direct evaluation and comparison of KPR performance across different biological contexts such as DNA replication and T cell receptor (TCR) signaling. We find that loss of fidelity for longer proofreading steps relies on the specific strategy of information extraction and show that in the first-passage time (FPT) discrimination strategy, longer proofreading steps can exponentially improve the accuracy of KPR at the cost of speed. Thus, KPR can still be an effective discrimination mechanism in the high noise regime. However, in a product concentration-based discrimination strategy, longer proofreading steps do not necessarily lead to an increase in performance. However, by introducing activation thresholds on product concentrations, can we decompose the product-based strategy into a series of FPT-based strategies to better resolve the subtleties of KPR-mediated product discrimination. Our findings underscore the importance of understanding KPR in the context of how information is extracted and processed in the cell.

Co-imaging of RelA and c-Rel reveals features of NF-κB signaling for ligand discrimination

Individual cell sensing of external cues has evolved through the temporal patterns in signaling. Since nuclear factor κB (NF-κB) signaling dynamics have been examined using a single subunit, RelA, it remains unclear whether more information might be transmitted via other subunits. Using NF-κB double-knockin reporter mice, we monitored both canonical NF-κB subunits, RelA and c-Rel, simultaneously in single macrophages by quantitative live-cell imaging. We show that signaling features of RelA and c-Rel convey more information about the stimuli than those of either subunit alone. Machine learning is used to predict the ligand identity accurately based on RelA and c-Rel signaling features without considering the co-activated factors. Ligand discrimination is achieved through selective non-redundancy of RelA and c-Rel signaling dynamics, as well as their temporal coordination. These results suggest a potential role of c-Rel in fine-tuning immune responses and highlight the need for approaches that will elucidate the mechanisms regulating NF-κB subunit specificity.

Unraveling IFN-I response dynamics and TNF crosstalk in the pathophysiology of systemic lupus erythematosus

Introduction

The innate immune system serves the crucial first line of defense against a wide variety of potential threats, during which the production of pro-inflammatory cytokines IFN-I and TNFα are key. This astonishing power to fight invaders, however, comes at the cost of risking IFN-I-related pathologies, such as observed during autoimmune diseases, during which IFN-I and TNFα response dynamics are dysregulated. Therefore, these response dynamics must be tightly regulated, and precisely matched with the potential threat. This regulation is currently far from understood.

Methods

Using droplet-based microfluidics and ODE modeling, we studied the fundamentals of single-cell decision-making upon TLR signaling in human primary immune cells (n = 23). Next, using biologicals used for treating autoimmune diseases [i.e., anti-TNFα, and JAK inhibitors], we unraveled the crosstalk between IFN-I and TNFα signaling dynamics. Finally, we studied primary immune cells isolated from SLE patients (n = 8) to provide insights into SLE pathophysiology.

Results

single-cell IFN-I and TNFα response dynamics display remarkable differences, yet both being highly heterogeneous. Blocking TNFα signaling increases the percentage of IFN-I-producing cells, while blocking IFN-I signaling decreases the percentage of TNFα-producing cells. Single-cell decision-making in SLE patients is dysregulated, pointing towards a dysregulated crosstalk between IFN-I and TNFα response dynamics.

Discussion

We provide a solid droplet-based microfluidic platform to study inherent immune secretory behaviors, substantiated by ODE modeling, which can challenge the conceptualization within and between different immune signaling systems. These insights will build towards an improved fundamental understanding on single-cell decision-making in health and disease.

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.

Information Transmission in G Protein-Coupled Receptors

G protein-coupled receptors (GPCRs) are the largest class of receptors in the human genome and constitute about 30% of all drug targets. In this article, intended for a non-mathematical audience, both experimental observations and new theoretical results are compared in the context of information transmission across the cell membrane. The amount of information actually currently used or projected to be used in clinical settings is a small fraction of the information transmission capacity of the GPCR. This indicates that the number of yet undiscovered drug targets within GPCRs is much larger than what is currently known. Theoretical studies with some experimental validation indicate that localized heat deposition and dissipation are key to the identification of sites and mechanisms for drug action.

Information Transmission via Molecular Communication in Astrobiological Environments

The ubiquity of information transmission via molecular communication between cells is comprehensively documented on Earth; this phenomenon might even have played a vital role in the origin(s) and early evolution of life. Motivated by these considerations, a simple model for molecular communication entailing the diffusion of signaling molecules from transmitter to receiver is elucidated. The channel capacity C (maximal rate of information transmission) and an optimistic heuristic estimate of the actual information transmission rate ℐ are derived for this communication system; the two quantities, especially the latter, are demonstrated to be broadly consistent with laboratory experiments and more sophisticated theoretical models. The channel capacity exhibits a potentially weak dependence on environmental parameters, whereas the actual information transmission rate may scale with the intercellular distance d as ℐ ∝ d-4 and could vary substantially across settings. These two variables are roughly calculated for diverse astrobiological environments, ranging from Earth's upper oceans (C ∼ 3.1 × 103 bits/s; ℐ ∼ 4.7 × 10-2 bits/s) and deep sea hydrothermal vents (C ∼ 4.2 × 103 bits/s; ℐ ∼ 1.2 × 10-1 bits/s) to the hydrocarbon lakes and seas of Titan (C ∼ 3.8 × 103 bits/s; ℐ ∼ 2.6 × 10-1 bits/s).

Small transcriptional differences among cell clones lead to distinct NF-κB dynamics

Transcription factor dynamics is fundamental to determine the activation of accurate transcriptional programs and yet is heterogeneous at a single-cell level, even within homogeneous populations. We asked how such heterogeneity emerges for the nuclear factor κB (NF-κB). We found that clonal populations of immortalized fibroblasts derived from a single mouse embryo display robustly distinct NF-κB dynamics upon tumor necrosis factor ɑ (TNF-ɑ) stimulation including persistent, oscillatory, and weak activation, giving rise to differences in the transcription of its targets. By combining transcriptomics and simulations we show how less than two-fold differences in the expression levels of genes coding for key proteins of the signaling cascade and feedback system are predictive of the differences of the NF-κB dynamic response of the clones to TNF-ɑ and IL-1β. We propose that small transcriptional differences in the regulatory circuit of a transcription factor can lead to distinct signaling dynamics in cells within homogeneous cell populations and among different cell types.

A guide to ERK dynamics, part 1: mechanisms and models

Extracellular signal-regulated kinase (ERK) has long been studied as a key driver of both essential cellular processes and disease. A persistent question has been how this single pathway is able to direct multiple cell behaviors, including growth, proliferation, and death. Modern biosensor studies have revealed that the temporal pattern of ERK activity is highly variable and heterogeneous, and critically, that these dynamic differences modulate cell fate. This two-part review discusses the current understanding of dynamic activity in the ERK pathway, how it regulates cellular decisions, and how these cell fates lead to tissue regulation and pathology. In part 1, we cover the optogenetic and live-cell imaging technologies that first revealed the dynamic nature of ERK, as well as current challenges in biosensor data analysis. We also discuss advances in mathematical models for the mechanisms of ERK dynamics, including receptor-level regulation, negative feedback, cooperativity, and paracrine signaling. While hurdles still remain, it is clear that higher temporal and spatial resolution provide mechanistic insights into pathway circuitry. Exciting new algorithms and advanced computational tools enable quantitative measurements of single-cell ERK activation, which in turn inform better models of pathway behavior. However, the fact that current models still cannot fully recapitulate the diversity of ERK responses calls for a deeper understanding of network structure and signal transduction in general.

Information-Theoretic Analysis of a Model of CAR-4-1BB-Mediated NFκB Activation

Systems biology utilizes computational approaches to examine an array of biological processes, such as cell signaling, metabolomics and pharmacology. This includes mathematical modeling of CAR T cells, a modality of cancer therapy by which genetically engineered immune cells recognize and combat a cancerous target. While successful against hematologic malignancies, CAR T cells have shown limited success against other cancer types. Thus, more research is needed to understand their mechanisms of action and leverage their full potential. In our work, we set out to apply information theory on a mathematical model of NFκB signaling initiated by the CAR following antigen encounter. First, we estimated channel capacity for CAR-4-1BB-mediated NFκB signal transduction. Next, we evaluated the pathway's ability to distinguish contrasting "low" and "high" antigen concentration levels, depending on the amount of variability in protein concentrations. Finally, we assessed the fidelity by which NFκB activation reflects the encountered antigen concentration, depending on the prevalence of antigen-positive targets in tumor population. We found that in most scenarios, fold change in the nuclear concentration of NFκB carries a higher channel capacity for the pathway than NFκB's absolute response. Additionally, we found that most errors in transducing the antigen signal through the pathway skew towards underestimating the concentration of encountered antigen. Finally, we found that disabling IKKβ deactivation could increase signaling fidelity against targets with antigen-negative cells. Our information-theoretic analysis of signal transduction can provide novel perspectives on biological signaling, as well as enable a more informed path to cell engineering.Kindly check and confirm whether the corresponding affiliation is correctly identified.this is correct.

CD3ζ ITAMs enable ligand discrimination and antagonism by inhibiting TCR signaling in response to low-affinity peptides

The T cell antigen receptor (TCR) contains ten immunoreceptor tyrosine-based activation motif (ITAM) signaling sequences distributed within six CD3 subunits; however, the reason for such structural complexity and multiplicity is unclear. Here we evaluated the effect of inactivating the three CD3ζ chain ITAMs on TCR signaling and T cell effector responses using a conditional 'switch' mouse model. Unexpectedly, we found that T cells expressing TCRs containing inactivated (non-signaling) CD3ζ ITAMs (6F-CD3ζ) exhibited reduced ability to discriminate between low- and high-affinity ligands, resulting in enhanced signaling and cytokine responses to low-affinity ligands because of a previously undetected inhibitory function of CD3ζ ITAMs. Also, 6F-CD3ζ TCRs were refractory to antagonism, as predicted by a new in silico adaptive kinetic proofreading model that revises the role of ITAM multiplicity in TCR signaling. Finally, T cells expressing 6F-CD3ζ displayed enhanced cytolytic activity against solid tumors expressing low-affinity ligands, identifying a new counterintuitive approach to TCR-mediated cancer immunotherapy.

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.

Heterogeneous off-target impact of ion-channel deletion on intrinsic properties of hippocampal model neurons that self-regulate calcium

How do neurons that implement cell-autonomous self-regulation of calcium react to knockout of individual ion-channel conductances? To address this question, we used a heterogeneous population of 78 conductance-based models of hippocampal pyramidal neurons that maintained cell-autonomous calcium homeostasis while receiving theta-frequency inputs. At calcium steady-state, we individually deleted each of the 11 active ion-channel conductances from each model. We measured the acute impact of deleting each conductance (one at a time) by comparing intrinsic electrophysiological properties before and immediately after channel deletion. The acute impact of deleting individual conductances on physiological properties (including calcium homeostasis) was heterogeneous, depending on the property, the specific model, and the deleted channel. The underlying many-to-many mapping between ion channels and properties pointed to ion-channel degeneracy. Next, we allowed the other conductances (barring the deleted conductance) to evolve towards achieving calcium homeostasis during theta-frequency activity. When calcium homeostasis was perturbed by ion-channel deletion, post-knockout plasticity in other conductances ensured resilience of calcium homeostasis to ion-channel deletion. These results demonstrate degeneracy in calcium homeostasis, as calcium homeostasis in knockout models was implemented in the absence of a channel that was earlier involved in the homeostatic process. Importantly, in reacquiring homeostasis, ion-channel conductances and physiological properties underwent heterogenous plasticity (dependent on the model, the property, and the deleted channel), even introducing changes in properties that were not directly connected to the deleted channel. Together, post-knockout plasticity geared towards maintaining homeostasis introduced heterogenous off-target effects on several channels and properties, suggesting that extreme caution be exercised in interpreting experimental outcomes involving channel knockouts.

Dynamical Information Synergy in Biochemical Signaling Networks

Biological cells encode information about their environment through biochemical signaling networks that control their internal state and response. This information is often encoded in the dynamical patterns of the signaling molecules, rather than just their instantaneous concentrations. Here, we analytically calculate the information contained in these dynamics for a number of paradigmatic cases in the linear regime, for both static and time-dependent input signals. When considering oscillatory output dynamics, we report on the emergence of synergy between successive measurements, meaning that the joint information in two measurements exceeds the sum of the individual information. We extend our analysis numerically beyond the scope of linear input encoding to reveal synergetic effects in the cases of frequency or damping modulation, both of which are relevant to classical biochemical signaling systems.

Integrating glycolysis, citric acid cycle, pentose phosphate pathway, and fatty acid beta-oxidation into a single computational model

The metabolic network of a living cell is highly intricate and involves complex interactions between various pathways. In this study, we propose a computational model that integrates glycolysis, the pentose phosphate pathway (PPP), the fatty acids beta-oxidation, and the tricarboxylic acid cycle (TCA cycle) using queueing theory. The model utilizes literature data on metabolite concentrations and enzyme kinetic constants to calculate the probabilities of individual reactions occurring on a microscopic scale, which can be viewed as the reaction rates on a macroscopic scale. However, it should be noted that the model has some limitations, including not accounting for all the reactions in which the metabolites are involved. Therefore, a genetic algorithm (GA) was used to estimate the impact of these external processes. Despite these limitations, our model achieved high accuracy and stability, providing real-time observation of changes in metabolite concentrations. This type of model can help in better understanding the mechanisms of biochemical reactions in cells, which can ultimately contribute to the prevention and treatment of aging, cancer, metabolic diseases, and neurodegenerative disorders.

Improved approaches to channel capacity estimation discover compromised GPCR signaling in diverse cancer cells

Intracellular signaling orchestrates an organism's development and functioning and underlies various pathologies, such as cancer, when aberrant. A universal cell signaling characteristic is channel capacity - the measure of how much information a given transmitting system can reliably transduce. Here, we describe improved approaches to quantify GPCR signaling channel capacity in single cells, averaged across cell population. We assess the channel capacity based on distribution of residuals by the cellular response amplitude. We further develop means to handle irregularly responding cancer cells using the integral values of their response to different agonist concentrations. These approaches enabled us to analyze, for the first time, channel capacity in single cancer cells. A universal feature emerging for different cancer cell types is a decreased channel capacity of their GPCR signaling. These findings provide experimental validation to the hypothesis that cancer is an information disease, bearing importance for basic cancer biology and anticancer drug discovery.

NFAT and NF-κB dynamically co-regulate TCR and CAR signaling responses in human T cells

While it has been established that the responses of T cells to antigens are combinatorially regulated by multiple signaling pathways, it remains elusive what mechanisms cells utilize to quantitatively modulate T cell responses during pathway integration. Here, we show that two key pathways in T cell signaling, calcium/nuclear factor of activated T cells (NFAT) and protein kinase C (PKC)/nuclear factor κB (NF-κB), integrate through a dynamic and combinatorial strategy to fine-tune T cell response genes. At the cis-regulatory level, the two pathways integrate through co-binding of NFAT and NF-κB to immune response genes. Pathway integration is further regulated temporally, where T cell receptor (TCR) and chimeric antigen receptor (CAR) activation signals modulate the temporal relationships between the nuclear localization dynamics of NFAT and NF-κB. Such physical and temporal integrations together contribute to distinct modes of expression modulation for genes. Thus, the temporal relationships between regulators can be modulated to affect their co-targets during immune responses, underscoring the importance of dynamic combinatorial regulation in cellular signaling.

Nearly maximal information gain due to time integration in central dogma reactions

Living cells process information about their environment through the central dogma processes of transcription and translation, which drive the cellular response to stimuli. Here, we study the transfer of information from environmental input to the transcript and protein expression levels. Evaluation of both experimental and analogous simulation data reveals that transcription and translation are not two simple information channels connected in series. Instead, we demonstrate that the central dogma reactions often create a time-integrating information channel, where the translation channel receives and integrates multiple outputs from the transcription channel. This information channel model of the central dogma provides new information-theoretic selection criteria for the central dogma rate constants. Using the data for four well-studied species we show that their central dogma rate constants achieve information gain because of time integration while also keeping the loss because of stochasticity in translation relatively low (<0.5 bits).

Model of ligand-triggered information transmission in G-protein coupled receptor complexes

We present a model for the effects of ligands on information transmission in G-Protein Coupled Receptor (GPCR) complexes. The model is built ab initio entirely on principles of statistical mechanics and tenets of information transmission theory and was validated in part using agonist-induced effector activity and signaling bias for the angiotensin- and adrenergic-mediated signaling pathways, with in vitro observations of phosphorylation sites on the C tail of the GPCR complex, and single-cell information-transmission experiments. The model extends traditional kinetic models that form the basis for many existing models of GPCR signaling. It is based on maximizing the rates of entropy production and information transmission through the GPCR complex. The model predicts that (1) phosphatase-catalyzed reactions, as opposed to kinase-catalyzed reactions, on the C-tail and internal loops of the GPCR are responsible for controlling the signaling activity, (2) signaling favors the statistical balance of the number of switches in the ON state and the number in the OFF state, and (3) biased-signaling response depends discontinuously on ligand concentration.

The MAPK/ERK channel capacity exceeds 6 bit/hour

Living cells utilize signaling pathways to sense, transduce, and process information. As the extracellular stimulation often has rich temporal characteristics which may govern dynamic cellular responses, it is important to quantify the rate of information flow through the signaling pathways. In this study, we used an epithelial cell line expressing a light-activatable FGF receptor and an ERK activity reporter to assess the ability of the MAPK/ERK pathway to transduce signal encoded in a sequence of pulses. By stimulating the cells with random light pulse trains, we demonstrated that the MAPK/ERK channel capacity is at least 6 bits per hour. The input reconstruction algorithm detects the light pulses with 1-min accuracy 5 min after their occurrence. The high information transmission rate may enable the pathway to coordinate multiple processes including cell movement and respond to rapidly varying stimuli such as chemoattracting gradients created by other cells.

How Well Can We Infer Selection Benefits and Mutation Rates from Allele Frequencies?

Experimentalists observe allele frequency distributions and try to infer mutation rates and selection coefficients. How easy is this? We calculate limits to their ability in the context of the Wright-Fisher model by first finding the maximal amount of information that can be acquired using allele frequencies about the mutation rate and selection coefficient- at least 2 bits per allele- and then by finding how the organisms would have shaped their mutation rates and selection coefficients so as to maximize the information transfer.

An optimized reporter of the transcription factor hypoxia-inducible factor 1α reveals complex HIF-1α activation dynamics in single cells

Immune cells adopt a variety of metabolic states to support their many biological functions, which include fighting pathogens, removing tissue debris, and tissue remodeling. One of the key mediators of these metabolic changes is the transcription factor hypoxia-inducible factor 1α (HIF-1α). Single-cell dynamics have been shown to be an important determinant of cell behavior; however, despite the importance of HIF-1α, little is known about its single-cell dynamics or their effect on metabolism. To address this knowledge gap, here we optimized a HIF-1α fluorescent reporter and applied it to study single-cell dynamics. First, we showed that single cells are likely able to differentiate multiple levels of prolyl hydroxylase inhibition, a marker of metabolic change, via HIF-1α activity. We then applied a physiological stimulus known to trigger metabolic change, interferon-γ, and observed heterogeneous, oscillatory HIF-1α responses in single cells. Finally, we input these dynamics into a mathematical model of HIF-1α-regulated metabolism and discovered a profound difference between cells exhibiting high versus low HIF-1α activation. Specifically, we found cells with high HIF-1α activation are able to meaningfully reduce flux through the tricarboxylic acid cycle and show a notable increase in the NAD+/NADH ratio compared with cells displaying low HIF-1α activation. Altogether, this work demonstrates an optimized reporter for studying HIF-1α in single cells and reveals previously unknown principles of HIF-1α activation.

Quantifying stimulus-response specificity to probe the functional state of macrophages

Immune sentinel macrophages initiate responses to pathogens via hundreds of immune response genes. Each immune threat demands a tailored response, suggesting that the capacity for stimulus-specific gene expression is a key functional hallmark of healthy macrophages. To quantify this property, termed "stimulus-response specificity" (SRS), we developed a single-cell experimental workflow and analytical approaches based on information theory and machine learning. We found that the response specificity of macrophages is driven by combinations of specific immune genes that show low cell-to-cell heterogeneity and are targets of separate signaling pathways. The "response specificity profile," a systematic comparison of multiple stimulus-response distributions, was distinctly altered by polarizing cytokines, and it enabled an assessment of the functional state of macrophages. Indeed, the response specificity profile of peritoneal macrophages from old and obese mice showed characteristic differences, suggesting that SRS may be a basis for measuring the functional state of innate immune cells. A record of this paper's transparent peer review process is included in the supplemental information.

Information transfer in mammalian glycan-based communication

Glycan-binding proteins, so-called lectins, are exposed on mammalian cell surfaces and decipher the information encoded within glycans translating it into biochemical signal transduction pathways in the cell. These glycan-lectin communication pathways are complex and difficult to analyze. However, quantitative data with single-cell resolution provide means to disentangle the associated signaling cascades. We chose C-type lectin receptors (CTLs) expressed on immune cells as a model system to study their capacity to transmit information encoded in glycans of incoming particles. In particular, we used nuclear factor kappa-B-reporter cell lines expressing DC-specific ICAM-3-grabbing nonintegrin (DC-SIGN), macrophage C-type lectin (MCL), dectin-1, dectin-2, and macrophage-inducible C-type lectin (MINCLE), as well as TNFαR and TLR-1&2 in monocytic cell lines and compared their transmission of glycan-encoded information. All receptors transmit information with similar signaling capacity, except dectin-2. This lectin was identified to be less efficient in information transmission compared to the other CTLs, and even when the sensitivity of the dectin-2 pathway was enhanced by overexpression of its co-receptor FcRγ, its transmitted information was not. Next, we expanded our investigation toward the integration of multiple signal transduction pathways including synergistic lectins, which is crucial during pathogen recognition. We show how the signaling capacity of lectin receptors using a similar signal transduction pathway (dectin-1 and dectin-2) is being integrated by compromising between the lectins. In contrast, co-expression of MCL synergistically enhanced the dectin-2 signaling capacity, particularly at low-glycan stimulant concentration. By using dectin-2 and other lectins as examples, we demonstrate how signaling capacity of dectin-2 is modulated in the presence of other lectins, and therefore, the findings provide insight into how immune cells translate glycan information using multivalent interactions.

Epigenetics as a mediator of plasticity in cancer

The concept of an epigenetic landscape describing potential cellular fates arising from pluripotent cells, first advanced by Conrad Waddington, has evolved in light of experiments showing nondeterministic outcomes of regulatory processes and mathematical methods for quantifying stochasticity. In this Review, we discuss modern approaches to epigenetic and gene regulation landscapes and the associated ideas of entropy and attractor states, illustrating how their definitions are both more precise and relevant to understanding cancer etiology and the plasticity of cancerous states. We address the interplay between different types of regulatory landscapes and how their changes underlie cancer progression. We also consider the roles of cellular aging and intrinsic and extrinsic stimuli in modulating cellular states and how landscape alterations can be quantitatively mapped onto phenotypic outcomes and thereby used in therapy development.

Radioresistance and radiosensitivity: a biophysical approach on bacterial cells robustness

The study of radiosensitivity and radioresistance of organisms exposed to ionizing radiation has acquired additional relevance since a new bio-concept, coined as The primacy of Proteome over Genome, was proposed and demonstrated elsewhere a few years ago. According to that finding, genome integrity would require an actively functioning Proteome. However, when exposure to radiation takes place, Reactive Oxygen Species (ROS) from water radiolysis induce protein carbonylation (PC), an irreversible oxidative Proteome damage. The bio-models used in that study were the radiosensitive Escherichia coli and the extraordinarily robust Deinococcus radiodurans. The production of ROS induces protective reactions rendering them non-reactive forms. Protective entities present in the cytosol, moieties smaller than 3 kDa, shield the Proteome against ROS, yielding protection against carbonylation. Shown in the present study is the fact that the fate of proteins functionality is determined by the magnitude of the Protein Carbonylation Yield (Y_PC), a quantity here analytically defined using published Y_PC numerical results. Analytical Y_PC expressions for E. coli and D. radiodurans were the input for a phenomenological approach, where the radiobiological magnitudes P_P and P_N, the probabilities for production of protein damage and ROS neutralization, respectively, were also analytically deduced. These highly relevant magnitudes, associated with key radiosensitivity and radioresistance issues, are addressed and discussed in this study. Among the plethora of information and conclusions derived from the present study, those endowed with higher conceptual degree, vis-à-vis the "Primacy of Proteome over Genome" concept, are as follows: (1) the ROS neutralization process in D. radiodurans reaches a maximum at a dose interval corresponding to the repairing shoulder. Therefore, it is a signature of the higher efficiency of the PC neutralization process. (2) ROS neutralization in D. radiodurans is nearly one order of magnitude higher than in E. coli, thus accounting for its extraordinary radioresistance. (3) Both physical (ROS-induced carbonyl radicals) and biological (protein modifications) processes are imbedded in the Protein Carbonylation Yield. The amalgamation of these two processes was accomplished by means of a statistical formalism.

Design principles of improving the dose-response alignment in coupled GTPase switches

"Dose-response alignment" (DoRA), where the downstream response of cellular signaling pathways closely matches the fraction of activated receptor, can improve the fidelity of dose information transmission. The negative feedback has been experimentally identified as a key component for DoRA, but numerical simulations indicate that negative feedback is not sufficient to achieve perfect DoRA, i.e., perfect match of downstream response and receptor activation level. Thus a natural question is whether there exist design principles for signaling motifs within only negative feedback loops to improve DoRA to near-perfect DoRA. Here, we investigated several model formulations of an experimentally validated circuit that couples two molecular switches-mGTPase (monomeric GTPase) and tGTPase (heterotrimeric GTPases) - with negative feedback loops. In the absence of feedback, the low and intermediate mGTPase activation levels benefit DoRA in mass action and Hill-function models, respectively. Adding negative feedback has versatile roles on DoRA: it may impair DoRA in the mass action model with low mGTPase activation level and Hill-function model with intermediate mGTPase activation level; in other cases, i.e., the mass action model with a high mGTPase activation level or the Hill-function model with a non-intermediate mGTPase activation level, it improves DoRA. Furthermore, we found that DoRA in a longer cascade (i.e., tGTPase) can be obtained using Hill-function kinetics under certain conditions. In summary, we show how ranges of activity of mGTPase, reaction kinetics, the negative feedback, and the cascade length affect DoRA. This work provides a framework for improving the DoRA performance in signaling motifs with negative feedback.

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.

NF-κB dynamics in the language of immune cells

Biological discovery has been driven by advances in throughput and resolution of analysis technologies. They have also created an indelible bias for snapshot-based knowledge. Even though recent methods such as multi-omics single-cell assays have empowered immunological investigations, they still provide snapshots of cellular behaviors and thus, have inherent limitations in reconstructing unsynchronized dynamic events across individual cells. Here, we present a rationale for how NF-κB may convey specificity of contextual information through subtle quantitative features of its signaling dynamics. The next frontier of predictive understanding should involve functional characterization of NF-κB signaling dynamics and their immunological implications. This may help solve the apparent paradox that a ubiquitously activated transcription factor can shape accurate responses to different immune challenges.

Understanding the physical mechanism of intrinsic noise inside viscous isotropic solids

We report acoustic impulse-response and system function of particle vibration inside viscous, dense solids and explain the physical mechanism of intrinsic-noise generation. With an external disturbance of a harmonic force acting on particles inside viscous solid media, the system of particle vibration goes through a gradual transition from a static state to a steady harmonic vibrational state. Based on the damped oscillator model, the transition frequency spectrum resembles the intrinsic noise generated by vibrating particles in viscous isotropic solids, which delivers a crucial understanding for applications to invert stratum characteristics around the drilled oil well and its abnormal geological structure.

Spatiotemporal NF-κB dynamics encodes the position, amplitude, and duration of local immune inputs

Infected cells communicate through secreted signaling molecules like cytokines, which carry information about pathogens. How differences in cytokine secretion affect inflammatory signaling over space and how responding cells decode information from propagating cytokines are not understood. By computationally and experimentally studying NF-κB dynamics in cocultures of signal-sending cells (macrophages) and signal-receiving cells (fibroblasts), we find that cytokine signals are transmitted by wave-like propagation of NF-κB activity and create well-defined activation zones in responding cells. NF-κB dynamics in responding cells can simultaneously encode information about cytokine dose, duration, and distance to the cytokine source. Spatially resolved transcriptional analysis reveals that responding cells transmit local cytokine information to distance-specific proinflammatory gene expression patterns, creating "gene expression zones." Despite single-cell variability, the size and duration of the signaling zone are tightly controlled by the macrophage secretion profile. Our results highlight how macrophages tune cytokine secretion to control signal transmission distance and how inflammatory signaling interprets these signals in space and time.

Efficient information coding and degeneracy in the nervous system

Efficient information coding (EIC) is a universal biological framework rooted in the fundamental principle that system responses should match their natural stimulus statistics for maximizing environmental information. Quantitatively assessed through information theory, such adaptation to the environment occurs at all biological levels and timescales. The context dependence of environmental stimuli and the need for stable adaptations make EIC a daunting task. We argue that biological complexity is the principal architect that subserves deft execution of stable EIC. Complexity in a system is characterized by several functionally segregated subsystems that show a high degree of functional integration when they interact with each other. Complex biological systems manifest heterogeneities and degeneracy, wherein structurally different subsystems could interact to yield the same functional outcome. We argue that complex systems offer several choices that effectively implement EIC and homeostasis for each of the different contexts encountered by the system.

NF-κB memory coordinates transcriptional responses to dynamic inflammatory stimuli

Many scenarios in cellular communication require cells to interpret multiple dynamic signals. It is unclear how exposure to inflammatory stimuli alters transcriptional responses to subsequent stimulus. Using high-throughput microfluidic live-cell analysis, we systematically profile the NF-κB response to different signal sequences in single cells. We find that NF-κB dynamics store the short-term history of received signals: depending on the prior pathogenic or cytokine signal, the NF-κB response to subsequent stimuli varies from no response to full activation. Using information theory, we reveal that these stimulus-dependent changes in the NF-κB response encode and reflect information about the identity and dose of the prior stimulus. Small-molecule inhibition, computational modeling, and gene expression profiling show that this encoding is driven by stimulus-dependent engagement of negative feedback modules. These results provide a model for how signal transduction networks process sequences of inflammatory stimuli to coordinate cellular responses in complex dynamic environments.

Quantifying the phenotypic information in mRNA abundance

Quantifying the dependency between mRNA abundance and downstream cellular phenotypes is a fundamental open problem in biology. Advances in multimodal single-cell measurement technologies provide an opportunity to apply new computational frameworks to dissect the contribution of individual genes and gene combinations to a given phenotype. Using an information theory approach, we analyzed multimodal data of the expression of 83 genes in the Ca2+ signaling network and the dynamic Ca2+ response in the same cell. We found that the overall expression levels of these 83 genes explain approximately 60% of Ca2+ signal entropy. The average contribution of each single gene was 17%, revealing a large degree of redundancy between genes. Using different heuristics, we estimated the dependency between the size of a gene set and its information content, revealing that on average, a set of 53 genes contains 54% of the information about Ca2+ signaling. Our results provide the first direct quantification of information content about complex cellular phenotype that exists in mRNA abundance measurements.

Information theory entering soils and tissues

For a long time, application of information theory for characterizing biological signaling fidelity was limited to "telegraph-like" situations, copying the classical scenario for which it was developed. A study published in this issue of Cell Systems applies this powerful framework to several distinct cases of dynamic signal sensing in complex geometries.

Multimodal perception links cellular state to decision making in single cells

Individual cells take decisions that are adapted to their internal state and surroundings, but how cells can reliably do this remains unclear. Using multiplexed quantification of signaling responses and markers of the cellular state, we find that signaling nodes in a network display adaptive information processing, which leads to heterogeneous growth factor responses and enables nodes to capture partially non-redundant information about the cellular state. Collectively, as a multimodal percept, this gives individual cells a large information processing capacity to accurately place growth factor concentration within the context of their cellular state and make cellular state-dependent decisions. We propose that heterogeneity and complexity in signaling networks have co-evolved to enable specific and context-aware cellular decision making in a multicellular setting.

Quantifying information of intracellular signaling: progress with machine learning

Cells convey information about their extracellular environment to their core functional machineries. Studying the capacity of intracellular signaling pathways to transmit information addresses fundamental questions about living systems. Here, we review how information-theoretic approaches have been used to quantify information transmission by signaling pathways that are functionally pleiotropic and subject to molecular stochasticity. We describe how recent advances in machine learning have been leveraged to address the challenges of complex temporal trajectory datasets and how these have contributed to our understanding of how cells employ temporal coding to appropriately adapt to environmental perturbations.

Insights on the NF-κB System Using Live Cell Imaging: Recent Developments and Future Perspectives

The transcription factor family of nuclear factor kappa B (NF-κB) proteins is widely recognized as a key player in inflammation and the immune responses, where it plays a fundamental role in translating external inflammatory cues into precise transcriptional programs, including the timely expression of a wide variety of cytokines/chemokines. Live cell imaging in single cells showed approximately 15 years ago that the canonical activation of NF-κB upon stimulus is very dynamic, including oscillations of its nuclear localization with a period close to 1.5 hours. This observation has triggered a fruitful interdisciplinary research line that has provided novel insights on the NF-κB system: how its heterogeneous response differs between cell types but also within homogeneous populations; how NF-κB dynamics translate external cues into intracellular signals and how NF-κB dynamics affects gene expression. Here we review the main features of this live cell imaging approach to the study of NF-κB, highlighting the key findings, the existing gaps of knowledge and hinting towards some of the potential future steps of this thriving research field.

Towards 'end-to-end' analysis and understanding of biological timecourse data

Petabytes of increasingly complex and multidimensional live cell and tissue imaging data are generated every year. These videos hold large promise for understanding biology at a deep and fundamental level, as they capture single-cell and multicellular events occurring over time and space. However, the current modalities for analysis and mining of these data are scattered and user-specific, preventing more unified analyses from being performed over different datasets and obscuring possible scientific insights. Here, we propose a unified pipeline for storage, segmentation, analysis, and statistical parametrization of live cell imaging datasets.