Shape, Size, and Robustness: Feasible Regions in the Parameter Space of Biochemical Networks

The concept of robustness of regulatory networks has received much attention in the last decade. One measure of robustness has been associated with the volume of the feasible region, namely, the region in the parameter space in which the system is functional. In this paper, we show that, in addition to volume, the geometry of this region has important consequences for the robustness and the fragility of a network. We develop an approximation within which we could algebraically specify the feasible region. We analyze the segment polarity gene network to illustrate our approach. The study of random walks in the parameter space and how they exit the feasible region provide us with a rich perspective on the different modes of failure of this network model. In particular, we found that, between two alternative ways of activating Wingless, one is more robust than the other. Our method provides a more complete measure of robustness to parameter variation. As a general modeling strategy, our approach is an interesting alternative to Boolean representation of biochemical networks

Metabolomics demonstrates divergent responses of two Eucalyptus species to water stress

Past studies of water stress in Eucalyptus spp. generally highlighted the role of fewer than five “important” metabolites, whereas recent metabolomic studies on other genera have shown tens of compounds are affected. There are currently no metabolite profiling data for responses of stress-tolerant species to water stress. We used GC–MS metabolite profiling to examine the response of leaf metabolites to a long (2 month) and severe (Ψpredawn < −2 MPa) water stress in two species of the perennial tree genus Eucalyptus (the mesic Eucalyptus pauciflora and the semi-arid Eucalyptus dumosa). Polar metabolites in leaves were analysed by GC–MS and inorganic ions by capillary electrophoresis. Pressure–volume curves and metabolite measurements showed that water stress led to more negative osmotic potential and increased total osmotically active solutes in leaves of both species. Water stress affected around 30–40% of measured metabolites in E. dumosa and 10–15% in E. pauciflora. There were many metabolites that were affected in E. dumosa but not E. pauciflora, and some that had opposite responses in the two species. For example, in E. dumosa there were increases in five acyclic sugar alcohols and four low-abundance carbohydrates that were unaffected by water stress in E. pauciflora. Re-watering increased osmotic potential and decreased total osmotically active solutes in E. pauciflora, whereas in E. dumosa re-watering led to further decreases in osmotic potential and increases in total osmotically active solutes. This experiment has added several extra dimensions to previous targeted analyses of water stress responses in Eucalyptus, and highlights that even species that are closely related (e.g. congeners) may respond differently to water stress and re-waterin

Robustness Analysis and Behavior Discrimination in Enzymatic Reaction Networks

Characterizing the behavior and robustness of enzymatic networks with numerous variables and unknown parameter values is a major challenge in biology, especially when some enzymes have counter-intuitive properties or switch-like behavior between activation and inhibition. In this paper, we propose new methodological and tool-supported contributions, based on the intuitive formalism of temporal logic, to express in a rigorous manner arbitrarily complex dynamical properties. Our multi-step analysis allows efficient sampling of the parameter space in order to define feasible regions in which the model exhibits imposed or experimentally observed behaviors. In a first step, an algorithmic methodology involving sensitivity analysis is conducted to determine bifurcation thresholds for a limited number of model parameters or initial conditions. In a second step, this boundary detection is supplemented by a global robustness analysis, based on quasi-Monte Carlo approach that takes into account all model parameters. We apply this method to a well-documented enzymatic reaction network describing collagen proteolysis by matrix metalloproteinase MMP2 and membrane type 1 metalloproteinase (MT1-MMP) in the presence of tissue inhibitor of metalloproteinase TIMP2. For this model, our method provides an extended analysis and quantification of network robustness toward paradoxical TIMP2 switching activity between activation or inhibition of MMP2 production. Further implication of our approach is illustrated by demonstrating and analyzing the possible existence of oscillatory behaviors when considering an extended open configuration of the enzymatic network. Notably, we construct bifurcation diagrams that specify key parameters values controlling the co-existence of stable steady and non-steady oscillatory proteolytic dynamics

Critical Dynamics in Genetic Regulatory Networks: Examples from Four Kingdoms

The coordinated expression of the different genes in an organism is essential to sustain functionality under the random external perturbations to which the organism might be subjected. To cope with such external variability, the global dynamics of the genetic network must possess two central properties. (a) It must be robust enough as to guarantee stability under a broad range of external conditions, and (b) it must be flexible enough to recognize and integrate specific external signals that may help the organism to change and adapt to different environments. This compromise between robustness and adaptability has been observed in dynamical systems operating at the brink of a phase transition between order and chaos. Such systems are termed critical. Thus, criticality, a precise, measurable, and well characterized property of dynamical systems, makes it possible for robustness and adaptability to coexist in living organisms. In this work we investigate the dynamical properties of the gene transcription networks reported for S. cerevisiae, E. coli, and B. subtilis, as well as the network of segment polarity genes of D. melanogaster, and the network of flower development of A. thaliana. We use hundreds of microarray experiments to infer the nature of the regulatory interactions among genes, and implement these data into the Boolean models of the genetic networks. Our results show that, to the best of the current experimental data available, the five networks under study indeed operate close to criticality. The generality of this result suggests that criticality at the genetic level might constitute a fundamental evolutionary mechanism that generates the great diversity of dynamically robust living forms that we observe around us

Using serological measures to monitor changes in malaria transmission in Vanuatu

BACKGROUND: With renewed interest in malaria elimination, island environments present unique opportunities to achieve this goal. However, as transmission decreases, monitoring and evaluation programmes need increasingly sensitive tools to assess Plasmodium falciparum and Plasmodium vivax exposure. In 2009, to assess the role of serological markers in evaluating malaria transmission, a cross-sectional seroprevalence study was carried out in Tanna and Aneityum, two of the southernmost islands of the Vanuatu archipelago, areas where malaria transmission has been variably reduced over the past few decades. METHODS: Malaria transmission was assessed using serological markers for exposure to P. falciparum and P. vivax. Filter blood spot papers were collected from 1,249 people from Tanna, and 517 people from Aneityum to assess the prevalence of antibodies to two P. falciparum antigens (MSP-119 and AMA-1) and two P. vivax antigens (MSP-119 and AMA-1). Age-specific prevalence was modelled using a simple catalytic conversion model based on maximum likelihood to generate a community seroconversion rate (SCR). RESULTS: Overall seropositivity in Tanna was 9.4%, 12.4% and 16.6% to P. falciparum MSP-119, AMA-1 and Schizont Extract respectively and 12.6% and 15.0% to P. vivax MSP-119 and AMA-1 respectively. Serological results distinguished between areas of differential dominance of either P. vivax or P. falciparum and analysis of age-stratified results showed a step in seroprevalence occurring approximately 30 years ago on both islands, indicative of a change in transmission intensity at this time. Results from Aneityum suggest that several children may have been exposed to malaria since the 2002 P. vivax epidemic. CONCLUSION: Seroepidemiology can provide key information on malaria transmission for control programmes, when parasite rates are low. As Vanuatu moves closer to malaria elimination, monitoring changes in transmission intensity and identification of residual malaria foci is paramount in order to concentrate intervention efforts

Characterization of a recent malaria outbreak in the autonomous indigenous region of Guna Yala, Panama

BackgroundThis study aims to describe the epidemiological and entomological factors associated with a recent malaria outbreak that occurred in 2012 in a socially marginalized population from Guna Yala Comarca in Panama.MethodsA descriptive and observational study was conducted by analysing demographic and epidemiological data from all malaria cases registered during 2012 in the Comarca Guna Yala, Panama. Malaria intensity indicators were calculated during the study period. Entomological evaluations were performed monthly, from October to December 2012, in the three communities that presented the most intense malaria transmission during the first semester of 2012. Anopheles breeding habitats were also characterized.ResultsDuring the studied period, 6754 blood smears were examined (17.8 % of the total population), and 143 were confirmed as positive for Plasmodium vivax. A significant increase of malaria transmission risk indicators (API: 3.8/1000, SPR: 2.1 %) was observed in Guna Yula, when compared with previous years, and also in comparison with estimates from the whole country. Anopheles albimanus was the most abundant and widespread (877; 72.0 %) vector species found in the three localities, followed by Anopheles punctimacula (231; 19.0 %) and Anopheles aquasalis (110; 9.0 %). Three An. albimanus pools were positive for P. vivax, showing an overall pooled prevalence estimate of 0.014.ConclusionsData analysis confirmed that during 2012 a malaria epidemic occurred in Guna Yala. Panama. This study provides baseline data on the local epidemiology of malaria in this vulnerable region of Panama. This information will be useful for targeting control strategies by the National Malaria Control Programme

Nutritional Systems Biology Modeling: From Molecular Mechanisms to Physiology

The use of computational modeling and simulation has increased in many biological fields, but despite their potential these techniques are only marginally applied in nutritional sciences. Nevertheless, recent applications of modeling have been instrumental in answering important nutritional questions from the cellular up to the physiological levels. Capturing the complexity of today's important nutritional research questions poses a challenge for modeling to become truly integrative in the consideration and interpretation of experimental data at widely differing scales of space and time. In this review, we discuss a selection of available modeling approaches and applications relevant for nutrition. We then put these models into perspective by categorizing them according to their space and time domain. Through this categorization process, we identified a dearth of models that consider processes occurring between the microscopic and macroscopic scale. We propose a “middle-out” strategy to develop the required full-scale, multilevel computational models. Exhaustive and accurate phenotyping, the use of the virtual patient concept, and the development of biomarkers from “-omics” signatures are identified as key elements of a successful systems biology modeling approach in nutrition research—one that integrates physiological mechanisms and data at multiple space and time scales

Propagation of kinetic uncertainties through a canonical topology of the TLR4 signaling network in different regions of biochemical reaction space

<p>Abstract</p> <p>Background</p> <p>Signal transduction networks represent the information processing systems that dictate which dynamical regimes of biochemical activity can be accessible to a cell under certain circumstances. One of the major concerns in molecular systems biology is centered on the elucidation of the robustness properties and information processing capabilities of signal transduction networks. Achieving this goal requires the establishment of causal relations between the design principle of biochemical reaction systems and their emergent dynamical behaviors.</p> <p>Methods</p> <p>In this study, efforts were focused in the construction of a relatively well informed, deterministic, non-linear dynamic model, accounting for reaction mechanisms grounded on standard mass action and Hill saturation kinetics, of the canonical reaction topology underlying Toll-like receptor 4 (TLR4)-mediated signaling events. This signaling mechanism has been shown to be deployed in macrophages during a relatively short time window in response to lypopolysaccharyde (LPS) stimulation, which leads to a rapidly mounted innate immune response. An extensive computational exploration of the biochemical reaction space inhabited by this signal transduction network was performed via local and global perturbation strategies. Importantly, a broad spectrum of biologically plausible dynamical regimes accessible to the network in widely scattered regions of parameter space was reconstructed computationally. Additionally, experimentally reported transcriptional readouts of target pro-inflammatory genes, which are actively modulated by the network in response to LPS stimulation, were also simulated. This was done with the main goal of carrying out an unbiased statistical assessment of the intrinsic robustness properties of this canonical reaction topology.</p> <p>Results</p> <p>Our simulation results provide convincing numerical evidence supporting the idea that a canonical reaction mechanism of the TLR4 signaling network is capable of performing information processing in a robust manner, a functional property that is independent of the signaling task required to be executed. Nevertheless, it was found that the robust performance of the network is not solely determined by its design principle (topology), but this may be heavily dependent on the network's current position in biochemical reaction space. Ultimately, our results enabled us the identification of key rate limiting steps which most effectively control the performance of the system under diverse dynamical regimes.</p> <p>Conclusions</p> <p>Overall, our <it>in silico </it>study suggests that biologically relevant and non-intuitive aspects on the general behavior of a complex biomolecular network can be elucidated only when taking into account a wide spectrum of dynamical regimes attainable by the system. Most importantly, this strategy provides the means for a suitable assessment of the inherent variational constraints imposed by the structure of the system when systematically probing its parameter space.</p