An official website of the United States government
The synchronization dynamics for the circadian gene expression in the suprachiasmatic nucleus is investigated using a transcriptional circadian clock gene oscillator model. With global coupling in constant dark (DD) conditions, the model exhibits a one-cluster phase synchronized state, in dim light (dim LL), bistability between one- and two-cluster states and in bright LL, a two-cluster state. The two-cluster phase synchronized state, where some oscillator pairs synchronize in-phase, and some anti-phase, can explain the splitting of the circadian clock, i.e., generation of two bouts of daily activities with certain species, e.g., with hamsters. The one- and two-cluster states can be reached by transferring the animal from DD or bright LL to dim LL, i.e., the circadian synchrony has a memory effect. The stability of the one- and two-cluster states was interpreted analytically by extracting phase models from the ordinary differential equation models. In a modular network with two strongly coupled oscillator populations with weak intragroup coupling, with appropriate initial conditions, one group is synchronized to the one-cluster state and the other group to the two-cluster state, resulting in a weak-chimera state. Computational modeling suggests that the daily rhythms in sleep–wake depend on light intensity acting on bilateral networks of suprachiasmatic nucleus (SCN) oscillators. Addition of a network heterogeneity (coupling between the left and right SCN) allowed the system to exhibit chimera states. The simulations can guide experiments in the circadian rhythm research to explore the effect of light intensity on the complexities of circadian desynchronization.
more »
« less
Non-monotonic transients to synchrony in Kuramoto networks and electrochemical oscillators
Abstract We performed numerical simulations with the Kuramoto model and experiments with oscillatory nickel electrodissolution to explore the dynamical features of the transients from random initial conditions to a fully synchronized (one-cluster) state. The numerical simulations revealed that certain networks (e.g., globally coupled or dense Erdős–Rényi random networks) showed relatively simple behavior with monotonic increase of the Kuramoto order parameter from the random initial condition to the fully synchronized state and that the transient times exhibited a unimodal distribution. However, some modular networks with bridge elements were identified which exhibited non-monotonic variation of the order parameter with local maximum and/or minimum. In these networks, the histogram of the transients times became bimodal and the mean transient time scaled well with inverse of the magnitude of the second largest eigenvalue of the network Laplacian matrix. The non-monotonic transients increase the relative standard deviations from about 0.3 to 0.5, i.e., the transient times became more diverse. The non-monotonic transients are related to generation of phase patterns where the modules are synchronized but approximately anti-phase to each other. The predictions of the numerical simulations were demonstrated in a population of coupled oscillatory electrochemical reactions in global, modular, and irregular tree networks. The findings clarify the role of network structure in generation of complex transients that can, for example, play a role in intermittent desynchronization of the circadian clock due to external cues or in deep brain stimulations where long transients are required after a desynchronization stimulus.
more »
« less
- Award ID(s):
- 1900011
- PAR ID:
- 10361405
- Publisher / Repository:
- IOP Publishing
- Date Published:
- Journal Name:
- Journal of Physics: Complexity
- Volume:
- 2
- Issue:
- 1
- ISSN:
- 2632-072X
- Page Range / eLocation ID:
- Article No. 015010
- Format(s):
- Medium: X
- Sponsoring Org:
- National Science Foundation
More Like this
-
-
To maintain normal functionality, it is necessary for a multicellular organism to generate robust responses to external temporal signals. However, the underlying mechanisms to coordinate the collective dynamics of cells remain poorly understood. Here, we study the calcium activity of biological neuron networks excited by periodic ATP stimuli. We use micropatterning to control the cells' physical connectivity. We find that whereas isolated cells become more synchronized in their calcium activity at long driving periods, connected cells become less synchronized, despite expressing more gap junctions which enable calcium exchange. To understand this result, we use a mathematical model in which a bifurcation analysis has previously shown coupling-induced desynchronization in an oscillatory network. Using parameters close to this bifurcation but in the excitable regime, we find that this desynchronization persists and can explain the experimental observations. The model further predicts that co-culturing with gap-junction-deficient cells should restore synchronization, which experiments confirm. Combining quantitative experiments, the physical and biological manipulation of cells, and mathematical modeling, our results suggest that cell-to-cell connectivity significantly affects how populations encode an external temporal signal as it slows down: Sparse networks synchronize due to longer entrainment, whereas highly connected networks can desynchronize due to dynamic frustration. Published by the American Physical Society2025more » « less
-
Csikász-Nagy, Attila (Ed.)The cell cycle consists of a series of orchestrated events controlled by molecular sensing and feedback networks that ultimately drive the duplication of total DNA and the subsequent division of a single parent cell into two daughter cells. The ability to block the cell cycle and synchronize cells within the same phase has helped understand factors that control cell cycle progression and the properties of each individual phase. Intriguingly, when cells are released from a synchronized state, they do not maintain synchronized cell division and rapidly become asynchronous. The rate and factors that control cellular desynchronization remain largely unknown. In this study, using a combination of experiments and simulations, we investigate the desynchronization properties in cervical cancer cells (HeLa) starting from the G1/S boundary following double-thymidine block. Propidium iodide (PI) DNA staining was used to perform flow cytometry cell cycle analysis at regular 8 hour intervals, and a custom auto-similarity function to assess the desynchronization and quantify the convergence to an asynchronous state. In parallel, we developed a single-cell phenomenological model the returns the DNA amount across the cell cycle stages and fitted the parameters using experimental data. Simulations of population of cells reveal that the cell cycle desynchronization rate is primarily sensitive to the variability of cell cycle duration within a population. To validate the model prediction, we introduced lipopolysaccharide (LPS) to increase cell cycle noise. Indeed, we observed an increase in cell cycle variability under LPS stimulation in HeLa cells, accompanied with an enhanced rate of cell cycle desynchronization. Our results show that the desynchronization rate of artificially synchronized in-phase cell populations can be used a proxy of the degree of variance in cell cycle periodicity, an underexplored axis in cell cycle research.more » « less
-
In this article, we devise two related control algorithms to change the degree of synchrony of a population of noise-free, identical, uncoupled neural oscillators using a single control input. The algorithms are based on phase reduction, and use a population-level partial differential equation formulation to change the phase distribution of the neurons as desired. Motivated by the pathological neural synchronization hypothesized to be present in patients suffering from essential and parkinsonian tremor, we take our control objective to be the desynchronization of an initially synchronized neural population. Through numerical simulations, we are able to show that our algorithms work for both Type I and Type II neural populations. To demonstrate the versatility of our control algorithms, we also show that they can be applied to synchronize an initially desynchronized neural population as well. For the systems considered in this paper, the control algorithms can be applied to achieve any desired traveling-wave neural phase distribution, as long as the combination of initial and desired phase distributions is non-degenerate.more » « less
-
Rasmussen, M.; Reich, S.; Zaslavski, A. J. (Ed.)We study patterns observed right after the loss of stability of mixing in the Kuramoto model of coupled phase oscillators with random intrinsic frequencies on large graphs, which can also be random. We show that the emergent patterns are formed via two independent mechanisms determined by the shape of the frequency distribution and the limiting structure of the underlying graph sequence. Specifically, we identify two nested eigenvalue problems whose eigenvectors (unstable modes) determine the structure of the nascent patterns. The analysis is illustrated with the results of the numerical experiments with the Kuramoto model with unimodal and bimodal frequency distributions on certain graphs.more » « less
