Saeed Shoaie

The RAVEN Toolbox and Its Use for Generating a Genome-scale Metabolic Model for Penicillium chrysogenum

We present the RAVEN (Reconstruction, Analysis and Visualization of Metabolic Networks) Toolbox: a software suite that allows for semi-automated reconstruction of genome-scale models. It makes use of published models and/or the KEGG database, coupled with extensive gap-filling and quality control features. The software suite also contains methods for visualizing simulation results and omics data, as well as a range of methods for performing simulations and analyzing the results. The software is a useful tool for system-wide data analysis in a metabolic context and for streamlined reconstruction of metabolic networks based on protein homology. The RAVEN Toolbox workflow was applied in order to reconstruct a genome-scale metabolic model for the important microbial cell factory Penicillium chrysogenum Wisconsin54-1255. The model was validated in a bibliomic study of in total 440 references, and it comprises 1471 unique biochemical reactions and 1006 ORFs. It was then used to study the roles of ATP and NADPH in the biosynthesis of penicillin, and to identify potential metabolic engineering targets for maximization of penicillin production.

Quantifying diet-induced metabolic changes of the human gut microbiome

The human gut microbiome is known to be associated with various human disorders, but a major challenge is to go beyond association studies and elucidate causalities. Mathematical modeling of the human gut microbiome at a genome scale is a useful tool to decipher microbe-microbe, diet-microbe and microbe-host interactions. Here, we describe the CASINO (Community And Systems-level INteractive Optimization) toolbox, a comprehensive computational platform for analysis of microbial communities through metabolic modeling. We first validated the toolbox by simulating and testing the performance of single bacteria and whole communities in vitro. Focusing on metabolic interactions between the diet, gut microbiota, and host metabolism, we demonstrated the predictive power of the toolbox in a diet-intervention study of 45 obese and overweight individuals and validated our predictions by fecal and blood metabolomics data. Thus, modeling could quantitatively describe altered fecal and serum amino acid levels in response to diet intervention.

Understanding the interactions between bacteria in the human gut through metabolic modeling

The human gut microbiome plays an influential role in maintaining human health, and it is a potential target for prevention and treatment of disease. Genome-scale metabolic models (GEMs) can provide an increased understanding of the mechanisms behind the effects of diet, the genotype-phenotype relationship and microbial robustness. Here we reconstructed GEMs for three key species, (Bacteroides thetaiotamicron, Eubacterium rectale and Methanobrevibacter smithii) as relevant representatives of three main phyla in the human gut (Bacteroidetes, Firmicutes and Euryarchaeota). We simulated the interactions between these three bacteria in different combinations of gut ecosystems and compared the predictions with the experimental results obtained from colonization of germ free mice. Furthermore, we used our GEMs for analyzing the contribution of each species to the overall metabolism of the gut microbiota based on transcriptome data and demonstrated that these models can be used as a scaffold for understanding bacterial interactions in the gut.

Genome-scale metabolic reconstructions of Pichia stipitis and Pichia pastoris and in silico evaluation of their potentials

BackgroundPichia stipitis and Pichia pastoris have long been investigated due to their native abilities to metabolize every sugar from lignocellulose and to modulate methanol consumption, respectively. The latter has been driving the production of several recombinant proteins. As a result, significant advances in their biochemical knowledge, as well as in genetic engineering and fermentation methods have been generated. The release of their genome sequences has allowed systems level research.ResultsIn this work, genome-scale metabolic models (GEMs) of P. stipitis (iSS884) and P. pastoris (iLC915) were reconstructed. iSS884 includes 1332 reactions, 922 metabolites, and 4 compartments. iLC915 contains 1423 reactions, 899 metabolites, and 7 compartments. Compared with the previous GEMs of P. pastoris, PpaMBEL1254 and iPP668, iLC915 contains more genes and metabolic functions, as well as improved predictive capabilities. Simulations of physiological responses for the growth of both yeasts on selected carbon sources using iSS884 and iLC915 closely reproduced the experimental data. Additionally, the iSS884 model was used to predict ethanol production from xylose at different oxygen uptake rates. Simulations with iLC915 closely reproduced the effect of oxygen uptake rate on physiological states of P. pastoris expressing a recombinant protein. The potential of P. stipitis for the conversion of xylose and glucose into ethanol using reactors in series, and of P. pastoris to produce recombinant proteins using mixtures of methanol and glycerol or sorbitol are also discussed.ConclusionsIn conclusion the first GEM of P. stipitis (iSS884) was reconstructed and validated. The expanded version of the P. pastoris GEM, iLC915, is more complete and has improved capabilities over the existing models. Both GEMs are useful frameworks to explore the versatility of these yeasts and to capitalize on their biotechnological potentials.

The gut microbiota modulates host amino acid and glutathione metabolism in mice

The gut microbiota has been proposed as an environmental factor that promotes the progression of metabolic diseases. Here, we investigated how the gut microbiota modulates the global metabolic differences in duodenum, jejunum, ileum, colon, liver, and two white adipose tissue depots obtained from conventionally raised (CONV‐R) and germ‐free (GF) mice using gene expression data and tissue‐specific genome‐scale metabolic models (GEMs). We created a generic mouse metabolic reaction (MMR) GEM, reconstructed 28 tissue‐specific GEMs based on proteomics data, and manually curated GEMs for small intestine, colon, liver, and adipose tissues. We used these functional models to determine the global metabolic differences between CONV‐R and GF mice. Based on gene expression data, we found that the gut microbiota affects the host amino acid (AA) metabolism, which leads to modifications in glutathione metabolism. To validate our predictions, we measured the level of AAs and N‐acetylated AAs in the hepatic portal vein of CONV‐R and GF mice. Finally, we simulated the metabolic differences between the small intestine of the CONV‐R and GF mice accounting for the content of the diet and relative gene expression differences. Our analyses revealed that the gut microbiota influences host amino acid and glutathione metabolism in mice.

Elucidating the interactions between the human gut microbiota and its host through metabolic modeling

Increased understanding of the interactions between the gut microbiota, diet and environmental effects may allow us to design efficient treatment strategies for addressing global health problems. Existence of symbiotic microorganisms in the human gut provides different functions for the host such as conversion of nutrients, training of the immune system, and resistance to pathogens. The gut microbiome also plays an influential role in maintaining human health, and it is a potential target for prevention and treatment of common disorders including obesity, type 2 diabetes, and atherosclerosis. Due to the extreme complexity of such disorders, it is necessary to develop mathematical models for deciphering the role of its individual elements as well as the entire system and such models may assist in better understanding of the interactions between the bacteria in the human gut and the host by use of genome-scale metabolic models (GEMs). Recently, GEMs have been employed to explore the interactions between predominant bacteria in the gut ecosystems. Additionally, these models enabled analysis of the contribution of each species to the overall metabolism of the microbiota through the integration of omics data. The outcome of these studies can be used for proposing optimal conditions for desired microbiome phenotypes. Here, we review the recent progress and challenges for elucidating the interactions between the human gut microbiota and host through metabolic modeling. We discuss how these models may provide scaffolds for analyzing high-throughput data, developing probiotics and prebiotics, evaluating the effects of probiotics and prebiotics and eventually designing clinical interventions.

Genome-scale metabolic reconstructions of Bifidobacterium adolescentis L2-32 and Faecalibacterium prausnitzii A2-165 and their interaction

BackgroundThe gut microbiota plays an important role in human health and disease by acting as a metabolic organ. Metagenomic sequencing has shown how dysbiosis in the gut microbiota is associated with human metabolic diseases such as obesity and diabetes. Modeling may assist to gain insight into the metabolic implication of an altered microbiota. Fast and accurate reconstruction of metabolic models for members of the gut microbiota, as well as methods to simulate a community of microorganisms, are therefore needed. The Integrated Microbial Genomes (IMG) database contains functional annotation for nearly 4,650 bacterial genomes. This tremendous new genomic information adds new opportunities for systems biology to reconstruct accurate genome scale metabolic models (GEMs).ResultsHere we assembled a reaction data set containing 2,340 reactions obtained from existing genome-scale metabolic models, where each reaction is assigned with KEGG Orthology. The reaction data set was then used to reconstruct two genome scale metabolic models for gut microorganisms available in the IMG database Bifidobacterium adolescentis L2-32, which produces acetate during fermentation, and Faecalibacterium prausnitzii A2-165, which consumes acetate and produces butyrate. F. prausnitzii is less abundant in patients with Crohn’s disease and has been suggested to play an anti-inflammatory role in the gut ecosystem. The B. adolescentis model, iBif452, comprises 699 reactions and 611 unique metabolites. The F. prausnitzii model, iFap484, comprises 713 reactions and 621 unique metabolites. Each model was validated with in vivo data. We used OptCom and Flux Balance Analysis to simulate how both organisms interact.ConclusionsThe consortium of iBif452 and iFap484 was applied to predict F. prausnitzii’s demand for acetate and production of butyrate which plays an essential role in colonic homeostasis and cancer prevention. The assembled reaction set is a useful tool to generate bacterial draft models from KEGG Orthology.

Identifying anti-growth factors for human cancer cell lines through genome-scale metabolic modeling

Human cancer cell lines are used as important model systems to study molecular mechanisms associated with tumor growth, hereunder how genomic and biological heterogeneity found in primary tumors affect cellular phenotypes. We reconstructed Genome scale metabolic models (GEMs) for eleven cell lines based on RNA-Seq data and validated the functionality of these models with data from metabolite profiling. We used cell line-specific GEMs to analyze the differences in the metabolism of cancer cell lines, and to explore the heterogeneous expression of the metabolic subsystems. Furthermore, we predicted 85 antimetabolites that can inhibit growth of, or even kill, any of the cell lines, while at the same time not being toxic for 83 different healthy human cell types. 60 of these antimetabolites were found to inhibit growth in all cell lines. Finally, we experimentally validated one of the predicted antimetabolites using two cell lines with different phenotypic origins, and found that it is effective in inhibiting the growth of these cell lines. Using immunohistochemistry, we also showed high or moderate expression levels of proteins targeted by the validated antimetabolite. Identified anti-growth factors for inhibition of cell growth may provide leads for the development of efficient cancer treatment strategies.

Candidalysin drives epithelial signaling, neutrophil recruitment, and immunopathology at the vaginal mucosa.

ABSTRACT Unlike other forms of candidiasis, vulvovaginal candidiasis, caused primarily by the fungal pathogen Candida albicans, is a disease of immunocompetent and otherwise healthy women. Despite its prevalence, the fungal factors responsible for initiating symptomatic infection remain poorly understood. One of the hallmarks of vaginal candidiasis is the robust recruitment of neutrophils to the site of infection, which seemingly do not clear the fungus, but rather exacerbate disease symptomatology. Candidalysin, a newly discovered peptide toxin secreted by C. albicans hyphae during invasion, drives epithelial damage, immune activation, and phagocyte attraction. Therefore, we hypothesized that Candidalysin is crucial for vulvovaginal candidiasis immunopathology. Anti-Candida immune responses are anatomical-site specific, as effective gastrointestinal, oral, and vaginal immunities are uniquely compartmentalized. Thus, we aimed to identify the immunopathologic role of Candidalysin and downstream signaling events at the vaginal mucosa. Microarray analysis of C. albicans-infected human vaginal epithelium in vitro revealed signaling pathways involved in epithelial damage responses, barrier repair, and leukocyte activation. Moreover, treatment of A431 vaginal epithelial cells with Candidalysin induced dose-dependent proinflammatory cytokine responses (including interleukin 1α [IL-1α], IL-1β, and IL-8), damage, and activation of c-Fos and mitogen-activated protein kinase (MAPK) signaling, consistent with fungal challenge. Mice intravaginally challenged with C. albicans strains deficient in Candidalysin exhibited no differences in colonization compared to isogenic controls. However, significant decreases in neutrophil recruitment, damage, and proinflammatory cytokine expression were observed with these strains. Our findings demonstrate that Candidalysin is a key hypha-associated virulence determinant responsible for the immunopathogenesis of C. albicans vaginitis.

Systems biology analysis of hepatitis C virus infection reveals the role of copy number increases in regions of chromosome 1q in hepatocellular carcinoma metabolism

Hepatitis C virus (HCV) infection is a worldwide healthcare problem; however, traditional treatment methods have failed to cure all patients, and HCV has developed resistance to new drugs. Systems biology-based analyses could play an important role in the holistic analysis of the impact of HCV on hepatocellular metabolism. Here, we integrated HCV assembly reactions with a genome-scale hepatocyte metabolic model to identify metabolic targets for HCV assembly and metabolic alterations that occur between different HCV progression states (cirrhosis, dysplastic nodule, and early and advanced hepatocellular carcinoma (HCC)) and healthy liver tissue. We found that diacylglycerolipids were essential for HCV assembly. In addition, the metabolism of keratan sulfate and chondroitin sulfate was significantly changed in the cirrhosis stage, whereas the metabolism of acyl-carnitine was significantly changed in the dysplastic nodule and early HCC stages. Our results explained the role of the upregulated expression of BCAT1, PLOD3 and six other methyltransferase genes involved in carnitine biosynthesis and S-adenosylmethionine metabolism in the early and advanced HCC stages. Moreover, GNPAT and BCAP31 expression was upregulated in the early and advanced HCC stages and could lead to increased acyl-CoA consumption. By integrating our results with copy number variation analyses, we observed that GNPAT, PPOX and five of the methyltransferase genes (ASH1L, METTL13, SMYD2, TARBP1 and SMYD3), which are all located on chromosome 1q, had increased copy numbers in the cancer samples relative to the normal samples. Finally, we confirmed our predictions with the results of metabolomics studies and proposed that inhibiting the identified targets has the potential to provide an effective treatment strategy for HCV-associated liver disorders.