Disentangling regular and chaotic motion in the standard map using complex network analysis of recurrences in phase space.

Recurrence in the phase space of complex systems is a well-studied phenomenon, which has provided deep insights into the nonlinear dynamics of such systems. For dissipative systems, characteristics based on recurrence plots have recently attracted much interest for discriminating qualitatively different types of dynamics in terms of measures of complexity, dynamical invariants, or even structural characteristics of the underlying attractor's geometry in phase space. Here, we demonstrate that the latter approach also provides a corresponding distinction between different co-existing dynamical regimes of the standard map, a paradigmatic example of a low-dimensional conservative system. Specifically, we show that the recently developed approach of recurrence network analysis provides potentially useful geometric characteristics distinguishing between regular and chaotic orbits. We find that chaotic orbits in an intermittent laminar phase (commonly referred to as sticky orbits) have a distinct geometric structure possibly differing in a subtle way from those of regular orbits, which is highlighted by different recurrence network properties obtained from relatively short time series. Thus, this approach can help discriminating regular orbits from laminar phases of chaotic ones, which presents a persistent challenge to many existing chaos detection techniques.

A max-plus model of ribosome dynamics during mRNA translation.

We examine the dynamics of the translation stage of cellular protein production, in which ribosomes move uni-directionally along an mRNA strand, building amino acid chains as they go. We describe the system using a timed event graph-a class of Petri net useful for studying discrete events, which have to satisfy constraints. We use max-plus algebra to describe a deterministic version of the model, where the constraints represent steric effects which prevent more than one ribosome reading a given codon at a given time and delays associated with the availability of the different tRNAs. We calculate the protein production rate and density of ribosomes on the mRNA and find exact agreement between these analytical results and numerical simulations of the deterministic model, even in the case of heterogeneous mRNAs.

The dynamics of supply and demand in mRNA translation.

We study the elongation stage of mRNA translation in eukaryotes and find that, in contrast to the assumptions of previous models, both the supply and the demand for tRNA resources are important for determining elongation rates. We find that increasing the initiation rate of translation can lead to the depletion of some species of aa-tRNA, which in turn can lead to slow codons and queueing. Particularly striking "competition" effects are observed in simulations of multiple species of mRNA which are reliant on the same pool of tRNA resources. These simulations are based on a recent model of elongation which we use to study the translation of mRNA sequences from the Saccharomyces cerevisiae genome. This model includes the dynamics of the use and recharging of amino acid tRNA complexes, and we show via Monte Carlo simulation that this has a dramatic effect on the protein production behaviour of the system.

Combinatorial stresses kill pathogenic Candida species.

Pathogenic microbes exist in dynamic niches and have evolved robust adaptive responses to promote survival in their hosts. The major fungal pathogens of humans, Candida albicans and Candida glabrata, are exposed to a range of environmental stresses in their hosts including osmotic, oxidative and nitrosative stresses. Significant efforts have been devoted to the characterization of the adaptive responses to each of these stresses. In the wild, cells are frequently exposed simultaneously to combinations of these stresses and yet the effects of such combinatorial stresses have not been explored. We have developed a common experimental platform to facilitate the comparison of combinatorial stress responses in C. glabrata and C. albicans. This platform is based on the growth of cells in buffered rich medium at 30°C, and was used to define relatively low, medium and high doses of osmotic (NaCl), oxidative (H(2)O(2)) and nitrosative stresses (e.g., dipropylenetriamine (DPTA)-NONOate). The effects of combinatorial stresses were compared with the corresponding individual stresses under these growth conditions. We show for the first time that certain combinations of combinatorial stress are especially potent in terms of their ability to kill C. albicans and C. glabrata and/or inhibit their growth. This was the case for combinations of osmotic plus oxidative stress and for oxidative plus nitrosative stress. We predict that combinatorial stresses may be highly significant in host defences against these pathogenic yeasts.

Limited resources in a driven diffusion process.

The advance of particles in many driven diffusion systems depends on the availability of resources in the surrounding environment. In the balance between supply and demand of such resources we are confronted with a regime in which, under limited resource availability, the flow is markedly reduced. In the context of mRNA translation this represents the finite availability of amino acid-tRNA molecules. In this limited resources regime a severe depletion of amino acid tRNAs is also observed. These dramatic effects are vital to our understanding of translation, and are likely to also be important for the many other applications of driven diffusion models.

Distinguishing direct from indirect interactions in oscillatory networks with multiple time scales.

We propose a method to infer the coupling structure in networks of nonlinear oscillatory systems with multiple time scales. The method of partial phase synchronization allows us to infer the coupling structure for coupled nonlinear oscillators with one well-defined time scale. The case of oscillators with multiple time scales has remained a challenge until now. Here, we introduce partial recurrence based synchronization analysis to tackle this challenge. We successfully apply the proposed method to model systems and experimental data from coupled electrochemical oscillators. The statistical significance of the results is evaluated based on a surrogate hypothesis test.

Abundance of stable periodic behavior in a Red Grouse population model with delay: a consequence of homoclinicity.

Shrimp-shaped periodic regions embedded in chaotic regions in two-dimensional parameter spaces are of specific interest for physical and biological systems. We provide the first observation of these shrimp-shaped stability regions in a parameter space of a continuous time-delayed population model, obtained by taking the delays as bifurcation parameters. The parameter space organization is governed by the presence of infinitely many periodicity hubs, which trigger the spiraling organization of these shrimp-shaped periodic regions around them. We provide evidence that this spiraling organization in the parameter space is a consequence of the existence of homoclinic orbits in the phase space.

Introduction to focus issue: dynamics in systems biology.

The methods of nonlinear systems form an extensive toolbox for the study of biology, and systems biology provides a rich source of motivation for the development of new mathematical techniques and the furthering of understanding of dynamical systems. This Focus Issue collects together a large variety of work which highlights the complementary nature of these two fields, showing what each has to offer the other. While a wide range of subjects is covered, the papers often have common themes such as "rhythms and oscillations," "networks and graph theory," and "switches and decision making." There is a particular emphasis on the links between experimental data and modeling and mathematical analysis.

Recurrences determine the dynamics.

We show that under suitable assumptions, Poincare recurrences of a dynamical system determine its topology in phase space. Therefore, dynamical systems with the same recurrences are dynamically equivalent. This conclusion can be drawn from a theorem proved in this paper which states that the recurrence matrix determines the topology of closed sets. The theorem states that if a set of points M is mapped onto another set N, such that two points in N are closer than some prescribed fixed distance if and only if the corresponding points in M are closer than some, in general different, prescribed fixed distance, then both sets are homeomorphic, i.e., identical up to a continuous change in the coordinate system. The theorem justifies a range of methods in nonlinear dynamics which are based on recurrence properties.

Hypothesis test for synchronization: twin surrogates revisited.

The method of twin surrogates has been introduced to test for phase synchronization of complex systems in the case of passive experiments. In this paper we derive new analytical expressions for the number of twins depending on the size of the neighborhood, as well as on the length of the trajectory. This allows us to determine the optimal parameters for the generation of twin surrogates. Furthermore, we determine the quality of the twin surrogates with respect to several linear and nonlinear statistics depending on the parameters of the method. In the second part of the paper we perform a hypothesis test for phase synchronization in the case of experimental data from fixational eye movements. These miniature eye movements have been shown to play a central role in neural information processing underlying the perception of static visual scenes. The high number of data sets (21 subjects and 30 trials per person) allows us to compare the generated twin surrogates with the "natural" surrogates that correspond to the different trials. We show that the generated twin surrogates reproduce very well all linear and nonlinear characteristics of the underlying experimental system. The synchronization analysis of fixational eye movements by means of twin surrogates reveals that the synchronization between the left and right eye is significant, indicating that either the centers in the brain stem generating fixational eye movements are closely linked, or, alternatively that there is only one center controlling both eyes.

Improving network inference: The impact of false positive and false negative conclusions about the presence or absence of links.

BACKGROUND: A reliable inference of networks from data is of key interest in the Neurosciences. Several methods have been suggested in the literature to reliably determine links in a network. To decide about the presence of links, these techniques rely on statistical inference, typically controlling the number of false positives, paying little attention to false negatives. NEW METHOD: In this paper, by means of a comprehensive simulation study, we analyse the influence of false positive and false negative conclusions about the presence or absence of links in a network on the network topology. We show that different values to balance false positive and false negative conclusions about links should be used in order to reliably estimate network characteristics. We propose to run careful simulation studies prior to making potentially erroneous conclusion about the network topology. RESULTS: Our analysis shows that optimal values to balance false positive and false negative conclusions about links depend on the network topology and characteristic of interest. COMPARISON WITH EXISTING METHODS: Existing methods rely on a choice of the rate for false positive conclusions. They aim to be sure about individual links rather than the entire network. The rate of false negative conclusions is typically not investigated. CONCLUSIONS: Our investigation shows that the balance of false positive and false negative conclusions about links in a network has to be tuned for any network topology that is to be estimated. Moreover, within the same network topology, the results are qualitatively the same for each network characteristic, but the actual values leading to reliable estimates of the characteristics are different.

Estimation of the direction of the coupling by conditional probabilities of recurrence.

We introduce a method to detect and quantify the asymmetry of the coupling between two interacting systems based on their recurrence properties. This method can detect the direction of the coupling in weakly as well as strongly coupled systems. It even allows detecting the asymmetry of the coupling in the more challenging case of structurally different systems and it is very robust against noise. We also address the problem of detecting the asymmetry of the coupling in passive experiments, i.e., when the strength of the coupling cannot be systematically changed, which is of great relevance for the analysis of experimental time series.

Network inference in the presence of latent confounders: the role of instantaneous causalities.

BACKGROUND: Detecting causal interactions in multivariate systems, in terms of Granger-causality, is of major interest in the Neurosciences. Typically, it is almost impossible to observe all components of the system. Missing certain components can lead to the appearance of spurious interactions. The aim of this study is to demonstrate the effect of this and to demonstrate that distinction between latent confounders and volume conduction is possible in some cases. NEW METHOD: Our new method uses a combination of renormalised partial directed coherence and analysis of the (partial) covariance matrix of residual noise process to detect instantaneous, spurious interactions. Sub-network analyses are performed to infer the true network structure of the underlying system. RESULTS: We provide evidence that it is possible to distinguish between instantaneous interactions that occur as a result of a latent confounder and those that occur as a result of volume conduction. COMPARISON WITH EXISTING METHODS: Our novel approach demonstrates to what extent inference of unobserved important processes as well as the distinction between latent confounders and volume conduction is possible. We suggest a combination of measures of Granger-causality and covariance selection models to achieve this numerically. CONCLUSIONS: Sub-network analyses enable a much more precise and correct inference of the true underlying network structure in some cases. From this it is possible to distinguish between unobserved processes and volume conduction. Our approach is straightforwardly adaptable to various measures of Granger-causality emphasising its ubiquitous successful applicability.

Community control in cellular protein production: consequences for amino acid starvation.

Deprivation of essential nutrients can have stark consequences for many processes in a cell. We consider amino acid starvation, which can result in bottlenecks in mRNA translation when ribosomes stall due to lack of resources, i.e. tRNAs charged with the missing amino acid. Recent experiments also show less obvious effects such as increased charging of other (non-starved) tRNA species and selective charging of isoaccepting tRNAs. We present a mechanism which accounts for these observations and shows that production of some proteins can actually increase under starvation. One might assume that such responses could only be a result of sophisticated control pathways, but here we show that these effects can occur naturally due to changes in the supply and demand for different resources, and that control can be accomplished through selective use of rare codons. We develop a model for translation which includes the dynamics of the charging and use of aminoacylated tRNAs, explicitly taking into account the effect of specific codon sequences. This constitutes a new control mechanism in gene regulation which emerges at the community level, i.e. via resources used by all ribosomes.

Networks: On the relation of bi- and multivariate measures.

A reliable inference of networks from observations of the nodes' dynamics is a major challenge in physics. Interdependence measures such as a the correlation coefficient or more advanced methods based on, e.g., analytic phases of signals are employed. For several of these interdependence measures, multivariate counterparts exist that promise to enable distinguishing direct and indirect connections. Here, we demonstrate analytically how bivariate measures relate to the respective multivariate ones; this knowledge will in turn be used to demonstrate the implications of thresholded bivariate measures for network inference. Particularly, we show, that random networks are falsely identified as small-world networks if observations thereof are treated by bivariate methods. We will employ the correlation coefficient as an example for such an interdependence measure. The results can be readily transferred to all interdependence measures partializing for information of thirds in their multivariate counterparts.

Assessing the strength of directed influences among neural signals: an approach to noisy data.

BACKGROUND: Measurements in the neurosciences are afflicted with observational noise. Granger-causality inference typically does not take this effect into account. We demonstrate that this leads to false positives conclusions and spurious causalities. NEW METHOD: State space modelling provides a convenient framework to obtain reliable estimates for Granger-causality. Despite its previous application in several studies, the analytical derivation of the statistics for parameter estimation in the state space model was missing. This prevented a rigorous evaluation of the results. RESULTS: In this manuscript we derive the statistics for parameter estimation in the state space model. We demonstrate in an extensive simulation study that our novel approach outperforms standard approaches and avoids false positive conclusions about Granger-causality. COMPARISON WITH EXISTING METHODS: In comparison with the naive application of Granger-causality inference, we demonstrate the superiority of our novel approach. The wide-spread applicability of our procedure provides a statistical framework for future studies. The application to mice electroencephalogram data demonstrates the immediate applicability of our approach. CONCLUSIONS: The analytical derivation of the statistics presented in this manuscript enables a rigorous evaluation of the results of Granger causal network inference. It is noteworthy that the statistics can be readily applied to various measures for Granger causality and other approaches that are based on vector autoregressive models.

Contribution of Fdh3 and Glr1 to Glutathione Redox State, Stress Adaptation and Virulence in Candida albicans.

The major fungal pathogen of humans, Candida albicans, is exposed to reactive nitrogen and oxygen species following phagocytosis by host immune cells. In response to these toxins, this fungus activates potent anti-stress responses that include scavenging of reactive nitrosative and oxidative species via the glutathione system. Here we examine the differential roles of two glutathione recycling enzymes in redox homeostasis, stress adaptation and virulence in C. albicans: glutathione reductase (Glr1) and the S-nitrosoglutathione reductase (GSNOR), Fdh3. We show that the NADPH-dependent Glr1 recycles GSSG to GSH, is induced in response to oxidative stress and is required for resistance to macrophage killing. GLR1 deletion increases the sensitivity of C. albicans cells to H2O2, but not to formaldehyde or NO. In contrast, Fdh3 detoxifies GSNO to GSSG and NH3, and FDH3 inactivation delays NO adaptation and increases NO sensitivity. C. albicans fdh3⎔ cells are also sensitive to formaldehyde, suggesting that Fdh3 also contributes to formaldehyde detoxification. FDH3 is induced in response to nitrosative, oxidative and formaldehyde stress, and fdh3Δ cells are more sensitive to killing by macrophages. Both Glr1 and Fdh3 contribute to virulence in the Galleria mellonella and mouse models of systemic infection. We conclude that Glr1 and Fdh3 play differential roles during the adaptation of C. albicans cells to oxidative, nitrosative and formaldehyde stress, and hence during the colonisation of the host. Our findings emphasise the importance of the glutathione system and the maintenance of intracellular redox homeostasis in this major pathogen.

Integrative Model of Oxidative Stress Adaptation in the Fungal Pathogen Candida albicans.

The major fungal pathogen of humans, Candida albicans, mounts robust responses to oxidative stress that are critical for its virulence. These responses counteract the reactive oxygen species (ROS) that are generated by host immune cells in an attempt to kill the invading fungus. Knowledge of the dynamical processes that instigate C. albicans oxidative stress responses is required for a proper understanding of fungus-host interactions. Therefore, we have adopted an interdisciplinary approach to explore the dynamical responses of C. albicans to hydrogen peroxide (H2O2). Our deterministic mathematical model integrates two major oxidative stress signalling pathways (Cap1 and Hog1 pathways) with the three major antioxidant systems (catalase, glutathione and thioredoxin systems) and the pentose phosphate pathway, which provides reducing equivalents required for oxidative stress adaptation. The model encapsulates existing knowledge of these systems with new genomic, proteomic, transcriptomic, molecular and cellular datasets. Our integrative approach predicts the existence of alternative states for the key regulators Cap1 and Hog1, thereby suggesting novel regulatory behaviours during oxidative stress. The model reproduces both existing and new experimental observations under a variety of scenarios. Time- and dose-dependent predictions of the oxidative stress responses for both wild type and mutant cells have highlighted the different temporal contributions of the various antioxidant systems during oxidative stress adaptation, indicating that catalase plays a critical role immediately following stress imposition. This is the first model to encapsulate the dynamics of the transcriptional response alongside the redox kinetics of the major antioxidant systems during H2O2 stress in C. albicans.

Optimized spectral estimation for nonlinear synchronizing systems.

In many fields of research nonlinear dynamical systems are investigated. When more than one process is measured, besides the distinct properties of the individual processes, their interactions are of interest. Often linear methods such as coherence are used for the analysis. The estimation of coherence can lead to false conclusions when applied without fulfilling several key assumptions. We introduce a data driven method to optimize the choice of the parameters for spectral estimation. Its applicability is demonstrated based on analytical calculations and exemplified in a simulation study. We complete our investigation with an application to nonlinear tremor signals in Parkinson's disease. In particular, we analyze electroencephalogram and electromyogram data.

The impact of latent confounders in directed network analysis in neuroscience.

In the analysis of neuroscience data, the identification of task-related causal relationships between various areas of the brain gives insights about the network of physiological pathways that are active during the task. One increasingly used approach to identify causal connectivity uses the concept of Granger causality that exploits predictability of activity in one region by past activity in other regions of the brain. Owing to the complexity of the data, selecting components for the analysis of causality as a preprocessing step has to be performed. This includes predetermined-and often arbitrary-exclusion of information. Therefore, the system is confounded by latent sources. In this paper, the effect of latent confounders is demonstrated, and paths of influence among three components are studied. While methods for analysing Granger causality are commonly based on linear vector autoregressive models, the effects of latent confounders are expected to be present also in nonlinear systems. Therefore, all analyses are also performed for a simulated nonlinear system and discussed with regard to applications in neuroscience.