RESUMO
Motivation: INCA is a powerful tool for metabolic flux analysis, however, import and export of data and results can be tedious and limit the use of INCA in automated workflows. Results: The INCAWrapper enables the use of INCA purely through Python, which allows the use of INCA in common data science workflows. Availability and implementation: The INCAWrapper is implemented in Python and can be found at https://github.com/biosustain/incawrapper. It is freely available under an MIT License. To run INCA, the user needs their own MATLAB and INCA licenses. INCA is freely available for noncommercial use at mfa.vueinnovations.com.
RESUMO
A biochemical pathway consists of a series of interconnected biochemical reactions to accomplish specific life activities. The participating reactants and resultant products of a pathway, including gene fragments, proteins, and small molecules, coalesce to form a complex reaction network. Biochemical pathways play a critical role in the biochemical domain as they can reveal the flow of biochemical reactions in living organisms, making them essential for understanding life processes. Existing studies of biochemical pathway networks are mainly based on experimentation and pathway database analysis methods, which are plagued by substantial cost constraints. Inspired by the success of representation learning approaches in biomedicine, we develop the biochemical pathway prediction (BPP) platform, which is an automatic BPP platform to identify potential links or attributes within biochemical pathway networks. Our BPP platform incorporates a variety of representation learning models, including the latest hypergraph neural networks technology to model biochemical reactions in pathways. In particular, BPP contains the latest biochemical pathway-based datasets and enables the prediction of potential participants or products of biochemical reactions in biochemical pathways. Additionally, BPP is equipped with an SHAP explainer to explain the predicted results and to calculate the contributions of each participating element. We conduct extensive experiments on our collected biochemical pathway dataset to benchmark the effectiveness of all models available on BPP. Furthermore, our detailed case studies based on the chronological pattern of our dataset demonstrate the effectiveness of our platform. Our BPP web portal, source code and datasets are freely accessible at https://github.com/Glasgow-AI4BioMed/BPP.
Assuntos
Biologia Computacional , Redes Neurais de Computação , Biologia Computacional/métodos , Redes e Vias Metabólicas , Software , Algoritmos , HumanosRESUMO
Generative modeling and representation learning of tandem mass spectrometry data aim to learn an interpretable and instrument-agnostic digital representation of metabolites directly from MS/MS spectra. Interpretable and instrument-agnostic digital representations would facilitate comparisons of MS/MS spectra between instrument vendors and enable better and more accurate queries of large MS/MS spectra databases for metabolite identification. In this study, we apply generative modeling and representation learning using variational autoencoders to understand the extent to which tandem mass spectra can be disentangled into their factors of generation (e.g., collision energy, ionization mode, instrument type, etc.) with minimal prior knowledge of the factors. We find that variational autoencoders can disentangle tandem mass spectra data with the proper choice of hyperparameters into meaningful latent representations aligned with known factors of variation. We develop a two-step approach to facilitate the selection of models that are disentangled, which could be applied to other complex and high-dimensional data sets.
Assuntos
Aprendizagem , Espectrometria de Massas em Tandem , Bases de Dados FactuaisRESUMO
Multi-omics datasets are becoming of key importance to drive discovery in fundamental research as much as generating knowledge for applied biotechnology. However, the construction of such large datasets is usually time-consuming and expensive. Automation might enable to overcome these issues by streamlining workflows from sample generation to data analysis. Here, we describe the construction of a complex workflow for the generation of high-throughput microbial multi-omics datasets. The workflow comprises a custom-built platform for automated cultivation and sampling of microbes, sample preparation protocols, analytical methods for sample analysis and automated scripts for raw data processing. We demonstrate possibilities and limitations of such workflow in generating data for three biotechnologically relevant model organisms, namely Escherichia coli, Saccharomyces cerevisiae, and Pseudomonas putida.
Assuntos
Multiômica , Fluxo de TrabalhoRESUMO
The gut microbiota impacts systemic levels of multiple metabolites including NAD+ precursors through diverse pathways. Nicotinamide riboside (NR) is an NAD+ precursor capable of regulating mammalian cellular metabolism. Some bacterial families express the NR-specific transporter, PnuC. We hypothesized that dietary NR supplementation would modify the gut microbiota across intestinal sections. We determined the effects of 12 weeks of NR supplementation on the microbiota composition of intestinal segments of high-fat diet-fed (HFD) rats. We also explored the effects of 12 weeks of NR supplementation on the gut microbiota in humans and mice. In rats, NR reduced fat mass and tended to decrease body weight. Interestingly, NR increased fat and energy absorption but only in HFD-fed rats. Moreover, 16S rRNA gene sequencing analysis of intestinal and fecal samples revealed an increased abundance of species within Erysipelotrichaceae and Ruminococcaceae families in response to NR. PnuC-positive bacterial strains within these families showed an increased growth rate when supplemented with NR. The abundance of species within the Lachnospiraceae family decreased in response to HFD irrespective of NR. Alpha and beta diversity and bacterial composition of the human fecal microbiota were unaltered by NR, but in mice, the fecal abundance of species within Lachnospiraceae increased while abundances of Parasutterella and Bacteroides dorei species decreased in response to NR. In conclusion, oral NR altered the gut microbiota in rats and mice, but not in humans. In addition, NR attenuated body fat mass gain in rats, and increased fat and energy absorption in the HFD context.
RESUMO
Automation is playing an increasingly significant role in synthetic biology. Groundbreaking technologies, developed over the past 20 years, have enormously accelerated the construction of efficient microbial cell factories. Integrating state-of-the-art tools (e.g. for genome engineering and analytical techniques) into the design-build-test-learn cycle (DBTLc) will shift the metabolic engineering paradigm from an almost artisanal labor towards a fully automated workflow. Here, we provide a perspective on how a fully automated DBTLc could be harnessed to construct the next-generation bacterial cell factories in a fast, high-throughput fashion. Innovative toolsets and approaches that pushed the boundaries in each segment of the cycle are reviewed to this end. We also present the most recent efforts on automation of the DBTLc, which heralds a fully autonomous pipeline for synthetic biology in the near future.
Assuntos
Engenharia Metabólica , Biologia Sintética , Engenharia Metabólica/métodosRESUMO
BACKGROUND: Individuals with long standing diabetes duration can experience damage to small microvascular blood vessels leading to diabetes complications (DCs) and increased mortality. Precision diagnostic tailors a diagnosis to an individual by using biomedical information. Blood small molecule profiling coupled with machine learning (ML) can facilitate the goals of precision diagnostics, including earlier diagnosis and individualized risk scoring. METHODS: Using data in a cohort of 537 adults with type 1 diabetes (T1D) we predicted five-year progression to DCs. Prediction models were computed first with clinical risk factors at baseline and then with clinical risk factors and blood-derived molecular data at baseline. Progression of diabetic kidney disease and diabetic retinopathy were predicted in two complication-specific models. FINDINGS: The model predicts the progression to diabetic kidney disease with accuracy: 0.96 ± 0.25 and 0.96 ± 0.06 area under curve, AUC, with clinical measurements and with small molecule predictors respectively and highlighted main predictors to be albuminuria, glomerular filtration rate, retinopathy status at baseline, sugar derivatives and ketones. For diabetic retinopathy, AUC 0.75 ± 0.14 and 0.79 ± 0.16 with clinical measurements and with small molecule predictors respectively and highlighted key predictors, albuminuria, glomerular filtration rate and retinopathy status at baseline. Individual risk scores were built to visualize results. INTERPRETATION: With further validation ML tools could facilitate the implementation of precision diagnosis in the clinic. It is envisaged that patients could be screened for complications, before these occur, thus preserving healthy life-years for persons with diabetes. FUNDING: This study has been financially supported by Novo Nordisk Foundation grant NNF14OC0013659.
Assuntos
Diabetes Mellitus Tipo 1 , Nefropatias Diabéticas , Retinopatia Diabética , Adulto , Albuminúria , Diabetes Mellitus Tipo 1/complicações , Diabetes Mellitus Tipo 1/diagnóstico , Nefropatias Diabéticas/diagnóstico , Nefropatias Diabéticas/etiologia , Retinopatia Diabética/diagnóstico , Retinopatia Diabética/etiologia , Taxa de Filtração Glomerular , Humanos , Fatores de RiscoRESUMO
Machine learning has greatly advanced over the past decade, owing to advances in algorithmic innovations, hardware acceleration, and benchmark datasets to train on domains such as computer vision, natural-language processing, and more recently the life sciences. In particular, the subfield of machine learning known as deep learning has found applications in genomics, proteomics, and metabolomics. However, a thorough assessment of how the data preprocessing methods required for the analysis of life science data affect the performance of deep learning is lacking. This work contributes to filling that gap by assessing the impact of commonly used as well as newly developed methods employed in data preprocessing workflows for metabolomics that span from raw data to processed data. The results from these analyses are summarized into a set of best practices that can be used by researchers as a starting point for downstream classification and reconstruction tasks using deep learning.
RESUMO
Technological advances in high-resolution mass spectrometry (MS) vastly increased the number of samples that can be processed in a life science experiment, as well as volume and complexity of the generated data. To address the bottleneck of high-throughput data processing, we present SmartPeak (https://github.com/AutoFlowResearch/SmartPeak), an application that encapsulates advanced algorithms to enable fast, accurate, and automated processing of capillary electrophoresis-, gas chromatography-, and liquid chromatography (LC)-MS(/MS) data and high-pressure LC data for targeted and semitargeted metabolomics, lipidomics, and fluxomics experiments. The application allows for an approximate 100-fold reduction in the data processing time compared to manual processing while enhancing quality and reproducibility of the results.
Assuntos
Processamento Eletrônico de Dados/métodos , Metabolômica/métodos , Automação , Cromatografia Líquida , Eletroforese Capilar , Espectrometria de Massas em Tandem , Fatores de TempoRESUMO
Small-molecule binding allosteric transcription factors (aTFs) derived from bacteria enable real-time monitoring of metabolite abundances, high-throughput screening of genetic designs, and dynamic control of metabolism. Yet, engineering of reporter promoter designs of prokaryotic aTF biosensors in eukaryotic cells is complex. Here we investigate the impact of aTF binding site positions at single-nucleotide resolution in >300 reporter promoter designs in Saccharomyces cerevisiae. From this we identify biosensor output landscapes with transient and distinct aTF binding site position effects for aTF repressors and activators, respectively. Next, we present positions for tunable reporter promoter outputs enabling metabolite-responsive designs for a total of four repressor-type and three activator-type aTF biosensors with dynamic output ranges up to 8- and 26-fold, respectively. This study highlights aTF binding site positions in reporter promoters as key for successful biosensor engineering and that repressor-type aTF biosensors allows for more flexibility in terms of choice of binding site positioning compared to activator-type aTF biosensors.
Assuntos
Técnicas Biossensoriais/métodos , Genes Reporter/genética , Saccharomyces cerevisiae/metabolismo , Sítios de Ligação , Plasmídeos/genética , Plasmídeos/metabolismo , Regiões Promotoras Genéticas , Engenharia de Proteínas , Saccharomyces cerevisiae/genética , Fatores de Transcrição/genética , Fatores de Transcrição/metabolismoRESUMO
[This corrects the article DOI: 10.3389/fgene.2019.00747.].
RESUMO
Fatty alcohols are widely used in various applications within a diverse set of industries, such as the soap and detergent industry, the personal care, and cosmetics industry, as well as the food industry. The total world production of fatty alcohols is over 2 million tons with approximately equal parts derived from fossil oil and from plant oils or animal fats. Due to the environmental impact of these production methods, there is an interest in alternative methods for fatty alcohol production via microbial fermentation using cheap renewable feedstocks. In this study, we aimed to obtain a better understanding of how fatty alcohol biosynthesis impacts the host organism, baker's yeast Saccharomyces cerevisiae or oleaginous yeast Yarrowia lipolytica. Producing and non-producing strains were compared in growth and nitrogen-depletion cultivation phases. The multi-omics analysis included physiological characterization, transcriptome analysis by RNAseq, 13Cmetabolic flux analysis, and intracellular metabolomics. Both species accumulated fatty alcohols under nitrogen-depletion conditions but not during growth. The fatty alcohol-producing Y. lipolytica strain had a higher fatty alcohol production rate than an analogous S. cerevisiae strain. Nitrogen-depletion phase was associated with lower glucose uptake rates and a decrease in the intracellular concentration of acetyl-CoA in both yeast species, as well as increased organic acid secretion rates in Y. lipolytica. Expression of the fatty alcohol-producing enzyme fatty acyl-CoA reductase alleviated the growth defect caused by deletion of hexadecenal dehydrogenase encoding genes (HFD1 and HFD4) in Y. lipolytica. RNAseq analysis showed that fatty alcohol production triggered a cell wall stress response in S. cerevisiae. RNAseq analysis also showed that both nitrogen-depletion and fatty alcohol production have substantial effects on the expression of transporter encoding genes in Y. lipolytica. In conclusion, through this multi-omics study, we uncovered some effects of fatty alcohol production on the host metabolism. This knowledge can be used as guidance for further strain improvement towards the production of fatty alcohols.
RESUMO
Growth rate and yield are fundamental features of microbial growth. However, we lack a mechanistic and quantitative understanding of the rate-yield relationship. Studies pairing computational predictions with experiments have shown the importance of maintenance energy and proteome allocation in explaining rate-yield tradeoffs and overflow metabolism. Recently, adaptive evolution experiments of Escherichia coli reveal a phenotypic diversity beyond what has been explained using simple models of growth rate versus yield. Here, we identify a two-dimensional rate-yield tradeoff in adapted E. coli strains where the dimensions are (A) a tradeoff between growth rate and yield and (B) a tradeoff between substrate (glucose) uptake rate and growth yield. We employ a multi-scale modeling approach, combining a previously reported coarse-grained small-scale proteome allocation model with a fine-grained genome-scale model of metabolism and gene expression (ME-model), to develop a quantitative description of the full rate-yield relationship for E. coli K-12 MG1655. The multi-scale analysis resolves the complexity of ME-model which hindered its practical use in proteome complexity analysis, and provides a mechanistic explanation of the two-dimensional tradeoff. Further, the analysis identifies modifications to the P/O ratio and the flux allocation between glycolysis and pentose phosphate pathway (PPP) as potential mechanisms that enable the tradeoff between glucose uptake rate and growth yield. Thus, the rate-yield tradeoffs that govern microbial adaptation to new environments are more complex than previously reported, and they can be understood in mechanistic detail using a multi-scale modeling approach.
Assuntos
Proteínas de Bactérias/metabolismo , Escherichia coli/metabolismo , Evolução Molecular , Proteínas de Bactérias/genética , Escherichia coli/genética , Genoma Bacteriano/genética , Modelos Biológicos , Proteoma/genética , Proteoma/metabolismo , Biologia de SistemasRESUMO
Current machine learning techniques enable robust association of biological signals with measured phenotypes, but these approaches are incapable of identifying causal relationships. Here, we develop an integrated "white-box" biochemical screening, network modeling, and machine learning approach for revealing causal mechanisms and apply this approach to understanding antibiotic efficacy. We counter-screen diverse metabolites against bactericidal antibiotics in Escherichia coli and simulate their corresponding metabolic states using a genome-scale metabolic network model. Regression of the measured screening data on model simulations reveals that purine biosynthesis participates in antibiotic lethality, which we validate experimentally. We show that antibiotic-induced adenine limitation increases ATP demand, which elevates central carbon metabolism activity and oxygen consumption, enhancing the killing effects of antibiotics. This work demonstrates how prospective network modeling can couple with machine learning to identify complex causal mechanisms underlying drug efficacy.
Assuntos
Antibacterianos/metabolismo , Antibacterianos/farmacologia , Redes e Vias Metabólicas/efeitos dos fármacos , Adenina/metabolismo , Biologia Computacional/métodos , Avaliação Pré-Clínica de Medicamentos/métodos , Escherichia coli/metabolismo , Aprendizado de Máquina , Redes e Vias Metabólicas/imunologia , Modelos Teóricos , Purinas/metabolismoRESUMO
Biological regulatory network architectures are multi-scale in their function and can adaptively acquire new functions. Gene knockout (KO) experiments provide an established experimental approach not just for studying gene function, but also for unraveling regulatory networks in which a gene and its gene product are involved. Here we study the regulatory architecture of Escherichia coli K-12 MG1655 by applying adaptive laboratory evolution (ALE) to metabolic gene KO strains. Multi-omic analysis reveal a common overall schema describing the process of adaptation whereby perturbations in metabolite concentrations lead regulatory networks to produce suboptimal states, whose function is subsequently altered and re-optimized through acquisition of mutations during ALE. These results indicate that metabolite levels, through metabolite-transcription factor interactions, have a dominant role in determining the function of a multi-scale regulatory architecture that has been molded by evolution.
Assuntos
Escherichia coli K12/fisiologia , Evolução Molecular , Redes Reguladoras de Genes/fisiologia , Redes e Vias Metabólicas/genética , Técnicas de Inativação de Genes , Mutação , FenótipoRESUMO
Adaptive laboratory evolution (ALE) has emerged as a new approach with which to pursue fundamental biological inquiries and, in particular, new insights into the systemic function of a gene product. Two E. coli knockout strains were constructed: one that blocked the Pentose Phosphate Pathway (gnd KO) and one that decoupled the TCA cycle from electron transport (sdhCDAB KO). Despite major perturbations in central metabolism, minimal growth rate changes were found in the two knockout strains. More surprisingly, many similarities were found in their initial transcriptomic states that could be traced to similarly perturbed metabolites despite the differences in the network location of the gene perturbations and concomitant re-routing of pathway fluxes around these perturbations. However, following ALE, distinct metabolomic and transcriptomic states were realized. These included divergent flux and gene expression profiles in the gnd and sdhCDAB KOs to overcome imbalances in NADPH production and nitrogen/sulfur assimilation, respectively, that were not obvious limitations of growth in the unevolved knockouts. Therefore, this work demonstrates that ALE provides a productive approach to reveal novel insights of gene function at a systems level that cannot be found by observing the fresh knockout alone.
RESUMO
A mechanistic understanding of how new phenotypes develop to overcome the loss of a gene product provides valuable insight on both the metabolic and regulatory functions of the lost gene. The pgi gene, whose product catalyzes the second step in glycolysis, was deleted in a growth-optimized Escherichia coli K-12 MG1655 strain. The initial knockout (KO) strain exhibited an 80% drop in growth rate that was largely recovered in eight replicate, but phenotypically distinct, cultures after undergoing adaptive laboratory evolution (ALE). Multi-omic data sets showed that the loss of pgi substantially shifted pathway usage, leading to a redox and sugar phosphate stress response. These stress responses were overcome by unique combinations of innovative mutations selected for by ALE. Thus, the coordinated mechanisms from genome to metabolome that lead to multiple optimal phenotypes after the loss of a major gene product were revealed.IMPORTANCE A mechanistic understanding of how microbes are able to overcome the loss of a gene through regulatory and metabolic changes is not well understood. Eight independent adaptive laboratory evolution (ALE) experiments with pgi knockout strains resulted in eight phenotypically distinct endpoints that were able to overcome the gene loss. Utilizing multi-omics analysis, the coordinated mechanisms from genome to metabolome that lead to multiple optimal phenotypes after the loss of a major gene product were revealed.
Assuntos
Escherichia coli K12/enzimologia , Escherichia coli K12/genética , Proteínas de Escherichia coli/genética , Glucose-6-Fosfato Isomerase/genética , Escherichia coli K12/metabolismo , Proteínas de Escherichia coli/metabolismo , Técnicas de Inativação de Genes , Glucose-6-Fosfato Isomerase/metabolismo , Glicólise , Mutação , Oxirredução , FenótipoRESUMO
Aromatic metabolites provide the backbone for numerous industrial and pharmaceutical compounds of high value. The Phosphotransferase System (PTS) is common to many bacteria, and is the primary mechanism for glucose uptake by Escherichia coli. The PTS was removed to conserve phosphoenolpyruvate (pep), which is a precursor for aromatic metabolites and consumed by the PTS, for aromatic metabolite production. Replicate adaptive laboratory evolution (ALE) of PTS and detailed omics data sets collected revealed that the PTS bridged the gap between respiration and fermentation, leading to distinct high fermentative and high respiratory rate phenotypes. It was also found that while all strains retained high levels of aromatic amino acid (AAA) biosynthetic precursors, only one replicate from the high glycolytic clade retained high levels of intracellular AAAs. The fast growth and high AAA precursor phenotypes could provide a starting host for cell factories targeting the overproduction aromatic metabolites.
Assuntos
Aminoácidos Aromáticos , Evolução Molecular Direcionada , Metabolismo Energético , Escherichia coli , Consumo de Oxigênio , Sistema Fosfotransferase de Açúcar do Fosfoenolpiruvato/genética , Aminoácidos Aromáticos/biossíntese , Aminoácidos Aromáticos/genética , Escherichia coli/genética , Escherichia coli/metabolismoRESUMO
Methylglyoxal is a highly toxic metabolite that can be produced in all living organisms. Methylglyoxal was artificially elevated by removal of the tpiA gene from a growth optimized Escherichia coli strain. The initial response to elevated methylglyoxal and its toxicity was characterized, and detoxification mechanisms were studied using adaptive laboratory evolution. We found that: 1) Multi-omics analysis revealed biological consequences of methylglyoxal toxicity, which included attack on macromolecules including DNA and RNA and perturbation of nucleotide levels; 2) Counter-intuitive cross-talk between carbon starvation and inorganic phosphate signalling was revealed in the tpiA deletion strain that required mutations in inorganic phosphate signalling mechanisms to alleviate; and 3) The split flux through lower glycolysis depleted glycolytic intermediates requiring a host of synchronized and coordinated mutations in non-intuitive network locations in order to re-adjust the metabolic flux map to achieve optimal growth. Such mutations included a systematic inactivation of the Phosphotransferase System (PTS) and alterations in cell wall biosynthesis enzyme activity. This study demonstrated that deletion of major metabolic genes followed by ALE was a productive approach to gain novel insight into the systems biology underlying optimal phenotypic states.