Stability and steady state of complex cooperative systems: a diakoptic approach.
Ontology highlight
ABSTRACT:
SUBMITTER: Greulich P
PROVIDER: S-EPMC6936286 | biostudies-literature | 2019 Dec
REPOSITORIES: biostudies-literature
Publications
Stability and steady state of complex cooperative systems: a diakoptic approach.
Greulich Philip P MacArthur Ben D BD Parigini Cristina C Sánchez-García Rubén J RJ
Royal Society open science 20191204 12
Cooperative dynamics are common in ecology and population dynamics. However, their commonly high degree of complexity with a large number of coupled degrees of freedom renders them difficult to analyse. Here, we present a graph-theoretical criterion, via a diakoptic approach (divide-and-conquer) to determine a cooperative system's stability by decomposing the system's dependence graph into its strongly connected components (SCCs). In particular, we show that a linear cooperative system is Lyapun ...[more]
Similar Datasets
Project description:In pharmacology and systems biology, it is a fundamental problem to determine how biological systems change their dose-response behavior upon perturbations. In particular, it is unclear how topologies, reactions, and parameters (differentially) affect the dose response. Because parameters are often unknown, systematic approaches should directly relate network structure and function. Here, we outline a procedure to compare general non-monotone dose-response curves and subsequently develop a comprehensive theory for differential dose responses of biochemical networks captured by non-equilibrium steady-state linear framework models. Although these models are amenable to analytical derivations of non-equilibrium steady states in principle, their size frequently increases (super) exponentially with model size. We extract general principles of differential responses based on a model's graph structure and thereby alleviate the combinatorial explosion. This allows us, for example, to determine reactions that affect differential responses, to identify classes of networks with equivalent differential, and to reject hypothetical models reliably without needing to know parameter values. We exemplify such applications for models of insulin signaling.
| S-EPMC5985043 | biostudies-literature
Project description:The â??small perturbationâ?? approach is critical in studying the â??steady stateâ?? of a biological system. In our experiments, small perturbations were generated by applying a series of repeating intermittent small doses of ultraviolet radiation to a human keratinocyte cell line, HaCaT. The biological response was assessed by monitoring the gene expression profiles using a high reliability and high resolution cDNA microarray system. Following intermittent 10 J/m2 UVB small perturbations, two opposite classes of genes, down-regulated and up-regulated, exhibited an immediate response followed by relaxation between each small perturbation, but were prolonged down- or up-regulated without relaxation while larger doses (233 or 582.5 J/m2) of UVB were applied. A repeated cycle pattern of gene expression following small perturbations is an indication of the existence of steady states. This cycle pattern is suppressed when large perturbations are applied. We believe that this is a universal phenomenon. In our experiments, the functions of up-regulated genes were mainly associated with anti-proliferation, anti-mitogenesis, and apoptosis. On the other hand, down-regulated genes were mainly related to proliferation, mitogenesis, and anti-apoptosis. In conclusion, this study experimentally proves the concept of steady state at the transcription level and demonstrates the feasibility of using small perturbation approaches for investigating steady states. This study could also set a foundation of computational systems biology, which has implicitly used the concept of steady state. Keywords: time course Three UVB exposures are indicated by UV at time points of 0, 8 and 16 hours. T1, T2, T3, T4, T5, and T6 denote the sampling time points. T1, T3 and T5 are allocated 30 minutes after the corresponding UVB irradiation. T2, T4 and T6 are allocated 8 hours after each UVB irradiation. At each sampling time point, two samples (control and UV-irradiated) are collected. See supplementary file Loop_design.pdf for further explanation.
2010-06-24 | E-GEOD-7060 | biostudies-arrayexpress
Project description:The ‘small perturbation’ approach is critical in studying the ‘steady state’ of a biological system. In our experiments, small perturbations were generated by applying a series of repeating intermittent small doses of ultraviolet radiation to a human keratinocyte cell line, HaCaT. The biological response was assessed by monitoring the gene expression profiles using a high reliability and high resolution cDNA microarray system. Following intermittent 10 J/m2 UVB small perturbations, two opposite classes of genes, down-regulated and up-regulated, exhibited an immediate response followed by relaxation between each small perturbation, but were prolonged down- or up-regulated without relaxation while larger doses (233 or 582.5 J/m2) of UVB were applied. A repeated cycle pattern of gene expression following small perturbations is an indication of the existence of steady states. This cycle pattern is suppressed when large perturbations are applied. We believe that this is a universal phenomenon. In our experiments, the functions of up-regulated genes were mainly associated with anti-proliferation, anti-mitogenesis, and apoptosis. On the other hand, down-regulated genes were mainly related to proliferation, mitogenesis, and anti-apoptosis. In conclusion, this study experimentally proves the concept of steady state at the transcription level and demonstrates the feasibility of using small perturbation approaches for investigating steady states. This study could also set a foundation of computational systems biology, which has implicitly used the concept of steady state. Keywords: time course
2009-01-30 | GSE7060 | GEO
Project description:Understanding how cells change their identity and behaviour in living systems is an important question in many fields of biology. The problem of inferring cell trajectories from single-cell measurements has been a major topic in the single-cell analysis community, with different methods developed for equilibrium and non-equilibrium systems (e.g. haematopoeisis vs. embryonic development). We show that optimal transport analysis, a technique originally designed for analysing time-courses, may also be applied to infer cellular trajectories from a single snapshot of a population in equilibrium. Therefore, optimal transport provides a unified approach to inferring trajectories that is applicable to both stationary and non-stationary systems. Our method, StationaryOT, is mathematically motivated in a natural way from the hypothesis of a Waddington's epigenetic landscape. We implement StationaryOT as a software package and demonstrate its efficacy in applications to simulated data as well as single-cell data from Arabidopsis thaliana root development.
| S-EPMC8691649 | biostudies-literature
Project description:Hundreds of millions of smallholders in emerging countries substantially overuse nitrogen (N) fertilizers, driving local environmental pollution and global climate change. Despite local demonstration-scale successes, widespread mobilization of smallholders to adopt precise N management practices remains a challenge, largely due to associated high costs and complicated sampling and calculations. Here, we propose a long-term steady-state N balance (SSNB) approach without these complications that is suitable for sustainable smallholder farming. The hypothesis underpinning the concept of SSNB is that an intensively cultivated soil-crop system with excessive N inputs and high N losses can be transformed into a steady-state system with minimal losses while maintaining high yields. Based on SSNB, we estimate the optimized N application range across 3,824 crop counties for the three staple crops in China. We evaluated SSNB first in ca. 18,000 researcher-managed on-farm trials followed by testing in on-farm trials with 13,760 smallholders who applied SSNB-optimized N rates under the guidance of local extension staff. Results showed that SSNB could significantly reduce N fertilizer use by 21 to 28% while maintaining or increasing yields by 6 to 7%, compared to current smallholder practices. The SSNB approach could become an effective tool contributing to the global N sustainability of smallholder agriculture.
| S-EPMC8488683 | biostudies-literature
Project description:Empirically based rate-and-state friction laws (RSFLs) have been proposed to model the dependence of friction forces with slip and time. The relevance of the RSFL for earthquake mechanics is that few constitutive parameters define critical conditions for fault stability (i.e., critical stiffness and frictional fault behavior). However, the RSFLs were determined from experiments conducted at subseismic slip rates (V < 1 cm/s), and their extrapolation to earthquake deformation conditions (V > 0.1 m/s) remains questionable on the basis of the experimental evidence of (1) large dynamic weakening and (2) activation of particular fault lubrication processes at seismic slip rates. Here we propose a modified RSFL (MFL) based on the review of a large published and unpublished data set of rock friction experiments performed with different testing machines. The MFL, valid at steady state conditions from subseismic to seismic slip rates (0.1 µm/s < V < 3 m/s), describes the initiation of a substantial velocity weakening in the 1-20 cm/s range resulting in a critical stiffness increase that creates a peak of potential instability in that velocity regime. The MFL leads to a new definition of fault frictional stability with implications for slip event styles and relevance for models of seismic rupture nucleation, propagation, and arrest.
| S-EPMC5021208 | biostudies-literature
Project description:Under stress, cells need to optimize the activity of a wide range of gene products during the response phases: shock, adaptation, and recovery. This requires coordination of several levels of regulation, including turnover and translation efficiencies of mRNAs. Mitogen-activated protein (MAP) kinase pathways are implicated in many aspects of the environmental stress response, including initiation of transcription, translation efficiency, and mRNA turnover. In this study, we analyze mRNA turnover rates and mRNA steady-state levels at different time points following mild hyperosmotic shock in Saccharomyces cerevisiae cells. The regulation of mRNA stability is transient and affects most genes for which there is a change in transcript level. These changes precede and prepare for the changes in steady-state levels, both regarding the initial increase and the later decline of stress-induced mRNAs. The inverse is true for stress-repressed genes, which become stabilized during hyperosmotic stress in preparation of an increase as the cells recover. The MAP kinase Hog1 affects both steady-state levels and stability of stress-responsive transcripts, whereas the Hog1-activated kinase Rck2 influences steady-state levels without a major effect on stability. Regulation of mRNA stability is a wide-spread, but not universal, effect on stress-responsive transcripts during transient hyperosmotic stress. By destabilizing stress-induced mRNAs when their steady-state levels have reached a maximum, the cell prepares for the subsequent recovery phase when these transcripts are to return to normal levels. Conversely, stabilization of stress-repressed mRNAs permits their rapid accumulation in the recovery phase. Our results show that mRNA turnover is coordinated with transcriptional induction.
| S-EPMC2661839 | biostudies-literature
Project description:Stochastic population models are widely used to model phenomena in different areas such as cyber-physical systems, chemical kinetics, collective animal behaviour, and beyond. Quantitative analysis of stochastic population models easily becomes challenging due to the combinatorial number of possible states of the population. Moreover, while the modeller easily hypothesises the mechanistic aspects of the model, the quantitative parameters associated to these mechanistic transitions are difficult or impossible to measure directly. In this paper, we investigate how formal verification methods can aid parameter inference for population discrete-time Markov chains in a scenario where only a limited sample of population-level data measurements-sample distributions among terminal states-are available. We first discuss the parameter identifiability and uncertainty quantification in this setup, as well as how the existing techniques of formal parameter synthesis and Bayesian inference apply. Then, we propose and implement four different methods, three of which incorporate formal parameter synthesis as a pre-computation step. We empirically evaluate the performance of the proposed methods over four representative case studies. We find that our proposed methods incorporating formal parameter synthesis as a pre-computation step allow us to significantly enhance the accuracy, precision, and scalability of inference. Specifically, in the case of unidentifiable parameters, we accurately capture the subspace of parameters which is data-compliant at a desired confidence level.
| S-EPMC10642793 | biostudies-literature
Project description:Microphysiological systems (MPS), consisting of tissue constructs, biomaterials, and culture media, aim to recapitulate relevant organ functions in vitro. MPS components are housed in fluidic hardware with operational protocols, such as periodic complete media replacement. Such batch-like operations provide relevant nutrients and remove waste products but also reset cell-secreted mediators (e.g. cytokines, hormones) and potentially limit exposure to drugs (and metabolites). While each component plays an essential role for tissue functionality, MPS-specific nutrient needs are not yet well-characterized nor utilized to operate MPSs at more physiologically-relevant conditions. MPS-specific nutrient needs for gut (immortalized cancer cells), liver (human primary hepatocytes) and cardiac (iPSC-derived cardiomyocytes) MPSs were experimentally quantified. In a long-term study of the gut MPS (10 days), this knowledge was used to design operational protocols to maintain glucose and lactate at desired levels. This quasi-steady state operation was experimentally validated by monitoring glucose and lactate as well as MPS functionality. In a theoretical study, nutrient needs of an integrated multi-MPS platform (gut, liver, cardiac MPSs) were computationally simulated to identify long-term quasi-steady state operations. This integrative experimental and computational approach demonstrates the utilization of quantitative multi-scale characterization of MPSs and incorporating MPS-specific information to establish more physiologically-relevant experimental operations.
| S-EPMC5964119 | biostudies-literature
Project description:Mycobacterium tuberculosis and many other members of the Actinomycetes family produce mycothiol, i.e., 1-d-myo-inosityl-2-(N-acetyl-l-cysteinyl)amido-2-deoxy-alpha-d-glucopyranoside (MSH or AcCys-GlcN-Ins), to act against oxidative and antibiotic stress. The biosynthesis of MSH is essential for cell growth and has been proposed to proceed via a biosynthetic pathway involving four key enzymes, MshA-MshD. The MSH biosynthetic enzymes present potential targets for inhibitor design. With this as a long-term goal, we have carried out a kinetic and mechanistic characterization, using steady-state and pre-steady-state approaches, of the recombinant Mycobacterium smegmatis MshC. MshC catalyzes the ATP-dependent condensation of GlcN-Ins and cysteine to form Cys-GlcN-Ins. Initial velocity and inhibition studies show that the steady-state kinetic mechanism of MshC is a Bi Uni Uni Bi Ping Pong mechanism, with ATP binding followed by cysteine binding, release of PPi, binding of GlcN-Ins, followed by the release of Cys-GlcN-Ins and AMP. The steady-state kinetic parameters were determined to be kcat equal to 3.15 s-1, and Km values of 1.8, 0.1, and 0.16 mM for ATP, cysteine, and GlcN-Ins, respectively. A stable bisubstrate analogue, 5'-O-[N-(l-cysteinyl)sulfamonyl]adenosine, exhibits competitive inhibition versus ATP and noncompetitive inhibition versus cysteine, with an inhibition constant of approximately 306 nM versus ATP. Single-turnover reactions of the first and second half reactions were determined using rapid-quench techniques, giving rates of approximately 9.4 and approximately 5.2 s-1, respectively, consistent with the cysteinyl adenylate being a kinetically competent intermediate in the reaction by MshC.
| S-EPMC2526253 | biostudies-literature