Genome-Scale Modeling of Light-Driven Reductant Partitioning and Carbon Fluxes in Diazotrophic Unicellular Cyanobacterium Cyanothece sp. ATCC 51142

Genome-scale metabolic models have proven useful for answering fundamental questions about metabolic capabilities of a variety of microorganisms, as well as informing their metabolic engineering. However, only a few models are available for oxygenic photosynthetic microorganisms, particularly in cyanobacteria in which photosynthetic and respiratory electron transport chains (ETC) share components. We addressed the complexity of cyanobacterial ETC by developing a genome-scale model for the diazotrophic cyanobacterium, Cyanothece sp. ATCC 51142. The resulting metabolic reconstruction, iCce806, consists of 806 genes associated with 667 metabolic reactions and includes a detailed representation of the ETC and a biomass equation based on experimental measurements. Both computational and experimental approaches were used to investigate light-driven metabolism in Cyanothece sp. ATCC 51142, with a particular focus on reductant production and partitioning within the ETC. The simulation results suggest that growth and metabolic flux distributions are substantially impacted by the relative amounts of light going into the individual photosystems. When growth is limited by the flux through photosystem I, terminal respiratory oxidases are predicted to be an important mechanism for removing excess reductant. Similarly, under photosystem II flux limitation, excess electron carriers must be removed via cyclic electron transport. Furthermore, in silico calculations were in good quantitative agreement with the measured growth rates whereas predictions of reaction usage were qualitatively consistent with protein and mRNA expression data, which we used to further improve the resolution of intracellular flux values

Identification of Functional Differences in Metabolic Networks Using Comparative Genomics and Constraint-Based Models

Genome-scale network reconstructions are useful tools for understanding cellular metabolism, and comparisons of such reconstructions can provide insight into metabolic differences between organisms. Recent efforts toward comparing genome-scale models have focused primarily on aligning metabolic networks at the reaction level and then looking at differences and similarities in reaction and gene content. However, these reaction comparison approaches are time-consuming and do not identify the effect network differences have on the functional states of the network. We have developed a bilevel mixed-integer programming approach, CONGA, to identify functional differences between metabolic networks by comparing network reconstructions aligned at the gene level. We first identify orthologous genes across two reconstructions and then use CONGA to identify conditions under which differences in gene content give rise to differences in metabolic capabilities. By seeking genes whose deletion in one or both models disproportionately changes flux through a selected reaction (e.g., growth or by-product secretion) in one model over another, we are able to identify structural metabolic network differences enabling unique metabolic capabilities. Using CONGA, we explore functional differences between two metabolic reconstructions of Escherichia coli and identify a set of reactions responsible for chemical production differences between the two models. We also use this approach to aid in the development of a genome-scale model of Synechococcus sp. PCC 7002. Finally, we propose potential antimicrobial targets in Mycobacterium tuberculosis and Staphylococcus aureus based on differences in their metabolic capabilities. Through these examples, we demonstrate that a gene-centric approach to comparing metabolic networks allows for a rapid comparison of metabolic models at a functional level. Using CONGA, we can identify differences in reaction and gene content which give rise to different functional predictions. Because CONGA provides a general framework, it can be applied to find functional differences across models and biological systems beyond those presented here