基因组规模代谢网络模型构建及其应用

微生物制造产业的发展迫切需要进一步提高认识、设计和改造微生物细胞代谢的能力,以推动工业生物技术快速发展。随着微生物全基因组序列等高通量数据的不断积聚和生物信息学策略的持续涌现,使全局性、系统化地解析、设计、调控微生物生理代谢功能成为可能。而基于基因组序列注释和详细生化信息整合的基因组规模代谢网络模型(GSMM)构建为全局理解和理性调控微生物生理代谢功能提供了最佳平台。以下在详述GSMM的应用基础上,描述了如何构建一个高精确度的GSMM,并展望了未来的发展方向。     

面向工业生物技术的系统生物学

工业生物技术是以微生物细胞工厂利用可再生的生物原料来生产能源、材料与化学品等的生物技术,在解决资源、能源与环境等问题方面起着越来越重要的作用。系统生物学是全面解析微生物细胞工厂及其发酵过程从"黑箱"到"白箱"的重要研究方法。系统生物学借助基因组、转录组、蛋白质组、代谢组以及代谢流组等多组学数据,可解析微生物细胞工厂在RNA、蛋白与代谢物等不同水平上的变化规律与调控机制。目前,系统生物学在微生物细胞工厂的设计创建与发酵工艺优化中起着越来越重要的指导作用,许多成功应用实例不断涌现,推动着工业生物技术的快速发展。文中重点综述基因组、转录组、蛋白质组、代谢组与代谢流组以及基因组规模的网络模型等各组学技术的最新发展及其在工业生物技术尤其是菌株改造与发酵优化中的应用,并就工业生物技术中系统生物学的未来发展方向进行展望。     

工业微生物：创新与突破专刊序言(2021)

工业微生物及其产品广泛用于工业、农业、医药等诸多领域,相关产业在国民经济中具有举足轻重的地位。高效的菌株是提高生产效率的核心,而先进发酵技术和仪器平台对充分开发菌株代谢潜能也很重要。近年来,工业微生物领域的研究取得了快速进展,人工智能、高效基因组编辑技术和合成生物学技术逐渐广泛使用,相关产业应用也在不断扩展。为进一步促进工业微生物在生物制造等领域的应用,《生物工程学报》特组织出版专刊,从微生物菌株的多样性和生理代谢、菌株改造技术、发酵过程优化和放大,高通量微液滴培养装备开发以及工业微生物应用等方面,分别阐述目前的研究进展,并展望未来的发展趋势,为促进工业微生物及生物制造等产业的发展奠定基础。     

乳酸乳球菌NZ9000基因组规模代谢网络模型的构建与验证

随着后基因组时代的到来,工业微生物的代谢工程改造在工业生产上发挥着越来越重要的作用。而基因组规模代谢网络模型(Genome-scalemetabolicmodel,GSMM)将生物体体内所有已知代谢信息进行整合,为全局理解生物体的代谢状态、理性指导代谢工程改造提供了最佳的平台。乳酸乳球菌NZ9000(Lactococcuslactis NZ9000)作为工业发酵领域的重要菌株之一,由于其遗传背景清晰且几乎不分泌蛋白,是基因工程改造和外源蛋白表达的理想模式菌株。文中基于基因组功能注释和比较基因组学构建了L.lactisNZ9000的首个基因组规模代谢网络模型iWK557,包含557个基因、668个代谢物、840个反应,并进一步在定性和定量两个层次验证了iWK557的准确性,以期为理性指导L. lactis NZ9000代谢工程改造提供良好工具。     

微生物基因组尺度代谢网络模型构建方法的比较分析

随着后基因时代的到来,微生物的定向菌种改造在生产中发挥着越来越重要的作用.基因组尺度代谢网络模型是微生物定向改造中一种不可缺少的指导性工具,可降低菌种改造的盲目性,增加目的性和成功率.随着研究的深入,基因组尺度代谢网络模型的构建方法也越来越多.究竟选择什么样的方法才能构建出全面准确的基因组尺度代谢网络模型,对于初学者来说是一个大难题.论文结合本课题组的研究,将近年文献报导中出现过的模型构建方法进行了分类和分析,并评述了各种方法的优、缺点,以期为初学者提供参考.具体介绍的方法有:基于基因组注释构建代谢网络模型,基于蛋白组构建代谢网络模型,基于文献挖掘构建代谢网络模型,通过软件和网络平台构建代谢网络模型,基于京都基因与基因组百科全书(KEGG)构建代谢网络模型.五种构建代谢网络模型方法都有其优点,但也有不可避免的缺点,要构建较为准确全面的基因组尺度代谢网络模型,需要将各种方法结合,弥补彼此的不足.图4表0参48     

基于Matlab对解脂耶氏酵母基因组规模代谢网络模型仿真结果的可视化分析

【目的】近十年来,基因组代谢网络模型迅速发展。通过构建基因组代谢网络模型进行计算机仿真模拟已成为研究生物体复杂的生理代谢不可或缺的工具。实现对仿真结果的可视化分析,可以直观地追踪模型中的代谢流向,从而更好地对仿真结果进行分析。【方法】在简要概述目前可视化方法的基础上,提出了一种基于Matlab实现基因组规模代谢网络模型仿真结果可视化的方法：通过CellDesigner预先绘制与模型相匹配的图,通过RAVEN toolbox中的函数于Matlab进行读图、并实现仿真结果的可视化。【结果】以解脂耶氏酵母基因组规模代谢网络模型iYL619_PCP v1.7为对象,实现并阐明其仿真结果的可视化。【结论】通过该方法可以清晰地监测模型中的流量和流向变化,提高仿真结果的分析效率。     

应用代谢网络模型解析工业微生物胞内代谢

为了快速、高效地理解工业微生物胞内代谢特征,寻找潜在的代谢工程改造靶点,基因组规模代谢网络模型(GSMM)作为一种系统生物学工具越来越受到人们的关注。文中在回顾GSMM 20年发展历程的基础上,分析了当前GSMM的研究现状,总结了GSMM的构建及分析方法,从预测细胞表型和指导代谢工程两个方面阐述了GSMM在解析工业微生物胞内代谢中的应用,并展望了GSMM未来的发展趋势。     

工业微生物代谢网络模型的研究进展及应用

基因组规模代谢网络(Genome-scale metabolic network model,GSMM)是工业微生物菌株定向改造过程中一种极为重要的指导性工具,有助于研究者快速获取特定性状的工业微生物,因此越来越受到人们的关注。文中回顾了GSMM的发展历程,总结并评述了GSMM的构建方法,以4种重要工业微生物(枯草芽孢杆菌Bacillus subtilis、大肠杆菌Escherichia coli、谷氨酸棒杆菌Corynebacterium glutamicum和酿酒酵母Saccharomyces cerevisiae)为例,阐述了GSMM在工业微生物中的发展与应用。此外,还对GSMM未来的发展趋势进行了展望。     

基元模式分析在生物网络和途径分析中的应用

基元模式分析是应用最广泛的代谢途径分析方法。基元模式分析的研究对象从代谢网络发展到信号传导网络;研究尺度从细胞到生物反应器,甚至生态系统;数学描述从稳态分解到动态解析;研究领域从微生物代谢到人类疾病。以下综述了基元模式分析的算法和软件开发现状,以及其在代谢途径与鲁棒性、代谢通量分解、稳态代谢通量分析、动态模型与生物过程模拟、网络结构与调控、菌株设计和信号传导网络等方面的应用。开发新的算法解决组合爆炸问题,探索基元模式与代谢调控的关系以及提高菌株设计算法效率是今后基元模式的重要发展方向。     

海洋氮循环过程及基于基因组代谢网络模型的预测

海洋氮循环在地球元素循环中充当着必不可少的角色。海洋氮循环是由一系列氧化还原反应构成的生物化学过程。固氮作用和氮同化作用为生态系统提供了生物可用氮(铵盐)。硝化作用可进一步将铵盐氧化为硝酸盐,硝酸盐又可以通过反硝化作用转化为氮气。整个氮循环实现了海洋中不同含氮无机盐间的转换。微生物是海洋氮循环的重要驱动者,海洋氮循环的研究可以帮助理解海洋生物与地球环境相互作用及协同演化的机制,从而更好地保护地球生态环境。随着氮循环关键微生物基因组尺度代谢网络模型的发表,研究者可以利用代谢网络模型来研究不同氮循环过程的效率、环境因子对氮循环过程的影响以及解析氮循环及生物网络的内在机理等,从而帮助人们更深入地研究海洋氮转化机制。本文主要综述了海洋氮循环过程中各个转化过程的主要微生物,以及基因组尺度代谢网络模型在分析氮循环中的应用。     

Use of genome-scale metabolic models for understanding microbial physiology

The exploitation of microorganisms in industrial, medical, food and environmental biotechnology requires a comprehensive understanding of their physiology. The availability of genome sequences and accumulation of high-throughput data allows gaining understanding of microbial physiology at the systems level, and genome-scale metabolic models represent a valuable framework for integrative analysis of metabolism of microorganisms. Genome-scale metabolic models are reconstructed based on a combination of genome sequence information and detailed biochemical information, and these reconstructed models can be used for analyzing and simulating the operation of metabolism in response to different stimuli. Here we discuss the requirement for having detailed physiological insight in order to exploit microorganisms for production of fuels, chemicals and pharmaceuticals. We further describe the reconstruction process of genome-scale metabolic models and different algorithms that can be used to apply these models to gain improved insight into microbial physiology.     

Constraints-based genome-scale metabolic simulation for systems metabolic engineering

Random mutagenesis and selection approaches used traditionally for the development of industrial strains have largely been complemented by metabolic engineering, which allows purposeful modification of metabolic and cellular characteristics by using recombinant DNA and other molecular biological techniques. As systems biology advances as a new paradigm of research thanks to the development of genome-scale computational tools and high-throughput experimental technologies including omics, systems metabolic engineering allowing modification of metabolic, regulatory and signaling networks of the cell at the systems-level is becoming possible. In silico genome-scale metabolic model and its simulation play increasingly important role in providing systematic strategies for metabolic engineering. The in silico genome-scale metabolic model is developed using genomic annotation, metabolic reactions, literature information, and experimental data. The advent of in silico genome-scale metabolic model brought about the development of various algorithms to simulate the metabolic status of the cell as a whole. In this paper, we review the algorithms developed for the system-wide simulation and perturbation of cellular metabolism, discuss the characteristics of these algorithms, and suggest future research direction.     

Genome-scale metabolic models as platforms for strain design and biological discovery

Genome-scale metabolic models (GEMs) have been developed and used in guiding systems’ metabolic engineering strategies for strain design and development. This strategy has been used in fermentative production of bio-based industrial chemicals and fuels from alternative carbon sources. However, computer-aided hypotheses building using established algorithms and software platforms for biological discovery can be integrated into the pipeline for strain design strategy to create superior strains of microorganisms for targeted biosynthetic goals. Here, I described an integrated workflow strategy using GEMs for strain design and biological discovery. Specific case studies of strain design and biological discovery using Escherichia coli genome-scale model are presented and discussed. The integrated workflow presented herein, when applied carefully would help guide future design strategies for high-performance microbial strains that have existing and forthcoming genome-scale metabolic models.     

Evolutionary programming as a platform for <Emphasis Type="Italic">in silico</Emphasis> metabolic engineering

Background  

Through genetic engineering it is possible to introduce targeted genetic changes and hereby engineer the metabolism of microbial cells with the objective to obtain desirable phenotypes. However, owing to the complexity of metabolic networks, both in terms of structure and regulation, it is often difficult to predict the effects of genetic modifications on the resulting phenotype. Recently genome-scale metabolic models have been compiled for several different microorganisms where structural and stoichiometric complexity is inherently accounted for. New algorithms are being developed by using genome-scale metabolic models that enable identification of gene knockout strategies for obtaining improved phenotypes. However, the problem of finding optimal gene deletion strategy is combinatorial and consequently the computational time increases exponentially with the size of the problem, and it is therefore interesting to develop new faster algorithms.     

Semi-automated Curation of Metabolic Models via Flux Balance Analysis: A Case Study with Mycoplasma gallisepticum

Primarily used for metabolic engineering and synthetic biology, genome-scale metabolic modeling shows tremendous potential as a tool for fundamental research and curation of metabolism. Through a novel integration of flux balance analysis and genetic algorithms, a strategy to curate metabolic networks and facilitate identification of metabolic pathways that may not be directly inferable solely from genome annotation was developed. Specifically, metabolites involved in unknown reactions can be determined, and potentially erroneous pathways can be identified. The procedure developed allows for new fundamental insight into metabolism, as well as acting as a semi-automated curation methodology for genome-scale metabolic modeling. To validate the methodology, a genome-scale metabolic model for the bacterium Mycoplasma gallisepticum was created. Several reactions not predicted by the genome annotation were postulated and validated via the literature. The model predicted an average growth rate of 0.358±0.12, closely matching the experimentally determined growth rate of M. gallisepticum of 0.244±0.03. This work presents a powerful algorithm for facilitating the identification and curation of previously known and new metabolic pathways, as well as presenting the first genome-scale reconstruction of M. gallisepticum.     

Genome-scale metabolic model in guiding metabolic engineering of microbial improvement

In the past few decades, despite all the significant achievements in industrial microbial improvement, the approaches of traditional random mutation and selection as well as the rational metabolic engineering based on the local knowledge cannot meet today’s needs. With rapid reconstructions and accurate in silico simulations, genome-scale metabolic model (GSMM) has become an indispensable tool to study the microbial metabolism and design strain improvements. In this review, we highlight the application of GSMM in guiding microbial improvements focusing on a systematic strategy and its achievements in different industrial fields. This strategy includes a repetitive process with four steps: essential data acquisition, GSMM reconstruction, constraints-based optimizing simulation, and experimental validation, in which the second and third steps are the centerpiece. The achievements presented here belong to different industrial application fields, including food and nutrients, biopharmaceuticals, biopolymers, microbial biofuel, and bioremediation. This strategy and its achievements demonstrate a momentous guidance of GSMM for metabolic engineering breeding of industrial microbes. More efforts are required to extend this kind of study in the meantime.     

Genome-scale metabolic representation of Amycolatopsis balhimycina

Infection caused by methicillin-resistant Staphylococcus aureus (MRSA) is an increasing societal problem. Typically, glycopeptide antibiotics are used in the treatment of these infections. The most comprehensively studied glycopeptide antibiotic biosynthetic pathway is that of balhimycin biosynthesis in Amycolatopsis balhimycina. The balhimycin yield obtained by A. balhimycina is, however, low and there is therefore a need to improve balhimycin production. In this study, we performed genome sequencing, assembly and annotation analysis of A. balhimycina and further used these annotated data to reconstruct a genome-scale metabolic model for the organism. Here we generated an almost complete A. balhimycina genome sequence comprising 10,562,587 base pairs assembled into 2,153 contigs. The high GC-genome (～ 69%) includes 8,585 open reading frames (ORFs). We used our integrative toolbox called SEQTOR for functional annotation and then integrated annotated data with biochemical and physiological information available for this organism to reconstruct a genome-scale metabolic model of A. balhimycina. The resulting metabolic model contains 583 ORFs as protein encoding genes (7% of the predicted 8,585 ORFs), 407 EC numbers, 647 metabolites and 1,363 metabolic reactions. During the analysis of the metabolic model, linear, quadratic and evolutionary programming algorithms using flux balance analysis (FBA), minimization of metabolic adjustment (MOMA), and OptGene, respectively were applied as well as phenotypic behavior and improved balhimycin production were simulated. The A. balhimycina model shows a good agreement between in silico data and experimental data and also identifies key reactions associated with increased balhimycin production. The reconstruction of the genome-scale metabolic model of A. balhimycina serves as a basis for physiological characterization. The model allows a rational design of engineering strategies for increasing balhimycin production in A. balhimycina and glycopeptide production in general.     

Large-scale bi-level strain design approaches and mixed-integer programming solution techniques

The use of computational models in metabolic engineering has been increasing as more genome-scale metabolic models and computational approaches become available. Various computational approaches have been developed to predict how genetic perturbations affect metabolic behavior at a systems level, and have been successfully used to engineer microbial strains with improved primary or secondary metabolite production. However, identification of metabolic engineering strategies involving a large number of perturbations is currently limited by computational resources due to the size of genome-scale models and the combinatorial nature of the problem. In this study, we present (i) two new bi-level strain design approaches using mixed-integer programming (MIP), and (ii) general solution techniques that improve the performance of MIP-based bi-level approaches. The first approach (SimOptStrain) simultaneously considers gene deletion and non-native reaction addition, while the second approach (BiMOMA) uses minimization of metabolic adjustment to predict knockout behavior in a MIP-based bi-level problem for the first time. Our general MIP solution techniques significantly reduced the CPU times needed to find optimal strategies when applied to an existing strain design approach (OptORF) (e.g., from ～10 days to ～5 minutes for metabolic engineering strategies with 4 gene deletions), and identified strategies for producing compounds where previous studies could not (e.g., malate and serine). Additionally, we found novel strategies using SimOptStrain with higher predicted production levels (for succinate and glycerol) than could have been found using an existing approach that considers network additions and deletions in sequential steps rather than simultaneously. Finally, using BiMOMA we found novel strategies involving large numbers of modifications (for pyruvate and glutamate), which sequential search and genetic algorithms were unable to find. The approaches and solution techniques developed here will facilitate the strain design process and extend the scope of its application to metabolic engineering.     

Genome-scale analysis of Mannheimia succiniciproducens metabolism

Mannheimia succiniciproducens MBEL55E isolated from bovine rumen is a capnophilic gram-negative bacterium that efficiently produces succinic acid, an industrially important four carbon dicarboxylic acid. In order to design a metabolically engineered strain which is capable of producing succinic acid with high yield and productivity, it is essential to optimize the whole metabolism at the systems level. Consequently, in silico modeling and simulation of the genome-scale metabolic network was employed for genome-scale analysis and efficient design of metabolic engineering experiments. The genome-scale metabolic network of M. succiniciproducens consisting of 686 reactions and 519 metabolites was constructed based on reannotation and validation experiments. With the reconstructed model, the network structure and key metabolic characteristics allowing highly efficient production of succinic acid were deciphered; these include strong PEP carboxylation, branched TCA cycle, relative weak pyruvate formation, the lack of glyoxylate shunt, and non-PTS for glucose uptake. Constraints-based flux analyses were then carried out under various environmental and genetic conditions to validate the genome-scale metabolic model and to decipher the altered metabolic characteristics. Predictions based on constraints-based flux analysis were mostly in excellent agreement with the experimental data. In silico knockout studies allowed prediction of new metabolic engineering strategies for the enhanced production of succinic acid. This genome-scale in silico model can serve as a platform for the systematic prediction of physiological responses of M. succiniciproducens to various environmental and genetic perturbations and consequently for designing rational strategies for strain improvement.