Identification of functional differences in metabolic networks using comparative genomics and constraint-based models

Community-scale models of microbiomes: Articulating metabolic modelling and metagenome sequencing

Building models is essential for understanding the functions and dynamics of microbial communities. Metabolic models built on genome-scale metabolic network reconstructions (GENREs) are especially relevant as a means to decipher the complex interactions occurring among species. Model reconstruction increasingly relies on metagenomics, which permits direct characterisation of naturally occurring communities that may contain organisms that cannot be isolated or cultured. In this review, we provide an overview of the field of metabolic modelling and its increasing reliance on and synergy with metagenomics and bioinformatics. We survey the means of assigning functions and reconstructing metabolic networks from (meta-)genomes, and present the variety and mathematical fundamentals of metabolic models that foster the understanding of microbial dynamics. We emphasise the characterisation of interactions and the scaling of model construction to large communities, two important bottlenecks in the applicability of these models. We give an overview of the current state of the art in metagenome sequencing and bioinformatics analysis, focusing on the reconstruction of genomes in microbial communities. Metagenomics benefits tremendously from third-generation sequencing, and we discuss the opportunities of long-read sequencing, strain-level characterisation and eukaryotic metagenomics. We aim at providing algorithmic and mathematical support, together with tool and application resources, that permit bridging the gap between metagenomics and metabolic modelling.

A Multi-Level Systems Biology Analysis of Aldrin's Metabolic Effects on Prostate Cancer Cells

Although numerous studies support a dose-effect relationship between Endocrine disruptors (EDs) and the progression and malignancy of tumors, the impact of a chronic exposure to non-lethal concentrations of EDs in cancer remains unknown. More specifically, a number of studies have reported the impact of Aldrin on a variety of cancer types, including prostate cancer. In previous studies, we demonstrated the induction of the malignant phenotype in DU145 prostate cancer (PCa) cells after a chronic exposure to Aldrin (an ED). Proteins are pivotal in the regulation and control of a variety of cellular processes. However, the mechanisms responsible for the impact of ED on PCa and the role of proteins in this process are not yet well understood. Here, two complementary computational approaches have been employed to investigate the molecular processes underlying the acquisition of malignancy in prostate cancer. First, the metabolic reprogramming associated with the chronic exposure to Aldrin in DU145 cells was studied by integrating transcriptomics and metabolomics via constraint-based metabolic modeling. Second, gene set enrichment analysis was applied to determine (i) altered regulatory pathways and (ii) the correlation between changes in the transcriptomic profile of Aldrin-exposed cells and tumor progression in various types of cancer. Experimental validation confirmed predictions revealing a disruption in metabolic and regulatory pathways. This alteration results in the modification of protein levels crucial in regulating triacylglyceride/cholesterol, linked to the malignant phenotype observed in Aldrin-exposed cells.

Music of metagenomics-a review of its applications, analysis pipeline, and associated tools

This humble effort highlights the intricate details of metagenomics in a simple, poetic, and rhythmic way. The paper enforces the significance of the research area, provides details about major analytical methods, examines the taxonomy and assembly of genomes, emphasizes some tools, and concludes by celebrating the richness of the ecosystem populated by the "metagenome."

Prospects for multi-omics in the microbial ecology of water engineering

Advances in high-throughput sequencing technologies and bioinformatics approaches over almost the last three decades have substantially increased our ability to explore microorganisms and their functions - including those that have yet to be cultivated in pure isolation. Genome-resolved metagenomic approaches have enabled linking powerful functional predictions to specific taxonomical groups with increasing fidelity. Additionally, related developments in both whole community gene expression surveys and metabolite profiling have permitted for direct surveys of community-scale functions in specific environmental settings. These advances have allowed for a shift in microbiome science away from descriptive studies and towards mechanistic and predictive frameworks for designing and harnessing microbial communities for desired beneficial outcomes. Water engineers, microbiologists, and microbial ecologists studying activated sludge, anaerobic digestion, and drinking water distribution systems have applied various (meta)omics techniques for connecting microbial community dynamics and physiologies to overall process parameters and system performance. However, the rapid pace at which new omics-based approaches are developed can appear daunting to those looking to apply these state-of-the-art practices for the first time. Here, we review how modern genome-resolved metagenomic approaches have been applied to a variety of water engineering applications from lab-scale bioreactors to full-scale systems. We describe integrated omics analysis across engineered water systems and the foundations for pairing these insights with modeling approaches. Lastly, we summarize emerging omics-based technologies that we believe will be powerful tools for water engineering applications. Overall, we provide a framework for microbial ecologists specializing in water engineering to apply cutting-edge omics approaches to their research questions to achieve novel functional insights. Successful adoption of predictive frameworks in engineered water systems could enable more economically and environmentally sustainable bioprocesses as demand for water and energy resources increases.

Towards model-driven characterization and manipulation of plant lipid metabolism

Plant lipids have versatile applications and provide essential fatty acids in human diet. Therefore, there has been a growing interest to better characterize the genetic basis, regulatory networks, and metabolic pathways that shape lipid quantity and composition. Addressing these issues is challenging due to context-specificity of lipid metabolism integrating environmental, developmental, and tissue-specific cues. Here we systematically review the known metabolic pathways and regulatory interactions that modulate the levels of storage lipids in oilseeds. We argue that the current understanding of lipid metabolism provides the basis for its study in the context of genome-wide plant metabolic networks with the help of approaches from constraint-based modeling and metabolic flux analysis. The focus is on providing a comprehensive summary of the state-of-the-art of modeling plant lipid metabolic pathways, which we then contrast with the existing modeling efforts in yeast and microalgae. We then point out the gaps in knowledge of lipid metabolism, and enumerate the recent advances of using genome-wide association and quantitative trait loci mapping studies to unravel the genetic regulations of lipid metabolism. Finally, we offer a perspective on how advances in the constraint-based modeling framework can propel further characterization of plant lipid metabolism and its rational manipulation.

Biochemical Characteristics and a Genome-Scale Metabolic Model of an Indian Euryhaline Cyanobacterium with High Polyglucan Content

Marine cyanobacteria are promising microbes to capture and convert atmospheric CO₂ and light into biomass and valuable industrial bio-products. Yet, reports on metabolic characteristics of non-model cyanobacteria are scarce. In this report, we show that an Indian euryhaline Synechococcus sp. BDU 130192 has biomass accumulation comparable to a model marine cyanobacterium and contains approximately double the amount of total carbohydrates, but significantly lower protein levels compared to Synechococcus sp. PCC 7002 cells. Based on its annotated chromosomal genome sequence, we present a genome scale metabolic model (GSMM) of this cyanobacterium, which we have named as iSyn706. The model includes 706 genes, 908 reactions, and 900 metabolites. The difference in the flux balance analysis (FBA) predicted flux distributions between Synechococcus sp. PCC 7002 and Synechococcus sp. BDU130192 strains mimicked the differences in their biomass compositions. Model-predicted oxygen evolution rate for Synechococcus sp. BDU130192 was found to be close to the experimentally-measured value. The model was analyzed to determine the potential of the strain for the production of various industrially-useful products without affecting growth significantly. This model will be helpful to researchers interested in understanding the metabolism as well as to design metabolic engineering strategies for the production of industrially-relevant compounds.

Integrating proteomic or transcriptomic data into metabolic models using linear bound flux balance analysis

Motivation

Transcriptomics and proteomics data have been integrated into constraint-based models to influence flux predictions. However, it has been reported recently for Escherichia coli and Saccharomyces cerevisiae, that model predictions from parsimonious flux balance analysis (pFBA), which does not use expression data, are as good or better than predictions from various algorithms that integrate transcriptomics or proteomics data into constraint-based models.

Results

In this paper, we describe a novel constraint-based method called Linear Bound Flux Balance Analysis (LBFBA), which uses expression data (either transcriptomic or proteomic) to predict metabolic fluxes. The method uses expression data to place soft constraints on individual fluxes, which can be violated. Parameters in the soft constraints are first estimated from a training expression and flux dataset before being used to predict fluxes from expression data in other conditions. We applied LBFBA to E.coli and S.cerevisiae datasets and found that LBFBA predictions were more accurate than pFBA predictions, with average normalized errors roughly half of those from pFBA. For the first time, we demonstrate a computational method that integrates expression data into constraint-based models and improves quantitative flux predictions over pFBA.

Availability and implementation

Code is available in the Supplementary data available at Bioinformatics online.

Supplementary information

Supplementary data are available at Bioinformatics online.

New Applications of Synthetic Biology Tools for Cyanobacterial Metabolic Engineering

Cyanobacteria are promising microorganisms for sustainable biotechnologies, yet unlocking their potential requires radical re-engineering and application of cutting-edge synthetic biology techniques. In recent years, the available devices and strategies for modifying cyanobacteria have been increasing, including advances in the design of genetic promoters, ribosome binding sites, riboswitches, reporter proteins, modular vector systems, and markerless selection systems. Because of these new toolkits, cyanobacteria have been successfully engineered to express heterologous pathways for the production of a wide variety of valuable compounds. Cyanobacterial strains with the potential to be used in real-world applications will require the refinement of genetic circuits used to express the heterologous pathways and development of accurate models that predict how these pathways can be best integrated into the larger cellular metabolic network. Herein, we review advances that have been made to translate synthetic biology tools into cyanobacterial model organisms and summarize experimental and in silico strategies that have been employed to increase their bioproduction potential. Despite the advances in synthetic biology and metabolic engineering during the last years, it is clear that still further improvements are required if cyanobacteria are to be competitive with heterotrophic microorganisms for the bioproduction of added-value compounds.

Advances in metabolic modeling of oleaginous microalgae

Production of biofuels and bioenergy precursors by phototrophic microorganisms, such as microalgae and cyanobacteria, is a promising alternative to conventional fuels obtained from non-renewable resources. Several species of microalgae have been investigated as potential candidates for the production of biofuels, for the most part due to their exceptional metabolic capability to accumulate large quantities of lipids. Constraint-based modeling, a systems biology approach that accurately predicts the metabolic phenotype of phototrophs, has been deployed to identify suitable culture conditions as well as to explore genetic enhancement strategies for bioproduction. Core metabolic models were employed to gain insight into the central carbon metabolism in photosynthetic microorganisms. More recently, comprehensive genome-scale models, including organelle-specific information at high resolution, have been developed to gain new insight into the metabolism of phototrophic cell factories. Here, we review the current state of the art of constraint-based modeling and computational method development and discuss how advanced models led to increased prediction accuracy and thus improved lipid production in microalgae.

Diversity, Function, and Application of Clostridium in Chinese Strong Flavor Baijiu Ecosystem: A Review

Baijiu is a Chinese traditional distilled liquor with an annual yield over 13.12 million tons. Strong flavor baijiu (SFB) also called Luzhou-flavor liquor, takes account for > 70% of the total baijiu produced. SFB is produced by an open solid fermentation process with a complex microbial ecosystem. Clostridium is one of the most important microorganisms for the formation of the main flavor compounds of SFB, such as ethyl caproate. In this paper, we review current research progress on the Clostridium in the SFB ecosystem, focusing on the species diversity, physiological and metabolic features along with interspecies interactions. Systems biology approaches for the study of Clostridium from SFB ecosystems were discussed and explored. Furthermore, current applications of Clostridium in SFB production were discussed.

PRACTICAL APPLICATION

Strong flavor baijiu (SFB) accounts for more than 70% of total yield of Chinese baijiu, which exists for hundreds of years. Clostridium is common in SFB ecosystem and identified to be one of main contributors of flavor compounds in SFB. Study on the Clostridium from SFB ecosystem is not only helpful for the understanding of flavor compounds formation mechanism, but also the improvement of SFB quality. This study focuses on the current researches on the Clostridium species in SFB ecosystem, including the species diversity, physiological and metabolic features, and applications.

Model-driven discovery of long-chain fatty acid metabolic reprogramming in heterogeneous prostate cancer cells

Epithelial-mesenchymal-transition promotes intra-tumoral heterogeneity, by enhancing tumor cell invasiveness and promoting drug resistance. We integrated transcriptomic data for two clonal subpopulations from a prostate cancer cell line (PC-3) into a genome-scale metabolic network model to explore their metabolic differences and potential vulnerabilities. In this dual cell model, PC-3/S cells express Epithelial-mesenchymal-transition markers and display high invasiveness and low metastatic potential, while PC-3/M cells present the opposite phenotype and higher proliferative rate. Model-driven analysis and experimental validations unveiled a marked metabolic reprogramming in long-chain fatty acids metabolism. While PC-3/M cells showed an enhanced entry of long-chain fatty acids into the mitochondria, PC-3/S cells used long-chain fatty acids as precursors of eicosanoid metabolism. We suggest that this metabolic reprogramming endows PC-3/M cells with augmented energy metabolism for fast proliferation and PC-3/S cells with increased eicosanoid production impacting angiogenesis, cell adhesion and invasion. PC-3/S metabolism also promotes the accumulation of docosahexaenoic acid, a long-chain fatty acid with antiproliferative effects. The potential therapeutic significance of our model was supported by a differential sensitivity of PC-3/M cells to etomoxir, an inhibitor of long-chain fatty acid transport to the mitochondria.

The coexistence within the same tumor of a variety of subpopulations, featuring different phenotypes (intra-tumoral heterogeneity) represents a challenge for diagnosis, prognosis and targeted therapies. In this work, we have explored the metabolic differences underlying tumor heterogeneity by building cell-type-specific genome-scale metabolic models that integrate transcriptome and metabolome data of two clonal subpopulations derived from the same prostate cancer cell line (PC-3). These subpopulations display either highly proliferative, cancer stem cell (PC-3/M) or highly invasive, epithelial-mesenchymal-transition-like phenotypes (PC-3/S). Our model-driven analysis and experimental validations have unveiled a differential utilization of the long-chain fatty acids pool in both subpopulations. More specifically, our findings show an enhanced entry of long-chain fatty acids into the mitochondria in PC-3/M cells, while in PC-3/S cells, long-chain fatty acids are used as precursors of eicosanoid metabolism. The different utilization of long-chain fatty acids between subpopulations endows PC-3/M cells with a highly proliferative phenotype while enhances PC-3/S invasive phenotype. The present work provides a tool to unveil key metabolic nodes associated with tumor heterogeneity and highlights potential subpopulation-specific targets with important therapeutic implications.

Rewiring of Cyanobacterial Metabolism for Hydrogen Production: Synthetic Biology Approaches and Challenges

With the demand for renewable energy growing, hydrogen (H₂) is becoming an attractive energy carrier. Developing H₂ production technologies with near-net zero carbon emissions is a major challenge for the "H₂ economy." Certain cyanobacteria inherently possess enzymes, nitrogenases, and bidirectional hydrogenases that are capable of H₂ evolution using sunlight, making them ideal cell factories for photocatalytic conversion of water to H₂. With the advances in synthetic biology, cyanobacteria are currently being developed as a "plug and play" chassis to produce H₂. This chapter describes the metabolic pathways involved and the theoretical limits to cyanobacterial H₂ production and summarizes the metabolic engineering technologies pursued.

Constraint-based metabolic modelling of marine microbes and communities

Constraint-based metabolic modelling (CBMM) consists in the use of computational methods and tools to perform genome-scale simulations and predict metabolic features at the whole cellular level. This approach is rapidly expanding in microbiology, as it combines reliable predictive abilities with conceptually and technically simple frameworks. Among the possible outcomes of CBMM, the capability to i) guide a focused planning of metabolic engineering experiments and ii) provide a system-level understanding of (single or community-level) microbial metabolic circuits also represent primary aims in present-day marine microbiology. In this work we briefly introduce the theoretical formulation behind CBMM and then review the most recent and effective case studies of CBMM of marine microbes and communities. Also, the emerging challenges and possibilities in the use of such methodologies in the context of marine microbiology/biotechnology are discussed. As the potential applications of CBMM have a very broad range, the topics presented in this review span over a large plethora of fields such as ecology, biotechnology and evolution.

Community metabolic modeling approaches to understanding the gut microbiome: Bridging biochemistry and ecology

Interest in the human microbiome is at an all time high. The number of human microbiome studies is growing exponentially, as are reported associations between microbial communities and disease. However, we have not been able to translate the ever-growing amount of microbiome sequence data into better health. To do this, we need a practical means of transforming a disease-associated microbiome into a health-associated microbiome. This will require a framework that can be used to generate predictions about community dynamics within the microbiome under different conditions, predictions that can be tested and validated. In this review, using the gut microbiome to illustrate, we describe two classes of model that are currently being used to generate predictions about microbial community dynamics: ecological models and metabolic models. We outline the strengths and weaknesses of each approach and discuss the insights into the gut microbiome that have emerged from modeling thus far. We then argue that the two approaches can be combined to yield a community metabolic model, which will supply the framework needed to move from high-throughput omics data to testable predictions about how prebiotic, probiotic, and nutritional interventions affect the microbiome. We are confident that with a suitable model, researchers and clinicians will be able to harness the stream of sequence data and begin designing strategies to make targeted alterations to the microbiome and improve health.

Toward Multiscale Models of Cyanobacterial Growth: A Modular Approach

Oxygenic photosynthesis dominates global primary productivity ever since its evolution more than three billion years ago. While many aspects of phototrophic growth are well understood, it remains a considerable challenge to elucidate the manifold dependencies and interconnections between the diverse cellular processes that together facilitate the synthesis of new cells. Phototrophic growth involves the coordinated action of several layers of cellular functioning, ranging from the photosynthetic light reactions and the electron transport chain, to carbon-concentrating mechanisms and the assimilation of inorganic carbon. It requires the synthesis of new building blocks by cellular metabolism, protection against excessive light, as well as diurnal regulation by a circadian clock and the orchestration of gene expression and cell division. Computational modeling allows us to quantitatively describe these cellular functions and processes relevant for phototrophic growth. As yet, however, computational models are mostly confined to the inner workings of individual cellular processes, rather than describing the manifold interactions between them in the context of a living cell. Using cyanobacteria as model organisms, this contribution seeks to summarize existing computational models that are relevant to describe phototrophic growth and seeks to outline their interactions and dependencies. Our ultimate aim is to understand cellular functioning and growth as the outcome of a coordinated operation of diverse yet interconnected cellular processes.

Flux balance analysis of photoautotrophic metabolism: Uncovering new biological details of subsystems involved in cyanobacterial photosynthesis

We have constructed and experimentally tested a comprehensive genome-scale model of photoautotrophic growth, denoted iSyp821, for the cyanobacterium Synechococcus sp. PCC 7002. iSyp821 incorporates a variable biomass objective function (vBOF), in which stoichiometries of the major biomass components vary according to light intensity. The vBOF was constrained to fit the measured cellular carbohydrate/protein content under different light intensities. iSyp821 provides rigorous agreement with experimentally measured cell growth rates and inorganic carbon uptake rates as a function of light intensity. iSyp821 predicts two observed metabolic transitions that occur as light intensity increases: 1) from PSI-cyclic to linear electron flow (greater redox energy), and 2) from carbon allocation as proteins (growth) to carbohydrates (energy storage) mode. iSyp821 predicts photoautotrophic carbon flux into 1) a hybrid gluconeogenesis-pentose phosphate (PP) pathway that produces glycogen by an alternative pathway than conventional gluconeogenesis, and 2) the photorespiration pathway to synthesize the essential amino acid, glycine. Quantitative fluxes through both pathways were verified experimentally by following the kinetics of formation of ¹³C metabolites from ¹³CO₂ fixation. iSyp821 was modified to include changes in gene products (enzymes) from experimentally measured transcriptomic data and applied to estimate changes in concentrations of metabolites arising from nutrient stress. Using this strategy, we found that iSyp821 correctly predicts the observed redistribution pattern of carbon products under nitrogen depletion, including decreased rates of CO₂ uptake, amino acid synthesis, and increased rates of glycogen and lipid synthesis.

Cyanobacterial Biofuels: Strategies and Developments on Network and Modeling

Cyanobacteria, the phototrophic microorganisms, have attracted much attention recently as a promising source for environmentally sustainable biofuels production. However, barriers for commercial markets of cyanobacteria-based biofuels concern the economic feasibility. Miscellaneous strategies for improving the production performance of cyanobacteria have thus been developed. Among these, the simple ad hoc strategies resulting in failure to optimize fully cell growth coupled with desired product yield are explored. With the advancement of genomics and systems biology, a new paradigm toward systems metabolic engineering has been recognized. In particular, a genome-scale metabolic network reconstruction and modeling is a crucial systems-based tool for whole-cell-wide investigation and prediction. In this review, the cyanobacterial genome-scale metabolic models, which offer a system-level understanding of cyanobacterial metabolism, are described. The main process of metabolic network reconstruction and modeling of cyanobacteria are summarized. Strategies and developments on genome-scale network and modeling through the systems metabolic engineering approach are advanced and employed for efficient cyanobacterial-based biofuels production.

Integrated in silico Analyses of Regulatory and Metabolic Networks of Synechococcus sp. PCC 7002 Reveal Relationships between Gene Centrality and Essentiality

Cyanobacteria dynamically relay environmental inputs to intracellular adaptations through a coordinated adjustment of photosynthetic efficiency and carbon processing rates. The output of such adaptations is reflected through changes in transcriptional patterns and metabolic flux distributions that ultimately define growth strategy. To address interrelationships between metabolism and regulation, we performed integrative analyses of metabolic and gene co-expression networks in a model cyanobacterium, Synechococcus sp. PCC 7002. Centrality analyses using the gene co-expression network identified a set of key genes, which were defined here as "topologically important." Parallel in silico gene knock-out simulations, using the genome-scale metabolic network, classified what we termed as "functionally important" genes, deletion of which affected growth or metabolism. A strong positive correlation was observed between topologically and functionally important genes. Functionally important genes exhibited variable levels of topological centrality; however, the majority of topologically central genes were found to be functionally essential for growth. Subsequent functional enrichment analysis revealed that both functionally and topologically important genes in Synechococcus sp. PCC 7002 are predominantly associated with translation and energy metabolism, two cellular processes critical for growth. This research demonstrates how synergistic network-level analyses can be used for reconciliation of metabolic and gene expression data to uncover fundamental biological principles.

Expanding Metabolic Engineering Algorithms Using Feasible Space and Shadow Price Constraint Modules

While numerous computational methods have been developed that use genome-scale models to propose mutants for the purpose of metabolic engineering, they generally compare mutants based on a single criteria (e.g., production rate at a mutant׳s maximum growth rate). As such, these approaches remain limited in their ability to include multiple complex engineering constraints. To address this shortcoming, we have developed feasible space and shadow price constraint (FaceCon and ShadowCon) modules that can be added to existing mixed integer linear adaptive evolution metabolic engineering algorithms, such as OptKnock and OptORF. These modules allow strain designs to be identified amongst a set of multiple metabolic engineering algorithm solutions that are capable of high chemical production while also satisfying additional design criteria. We describe the various module implementations and their potential applications to the field of metabolic engineering. We then incorporated these modules into the OptORF metabolic engineering algorithm. Using an Escherichia coli genome-scale model (iJO1366), we generated different strain designs for the anaerobic production of ethanol from glucose, thus demonstrating the tractability and potential utility of these modules in metabolic engineering algorithms.

•

Added modules to impose multiple design criteria for engineering algorithms.

•

Examples are provided to eliminate by-product secretion.

•

Examples are provided to control coupling between product and biomass formation.

•

Modules are tractable for genome-scale design.

Genome-scale reconstruction of metabolic networks of Lactobacillus casei ATCC 334 and 12A

Lactobacillus casei strains are widely used in industry and the utility of this organism in these industrial applications is strain dependent. Hence, tools capable of predicting strain specific phenotypes would have utility in the selection of strains for specific industrial processes. Genome-scale metabolic models can be utilized to better understand genotype-phenotype relationships and to compare different organisms. To assist in the selection and development of strains with enhanced industrial utility, genome-scale models for L. casei ATCC 334, a well characterized strain, and strain 12A, a corn silage isolate, were constructed. Draft models were generated from RAST genome annotations using the Model SEED database and refined by evaluating ATP generating cycles, mass-and-charge-balances of reactions, and growth phenotypes. After the validation process was finished, we compared the metabolic networks of these two strains to identify metabolic, genetic and ortholog differences that may lead to different phenotypic behaviors. We conclude that the metabolic capabilities of the two networks are highly similar. The L. casei ATCC 334 model accounts for 1,040 reactions, 959 metabolites and 548 genes, while the L. casei 12A model accounts for 1,076 reactions, 979 metabolites and 640 genes. The developed L. casei ATCC 334 and 12A metabolic models will enable better understanding of the physiology of these organisms and be valuable tools in the development and selection of strains with enhanced utility in a variety of industrial applications.

Generation and Evaluation of a Genome-Scale Metabolic Network Model of Synechococcus elongatus PCC7942

The reconstruction of genome-scale metabolic models and their applications represent a great advantage of systems biology. Through their use as metabolic flux simulation models, production of industrially-interesting metabolites can be predicted. Due to the growing number of studies of metabolic models driven by the increasing genomic sequencing projects, it is important to conceptualize steps of reconstruction and analysis. We have focused our work in the cyanobacterium Synechococcus elongatus PCC7942, for which several analyses and insights are unveiled. A comprehensive approach has been used, which can be of interest to lead the process of manual curation and genome-scale metabolic analysis. The final model, iSyf715 includes 851 reactions and 838 metabolites. A biomass equation, which encompasses elementary building blocks to allow cell growth, is also included. The applicability of the model is finally demonstrated by simulating autotrophic growth conditions of Synechococcus elongatus PCC7942.

Comparisons of Shewanella strains based on genome annotations, modeling, and experiments

BACKGROUND

Shewanella is a genus of facultatively anaerobic, Gram-negative bacteria that have highly adaptable metabolism which allows them to thrive in diverse environments. This quality makes them an attractive bacterial target for research in bioremediation and microbial fuel cell applications. Constraint-based modeling is a useful tool for helping researchers gain insights into the metabolic capabilities of these bacteria. However, Shewanella oneidensis MR-1 is the only strain with a genome-scale metabolic model constructed out of 21 sequenced Shewanella strains.

RESULTS

In this work, we updated the model for Shewanella oneidensis MR-1 and constructed metabolic models for three other strains, namely Shewanella sp. MR-4, Shewanella sp. W3-18-1, and Shewanella denitrificans OS217 which span the genus based on the number of genes lost in comparison to MR-1. We also constructed a Shewanella core model that contains the genes shared by all 21 sequenced strains and a few non-conserved genes associated with essential reactions. Model comparisons between the five constructed models were done at two levels - for wildtype strains under different growth conditions and for knockout mutants under the same growth condition. In the first level, growth/no-growth phenotypes were predicted by the models on various carbon sources and electron acceptors. Cluster analysis of these results revealed that the MR-1 model is most similar to the W3-18-1 model, followed by the MR-4 and OS217 models when considering predicted growth phenotypes. However, a cluster analysis done based on metabolic gene content revealed that the MR-4 and W3-18-1 models are the most similar, with the MR-1 and OS217 models being more distinct from these latter two strains. As a second level of comparison, we identified differences in reaction and gene content which give rise to different functional predictions of single and double gene knockout mutants using Comparison of Networks by Gene Alignment (CONGA). Here, we showed how CONGA can be used to find biomass, metabolic, and genetic differences between models.

CONCLUSIONS

We developed four strain-specific models and a general core model that can be used to do various in silico studies of Shewanella metabolism. The developed models provide a platform for a systematic investigation of Shewanella metabolism to aid researchers using Shewanella in various biotechnology applications.

Rapid construction of metabolic models for a family of Cyanobacteria using a multiple source annotation workflow

Background

Cyanobacteria are photoautotrophic prokaryotes that exhibit robust growth under diverse environmental conditions with minimal nutritional requirements. They can use solar energy to convert CO₂ and other reduced carbon sources into biofuels and chemical products. The genus Cyanothece includes unicellular nitrogen-fixing cyanobacteria that have been shown to offer high levels of hydrogen production and nitrogen fixation. The reconstruction of quality genome-scale metabolic models for organisms with limited annotation resources remains a challenging task.

Results

Here we reconstruct and subsequently analyze and compare the metabolism of five Cyanothece strains, namely Cyanothece sp. PCC 7424, 7425, 7822, 8801 and 8802, as the genome-scale metabolic reconstructions iCyc792, iCyn731, iCyj826, iCyp752, and iCyh755 respectively. We compare these phylogenetically related Cyanothece strains to assess their bio-production potential. A systematic workflow is introduced for integrating and prioritizing annotation information from the Universal Protein Resource (Uniprot), NCBI Protein Clusters, and the Rapid Annotations using Subsystems Technology (RAST) method. The genome-scale metabolic models include fully traced photosynthesis reactions and respiratory chains, as well as balanced reactions and GPR associations. Metabolic differences between the organisms are highlighted such as the non-fermentative pathway for alcohol production found in only Cyanothece 7424, 8801, and 8802.

Conclusions

Our development workflow provides a path for constructing models using information from curated models of related organisms and reviewed gene annotations. This effort lays the foundation for the expedient construction of curated metabolic models for organisms that, while not being the target of comprehensive research, have a sequenced genome and are related to an organism with a curated metabolic model. Organism-specific models, such as the five presented in this paper, can be used to identify optimal genetic manipulations for targeted metabolite overproduction as well as to investigate the biology of diverse organisms.

Pathogenicity phenomena in three model systems: from network mining to emerging system-level properties

Understanding the interconnections of microbial pathogenicity phenomena, such as biofilm formation, quorum sensing and antimicrobial resistance, is a tremendous open challenge for biomedical research. Progress made by wet-lab researchers and bioinformaticians in understanding the underlying regulatory phenomena has been significant, with converging evidence from multiple high-throughput technologies. Notably, network reconstructions are already of considerable size and quality, tackling both intracellular regulation and signal mediation in microbial infection. Therefore, it stands to reason that in silico investigations would play a more active part in this research. Drug target identification and drug repurposing could take much advantage of the ability to simulate pathogen regulatory systems, host-pathogen interactions and pathogen cross-talking. Here, we review the bioinformatics resources and tools available for the study of the gram-negative bacterium Pseudomonas aeruginosa, the gram-positive bacterium Staphylococcus aureus and the fungal species Candida albicans. The choice of these three microorganisms fits the rationale of the review converging into pathogens of great clinical importance, which thrive in biofilm consortia and manifest growing antimicrobial resistance.

Synechocystis sp. PCC6803 metabolic models for the enhanced production of hydrogen

In the present economy, difficulties to access energy sources are real drawbacks to maintain our current lifestyle. In fact, increasing interests have been gathered around efficient strategies to use energy sources that do not generate high CO2 titers. Thus, science-funding agencies have invested more resources into research on hydrogen among other biofuels as interesting energy vectors. This article reviews present energy challenges and frames it into the present fuel usage landscape. Different strategies for hydrogen production are explained and evaluated. Focus is on biological hydrogen production; fermentation and photon-fuelled hydrogen production are compared. Mathematical models in biology can be used to assess, explore and design production strategies for industrially relevant metabolites, such as biofuels. We assess the diverse construction and uses of genome-scale metabolic models of cyanobacterium Synechocystis sp. PCC6803 to efficiently obtain biofuels. This organism has been studied as a potential photon-fuelled production platform for its ability to grow from carbon dioxide, water and photons, on simple culture media. Finally, we review studies that propose production strategies to weigh this organism's viability as a biofuel production platform. Overall, the work presented in this review unveils the industrial capabilities of cyanobacterium Synechocystis sp. PCC6803 to evolve interesting metabolites as a clean biofuel production platform.

Synthetic biology of cyanobacteria: unique challenges and opportunities

Photosynthetic organisms, and especially cyanobacteria, hold great promise as sources of renewably-produced fuels, bulk and specialty chemicals, and nutritional products. Synthetic biology tools can help unlock cyanobacteria's potential for these functions, but unfortunately tool development for these organisms has lagged behind that for S. cerevisiae and E. coli. While these organisms may in many cases be more difficult to work with as “chassis” strains for synthetic biology than certain heterotrophs, the unique advantages of autotrophs in biotechnology applications as well as the scientific importance of improved understanding of photosynthesis warrant the development of these systems into something akin to a “green E. coli.” In this review, we highlight unique challenges and opportunities for development of synthetic biology approaches in cyanobacteria. We review classical and recently developed methods for constructing targeted mutants in various cyanobacterial strains, and offer perspective on what genetic tools might most greatly expand the ability to engineer new functions in such strains. Similarly, we review what genetic parts are most needed for the development of cyanobacterial synthetic biology. Finally, we highlight recent methods to construct genome-scale models of cyanobacterial metabolism and to use those models to measure properties of autotrophic metabolism. Throughout this paper, we discuss some of the unique challenges of a diurnal, autotrophic lifestyle along with how the development of synthetic biology and biotechnology in cyanobacteria must fit within those constraints.

Quantitative assessment of thermodynamic constraints on the solution space of genome-scale metabolic models

Constraint-based methods provide powerful computational techniques to allow understanding and prediction of cellular behavior. These methods rely on physiochemical constraints to eliminate infeasible behaviors from the space of available behaviors. One such constraint is thermodynamic feasibility, the requirement that intracellular flux distributions obey the laws of thermodynamics. The past decade has seen several constraint-based methods that interpret this constraint in different ways, including those that are limited to small networks, rely on predefined reaction directions, and/or neglect the relationship between reaction free energies and metabolite concentrations. In this work, we utilize one such approach, thermodynamics-based metabolic flux analysis (TMFA), to make genome-scale, quantitative predictions about metabolite concentrations and reaction free energies in the absence of prior knowledge of reaction directions, while accounting for uncertainties in thermodynamic estimates. We applied TMFA to a genome-scale network reconstruction of Escherichia coli and examined the effect of thermodynamic constraints on the flux space. We also assessed the predictive performance of TMFA against gene essentiality and quantitative metabolomics data, under both aerobic and anaerobic, and optimal and suboptimal growth conditions. Based on these results, we propose that TMFA is a useful tool for validating phenotypes and generating hypotheses, and that additional types of data and constraints can improve predictions of metabolite concentrations.

Flux balance analysis of cyanobacterial metabolism: the metabolic network of Synechocystis sp. PCC 6803

Cyanobacteria are versatile unicellular phototrophic microorganisms that are highly abundant in many environments. Owing to their capability to utilize solar energy and atmospheric carbon dioxide for growth, cyanobacteria are increasingly recognized as a prolific resource for the synthesis of valuable chemicals and various biofuels. To fully harness the metabolic capabilities of cyanobacteria necessitates an in-depth understanding of the metabolic interconversions taking place during phototrophic growth, as provided by genome-scale reconstructions of microbial organisms. Here we present an extended reconstruction and analysis of the metabolic network of the unicellular cyanobacterium Synechocystis sp. PCC 6803. Building upon several recent reconstructions of cyanobacterial metabolism, unclear reaction steps are experimentally validated and the functional consequences of unknown or dissenting pathway topologies are discussed. The updated model integrates novel results with respect to the cyanobacterial TCA cycle, an alleged glyoxylate shunt, and the role of photorespiration in cellular growth. Going beyond conventional flux-balance analysis, we extend the computational analysis to diurnal light/dark cycles of cyanobacterial metabolism.

Computational evaluation of Synechococcus sp. PCC 7002 metabolism for chemical production

Cyanobacteria are ideal metabolic engineering platforms for carbon-neutral biotechnology because they directly convert CO2 to a range of valuable products. In this study, we present a computational assessment of biochemical production in Synechococcus sp. PCC 7002 (Synechococcus 7002), a fast growing cyanobacterium whose genome has been sequenced, and for which genetic modification methods have been developed. We evaluated the maximum theoretical yields (mol product per mol CO2 or mol photon) of producing various chemicals under photoautotrophic and dark conditions using a genome-scale metabolic model of Synechococcus 7002. We found that the yields were lower under dark conditions, compared to photoautotrophic conditions, due to the limited amount of energy and reductant generated from glycogen. We also examined the effects of photon and CO2 limitations on chemical production under photoautotrophic conditions. In addition, using various computational methods such as minimization of metabolic adjustment (MOMA), relative metabolic change (RELATCH), and OptORF, we identified gene-knockout mutants that are predicted to improve chemical production under photoautotrophic and/or dark anoxic conditions. These computational results are useful for metabolic engineering of cyanobacteria to synthesize value-added products.

FOCAL: an experimental design tool for systematizing metabolic discoveries and model development

Current computational tools can generate and improve genome-scale models based on existing data; however, for many organisms, the data needed to test and refine such models are not available. To facilitate model development, we created the forced coupling algorithm, FOCAL, to identify genetic and environmental conditions such that a reaction becomes essential for an experimentally measurable phenotype. This reaction's conditional essentiality can then be tested experimentally to evaluate whether network connections occur or to create strains with desirable phenotypes. FOCAL allows network connections to be queried, which improves our understanding of metabolism and accuracy of developed models.