Genome-scale Biochemical Networks Research Articles

CHRR: coordinate hit-and-run with rounding for uniform sampling of constraint-based models.

SummaryIn constraint-based metabolic modelling, physical and biochemical constraints define a polyhedral convex set of feasible flux vectors. Uniform sampling of this set provides an unbiased characterization of the metabolic capabilities of a biochemical network. However, reliable uniform sampling of genome-scale biochemical networks is challenging due to their high dimensionality and inherent anisotropy. Here, we present an implementation of a new sampling algorithm, coordinate hit-and-run with rounding (CHRR). This algorithm is based on the provably efficient hit-and-run random walk and crucially uses a preprocessing step to round the anisotropic flux set. CHRR provably converges to a uniform stationary sampling distribution. We apply it to metabolic networks of increasing dimensionality. We show that it converges several times faster than a popular artificial centering hit-and-run algorithm, enabling reliable and tractable sampling of genome-scale biochemical networks.Availability and Implementation https://github.com/opencobra/cobratoolbox.Supplementary information Supplementary data are available at Bioinformatics online.

Conditions for duality between fluxes and concentrations in biochemical networks

Mathematical and computational modelling of biochemical networks is often done in terms of either the concentrations of molecular species or the fluxes of biochemical reactions. When is mathematical modelling from either perspective equivalent to the other? Mathematical duality translates concepts, theorems or mathematical structures into other concepts, theorems or structures, in a one-to-one manner. We present a novel stoichiometric condition that is necessary and sufficient for duality between unidirectional fluxes and concentrations. Our numerical experiments, with computational models derived from a range of genome-scale biochemical networks, suggest that this flux-concentration duality is a pervasive property of biochemical networks. We also provide a combinatorial characterisation that is sufficient to ensure flux-concentration duality.The condition prescribes that, for every two disjoint sets of molecular species, there is at least one reaction complex that involves species from only one of the two sets. When unidirectional fluxes and molecular species concentrations are dual vectors, this implies that the behaviour of the corresponding biochemical network can be described entirely in terms of either concentrations or unidirectional fluxes.

Systematic evaluation of methods for integration of transcriptomic data into constraint-based models of metabolism.

Constraint-based models of metabolism are a widely used framework for predicting flux distributions in genome-scale biochemical networks. The number of published methods for integration of transcriptomic data into constraint-based models has been rapidly increasing. So far the predictive capability of these methods has not been critically evaluated and compared. This work presents a survey of recently published methods that use transcript levels to try to improve metabolic flux predictions either by generating flux distributions or by creating context-specific models. A subset of these methods is then systematically evaluated using published data from three different case studies in E. coli and S. cerevisiae. The flux predictions made by different methods using transcriptomic data are compared against experimentally determined extracellular and intracellular fluxes (from 13C-labeling data). The sensitivity of the results to method-specific parameters is also evaluated, as well as their robustness to noise in the data. The results show that none of the methods outperforms the others for all cases. Also, it is observed that for many conditions, the predictions obtained by simple flux balance analysis using growth maximization and parsimony criteria are as good or better than those obtained using methods that incorporate transcriptomic data. We further discuss the differences in the mathematical formulation of the methods, and their relation to the results we have obtained, as well as the connection to the underlying biological principles of metabolic regulation.

Design and strain selection criteria for bacterial communication networks

In this paper, we discuss data and methodological challenges for building bacterial communication networks using two examples: E. coli as a flagellate bacterium and G. sulfurreducens as a biofilm forming bacterium. We first highlight the link between the bacterial network communication design with respect to metabolic information processing design. The potentialities of designing routing network schemes described previously in the literature and based on bacteria motility and genetic message exchanges will depend on the genes coding for the intracellular and intercellular signalling pathways. In bacteria, the “mobilome” is related to horizontal gene transfer. Bacteria trade-off the acquisition of new genes which could improve their survival (and often their communication bandwidth), keeping their genome small enough to ensure quick DNA replication and fast increase of the biomass to speed up cell division. First, by using a multi-objective optimisation procedure, we search for the optimal trade-off between energy production, which is a requirement for the motility, and biomass growth, which is related to the overall survival and fitness of the bacterium. We use flux balance analysis of the genome-scale biochemical network of E. coli k-13 MG1655. Then, as a second case study, we analyse the electric properties and biomass trade-off of the bacterium G. sulfurreducens, which constructs an electric biofilm where electrons move across the nanowires. Here we discuss the potentialities of optimisation methodologies to design and select bacterial strains with desiderata properties. The optimisation methodologies establish also a relation between metabolism, network communication and computing. Moreover, we point to genetic design and synthetic biology as key areas to develop bacterial nano communication networks.

A variational principle for computing nonequilibrium fluxes and potentials in genome-scale biochemical networks

We derive a convex optimization problem on a steady-state nonequilibrium network of biochemical reactions, with the property that energy conservation and the second law of thermodynamics both hold at the problem solution. This suggests a new variational principle for biochemical networks that can be implemented in a computationally tractable manner. We derive the Lagrange dual of the optimization problem and use strong duality to demonstrate that a biochemical analogue of Tellegen׳s theorem holds at optimality. Each optimal flux is dependent on a free parameter that we relate to an elementary kinetic parameter when mass action kinetics is assumed.

Quantitative assignment of reaction directionality in constraint-based models of metabolism: Application to Escherichia coli

Constraint-based modeling is an approach for quantitative prediction of net reaction flux in genome-scale biochemical networks. In vivo, the second law of thermodynamics requires that net macroscopic flux be forward, when the transformed reaction Gibbs energy is negative. We calculate the latter by using (i) group contribution estimates of metabolite species Gibbs energy, combined with (ii) experimentally measured equilibrium constants. In an application to a genome-scale stoichiometric model of Escherichia coli metabolism, iAF1260, we demonstrate that quantitative prediction of reaction directionality is increased in scope and accuracy by integration of both data sources, transformed appropriately to in vivo pH, temperature and ionic strength. Comparison of quantitative versus qualitative assignment of reaction directionality in iAF1260, assuming an accommodating reactant concentration range of 0.02–20 mM, revealed that quantitative assignment leads to a low false positive, but high false negative, prediction of effectively irreversible reactions. The latter is partly due to the uncertainty associated with group contribution estimates. We also uncovered evidence that the high intracellular concentration of glutamate in E. coli may be essential to direct otherwise thermodynamically unfavorable essential reactions, such as the leucine transaminase reaction, in an anabolic direction.

Properties of Metabolic Networks: Structure versus Function

Biological data from high-throughput technologies describing the network components (genes, proteins, metabolites) and their associated interactions have driven the reconstruction and study of structural (topological) properties of large-scale biological networks. In this article, we address the relation of the functional and structural properties by using extensively experimentally validated genome-scale metabolic network models to compute observable functional states of a microorganism and compare the “structure versus function” attributes of metabolic networks. It is observed that, functionally speaking, the essentiality of reactions in a node is not correlated with node connectivity as structural analyses of other biological networks have suggested. These findings are illustrated with the analysis of the genome-scale biochemical networks of three species with distinct modes of metabolism. These results also suggest fundamental differences among different biological networks arising out of their representation and functional constraints.