Dynamical properties induced by state-dependent delays in photonic systems.
Ontology highlight
ABSTRACT:
SUBMITTER: Martinez-Llinas J
PROVIDER: S-EPMC4557356 | biostudies-other | 2015
REPOSITORIES: biostudies-other
Similar Datasets
Project description:Farey graphs are simultaneously small-world, uniquely Hamiltonian, minimally 3-colorable, maximally outerplanar and perfect. Farey graphs are therefore famous in deterministic models for complex networks. By lacking of the most important characteristics of scale-free, Farey graphs are not a good model for networks associated with some empirical complex systems. We discuss here a category of graphs which are extension of the well-known Farey graphs. These new models are named generalized Farey graphs here. We focus on the analysis of the topological characteristics of the new models and deduce the complicated and graceful analytical results from the growth mechanism used in generalized Farey graphs. The conclusions show that the new models not only possess the properties of being small-world and highly clustered, but also possess the quality of being scale-free. We also find that it is precisely because of the exponential increase of nodes' degrees in generalized Farey graphs as they grow that caused the new networks to have scale-free characteristics. In contrast, the linear incrementation of nodes' degrees in Farey graphs can only cause an exponential degree distribution.
| S-EPMC6093942 | biostudies-literature
Project description:BACKGROUND: Parameter estimation for differential equation models of intracellular processes is a highly relevant bu challenging task. The available experimental data do not usually contain enough information to identify all parameters uniquely, resulting in ill-posed estimation problems with often highly correlated parameters. Sampling-based Bayesian statistical approaches are appropriate for tackling this problem. The samples are typically generated via Markov chain Monte Carlo, however such methods are computationally expensive and their convergence may be slow, especially if there are strong correlations between parameters. Monte Carlo methods based on Euclidean or Riemannian Hamiltonian dynamics have been shown to outperform other samplers by making proposal moves that take the local sensitivities of the system's states into account and accepting these moves with high probability. However, the high computational cost involved with calculating the Hamiltonian trajectories prevents their widespread use for all but the smallest differential equation models. The further development of efficient sampling algorithms is therefore an important step towards improving the statistical analysis of predictive models of intracellular processes. RESULTS: We show how state of the art Hamiltonian Monte Carlo methods may be significantly improved for steady state dynamical models. We present a novel approach for efficiently calculating the required geometric quantities by tracking steady states across the Hamiltonian trajectories using a Newton-Raphson method and employing local sensitivity information. Using our approach, we compare both Euclidean and Riemannian versions of Hamiltonian Monte Carlo on three models for intracellular processes with real data and demonstrate at least an order of magnitude improvement in the effective sampling speed. We further demonstrate the wider applicability of our approach to other gradient based MCMC methods, such as those based on Langevin diffusions. CONCLUSION: Our approach is strictly benefitial in all test cases. The Matlab sources implementing our MCMC methodology is available from https://github.com/a-kramer/ode_rmhmc.
| S-EPMC4262080 | biostudies-other
Project description:This Teaching Resource provides lecture notes, slides, and a problem set that can assist in teaching concepts related to dynamical systems tools for the analysis of ordinary differential equation (ODE)-based models. The concepts are applied to familiar biological problems, and the material is appropriate for graduate students or advanced undergraduates. The lecture explains how equations describing biochemical signaling networks can be derived from diagrams that illustrate the reactions in graphical form. Because such reactions are most frequently described using systems of ODEs, the lecture discusses and illustrates the principles underlying the numerical solution of ODEs. Methods for determining the stability of steady-state solutions of one or two-dimensional ODE systems are covered and illustrated using standard graphical methods. The concept of a bifurcation, a condition at which a system's behavior changes qualitatively, is also introduced. A problem set is included that (i) requires students to implement an ODE model of biochemical reactions using MATLAB and (ii) allows them to explore dynamical systems concepts.
| S-EPMC4118930 | biostudies-literature
Project description:Nonlinear stochastic dynamical systems are widely used to model systems across the sciences and engineering. Such models are natural to formulate and can be analyzed mathematically and numerically. However, difficulties associated with inference from time-series data about unknown parameters in these models have been a constraint on their application. We present a new method that makes maximum likelihood estimation feasible for partially-observed nonlinear stochastic dynamical systems (also known as state-space models) where this was not previously the case. The method is based on a sequence of filtering operations which are shown to converge to a maximum likelihood parameter estimate. We make use of recent advances in nonlinear filtering in the implementation of the algorithm. We apply the method to the study of cholera in Bangladesh. We construct confidence intervals, perform residual analysis, and apply other diagnostics. Our analysis, based upon a model capturing the intrinsic nonlinear dynamics of the system, reveals some effects overlooked by previous studies.
| S-EPMC3020138 | biostudies-literature
Project description:To perform nontrivial, real-time computations on a sensory input stream, biological systems must retain a short-term memory trace of their recent inputs. It has been proposed that generic high-dimensional dynamical systems could retain a memory trace for past inputs in their current state. This raises important questions about the fundamental limits of such memory traces and the properties required of dynamical systems to achieve these limits. We address these issues by applying Fisher information theory to dynamical systems driven by time-dependent signals corrupted by noise. We introduce the Fisher Memory Curve (FMC) as a measure of the signal-to-noise ratio (SNR) embedded in the dynamical state relative to the input SNR. The integrated FMC indicates the total memory capacity. We apply this theory to linear neuronal networks and show that the capacity of networks with normal connectivity matrices is exactly 1 and that of any network of N neurons is, at most, N. A nonnormal network achieving this bound is subject to stringent design constraints: It must have a hidden feedforward architecture that superlinearly amplifies its input for a time of order N, and the input connectivity must optimally match this architecture. The memory capacity of networks subject to saturating nonlinearities is further limited, and cannot exceed square root N. This limit can be realized by feedforward structures with divergent fan out that distributes the signal across neurons, thereby avoiding saturation. We illustrate the generality of the theory by showing that memory in fluid systems can be sustained by transient nonnormal amplification due to convective instability or the onset of turbulence.
| S-EPMC2596211 | biostudies-literature
Project description:The notion of a part of phase space containing desired (or allowed) states of a dynamical system is important in a wide range of complex systems research. It has been called the safe operating space, the viability kernel or the sunny region. In this paper we define the notion of survivability: Given a random initial condition, what is the likelihood that the transient behaviour of a deterministic system does not leave a region of desirable states. We demonstrate the utility of this novel stability measure by considering models from climate science, neuronal networks and power grids. We also show that a semi-analytic lower bound for the survivability of linear systems allows a numerically very efficient survivability analysis in realistic models of power grids. Our numerical and semi-analytic work underlines that the type of stability measured by survivability is not captured by common asymptotic stability measures.
| S-EPMC4942794 | biostudies-literature
Project description:We propose an optical simulation of dissipation-induced correlations in one-dimensional (1D) interacting bosonic systems, using a two-dimensional (2D) array of linear photonic waveguides and only classical light. We show that for the case of two bosons in a 1D lattice, one can simulate on-site two-body dissipative dynamics using a linear 2D waveguide array with lossy diagonal waveguides. The intensity distribution of the propagating light directly maps out the wave function, allowing one to observe the dissipation-induced correlations with simple measurements. Beyond the on-site model, we also show that a generalised model containing nearest-neighbour dissipative interaction can be engineered and probed in the proposed set-up.
| S-EPMC4338417 | biostudies-other
Project description:Enzyme catalysis is a complex process involving several steps along the reaction coordinates, including substrate recognition and binding, chemical transformation, and product release. Evidence continues to emerge linking the functional and evolutionary role of conformational exchange processes in optimal catalytic activity. Ligand binding changes the conformational landscape of enzymes, inducing long-range conformational rearrangements. Using functionally distinct members of the pancreatic ribonuclease superfamily as a model system, we characterized the structural and conformational changes associated with the binding of two mononucleotide ligands. By combining NMR chemical shift titration experiments with the chemical shift projection analysis (CHESPA) and relaxation dispersion experiments, we show that biologically distinct members of the RNase superfamily display discrete chemical shift perturbations upon ligand binding that are not conserved even in structurally related members. Amino acid networks exhibiting coordinated chemical shift displacements upon binding of the two ligands are unique to each of the RNases analyzed. Our results reveal the contribution of conformational rearrangements to the observed chemical shift perturbations. These observations provide important insights into the contribution of the different ligand binding specificities and effects of conformational exchange on the observed perturbations associated with ligand binding for functionally diverse members of the pancreatic RNase superfamily.
| S-EPMC6005897 | biostudies-literature
Project description:When simulating molecular systems using deterministic equations of motion (e.g., Newtonian dynamics), such equations are generally numerically integrated according to a well-developed set of algorithms that share commonly agreed-upon desirable properties. However, for stochastic equations of motion (e.g., Langevin dynamics), there is still broad disagreement over which integration algorithms are most appropriate. While multiple desiderata have been proposed throughout the literature, consensus on which criteria are important is absent, and no published integration scheme satisfies all desiderata simultaneously. Additional nontrivial complications stem from simulating systems driven out of equilibrium using existing stochastic integration schemes in conjunction with recently developed nonequilibrium fluctuation theorems. Here, we examine a family of discrete time integration schemes for Langevin dynamics, assessing how each member satisfies a variety of desiderata that have been enumerated in prior efforts to construct suitable Langevin integrators. We show that the incorporation of a novel time step rescaling in the deterministic updates of position and velocity can correct a number of dynamical defects in these integrators. Finally, we identify a particular splitting (related to the velocity Verlet discretization) that has essentially universally appropriate properties for the simulation of Langevin dynamics for molecular systems in equilibrium, nonequilibrium, and path sampling contexts.
| S-EPMC4065221 | biostudies-literature
Project description:The effect of p53-dependent cell-cycle arrest and senescence on Emu-myc-induced B-cell lymphoma development remains controversial. To address this question, we crossed Emu-myc mice with the p53(515C) mutant mouse, encoding the mutant p53R172P protein that retains the ability to activate the cell-cycle inhibitor and senescence activator p21. Importantly, this mutant lacks the ability to activate p53-dependent apoptotic genes. Hence, Emu-myc mice that harbor two p53(515C) alleles are completely defective for p53-dependent apoptosis. Both Emu-myc::p53(515C/515C) and Emu-myc::p53(515C/+) mice survive significantly longer than Emu-myc::p53(+/-) mice, indicating the importance of the p53-dependent non-apoptotic pathways in B-cell lymphomagenesis. In addition, the p53(515C) allele is deleted in several Emu-myc::p53(515C/+) lymphomas, further emphasizing the functionality of p53R172P in tumor inhibition. Lymphomas from both Emu-myc::p53(515C/515C) and Emu-myc::p53(515C/+) mice retain the ability to upregulate p21, resulting in cellular senescence. Senescence-associated beta-galactosidase (SA beta-gal) activity was observed in lymphomas from Emu-myc::p53(+/+), Emu-myc::p53(515C/515C) and Emu-myc::p53(515C /+) mice but not in lymphomas isolated from Emu-myc::p53(+/-) mice. Thus, in the absence of p53-dependent apoptosis, the ability of p53R172P to induce senescence leads to a significant delay in B-cell lymphoma development.
| S-EPMC2903442 | biostudies-literature