Automated interpretation of metabolic capacity from genome and metagenome sequences

The KEGG pathway maps are widely used as a reference data set for inferring high-level functions of the organism or the ecosystem from its genome or metagenome sequence data. The KEGG modules, which are tighter functional units often corresponding to subpathways in the KEGG pathway maps, are designed for better automation of genome interpretation. Each KEGG module is represented by a simple Boolean expression of KEGG Orthology (KO) identifiers (K numbers), enabling automatic evaluation of the completeness of genes in the genome. Here we focus on metabolic functions and introduce reaction modules for improving annotation and signature modules for inferring metabolic capacity. We also describe how genome annotation is performed in KEGG using the manually created KO database and the computationally generated SSDB database. The resulting KEGG GENES database with KO (K number) annotation is a reference sequence database to be compared for automated annotation and interpretation of newly determined genomes.     

不同生物土壤结皮微生物组跨膜转运蛋白基因多样性及差异

【背景】跨膜转运蛋白在微生物转运各种物质的过程中具有重要作用。【目的】通过比较原核微生物组磷酸转移酶(phosphotransferasesystem,PTS)系统和腺苷三磷酸结合盒(ATP-binding cassette,ABC)转运蛋白编码基因在两种不同生物土壤结皮中(藻结皮与藓结皮)的差异,以期揭示随着生物土壤结皮的发育演替,微生物组跨膜转运物质的生物学过程中的潜在变化趋势。【方法】对腾格里沙漠东南缘的藻结皮和藓结皮12个样品进行宏基因组测序,参照KEGG数据库PTS系统,与ABC转运蛋白代谢通路进行比较并筛选相关基因,分析其差异显著性。【结果】藻结皮和藓结皮PTS系统和ABC转运蛋白编码基因的多样性一致。在生物土壤结皮中共检测到16种PTS系统的转运蛋白的编码基因,具有显著性差异的有5种;检测到106种ABC转运蛋白的编码基因,具有显著性差异的有46种,并对这46种转运蛋白结合的底物以及变化趋势进行了详细的描述。【结论】生物土壤结皮发育演替过程中,微生物组从环境中摄取能够增加渗透势物质的潜力总体呈现降低趋势,转运氨基酸、细胞膜和细胞壁组分的潜力总体呈现增加趋势,对于矿物离子、...     

ABCdb: an online resource for ABC transporter repertories from sequenced archaeal and bacterial genomes

The ATP-binding cassette (ABC) transporters are one of the major classes of active transporters. They are widespread in archaea, bacteria, and eukaryota, indicating that they have arisen early in evolution. They are involved in many essential physiological processes, but the majority import or export a wide variety of compounds across cellular membranes. These systems share a common architecture composed of four (exporters) or five (importers) domains. To identify and reconstruct functional ABC transporters encoded by archaeal and bacterial genomes, we have developed a bioinformatic strategy. Cross-reference to the transport classification system is used to predict the type of compound transported. A high quality of annotation is achieved by manual verification of the predictions. However, in order to face the rapid increase in the number of published genomes, we also include analyses of genomes issuing directly from the automated strategy. Querying the database (http://www-abcdb.biotoul.fr) allows to easily retrieve ABC transporter repertories and related data. Additional query tools have been developed for the analysis of the ABC family from both functional and evolutionary perspectives.     

Towards inference of a biochemical ontology from a metabolic database

In order to predict the metabolic fate of an arbitrary compound based solely on structure, it is useful to be able to identify substructural 'functional groups' that are biochemically reactive. These functional groups are the substructural elements that can be removed and replaced to transform one compound into another. This problem of identifying functional groups is related to the problem of classifying compounds. The research presented here discusses the state of the art in biochemical databases and how these sources may be applied to the problem of classifying compounds based solely on structure. We describe a biochemical informatics system for processing molecular data and describe how 100 255 compositional (hasA) relationships are inferred between 835 abstractions and 9500 metabolites from the KEGG Ligand database. Specifically, we focus on the identification of amino acids and consider ways in which the inference of biochemical ontologies for metabolites will be improved in the future.     

KEGG: kyoto encyclopedia of genes and genomes

KEGG (Kyoto Encyclopedia of Genes and Genomes) is a knowledge base for systematic analysis of gene functions, linking genomic information with higher order functional information. The genomic information is stored in the GENES database, which is a collection of gene catalogs for all the completely sequenced genomes and some partial genomes with up-to-date annotation of gene functions. The higher order functional information is stored in the PATHWAY database, which contains graphical representations of cellular processes, such as metabolism, membrane transport, signal transduction and cell cycle. The PATHWAY database is supplemented by a set of ortholog group tables for the information about conserved subpathways (pathway motifs), which are often encoded by positionally coupled genes on the chromosome and which are especially useful in predicting gene functions. A third database in KEGG is LIGAND for the information about chemical compounds, enzyme molecules and enzymatic reactions. KEGG provides Java graphics tools for browsing genome maps, comparing two genome maps and manipulating expression maps, as well as computational tools for sequence comparison, graph comparison and path computation. The KEGG databases are daily updated and made freely available (http://www. genome.ad.jp/kegg/).     

A nearest neighbor approach for automated transporter prediction and categorization from protein sequences

MOTIVATION: Membrane transport proteins play a crucial role in the import and export of ions, small molecules or macromolecules across biological membranes. Currently, there are a limited number of published computational tools which enable the systematic discovery and categorization of transporters prior to costly experimental validation. To approach this problem, we utilized a nearest neighbor method which seamlessly integrates homologous search and topological analysis into a machine-learning framework. RESULTS: Our approach satisfactorily distinguished 484 transporter families in the Transporter Classification Database, a curated and representative database for transporters. A five-fold cross-validation on the database achieved a positive classification rate of 72.3% on average. Furthermore, this method successfully detected transporters in seven model and four non-model organisms, ranging from archaean to mammalian species. A preliminary literature-based validation has cross-validated 65.8% of our predictions on the 11 organisms, including 55.9% of our predictions overlapping with 83.6% of the predicted transporters in TransportDB.     

A functional genomics approach to establish the complement of carbohydrate transporters in Streptococcus pneumoniae

The aerotolerant anaerobe Streptococcus pneumoniae is part of the normal nasopharyngeal microbiota of humans and one of the most important invasive pathogens. A genomic survey allowed establishing the occurrence of twenty-one phosphotransferase systems, seven carbohydrate uptake ABC transporters, one sodium:solute symporter and a permease, underlining an exceptionally high capacity for uptake of carbohydrate substrates. Despite high genomic variability, combined phenotypic and genomic analysis of twenty sequenced strains did assign the substrate specificity only to two uptake systems. Systematic analysis of mutants for most carbohydrate transporters enabled us to assign a phenotype and substrate specificity to twenty-three transport systems. For five putative transporters for galactose, pentoses, ribonucleosides and sulphated glycans activity was inferred, but not experimentally confirmed and only one transport system remains with an unknown substrate and lack of any functional annotation. Using a metabolic approach, 80% of the thirty-two fermentable carbon substrates were assigned to the corresponding transporter. The complexity and robustness of sugar uptake is underlined by the finding that many transporters have multiple substrates, and many sugars are transported by more than one system. The present work permits to draw a functional map of the complete arsenal of carbohydrate utilisation proteins of pneumococci, allows re-annotation of genomic data and might serve as a reference for related species. These data provide tools for specific investigation of the roles of the different carbon substrates on pneumococcal physiology in the host during carriage and invasive infection.     

Evidence for simultaneous binding of dissimilar substrates by the Escherichia coli multidrug transporter MdfA

The mechanism by which multidrug transporters interact with structurally unrelated substrates remains enigmatic. Based on transport competition experiments, photoaffinity labeling, and effects on enzymatic activities, it was proposed in the past that multidrug transporters can interact simultaneously with a number of dissimilar substrate molecules. To study this phenomenon, we applied a direct binding approach and transport assays using the Escherichia coli multidrug transporter MdfA, which exports both positively charged (e.g., tetraphenylphosphonium, TPP(+)), zwitterionic (e.g., ciprofloxacin), and neutral (e.g., chloramphenicol) drugs. The interaction of MdfA with various substrates was examined by direct binding assays with the purified transporter. The immobilized MdfA binds TPP(+) in a specific manner, and all the tested positively charged substrates inhibit TPP(+) binding. Surprisingly, although TPP(+) binding is not affected by zwitterionic substrates, the neutral substrate chloramphenicol stimulates TPP(+) binding by enhancing its affinity to MdfA. In contrast, transport competition assays show inhibition of TPP(+) transport by chloramphenicol. We suggest that MdfA binds TPP(+) and chloramphenicol simultaneously to distinct but interacting binding sites, and the interaction between these two substrates during transport is discussed.     

Conformational changes in MetNI: steered molecular dynamic studies of the methionine ABC transporter with and without substrates

ATP-binding cassette (ABC) transporters are integral membrane proteins that utilised energy from ATP hydrolysis to translocate substrates across the membrane. In addition to the common nucleotide-binding domains (NBDs) and transmembrane domains (TMDs), the methionine ABC transporter has C-terminal regulatory domains (C2 domains) that belong to ACT protein family. When the amount of methionine in the cell is high, the transport stops. This phenomenon is called trans-inhibition. To understand how a trans-inhibited protein returns to an uninhibited, resting state, we performed steered molecular dynamic simulations with and without the substrates. From the simulations, we observed some important conformational changes in the whole ABC transporter, including the constriction in the translocation pathway in the TMDs and approach of the NBDs. However, the C2 domains behaved differently in two types of the simulations. These findings might help to explain the relationship of the conformational changes of the C2 domains with the rearrangements of the NBDs or TMDs, and provide a way to understand the trans-inhibition from an opposite direction.     

绿色杜氏藻转录组分析

为了深入了解绿色杜氏藻(Dunaliella viridis)基因信息及功能、耐盐相关通路(甘油脂代谢)及关键酶,本文首次通过Illumina HiSeq^TM 2000高通量测序技术对绿色杜氏藻转录组进行测序,利用Trinity软件将数据组装形成转录本,对所有转录本进行COG(Clusters of Orthologous Groups)、GO(Gene Ontology)和KEGG(Kyoto Encyclopedia of Genes and Genomes)分类和功能注释、Pathway注释以及蛋白编码区(Opening reading fragment,ORF)的预测,并对甘油脂代谢通路关键酶基因进行了分析。转录组测序共获得81 593个转录本,其中ORF共有77 117条,约占所有转录本的94.50%。COG分类结果表明,16 569条转录本被分为24个类别。GO分类结果表明,76 436条转录本被注释。在所有注释分类中,生物学过程转录本数量最多,为30 678条,占总转录本数的40.14%。KEGG分析结果表明,317个标准途径中包含26 428条转录本,含转录本最多的类别是代谢,为9949条(37.65%)。与代谢有关的途径为131条,占所有注释途径的41.32%。在甘油脂代谢通路中仅发现1条关键酶转录本(二羟丙酮激酶),该酶可能与绿色杜氏藻耐盐胁迫中甘油的合成有较大关系。本研究进一步完善了绿色杜氏藻的基因信息,为绿色杜氏藻代谢途径研究奠定了坚实的基础。     

EmrE, a multidrug transporter from Escherichia coli, transports monovalent and divalent substrates with the same stoichiometry

Multidrug transporters recognize and transport substrates with apparently little common structural features. At times these substrates are neutral, negatively, or positively charged, and only limited information is available as to how these proteins deal with the energetic consequences of transport of substrates with different charges. Multidrug transporters and drug-specific efflux systems are responsible for clinically significant resistance to chemotherapeutic agents in pathogenic bacteria, fungi, parasites, and human cancer cells. Understanding how these efflux systems handle different substrates may also have practical implications in the development of strategies to overcome the resistance mechanisms mediated by these proteins. Here, we compare transport of monovalent and divalent substrates by EmrE, a multidrug transporter from Escherichia coli, in intact cells and in proteoliposomes reconstituted with the purified protein. The results demonstrated that whereas the transport of monovalent substrates involves charge movement (i.e. electrogenic), the transport of divalent substrate does not (i.e. electroneutral). Together with previous results, these findings suggest that an EmrE dimer exchanges two protons per substrate molecule during each transport cycle. In intact cells, under conditions where the only driving force is the electrical potential, EmrE confers resistance to monovalent substrates but not to divalent ones. In the presence of proton gradients, resistance to both types of substrates is detected. The finding that under some conditions EmrE does not remove certain types of drugs points out the importance of an in-depth understanding of mechanisms of action of multidrug transporters to devise strategies for coping with the problem of multidrug resistance.     

基于转录组分析铜绿假单胞菌DN1降解荧蒽特性

[背景]铜绿假单胞菌DN1是一株从石油污染土壤中分离筛选到的具有广谱降解功能的菌株。[目的]深入了解荧蒽胁迫条件下铜绿假单胞菌DN1降解污染物过程中重要的降解相关基因信息。[方法]通过高通量测序技术对铜绿假单胞菌DN1进行转录组测序,对其所有的转录本进行KEGG (kyoto encyclopedia of genes and genomes)分类和Pathway注释、GO (gene ontology)分类和富集分析。[结果]转录组测序显示:与对照组相比,荧蒽诱导组检测到6 189个基因,其中1 919个基因上调表达,1 603个基因下调表达。KEGG注释分析显示差异上调表达基因匹配到了112个KEGG代谢途径,注释到"代谢途径"的1 408个基因(约占总差异基因的73.4%)中有317个基因参与了碳氢化合物代谢及含有苯环结构的异源生物质的生物降解,占"代谢途径"的16.53%,暗示了菌株DN1降解荧蒽可能与这些途径有密切关系。另外,主要代谢途径中的差异表达基因主要集中在ABC转运系统、氨基酸生物合成、双组分系统及碳代谢,这些途径大多数参与了底物的识别转运、信号转导及基因表达调控。[结论]进一步拓展了铜绿假单胞菌DN1在荧蒽胁迫条件下的代谢途径和逆境反应,也为微生物修复环境污染物研究夯实了理论基础。     

榆树叶片形成虫瘿过程中转录组学分析

榆瘿蚜取食侵染榆树叶片形成了榆树虫瘿,本研究采用新一代的高通量Illumina Hi Seq^TM 2000技术测序平台对榆瘿蚜取食刺激的榆树叶片进行转录组测序和功能注释,利用生物学方法对基因表达和功能进行研究。测序获得23.19 Gb碱基序列信息,通过对测序数据进行序列过滤、拼接和去冗余,共获得102 017个Unigenes,通过NR与BLAST等数据库比对,其中有37 899个(37.15%) Unigense被注释。利用KOG、GO、KEGG等数据库对榆树虫瘿叶片的Unigense进行比对,按功其能将匹配的Unigenes基因划分25大类;GO注释将信息归纳为基因的3大主类,57个亚类;以KEGG数据库为参考,将Unigene定位到110个不同的代谢通路,包括氧化应激防御、植物激素信号转导、碳水化合物以及次生物代谢等代谢相关的Unigenes通路。本研究通过二代高通量转录组测序技术研究榆瘿蚜侵染下榆树虫瘿的相关基因,为今后研究榆瘿蚜侵染榆树叶片形成虫瘿的分子机理提供了基础资料。     

Multidrug efflux transporter, AcrB--the pumping mechanism

Resistance nodulation cell division (RND) transporters are one of the main causes of the bacterial multidrug resistance. They pump a wide range of antibiotics out of the cell by proton motive force. AcrB is the major RND transporter in Escherichia coli. Recently, the crystal structures of AcrB have been determined by different space groups. All these structures are consistent with asymmetric trimer. Each monomer has different conformation corresponding to one of the three functional states of the transport cycle. Transporting hydrophobic drug was bound in the periplasmic domain on one of the three monomers. The transport pathway with alternating access mechanism is located at the hydrophilic domain protruded into the periplasmic space while this mechanism of other transporter families like ATP binding cassette (ABC) and major facilitator superfamily (MFS) transporter is located in the membrane-embedded region. For the RND, protonation might also take place asymmetrically at the functionally important charged residues in the transmembrane (TM) region. The structures indicate that drugs are transported by a three-step functional rotation in which substrates undergo ordered binding change.     

Prediction of transporter targets using efficient RBF networks with PSSM profiles and biochemical properties

SUMMARY: Transporters are proteins that are involved in the movement of ions or molecules across biological membranes. Currently, our knowledge about the functions of transporters is limited due to the paucity of their 3D structures. Hence, computational techniques are necessary to annotate the functions of transporters. In this work, we focused on an important functional aspect of transporters, namely annotation of targets for transport proteins. We have systematically analyzed four major classes of transporters with different transporter targets: (i) electron, (ii) protein/mRNA, (iii) ion and (iv) others, using amino acid properties. We have developed a radial basis function network-based method for predicting transport targets with amino acid properties and position specific scoring matrix profiles. Our method showed a 10-fold cross-validation accuracy of 90.1, 80.1, 70.3 and 82.3% for electron transporters, protein/mRNA transporters, ion transporters and others, respectively, in a dataset of 543 transporters. We have also evaluated the performance of the method with an independent dataset of 108 proteins and we obtained similar accuracy. We suggest that our method could be an effective tool for functional annotation of transport proteins. AVAILABILITY: http://rbf.bioinfo.tw/~sachen/ttrbf.html     

Transport of multidrug resistance substrates by the Streptococcus agalactiae hemolysin transporter

Streptococcus agalactiae (group B streptococcus [GBS]) causes neonatal sepsis, pneumonia, and meningitis, as well as infections of the bovine udder. The S. agalactiae hemolysin is regarded as an important virulence factor, and hemolysin expression is dependent on the cyl gene cluster. cylA and cylB encode the ATP binding and transmembrane domains of a typical ATP binding cassette (ABC) transporter. The deduced proteins contain the signature sequence of a multidrug resistance (MDR) transporter, and mutation of the genes results in a nonhemolytic and nonpigmented phenotype. To further elucidate the function of the putative transporter, nonpolar deletion mutants of cylA were constructed. These mutants are nonhemolytic and can be complemented by the transporter genes. Wild-type strain and nonhemolytic cylA and cylK deletion mutants were exposed to known substrates of MDR transporters. Mutation of cylA significantly impaired growth in the presence of daunorubicin, doxorubicin, and rhodamine 6G and resulted in a decreased export of doxorubicin from the cells. The mutation of cylK, a gene of unknown function located downstream from cylA, caused a loss of hemolysis but had no effect on the transport of MDR substrates. Furthermore, the hemolytic activity of the wild-type strain was inhibited by reserpine in a dose-dependent manner. We conclude that CylAB closely resembles an ABC-type MDR transporter and propose that the GBS hemolysin molecule represents a natural substrate of the transporter.     

MedicCyc: a biochemical pathway database for Medicago truncatula

MOTIVATION: There is an imperative need to integrate functional genomics data to obtain a more comprehensive systems-biology view of the results. We believe that this is best achieved through the visualization of data within the biological context of metabolic pathways. Accordingly, metabolic pathway reconstruction was used to predict the metabolic composition for Medicago truncatula and these pathways were engineered to enable the correlated visualization of integrated functional genomics data. Results: Metabolic pathway reconstruction was used to generate a pathway database for M. truncatula (MedicCyc), which currently features more than 250 pathways with related genes, enzymes and metabolites. MedicCyc was assembled from more than 225,000 M. truncatula ESTs (MtGI Release 8.0) and available genomic sequences using the Pathway Tools software and the MetaCyc database. The predicted pathways in MedicCyc were verified through comparison with other plant databases such as AraCyc and RiceCyc. The comparison with other plant databases provided crucial information concerning enzymes still missing from the ongoing, but currently incomplete M. truncatula genome sequencing project. MedicCyc was further manually curated to remove non-plant pathways, and Medicago-specific pathways including isoflavonoid, lignin and triterpene saponin biosynthesis were modified or added based upon available literature and in-house expertise. Additional metabolites identified in metabolic profiling experiments were also used for pathway predictions. Once the metabolic reconstruction was completed, MedicCyc was engineered to visualize M. truncatula functional genomics datasets within the biological context of metabolic pathways. Availability: freely accessible at http://www.noble.org/MedicCyc/