Modeling Cell-to-Cell Communication Networks Using Response-Time Distributions

Dynamic Time-Programming Circuit for Encoding Information, Programming Dissipative Systems, and Delaying Release of Cargo

Living systems have some of the most sophisticated reaction circuits in the world, realizing many incredibly complex functions through a variety of simple molecular reactions, in which the most notable feature that distinguishes them from artificial molecular reaction networks is the precise control of reaction times and programmable expression. Here, we exploit the hydrolysis-directed nature of λ exonuclease and the programmed responses of the dynamic nanotechnology of nucleic acids to construct a simple, complete, and powerful set of temporally programmed circuits. This system can arbitrarily regulate the degradation rate of the blocker, thereby delaying the nucleic acid chain substitution reaction with less signal leakage. In addition, the powerful dynamic reaction network of nucleic acids enabled us to control the programmed execution of a wide range of reactions in different fields. We have developed a simple strategy to introduce precise control of the time dimension into nucleic acid reaction circuits, which greatly enriches the functionality and applicability of the reaction programs, which can be easily used as timers, compilers, converters, etc. The simplicity, precision, stability, and versatility of such dynamic temporal programming circuits greatly expand the potential of artificial molecular reaction networks for more complex practical applications in biochemistry and molecular biology.

Diffusion-limited cytokine signaling in T cell populations

Effective immune-cell responses depend on collective decision-making mediated by diffusible intercellular signaling proteins called cytokines. Here, we designed a three-dimensional spatiotemporal modeling framework and a precise finite-element simulation setup to systematically investigate the origin and consequences of spatially inhomogeneous cytokine distributions in lymph nodes. We found that such inhomogeneities are critical for effective paracrine signaling, and they do not arise by diffusion and uptake alone, but rather depend on properties of the cell population such as an all-or-none behavior of cytokine secreting cells. Furthermore, we assessed the regulatory properties of negative and positive feedback in combination with diffusion-limited signaling dynamics, and we derived statistical quantities to characterize the spatiotemporal signaling landscape in the context of specific tissue architectures. Overall, our simulations highlight the complex spatiotemporal dynamics imposed by cell-cell signaling with diffusible ligands, which entails a large potential for fine-tuned biological control especially if combined with feedback mechanisms.

Bridging systems biology and tissue engineering: Unleashing the full potential of complex 3D in vitro tissue models of disease

Rapid advances in tissue engineering have resulted in more complex and physiologically relevant 3D in vitro tissue models with applications in fundamental biology and therapeutic development. However, the complexity provided by these models is often not leveraged fully due to the reductionist methods used to analyze them. Computational and mathematical models developed in the field of systems biology can address this issue. Yet, traditional systems biology has been mostly applied to simpler in vitro models with little physiological relevance and limited cellular complexity. Therefore, integrating these two inherently interdisciplinary fields can result in new insights and move both disciplines forward. In this review, we provide a systematic overview of how systems biology has been integrated with 3D in vitro tissue models and discuss key application areas where the synergies between both fields have led to important advances with potential translational impact. We then outline key directions for future research and discuss a framework for further integration between fields.

A practical guide to bioimaging research data management in core facilities

Bioimage data are generated in diverse research fields throughout the life and biomedical sciences. Its potential for advancing scientific progress via modern, data-driven discovery approaches reaches beyond disciplinary borders. To fully exploit this potential, it is necessary to make bioimaging data, in general, multidimensional microscopy images and image series, FAIR, that is, findable, accessible, interoperable and reusable. These FAIR principles for research data management are now widely accepted in the scientific community and have been adopted by funding agencies, policymakers and publishers. To remain competitive and at the forefront of research, implementing the FAIR principles into daily routines is an essential but challenging task for researchers and research infrastructures. Imaging core facilities, well-established providers of access to imaging equipment and expertise, are in an excellent position to lead this transformation in bioimaging research data management. They are positioned at the intersection of research groups, IT infrastructure providers, the institution´s administration, and microscope vendors. In the frame of German BioImaging - Society for Microscopy and Image Analysis (GerBI-GMB), cross-institutional working groups and third-party funded projects were initiated in recent years to advance the bioimaging community's capability and capacity for FAIR bioimage data management. Here, we provide an imaging-core-facility-centric perspective outlining the experience and current strategies in Germany to facilitate the practical adoption of the FAIR principles closely aligned with the international bioimaging community. We highlight which tools and services are ready to be implemented and what the future directions for FAIR bioimage data have to offer.

Complementing Cell Taxonomies with a Multicellular Analysis of Tissues

The application of single-cell molecular profiling coupled with spatial technologies has enabled charting of cellular heterogeneity in reference tissues and in disease. This new wave of molecular data has highlighted the expected diversity of single-cell dynamics upon shared external queues and spatial organizations. However, little is known about the relationship between single-cell heterogeneity and the emergence and maintenance of robust multicellular processes in developed tissues and its role in (patho)physiology. Here, we present emerging computational modeling strategies that use increasingly available large-scale cross-condition single-cell and spatial datasets to study multicellular organization in tissues and complement cell taxonomies. This perspective should enable us to better understand how cells within tissues collectively process information and adapt synchronized responses in disease contexts and to bridge the gap between structural changes and functions in tissues.

Noisy Delay Denoises Biochemical Oscillators

Genetic oscillations are generated by delayed transcriptional negative feedback loops, wherein repressor proteins inhibit their own synthesis after a temporal production delay. This delay is distributed because it arises from a sequence of noisy processes, including transcription, translocation, translation, and folding. Because the delay determines repression timing and, therefore, oscillation period, it has been commonly believed that delay noise weakens oscillatory dynamics. Here, we demonstrate that noisy delay can surprisingly denoise genetic oscillators. Specifically, moderate delay noise improves the signal-to-noise ratio and sharpens oscillation peaks, all without impacting period and amplitude. We show that this denoising phenomenon occurs in a variety of well-studied genetic oscillators, and we use queueing theory to uncover the universal mechanisms that produce it.

Density physics-informed neural networks reveal sources of cell heterogeneity in signal transduction

The transduction time between signal initiation and final response provides valuable information on the underlying signaling pathway, including its speed and precision. Furthermore, multi-modality in a transduction-time distribution indicates that the response is regulated by multiple pathways with different transduction speeds. Here, we developed a method called density physics-informed neural networks (Density-PINNs) to infer the transduction-time distribution from measurable final stress response time traces. We applied Density-PINNs to single-cell gene expression data from sixteen promoters regulated by unknown pathways in response to antibiotic stresses. We found that promoters with slower signaling initiation and transduction exhibit larger cell-to-cell heterogeneity in response intensity. However, this heterogeneity was greatly reduced when the response was regulated by slow and fast pathways together. This suggests a strategy for identifying effective signaling pathways for consistent cellular responses to disease treatments. Density-PINNs can also be applied to understand other time delay systems, including infectious diseases.

Timed Pulses in DNA Strand Displacement Reactions

Inspired by naturally occurring regulatory mechanisms that allow complex temporal pulse features with programmable delays, we demonstrate here a strategy to achieve temporally programmed pulse output signals in DNA-based strand displacement reactions (SDRs). To achieve this, we rationally designed input strands that, once bound to their target duplex, can be gradually degraded, resulting in a pulse output signal. We also designed blocker strands that suppress strand displacement and determine the time at which the pulse reaction is generated. We show that by controlling the degradation rate of blocker and input strands, we can finely control the delayed pulse output over a range of 10 h. We also prove that it is possible to orthogonally delay two different pulse reactions in the same solution by taking advantage of the specificity of the degradation reactions for the input and blocker strands. Finally, we show here two possible applications of such delayed pulse SDRs: the time-programmed pulse decoration of DNA nanostructures and the sequentially appearing and self-erasing formation of DNA-based patterns.

Distribution modeling quantifies collective TH cell decision circuits in chronic inflammation

Immune responses are tightly regulated by a diverse set of interacting immune cell populations. Alongside decision-making processes such as differentiation into specific effector cell types, immune cells initiate proliferation at the beginning of an inflammation, forming two layers of complexity. Here, we developed a general mathematical framework for the data-driven analysis of collective immune cell dynamics. We identified qualitative and quantitative properties of generic network motifs, and we specified differentiation dynamics by analysis of kinetic transcriptome data. Furthermore, we derived a specific, data-driven mathematical model for T helper 1 versus T follicular helper cell-fate decision dynamics in acute and chronic lymphocytic choriomeningitis virus infections in mice. The model recapitulates important dynamical properties without model fitting and solely by using measured response-time distributions. Model simulations predict different windows of opportunity for perturbation in acute and chronic infection scenarios, with potential implications for optimization of targeted immunotherapy.

A Review of Single-Cell RNA-Seq Annotation, Integration, and Cell-Cell Communication

Single-cell RNA sequencing (scRNA-seq) has emerged as a powerful tool for investigating cellular biology at an unprecedented resolution, enabling the characterization of cellular heterogeneity, identification of rare but significant cell types, and exploration of cell-cell communications and interactions. Its broad applications span both basic and clinical research domains. In this comprehensive review, we survey the current landscape of scRNA-seq analysis methods and tools, focusing on count modeling, cell-type annotation, data integration, including spatial transcriptomics, and the inference of cell-cell communication. We review the challenges encountered in scRNA-seq analysis, including issues of sparsity or low expression, reliability of cell annotation, and assumptions in data integration, and discuss the potential impact of suboptimal clustering and differential expression analysis tools on downstream analyses, particularly in identifying cell subpopulations. Finally, we discuss recent advancements and future directions for enhancing scRNA-seq analysis. Specifically, we highlight the development of novel tools for annotating single-cell data, integrating and interpreting multimodal datasets covering transcriptomics, epigenomics, and proteomics, and inferring cellular communication networks. By elucidating the latest progress and innovation, we provide a comprehensive overview of the rapidly advancing field of scRNA-seq analysis.

Cell-cell interactions and fluctuations in the direction of motility promote directed migration of osteoblasts in direct current electrotaxis

Under both physiological (development, regeneration) and pathological conditions (cancer metastasis), cells migrate while sensing environmental cues in the form of mechanical, chemical or electrical stimuli. In the case of bone tissue, osteoblast migration is essential in bone regeneration. Although it is known that osteoblasts respond to exogenous electric fields, the underlying mechanism of electrotactic collective movement of human osteoblasts is unclear. Here, we present a computational model that describes the osteoblast cell migration in a direct current electric field as the motion of a collection of active self-propelled particles and takes into account fluctuations in the direction of single-cell migration, finite-range cell-cell interactions, and the interaction of a cell with the external electric field. By comparing this model with in vitro experiments in which human primary osteoblasts are exposed to a direct current electric field of different field strengths, we show that cell-cell interactions and fluctuations in the migration direction promote anode-directed collective migration of osteoblasts.

Dynamically modeling the effective range of IL-2 dosage in the treatment of systemic lupus erythematosus

Systemic lupus erythematosus (SLE) is a complex systemic autoimmune disease characterized by an overactive immune response to self-antigen. The overactivation of CD4+ Foxp3- conventional T cells (Tcons) and the inactivation of CD4+ CD25+ Foxp3+ regulatory T cells (Tregs) play important roles in the progression of SLE. Clinical trials showed that low-dose interleukin-2 (IL-2) is effective in treating SLE. Here, we developed a mathematical model involving Tcons, Tregs, natural killer (NK) cells, and IL-2 to simulate the dynamic processes involved in the treatment of SLE. We found an effective range of IL-2 dosage defined by the Tcon/Treg ratio in SLE treatment, termed the IL-2 dosage therapeutic window (IDTW). Our results showed that high levels of self-antigen result in a narrow IDTW and high post-treatment Tcon/Treg ratio. Furthermore, we proposed a classification method based on the ratio of pre-treatment Treg to CD4+ T cells to predict the treatment outcome of SLE patients.

Dissecting the dynamic transcriptional landscape of early T helper cell differentiation into Th1, Th2, and Th1/2 hybrid cells

Selective differentiation of CD4+ T helper (Th) cells into specialized subsets such as Th1 and Th2 cells is a key element of the adaptive immune system driving appropriate immune responses. Besides those canonical Th-cell lineages, hybrid phenotypes such as Th1/2 cells arise in vivo, and their generation could be reproduced in vitro. While master-regulator transcription factors like T-bet for Th1 and GATA-3 for Th2 cells drive and maintain differentiation into the canonical lineages, the transcriptional architecture of hybrid phenotypes is less well understood. In particular, it has remained unclear whether a hybrid phenotype implies a mixture of the effects of several canonical lineages for each gene, or rather a bimodal behavior across genes. Th-cell differentiation is a dynamic process in which the regulatory factors are modulated over time, but longitudinal studies of Th-cell differentiation are sparse. Here, we present a dynamic transcriptome analysis following Th-cell differentiation into Th1, Th2, and Th1/2 hybrid cells at 3-h time intervals in the first hours after stimulation. We identified an early bifurcation point in gene expression programs, and we found that only a minority of ~20% of Th cell-specific genes showed mixed effects from both Th1 and Th2 cells on Th1/2 hybrid cells. While most genes followed either Th1- or Th2-cell gene expression, another fraction of ~20% of genes followed a Th1 and Th2 cell-independent transcriptional program associated with the transcription factors STAT1 and STAT4. Overall, our results emphasize the key role of high-resolution longitudinal data for the characterization of cellular phenotypes.

Modular, robust and extendible multicellular circuit design in yeast

Division of labor between cells is ubiquitous in biology but the use of multi-cellular consortia for engineering applications is only beginning to be explored. A significant advantage of multi-cellular circuits is their potential to be modular with respect to composition but this claim has not yet been extensively tested using experiments and quantitative modeling. Here, we construct a library of 24 yeast strains capable of sending, receiving or responding to three molecular signals, characterize them experimentally and build quantitative models of their input-output relationships. We then compose these strains into two- and three-strain cascades as well as a four-strain bistable switch and show that experimentally measured consortia dynamics can be predicted from the models of the constituent parts. To further explore the achievable range of behaviors, we perform a fully automated computational search over all two-, three- and four-strain consortia to identify combinations that realize target behaviors including logic gates, band-pass filters and time pulses. Strain combinations that are predicted to map onto a target behavior are further computationally optimized and then experimentally tested. Experiments closely track computational predictions. The high reliability of these model descriptions further strengthens the feasibility and highlights the potential for distributed computing in synthetic biology.

Protease-controlled secretion and display of intercellular signals

To program intercellular communication for biomedicine, it is crucial to regulate the secretion and surface display of signaling proteins. If such regulations are at the protein level, there are additional advantages, including compact delivery and direct interactions with endogenous signaling pathways. Here we create a modular, generalizable design called Retained Endoplasmic Cleavable Secretion (RELEASE), with engineered proteins retained in the endoplasmic reticulum and displayed/secreted in response to specific proteases. The design allows functional regulation of multiple synthetic and natural proteins by synthetic protease circuits to realize diverse signal processing capabilities, including logic operation and threshold tuning. By linking RELEASE to additional sensing and processing circuits, we can achieve elevated protein secretion in response to "undruggable" oncogene KRAS mutants. RELEASE should enable the local, programmable delivery of intercellular cues for a broad variety of fields such as neurobiology, cancer immunotherapy and cell transplantation.

A single-cell resolved cell-cell communication model explains lineage commitment in hematopoiesis

Cells do not make fate decisions independently. Arguably, every cell-fate decision occurs in response to environmental signals. In many cases, cell-cell communication alters the dynamics of the internal gene regulatory network of a cell to initiate cell-fate transitions, yet models rarely take this into account. Here, we have developed a multiscale perspective to study the granulocyte-monocyte versus megakaryocyte-erythrocyte fate decisions. This transition is dictated by the GATA1-PU.1 network: a classical example of a bistable cell-fate system. We show that, for a wide range of cell communication topologies, even subtle changes in signaling can have pronounced effects on cell-fate decisions. We go on to show how cell-cell coupling through signaling can spontaneously break the symmetry of a homogenous cell population. Noise, both intrinsic and extrinsic, shapes the decision landscape profoundly, and affects the transcriptional dynamics underlying this important hematopoietic cell-fate decision-making system.

This article has an associated ‘The people behind the papers’ interview.

Summary: Through theory and computational modeling, cell-cell communication is revealed to be a crucial and under-appreciated determinant of cell-fate decision-making during hematopoiesis.

Cell-Cell Communication Networks in Tissue: Toward Quantitatively Linking Structure with Function

Forefront techniques for molecular interrogation of mammalian tissues, such as multiplexed tissue imaging, intravital microscopy, and single-cell RNA sequencing (scRNAseq), can combine to quantify cell-type abundance, co-localization, and global levels of receptors and their ligands. Nonetheless, it remains challenging to translate these various quantities into a more comprehensive understanding of how cell-cell communication networks dynamically operate. Therefore, construction of computational models for network-level functions - including niche-dependent actions, homeostasis, and multi-scale coordination - will be valuable for productively integrating the battery of experimental approaches. Here, we review recent progress in understanding cell-cell communication networks in tissue. Featured examples include ligand-receptor dissection of immunosuppressive and mitogenic signaling in the tumor microenvironment. As a future direction, we highlight an unmet potential to bridge high-level statistical approaches with low-level physicochemical mechanisms.

Modeling the Th17 and Tregs Paradigm: Implications for Cancer Immunotherapy

CD4 + T cell differentiation is governed by gene regulatory and metabolic networks, with both networks being highly interconnected and able to adapt to external stimuli. Th17 and Tregs differentiation networks play a critical role in cancer, and their balance is affected by the tumor microenvironment (TME). Factors from the TME mediate recruitment and expansion of Th17 cells, but these cells can act with pro or anti-tumor immunity. Tregs cells are also involved in tumor development and progression by inhibiting antitumor immunity and promoting immunoevasion. Due to the complexity of the underlying molecular pathways, the modeling of biological systems has emerged as a promising solution for better understanding both CD4 + T cell differentiation and cancer cell behavior. In this review, we present a context-dependent vision of CD4 + T cell transcriptomic and metabolic network adaptability. We then discuss CD4 + T cell knowledge-based models to extract the regulatory elements of Th17 and Tregs differentiation in multiple CD4 + T cell levels. We highlight the importance of complementing these models with data from omics technologies such as transcriptomics and metabolomics, in order to better delineate existing Th17 and Tregs bifurcation mechanisms. We were able to recompilate promising regulatory components and mechanisms of Th17 and Tregs differentiation under normal conditions, which we then connected with biological evidence in the context of the TME to better understand CD4 + T cell behavior in cancer. From the integration of mechanistic models with omics data, the transcriptomic and metabolomic reprograming of Th17 and Tregs cells can be predicted in new models with potential clinical applications, with special relevance to cancer immunotherapy.

Memorizing environmental signals through feedback and feedforward loops

Cells in diverse organisms can store the information of previous environmental conditions for long periods of time. This form of cellular memory adjusts the cell's responses to future challenges, providing fitness advantages in fluctuating environments. Many biological functions, including cellular memory, are mediated by specific recurring patterns of interactions among proteins and genes, known as 'network motifs.' In this review, we focus on three well-characterized network motifs - negative feedback loops, positive feedback loops, and feedforward loops, which underlie different types of cellular memories. We describe the latest studies identifying these motifs in various molecular processes and discuss how the topologies and dynamics of these motifs can enable memory encoding and storage.

Reconstruction of Intercellular Signaling Network by Cytokine-Receptor Interactions

Background:

The immune system function depends on the coordination activity of the components of system and communications between them which leads to the formation of a complex communication network between immune cells. In this network, cytokines have an important role in the communication between immune cells through the interaction to their specific receptors. These molecules cause to cellular communications and normal function of a tissue. Reconstruction of such a complex network can be a way to provide a better understanding of cytokines’ function.

Objective:

Our main goal from reconstructing such a network was investigation of expressed cytokines and cytokines receptors in various lineage and tissues of immune cells and identifying the lineage and tissue with the highest expression of cytokines and their receptors.

Materials and Methods:

In this study, gene expression data related to part of the Immunological Genome Project (ImmGen) and receptor-ligand interactions dataset were used to reconstruct the immune network in mouse. In next step, the topological properties of reconstructed network, expression specificity of cytokines and their receptors and interactions specificity were analyzed.

Results:

The results of the network analysis were indicated that non- hematopoietic stromal cells have the highest expression of cytokines and cytokine receptors and interactions specificity is very high. Our results show that chemokine receptor of Ccr1 receives the largest number of signals between receptors and only expressed in three hematopoietic lineages.

Conclusions:

The most of the network communications belonged to non-hematopoietic stromal and macrophage cells. The relationships between stromal cells and macrophages are necessary to create an appropriate environment for differentiation of immune cells. Studying the cellular expression specificity of receptor and ligand genes reveal the high degree of specificity of these genes that indicate non-random transfer of information between cells in multicellular organisms.

Phenotypic Heterogeneity in Tumor Progression, and Its Possible Role in the Onset of Cancer

Heterogeneity among isogenic cells/individuals has been known for at least 150 years. Even Mendel, working on pea plants, realized that not all tall plants were identical. However, Mendel was more interested in the discontinuous variation between genetically distinct individuals. The concept of environment dictating distinct phenotypes among isogenic individuals has since been shown to impact the evolution of populations in numerous examples at different scales of life. In this review, we discuss how phenotypic heterogeneity and its evolutionary implications exist at all levels of life, from viruses to mammals. In particular, we discuss how a particular disease condition (cancer) is impacted by heterogeneity among isogenic cells, and propose a potential role that phenotypic heterogeneity might play toward the onset of the disease.

Multiple random phosphorylations in clock proteins provide long delays and switches

Theory predicts that self-sustained oscillations require robust delays and nonlinearities (ultrasensitivity). Delayed negative feedback loops with switch-like inhibition of transcription constitute the core of eukaryotic circadian clocks. The kinetics of core clock proteins such as PER2 in mammals and FRQ in Neurospora crassa is governed by multiple phosphorylations. We investigate how multiple, slow and random phosphorylations control delay and molecular switches. We model phosphorylations of intrinsically disordered clock proteins (IDPs) using conceptual models of sequential and distributive phosphorylations. Our models help to understand the underlying mechanisms leading to delays and ultrasensitivity. The model shows temporal and steady state switches for the free kinase and the phosphoprotein. We show that random phosphorylations and sequestration mechanisms allow high Hill coefficients required for self-sustained oscillations.

Cell-cycle-gated feedback control mediates desensitization to interferon stimulation

Cells use molecular circuits to interpret and respond to extracellular cues, such as hormones and cytokines, which are often released in a temporally varying fashion. In this study, we combine microfluidics, time-lapse microscopy, and computational modeling to investigate how the type I interferon (IFN)-responsive regulatory network operates in single human cells to process repetitive IFN stimulation. We found that IFN-α pretreatments lead to opposite effects, priming versus desensitization, depending on input durations. These effects are governed by a regulatory network composed of a fast-acting positive feedback loop and a delayed negative feedback loop, mediated by upregulation of ubiquitin-specific peptidase 18 (USP18). We further revealed that USP18 upregulation can only be initiated at the G1/early S phases of cell cycle upon the treatment onset, resulting in heterogeneous and delayed induction kinetics in single cells. This cell cycle gating provides a temporal compartmentalization of feedback loops, enabling duration-dependent desensitization to repetitive stimulations.

Toward Systems Biomarkers of Response to Immune Checkpoint Blockers

Immunotherapy with checkpoint blockers (ICBs), aimed at unleashing the immune response toward tumor cells, has shown a great improvement in overall patient survival compared to standard therapy, but only in a subset of patients. While a number of recent studies have significantly improved our understanding of mechanisms playing an important role in the tumor microenvironment (TME), we still have an incomplete view of how the TME works as a whole. This hampers our ability to effectively predict the large heterogeneity of patients' response to ICBs. Systems approaches could overcome this limitation by adopting a holistic perspective to analyze the complexity of tumors. In this Mini Review, we focus on how an integrative view of the increasingly available multi-omics experimental data and computational approaches enables the definition of new systems-based predictive biomarkers. In particular, we will focus on three facets of the TME toward the definition of new systems biomarkers. First, we will review how different types of immune cells influence the efficacy of ICBs, not only in terms of their quantification, but also considering their localization and functional state. Second, we will focus on how different cells in the TME interact, analyzing how inter- and intra-cellular networks play an important role in shaping the immune response and are responsible for resistance to immunotherapy. Finally, we will describe the potential of looking at these networks as dynamic systems and how mathematical models can be used to study the rewiring of the complex interactions taking place in the TME.

Molecular-memory-driven phenotypic switching in a genetic toggle switch without cooperative binding

A genetic toggle switch would involve multistep reaction processes (e.g., complex promoter activation), creating memories between individual reaction events. Revealing the effect of this molecular memory is important for understanding intracellular processes such as cellular decision making. We propose a generalized genetic toggle switch model and use a generalized chemical master equation theory to account for the memory effect. Interestingly, we find that molecular memory can induce bimodality in this memory system although the corresponding memoryless counterpart is not bimodal. This finding implies a plausible alternative mechanism for phenotypic switching that is driven by molecular memory rather than by ultrasensitivity or cooperative binding as shown in previous studies. We also find that unbalanced memories arising from the processes by which mutually inhibiting transcription factors are produced can give rise to asymmetric bimodality without changing the positions of two peaks in the bimodal protein distribution. Given the prevalence of molecular memory in gene regulation, our findings would provide insights into cell fate decisions in growth and development.

Agent-based model of multicellular tumor spheroid evolution including cell metabolism

Computational models aiming at the spatio-temporal description of cancer evolution are a suitable framework for testing biological hypotheses from experimental data, and generating new ones. Building on our recent work (J. Theor. Biol. 389, 146 (2016)) we develop a 3D agent-based model, capable of tracking hundreds of thousands of interacting cells, over time scales ranging from seconds to years. Cell dynamics is driven by a Monte Carlo solver, incorporating partial differential equations to describe chemical pathways and the activation/repression of "genes", leading to the up- or down-regulation of specific cell markers. Each cell-agent of different kind (stem, cancer, stromal etc.) runs through its cycle, undergoes division, can exit to a dormant, senescent, necrotic state, or apoptosis, according to the inputs from its systemic network. The basic network at this stage describes glucose/oxygen/ATP cycling, and can be readily extended to cancer-cell specific markers. Eventual accumulation of chemical/radiation damage to each cell's DNA is described by a Markov chain of internal states, and by a damage-repair network, whose evolution is linked to the cell systemic network. Aimed at a direct comparison with experiments of tumorsphere growth from stem cells, the present model will allow to quantitatively study the role of transcription factors involved in the reprogramming and variable radio-resistance of simulated cancer-stem cells, evolving in a realistic computer simulation of a growing multicellular tumorsphere.

Network modeling of single-cell omics data: challenges, opportunities, and progresses

Single-cell multi-omics technologies are rapidly evolving, prompting both methodological advances and biological discoveries at an unprecedented speed. Gene regulatory network modeling has been used as a powerful approach to elucidate the complex molecular interactions underlying biological processes and systems, yet its application in single-cell omics data modeling has been met with unique challenges and opportunities. In this review, we discuss these challenges and opportunities, and offer an overview of the recent development of network modeling approaches designed to capture dynamic networks, within-cell networks, and cell-cell interaction or communication networks. Finally, we outline the remaining gaps in single-cell gene network modeling and the outlooks of the field moving forward.

Transmission of Mechanical Information by Purinergic Signaling

The skeleton constantly interacts and adapts to the physical world. We have previously reported that physiologically relevant mechanical forces lead to small repairable membrane injuries in bone-forming osteoblasts, resulting in release of ATP and stimulation of purinergic (P2) calcium responses in neighboring cells. The goal of this study was to develop a theoretical model describing injury-related ATP and ADP release, their extracellular diffusion and degradation, and purinergic responses in neighboring cells. After validation using experimental data for intracellular free calcium elevations, ATP, and vesicular release after mechanical stimulation of a single osteoblast, the model was scaled to a tissue-level injury to investigate how purinergic signaling communicates information about injuries with varying geometries. We found that total ATP released, peak extracellular ATP concentration, and the ADP-mediated signaling component contributed complementary information regarding the mechanical stimulation event. The total amount of ATP released governed spatial factors, such as the maximal distance from the injury at which purinergic responses were stimulated. The peak ATP concentration reflected the severity of an individual cell injury, allowing to discriminate between minor and severe injuries that released similar amounts of ATP because of differences in injury repair, and determined temporal aspects of the response, such as signal propagation velocity. ADP-mediated signaling became relevant only in larger tissue-level injuries, conveying information about the distance to the injury site and its geometry. Thus, we identified specific features of extracellular ATP and ADP spatiotemporal signals that depend on tissue mechanoresilience and encode the severity, scope, and proximity of the mechanical stimulus.

Synthetic Switches and Regulatory Circuits in Plants

Synthetic biology is an established but ever-growing interdisciplinary field of research currently revolutionizing biomedicine studies and the biotech industry. The engineering of synthetic circuitry in bacterial, yeast, and animal systems prompted considerable advances for the understanding and manipulation of genetic and metabolic networks; however, their implementation in the plant field lags behind. Here, we review theoretical-experimental approaches to the engineering of synthetic chemical- and light-regulated (optogenetic) switches for the targeted interrogation and control of cellular processes, including existing applications in the plant field. We highlight the strategies for the modular assembly of genetic parts into synthetic circuits of different complexity, ranging from Boolean logic gates and oscillatory devices up to semi- and fully synthetic open- and closed-loop molecular and cellular circuits. Finally, we explore potential applications of these approaches for the engineering of novel functionalities in plants, including understanding complex signaling networks, improving crop productivity, and the production of biopharmaceuticals.

Computational Structural Biology: Successes, Future Directions, and Challenges

Computational biology has made powerful advances. Among these, trends in human health have been uncovered through heterogeneous 'big data' integration, and disease-associated genes were identified and classified. Along a different front, the dynamic organization of chromatin is being elucidated to gain insight into the fundamental question of genome regulation. Powerful conformational sampling methods have also been developed to yield a detailed molecular view of cellular processes. when combining these methods with the advancements in the modeling of supramolecular assemblies, including those at the membrane, we are finally able to get a glimpse into how cells' actions are regulated. Perhaps most intriguingly, a major thrust is on to decipher the mystery of how the brain is coded. Here, we aim to provide a broad, yet concise, sketch of modern aspects of computational biology, with a special focus on computational structural biology. We attempt to forecast the areas that computational structural biology will embrace in the future and the challenges that it may face. We skirt details, highlight successes, note failures, and map directions.

Multi-Omics Profiling of the Tumor Microenvironment: Paving the Way to Precision Immuno-Oncology

The tumor microenvironment (TME) is a multifaceted ecosystem characterized by profound cellular heterogeneity, dynamicity, and complex intercellular cross-talk. The striking responses obtained with immune checkpoint blockers, i.e., antibodies targeting immune-cell regulators to boost antitumor immunity, have demonstrated the enormous potential of anticancer treatments that target TME components other than tumor cells. However, as checkpoint blockade is currently beneficial only to a limited fraction of patients, there is an urgent need to understand the mechanisms orchestrating the immune response in the TME to guide the rational design of more effective anticancer therapies. In this Mini Review, we give an overview of the methodologies that allow studying the heterogeneity of the TME from multi-omics data generated from bulk samples, single cells, or images of tumor-tissue slides. These include approaches for the characterization of the different cell phenotypes and for the reconstruction of their spatial organization and inter-cellular cross-talk. We discuss how this broader vision of the cellular heterogeneity and plasticity of tumors, which is emerging thanks to these methodologies, offers the opportunity to rationally design precision immuno-oncology treatments. These developments are fundamental to overcome the current limitations of targeted agents and checkpoint blockers and to bring long-term clinical benefits to a larger fraction of cancer patients.

Sharing and reusing cell image data

The rapid growth in content and complexity of cell image data creates an opportunity for synergy between experimental and computational scientists. Sharing microscopy data enables computational scientists to develop algorithms and tools for data analysis, integration, and mining. These tools can be applied by experimentalists to promote hypothesis-generation and discovery. We are now at the dawn of this revolution: infrastructure is being developed for data standardization, deposition, sharing, and analysis; some journals and funding agencies mandate data deposition; data journals publish high-content microscopy data sets; quantification becomes standard in scientific publications; new analytic tools are being developed and dispatched to the community; and huge data sets are being generated by individual labs and philanthropic initiatives. In this Perspective, I reflect on sharing and reusing cell image data and the opportunities that will come along with it.

The stretch to stray on time: Resonant length of random walks in a transient

First-passage times in random walks have a vast number of diverse applications in physics, chemistry, biology, and finance. In general, environmental conditions for a stochastic process are not constant on the time scale of the average first-passage time or control might be applied to reduce noise. We investigate moments of the first-passage time distribution under an exponential transient describing relaxation of environmental conditions. We solve the Laplace-transformed (generalized) master equation analytically using a novel method that is applicable to general state schemes. The first-passage time from one end to the other of a linear chain of states is our application for the solutions. The dependence of its average on the relaxation rate obeys a power law for slow transients. The exponent ν depends on the chain length N like ν=-N/(N+1) to leading order. Slow transients substantially reduce the noise of first-passage times expressed as the coefficient of variation (CV), even if the average first-passage time is much longer than the transient. The CV has a pronounced minimum for some lengths, which we call resonant lengths. These results also suggest a simple and efficient noise control strategy and are closely related to the timing of repetitive excitations, coherence resonance, and information transmission by noisy excitable systems. A resonant number of steps from the inhibited state to the excitation threshold and slow recovery from negative feedback provide optimal timing noise reduction and information transmission.