From biological data to oscillator models using SINDy.

Periodic changes in the concentration or activity of different molecules regulate vital cellular processes such as cell division and circadian rhythms. Developing mathematical models is essential to better understand the mechanisms underlying these oscillations. Recent data-driven methods like SINDy have fundamentally changed model identification, yet their application to experimental biological data remains limited. This study investigates SINDy's constraints by directly applying it to biological oscillatory data. We identify insufficient resolution, noise, dimensionality, and limited prior knowledge as primary limitations. Using various generic oscillator models of different complexity and/or dimensionality, we systematically analyze these factors. We then propose a comprehensive guide for inferring models from biological data, addressing these challenges step by step. Our approach is validated using glycolytic oscillation data from yeast.

Mitotic waves in frog egg extracts: Transition from phase waves to trigger waves.

Cyclin-dependent kinase 1 (Cdk1) activity rises and falls throughout the cell cycle, a cell-autonomous process known as mitotic oscillations. These oscillators can synchronize when spatially coupled, providing a crucial foundation for rapid synchronous divisions in large early embryos like Drosophila (~ 0.5 mm) and Xenopus (~ 1.2 mm). While diffusion alone cannot achieve such long-range coordination, recent studies have proposed two types of mitotic waves, phase and trigger waves, to explain the phenomena. How the waves establish over time for efficient spatial coordination remains unclear. Using Xenopus laevis egg extracts and a Cdk1 FRET sensor, we observe a transition from phase waves to a trigger wave regime in an initially homogeneous cytosol. Adding nuclei accelerates such transition. Moreover, the system transitions almost immediately to this regime when externally driven by metaphase-arrested extracts from the boundary. Employing computational modeling, we pinpoint how wave nature, including speed-period relation, depends on transient dynamics and oscillator properties, suggesting that phase waves appear transiently due to the time required for trigger waves to entrain the system and that spatial heterogeneity promotes entrainment. Therefore, we show that both waves belong to a single biological process capable of coordinating the cell cycle over long distances.

Challenges in identifying simple pattern-forming mechanisms in the development of settlements using demographic data.

The rapid increase of population and settlement structures in the Global South during recent decades has motivated the development of suitable models to describe their formation and evolution. Such settlement formation has been previously suggested to be dynamically driven by simple pattern-forming mechanisms. Here, we explore the use of a data-driven white-box approach, called SINDy, to discover differential equation models directly from available spatiotemporal demographic data for three representative regions of the Global South. We show that the current resolution and observation time of the available data are insufficient to uncover relevant pattern-forming mechanisms in settlement development. Using synthetic data generated with a generic pattern-forming model, the Allen-Cahn equation, we characterize what the requirements are for spatial and temporal resolution, as well as observation time, to successfully identify possible model system equations. Overall, the study provides a theoretical framework for the analysis of large-scale geographical and/or ecological systems, and it motivates further improvements in optimization approaches and data collection.

The ups and downs of biological oscillators: a comparison of time-delayed negative feedback mechanisms.

Many biochemical oscillators are driven by the periodic rise and fall of protein concentrations or activities. A negative feedback loop underlies such oscillations. The feedback can act on different parts of the biochemical network. Here, we mathematically compare time-delay models where the feedback affects production and degradation. We show a mathematical connection between the linear stability of the two models, and derive how both mechanisms impose different constraints on the production and degradation rates that allow oscillations. We show how oscillations are affected by the inclusion of a distributed delay, of double regulation (acting on production and degradation) and of enzymatic degradation.

Multiomics and spatial mapping characterizes human CD8+ T cell states in cancer.

Clinically relevant immunological biomarkers that discriminate between diverse hypofunctional states of tumor-associated CD8+ T cells remain disputed. Using multiomics analysis of CD8+ T cell features across multiple patient cohorts and tumor types, we identified tumor niche-dependent exhausted and other types of hypofunctional CD8+ T cell states. CD8+ T cells in "supportive" niches, like melanoma or lung cancer, exhibited features of tumor reactivity-driven exhaustion (CD8+ TEX). These included a proficient effector memory phenotype, an expanded T cell receptor (TCR) repertoire linked to effector exhaustion signaling, and a cancer-relevant T cell-activating immunopeptidome composed of largely shared cancer antigens or neoantigens. In contrast, "nonsupportive" niches, like glioblastoma, were enriched for features of hypofunctionality distinct from canonical exhaustion. This included immature or insufficiently activated T cell states, high wound healing signatures, nonexpanded TCR repertoires linked to anti-inflammatory signaling, high T cell-recognizable self-epitopes, and an antiproliferative state linked to stress or prodeath responses. In situ spatial mapping of glioblastoma highlighted the prevalence of dysfunctional CD4+:CD8+ T cell interactions, whereas ex vivo single-cell secretome mapping of glioblastoma CD8+ T cells confirmed negligible effector functionality and a promyeloid, wound healing-like chemokine profile. Within immuno-oncology clinical trials, anti-programmed cell death protein 1 (PD-1) immunotherapy facilitated glioblastoma's tolerogenic disparities, whereas dendritic cell (DC) vaccines partly corrected them. Accordingly, recipients of a DC vaccine for glioblastoma had high effector memory CD8+ T cells and evidence of antigen-specific immunity. Collectively, we provide an atlas for assessing different CD8+ T cell hypofunctional states in immunogenic versus nonimmunogenic cancers.

Cell cycle oscillations driven by two interlinked bistable switches.

Regular transitions between interphase and mitosis during the cell cycle are driven by changes in the activity of the enzymatic protein complex cyclin B with cyclin-dependent kinase 1 (Cdk1). At the most basic level, this cell cycle oscillator is driven by negative feedback: active cyclin B-Cdk1 activates the anaphase-promoting complex/cyclosome, which triggers the degradation of cyclin B. Such cell cycle oscillations occur fast and periodically in the early embryos of the frog Xenopus laevis, where several positive-feedback loops leading to bistable switches in parts of the regulatory network have been experimentally identified. Here, we build cell cycle oscillator models to show how single and multiple bistable switches in parts of the underlying regulatory network change the properties of the oscillations and how they can confer robustness to the oscillator. We present a detailed bifurcation analysis of these models.

Analytical approximations for the speed of pacemaker-generated waves.

In oscillatory media, waves can be generated by pacemaker regions which oscillate faster than their surroundings. In many chemical and biological systems, such waves can synchronize the whole medium and as such they are a means of transmitting information at a fixed speed over large distances. In this paper, we apply analytical tools to investigate the factors that determine the speed of these waves. More precisely, we apply singular perturbation and phase reduction methods to two types of negative-feedback oscillators, one built on underlying bistability and one including a time delay in the negative feedback. In both systems, we investigate the influence of time-scale separation on the resulting wave speed, as well as the effect of size and frequency of the pacemaker region. We compare our analytical estimates to numerical simulations which we described previously [J. Rombouts and L. Gelens, Phys. Rev. Research 2, 043038 (2020)2643-156410.1103/PhysRevResearch.2.043038].

A modular approach for modeling the cell cycle based on functional response curves.

Modeling biochemical reactions by means of differential equations often results in systems with a large number of variables and parameters. As this might complicate the interpretation and generalization of the obtained results, it is often desirable to reduce the complexity of the model. One way to accomplish this is by replacing the detailed reaction mechanisms of certain modules in the model by a mathematical expression that qualitatively describes the dynamical behavior of these modules. Such an approach has been widely adopted for ultrasensitive responses, for which underlying reaction mechanisms are often replaced by a single Hill function. Also time delays are usually accounted for by using an explicit delay in delay differential equations. In contrast, however, S-shaped response curves, which by definition have multiple output values for certain input values and are often encountered in bistable systems, are not easily modeled in such an explicit way. Here, we extend the classical Hill function into a mathematical expression that can be used to describe both ultrasensitive and S-shaped responses. We show how three ubiquitous modules (ultrasensitive responses, S-shaped responses and time delays) can be combined in different configurations and explore the dynamics of these systems. As an example, we apply our strategy to set up a model of the cell cycle consisting of multiple bistable switches, which can incorporate events such as DNA damage and coupling to the circadian clock in a phenomenological way.

Dynamic bistable switches enhance robustness and accuracy of cell cycle transitions.

Bistability is a common mechanism to ensure robust and irreversible cell cycle transitions. Whenever biological parameters or external conditions change such that a threshold is crossed, the system abruptly switches between different cell cycle states. Experimental studies have uncovered mechanisms that can make the shape of the bistable response curve change dynamically in time. Here, we show how such a dynamically changing bistable switch can provide a cell with better control over the timing of cell cycle transitions. Moreover, cell cycle oscillations built on bistable switches are more robust when the bistability is modulated in time. Our results are not specific to cell cycle models and may apply to other bistable systems in which the bistable response curve is time-dependent.

Bistable, Biphasic Regulation of PP2A-B55 Accounts for the Dynamics of Mitotic Substrate Phosphorylation.

The phosphorylation of mitotic proteins is bistable, which contributes to the decisiveness of the transitions into and out of M phase. The bistability in substrate phosphorylation has been attributed to bistability in the activation of the cyclin-dependent kinase Cdk1. However, more recently it has been suggested that bistability also arises from positive feedback in the regulation of the Cdk1-counteracting phosphatase PP2A-B55. Here, we demonstrate biochemically using Xenopus laevis egg extracts that the Cdk1-counteracting phosphatase PP2A-B55 functions as a bistable switch, even when the bistability of Cdk1 activation is suppressed. In addition, Cdk1 regulates PP2A-B55 in a biphasic manner; low concentrations of Cdk1 activate PP2A-B55 and high concentrations inactivate it. As a consequence of this incoherent feedforward regulation, PP2A-B55 activity rises concurrently with Cdk1 activity during interphase and suppresses substrate phosphorylation. PP2A-B55 activity is then sharply downregulated at the onset of mitosis. During mitotic exit, Cdk1 activity initially falls with no obvious change in substrate phosphorylation; dephosphorylation then commences once PP2A-B55 spikes in activity. These findings suggest that changes in Cdk1 activity are permissive for mitotic entry and exit but that the changes in PP2A-B55 activity are the ultimate trigger.

Nuclei determine the spatial origin of mitotic waves.

Traveling waves play an essential role in coordinating mitosis over large distances, but what determines the spatial origin of mitotic waves remains unclear. Here, we show that such waves initiate at pacemakers, regions that oscillate faster than their surroundings. In cell-free extracts of Xenopus laevis eggs, we find that nuclei define such pacemakers by concentrating cell cycle regulators. In computational models of diffusively coupled oscillators that account for nuclear import, nuclear positioning determines the pacemaker location. Furthermore, we find that the spatial dimensions of the oscillatory medium change the nuclear positioning and strongly influence whether a pacemaker is more likely to be at a boundary or an internal region. Finally, we confirm experimentally that increasing the system width increases the proportion of pacemakers at the boundary. Our work provides insight into how nuclei and spatial system dimensions can control local concentrations of regulators and influence the emergent behavior of mitotic waves.

Mutualistic cross-feeding in microbial systems generates bistability via an Allee effect.

In microbial ecosystems, species not only compete for common resources but may also display mutualistic interactions as a result from metabolic cross-feeding. Such mutualism can lead to bistability. Depending on the initial population sizes, species will either survive or go extinct. Various phenomenological models have been suggested to describe bistability in mutualistic systems. However, these models do not account for interaction mediators such as nutrients. In contrast, nutrient-explicit models do not provide an intuitive understanding of what causes bistability. Here, we reduce a theoretical nutrient-explicit model of two mutualistic cross-feeders in a chemostat, uncovering an explicit relation to a growth model with an Allee effect. We show that the dilution rate in the chemostat leads to bistability by turning a weak Allee effect into a strong Allee effect. This happens as long as there is more production than consumption of cross-fed nutrients. Thanks to the explicit relationship of the reduced model with the underlying experimental parameters, these results allow to predict the biological conditions that sustain or prevent the survival of mutualistic species.

Co-regulation of the antagonistic RepoMan:Aurora-B pair in proliferating cells.

Chromosome segregation during mitosis is antagonistically regulated by the Aurora-B kinase and RepoMan (recruits PP1 onto mitotic chromatin at anaphase)-associated phosphatases PP1/PP2A. Aurora B is overexpressed in many cancers but, surprisingly, this only rarely causes lethal aneuploidy. Here we show that RepoMan abundance is regulated by the same mechanisms that control Aurora B, including FOXM1-regulated expression and proteasomal degradation following ubiquitination by APC/C-CDH1 or SCFFBXW7. The deregulation of these mechanisms can account for the balanced co-overexpression of Aurora B and RepoMan in many cancers, which limits chromosome segregation errors. In addition, Aurora B and RepoMan independently promote cancer cell proliferation by reducing checkpoint--induced cell-cycle arrest during interphase. The co-up-regulation of RepoMan and Aurora B in tumors is inversely correlated with patient survival, underscoring its potential importance for tumor progression. Finally, we demonstrate that high RepoMan levels sensitize cancer cells to Aurora-B inhibitors. Hence, the co-up-regulation of RepoMan and Aurora B is associated with tumor aggressiveness but also exposes a vulnerable target for therapeutic intervention.

Travelling fronts in time-delayed reaction-diffusion systems.

We review a series of key travelling front problems in reaction-diffusion systems with a time-delayed feedback, appearing in ecology, nonlinear optics and neurobiology. For each problem, we determine asymptotic approximations for the wave shape and its speed. Particular attention is devoted to their validity and all analytical solutions are compared to solutions obtained numerically. We also extend the work by Erneux et al. (Erneux et al. 2010 Phil. Trans. R. Soc. A 368, 483-493 (doi:10.1098/rsta.2009.0228)) by considering the case of a slowly propagating front subject to a weak delayed feedback. The delay may either speed up the front in the same direction or reverse its direction. This article is part of the theme issue 'Nonlinear dynamics of delay systems'.

Frequency comb generation through the locking of domain walls in doubly resonant dispersive optical parametric oscillators.

In this Letter we theoretically investigate the formation of localized temporal dissipative structures, and their corresponding frequency combs in doubly resonant dispersive optical parametric oscillators. We derive a nonlocal mean field model, and show that domain wall locking allows for the formation of stable coherent optical frequency combs.

Excitable dynamics through toxin-induced mRNA cleavage in bacteria.

Toxin-antitoxin (TA) systems in bacteria and archaea are small genetic elements consisting of the genes coding for an intracellular toxin and an antitoxin that can neutralize this toxin. In various cases, the toxins cleave the mRNA. In this theoretical work we use deterministic and stochastic modeling to explain how toxin-induced cleavage of mRNA in TA systems can lead to excitability, allowing large transient spikes in toxin levels to be triggered. By using a simplified network where secondary complex formation and transcriptional regulation are not included, we show that a two-dimensional, deterministic model captures the origin of such toxin excitations. Moreover, it allows to increase our understanding by examining the dynamics in the phase plane. By systematically comparing the deterministic results with Gillespie simulations we demonstrate that even though the real TA system is intrinsically stochastic, toxin excitations can be accurately described deterministically. A bifurcation analysis of the system shows that the excitable behavior is due to a nearby Hopf bifurcation in the parameter space, where the system becomes oscillatory. The influence of stress is modeled by varying the degradation rate of the antitoxin and the translation rate of the toxin. We find that stress increases the frequency of toxin excitations. The inclusion of secondary complex formation and transcriptional regulation does not fundamentally change the mechanism of toxin excitations. Finally, we show that including growth rate suppression and translational inhibition can lead to longer excitations, and even cause excitations in cases when the system would otherwise be non-excitable. To conclude, the deterministic model used in this work provides a simple and intuitive explanation of toxin excitations in TA systems.

Coordination of Timers and Sensors in Cell Signaling.

Timers and sensors are common devices that make our daily life safer, more convenient, and more efficient. In a cellular context, they arguably play an even more crucial role as they ensure the survival of cells in the presence of various extrinsic and intrinsic stresses. Biological timers and sensors generate distinct signaling profiles, enabling them to produce different types of cellular responses. Recent data suggest that they can work together to guarantee correct timing and responsiveness. By exploring examples of cellular stress signaling from mitosis, DNA damage, and hypoxia, the authors discuss the common architecture of timer-sensor integration, and how its added features contribute to the generation of desired signaling profiles when dealing with stresses of variable duration and strength. The authors propose timer-sensor integration as a widespread mechanism with profound biological implications and therapeutic potential.

Quadratic soliton combs in doubly resonant second-harmonic generation.

We report a theoretical investigation of quadratic frequency combs in a dispersive second-harmonic generation cavity system. We identify different dynamical regimes and demonstrate that the same system can exhibit both bright and dark localized cavity solitons in the absence of a temporal walk-off.