Dynamical quorum sensing: Population density encoded in cellular dynamics

Kinetic modelling of glycolytic oscillations

Glycolytic oscillations have been studied for well over 60 years, but aspects of their function, and mechanisms of regulation and synchronisation remain unclear. Glycolysis is amenable to mechanistic mathematical modelling, as its components have been well characterised, and the system can be studied at many organisational levels: in vitro reconstituted enzymes, cell free extracts, individual cells, and cell populations. In recent years, the emergence of individual cell analysis has opened new ways of studying this intriguing system.

Glycolytic oscillations under periodic drivings

In many living organisms displaying circadian rhythms, the intake of energy often occurs in a periodic manner. Glycolysis is a prototypical biochemical reaction that exhibits a self-sustained oscillation under continuous injection of glucose. Here we study the effect of periodic injection of glucose on the glycolytic oscillation from a dynamical systems perspective. In particular, we employ Goldbeter's allosteric model of phosphofructokinase as a model system for glycolytic oscillations, and explore the effect of periodic substrate influx of varying frequencies and amplitudes by building the phase diagrams of Lyapunov exponents and oscillatory periods. When the frequency of driving is tuned around the harmonic and sub/super-harmonic conditions of the natural frequency, the system is entrained to a frequency-locked state, forming an entrainment band that broadens with an increasing amplitude of driving. On the other hand, if the amplitude is substantial, the system may transition, albeit infrequent, to a chaotic state which defies prediction of dynamical behaviour. Our study offers in-depth understandings into the controllability of glycolytic oscillation as well as explaining physical underpinnings that enable the synchronous oscillations among a dense population of cells.

Mobility-induced order in active XY spins on a substrate

We elucidate that the nearly phase-ordered active XY spins in contact with a conserved, diffusing species on a substrate can be stable. For wide-ranging model parameters, it has stable uniform phases robust against noises. These are distinguished by generalized quasi-long-range (QLRO) orientational order logarithmically stronger or weaker than the well-known QLRO in equilibrium, together with miniscule (i.e., hyperuniform) or giant number fluctuations, respectively. This illustrates a direct correspondence between the two. The scaling of both phase and density fluctuations in the stable phase-ordered states is nonuniversal: they depend on the nonlinear dynamical couplings. For other parameters, it has no stable uniformly ordered phase. Our model, a theory for active spinners, provides a minimal framework for wide-ranging systems, e.g., active superfluids on substrates, synchronization of oscillators, active carpets of cilia and bacterial flagella, and active membranes.

Active XY model on a substrate: Density fluctuations and phase ordering

We explore the generic long-wavelength properties of an active XY model on a substrate, consisting of a collection of nearly phase-ordered active XY spins in contact with a diffusing, conserved species, as a representative system of active spinners with a conservation law. The spins rotate actively in response to the local density fluctuations and local phase differences, on a solid substrate. We investigate this system by Monte Carlo simulations of an agent-based model, which we set up, complemented by the hydrodynamic theory for the system. We demonstrate that this system can phase-synchronize without any hydrodynamic interactions. Our combined numerical and analytical studies show that this model, when stable, displays hitherto unstudied scaling behavior: As a consequence of the interplay between the mobility, active rotation, and number conservation, such a system can be stable over a wide range of the model parameters characterized by a novel correspondence between the phase and density fluctuations. In different regions of the phase space where the phase-ordered system is stable, it displays generalized quasi-long-range order (QLRO): It shows phase ordering which is generically either logarithmically stronger than the conventional QLRO found in its equilibrium limit, together with "miniscule number fluctuations," or logarithmically weaker than QLRO along with "giant number fluctuations," showing a novel one-to-one correspondence between phase ordering and density fluctuations in the ordered states. Intriguingly, these scaling exponents are found to depend explicitly on the model parameters. We further show that in other parameter regimes there are no stable, ordered phases. Instead, two distinct types of disordered states with short-range phase order are found, characterized by the presence or absence of stable clusters of finite sizes. In a surprising connection, the hydrodynamic theory for this model also describes the fluctuations in a Kardar-Parisi-Zhang (KPZ) surface with a conserved species on it, or an active fluid membrane with a finite tension, without momentum conservation and a conserved species living on it. This implies the existence of stable fluctuating surfaces that are only logarithmically smoother or rougher than the Edward-Wilkinson surface at two dimensions (2D) can exist, in contrast to the 2D pure KPZ-like "rough" surfaces.

Coevolutionary Dynamics with Global Fields

We investigate the effects of external and autonomous global interaction fields on an adaptive network of social agents with an opinion formation dynamics based on a simple imitation rule. We study the competition between global fields and adaptive rewiring on the space of parameters of the system. The model represents an adaptive society subject to global mass media such as a directed opinion influence or feedback of endogenous cultural trends. We show that, in both situations, global mass media contribute to consensus and to prevent the fragmentation of the social network induced by the coevolutionary dynamics. We present a discussion of these results in the context of dynamical systems and opinion formation dynamics.

The macroscopic limit to synchronization of cellular clocks in single cells of Neurospora crassa

We determined the macroscopic limit for phase synchronization of cellular clocks in an artificial tissue created by a “big chamber” microfluidic device to be about 150,000 cells or less. The dimensions of the microfluidic chamber allowed us to calculate an upper limit on the radius of a hypothesized quorum sensing signal molecule of 13.05 nm using a diffusion approximation for signal travel within the device. The use of a second microwell microfluidic device allowed the refinement of the macroscopic limit to a cell density of 2166 cells per fixed area of the device for phase synchronization. The measurement of averages over single cell trajectories in the microwell device supported a deterministic quorum sensing model identified by ensemble methods for clock phase synchronization. A strong inference framework was used to test the communication mechanism in phase synchronization of quorum sensing versus cell-to-cell contact, suggesting support for quorum sensing. Further evidence came from showing phase synchronization was density-dependent.

Principles, mechanisms and functions of entrainment in biological oscillators

Entrainment is a phenomenon in which two oscillators interact with each other, typically through physical or chemical means, to synchronize their oscillations. This phenomenon occurs in biology to coordinate processes from the molecular to organismal scale. Biological oscillators can be entrained within a single cell, between cells or to an external input. Using six illustrative examples of entrainable biological oscillators, we discuss the distinctions between entrainment and synchrony and explore features that contribute to a system's propensity to entrain. Entrainment can either enhance or reduce the heterogeneity of oscillations within a cell population, and we provide examples and mechanisms of each case. Finally, we discuss the known functions of entrainment and discuss potential functions from an evolutionary perspective.

Collective Behavior of Urease pH Clocks in Nano- and Microvesicles Controlled by Fast Ammonia Transport

The transmission of chemical signals via an extracellular solution plays a vital role in collective behavior in cellular biological systems and may be exploited in applications of lipid vesicles such as drug delivery. Here, we investigated chemical communication in synthetic micro- and nanovesicles containing urease in a solution of urea and acid. We combined experiments with simulations to demonstrate that the fast transport of ammonia to the external solution governs the pH-time profile and synchronizes the timing of the pH clock reaction in a heterogeneous population of vesicles. This study shows how the rate of production and emission of a small basic product controls pH changes in active vesicles with a distribution of sizes and enzyme amounts, which may be useful in bioreactor or healthcare applications.

Pattern forming systems coupling linear bulk diffusion to dynamically active membranes or cells

Some analytical and numerical results are presented for pattern formation properties associated with novel types of reaction-diffusion (RD) systems that involve the coupling of bulk diffusion in the interior of a multi-dimensional spatial domain to nonlinear processes that occur either on the domain boundary or within localized compartments that are confined within the domain. The class of bulk-membrane system considered herein is derived from an asymptotic analysis in the limit of small thickness of a thin domain that surrounds the bulk medium. When the bulk domain is a two-dimensional disk, a weakly nonlinear analysis is used to characterize Turing and Hopf bifurcations that can arise from the linearization around a radially symmetric, but spatially non-uniform, steady-state of the bulk-membrane system. In a singularly perturbed limit, the existence and linear stability of localized membrane-bound spike patterns is analysed for a Gierer-Meinhardt activator-inhibitor model that includes bulk coupling. Finally, the emergence of collective intracellular oscillations is studied for a class of PDE-ODE bulk-cell model in a bounded two-dimensional domain that contains spatially localized, but dynamically active, circular cells that are coupled through a linear bulk diffusion field. Applications of such coupled bulk-membrane or bulk-cell systems to some biological systems are outlined, and some open problems in this area are discussed. This article is part of the theme issue 'Recent progress and open frontiers in Turing's theory of morphogenesis'.

Disorder in cellular packing can alter proliferation dynamics to regulate growth

The mechanisms by which an organ regulates its growth are not yet fully understood, especially when the cells are closely packed as in epithelial tissues. We explain growth arrest as a collective dynamical transition in coupled oscillators on disordered lattices. As the cellular morphologies become homogeneous over the course of development, the signals induced by cell-cell contact increase beyond a critical value that triggers coordinated cessation of the cell-cycle oscillators driving cell division. Thus, control of cell proliferation is causally related to the geometry of cellular packing.

Quorum-sensing induced transitions between bistable steady-states for a cell-bulk ODE-PDE model with lux intracellular kinetics

Intercellular signaling and communication are used by bacteria to regulate a variety of behaviors. In a type of cell-cell communication known as quorum sensing (QS), which is mediated by a diffusible signaling molecule called an autoinducer, bacteria can undergo sudden changes in their behavior at a colony-wide level when the density of cells exceeds a critical threshold. In mathematical models of QS behavior, these changes can include the switch-like emergence of intracellular oscillations through a Hopf bifurcation, or sudden transitions between bistable steady-states as a result of a saddle-node bifurcation of equilibria. As an example of this latter type of QS transition, we formulate and analyze a cell-bulk ODE-PDE model in a 2-D spatial domain that incorporates the prototypical LuxI/LuxR QS system for a collection of stationary bacterial cells, as modeled by small circular disks of a common radius with a cell membrane that is permeable only to the autoinducer. By using the method of matched asymptotic expansions, it is shown that the steady-state solutions for the cell-bulk model exhibit a saddle-node bifurcation structure. The linear stability of these branches of equilibria are determined from the analysis of a nonlinear matrix eigenvalue problem, called the globally coupled eigenvalue problem. The key role on QS behavior of a bulk degradation of the autoinducer field, which arises from either a Robin boundary condition on the domain boundary or from a constant bulk decay, is highlighted. With bulk degradation, it is shown analytically that the effect of coupling identical bacterial cells to the bulk autoinducer diffusion field is to create an effective bifurcation parameter that depends on the population of the colony, the bulk diffusivity, the membrane permeabilities, and the cell radius. QS transitions occur when this effective parameter passes through a saddle-node bifurcation point of the Lux ODE kinetics for an isolated cell. In the limit of a large but finite bulk diffusivity, it is shown that the cell-bulk system is well-approximated by a simpler ODE-DAE system. This reduced system, which is used to study the effect of cell location on QS behavior, is easily implemented for a large number of cells. Predictions from the asymptotic theory for QS transitions between bistable states are favorably compared with full numerical solutions of the cell-bulk ODE-PDE system.

Synchronisation of glycolytic activity in yeast cells

Glycolysis is the central metabolic pathway of almost every cell and organism. Under appropriate conditions, glycolytic oscillations may occur in individual cells as well as in entire cell populations or tissues. In many biological systems, glycolytic oscillations drive coherent oscillations of other metabolites, for instance in cardiomyocytes near anorexia, or in pancreas where they lead to a pulsatile release of insulin. Oscillations at the population or tissue level require the cells to synchronize their metabolism. We review the progress achieved in studying a model organism for glycolytic oscillations, namely yeast. Oscillations may occur on the level of individual cells as well as on the level of the cell population. In yeast, the cell-to-cell interaction is realized by diffusion-mediated intercellular communication via a messenger molecule. The present mini-review focuses on the synchronisation of glycolytic oscillations in yeast. Synchronisation is a quorum-sensing phenomenon because the collective oscillatory behaviour of a yeast cell population ceases when the cell density falls below a threshold. We review the question, under which conditions individual cells in a sparse population continue or cease to oscillate. Furthermore, we provide an overview of the pathway leading to the onset of synchronized oscillations. We also address the effects of spatial inhomogeneities (e.g., the formation of spatial clusters) on the collective dynamics, and also review the emergence of travelling waves of glycolytic activity. Finally, we briefly review the approaches used in numerical modelling of synchronized cell populations.

Arnold tongues in oscillator systems with nonuniform spatial driving

Nonlinear oscillator systems are ubiquitous in biology and physics, and their control is a practical problem in many experimental systems. Here we study this problem in the context of the two models of spatially coupled oscillators: the complex Ginzburg-Landau equation (CGLE) and a generalization of the CGLE in which oscillators are coupled through an external medium (emCGLE). We focus on external control drives that vary in both space and time. We find that the spatial distribution of the drive signal controls the frequency ranges over which oscillators synchronize to the drive and that boundary conditions strongly influence synchronization to external drives for the CGLE. Our calculations also show that the emCGLE has a low density regime in which a broad range of frequencies can be synchronized for low drive amplitudes. We study the bifurcation structure of these models and find that they are very similar to results for the driven Kuramoto model, a system with no spatial structure. We conclude by discussing qualitative implications of our results for controlling coupled oscillator systems such as the social amoebae Dictyostelium and populations of Belousov Zhabotinsky (BZ) catalytic particles using spatially structured external drives.

Emergence of multistability and strongly asymmetric collective modes in two quorum sensing coupled identical ring oscillators

The simplest ring oscillator is made from three strongly nonlinear elements repressing each other unidirectionally, resulting in the emergence of a limit cycle. A popular implementation of this scheme uses repressor genes in bacteria, creating the synthetic genetic oscillator known as the Repressilator. We consider the main collective modes produced when two identical Repressilators are mean-field-coupled via the quorum-sensing mechanism. In-phase and anti-phase oscillations of the coupled oscillators emerge from two Andronov-Hopf bifurcations of the homogeneous steady state. Using the rate of the repressor's production and the value of coupling strength as the bifurcation parameters, we performed one-parameter continuations of limit cycles and two-parameter continuations of their bifurcations to show how bifurcations of the in-phase and anti-phase oscillations influence the dynamical behaviors for this system. Pitchfork bifurcation of the unstable in-phase cycle leads to the creation of novel inhomogeneous limit cycles with very different amplitudes, in contrast to the well-known asymmetrical limit cycles arising from oscillation death. The Neimark-Sacker bifurcation of the anti-phase cycle determines the border of an island in two-parameter space containing almost all the interesting regimes including the set of resonant limit cycles, the area with stable inhomogeneous cycle, and very large areas with chaotic regimes resulting from torus destruction and period doubling of resonant cycles and inhomogeneous cycles. We discuss the structure of the chaos skeleton to show the role of inhomogeneous cycles in its formation. Many regions of multistability and transitions between regimes are presented. These results provide new insights into the coupling-dependent mechanisms of multistability and collective regime symmetry breaking in populations of identical multidimensional oscillators.

Partial synchronisation of glycolytic oscillations in yeast cell populations

The transition between synchronized and asynchronous behaviour of immobilized yeast cells of the strain Saccharomyces carlsbergensis was investigated by monitoring the autofluorescence of the coenzyme NADH. In populations of intermediate cell densities the individual cells remained oscillatory, whereas on the level of the cell population both a partially synchronized and an asynchronous state were accessible for experimental studies. In the partially synchronized state, the mean oscillatory frequency was larger than that of the cells in the asynchronous state. This suggests that synchronisation occurred due to entrainment by the cells that oscillated more rapidly. This is typical for synchronisation due to phase advancement. Furthermore, the synchronisation of the frequency of the glycolytic oscillations preceded the synchronisation of their phases. However, the cells did not synchronize completely, as the distribution of the oscillatory frequencies only narrowed but did not collapse to a unique frequency. Cells belonging to spatially denser clusters showed a slightly enhanced local synchronisation during the episode of partial synchronisation. Neither the clusters nor a transition from partially synchronized glycolytic oscillations to travelling glycolytic waves did substantially affect the degree of partial synchronisation. Chimera states, i.e., the coexistence of a synchronized and an asynchronous part of the population, could not be found.

Phase synchronization of fluid-fluid interfaces as hydrodynamically coupled oscillators

Hydrodynamic interactions play a role in synchronized motions of coupled oscillators in fluids, and understanding the mechanism will facilitate development of applications in fluid mechanics. For example, synchronization phenomenon in two-phase flow will benefit the design of future microfluidic devices, allowing spatiotemporal control of microdroplet generation without additional integration of control elements. In this work, utilizing a characteristic oscillation of adjacent interfaces between two immiscible fluids in a microfluidic platform, we discover that the system can act as a coupled oscillator, notably showing spontaneous in-phase synchronization of droplet breakup. With this observation of in-phase synchronization, the coupled droplet generator exhibits a complete set of modes of coupled oscillators, including out-of-phase synchronization and nonsynchronous modes. We present a theoretical model to elucidate how a negative feedback mechanism, tied to the distance between the interfaces, induces the in-phase synchronization. We also identify the criterion for the transition from in-phase to out-of-phase oscillations.

The robust production of droplets by microfluidic T-junctions is a well-established technique. Um et al. demonstrate how the mutual interaction between droplets can be used to achieve additional control including the simultaneous release of droplets caused by synchronization phenomena.

Cellular metabolism and colloids: Realistically linking physiology and biological physical chemistry

Important concepts from colloidal physical chemistry such as coacervation, phase transitions, emergent properties and ionic association, are currently emerging in the lexicon of cellular biology, prompted mostly by recent experimental observations of liquid phase coexistence in the cell cytosol. Nevertheless, from an historical point of view, the application of these concepts in cell biology is not new. They were key concepts into the so-called protoplasmic doctrine, an alternative (and largely forgotten) approach to cell physiology. The most complete theory originating from this line of thinking was the Association-Induction Hypothesis (AIH), introduced by Gilbert N. Ling in 1962. The AIH, which envisions living cells as complex dynamical colloidal systems, provides ample theory and experimental evidence to call into question the now dominant view of living cells as fluid-filled vesicles. This review attempts to present and discuss the usefulness of the AIH to understand a series of experimental observations from our laboratory from living suspensions of the yeast Saccharomyces cerevisiae exhibiting glycolytic oscillations. Particularly, the AIH helped us integrate, in a mechanistic sense, the basis of a strong temporal coupling observed between ATP and a series of cellular properties such as intracellular water dipolar relaxation, intracellular K⁺ concentration, among many others, where the colloidal physical chemistry of the cell interior plays a fundamental role.

Coordinating an Ensemble of Chemical Micromotors via Spontaneous Synchronization

Spatiotemporal coordination of a nanorobot ensemble is critical for their operation in complex environments, such as tissue removal or drug delivery. Current strategies of achieving this task, however, rely heavily on sophisticated, external manipulation. We here present an alternative, biomimetic strategy by which oscillating Ag Janus micromotors spontaneously synchronize their dynamics as chemically coupled oscillators. By quantitatively tracking the kinetics at both an individual and cluster level, we find that synchronization emerges as the oscillating entities are increasingly coupled as they approach each other. In addition, the synchronized beating of a cluster of these oscillating colloids was found to be dominated by substrate electroosmosis, revealed with the help of an acoustic trapping technique. This quantitative, systematic study of synchronizing micromotors could facilitate the design of biomimetic nanorobots that spontaneously communicate and organize at micro- and nanoscales. It also serves as a model system for nonlinear active matter.

Symmetry breakings in two populations of oscillators coupled via diffusive environments: Chimera and heterosynchrony

We consider two diffusively coupled populations of identical oscillators, where the oscillators in each population are coupled with a common dynamic environment. Existence and stability of a variety of stationary states are analyzed on the basis of the Ott-Antonsen reduction method, which reveals that the chimera state occurs under the diffusive coupling scheme. Furthermore, we find an exotic symmetry-breaking behavior, the so-called the heterosynchronous state, in which each population exhibits in-phase coherence, while the order parameters of two populations rotate at different phase velocities. The chimera and heterosynchronous states emerge from bistabilities of distinct states for decoupled population and occur as a unique continuation for weak diffusive couplings. The heterosynchronous state is caused by an indirect coupling scheme via dynamic environments and could occur for a finite-size system as well, even for the system that consists of one oscillator per population.

Evolution of quasiperiodicity in quorum-sensing coupled identical repressilators

The dynamics of three three-dimensional repressilators globally coupled by a quorum sensing mechanism was numerically studied. This number (three) of coupled repressilators is sufficient to obtain such a set of self-consistent oscillation frequencies of signal molecules in the mean field that results in the appearance of self-organized quasiperiodicity and its complex evolution over wide areas of model parameters. Numerically analyzing the invariant curves as a function of coupling strength, we observed torus doubling, three torus arising via quasiperiodic Hopf bifurcation, the emergence of resonant cycles, and secondary Neimark-Sacker bifurcation. A gradual increase in the oscillation amplitude leads to chaotizations of the tori and to the birth of weak, but multidimensional chaos.

Emergence of collective oscillations in adaptive cells

Collective oscillations of cells in a population appear under diverse biological contexts. Here, we establish a set of common principles by categorising the response of individual cells against a time-varying signal. A positive intracellular signal relay of sufficient gain from participating cells is required to sustain the oscillations, together with phase matching. The two conditions yield quantitative predictions for the onset cell density and frequency in terms of measured single-cell and signal response functions. Through mathematical constructions, we show that cells that adapt to a constant stimulus fulfil the phase requirement by developing a leading phase in an active frequency window that enables cell-to-signal energy flow. Analysis of dynamical quorum sensing in several cellular systems with increasing biological complexity reaffirms the pivotal role of adaptation in powering oscillations in an otherwise dissipative cell-to-cell communication channel. The physical conditions identified also apply to synthetic oscillatory systems.

There are many examples of cell populations exhibiting density-dependent collective oscillatory behaviour. Here, the authors show that sustained collective oscillations emerge when cells anticipate variation in signal and attempt to amplify it, a property that can be linked to adaptation.

Explosive death in nonlinear oscillators coupled by quorum sensing

Many biological and chemical systems exhibit collective behavior in response to the change in their population density. These elements or cells communicate with each other via dynamical agents or signaling molecules. In this work, we explore the dynamics of nonlinear oscillators, specifically Stuart-Landau oscillators and Rayleigh oscillators, interacting globally through dynamical agents in the surrounding environment modeled as a quorum sensing interaction. The system exhibits the typical continuous second-order transition from oscillatory state to death state, when the oscillation amplitude is small. However, interestingly, when the amplitude of oscillations is large we find that the system shows an abrupt transition from oscillatory to death state, a transition termed "explosive death." So the quorum-sensing form of interaction can induce the usual second-order transition, as well as sudden first-order transitions. Further, in the case of the explosive death transitions, the oscillatory state and the death state coexist over a range of coupling strengths near the transition point. This emergent regime of hysteresis widens with increasing strength of the mean-field feedback, and is relevant to hysteresis that is widely observed in biological, chemical, and physical processes.

Phase transition to synchronization in generalized Kuramoto model with low-pass filter

A second-order continuous synchronization has been well documented for the classic Kuramoto model. Here we generalize the classic Kuramoto model by incorporating a low-pass filter (LPF) in the coupling, which serves as a simple form of indirect coupling through a common external dynamic environment. We uncover that a first-order explosive synchronization turns out to be a very generic phenomenon in this generalized Kuramoto model with LPF. We establish theoretical results by providing a rigorous analytical treatment, which is validated by conducting extensive numerical simulations. Our study provides a new root for the emergence of first-order explosive synchronization, which could substantially deepen the understanding of the underlying mechanism of a first-order phase transition towards synchronization in coupled dynamical networks.

Selection of spatiotemporal patterns in arrays of spatially distributed oscillators indirectly coupled via a diffusive environment

Emergence of self-organized behaviors in diverse living systems often depends on population density. In these systems, cell-cell communications are usually mediated by the surrounding environment. Collective behaviors (e.g., synchrony and dynamical quorum sensing) of such systems with stirred environment have been extensively studied, but the spatiotemporal dynamics of the oscillators coupled via a diffusive environment (without stirring) is rather understudied. We here perform a computational study on the selection and competition of wave patterns in arrays of spatially distributed oscillators immersed in a diffusive medium. We find that population density plays a crucial role in the selection of wave patterns: (i) for a single spiral in the system, its rotation either inward or outward could be controlled by population density, and (ii) for spiral and target waves coexisting initially in the system, wave competition happens and population density decides which type of wave will finally survive. The latter phenomenon is further confirmed in a system whose individual element is excitable rather than self-sustained oscillatory. The mechanism underlying all these observations is attributed to the frequency competition. Our results in the excitable case may have implications on the experimental results.

Oscillatory activity regulation in an ensemble of autonomous mercury beating heart oscillators

Collective behavior of an ensemble of directly or indirectly coupled oscillators can be a function of population density. Experiments using autonomous mercury beating heart (MBH) oscillators coupled through their surroundings are employed, to study the existence of quorum-like (population dependent) phenomena. Two coupling mechanisms are used, namely, static and dynamic coupling. For the static coupling scheme, the transitions of a subset of the coupled oscillators occur from active (oscillatory) to inactive (quiescent) state and vice versa. A continuous variation of collective dynamics was observed as the population of the oscillators increased. For the dynamic coupling scheme, the time for which the coupled oscillators are active changes sharply as the population increases beyond a certain threshold.

Modeling studies of heterogeneities in glycolytic oscillations in HeLa cervical cancer cells

Previous experiments demonstrated that a population of HeLa cells starved of glucose or both glucose and serum exhibited a strong heterogeneity in the glycolytic oscillations in terms of the number of oscillatory cells, periods of oscillations, and duration of oscillations. Here, we report numerical simulations of this heterogeneous oscillatory behavior in HeLa cells by using a newly developed mathematical model. It is simple enough that we can apply a mathematical analysis, but capture the core of the glycolytic pathway and the activity of the glucose transporter (GLUT). Lognormal distributions of the values of the four rate constants in the model were obtained from the experimental distributions in the periods of oscillations. Thus, the heterogeneity in the periods of oscillations can be attributed to the difference in the rate constants of the enzymatic reactions. The activity of GLUT is found to determine whether the HeLa cells were oscillatory or non-oscillatory under the same experimental conditions. Simulation with the log-normal distribution of the maximum uptake velocity of glucose and the four randomized rate constants based on the log-normal distributions successfully reproduced the time-dependent number of oscillatory cells (oscillatory ratios) under the two starving conditions. The difference in the initial values of the metabolites has little effect on the simulated results.

Abrupt and gradual onset of synchronized oscillations due to dynamical quorum sensing in the single-cathode multi-anode nickel electrodissolution system

The nonlinear dynamics of an oscillatory Ni electrodissolution-hydrogen ion reduction system are explored in a multi-electrode anode-single cathode system. A mathematical analysis of the charge balance equations reveals that the coupling scheme is similar to dynamical quorum sensing, where the number of anode wires affects a parameter related to the population density. In a parameter region where the large population exhibits stationary behavior, with sufficiently strong coupling (with small individual resistances attached to the anode wires), synchronized oscillations emerge abruptly with decreasing the number of anodes. Therefore, an "inverse" dynamical quorum sensing takes place. With weak coupling the transition is gradual. The experiments are supported by numerical simulation of a kinetic model of the process. The results thus show that the description of nontrivial cathode-anode interactions in the form of dynamical quorum sensing provides an efficient way of analyzing the dynamical response of complex, interacting electrochemical reactions.

Influence of reaction-induced convection on quorum sensing in enzyme-loaded agarose beads

In theory, groups of enzyme-loaded particles producing an acid or base may show complex behavior including dynamical quorum sensing, the appearance of synchronized oscillations above a critical number or density of particles. Here, experiments were performed with the enzyme urease loaded into mm-sized agarose beads and placed in a solution of urea, resulting in an increase in pH. This behavior was found to be dependent upon the number of beads present in the array; however, reaction-induced convection occurred and plumes of high pH developed that extended to the walls of the reactor. The convection resulted in the motion of the mm-sized particles and conversion of the solution to high pH. Simulations in a simple model of the beads demonstrated the suppression of dynamical quorum sensing in the presence of flow.

Measuring Lyapunov exponents of large chaotic systems with global coupling by time series analysis

Despite the prominent importance of the Lyapunov exponents for characterizing chaos, it still remains a challenge to measure them for large experimental systems, mainly because of the lack of recurrences in time series analysis. Here, we develop a method to overcome this difficulty, valid for highly symmetric systems such as systems with global coupling, for which the dimensionality of recurrence analysis can be reduced drastically. We test our method numerically with two globally coupled systems, namely, logistic maps and limit-cycle oscillators with global coupling. The evaluated exponent values are successfully compared with the true ones obtained by the standard numerical method. We also describe a few techniques to improve the accuracy of the proposed method.

Collective excitability in a mesoscopic neuronal model of epileptic activity

At the mesoscopic scale, the brain can be understood as a collection of interacting neuronal oscillators, but the extent to which its sustained activity is due to coupling among brain areas is still unclear. Here we address this issue in a simplified situation by examining the effect of coupling between two cortical columns described via Jansen-Rit neural mass models. Our results show that coupling between the two neuronal populations gives rise to stochastic initiations of sustained collective activity, which can be interpreted as epileptic events. For large enough coupling strengths, termination of these events results mainly from the emergence of synchronization between the columns, and thus it is controlled by coupling instead of noise. Stochastic triggering and noise-independent durations are characteristic of excitable dynamics, and thus we interpret our results in terms of collective excitability.

Bistable emergence of oscillations in growing Bacillus subtilis biofilms

Biofilm communities of Bacillus subtilis bacteria have recently been shown to exhibit collective growth-rate oscillations mediated by electrochemical signaling to cope with nutrient starvation. These oscillations emerge once the colony reaches a large enough number of cells. However, it remains unclear whether the amplitude of the oscillations, and thus their effectiveness, builds up over time gradually or if they can emerge instantly with a nonzero amplitude. Here we address this question by combining microfluidics-based time-lapse microscopy experiments with a minimal theoretical description of the system in the form of a delay-differential equation model. Analytical and numerical methods reveal that oscillations arise through a subcritical Hopf bifurcation, which enables instant high-amplitude oscillations. Consequently, the model predicts a bistable regime where an oscillating and a nonoscillating attractor coexist in phase space. We experimentally validate this prediction by showing that oscillations can be triggered by perturbing the media conditions, provided the biofilm size lies within an appropriate range. The model also predicts that the minimum size at which oscillations start decreases with stress, a fact that we also verify experimentally. Taken together, our results show that collective oscillations in cell populations can emerge suddenly with nonzero amplitude via a discontinuous transition.

Dissipative structures in biological systems: bistability, oscillations, spatial patterns and waves

The goal of this review article is to assess how relevant is the concept of dissipative structure for understanding the dynamical bases of non-equilibrium self-organization in biological systems, and to see where it has been applied in the five decades since it was initially proposed by Ilya Prigogine. Dissipative structures can be classified into four types, which will be considered, in turn, and illustrated by biological examples: (i) multistability, in the form of bistability and tristability, which involve the coexistence of two or three stable steady states, or in the form of birhythmicity, which involves the coexistence between two stable rhythms; (ii) temporal dissipative structures in the form of sustained oscillations, illustrated by biological rhythms; (iii) spatial dissipative structures, known as Turing patterns; and (iv) spatio-temporal structures in the form of propagating waves. Rhythms occur with widely different periods at all levels of biological organization, from neural, cardiac and metabolic oscillations to circadian clocks and the cell cycle; they play key roles in physiology and in many disorders. New rhythms are being uncovered while artificial ones are produced by synthetic biology. Rhythms provide the richest source of examples of dissipative structures in biological systems. Bistability has been observed experimentally, but has primarily been investigated in theoretical models in an increasingly wide range of biological contexts, from the genetic to the cell and animal population levels, both in physiological conditions and in disease. Bistable transitions have been implicated in the progression between the different phases of the cell cycle and, more generally, in the process of cell fate specification in the developing embryo. Turing patterns are exemplified by the formation of some periodic structures in the course of development and by skin stripe patterns in animals. Spatio-temporal patterns in the form of propagating waves are observed within cells as well as in intercellular communication. This review illustrates how dissipative structures of all sorts abound in biological systems.This article is part of the theme issue 'Dissipative structures in matter out of equilibrium: from chemistry, photonics and biology (part 1)'.

Universal relations of local order parameters for partially synchronized oscillators

Interactions among discrete oscillatory units (e.g., cells) can result in partially synchronized states when some of the units exhibit phase locking and others phase slipping. Such states are typically characterized by a global order parameter that expresses the extent of synchrony in the system. Here we show that such states carry data-rich information of the system behavior, and a local order parameter analysis reveals universal relations through a semicircle representation. The universal relations are derived from thermodynamic limit analysis of a globally coupled Kuramoto-type phase oscillator model. The relations are confirmed with the partially synchronized states in numerical simulations with a model of circadian cells and in laboratory experiments with chemical oscillators. The application of the theory allows direct approximation of coupling strength, the natural frequency of oscillations, and the phase lag parameter without extensive nonlinear fits as well as a self-consistency check for presence of network interactions and higher harmonic components in the phase model.

Stopping Autoinducer-2 Chatter by Means of an Indigenous Bacterium ( Acinetobacter sp. DKY-1): A New Antibiofouling Strategy in a Membrane Bioreactor for Wastewater Treatment

Bacterial quorum quenching (QQ) by means of degrading signaling molecules has been applied to antibiofouling strategies in a membrane bioreactor (MBR) for wastewater treatment. However, the target signaling molecules have been limited to N-acyl homoserine lactones participating in intraspecies quorum sensing. Here, an approach to disrupting autoinducer-2 (AI-2) signaling molecules participating in interspecies quorum sensing was pursued as a next-generation antibiofouling strategy in an MBR for wastewater treatment. We isolated an indigenous QQ bacterium ( Acinetobacter sp. DKY-1) that can attenuate the expression of the quorum-sensing (QS) response through the inactivation of an autoinducer-2 signaling molecule, 4,5-dihydroxy-2,3-pentanedione (DPD), among four kinds of autoinducer-2 QS bacteria. DKY-1 released AI-2 QQ compounds, which were verified to be hydrophilic with a molecular weight of <400 Da. The addition of DKY-1 entrapping beads into an MBR significantly decreased DPD concentration and remarkably reduced membrane biofouling. This new approach, combining molecular biology with wastewater engineering, could enlarge the range of QQ-MBR for antibiofouling and energy savings in the field of wastewater treatment.

Is a constant low-entropy process at the root of glycolytic oscillations?

We measured temporal oscillations in thermodynamic variables such as temperature, heat flux, and cellular volume in suspensions of non-dividing yeast cells which exhibit temporal glycolytic oscillations. Oscillations in these variables have the same frequency as oscillations in the activity of intracellular metabolites, suggesting strong coupling between them. These results can be interpreted in light of a recently proposed theoretical formalism in which isentropic thermodynamic systems can display coupled oscillations in all extensive and intensive variables, reminiscent of adiabatic waves. This interpretation suggests that oscillations may be a consequence of the requirement of living cells for a constant low-entropy state while simultaneously performing biochemical transformations, i.e., remaining metabolically active. This hypothesis, which is in line with the view of the cellular interior as a highly structured and near equilibrium system where energy inputs can be low and sustain regular oscillatory regimes, calls into question the notion that metabolic processes are essentially dissipative.

Promotion and inhibition of synchronous glycolytic oscillations in yeast by chitosan

Synchronous rhythmic activities play crucial roles in diverse biological systems. Glycolytic oscillations in yeast cells have been studied for 50 years with the aim of elucidating the mechanisms underlying the intracellular oscillations and their synchronization. We investigated the effects of chemical disturbances on the individual and collective glycolytic oscillations in yeast cells encapsulated in alginate microparticles, and demonstrated that the addition of chitosan, an antimicrobial agent, decreased the duration of these oscillations. In contrast, the periods and the synchronicity states showed two different responses to the chitosan treatments. The periods were shown to be prolonged following the treatment with 5-50 mg·L^-1 and shortened at 75 mg·L^-1 of chitosan. Collective oscillations became more synchronized at 5 mg·L^-1 of chitosan, and desynchronized at 25-75 mg·L^-1 of this compound. These findings can be explained by the balance between two chitosan features, increasing cell membrane permeability and acetaldehyde scavenging. At low concentrations, chitosan presumably acts as a synchronization promoter that does not mediate the synchronization itself but induces an increase in intercellular coupling. We believe that our findings may provide new insights into the synchronous rhythmic activities in biological systems.

Statistical Dynamics of Spatial-Order Formation by Communicating Cells

Communicating cells can coordinate their gene expressions to form spatial patterns, generating order from disorder. Ubiquitous "secrete-and-sense cells" secrete and sense the same molecule to do so. Here we present a modeling framework-based on cellular automata and mimicking approaches of statistical mechanics-for understanding how secrete-and-sense cells with bistable gene expression, from disordered beginnings, can become spatially ordered by communicating through rapidly diffusing molecules. Classifying lattices of cells by two "macrostate" variables-"spatial index," measuring degree of order, and average gene-expression level-reveals a conceptual picture: a group of cells behaves as a single particle, in an abstract space, that rolls down on an adhesive "pseudo-energy landscape" whose shape is determined by cell-cell communication and an intracellular gene-regulatory circuit. Particles rolling down the landscape represent cells becoming more spatially ordered. We show how to extend this framework to more complex forms of cellular communication.

Switches induced by quorum sensing in a model of enzyme-loaded microparticles

Quorum sensing refers to the ability of bacteria and other single-celled organisms to respond to changes in cell density or number with population-wide changes in behaviour. Here, simulations were performed to investigate quorum sensing in groups of diffusively coupled enzyme microparticles using a well-characterized autocatalytic reaction which raises the pH of the medium: hydrolysis of urea by urease. The enzyme urease is found in both plants and microorganisms, and has been widely exploited in engineering processes. We demonstrate how increases in group size can be used to achieve a sigmoidal switch in pH at high enzyme loading, oscillations in pH at intermediate enzyme loading and a bistable, hysteretic switch at low enzyme loading. Thus, quorum sensing can be exploited to obtain different types of response in the same system, depending on the enzyme concentration. The implications for microorganisms in colonies are discussed, and the results could help in the design of synthetic quorum sensing for biotechnology applications such as drug delivery.

The dynamics of intracellular water constrains glycolytic oscillations in Saccharomyces cerevisiae

We explored the dynamic coupling of intracellular water with metabolism in yeast cells. Using the polarity-sensitive probe 6-acetyl-2-dimethylaminonaphthalene (ACDAN), we show that glycolytic oscillations in the yeast S. cerevisiae BY4743 wild-type strain are coupled to the generalized polarization (GP) function of ACDAN, which measures the physical state of intracellular water. We analysed the oscillatory dynamics in wild type and 24 mutant strains with mutations in many different enzymes and proteins. Using fluorescence spectroscopy, we measured the amplitude and frequency of the metabolic oscillations and ACDAN GP in the resting state of all 25 strains. The results showed that there is a lower and an upper threshold of ACDAN GP, beyond which oscillations do not occur. This critical GP range is also phenomenologically linked to the occurrence of oscillations when cells are grown at different temperatures. Furthermore, the link between glycolytic oscillations and the ACDAN GP value also holds when ATP synthesis or the integrity of the cell cytoskeleton is perturbed. Our results represent the first demonstration that the dynamic behaviour of a metabolic process can be regulated by a cell-wide physical property: the dynamic state of intracellular water, which represents an emergent property.

Dissipative structures and biological rhythms

Sustained oscillations abound in biological systems. They occur at all levels of biological organization over a wide range of periods, from a fraction of a second to years, and with a variety of underlying mechanisms. They control major physiological functions, and their dysfunction is associated with a variety of physiological disorders. The goal of this review is (i) to give an overview of the main rhythms observed at the cellular and supracellular levels, (ii) to briefly describe how the study of biological rhythms unfolded in the course of time, in parallel with studies on chemical oscillations, (iii) to present the major roles of biological rhythms in the control of physiological functions, and (iv) the pathologies associated with the alteration, disappearance, or spurious occurrence of biological rhythms. Two tables present the main examples of cellular and supracellular rhythms ordered according to their period, and their role in physiology and pathophysiology. Among the rhythms discussed are neural and cardiac rhythms, metabolic oscillations such as those occurring in glycolysis in yeast, intracellular Ca⁺⁺ oscillations, cyclic AMP oscillations in Dictyostelium amoebae, the segmentation clock that controls somitogenesis, pulsatile hormone secretion, circadian rhythms which occur in all eukaryotes and some bacteria with a period close to 24 h, the oscillatory dynamics of the enzymatic network driving the cell cycle, and oscillations in transcription factors such as NF-ΚB and tumor suppressors such as p53. Ilya Prigogine's concept of dissipative structures applies to temporal oscillations and allows us to unify within a common framework the various rhythms observed at different levels of biological organization, regardless of their period and underlying mechanism.

Primordial oscillations in life: Direct observation of glycolytic oscillations in individual HeLa cervical cancer cells

We report the first direct observation of glycolytic oscillations in HeLa cervical cancer cells, which we regard as primordial oscillations preserved in living cells. HeLa cells starved of glucose or both glucose and serum exhibited glycolytic oscillations in nicotinamide adenine dinucleotide (NADH), exhibiting asynchronous intercellular behaviors. Also found were spatially homogeneous and inhomogeneous intracellular NADH oscillations in the individual cells. Our results demonstrate that starved HeLa cells may be induced to exhibit glycolytic oscillations by either high-uptake of glucose or the enhancement of a glycolytic pathway (Crabtree effect or the Warburg effect), or both. Their asynchronous collective behaviors in the oscillations were probably due to a weak intercellular coupling. Elucidation of the relationship between the mechanism of glycolytic dynamics in cancer cells and their pathophysiological characteristics remains a challenge in future.

Emergence of localized patterns in globally coupled networks of relaxation oscillators with heterogeneous connectivity

Oscillations in far-from-equilibrium systems (e.g., chemical, biochemical, biological) are generated by the nonlinear interplay of positive and negative feedback effects operating at different time scales. Relaxation oscillations emerge when the time scales between the activators and the inhibitors are well separated. In addition to the large-amplitude oscillations (LAOs) or relaxation type, these systems exhibit small-amplitude oscillations (SAOs) as well as abrupt transitions between them (canard phenomenon). Localized cluster patterns in networks of relaxation oscillators consist of one cluster oscillating in the LAO regime or exhibiting mixed-mode oscillations (LAOs interspersed with SAOs), while the other oscillates in the SAO regime. Because the individual oscillators are monostable, localized patterns are a network phenomenon that involves the interplay of the connectivity and the intrinsic dynamic properties of the individual nodes. Motivated by experimental and theoretical results on the Belousov-Zhabotinsky reaction, we investigate the mechanisms underlying the generation of localized patterns in globally coupled networks of piecewise-linear relaxation oscillators where the global feedback term affects the rate of change of the activator (fast variable) and depends on the weighted sum of the inhibitor (slow variable) at any given time. We also investigate whether these patterns are affected by the presence of a diffusive type of coupling whose synchronizing effects compete with the symmetry-breaking global feedback effects.

Synthetic quorum sensing in model microcapsule colonies

Biological quorum sensing refers to the ability of cells to gauge their population density and collectively initiate a new behavior once a critical density is reached. Designing synthetic materials systems that exhibit quorum sensing-like behavior could enable the fabrication of devices with both self-recognition and self-regulating functionality. Herein, we develop models for a colony of synthetic microcapsules that communicate by producing and releasing signaling molecules. Production of the chemicals is regulated by a biomimetic negative feedback loop, the "repressilator" network. Through theory and simulation, we show that the chemical behavior of such capsules is sensitive to both the density and number of capsules in the colony. For example, decreasing the spacing between a fixed number of capsules can trigger a transition in chemical activity from the steady, repressed state to large-amplitude oscillations in chemical production. Alternatively, for a fixed density, an increase in the number of capsules in the colony can also promote a transition into the oscillatory state. This configuration-dependent behavior of the capsule colony exemplifies quorum-sensing behavior. Using our theoretical model, we predict the transitions from the steady state to oscillatory behavior as a function of the colony size and capsule density.

Frequency and phase synchronization in large groups: Low dimensional description of synchronized clapping, firefly flashing, and cricket chirping

A common observation is that large groups of oscillatory biological units often have the ability to synchronize. A paradigmatic model of such behavior is provided by the Kuramoto model, which achieves synchronization through coupling of the phase dynamics of individual oscillators, while each oscillator maintains a different constant inherent natural frequency. Here we consider the biologically likely possibility that the oscillatory units may be capable of enhancing their synchronization ability by adaptive frequency dynamics. We propose a simple augmentation of the Kuramoto model which does this. We also show that, by the use of a previously developed technique [Ott and Antonsen, Chaos 18, 037113 (2008)], it is possible to reduce the resulting dynamics to a lower dimensional system for the macroscopic evolution of the oscillator ensemble. By employing this reduction, we investigate the dynamics of our system, finding a characteristic hysteretic behavior and enhancement of the quality of the achieved synchronization.

Noise-Induced and Control of Collective Behavior in a Population of Coupled Chemical Oscillators

Synchronization of intercommunicating individual oscillators is an important form of collective behavior used in nature as a mechanism to face dangers, act collectively, and communicate. The involvement of the medium where oscillators exist is an important ingredient. Because of their nature and their multiple different components, the medium and the environment are often perceived as stochastic relative to the deterministic nature of the individuals on some scale. This injects energy/matter into the system in ways that can enhance or de-enhance communication in a stochastic manner. Here we experimentally consider a large number of coupled nonlinear-chemical oscillators under the effect of a controlled normally distributed noise. Experiments show that the collective behavior of the oscillator is triggered by this stochastic perturbation, and we observe the dependence on the noise parameters. Our results point to the potential use of environmental fluctuations in determining the emergence and properties of collective behaviors in complex systems.

Weak synchronization and large-scale collective oscillation in dense bacterial suspensions

Collective oscillatory behaviour is ubiquitous in nature, having a vital role in many biological processes from embryogenesis and organ development to pace-making in neuron networks. Elucidating the mechanisms that give rise to synchronization is essential to the understanding of biological self-organization. Collective oscillations in biological multicellular systems often arise from long-range coupling mediated by diffusive chemicals, by electrochemical mechanisms, or by biomechanical interaction between cells and their physical environment. In these examples, the phase of some oscillatory intracellular degree of freedom is synchronized. Here, in contrast, we report the discovery of a weak synchronization mechanism that does not require long-range coupling or inherent oscillation of individual cells. We find that millions of motile cells in dense bacterial suspensions can self-organize into highly robust collective oscillatory motion, while individual cells move in an erratic manner, without obvious periodic motion but with frequent, abrupt and random directional changes. So erratic are individual trajectories that uncovering the collective oscillations of our micrometre-sized cells requires individual velocities to be averaged over tens or hundreds of micrometres. On such large scales, the oscillations appear to be in phase and the mean position of cells typically describes a regular elliptic trajectory. We found that the phase of the oscillations is organized into a centimetre-scale travelling wave. We present a model of noisy self-propelled particles with strictly local interactions that accounts faithfully for our observations, suggesting that self-organized collective oscillatory motion results from spontaneous chiral and rotational symmetry breaking. These findings reveal a previously unseen type of long-range order in active matter systems (those in which energy is spent locally to produce non-random motion). This mechanism of collective oscillation may inspire new strategies to control the self-organization of active matter and swarming robots.

Communication shapes sensory response in multicellular networks

Collective sensing by interacting cells is observed in a variety of biological systems, and yet, a quantitative understanding of how sensory information is collectively encoded is lacking. Here, we investigate the ATP-induced calcium dynamics of monolayers of fibroblast cells that communicate via gap junctions. Combining experiments and stochastic modeling, we find that increasing the ATP stimulus increases the propensity for calcium oscillations, despite large cell-to-cell variability. The model further predicts that the oscillation propensity increases with not only the stimulus, but also the cell density due to increased communication. Experiments confirm this prediction, showing that cell density modulates the collective sensory response. We further implicate cell-cell communication by coculturing the fibroblasts with cancer cells, which we show act as "defects" in the communication network, thereby reducing the oscillation propensity. These results suggest that multicellular networks sit at a point in parameter space where cell-cell communication has a significant effect on the sensory response, allowing cells to simultaneously respond to a sensory input and the presence of neighbors.