Consistent robustness analysis (CRA) identifies biologically relevant properties of regulatory network models.
Ontology highlight
ABSTRACT: Background
A number of studies have previously demonstrated that "goodness of fit" is insufficient in reliably classifying the credibility of a biological model. Robustness and/or sensitivity analysis is commonly employed as a secondary method for evaluating the suitability of a particular model. The results of such analyses invariably depend on the particular parameter set tested, yet many parameter values for biological models are uncertain.Results
Here, we propose a novel robustness analysis that aims to determine the "common robustness" of the model with multiple, biologically plausible parameter sets, rather than the local robustness for a particular parameter set. Our method is applied to two published models of the Arabidopsis circadian clock (the one-loop [1] and two-loop [2] models). The results reinforce current findings suggesting the greater reliability of the two-loop model and pinpoint the crucial role of TOC1 in the circadian network.Conclusions
Consistent Robustness Analysis can indicate both the relative plausibility of different models and also the critical components and processes controlling each model.
SUBMITTER: Saithong T
PROVIDER: S-EPMC3002950 | biostudies-literature |
REPOSITORIES: biostudies-literature
Similar Datasets
Project description:BackgroundIn cancer research, robustness of a complex biochemical network is one of the most relevant properties to investigate for the development of novel targeted therapies. In cancer systems biology, biological networks are typically modeled through Ordinary Differential Equation (ODE) models. Hence, robustness analysis consists in quantifying how much the temporal behavior of a specific node is influenced by the perturbation of model parameters. The Conditional Robustness Algorithm (CRA) is a valuable methodology to perform robustness analysis on a selected output variable, representative of the proliferation activity of cancer disease.ResultsHere we introduce our new freely downloadable software, the CRA Toolbox. The CRA Toolbox is an Object-Oriented MATLAB package which implements the features of CRA for ODE models. It offers the users the ability to import a mathematical model in Systems Biology Markup Language (SBML), to perturb the model parameter space and to choose the reference node for the robustness analysis. The CRA Toolbox allows the users to visualize and save all the generated results through a user-friendly Graphical User Interface (GUI). The CRA Toolbox has a modular and flexible architecture since it is designed according to some engineering design patterns. This tool has been successfully applied in three nonlinear ODE models: the Prostate-specific Pten-/- mouse model, the Pulse Generator Network and the EGFR-IGF1R pathway.ConclusionsThe CRA Toolbox for MATLAB is an open-source tool implementing the CRA to perform conditional robustness analysis. With its unique set of functions, the CRA Toolbox is a remarkable software for the topological study of biological networks. The source and example code and the corresponding documentation are freely available at the web site: http://gitlab.ict4life.com/SysBiOThe/CRA-Matlab .
| S-EPMC6617887 | biostudies-literature
Project description:Gene regulatory networks are powerful tools which facilitate hypothesis generation and candidate gene discovery. However, the extent to which the network predictions are biologically relevant is often unclear. Recently a GENIE3 network which predicted targets of wheat transcription factors was produced. Here we used an independent RNA-Seq dataset to test the predictions of the wheat GENIE3 network for the senescence-regulating transcription factor NAM-A1 (TraesCS6A02G108300). We re-analyzed the RNA-Seq data against the RefSeqv1.0 genome and identified a set of differentially expressed genes (DEGs) between the wild-type and nam-a1 mutant which recapitulated the known role of NAM-A1 in senescence and nutrient remobilisation. We found that the GENIE3-predicted target genes of NAM-A1 overlap significantly with the DEGs, more than would be expected by chance. Based on high levels of overlap between GENIE3-predicted target genes and the DEGs, we identified candidate senescence regulators. We then explored genome-wide trends in the network related to polyploidy and found that only homeologous transcription factors are likely to share predicted targets in common. However, homeologs which vary in expression levels across tissues are less likely to share predicted targets than those that do not, suggesting that they may be more likely to act in distinct pathways. This work demonstrates that the wheat GENIE3 network can provide biologically-relevant predictions of transcription factor targets, which can be used for candidate gene prediction and for global analyses of transcription factor function. The GENIE3 network has now been integrated into the KnetMiner web application, facilitating its use in future studies.
| S-EPMC7534433 | biostudies-literature
Project description:Pseudomonas aeruginosa (P. aeruginosa) is an opportunistic pathogen and the leading cause of infection in patients with cystic fibrosis (CF). The ability of P. aeruginosa to evade host responses and develop into chronic infection causes significant morbidity and mortality. Several mouse models have been developed to study chronic respiratory infections induced by P. aeruginosa, with the bead agar model being the most widely used. However, this model has several limitations, including the requirement for surgical procedures and high mortality rates. Herein, we describe novel and adapted biologically relevant models of chronic lung infection caused by P. aeruginosa. Three methods are described: a clinical isolate infection model, utilising isolates obtained from patients with CF; an incomplete antibiotic clearance model, leading to bacterial bounce-back; and the establishment of chronic infection; and an adapted water bottle chronic infection model. These models circumvent the requirement for a surgical procedure and, importantly, can be induced with clinical isolates of P. aeruginosa and in wild-type mice. We also demonstrate successful induction of chronic infection in the transgenic βENaC murine model of CF. We envisage that the models described will facilitate the investigations of host and microbial factors, and the efficacy of novel antimicrobials, during chronic P. aeruginosa respiratory infections.
| S-EPMC10458525 | biostudies-literature
Project description:Genome-scale metabolic models comprehensively describe an organism's metabolism and can be tailored using omics data to model condition-specific physiology. The quality of context-specific models is impacted by (i) choice of algorithm and parameters and (ii) alternate context-specific models that equally explain the -omics data. Here we quantify the influence of alternate optima on microbial and mammalian model extraction using GIMME, iMAT, MBA, and mCADRE. We find that metabolic tasks defining an organism's phenotype must be explicitly and quantitatively protected. The scope of alternate models is strongly influenced by algorithm choice and the topological properties of the parent genome-scale model with fatty acid metabolism and intracellular metabolite transport contributing much to alternate solutions in all models. mCADRE extracted the most reproducible context-specific models and models generated using MBA had the most alternate solutions. There were fewer qualitatively different solutions generated by GIMME in E. coli, but these increased substantially in the mammalian models. Screening ensembles using a receiver operating characteristic plot identified the best-performing models. A comprehensive evaluation of models extracted using combinations of extraction methods and expression thresholds revealed that GIMME generated the best-performing models in E. coli, whereas mCADRE is better suited for complex mammalian models. These findings suggest guidelines for benchmarking -omics integration algorithms and motivate the development of a systematic workflow to enumerate alternate models and extract biologically relevant context-specific models.
| S-EPMC10258867 | biostudies-literature
Project description:Liebau pump is a tubular, non-peristaltic, pulsatile pump capable of creating unidirectional flow in the absence of valves. It requires asymmetrical positioning of the pincher relative to the attachment sites of its elastic segment to the rest of the circuit. Biological feasibility of such valveless pumps remains a hotly debated topic. To test the feasibility of the Liebau-based pumping in vessels with biologically relevant properties we quantified the output of Liebau pumps with their compliant segments made of a silicone rubber that mimicked the Young modulus of soft tissues. The lengths, the inner diameters, thicknesses of the tested compliant segments ranged from 1 to 5 cm, 3 to 8 mm and 0.3 to 1 mm, respectively. The compliant segment of the setup was compressed at 0.5-2.5 Hz frequencies using a 3.5-mm-wide rectangular piston. A nearest-neighbor tracking algorithm was used to track movements of 0.5-mm carbon particles within the system. The viscosity of the aqueous solution was varied by increased percentage of glycerin. Measurements yielded quantitative relationships between viscosity, frequency of compression and the net flowrate. The use of the Liebau principle of valveless pumping in conjunction with physiologically sized vessel and contraction frequencies yields flowrates comparable to peristaltic pumps of the same dimensions. We conclude that the data confirm physiological feasibility of Liebau-based pumping and warrant further testing of its mechanism using excised biological conduits or tissue engineered components. Such biomimetic pumps can serve as energy-efficient flow generators in microdevices or to study the function of embryonic heart during its normal development or in diseased states.
| S-EPMC8169938 | biostudies-literature
Project description:Gene regulatory networks are useful to understand the activity behind the complex mechanisms in transcriptional regulation. A main goal in contemporary biology is using such networks to understand the systemic regulation of gene expression. In this work, we carried out a systematic study of a transcriptional regulatory network derived from a comprehensive selection of all potential transcription factor interactions downstream from MEF2C, a human transcription factor master regulator. By analyzing the connectivity structure of such network, we were able to find different biologically functional processes and specific biochemical pathways statistically enriched in communities of genes into the network, such processes are related to cell signaling, cell cycle and metabolism. In this way we further support the hypothesis that structural properties of biological networks encode an important part of their functional behavior in eukaryotic cells.
| S-EPMC4878384 | biostudies-literature
Project description:BackgroundBecause biological networks exhibit a high-degree of robustness, a systemic understanding of their architecture and function requires an appraisal of the network design principles that confer robustness. In this project, we conduct a computational study of the contribution of three degree-based topological properties (transcription factor-target ratio, degree distribution, cross-talk suppression) and their combinations on the robustness of transcriptional regulatory networks. We seek to quantify the relative degree of robustness conferred by each property (and combination) and also to determine the extent to which these properties alone can explain the robustness observed in transcriptional networks.ResultsTo study individual properties and their combinations, we generated synthetic, random networks that retained one or more of the three properties with values derived from either the yeast or E. coli gene regulatory networks. Robustness of these networks were estimated through simulation. Our results indicate that the combination of the three properties we considered explains the majority of the structural robustness observed in the real transcriptional networks. Surprisingly, scale-free degree distribution is, overall, a minor contributor to robustness. Instead, most robustness is gained through topological features that limit the complexity of the overall network and increase the transcription factor subnetwork sparsity.ConclusionsOur work demonstrates that (i) different types of robustness are implemented by different topological aspects of the network and (ii) size and sparsity of the transcription factor subnetwork play an important role for robustness induction. Our results are conserved across yeast and E Coli, which suggests that the design principles examined are present within an array of living systems.
| S-EPMC3541983 | biostudies-literature
Project description:Topology affects physical and biological properties of DNA and impacts fundamental cellular processes, such as gene expression, genome replication, chromosome structure and segregation. In all organisms DNA topology is carefully modulated and the supercoiling degree of defined genome regions may change according to physiological and environmental conditions. Elucidation of structural properties of DNA molecules with different topology may thus help to better understand genome functions. Whereas a number of structural studies have been published on highly negatively supercoiled DNA molecules, only preliminary observations of highly positively supercoiled are available, and a description of DNA structural properties over the full range of supercoiling degree is lacking. Atomic Force Microscopy (AFM) is a powerful tool to study DNA structure at single molecule level. We here report a comprehensive analysis by AFM of DNA plasmid molecules with defined supercoiling degree, covering the full spectrum of biologically relevant topologies, under different observation conditions. Our data, supported by statistical and biochemical analyses, revealed striking differences in the behavior of positive and negative plasmid molecules.
| S-EPMC5906655 | biostudies-literature
Project description:Frailty is a late-life syndrome of vulnerability to adverse health outcomes characterized by a phenotype that includes muscle weakness, fatigue, and inflammatory pathway activation. The identification of biologically relevant pathways that influence frailty is challenged by its biological complexity and the necessity in separating disease states from the syndrome of frailty. As with longevity research, genetic analyses may help to provide insights into biologically relevant pathways that contribute to frailty.Based on current understanding of the physiological basis of frailty, we hypothesize that variation in genes related to inflammation and muscle maintenance would associate with frailty. One thousand three hundred and fifty-four single-nucleotide polymorphisms were genotyped across 134 candidate genes using the Illumina Genotyping platform, and the rank order by strength of association between frailty and genotype was determined in a cross-sectional study.Although no single-nucleotide polymorphism reached study-wide significance after controlling family-wise false-discovery rate at 0.05, single-nucleotide polymorphisms within the 5-methyltetrahydrofolate-homocysteine methyltransferase (MTR), Caspase 8 (CASP8), CREB-binding protein (CREBBP), lysine acetyltransferase 2B (KAT2B), and beta-transducin repeat containing (BTRC) loci were among those strongly associated with frailty.The apoptosis- and transcription regulation-related pathways highlighted by this preliminary analysis were consistent with prior gene expression studies in a frail mouse model and provide useful etiological insights for future biological studies of frailty.
| S-EPMC3156628 | biostudies-literature
Project description:BackgroundThe use of 3-D similarity techniques in the analysis of biological data and virtual screening is pervasive, but what is a biologically meaningful 3-D similarity value? Can one find statistically significant separation between "active/active" and "active/inactive" spaces? These questions are explored using 734,486 biologically tested chemical structures, 1,389 biological assay data sets, and six different 3-D similarity types utilized by PubChem analysis tools.ResultsThe similarity value distributions of 269.7 billion unique conformer pairs from 734,486 biologically tested compounds (all-against-all) from PubChem were utilized to help work towards an answer to the question: what is a biologically meaningful 3-D similarity score? The average and standard deviation for the six similarity measures STST-opt, CTST-opt, ComboTST-opt, STCT-opt, CTCT-opt, and ComboTCT-opt were 0.54 ± 0.10, 0.07 ± 0.05, 0.62 ± 0.13, 0.41 ± 0.11, 0.18 ± 0.06, and 0.59 ± 0.14, respectively. Considering that this random distribution of biologically tested compounds was constructed using a single theoretical conformer per compound (the "default" conformer provided by PubChem), further study may be necessary using multiple diverse conformers per compound; however, given the breadth of the compound set, the single conformer per compound results may still apply to the case of multi-conformer per compound 3-D similarity value distributions. As such, this work is a critical step, covering a very wide corpus of chemical structures and biological assays, creating a statistical framework to build upon.The second part of this study explored the question of whether it was possible to realize a statistically meaningful 3-D similarity value separation between reputed biological assay "inactives" and "actives". Using the terminology of noninactive-noninactive (NN) pairs and the noninactive-inactive (NI) pairs to represent comparison of the "active/active" and "active/inactive" spaces, respectively, each of the 1,389 biological assays was examined by their 3-D similarity score differences between the NN and NI pairs and analyzed across all assays and by assay category types. While a consistent trend of separation was observed, this result was not statistically unambiguous after considering the respective standard deviations. While not all "actives" in a biological assay are amenable to this type of analysis, e.g., due to different mechanisms of action or binding configurations, the ambiguous separation may also be due to employing a single conformer per compound in this study. With that said, there were a subset of biological assays where a clear separation between the NN and NI pairs found. In addition, use of combo Tanimoto (ComboT) alone, independent of superposition optimization type, appears to be the most efficient 3-D score type in identifying these cases.ConclusionThis study provides a statistical guideline for analyzing biological assay data in terms of 3-D similarity and PubChem structure-activity analysis tools. When using a single conformer per compound, a relatively small number of assays appear to be able to separate "active/active" space from "active/inactive" space.
| S-EPMC3223603 | biostudies-literature