Applications and trends in systems biology in biochemistry

A Fermentation State Marker Rule Design Task in Metabolic Engineering

There are several ways in which mathematical modeling is used in fermentation control, but mechanistic mathematical genome-scale models of metabolism within the cell have not been applied or implemented so far. As part of the metabolic engineering task setting, we propose that metabolite fluxes and/or biomass growth rate be used to search for a fermentation steady state marker rule. During fermentation, the bioreactor control system can automatically detect the desired steady state using a logical marker rule. The marker rule identification can be also integrated with the production growth coupling approach, as presented in this study. A design of strain with marker rule is demonstrated on genome scale metabolic model iML1515 of Escherichia coli MG1655 proposing two gene deletions enabling a measurable marker rule for succinate production using glucose as a substrate. The marker rule example at glucose consumption 10.0 is: IF (specific growth rate μ is above 0.060 h-1, AND CO2 production under 1.0, AND ethanol production above 5.5), THEN succinate production is within the range 8.2-10, where all metabolic fluxes units are mmol ∗ gDW-1 ∗ h-1. An objective function for application in metabolic engineering, including productivity features and rule detecting sensor set characterizing parameters, is proposed. Two-phase approach to implementing marker rules in the cultivation control system is presented to avoid the need for a modeler during production.

A comprehensive mechanistic model of adipocyte signaling with layers of confidence

Adipocyte signaling, normally and in type 2 diabetes, is far from fully understood. We have earlier developed detailed dynamic mathematical models for several well-studied, partially overlapping, signaling pathways in adipocytes. Still, these models only cover a fraction of the total cellular response. For a broader coverage of the response, large-scale phosphoproteomic data and systems level knowledge on protein interactions are key. However, methods to combine detailed dynamic models with large-scale data, using information about the confidence of included interactions, are lacking. We have developed a method to first establish a core model by connecting existing models of adipocyte cellular signaling for: (1) lipolysis and fatty acid release, (2) glucose uptake, and (3) the release of adiponectin. Next, we use publicly available phosphoproteome data for the insulin response in adipocytes together with prior knowledge on protein interactions, to identify phosphosites downstream of the core model. In a parallel pairwise approach with low computation time, we test whether identified phosphosites can be added to the model. We iteratively collect accepted additions into layers and continue the search for phosphosites downstream of these added layers. For the first 30 layers with the highest confidence (311 added phosphosites), the model predicts independent data well (70-90% correct), and the predictive capability gradually decreases when we add layers of decreasing confidence. In total, 57 layers (3059 phosphosites) can be added to the model with predictive ability kept. Finally, our large-scale, layered model enables dynamic simulations of systems-wide alterations in adipocytes in type 2 diabetes.

Reproducibility and FAIR principles: the case of a segment polarity network model

The issue of reproducibility of computational models and the related FAIR principles (findable, accessible, interoperable, and reusable) are examined in a specific test case. I analyze a computational model of the segment polarity network in Drosophila embryos published in 2000. Despite the high number of citations to this publication, 23 years later the model is barely accessible, and consequently not interoperable. Following the text of the original publication allowed successfully encoding the model for the open source software COPASI. Subsequently saving the model in the SBML format allowed it to be reused in other open source software packages. Submission of this SBML encoding of the model to the BioModels database enables its findability and accessibility. This demonstrates how the FAIR principles can be successfully enabled by using open source software, widely adopted standards, and public repositories, facilitating reproducibility and reuse of computational cell biology models that will outlive the specific software used.

Computational modeling and quantitative physiology reveal central parameters for brassinosteroid-regulated early cell physiological processes linked to elongation growth of the Arabidopsis root

Brassinosteroids (BR) are key hormonal regulators of plant development. However, whereas the individual components of BR perception and signaling are well characterized experimentally, the question of how they can act and whether they are sufficient to carry out the critical function of cellular elongation remains open. Here, we combined computational modeling with quantitative cell physiology to understand the dynamics of the plasma membrane (PM)-localized BR response pathway during the initiation of cellular responses in the epidermis of the Arabidopsis root tip that are be linked to cell elongation. The model, consisting of ordinary differential equations, comprises the BR-induced hyperpolarization of the PM, the acidification of the apoplast and subsequent cell wall swelling. We demonstrate that the competence of the root epidermal cells for the BR response predominantly depends on the amount and activity of H+-ATPases in the PM. The model further predicts that an influx of cations is required to compensate for the shift of positive charges caused by the apoplastic acidification. A potassium channel was subsequently identified and experimentally characterized, fulfilling this function. Thus, we established the landscape of components and parameters for physiological processes potentially linked to cell elongation, a central process in plant development.

Dynamic publication media with the COPASI R Connector (CoRC)

In this article we show how dynamic publication media and the COPASI R Connector (CoRC) can be combined in a natural and synergistic way to communicate (biochemical) models. Dynamic publication media are becoming a popular tool for authors to effectively compose and publish their work. They are built from templates and the final documents are created dynamically. In addition, they can also be interactive. Working with dynamic publication media is made easy with the programming environment R via its integration with tools such as R Markdown, Jupyter and Shiny. Additionally, the COmplex PAthway SImulator COPASI (http://www.copasi.org), a widely used biochemical modelling toolkit, is available in R through the use of the COPASI R Connector (CoRC, https://jpahle.github.io/CoRC). Models are a common tool in the mathematical biosciences, in particular kinetic models of biochemical networks in (computational) systems biology. We focus on three application areas of dynamic publication media and CoRC: Documentation (reproducible workflows), Teaching (creating self-paced lessons) and Science Communication (immersive and engaging presentation). To illustrate these, we created six dynamic document examples in the form of R Markdown and Jupyter notebooks, hosted on the platforms GitHub, shinyapps.io, Google Colaboratory. Having code and output in one place, creating documents in template-form and the option of interactivity make the combination of dynamic documents and CoRC a versatile tool. All our example documents are freely available at https://jpahle.github.io/DynamiCoRC under the Creative Commons BY 4.0 licence.

Impact of explicit area scaling on kinetic models involving multiple compartments

BACKGROUND

Computational modelling of cell biological processes is a frequently used technique to analyse the underlying mechanisms and to generally understand the behaviour of these processes in the context of a pathway, network or even the whole cell. The most common technique in this context is the usage of ordinary differential equations that describe the kinetics of the relevant processes in mechanistic detail. Here, it is usually assumed that the content of the cell is well-stirred and thus homogeneous - which is of course an over-simplification, but often worked in the past. However, many processes happen at membranes and thus not in 3D, but in 2D. The scaling of the rates of these processes poses a special problem, if volumes of compartments are changed. They will typically scale with an area, but not with the volume of the involved compartment. However, commonly, this is neglected when setting up models and/or volume scaling also sometimes automatically happens when using modelling software in the field.

RESULTS

Here, we investigate generic as well as specific, realistic cases to find out, how strong the impact of the wrong scaling is for the outcome of simulations. We show that the importance of correct area scaling depends on the architecture of the reaction site and its changes upon volume alterations and it is hard to foresee, if it has a significant impact or not just by looking at the original model set-up. Moreover, scaled rates might exhibit more or less control over the behaviour of the system and therefore, accordingly, incorrect scaling will have more or less influence.

CONCLUSIONS

Working with multi-compartment reactions requires a careful consideration of the correct scaling of the rates when changing the volumes of the involved compartments. The error following incorrect scaling - often done by scaling with the volume of the respective compartments can lead to significant aberrations of model behaviour.

The Construction of Biological 'Inter-Identity' as the Outcome of a Complex Process of Protocell Development in Prebiotic Evolution

The concept of identity is used both (i) to distinguish a system as a particular material entity that is conserved as such in a given environment (token-identity: i.e., identity as permanence or endurance over time), and (ii) to relate a system with other members of a set (type-identity: i.e., identity as an equivalence relationship). Biological systems are characterized, in a minimal and universal sense, by a highly complex and dynamic, far-from-equilibrium organization of very diverse molecular components and transformation processes (i.e., 'genetically instructed cellular metabolisms') that maintain themselves in constant interaction with their corresponding environments, including other systems of similar nature. More precisely, all living entities depend on a deeply convoluted organization of molecules and processes (a naturalized von Neumann constructor architecture) that subsumes, in the form of current individuals (autonomous cells), a history of ecological and evolutionary interactions (across cell populations). So one can defend, on those grounds, that living beings have an identity of their own from both approximations: (i) and (ii). These transversal and trans-generational dimensions of biological phenomena, which unfold together with the actual process of biogenesis, must be carefully considered in order to understand the intricacies and metabolic robustness of the first living cells, their underlying uniformity (i.e., their common biochemical core) and the eradication of previous -or alternative- forms of complex natural phenomena. Therefore, a comprehensive approach to the origins of life requires conjugating the actual properties of the developing complex individuals (fusing and dividing protocells, at various stages) with other, population-level features, linked to their collective-evolutionary behavior, under much wider and longer-term parameters. On these lines, we will argue that life, in its most basic sense, here on Earth or anywhere else, demands crossing a high complexity threshold and that the concept of 'inter-identity' can help us realize the different aspects involved in the process. The article concludes by pointing out some of the challenges ahead if we are to integrate the corresponding explanatory frameworks, physiological and evolutionary, in the hope that a more general theory of biology is on its way.

Computational systems biology of cellular processes in Arabidopsis thaliana: an overview

Systems biology strives for gaining an understanding of biological phenomena by studying the interactions of different parts of a system and integrating the knowledge obtained into the current view of the underlying processes. This is achieved by a tight combination of quantitative experimentation and computational modeling. While there is already a large quantity of systems biology studies describing human, animal and especially microbial cell biological systems, plant biology has been lagging behind for many years. However, in the case of the model plant Arabidopsis thaliana, the steadily increasing amount of information on the levels of its genome, proteome and on a variety of its metabolic and signalling pathways is progressively enabling more researchers to construct models for cellular processes for the plant, which in turn encourages more experimental data to be generated, showing also for plant sciences how fruitful systems biology research can be. In this review, we provide an overview over some of these recent studies which use different systems biological approaches to get a better understanding of the cell biology of A. thaliana. The approaches used in these are genome-scale metabolic modeling, as well as kinetic modeling of metabolic and signalling pathways. Furthermore, we selected several cases to exemplify how the modeling approaches have led to significant advances or new perspectives in the field.

Approximate analytical solution of fractional order biochemical reaction model and its stability analysis

Approximate analytical solution of the system of coupled nonlinear Ordinary Differential Equations (ODEs) of a biochemical reaction model is much relevant due to its practical significance to biochemists. In this paper, an effective and powerful mathematical technique, viz. fractional homotopy analysis transform method (FHATM), is employed to get the numerical solutions of biochemical reaction model with time fractional derivatives. The adopted scheme is the beautiful copulation of homotopy analysis technique and Laplace transform algorithm. This paper shows that the adopted scheme is quite easy as well as computationally attractive in the context of a solution procedure. The Caputo-type fractional derivatives are considered in the present paper. Approximate results of the probability density functions of the time fractional biochemical reaction model are computed for miscellaneous fractional Brownian motions as well as for classical motion and are presented graphically. The time fractional biochemical reaction model with respect to stability analysis for various values of fractional order [Formula: see text] is also analyzed. In the context of stability discussion, we have used the fractional Routh–Hurwitz stability criterion to establish the local stability of the biochemical reaction model of fractional order.

Distribution of control in the sulfur assimilation in Arabidopsis thaliana depends on environmental conditions

Sulfur assimilation is central to the survival of plants and has been studied under different environmental conditions. Multiple studies have been published trying to determine rate-limiting or controlling steps in this pathway. However, the picture remains inconclusive with at least two different enzymes proposed to represent such rate-limiting steps. Here, we used computational modeling to gain an integrative understanding of the distribution of control in the sulfur assimilation pathway of Arabidopsis thaliana. For this purpose, we set up a new ordinary differential equation (ODE)-based, kinetic model of sulfur assimilation encompassing all biochemical reactions directly involved in this process. We fitted the model to published experimental data and produced a model ensemble to deal with parameter uncertainties. The ensemble was validated against additional published experimental data. We used the model ensemble to subsequently analyse the control pattern and robustly identified a set of processes that share the control in this pathway under standard conditions. Interestingly, the pattern of control is dynamic and not static, that is it changes with changing environmental conditions. Therefore, while adenosine-5'-phosphosulfate reductase (APR) and sulfite reductase (SiR) share control under standard laboratory conditions, APR takes over an even more dominant role under sulfur starvation conditions.

Reproducible Research Using Biomodels

Like other types of computational research, modeling and simulation of biological processes (biomodels) is still largely communicated without sufficient detail to allow independent reproduction of results. But reproducibility in this area of research could easily be achieved by making use of existing resources, such as supplying models in standard formats and depositing code, models, and results in public repositories.

MAMMOTh: A new database for curated mathematical models of biomolecular systems

MOTIVATION

Living systems have a complex hierarchical organization that can be viewed as a set of dynamically interacting subsystems. Thus, to simulate the internal nature and dynamics of the entire biological system, we should use the iterative way for a model reconstruction, which is a consistent composition and combination of its elementary subsystems. In accordance with this bottom-up approach, we have developed the MAthematical Models of bioMOlecular sysTems (MAMMOTh) tool that consists of the database containing manually curated MAMMOTh fitted to the experimental data and a software tool that provides their further integration.

RESULTS

The MAMMOTh database entries are organized as building blocks in a way that the model parts can be used in different combinations to describe systems with higher organizational level (metabolic pathways and/or transcription regulatory networks). The tool supports export of a single model or their combinations in SBML or Mathematica standards. The database currently contains 110 mathematical sub-models for Escherichia coli elementary subsystems (enzymatic reactions and gene expression regulatory processes) that can be combined in at least 5100 complex/sophisticated models concerning more complex biological processes as de novo nucleotide biosynthesis, aerobic/anaerobic respiration and nitrate/nitrite utilization in E. coli. All models are functionally interconnected and sufficiently complement public model resources.

AVAILABILITY

http://mammoth.biomodelsgroup.ru.

COPASI and its applications in biotechnology

COPASI is software used for the creation, modification, simulation and computational analysis of kinetic models in various fields. It is open-source, available for all major platforms and provides a user-friendly graphical user interface, but is also controllable via the command line and scripting languages. These are likely reasons for its wide acceptance. We begin this review with a short introduction describing the general approaches and techniques used in computational modeling in the biosciences. Next we introduce the COPASI package, and its capabilities, before looking at typical applications of COPASI in biotechnology.

Exploiting intrinsic fluctuations to identify model parameters

Parameterisation of kinetic models plays a central role in computational systems biology. Besides the lack of experimental data of high enough quality, some of the biggest challenges here are identification issues. Model parameters can be structurally non‐identifiable because of functional relationships. Noise in measured data is usually considered to be a nuisance for parameter estimation. However, it turns out that intrinsic fluctuations in particle numbers can make parameters identifiable that were previously non‐identifiable. The authors present a method to identify model parameters that are structurally non‐identifiable in a deterministic framework. The method takes time course recordings of biochemical systems in steady state or transient state as input. Often a functional relationship between parameters presents itself by a one‐dimensional manifold in parameter space containing parameter sets of optimal goodness. Although the system's behaviour cannot be distinguished on this manifold in a deterministic framework it might be distinguishable in a stochastic modelling framework. Their method exploits this by using an objective function that includes a measure for fluctuations in particle numbers. They show on three example models, immigration‐death, gene expression and Epo‐EpoReceptor interaction, that this resolves the non‐identifiability even in the case of measurement noise with known amplitude. The method is applied to partially observed recordings of biochemical systems with measurement noise. It is simple to implement and it is usually very fast to compute. This optimisation can be realised in a classical or Bayesian fashion.

The (Mathematical) Modeling Process in Biosciences

In this communication, we introduce a general framework and discussion on the role of models and the modeling process in the field of biosciences. The objective is to sum up the common procedures during the formalization and analysis of a biological problem from the perspective of Systems Biology, which approaches the study of biological systems as a whole. We begin by presenting the definitions of (biological) system and model. Particular attention is given to the meaning of mathematical model within the context of biology. Then, we present the process of modeling and analysis of biological systems. Three stages are described in detail: conceptualization of the biological system into a model, mathematical formalization of the previous conceptual model and optimization and system management derived from the analysis of the mathematical model. All along this work the main features and shortcomings of the process are analyzed and a set of rules that could help in the task of modeling any biological system are presented. Special regard is given to the formative requirements and the interdisciplinary nature of this approach. We conclude with some general considerations on the challenges that modeling is posing to current biology.

Nonlinear mixed-effects modelling for single cell estimation: when, why, and how to use it

Background

Studies of cell-to-cell variation have in recent years grown in interest, due to improved bioanalytical techniques which facilitates determination of small changes with high uncertainty. Like much high-quality data, single-cell data is best analysed using a systems biology approach. The most common systems biology approach to single-cell data is the standard two-stage (STS) approach. In STS, data from each cell is analysed in a separate sub-problem, meaning that only data from the same cell is used to calculate the parameter values within that cell. Because only parts of the data are considered, problems with parameter unidentifiability are exaggerated in STS. In contrast, a related approach to data analysis has been developed for the studies of patient-to-patient variations. This approach, called nonlinear mixed-effects modelling (NLME), makes use of all data, when estimating the patient-specific parameters. NLME would therefore be advantageous compared to STS also for the study of cell-to-cell variation. However, no such systematic evaluation of the two approaches exists.

Results

Herein, such a systematic comparison between STS and NLME has been performed. Different examples, both linear and nonlinear, and both simulated and real experimental data, have been examined. With informative data, there is no significant difference in the results for either parameter or noise estimation. However, when data becomes uninformative, NLME is significantly superior to STS. These results hold independently of whether the loss of information is due to a low signal-to-noise ratio, too few data points, or a bad input signal. The improvement is shown to come from both the consideration of a joint likelihood (JLH) function, describing all parameters and data, and from an a priori postulated form of the population parameters. Finally, we provide a small tutorial that shows how to use NLME for single-cell analysis, using the free and user-friendly software Monolix.

Conclusions

When considering uninformative single-cell data, NLME yields more accurate parameter and noise estimates, compared to more traditional approaches, such as STS and JLH.

Electronic supplementary material

The online version of this article (doi:10.1186/s12918-015-0203-x) contains supplementary material, which is available to authorized users.

Flux balance analysis predicts essential genes in clear cell renal cell carcinoma metabolism

Flux balance analysis is the only modelling approach that is capable of producing genome-wide predictions of gene essentiality that may aid to unveil metabolic liabilities in cancer. Nevertheless, a systemic validation of gene essentiality predictions by flux balance analysis is currently missing. Here, we critically evaluated the accuracy of flux balance analysis in two cancer types, clear cell renal cell carcinoma (ccRCC) and prostate adenocarcinoma, by comparison with large-scale experiments of gene essentiality in vitro. We found that in ccRCC, but not in prostate adenocarcinoma, flux balance analysis could predict essential metabolic genes beyond random expectation. Five of the identified metabolic genes, AGPAT6, GALT, GCLC, GSS, and RRM2B, were predicted to be dispensable in normal cell metabolism. Hence, targeting these genes may selectively prevent ccRCC growth. Based on our analysis, we discuss the benefits and limitations of flux balance analysis for gene essentiality predictions in cancer metabolism, and its use for exposing metabolic liabilities in ccRCC, whose emergent metabolic network enforces outstanding anabolic requirements for cellular proliferation.

Integrative analysis of metabolic models - from structure to dynamics

The characterization of biological systems with respect to their behavior and functionality based on versatile biochemical interactions is a major challenge. To understand these complex mechanisms at systems level modeling approaches are investigated. Different modeling formalisms allow metabolic models to be analyzed depending on the question to be solved, the biochemical knowledge and the availability of experimental data. Here, we describe a method for an integrative analysis of the structure and dynamics represented by qualitative and quantitative metabolic models. Using various formalisms, the metabolic model is analyzed from different perspectives. Determined structural and dynamic properties are visualized in the context of the metabolic model. Interaction techniques allow the exploration and visual analysis thereby leading to a broader understanding of the behavior and functionality of the underlying biological system. The System Biology Metabolic Model Framework (SBM (2) - Framework) implements the developed method and, as an example, is applied for the integrative analysis of the crop plant potato.

Computationally Modeling Lipid Metabolism and Aging: A Mini-review

One of the greatest challenges in biology is to improve the understanding of the mechanisms which underpin aging and how these affect health. The need to better understand aging is amplified by demographic changes, which have caused a gradual increase in the global population of older people. Aging western populations have resulted in a rise in the prevalence of age-related pathologies. Of these diseases, cardiovascular disease is the most common underlying condition in older people. The dysregulation of lipid metabolism due to aging impinges significantly on cardiovascular health. However, the multifaceted nature of lipid metabolism and the complexities of its interaction with aging make it challenging to understand by conventional means. To address this challenge computational modeling, a key component of the systems biology paradigm is being used to study the dynamics of lipid metabolism. This mini-review briefly outlines the key regulators of lipid metabolism, their dysregulation, and how computational modeling is being used to gain an increased insight into this system.

Combining test statistics and models in bootstrapped model rejection: it is a balancing act

BACKGROUND

Model rejections lie at the heart of systems biology, since they provide conclusive statements: that the corresponding mechanistic assumptions do not serve as valid explanations for the experimental data. Rejections are usually done using e.g. the chi-square test (χ2) or the Durbin-Watson test (DW). Analytical formulas for the corresponding distributions rely on assumptions that typically are not fulfilled. This problem is partly alleviated by the usage of bootstrapping, a computationally heavy approach to calculate an empirical distribution. Bootstrapping also allows for a natural extension to estimation of joint distributions, but this feature has so far been little exploited.

RESULTS

We herein show that simplistic combinations of bootstrapped tests, like the max or min of the individual p-values, give inconsistent, i.e. overly conservative or liberal, results. A new two-dimensional (2D) approach based on parametric bootstrapping, on the other hand, is found both consistent and with a higher power than the individual tests, when tested on static and dynamic examples where the truth is known. In the same examples, the most superior test is a 2D χ2vsχ2, where the second χ2-value comes from an additional help model, and its ability to describe bootstraps from the tested model. This superiority is lost if the help model is too simple, or too flexible. If a useful help model is found, the most powerful approach is the bootstrapped log-likelihood ratio (LHR). We show that this is because the LHR is one-dimensional, because the second dimension comes at a cost, and because LHR has retained most of the crucial information in the 2D distribution. These approaches statistically resolve a previously published rejection example for the first time.

CONCLUSIONS

We have shown how to, and how not to, combine tests in a bootstrap setting, when the combination is advantageous, and when it is advantageous to include a second model. These results also provide a deeper insight into the original motivation for formulating the LHR, for the more general setting of nonlinear and non-nested models. These insights are valuable in cases when accuracy and power, rather than computational speed, are prioritized.

What can we learn from global sensitivity analysis of biochemical systems?

Most biological models of intermediate size, and probably all large models, need to cope with the fact that many of their parameter values are unknown. In addition, it may not be possible to identify these values unambiguously on the basis of experimental data. This poses the question how reliable predictions made using such models are. Sensitivity analysis is commonly used to measure the impact of each model parameter on its variables. However, the results of such analyses can be dependent on an exact set of parameter values due to nonlinearity. To mitigate this problem, global sensitivity analysis techniques are used to calculate parameter sensitivities in a wider parameter space. We applied global sensitivity analysis to a selection of five signalling and metabolic models, several of which incorporate experimentally well-determined parameters. Assuming these models represent physiological reality, we explored how the results could change under increasing amounts of parameter uncertainty. Our results show that parameter sensitivities calculated with the physiological parameter values are not necessarily the most frequently observed under random sampling, even in a small interval around the physiological values. Often multimodal distributions were observed. Unsurprisingly, the range of possible sensitivity coefficient values increased with the level of parameter uncertainty, though the amount of parameter uncertainty at which the pattern of control was able to change differed among the models analysed. We suggest that this level of uncertainty can be used as a global measure of model robustness. Finally a comparison of different global sensitivity analysis techniques shows that, if high-throughput computing resources are available, then random sampling may actually be the most suitable technique.

Cardiovascular risk assessment in prediabetes and undiagnosed diabetes mellitus study: international collaboration research overview

The study aims to develop a screening protocol for the risk of future cardiovascular disease and diabetes mellitus in people with prediabetes and undiagnosed diabetes; and to establish a framework for early identification and intervention of prediabetes including strategies for holistic management and monitoring of progression. The first phase is to identify prediabetes and undiagnosed diabetes in volunteers who are ≥18-years-old for 5 years. Point-of-care testing and questionnaire will be used to screen for prediabetes and cardiovascular disease. We anticipate screening more than 2000 individuals of both genders by the end of first phase. The second and third phases which shall run for 5-10 years will be longitudinal study involving participants identified in the first phase as having prediabetes without dyslipidaemia, or clinically established cardiovascular disease. The second phase shall focus on preventive management of risk of progress to diabetes with explicit diagnosis of cardiovascular disease. Oxidative stress measurements will be performed cum evaluation of the use of antioxidants, exercise, and nutrition. The third phase will include probing the development of diabetes and cardiovascular disease. Binomial logistic regression would be performed to generate and propose a model chart for the assessment of cardiovascular disease risk in prediabetes.

Can Computer Meridian Diagnostic be useful in diagnosis and management of diseases?

The need for the development of appropriate guidelines for effective or safe use of antioxidants and herbs has always been a concern, especially for the alternative medicine practices. Computer Meridian Diagnostic (CMD) is one of emerging computer-based diagnostic technologies available to alternative medicine practitioners. However, case report of the agents monitored with CMD is uncommon; and concerted effort to bring this into conventional medical practice is yet to be. This hypothesis builds on an anecdotal observation of anti-stress effect monitored with CMD, with a view to highlight a potential tool that requires expatiation, as well as proof of concept and validation studies for possible integration in conventional and traditional medicine practices for therapeutic monitoring.

From top-down to bottom-up: computational modeling approaches for cellular redoxin networks

SIGNIFICANCE

Thioredoxin, glutaredoxin, and peroxiredoxin systems play critical roles in a large number of redox-sensitive cellular processes. These systems are linked to each other by coupled redox cycles and common reaction intermediates into a larger network. Given the scale and connectivity of this network, computational approaches are required to analyze its dynamics and organization.

RECENT ADVANCES

Theoretical advances, as well as new redox proteomic methods, have led to the development of both top-down and bottom-up systems biology approaches to analyze the these systems and the network as a whole. Top-down approaches have been based on modifications to the Nernst equation or on graph theoretical approaches, while bottom-up approaches have been based on kinetic or stoichiometric modeling techniques.

CRITICAL ISSUES

This review will consider the rationale behind these approaches and focus on their advantages and limitations. Further, the review will discuss modeling standards to ensure model accuracy and availability.

FUTURE DIRECTIONS

Top-down and bottom-up approaches have distinct strengths and limitations in describing cellular redoxin networks. The availability of methods to overcome these limitations, together with the adoption of common modeling standards, is expected to increase the pace of model-led discovery within the redox biology field.

Yeast 5 - an expanded reconstruction of the Saccharomyces cerevisiae metabolic network

Background

Efforts to improve the computational reconstruction of the Saccharomyces cerevisiae biochemical reaction network and to refine the stoichiometrically constrained metabolic models that can be derived from such a reconstruction have continued since the first stoichiometrically constrained yeast genome scale metabolic model was published in 2003. Continuing this ongoing process, we have constructed an update to the Yeast Consensus Reconstruction, Yeast 5. The Yeast Consensus Reconstruction is a product of efforts to forge a community-based reconstruction emphasizing standards compliance and biochemical accuracy via evidence-based selection of reactions. It draws upon models published by a variety of independent research groups as well as information obtained from biochemical databases and primary literature.

Results

Yeast 5 refines the biochemical reactions included in the reconstruction, particularly reactions involved in sphingolipid metabolism; updates gene-reaction annotations; and emphasizes the distinction between reconstruction and stoichiometrically constrained model. Although it was not a primary goal, this update also improves the accuracy of model prediction of viability and auxotrophy phenotypes and increases the number of epistatic interactions. This update maintains an emphasis on standards compliance, unambiguous metabolite naming, and computer-readable annotations available through a structured document format. Additionally, we have developed MATLAB scripts to evaluate the model’s predictive accuracy and to demonstrate basic model applications such as simulating aerobic and anaerobic growth. These scripts, which provide an independent tool for evaluating the performance of various stoichiometrically constrained yeast metabolic models using flux balance analysis, are included as Additional files 1, 2 and 3.

Conclusions

Yeast 5 expands and refines the computational reconstruction of yeast metabolism and improves the predictive accuracy of a stoichiometrically constrained yeast metabolic model. It differs from previous reconstructions and models by emphasizing the distinction between the yeast metabolic reconstruction and the stoichiometrically constrained model, and makes both available as Additional file 4 and Additional file 5 and at http://yeast.sf.net/ as separate systems biology markup language (SBML) files. Through this separation, we intend to make the modeling process more accessible, explicit, transparent, and reproducible.

Simplification of biochemical models: a general approach based on the analysis of the impact of individual species and reactions on the systems dynamics

Background

Given the complex mechanisms underlying biochemical processes systems biology researchers tend to build ever increasing computational models. However, dealing with complex systems entails a variety of problems, e.g. difficult intuitive understanding, variety of time scales or non-identifiable parameters. Therefore, methods are needed that, at least semi-automatically, help to elucidate how the complexity of a model can be reduced such that important behavior is maintained and the predictive capacity of the model is increased. The results should be easily accessible and interpretable. In the best case such methods may also provide insight into fundamental biochemical mechanisms.

Results

We have developed a strategy based on the Computational Singular Perturbation (CSP) method which can be used to perform a "biochemically-driven" model reduction of even large and complex kinetic ODE systems. We provide an implementation of the original CSP algorithm in COPASI (a COmplex PAthway SImulator) and applied the strategy to two example models of different degree of complexity - a simple one-enzyme system and a full-scale model of yeast glycolysis.

Conclusion

The results show the usefulness of the method for model simplification purposes as well as for analyzing fundamental biochemical mechanisms. COPASI is freely available at http://www.copasi.org.