Mass++: A Visualization and Analysis Tool for Mass Spectrometry.

We have developed Mass++, a plug-in style visualization and analysis tool for mass spectrometry. Its plug-in style enables users to customize it and to develop original functions. Mass++ has several kinds of plug-ins, including rich viewers and analysis methods for proteomics and metabolomics. Plug-ins for supporting vendors' raw data are currently available; hence, Mass++ can read several data formats. Mass++ is both a desktop tool and a software development platform. Original functions can be developed without editing the Mass++ source code. Here, we present this tool's capability to rapidly analyze MS data and develop functions by providing examples of label-free quantitation and implementing plug-ins or scripts. Mass++ is freely available at http://www.first-ms3d.jp/english/ .

A single sex pheromone receptor determines chemical response specificity of sexual behavior in the silkmoth Bombyx mori.

In insects and other animals, intraspecific communication between individuals of the opposite sex is mediated in part by chemical signals called sex pheromones. In most moth species, male moths rely heavily on species-specific sex pheromones emitted by female moths to identify and orient towards an appropriate mating partner among a large number of sympatric insect species. The silkmoth, Bombyx mori, utilizes the simplest possible pheromone system, in which a single pheromone component, (E, Z)-10,12-hexadecadienol (bombykol), is sufficient to elicit full sexual behavior. We have previously shown that the sex pheromone receptor BmOR1 mediates specific detection of bombykol in the antennae of male silkmoths. However, it is unclear whether the sex pheromone receptor is the minimally sufficient determination factor that triggers initiation of orientation behavior towards a potential mate. Using transgenic silkmoths expressing the sex pheromone receptor PxOR1 of the diamondback moth Plutella xylostella in BmOR1-expressing neurons, we show that the selectivity of the sex pheromone receptor determines the chemical response specificity of sexual behavior in the silkmoth. Bombykol receptor neurons expressing PxOR1 responded to its specific ligand, (Z)-11-hexadecenal (Z11-16:Ald), in a dose-dependent manner. Male moths expressing PxOR1 exhibited typical pheromone orientation behavior and copulation attempts in response to Z11-16:Ald and to females of P. xylostella. Transformation of the bombykol receptor neurons had no effect on their projections in the antennal lobe. These results indicate that activation of bombykol receptor neurons alone is sufficient to trigger full sexual behavior. Thus, a single gene defines behavioral selectivity in sex pheromone communication in the silkmoth. Our findings show that a single molecular determinant can not only function as a modulator of behavior but also as an all-or-nothing initiator of a complex species-specific behavioral sequence.

Pheromonal activities of the bombykol isomer, (10E,12E)-10,12-hexadecadien-1-ol, in the pheromone gland of the silkmoth Bombyx mori.

Bombykol (EZ) is the single component of the female sex pheromone in the silkmoth Bombyx mori. EZ alone evokes full courtship behaviors from conspecific males; however, its geometric isomer (EE) was consistently detected in the pheromone glands (PG) of 16 B. mori strains and a field population of the wild silkmoth Bombyx mandarina, which also uses EZ as the single pheromone component. We investigated the pheromonal activities of EE using a commercial hybrid strain of B. mori, Kinshu × Showa. The behavioral assay demonstrated that a 104-105-fold higher dose of EE than EZ was able to elicit behavioral responses from males. To elucidate whether the trace contaminant of EZ in the EE standard is responsible for these responses, we examined the responses of male antennae to EE using a gas chromatograph-electroantennographic detector system (GC-EAD). The EE, at high doses elicited marginal responses from the male antennae. We next examined antennal and behavioral responses of B. mori whose BmOR1 gene, which is responsible for the reception of bombykol, was knocked out. The knockout of BmOR1 resulted in the complete loss of antennal and behavioral responses to EE and EZ, demonstrating that if EE itself is active, it induces these responses via the incidental stimulation of BmOR1, not via the stimulation of EE-specific receptors. The existence of EE in the PG of B. mori and B. mandarina is discussed from the viewpoints of pheromone biosynthesis and the evolution of pheromone communication systems.

Characterization of the inhibitor binding site in mitochondrial NADH-ubiquinone oxidoreductase by photoaffinity labeling using a quinazoline-type inhibitor.

The diverse inhibitors of bovine heart mitochondrial complex I (NADH-ubiquinone oxidoreductase) are believed to share a common large binding domain with partially overlapping sites, though it remains unclear how these binding sites relate to each other. To obtain new insight into the inhibitor binding domain in complex I, we synthesized a photoreactive azidoquinazoline {[(125)I]-6-azido-4-(4-iodophenethylamino)quinazoline, [(125)I]AzQ}, in which a photolabile azido group was introduced into the toxophoric quinazoline ring to allow specific cross-linking, and carried out a photoaffinity labeling study using bovine heart submitochondrial particles. Analysis of the photo-cross-linked proteins by peptide mass fingerprinting and immunoblotting revealed that [(125)I]AzQ specifically binds to the 49 kDa and ND1 subunits with a frequency of approximately 4:1. The cross-linking was completely blocked by excess amounts of other inhibitors such as acetogenin and fenpyroximate. Considerable cross-linking was also detected in the ADP/ATP carrier and 3-hydroxybutyrate dehydrogenase, though it was not associated with dysfunction of the two proteins. The partial proteolysis of the [(125)I]AzQ-labeled 49 kDa subunit by V8-protease and N-terminal sequencing of the resulting peptides revealed that the amino acid residue cross-linked by [(125)I]AzQ is within the sequence region Thr25-Glu143 (118 amino acids). Furthermore, examination of fragment patterns generated by exhaustive digestion of the [(125)I]AzQ-labeled 49 kDa subunit by V8-protease, lysylendopeptidase, or trypsin strongly suggested that the cross-linked residue is located within the region Asp41-Arg63 (23 amino acids). The present study has revealed, for the first time, the inhibitor binding site in complex I at the sub-subunit level.

The tryptophan pathway is involved in the defense responses of rice against pathogenic infection via serotonin production.

The upregulation of the tryptophan (Trp) pathway in rice leaves infected by Bipolaris oryzae was indicated by: (i) enhanced enzyme activity of anthranilate synthase (AS), which regulates metabolic flux in the Trp pathway; (ii) elevated levels of the AS (OASA2, OASB1, and OASB2) transcripts; and (iii) increases in the contents of anthranilate, indole, and Trp. The measurement of the contents of Trp-derived metabolites by high-performance liquid chromatography coupled with tandem mass spectrometry revealed that serotonin and its hydroxycinnamic acid amides were accumulated in infected leaves. Serotonin accumulation was preceded by a transient increase in the tryptamine content and by marked activation of Trp decarboxylase, indicating that enhanced Trp production is linked to the formation of serotonin from Trp via tryptamine. Feeding of radiolabeled serotonin to inoculated leaves demonstrated that serotonin is incorporated into the cell walls of lesion tissue. The leaves of a propagating-type lesion mimic mutant (sl, Sekiguchi lesion) lacked both serotonin production and deposition of unextractable brown material at the infection sites, and showed increased susceptibility to B. oryzae infection. Treating the mutant with serotonin restored deposition of brown material at the lesion site. In addition, the serotonin treatment suppressed the growth of fungal hyphae in the leaf tissues of the sl mutant. These findings indicated that the activation of the Trp pathway is involved in the establishment of effective physical defenses by producing serotonin in rice leaves.

Exogenous polyamines elicit herbivore-induced volatiles in lima bean leaves: involvement of calcium, H2O2 and Jasmonic acid.

We investigated the role of polyamines (PAs) in lima bean (Phaseolus lunatus) leaves on the production of herbivorous mite (Tetranychus urticae)-induced plant volatiles that attract carnivorous natural enemies of the herbivores. To do this, we focused on the effects of the exogenous PAs [cadaverine, putrescine, spermidine and spermine (Spm)] on the production of volatiles, H(2)O(2) and jasmonic acid (JA) and the levels of defensive genes, cytosolic calcium and reactive oxygen species (ROS). Among the tested PAs, Spm was the most active in inducing the production of volatile terpenoids known to be induced by T. urticae. An increase in JA levels was also found after Spm treatment, indicating that Spm induces the biosynthesis of JA, which has been shown elsewhere to regulate the production of some volatile terpenoids. Further, treatment with JA and Spm together resulted in greater volatile emission than that with JA alone. In a Y-tube olfactometer, leaves treated with Spm + JA attracted more predatory mites (Phytoseiulus persimilis) than those treated with JA alone. After treatment with Spm + JA, no effects were found on the enzyme activity of polyamine oxidase and copper amine oxidase. However, induction of calcium influx and ROS production, and increased enzyme activities and gene expression for NADPH oxidase complex, superoxide dismutase, catalase, ascorbate peroxidase, glutathione reductase and glutathione peroxidase were found after treatment with Spm + JA. These results indicate that Spm plays an important role in the production of T. urticae-induced lima bean leaf volatiles.

Accumulation of hydroxycinnamic acid amides induced by pathogen infection and identification of agmatine coumaroyltransferase in Arabidopsis thaliana.

Hydroxycinnamic acid amides (HCAAs) are secondary metabolites involved in the defense of plants against pathogens. Here, we report the first identification of HCAAs, p-coumaroylagmatine, feruloylagmatine, p-coumaroylputrescine and feruloylputrescine, in Arabidopsis thaliana rosette leaves infected with Alternaria brassicicola and the assignment of At5g61160 as the agmatine coumaroyltransferase (AtACT) that catalyzes the last reaction in the biosynthesis of the HCAAs. Feeding experiments with putative labeled precursors revealed that the four HCAAs were synthesized from hydroxycinnamic acids and agmatine or putrescine. AtACT gene function was identified from an analysis of a mutant that did not accumulate HCAAs. In wild-type Arabidopsis, AtACT transcripts markedly increased in response to A. brassicicola infection. Enzymatic activity that catalyzes the synthesis of the HCAAs was confirmed in vitro by using a recombinant AtACT expressed in Escherichia coli. The Atact mutant was susceptible to infection by A. brassicicola, indicating that HCAAs are responsible for defense against pathogens in A. thaliana.

Synthesis and characterization of new piperazine-type inhibitors for mitochondrial NADH-ubiquinone oxidoreductase (complex I).

The mode of action of Deltalac-acetogenins, strong inhibitors of bovine heart mitochondrial complex I, is different from that of traditional inhibitors such as rotenone and piericidin A [Murai, M., et al. (2007) Biochemistry 46 , 6409-6416]. As further exploration of these unique inhibitors might provide new insights into the terminal electron transfer step of complex I, we drastically modified the structure of Deltalac-acetogenins and characterized their inhibitory action. In particular, on the basis of structural similarity between the bis-THF and the piperazine rings, we here synthesized a series of piperazine derivatives. Some of the derivatives exhibited very potent inhibition at nanomolar levels. The hydrophobicity of the side chains and their balance were important structural factors for the inhibition, as is the case for the original Deltalac-acetogenins. However, unlike in the case of the original Deltalac-acetogenins, (i) the presence of two hydroxy groups is not crucial for the activity, (ii) the level of superoxide production induced by the piperazines is relatively high, (iii) the inhibitory potency for the reverse electron transfer is remarkably weaker than that for the forward event, and (iv) the piperazines efficiently suppressed the specific binding of a photoaffinity probe of natural-type acetogenins ([ (125)I]TDA) to the ND1 subunit. We therefore conclude that the action mechanism of the piperazine series differs from that of the original Deltalac-acetogenins. The photoaffinity labeling study using a newly synthesized photoreactive piperazine ([ (125)I]AFP) revealed that this compound binds to the 49 kDa subunit and an unidentified subunit, not ND1, with a frequency of approximately 1:3. A variety of traditional complex I inhibitors as well as Deltalac-acetogenins suppressed the specific binding of [ (125)I]AFP to the subunits. The apparent competitive behavior of inhibitors that seem to bind to different sites may be due to structural changes at the binding site, rather than occupying the same site. The meaning of the occurrence of diverse inhibitors exhibiting different mechanisms of action is discussed in light of the functionality of the membrane arm of complex I.

Dynamic function of the spacer region of acetogenins in the inhibition of bovine mitochondrial NADH-ubiquinone oxidoreductase (complex I).

Studies of the action mechanism of acetogenins, the most potent and structurally unique inhibitors of bovine heart mitochondrial complex I (NADH-ubiquinone oxidoreductase), are valuable in characterizing the inhibitor binding site in this enzyme. Our previous study deepened our understanding of the dynamic function of the spacer region of bis-THF acetogenins [Abe, M., et al. (2005) Biochemistry 44, 14898-14906] but, at the same time, posed new important questions. First, while the two toxophores (i.e., the hydroxylated THF and the gamma-lactone rings) span a distance shorter than that of the extended 13 carbon atoms [-(CH 2) 13-], what is the apparent optimal length of the spacer for the inhibition of 13 carbon atoms? In other words, what is the functional role of the additional methylene groups? Second, why was the inhibitory potency of the mono-THF derivative, but not the bis-THF derivative, drastically reduced by hardening the spacer covering 10 carbon atoms into a rodlike shape [-CH 2-(C identical withC) 4-CH 2-]? This study was designed not only to answer these questions but also to further disclose the dynamic functions of the spacer. We here synthesized systematically designed acetogenins, including mono- and bis-THF derivatives, and evaluated their inhibitory effects on bovine complex I. With regard to the first question, we demonstrated that the additional methylenes enhance the hydrophobicity of the spacer region, which may be thermodynamically advantageous for bringing the polar gamma-lactone ring into the membrane-embedded segment of complex I. With regard to the second question, we observed that a decrease in the flexibility of the spacer region is more adverse to the action of the mono-THF series than that of the bis-THF series. As a cause of this difference, we suggest that for bis-THF derivatives, one of the two THF rings, being adjacent to the spacer, is capable of working as a pseudospacer to overcome the remarkable decrease in the conformational freedom and/or the length of the spacer. Moreover, using photoresponsive acetogenins that undergo drastic and reversible conformational changes with alternating UV-vis irradiation, we provided further evidence that the spacer region is free from steric congestion arising from the putative binding site probably because there is no receptor wall for the spacer region.

Identification of receptors of main sex-pheromone components of three Lepidopteran species.

Male moths discriminate conspecific female-emitted sex pheromones. Although the chemical components of sex pheromones have been identified in more than 500 moth species, only three components in Bombyx mori and Heliothis virescens have had their receptors identified. Here we report the identification of receptors for the main sex-pheromone components in three moth species, Plutella xylostella, Mythimna separata and Diaphania indica. We cloned putative sex-pheromone receptor genes PxOR1, MsOR1 and DiOR1 from P. xylostella, M. separata and D. indica, respectively. Each of the three genes was exclusively expressed with an Or83b orthologous gene in male olfactory receptor neurons (ORNs) that are surrounded by supporting cells expressing pheromone-binding-protein (PBP) genes. By two-electrode voltage-clamp recording, we tested the ligand specificity of Xenopus oocytes co-expressing PxOR1, MsOR1 or DiOR1 with an OR83b family protein. Among the seven sex-pheromone components of the three moth species, the oocytes dose-dependently responded only to the main sex-pheromone component of the corresponding moth species. In our study, PBPs were not essential for ligand specificity of the receptors. On the phylogenetic tree of insect olfactory receptors, the six sex-pheromone receptors identified in the present and previous studies are grouped in the same subfamily but have no relation with the taxonomy of moths. It is most likely that sex-pheromone receptors have randomly evolved from ancestral sex-pheromone receptors before the speciation of moths and that their ligand specificity was modified by mutations of local amino acid sequences after speciation.

Novel Approach to Classify Plants Based on Metabolite-Content Similarity.

Secondary metabolites are bioactive substances with diverse chemical structures. Depending on the ecological environment within which they are living, higher plants use different combinations of secondary metabolites for adaptation (e.g., defense against attacks by herbivores or pathogenic microbes). This suggests that the similarity in metabolite content is applicable to assess phylogenic similarity of higher plants. However, such a chemical taxonomic approach has limitations of incomplete metabolomics data. We propose an approach for successfully classifying 216 plants based on their known incomplete metabolite content. Structurally similar metabolites have been clustered using the network clustering algorithm DPClus. Plants have been represented as binary vectors, implying relations with structurally similar metabolite groups, and classified using Ward's method of hierarchical clustering. Despite incomplete data, the resulting plant clusters are consistent with the known evolutional relations of plants. This finding reveals the significance of metabolite content as a taxonomic marker. We also discuss the predictive power of metabolite content in exploring nutritional and medicinal properties in plants. As a byproduct of our analysis, we could predict some currently unknown species-metabolite relations.

Metabolomic Studies of Indonesian Jamu Medicines: Prediction of Jamu Efficacy and Identification of Important Metabolites.

In order to obtain a better understanding why some Jamu formulas can be used to treat a specific disease, we performed metabolomic studies of Jamu by taking into consideration the biologically active compounds existing in plants used as Jamu ingredients. A thorough integration of information from omics is expected to provide solid evidence-based scientific rationales for the development of modern phytomedicines. This study focused on prediction of Jamu efficacy based on its component metabolites and also identification of important metabolites related to each efficacy group. Initially, we compared the performance of Support Vector Machines and Random Forest to predict the Jamu efficacy with three different data pre-processing approaches, such as no filtering, Single Filtering algorithm, and a combination of Single Filtering algorithm and feature selection using Regularized Random Forest. Both classifiers performed very well and according to 5-fold cross-validation results, the mean accuracy of Support Vector Machine with linear kernel was slightly better than Random Forest. It can be concluded that machine learning methods can successfully relate Jamu efficacy with metabolites. In addition, we extended our analysis by identifying important metabolites from the Random Forest model. The inTrees framework was used to extract the rules and to select important metabolites for each efficacy group. Overall, we identified 94 significant metabolites associated to 12 efficacy groups and many of them were validated by published literature and KNApSAcK Metabolite Activity database.

MathDAMP: a package for differential analysis of metabolite profiles.

BACKGROUND: With the advent of metabolomics as a powerful tool for both functional and biomarker discovery, the identification of specific differences between complex metabolite profiles is becoming a major challenge in the data analysis pipeline. The task remains difficult, given the datasets' size, complexity, and common shifts in migration (elution/retention) times between samples analyzed by hyphenated mass spectrometry methods. RESULTS: We present a Mathematica (Wolfram Research, Inc.) package MathDAMP (Mathematica package for Differential Analysis of Metabolite Profiles), which highlights differences between raw datasets acquired by hyphenated mass spectrometry methods by applying arithmetic operations to all corresponding signal intensities on a datapoint-by-datapoint basis. Peak identification and integration is thus bypassed and the results are displayed graphically. To facilitate direct comparisons, the raw datasets are automatically preprocessed and normalized in terms of both migration times and signal intensities. A combination of dynamic programming and global optimization is used for the alignment of the datasets along the migration time dimension. The processed datasets and the results of direct comparisons between them are visualized using density plots (axes represent migration time and m/z values while peaks appear as color-coded spots) providing an intuitive overall view. Various forms of comparisons and statistical tests can be applied to highlight subtle differences. Overlaid electropherograms (chromatograms) corresponding to the vicinities of the candidate differences from any result may be generated in a descending order of significance for visual confirmation. Additionally, a standard library table (a list of m/z values and migration times for known compounds) may be aligned and overlaid on the plots to allow easier identification of metabolites. CONCLUSION: Our tool facilitates the visualization and identification of differences between complex metabolite profiles according to various criteria in an automated fashion and is useful for data-driven discovery of biomarkers and functional genomics.

Metabolic changes in Arabidopsis thaliana expressing the feedback-resistant anthranilate synthase alpha subunit gene OASA1D.

Anthranilate synthase (AS) is a key enzyme in tryptophan (Trp) biosynthesis. Metabolic changes in transgenic Arabidopsis plants expressing the feedback-resistant anthranilate synthase alpha subunit gene OASA1D were investigated with respect to Trp synthesis and effects on secondary metabolism. The Trp content varied depending on the transgenic line, with some lines showing an approximately 200-fold increase. The levels of AS activity in crude extracts from the transgenic lines were comparable to those in the wild type. On the other hand, the enzyme prepared from the lines accumulating high levels of Trp showed a relaxed feedback sensitivity. The AS activity, determined in the presence of 50 microM L-Trp, correlated well with the amount of free Trp in the transgenic lines, indicating the important role of feedback inhibition in control of Trp pool size. In Arabidopsis, Trp is a precursor of multiple secondary metabolites, including indole glucosinolates and camalexin. The amount of indol-3-ylmethyl glucosinolate (I3 M) in rosette leaves of the high-Trp accumulating lines was 1.5- to 2.1-fold greater than that in wild type. The treatment of the leaves with jasmonic acid resulted in a more pronounced accumulation of I3 M in the high-Trp accumulating lines than in wild type. The induction of camalexin formation after the inoculation of Alternaria brassicicola was not affected by the accumulation of a large amount of Trp. The accumulation of constitutive phenylpropanoids and flavonoids was suppressed in high-Trp accumulating lines, while the amounts of Phe and Tyr increased, thereby indicating an interaction between the Trp branch and the Phe and Tyr branch in the shikimate pathway.

LIGAND: database of chemical compounds and reactions in biological pathways.

LIGAND is a composite database comprising three sections: COMPOUND for the information about metabolites and other chemical compounds, REACTION for the collection of substrate-product relations representing metabolic and other reactions, and ENZYME for the information about enzyme molecules. The current release (as of September 7, 2001) includes 7298 compounds, 5166 reactions and 3829 enzymes. In addition to the keyword search provided by the DBGET/LinkDB system, a substructure search to the COMPOUND and REACTION sections is now available through the World Wide Web (http://www.genome.ad.jp/ligand/). LIGAND may be also downloaded by anonymous FTP (ftp://ftp.genome.ad.jp/pub/kegg/ligand/).

Characterization of inhibitor binding sites of mitochondrial complex I using fluorescent inhibitor.

Recent progress in complex I research suggests that a wide variety of complex I inhibitors share a common large binding domain with partially overlapping sites. To verify this concept, we carried out real-time displacement tests of a fluorescent ligand with various competitors using a novel quinazoline-type inhibitor (aminoquinazoline, AQ). In the presence of an excess amount of the competitors, the binding of AQ to the enzyme was completely suppressed, being in line with the concept mentioned above. However, AQ bound to the enzyme was not displaced by subsequent addition of an increasing amount of competitors in the concentration range expected from the relative magnitude of the K(d) values of AQ and competitors, rather, much higher concentrations of the competitors were needed to displace bound AQ. These results cannot be explained merely by the premise of a common or partially overlapping binding site(s) between AQ and competitors. On the other hand, double-inhibitor titration of steady state complex I activity suggested that additivity of inhibition is not necessarily observed for all pairs of complex I inhibitors. Our results are discussed in light of the cooperativity of the inhibitor binding sites.

Development and mining of a volatile organic compound database.

Volatile organic compounds (VOCs) are small molecules that exhibit high vapor pressure under ambient conditions and have low boiling points. Although VOCs contribute only a small proportion of the total metabolites produced by living organisms, they play an important role in chemical ecology specifically in the biological interactions between organisms and ecosystems. VOCs are also important in the health care field as they are presently used as a biomarker to detect various human diseases. Information on VOCs is scattered in the literature until now; however, there is still no available database describing VOCs and their biological activities. To attain this purpose, we have developed KNApSAcK Metabolite Ecology Database, which contains the information on the relationships between VOCs and their emitting organisms. The KNApSAcK Metabolite Ecology is also linked with the KNApSAcK Core and KNApSAcK Metabolite Activity Database to provide further information on the metabolites and their biological activities. The VOC database can be accessed online.

Metabolonote: a wiki-based database for managing hierarchical metadata of metabolome analyses.

Metabolomics - technology for comprehensive detection of small molecules in an organism - lags behind the other "omics" in terms of publication and dissemination of experimental data. Among the reasons for this are difficulty precisely recording information about complicated analytical experiments (metadata), existence of various databases with their own metadata descriptions, and low reusability of the published data, resulting in submitters (the researchers who generate the data) being insufficiently motivated. To tackle these issues, we developed Metabolonote, a Semantic MediaWiki-based database designed specifically for managing metabolomic metadata. We also defined a metadata and data description format, called "Togo Metabolome Data" (TogoMD), with an ID system that is required for unique access to each level of the tree-structured metadata such as study purpose, sample, analytical method, and data analysis. Separation of the management of metadata from that of data and permission to attach related information to the metadata provide advantages for submitters, readers, and database developers. The metadata are enriched with information such as links to comparable data, thereby functioning as a hub of related data resources. They also enhance not only readers' understanding and use of data but also submitters' motivation to publish the data. The metadata are computationally shared among other systems via APIs, which facilitate the construction of novel databases by database developers. A permission system that allows publication of immature metadata and feedback from readers also helps submitters to improve their metadata. Hence, this aspect of Metabolonote, as a metadata preparation tool, is complementary to high-quality and persistent data repositories such as MetaboLights. A total of 808 metadata for analyzed data obtained from 35 biological species are published currently. Metabolonote and related tools are available free of cost at http://metabolonote.kazusa.or.jp/.

Finding evolutionary relations beyond superfamilies: fold-based superfamilies.

Superfamily classifications are based variably on similarity of sequences, global folds, local structures, or functions. We have examined the possibility of defining superfamilies purely from the viewpoint of the global fold/function relationship. For this purpose, we first classified protein domains according to the beta-sheet topology. We then introduced the concept of kinship relations among the classified beta-sheet topology by assuming that the major elementary event leading to creation of a new beta-sheet topology is either an addition or deletion of one beta-strand at the edge of an existing beta-sheet during the molecular evolution. Based on this kinship relation, a network of protein domains was constructed so that the distance between a pair of domains represents the number of evolutionary events that lead one from the other domain. We then mapped on it all known domains with a specific core chemical function (here taken, as an example, that involving ATP or its analogs). Careful analyses revealed that the domains are found distributed on the network as >20 mutually disjointed clusters. The proteins in each cluster are defined to form a fold-based superfamily. The results indicate that >20 ATP-binding protein superfamilies have been invented independently in the process of molecular evolution, and the conservative evolutionary diffusion of global folds and functions is the origin of the relationship between them.

Analysis of carboxylic acid metabolites from the tricarboxylic acid cycle in Bacillus subtilis cell extract by capillary electrophoresis using an indirect photometric detection method.

With a growing interest in metabolome analysis, there is a need for developing robust methods for analysis of intracellular metabolites profiles in real samples like e.g., bacteria cell. Due to their weak absorbance properties, tri- and dicarboxylic acids from TCA cycle (citric, isocitric, 2-oxoglutaric, succinic, fumaric, malic) as well as carboxylic acid metabolites from glycolysis pathway, urea cycle and metabolism of amino compounds (formic, pyruvic, lactic, acetic, glutamic) were analyzed by capillary electrophoresis (CE) with indirect UV detection. Using 4 mM 2,6-pyridinedicarboxylic acid as a highly UV absorbing carrier electrolyte, 0.2 mM cetyltrimethylammonium bromide, 10% ethylene glycol and 10% acetonitrile, pH 3.5, carboxylic acids metabolites were analyzed in Bacillus subtilis cell extract from two different cultures: glucose and malate. CE with an electrokinetic injection mode achieved limits of detection in the range of 13-54 ppb (1.12-10(-7) - 5.96-10(-7) M). The reproducibility and linearity of method was investigated with RSD for migration time less than 1.3% and acceptable correlation coefficients. The optimized CE method was used to compare metabolome content of cell extract derived from two different culture media containing either glucose or malate as a carbon source. The changes in carboxylic acid metabolites profile were observed depending from used culture medium. Carboxylic acid concentrations ranged: in cell extract from malate culture from 59 to 0.5 microM for lactate and citrate, respectively, and in cell extract from glucose culture from 133 to 0.5 microM for glutamate and citrate, respectively. Appropriate concentrations of carboxylic acid in the single bacterium cell were estimated at mM and sub-mM levels.