Assessment of the metabolic capabilities of Haemophilus influenzae Rd through a genome-scale pathway analysis

The annotated full DNA sequence is becoming available for a growing number of organisms. This information along with additional biochemical and strain-specific data can be used to define metabolic genotypes and reconstruct cellular metabolic networks. The first free-living organism for which the entire genomic sequence was established was Haemophilus influenzae. Its metabolic network is reconstructed herein and contains 461 reactions operating on 367 intracellular and 84 extracellular metabolites. With the metabolic reaction network established, it becomes necessary to determine its underlying pathway structure as defined by the set of extreme pathways. The H. influenzae metabolic network was subdivided into six subsystems and the extreme pathways determined for each subsystem based on stoichiometric, thermodynamic, and systems-specific constraints. Positive linear combinations of these pathways can be taken to determine the extreme pathways for the complete system. Since these pathways span the capabilities of the full system, they could be used to address a number of important physiological questions. First, they were used to reconcile and curate the sequence annotation by identifying reactions whose function was not supported in any of the extreme pathways. Second, they were used to predict gene products that should be co-regulated and perhaps co-expressed. Third, they were used to determine the composition of the minimal substrate requirements needed to support the production of 51 required metabolic products such as amino acids, nucleotides, phospholipids, etc. Fourth, sets of critical gene deletions from core metabolism were determined in the presence of the minimal substrate conditions and in more complete conditions reflecting the environmental niche of H. influenzae in the human host. In the former case, 11 genes were determined to be critical while six remained critical under the latter conditions. This study represents an important milestone in theoretical biology, namely the establishment of the first extreme pathway structure of a whole genome.     

基因组尺度代谢网络研究进展

基因组尺度代谢网络从基因组序列出发,结合基因、蛋白质、代谢数据库和实验数据,从系统的角度定量研究生命体的代谢过程,了解各个组分之间的相互作用关系。这类网络模型对于生命活动理论研究和优良工程菌的构建都具有重要的理论和实践意义。以下结合作者的实际研究经验,对基因组尺度代谢网络从重构到模拟直至应用进行了较为详细的介绍,并讨论了一些目前存在的难题和未来的研究方向。     

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

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

Phylogenetic analysis based on genome-scale metabolic pathway reaction content

Phylogenetic classifications based on single genes such as rRNA genes do not provide a complete and accurate picture of evolution because they do not account for evolutionary leaps caused by gene transfer, duplication, deletion and functional replacement. Here, we present a whole-genome-scale phylogeny based on metabolic pathway reaction content. From the genome sequences of 42 microorganisms, we deduced the metabolic pathway reactions and used the relatedness of these contents to construct a phylogenetic tree that represents the similarity of metabolic profiles (relatedness) as well as the extent of metabolic pathway similarity (evolutionary distance). This method accounts for horizontal gene transfer and specific gene loss by comparison of whole metabolic subpathways, and allows evaluation of evolutionary relatedness and changes in metabolic pathways. Thus, a tree based on metabolic pathway content represents both the evolutionary time scale (changes in genetic content) and the evolutionary process (changes in metabolism).     

基于约束的基因组规模代谢网络模型构建方法研究进展

基因组规模代谢网络模型(Genome-scale metabolic network model,GSMM)正成为细胞代谢特性研究的重要工具,经过多年发展相关理论方法取得了诸多进展.近年来,在基础GSMM模型基础上,通过整合基因组、转录组、蛋白组和热力学数据,实现基于各种约束的GSMM构建,在基因靶点识别、系统代谢工程...     

基因组尺度代谢网络自动重构及分析工具研究进展

高通量数据的产出为基因组尺度代谢网络的构建提供了基础,但同时也对网络构建和分析方法的改进提出了挑战。随着数据量的不断增大,耗时耗力的人工构建及分析已经无法满足模型发展的需要,因而各种自动化的方法应运而生。模型构建和分析的自动化不仅能够大幅度提高模型构建和解析的速度,同时对于模型构建和分析方法的标准化和程序化也有着不可替代的作用。文中结合作者的实际研究经验,对基因组尺度代谢网络构建的自动化进程和主要的代谢网络分析工具进行了较为详细的介绍,总结了代谢网络自动重构的流程,并提出了目前面对的主要问题和未来的研究方向。     

Porin OmpP2 of Haemophilus influenzae shows specificity for nicotinamide-derived nucleotide substrates

Haemophilus influenzae has an absolute requirement for NAD (factor V) because it lacks all biosynthetic enzymes necessary for de novo synthesis of that cofactor. Therefore, growth in vitro requires the presence of NAD itself, NMN, or nicotinamide riboside (NR). To address uptake abilities of these compounds, we investigated outer membrane proteins. By analyzing ompP2 knockout mutants, we found that NAD and NMN uptake was prevented, whereas NR uptake was not. Through investigation of the properties of purified OmpP2 in artificial lipid membrane systems, the substrate specificity of OmpP2 for NAD and NMN was determined, with KS values of approximately 8 and 4mm, respectively, in 0.1 m KCl, whereas no interaction was detected for the nucleoside NR and other purine or pyrimidine nucleotide or nucleoside species. Based on our analysis, we assume that an intrinsic binding site within OmpP2 exists that facilitates diffusion of these compounds across the outer membrane, recognizing carbonyl and exposed phosphate groups. Because OmpP2 was formerly described as a general diffusion porin, an additional property of acting as a facilitator for nicotinamide-based nucleotide transport may have evolved to support and optimize utilization of the essential cofactor sources NAD and NMN in H. influenzae.     

Detecting and investigating substrate cycles in a genome-scale human metabolic network

Substrate cycles, also known as futile cycles, are cyclic metabolic routes that dissipate energy by hydrolysing cofactors such as ATP. They were first described to occur in the muscles of bumblebees and brown adipose tissue in the 1970s. A popular example is the conversion of fructose?6-phosphate to fructose?1,6-bisphosphate and back. In the present study, we analyze a large number of substrate cycles in human metabolism that consume ATP and discuss their statistics. For this purpose, we use two recently published methods (i.e. EFMEvolver and the K-shortest EFM method) to calculate samples of 100?000 and 15?000 substrate cycles, respectively. We find an unexpectedly high number of substrate cycles in human metabolism, with up to 100 reactions per cycle, utilizing reactions from up to six different compartments. An analysis of tissue-specific models of liver and brain metabolism shows that there is selective pressure that acts against the uncontrolled dissipation of energy by avoiding the coexpression of enzymes belonging to the same substrate cycle. This selective force is particularly strong against futile cycles that have a high flux as a result of thermodynamic principles.     

Systems properties of the Haemophilus influenzae Rd metabolic genotype.

Haemophilus influenzae Rd was the first free-living organism for which the complete genomic sequence was established. The annotated sequence and known biochemical information was used to define the H. influenzae Rd metabolic genotype. This genotype contains 488 metabolic reactions operating on 343 metabolites. The stoichiometric matrix was used to determine the systems characteristics of the metabolic genotype and to assess the metabolic capabilities of H. influenzae. The need to balance cofactor and biosynthetic precursor production during growth on mixed substrates led to the definition of six different optimal metabolic phenotypes arising from the same metabolic genotype, each with different constraining features. The effects of variations in the metabolic genotype were also studied, and it was shown that the H. influenzae Rd metabolic genotype contains redundant functions under defined conditions. We thus show that the synthesis of in silico metabolic genotypes from annotated genome sequences is possible and that systems analysis methods are available that can be used to analyze and interpret phenotypic behavior of such genotypes.     

The topology of genome-scale metabolic reconstructions unravels independent modules and high network flexibility

The topology of metabolic networks is recognisably modular with modules weakly connected apart from sharing a pool of currency metabolites. Here, we defined modules as sets of reversible reactions isolated from the rest of metabolism by irreversible reactions except for the exchange of currency metabolites. Our approach identifies topologically independent modules under specific conditions associated with different metabolic functions. As case studies, the E.coli iJO1366 and Human Recon 2.2 genome-scale metabolic models were split in 103 and 321 modules respectively, displaying significant correlation patterns in expression data. Finally, we addressed a fundamental question about the metabolic flexibility conferred by reversible reactions: “Of all Directed Topologies (DTs) defined by fixing directions to all reversible reactions, how many are capable of carrying flux through all reactions?”. Enumeration of the DTs for iJO1366 model was performed using an efficient depth-first search algorithm, rejecting infeasible DTs based on mass-imbalanced and loopy flux patterns. We found the direction of 79% of reversible reactions must be defined before all directions in the network can be fixed, granting a high degree of flexibility.     

Crystal structure of dephospho-coenzyme A kinase from Haemophilus influenzae

Dephospho-coenzyme A kinase catalyzes the final step in CoA biosynthesis, the phosphorylation of the 3'-hydroxyl group of ribose using ATP as a phosphate donor. The protein from Haemophilus influenzae was cloned and expressed, and its crystal structure was determined at 2.0-A resolution in complex with ATP. The protein molecule consists of three domains: the canonical nucleotide-binding domain with a five-stranded parallel beta-sheet, the substrate-binding alpha-helical domain, and the lid domain formed by a pair of alpha-helices. The overall topology of the protein resembles the structures of nucleotide kinases. ATP binds in the P-loop in a manner observed in other kinases. The CoA-binding site is located at the interface of all three domains. The double-pocket structure of the substrate-binding site is unusual for nucleotide kinases. Amino acid residues implicated in substrate binding and catalysis have been identified. The structure analysis suggests large domain movements during the catalytic cycle.     

Identification of genome-scale metabolic network models using experimentally measured flux profiles

Genome-scale metabolic network models can be reconstructed for well-characterized organisms using genomic annotation and literature information. However, there are many instances in which model predictions of metabolic fluxes are not entirely consistent with experimental data, indicating that the reactions in the model do not match the active reactions in the in vivo system. We introduce a method for determining the active reactions in a genome-scale metabolic network based on a limited number of experimentally measured fluxes. This method, called optimal metabolic network identification (OMNI), allows efficient identification of the set of reactions that results in the best agreement between in silico predicted and experimentally measured flux distributions. We applied the method to intracellular flux data for evolved Escherichia coli mutant strains with lower than predicted growth rates in order to identify reactions that act as flux bottlenecks in these strains. The expression of the genes corresponding to these bottleneck reactions was often found to be downregulated in the evolved strains relative to the wild-type strain. We also demonstrate the ability of the OMNI method to diagnose problems in E. coli strains engineered for metabolite overproduction that have not reached their predicted production potential. The OMNI method applied to flux data for evolved strains can be used to provide insights into mechanisms that limit the ability of microbial strains to evolve towards their predicted optimal growth phenotypes. When applied to industrial production strains, the OMNI method can also be used to suggest metabolic engineering strategies to improve byproduct secretion. In addition to these applications, the method should prove to be useful in general for reconstructing metabolic networks of ill-characterized microbial organisms based on limited amounts of experimental data.     

The structure of the capsular antigen from Haemophilus influenzae type A.

Structural investigation of the capsular antigen from Haemophilus influenzae type a has shown it to be composed of 4-O-beta-D-glucopyranosyl-D-ribitol residues joined through phosphoric diester linkages between O-4 of D-glucose and O-5 of D-ribitol. Chemical degradations and 13C-n.m.r. spectroscopy were the main methods used.     

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

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

Metabolite coupling in genome-scale metabolic networks

Background  

Biochemically detailed stoichiometric matrices have now been reconstructed for various bacteria, yeast, and for the human cardiac mitochondrion based on genomic and proteomic data. These networks have been manually curated based on legacy data and elementally and charge balanced. Comparative analysis of these well curated networks is now possible. Pairs of metabolites often appear together in several network reactions, linking them topologically. This co-occurrence of pairs of metabolites in metabolic reactions is termed herein "metabolite coupling." These metabolite pairs can be directly computed from the stoichiometric matrix, S. Metabolite coupling is derived from the matrix ŜŜ ^T, whose off-diagonal elements indicate the number of reactions in which any two metabolites participate together, where Ŝ is the binary form of S.     

Haemophilus influenzae Meningitis in Adults

Haemophilus influenzae meningitis though common in childhood is rarely seen in the adult. During the past four years eight cases of H. influenzae meningitis have been seen in St. Thomas''s Hospital and four of these were in patients over 20 years old. There was a possible predisposing condition in two patients. In each case there was difficulty in identification of the organism in the Gram-stained film of the cerebrospinal fluid deposit.