A structural property for reduction of biochemical networks

The effective deficiency of biochemical networks

The deficiency of a (bio)chemical reaction network can be conceptually interpreted as a measure of its ability to support exotic dynamical behavior and/or multistationarity. The classical definition of deficiency relates to the capacity of a network to permit variations of the complex formation rate vector at steady state, irrespective of the network kinetics. However, the deficiency is by definition completely insensitive to the fine details of the directionality of reactions as well as bounds on reaction fluxes. While the classical definition of deficiency can be readily applied in the analysis of unconstrained, weakly reversible networks, it only provides an upper bound in the cases where relevant constraints on reaction fluxes are imposed. Here we propose the concept of effective deficiency, which provides a more accurate assessment of the network's capacity to permit steady state variations at the complex level for constrained networks of any reversibility patterns. The effective deficiency relies on the concept of nonstoichiometric balanced complexes, which we have already shown to be present in real-world biochemical networks operating under flux constraints. Our results demonstrate that the effective deficiency of real-world biochemical networks is smaller than the classical deficiency, indicating the effects of reaction directionality and flux bounds on the variation of the complex formation rate vector at steady state.

Maximizing multi-reaction dependencies provides more accurate and precise predictions of intracellular fluxes than the principle of parsimony

Intracellular fluxes represent a joint outcome of cellular transcription and translation and reflect the availability and usage of nutrients from the environment. While approaches from the constraint-based metabolic framework can accurately predict cellular phenotypes, such as growth and exchange rates with the environment, accurate prediction of intracellular fluxes remains a pressing problem. Parsimonious flux balance analysis (pFBA) has become an approach of choice to predict intracellular fluxes by employing the principle of efficient usage of protein resources. Nevertheless, comparative analyses of intracellular flux predictions from pFBA against fluxes estimated from labeling experiments remain scarce. Here, we posited that steady-state flux distributions derived from the principle of maximizing multi-reaction dependencies are of improved accuracy and precision than those resulting from pFBA. To this end, we designed a constraint-based approach, termed complex-balanced FBA (cbFBA), to predict steady-state flux distributions that support the given specific growth rate and exchange fluxes. We showed that the steady-state flux distributions resulting from cbFBA in comparison to pFBA show better agreement with experimentally measured fluxes from 17 Escherichia coli strains and are more precise, due to the smaller space of alternative solutions. We also showed that the same principle holds in eukaryotes by comparing the predictions of pFBA and cbFBA against experimentally derived steady-state flux distributions from 26 knock-out mutants of Saccharomyces cerevisiae. Furthermore, our results showed that intracellular fluxes predicted by cbFBA provide better support for the principle of minimizing metabolic adjustment between mutants and wild types. Together, our findings point that other principles that consider the dynamics and coordination of steady states may govern the distribution of intracellular fluxes.

A universal dynamical metabolic model representing mixotrophic growth of Chlorella sp. on wastes

An emerging idea is to couple wastewater treatment and biofuel production using microalgae to achieve higher productivities and lower costs. This paper proposes a metabolic modeling of Chlorella sp. growing on fermentation wastes (blend of acetate, butyrate and other acids) in mixotrophic conditions, accounting also for the possible inhibitory substrates. This model extends previous works by modifying the metabolic network to include the consumption of glycerol and glucose by Chlorella sp., with the goal to test the addition of these substrates in order to overcome butyrate inhibition. The metabolic model was built using the DRUM framework and consists of 188 reactions and 173 metabolites. After a calibration phase, the model was successfully challenged with data from 122 experiments collected from scientific literature in autotrophic, heterotrophic and mixotrophic conditions. The optimal feeding strategy estimated with the model reduces the time to consume the volatile fatty acids from 16 days to 2 days. The high prediction capability of this model opens new routes for enhancing design and operation in waste valorization using microalgae.

The hidden simplicity of metabolic networks is revealed by multireaction dependencies

Understanding the complexity of metabolic networks has implications for manipulation of their functions. The complexity of metabolic networks can be characterized by identifying multireaction dependencies that are challenging to determine due to the sheer number of combinations to consider. Here, we propose the concept of concordant complexes that captures multireaction dependencies and can be efficiently determined from the algebraic structure and operational constraints of metabolic networks. The concordant complexes imply the existence of concordance modules based on which the apparent complexity of 12 metabolic networks of organisms from all kingdoms of life can be reduced by at least 78%. A comparative analysis against an ensemble of randomized metabolic networks shows that the metabolic network of Escherichia coli contains fewer concordance modules and is, therefore, more tightly coordinated than expected by chance. Together, our findings demonstrate that metabolic networks are considerably simpler than what can be perceived from their structure alone.

Multiscale kinetic analysis of proteins

The stochasticity of molecular motion results in the existence of multiple kinetically relevant pathways in many biomolecular mechanisms. Because it is highly demanding to characterize them for complex systems, mechanisms are often described with a single-pathway perspective. However, kinetic network analysis and sub-ensemble experimental insight are increasingly demonstrating not only the existence of competing pathways but also the importance of kinetic selection in biology. This review focuses on advances in multiscale kinetic analysis of proteins, which connects molecular level information from simulations to macroscopic data to characterize mechanistic reaction networks and the reactive flux through them. We describe a range of methods used and highlight several examples where kinetic modeling has revealed functional importance of pathway heterogeneity.