Diffusive coupling facilitates and impedes noise-induced escape in interacting bistable elements.

Diverse complex systems often undergo sudden changes in their states, such as epileptic seizures, climate changes, and social uprisings. Such behavior has been modeled by noise-induced escape of bistable elements, which is the escape from an attracting state driven by a fluctuation in the system's state. We consider a system of interacting bistable elements and investigate the effect of diffusive coupling among elements on the process of noise-induced escape. We focus on the influence of the coupling strength over the escape time, which is the time it takes for noise-induced escape to occur. We performed numerical simulations and observed that weak coupling reduced the mean escape time, whereas strong coupling impeded escape. We argue that, although diffusive coupling both facilitates and impedes escape, the facilitating effect is dominant when coupling is weak. For weak coupling cases, we develop an approximate theory that can predict the mean and variance of escape times. In contrast, strong coupling reduces the effective noise intensity to impede escape. Our results suggest that diffusive coupling among multistable elements contributes to regulating the rate of transitions among attracting states.

Weakly nonlinear analysis on synchronization and oscillation quenching of coupled mechanical oscillators.

Various oscillatory phenomena occur in the world. Because some are associated with abnormal states (e.g. epilepsy), it is important to establish ways to terminate oscillations by external stimuli. However, despite the prior development of techniques for stabilizing unstable oscillations, relatively few studies address the transition from oscillatory to resting state in nonlinear dynamics. This study mainly analyzes the oscillation-quenching of metronomes on a platform as an example of such transitions. To facilitate the analysis, we describe the impulsive force (escapement mechanism) of a metronome by a fifth-order polynomial. By performing both averaging approximation and numerical simulation, we obtain a phase diagram for synchronization and oscillation quenching. We find that quenching occurs when the feedback to the oscillator increases, which will help explore the general principle regarding the state transition from oscillatory to resting state. We also numerically investigate the bifurcation of out-of-phase synchronization and beat-like solution. Despite the simplicity, our model successfully reproduces essential phenomena in interacting mechanical clocks, such as the bistability of in-phase and anti-phase synchrony and oscillation quenching occurring for a large mass ratio between the oscillator and the platform. We believe that our simple model will contribute to future analyses of other dynamics of mechanical clocks.

Synchronization and stability analysis of an exponentially diverging solution in a mathematical model of asymmetrically interacting agents.

This study deals with an existing mathematical model of asymmetrically interacting agents. We analyze the following two previously unfocused features of the model: (i) synchronization of growth rates and (ii) initial value dependence of damped oscillation. By applying the techniques of variable transformation and timescale separation, we perform the stability analysis of a diverging solution. We find that (i) all growth rates synchronize to the same value that is as small as the smallest growth rate and (ii) oscillatory dynamics appear if the initial value of the slowest-growing agent is sufficiently small. Furthermore, our analytical method proposes a way to apply stability analysis to an exponentially diverging solution, which we believe is also a contribution of this study. Although the employed model is originally proposed as a model of infectious disease, we do not discuss its biological relevance but merely focus on the technical aspects.

An extended Hilbert transform method for reconstructing the phase from an oscillatory signal.

Rhythmic activity is ubiquitous in biological systems from the cellular to organism level. Reconstructing the instantaneous phase is the first step in analyzing the essential mechanism leading to a synchronization state from the observed signals. A popular method of phase reconstruction is based on the Hilbert transform, which can only reconstruct the interpretable phase from a limited class of signals, e.g., narrow band signals. To address this issue, we propose an extended Hilbert transform method that accurately reconstructs the phase from various oscillatory signals. The proposed method is developed by analyzing the reconstruction error of the Hilbert transform method with the aid of Bedrosian's theorem. We validate the proposed method using synthetic data and show its systematically improved performance compared with the conventional Hilbert transform method with respect to accurately reconstructing the phase. Finally, we demonstrate that the proposed method is potentially useful for detecting the phase shift in an observed signal. The proposed method is expected to facilitate the study of synchronization phenomena from experimental data.

Partial synchronization and community switching in phase-oscillator networks and its analysis based on a bidirectional, weighted chain of three oscillators.

Complex networks often possess communities defined based on network connectivity. When dynamics undergo in a network, one can also consider dynamical communities, i.e., a group of nodes displaying a similar dynamical process. We have investigated both analytically and numerically the development of a dynamical community structure, where the community is referred to as a group of nodes synchronized in frequency, in networks of phase oscillators. We first demonstrate that using a few example networks, the community structure changes when network connectivity or interaction strength is varied. In particular, we found that community switching, i.e., a portion of oscillators change the group to which they synchronize, occurs for a range of parameters. We then propose a three-oscillator model: a bidirectional, weighted chain of three Kuramoto phase oscillators, as a theoretical framework for understanding the community formation and its variation. Our analysis demonstrates that the model shows a variety of partially synchronized patterns: oscillators with similar natural frequencies tend to synchronize for weak coupling, while tightly connected oscillators tend to synchronize for strong coupling. We obtain approximate expressions for the critical coupling strengths by employing a perturbative approach in a weak coupling regime and a geometric approach in strong coupling regimes. Moreover, we elucidate the bifurcation types of transitions between different patterns. Our theory might be useful for understanding the development of partially synchronized patterns in a wider class of complex networks than community structured networks.

A light-induced small G-protein gem limits the circadian clock phase-shift magnitude by inhibiting voltage-dependent calcium channels.

Calcium signaling is pivotal to the circadian clockwork in the suprachiasmatic nucleus (SCN), particularly in rhythm entrainment to environmental light-dark cycles. Here, we show that a small G-protein Gem, an endogenous inhibitor of high-voltage-activated voltage-dependent calcium channels (VDCCs), is rapidly induced by light in SCN neurons via the calcium (Ca2+)-mediated CREB/CRE transcriptional pathway. Gem attenuates light-induced calcium signaling through its interaction with VDCCs. The phase-shift magnitude of locomotor activity rhythms by light, at night, increases in Gem-deficient (Gem-/-) mice. Similarly, in SCN slices from Gem-/- mice, depolarizing stimuli induce larger phase shifts of clock gene transcription rhythms that are normalized by the application of an L-type VDCC blocker, nifedipine. Voltage-clamp recordings from SCN neurons reveal that Ca2+ currents through L-type channels increase in Gem-/- mice. Our findings suggest that transcriptionally activated Gem feeds back to suppress excessive light-evoked L-type VDCC activation, adjusting the light-induced phase-shift magnitude to an appropriate level in mammals.

Phase and frequency linear response theory for hyperbolic chaotic oscillators.

We formulate a linear phase and frequency response theory for hyperbolic flows, which generalizes phase response theory for autonomous limit cycle oscillators to hyperbolic chaotic dynamics. The theory is based on a shadowing conjecture, stating the existence of a perturbed trajectory shadowing every unperturbed trajectory on the system attractor for any small enough perturbation of arbitrary duration and a corresponding unique time isomorphism, which we identify as phase such that phase shifts between the unperturbed trajectory and its perturbed shadow are well defined. The phase sensitivity function is the solution of an adjoint linear equation and can be used to estimate the average change of phase velocity to small time dependent or independent perturbations. These changes in frequency are experimentally accessible, giving a convenient way to define and measure phase response curves for chaotic oscillators. The shadowing trajectory and the phase can be constructed explicitly in the tangent space of an unperturbed trajectory using co-variant Lyapunov vectors. It can also be used to identify the limits of the regime of linear response.

Accuracy of a one-dimensional reduction of dynamical systems on networks.

Resilience is an ability of a system with which the system can adjust its activity to maintain its functionality when it is perturbed. To study resilience of dynamics on networks, Gao et al. [Nature (London) 530, 307 (2016)0028-083610.1038/nature16948] proposed a theoretical framework to reduce dynamical systems on networks, which are high dimensional in general, to one-dimensional dynamical systems. The accuracy of this one-dimensional reduction relies on three approximations in addition to the assumption that the network has a negligible degree correlation. In the present study, we analyze the accuracy of the one-dimensional reduction assuming networks without degree correlation. We do so mainly through examining the validity of the individual assumptions underlying the method. Across five dynamical system models, we find that the accuracy of the one-dimensional reduction hinges on the spread of the equilibrium value of the state variable across the nodes in most cases. Specifically, the one-dimensional reduction tends to be accurate when the dispersion of the node's state is small. We also find that the correlation between the node's state and the node's degree, which is common for various dynamical systems on networks, is unrelated to the accuracy of the one-dimensional reduction.

Noninvasive inference methods for interaction and noise intensities of coupled oscillators using only spike time data.

Measurements of interaction intensity are generally achieved by observing responses to perturbations. In biological and chemical systems, external stimuli tend to deteriorate their inherent nature, and thus, it is necessary to develop noninvasive inference methods. In this paper, we propose theoretical methods to infer coupling strength and noise intensity simultaneously in two well-synchronized noisy oscillators through observations of spontaneously fluctuating events such as neural spikes. A phase oscillator model is applied to derive formulae relating each of the parameters to spike time statistics. Using these formulae, each parameter is inferred from a specific set of statistics. We verify these methods using the FitzHugh-Nagumo model as well as the phase model. Our methods do not require external perturbations and thus can be applied to various experimental systems.

Feedback-induced desynchronization and oscillation quenching in a population of globally coupled oscillators.

Motivated from a wide range of applications, various methods to control synchronization in coupled oscillators have been proposed. Previous studies have demonstrated that global feedback typically induces three macroscopic behaviors: synchronization, desynchronization, and oscillation quenching. However, analyzing all of these transitions within a single theoretical framework is difficult, and thus the feedback effect is only partially understood in each framework. Herein, we analyze a model of globally coupled phase oscillators exposed to global feedback, which shows all of the typical macroscopic dynamical states. Analytical tractability of the model enables us to obtain detailed phase diagrams where transitions and bistabilities between different macroscopic states are identified. Additionally, we propose strategies to steer the oscillators into targeted states with minimal feedback strength. Our study provides a useful overview of the effect of global feedback and is expected to serve as a benchmark when more sophisticated feedback needs to be designed.

A Liposomal Gemcitabine, FF-10832, Improves Plasma Stability, Tumor Targeting, and Antitumor Efficacy of Gemcitabine in Pancreatic Cancer Xenograft Models.

PURPOSE: The clinical application of gemcitabine (GEM) is limited by its pharmacokinetic properties. The aim of this study was to characterize the stability in circulating plasma, tumor targeting, and payload release of liposome-encapsulated GEM, FF-10832. METHODS: Antitumor activity was assessed in xenograft mouse models of human pancreatic cancer. The pharmacokinetics of GEM and its active metabolite dFdCTP were also evaluated. RESULTS: In mice with Capan-1 tumors, the dose-normalized areas under the curve (AUCs) after FF-10832 administration in plasma and tumor were 672 and 1047 times higher, respectively, than after using unencapsulated GEM. The tumor-to-bone marrow AUC ratio of dFdCTP was approximately eight times higher after FF-10832 administration than after GEM administration. These results indicated that liposomal encapsulation produced long-term stability in circulating plasma and tumor-selective targeting of GEM. In mice with Capan-1, SUIT-2, and BxPC-3 tumors, FF-10832 had better antitumor activity and tolerability than GEM. Internalization of FF-10832 in tumor-associated macrophages (TAMs) was revealed by flow cytometry and confocal laser scanning microscopy, and GEM was efficiently released from isolated macrophages of mice treated with FF-10832. These results suggest that TAMs are one of the potential reservoirs of GEM in tumors. CONCLUSION: This study found that FF-10832 had favorable pharmacokinetic properties. The liposomal formulation was more effective and tolerable than unencapsulated GEM in mouse xenograft tumor models. Hence, FF-10832 is a promising candidate for the treatment of pancreatic cancer.

The effects of coating culture dishes with collagen on fibroblast cell shape and swirling pattern formation.

Motile human-skin fibroblasts form macroscopic swirling patterns when grown to confluence on a culture dish. In this paper, we investigate the effect of coating the culture-dish surface with collagen on the resulting pattern, using human-skin fibroblast NB1RGB cells as the model system. The presence of the collagen coating is expected to enhance the adherence of the fibroblasts to the dish surface, and thereby also enhance the traction that the fibroblasts have as they move. We find that, contrary to our initial expectation, the coating does not significantly affect the motility of the fibroblasts. Their eventual number density at confluence is also unaffected. However, the coherence length of cell orientation in the swirling pattern is diminished. We also find that the fibroblasts cultured in collagen-coated dishes are rounder in shape and shorter in perimeter, compared with those cultured in uncoated polystyrene or glass culture dishes. We hypothesise that the rounder cell-shape which weakens the cell-cell nematic contact interaction is responsible for the change in coherence length. A simple mathematical model of the migrating fibroblasts is constructed, which demonstrates that constant motility with weaker nematic interaction strength does indeed lead to the shortening of the coherence length.

Noise stability of synchronization and optimal network structures.

We provide a theoretical framework for quantifying the expected level of synchronization in a network of noisy oscillators. Through linearization around the synchronized state, we derive the following quantities as functions of the eigenvalues and eigenfunctions of the network Laplacian using a standard technique for dealing with multivariate Ornstein-Uhlenbeck processes: the magnitude of the fluctuations around a synchronized state and the disturbance coefficients αi that represent how strongly node i disturbs the synchronization. With this approach, we can quantify the effect of individual nodes and links on synchronization. Our theory can thus be utilized to find the optimal network structure for accomplishing the best synchronization. Furthermore, when the noise levels of the oscillators are heterogeneous, we can also find optimal oscillator configurations, i.e., where to place oscillators in a given network depending on their noise levels. We apply our theory to several example networks to elucidate optimal network structures and oscillator configurations.

Coupling delay controls synchronized oscillation in the segmentation clock.

Individual cellular activities fluctuate but are constantly coordinated at the population level via cell-cell coupling. A notable example is the somite segmentation clock, in which the expression of clock genes (such as Hes7) oscillates in synchrony between the cells that comprise the presomitic mesoderm (PSM)1,2. This synchronization depends on the Notch signalling pathway; inhibiting this pathway desynchronizes oscillations, leading to somite fusion3-7. However, how Notch signalling regulates the synchronicity of HES7 oscillations is unknown. Here we establish a live-imaging system using a new fluorescent reporter (Achilles), which we fuse with HES7 to monitor synchronous oscillations in HES7 expression in the mouse PSM at a single-cell resolution. Wild-type cells can rapidly correct for phase fluctuations in HES7 oscillations, whereas the absence of the Notch modulator gene lunatic fringe (Lfng) leads to a loss of synchrony between PSM cells. Furthermore, HES7 oscillations are severely dampened in individual cells of Lfng-null PSM. However, when Lfng-null PSM cells were completely dissociated, the amplitude and periodicity of HES7 oscillations were almost normal, which suggests that LFNG is involved mostly in cell-cell coupling. Mixed cultures of control and Lfng-null PSM cells, and an optogenetic Notch signalling reporter assay, revealed that LFNG delays the signal-sending process of intercellular Notch signalling transmission. These results-together with mathematical modelling-raised the possibility that Lfng-null PSM cells shorten the coupling delay, thereby approaching a condition known as the oscillation or amplitude death of coupled oscillators8. Indeed, a small compound that lengthens the coupling delay partially rescues the amplitude and synchrony of HES7 oscillations in Lfng-null PSM cells. Our study reveals a delay control mechanism of the oscillatory networks involved in somite segmentation, and indicates that intercellular coupling with the correct delay is essential for synchronized oscillation.

Enhanced lysosomal degradation maintains the quiescent state of neural stem cells.

Quiescence is important for sustaining neural stem cells (NSCs) in the adult brain over the lifespan. Lysosomes are digestive organelles that degrade membrane receptors after they undergo endolysosomal membrane trafficking. Enlarged lysosomes are present in quiescent NSCs (qNSCs) in the subventricular zone of the mouse brain, but it remains largely unknown how lysosomal function is involved in the quiescence. Here we show that qNSCs exhibit higher lysosomal activity and degrade activated EGF receptor by endolysosomal degradation more rapidly than proliferating NSCs. Chemical inhibition of lysosomal degradation in qNSCs prevents degradation of signaling receptors resulting in exit from quiescence. Furthermore, conditional knockout of TFEB, a lysosomal master regulator, delays NSCs quiescence in vitro and increases NSC proliferation in the dentate gyrus of mice. Taken together, our results demonstrate that enhanced lysosomal degradation is an important regulator of qNSC maintenance.

Partial synchronization of relaxation oscillators with repulsive coupling in autocatalytic integrate-and-fire model and electrochemical experiments.

Experiments and supporting theoretical analysis are presented to describe the synchronization patterns that can be observed with a population of globally coupled electrochemical oscillators close to a homoclinic, saddle-loop bifurcation, where the coupling is repulsive in the electrode potential. While attractive coupling generates phase clusters and desynchronized states, repulsive coupling results in synchronized oscillations. The experiments are interpreted with a phenomenological model that captures the waveform of the oscillations (exponential increase) followed by a refractory period. The globally coupled autocatalytic integrate-and-fire model predicts the development of partially synchronized states that occur through attracting heteroclinic cycles between out-of-phase two-cluster states. Similar behavior can be expected in many other systems where the oscillations occur close to a saddle-loop bifurcation, e.g., with Morris-Lecar neurons.

A reduced cell-based phase model for tissue polarity alignment through global anisotropic cues.

Ordered polarity alignment of cell populations plays vital roles in biology, such as in hair follicle alignment and asymmetric cell division. Although cell polarity is uniformly oriented along a tissue axis in many tissues, its mechanism is not well understood. In this paper, we propose a theoretical framework to understand the generic dynamical properties of polarity alignment in interacting cellular units, where each cell is described by a reaction-diffusion system, and the cells further interact with one another through the contacting surfaces between them. Using a perturbation method under the assumption of weak coupling between cells, we derive a reduced model in which polarity of each cell is described by only one variable. Essential dynamical properties including the effects of cell shape, coupling heterogeneity, external signal and noise can be clarified analytically. In particular, we show that the anisotropicity of the system, such as oriented cell elongation and axial asymmetry in the coupling strength, can serve as a global cue that drives the uniform orientation of cell polarity along a certain axis. Our study bridges the gap between detailed and phenomenological models, and it is expected to facilitate the study of polarity dynamics in various nonequilibrium systems.

Analyzing orientation patterns of human skin fibroblasts using image gradient information.

Fibroblasts are the cells within the dermis layer of skin which are responsible for generating connective tissue, where collagen is the main component. However, the effects of collagen on the motion behavior of human skin fibroblasts were not studied well. We therefore propose an image analysis method to quantify the effects of collagen on the motion characteristics of fibroblasts. Images are divided into disjoint blocks and three gradient based descriptors, namely histogram of oriented gradients (HOG), the modified HOG (modiHOG), and the gradient direction histogram (GDH) are then applied to each block. These features can extract the orientation and shape or orientation only information of the skin fibroblasts from images. An orientation similarity function and a quantitative index, called first minimum similarity distance, are defined to quantify the orientation behavior of fibroblasts. An experiment with two sets of fibroblasts images, obtained from cell culture with and without collagen coating, demonstrates that the skin fibroblasts produce smaller orientation patterns, when they are cultured in the presence of collagen proteins. A statistical t-test is then performed, which validates the significant differences in the orientation behavior of the skin fibroblasts in the two cases mentioned.

Segmentation Genes Enter an Excited State.

The expression of segmentation clock genes is self-oscillatory due to a delayed negative-feedback mechanism. Hubaud et al. (2017), reporting in Cell, now reveal that segmentation clock genes also exhibit properties of an excitable system, with Yap signaling setting the threshold for adopting an oscillatory state and Notch signaling triggering oscillations.