NetFlow: A tool for isolating carbon flows in genome-scale metabolic networks.

Genome-scale stoichiometric models (GSMs) have been widely utilized to predict and understand cellular metabolism. GSMs and the flux predictions resulting from them have proven indispensable to fields ranging from metabolic engineering to human disease. Nonetheless, it is challenging to parse these flux predictions due to the inherent size and complexity of the GSMs. Several previous approaches have reduced this complexity by identifying key pathways contained within the genome-scale flux predictions. However, a reduction method that overlays carbon atom transitions on stoichiometry and flux predictions is lacking. To fill this gap, we developed NetFlow, an algorithm that leverages genome-scale carbon mapping to extract and quantitatively distinguish biologically relevant metabolic pathways from a given genome-scale flux prediction. NetFlow extends prior approaches by utilizing both full carbon mapping and context-specific flux predictions. Thus, NetFlow is uniquely able to quantitatively distinguish between biologically relevant pathways of carbon flow within the given flux map. NetFlow simulates 13C isotope labeling experiments to calculate the extent of carbon exchange, or carbon yield, between every metabolite in the given GSM. Based on the carbon yield, the carbon flow to or from any metabolite or between any pair of metabolites of interest can be isolated and readily visualized. The resulting pathways are much easier to interpret, which enables an in-depth mechanistic understanding of the metabolic phenotype of interest. Here, we first demonstrate NetFlow with a simple network. We then depict the utility of NetFlow on a model of central carbon metabolism in E. coli. Specifically, we isolated the production pathway for succinate synthesis in this model and the metabolic mechanism driving the predicted increase in succinate yield in a double knockout of E. coli. Finally, we describe the application of NetFlow to a GSM of lycopene-producing E. coli, which enabled the rapid identification of the mechanisms behind the measured increases in lycopene production following single, double, and triple knockouts.

NetRed, an algorithm to reduce genome-scale metabolic networks and facilitate the analysis of flux predictions.

Flux balance analysis (FBA) of large, genome-scale stoichiometric models (GSMs) is a powerful and popular method to predict cell-wide metabolic activity. FBA typically generates a flux vector containing O(1,000) fluxes. The interpretation of such a flux vector is difficult, even for expert users, because of the large size and complex topology of the underlying metabolic network. This interpretation could be simplified by condensing the network to a reduced, yet fully representative version. Toward this goal we report NetRed, an algorithm that systematically reduces a stoichiometric matrix and a corresponding flux vector to a more easily interpretable form. The reduction offered by NetRed is transparent because it relies purely on matrix algebra and not on optimization. Uniquely, it involves zero information loss; therefore, the original unreduced network can be easily recovered from the reduced network. The inputs to NetRed are (i) a stoichiometric matrix, (ii) a flux vector with numerical flux values, and (iii) a list of "protected" metabolites recommended by the user to remain in the reduced network. NetRed outputs a reduced metabolic network containing a reduced number of metabolites, of which the protected metabolites are a subset. The algorithm also generates a corresponding reduced flux vector. Due to its simplified presentation and easier interpretability, the reduced network allows the user to quickly find fluxes through metabolites and reaction modes or pathways of interest. In this manuscript, we first demonstrate NetRed on a simple network consisting of glycolysis and the pentose phosphate pathway (PPP), wherein NetRed reduced the PPP to a single net reaction. We followed this with applications of NetRed to E. coli and yeast GSMs. NetRed reduced the size of an E. coli GSM by 20- to 30-fold and enabled a comprehensive comparison of aerobic and anaerobic metabolism. The application of NetRed to a yeast GSM allowed for easy mechanistic interpretation of a double-gene knockout that rerouted flux toward dihydroartemisinic acid. When applied to an E. coli strain engineered for enhanced valine production, NetRed allowed for a holistic interpretation of the metabolic rerouting resulting from multiple genetic interventions.

Achieving a Predictive Understanding of Antimicrobial Stress Physiology through Systems Biology.

The dramatic spread and diversity of antibiotic-resistant pathogens has significantly reduced the efficacy of essentially all antibiotic classes, bringing us ever closer to a postantibiotic era. Exacerbating this issue, our understanding of the multiscale physiological impact of antimicrobial challenge on bacterial pathogens remains incomplete. Concerns over resistance and the need for new antibiotics have motivated the collection of omics measurements to provide systems-level insights into antimicrobial stress responses for nearly 20 years. Although technological advances have markedly improved the types and resolution of such measurements, continued development of mathematical frameworks aimed at providing a predictive understanding of complex antimicrobial-associated phenotypes is critical to maximize the utility of multiscale data. Here we highlight recent efforts utilizing systems biology to enhance our knowledge of antimicrobial stress physiology. We provide a brief historical perspective of antibiotic-focused omics measurements, highlight new measurement discoveries and trends, discuss examples and opportunities for integrating measurements with mathematical models, and describe future challenges for the field.

Systems biology investigation of cAMP modulation to increase SMN levels for the treatment of spinal muscular atrophy.

Spinal muscular atrophy (SMA), a leading genetic cause of infant death worldwide, is an autosomal recessive disorder caused by the loss of SMN1 (survival motor neuron 1), which encodes the protein SMN. The loss of SMN1 causes a deficiency in SMN protein levels leading to motor neuron cell death in the anterior horn of the spinal cord. SMN2, however, can also produce some functional SMN to partially compensate for loss of SMN1 in SMA suggesting increasing transcription of SMN2 as a potential therapy to treat patients with SMA. A cAMP response element was identified on the SMN2 promoter, implicating cAMP activation as a step in the transcription of SMN2. Therefore, we investigated the effects of modulating the cAMP signaling cascade on SMN production in vitro and in silico. SMA patient fibroblasts were treated with the cAMP signaling modulators rolipram, salbutamol, dbcAMP, epinephrine and forskolin. All of the modulators tested were able to increase gem formation, a marker for SMN protein in the nucleus, in a dose-dependent manner. We then derived two possible mathematical models simulating the regulation of SMN2 expression by cAMP signaling. Both models fit well with our experimental data. In silico treatment of SMA fibroblasts simultaneously with two different cAMP modulators resulted in an additive increase in gem formation. This study shows how a systems biology approach can be used to develop potential therapeutic targets for treating SMA.