Malate/lactate dehydrogenase in mollicutes: evidence for a multienzyme protein

A "plus one" strategy impacts replication of felid alphaherpesvirus 1, Mycoplasma and Chlamydia, and the metabolism of coinfected feline cells

Coinfections are known to play an important role in disease progression and severity. Coinfections are common in cats, but no coinfection studies have investigated the in vitro dynamics between feline viral and bacterial pathogens. In this study, we performed co-culture and invasion assays to investigate the ability of common feline bacterial respiratory pathogens, Chlamydia felis and Mycoplasma felis, to replicate in and invade into Crandell-Rees feline kidney cells. We subsequently investigated how coinfection of these feline cells with each bacterium (C. felis or M. felis) and the common feline viral pathogen, felid alphaherpesvirus 1 (FHV-1), affects replication of each agent in this cell culture system. We also investigated the metabolic impact of each co-pathogen using metabolomic analysis of infected and coinfected cells. C. felis was able to invade and replicate in CRFKs, while M. felis had little capacity to invade. During coinfection, FHV-1 replication was minimally affected by the presence of either bacterial pathogen, but bacterial replication kinetics were more affected, particularly in M. felis. Both C. felis and M. felis replicated to higher levels in the presence of a secondary pathogen. Coinfections resulted in reprogramming of the glycolysis pathway, the pentose phosphate pathway, and the tricarboxylic acid cycle. The distinct metabolic profiles of coinfected cells compared to those of cells infected with just one of these three pathogens, as well as the impact of coinfections on viral or bacterial load, suggest strong interactions between these three pathogens and possible synergistic mechanisms enhancing virulence that need further investigation.IMPORTANCEIn the natural world, respiratory pathogens coexist within their hosts, but their dynamics and interactions remain largely unexplored. Herpesviruses, mycoplasmas, and chlamydias are common and significant causes of acute and chronic respiratory and system disease in animals and people, and these diseases are increasingly found to be polymicrobial. This study investigates how coinfection of feline cells between three respiratory pathogens of cats impact each other as well as the host innate metabolic response to infection. Each of these pathogens have been implicated in the induction of feline upper respiratory tract disease in cats, which is the leading cause of euthanasia in shelters. Understanding how coinfection impacts co-pathogenesis and host responses is critical for improving disease management.

Reversion of mutations in a live mycoplasma vaccine alters its metabolism

The live attenuated temperature sensitive vaccine strain MS-H (Vaxsafe® MS, Bioproperties Pty. Ltd., Australia) is widely used to control disease associated with M. synoviae infection in commercial poultry. MS-H was derived from a field strain (86079/7NS) through N-methyl-N'-nitro-N-nitrosoguanidine (NTG)-induced mutagenesis. Whole genomic sequence analysis of the MS-H and comparison with that of the 86079/7NS have found that MS-H contains 32 single nucleotide polymorphisms (SNPs). Three of these SNPs, found in the obgE, oppF and gapdh genes, have been shown to be prone to reversion under field condition, albeit at a low frequency. Three MS-H reisolates containing the 86079/7NS genotype in obgE (AS2), obgE and oppF (AB1), and obgE, oppF and gapdh (TS4), appeared to be more immunogenic and transmissible compared to MS-H in chickens. To investigate the influence of these reversions in the in vitro fitness of M. synoviae, the growth kinetics and steady state metabolite profiles of the MS-H reisolates, AS2, AB1 and TS4, were compared to those of the vaccine strain. Steady state metabolite profiling of the reisolates showed that changes in ObgE did not significantly influence the metabolism, while changes in OppF was associated with significant alterations in uptake of peptides and/or amino acids into the M. synoviae cell. It was also found that GAPDH plays a role in metabolism of the glycerophospholipids as well as an arginine deiminase (ADI) pathway. This study underscores the role of ObgE, OppF and GAPDH in M. synoviae metabolism, and suggests that the impaired fitness arising from variations in ObgE, OppF and GAPDH contributes to attenuation of MS-H.

Analysis of the Mycoplasma bovis lactate dehydrogenase reveals typical enzymatic activity despite the presence of an atypical catalytic site motif

The lactate dehydrogenase (LDH) of Mycoplasma genitalium has been predicted to also act as a malate dehydrogenase (MDH), but there has been no experimental validation of this hypothesized dual function for any mollicute. Our analysis of the metabolite profile of Mycoplasma bovis using gas chromatography/mass spectrometry (GC/MS) and liquid chromatography/mass spectrometry (LC/MS) detected malate, suggesting that there may be MDH activity in M. bovis. To investigate whether the putative l-LDH enzyme of M. bovis has a dual function (MDH and LDH), we performed bioinformatic and functional biochemical analyses. Although the amino acid sequence and predicted structural analysis of M. bovisl-LDH revealed unusual residues within the catalytic site, suggesting that it may have the flexibility to possess a dual function, our biochemical studies using recombinant M. bovis -LDH did not detect any MDH activity. However, we did show that the enzyme has typical LDH activity that could be inhibited by both MDH substrates oxaloacetate (OAA) and malate, suggesting that these substrates may be able to bind to M. bovis LDH. Inhibition of the conversion of pyruvate to lactate by OAA may be one method the mycoplasma cell uses to reduce the potential for accumulation of intracellular lactate.

Comparative Metabolomics of Mycoplasma bovis and Mycoplasma gallisepticum Reveals Fundamental Differences in Active Metabolic Pathways and Suggests Novel Gene Annotations

Mycoplasmas are pathogenic bacteria that cause serious chronic infections in production animals, resulting in considerable losses worldwide, as well as causing disease in humans. These bacteria have extremely reduced genomes and are thought to have limited metabolic flexibility, even though they are highly successful persistent parasites in a diverse number of species. The extent to which different Mycoplasma species are capable of catabolizing host carbon sources and nutrients, or synthesizing essential metabolites, remains poorly defined. We have used advanced metabolomic techniques to identify metabolic pathways that are active in two species of Mycoplasma that infect distinct hosts (poultry and cattle). We show that these species exhibit marked differences in metabolite steady-state levels and carbon source utilization. This information has been used to functionally characterize previously unknown genes in the genomes of these pathogens. These species-specific differences are likely to reflect important differences in host nutrient levels and pathogenic mechanisms.

Mycoplasmas are simple, but successful parasites that have the smallest genome of any free-living cell and are thought to have a highly streamlined cellular metabolism. Here, we have undertaken a detailed metabolomic analysis of two species, Mycoplasma bovis and Mycoplasma gallisepticum, which cause economically important diseases in cattle and poultry, respectively. Untargeted gas chromatography-mass spectrometry and liquid chromatography-mass spectrometry analyses of mycoplasma metabolite extracts revealed significant differences in the steady-state levels of many metabolites in central carbon metabolism, while ¹³C stable isotope labeling studies revealed marked differences in carbon source utilization. These data were mapped onto in silico metabolic networks predicted from genome wide annotations. The analyses elucidated distinct differences, including a clear difference in glucose utilization, with a marked decrease in glucose uptake and glycolysis in M. bovis compared to M. gallisepticum, which may reflect differing host nutrient availabilities. The ¹³C-labeling patterns also revealed several functional metabolic pathways that were previously unannotated in these species, allowing us to assign putative enzyme functions to the products of a number of genes of unknown function, especially in M. bovis. This study demonstrates the considerable potential of metabolomic analyses to assist in characterizing significant differences in the metabolism of different bacterial species and in improving genome annotation.

IMPORTANCE Mycoplasmas are pathogenic bacteria that cause serious chronic infections in production animals, resulting in considerable losses worldwide, as well as causing disease in humans. These bacteria have extremely reduced genomes and are thought to have limited metabolic flexibility, even though they are highly successful persistent parasites in a diverse number of species. The extent to which different Mycoplasma species are capable of catabolizing host carbon sources and nutrients, or synthesizing essential metabolites, remains poorly defined. We have used advanced metabolomic techniques to identify metabolic pathways that are active in two species of Mycoplasma that infect distinct hosts (poultry and cattle). We show that these species exhibit marked differences in metabolite steady-state levels and carbon source utilization. This information has been used to functionally characterize previously unknown genes in the genomes of these pathogens. These species-specific differences are likely to reflect important differences in host nutrient levels and pathogenic mechanisms.

Insights on the virulence of swine respiratory tract mycoplasmas through genome-scale metabolic modeling

Background

The respiratory tract of swine is colonized by several bacteria among which are three Mycoplasma species: Mycoplasma flocculare, Mycoplasma hyopneumoniae and Mycoplasma hyorhinis. While colonization by M. flocculare is virtually asymptomatic, M. hyopneumoniae is the causative agent of enzootic pneumonia and M. hyorhinis is present in cases of pneumonia, polyserositis and arthritis. The genomic resemblance among these three Mycoplasma species combined with their different levels of pathogenicity is an indication that they have unknown mechanisms of virulence and differential expression, as for most mycoplasmas.

Methods

In this work, we performed whole-genome metabolic network reconstructions for these three mycoplasmas. Cultivation tests and metabolomic experiments through nuclear magnetic resonance spectroscopy (NMR) were also performed to acquire experimental data and further refine the models reconstructed in silico.

Results

Even though the refined models have similar metabolic capabilities, interesting differences include a wider range of carbohydrate uptake in M. hyorhinis, which in turn may also explain why this species is a widely contaminant in cell cultures. In addition, the myo-inositol catabolism is exclusive to M. hyopneumoniae and may be an important trait for virulence. However, the most important difference seems to be related to glycerol conversion to dihydroxyacetone-phosphate, which produces toxic hydrogen peroxide. This activity, missing only in M. flocculare, may be directly involved in cytotoxicity, as already described for two lung pathogenic mycoplasmas, namely Mycoplasma pneumoniae in human and Mycoplasma mycoides subsp. mycoides in ruminants. Metabolomic data suggest that even though these mycoplasmas are extremely similar in terms of genome and metabolism, distinct products and reaction rates may be the result of differential expression throughout the species.

Conclusions

We were able to infer from the reconstructed networks that the lack of pathogenicity of M. flocculare if compared to the highly pathogenic M. hyopneumoniae may be related to its incapacity to produce cytotoxic hydrogen peroxide. Moreover, the ability of M. hyorhinis to grow in diverse sites and even in different hosts may be a reflection of its enhanced and wider carbohydrate uptake. Altogether, the metabolic differences highlighted in silico and in vitro provide important insights to the different levels of pathogenicity observed in each of the studied species.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-016-2644-z) contains supplementary material, which is available to authorized users.

Characterization of triosephosphate isomerase from Mycoplasma gallisepticum

Triosephosphate isomerase (Tpi) is a glycolytic enzyme that is essential for efficient energy production in many pathogens. However, its function in Mycoplasma gallisepticum has not been fully elucidated. In this study, the mga0357 gene of M. gallisepticum, which encodes TpiA (MGTpiA), was amplified and expressed in Escherichia coli by IPTG induction. The purified recombinant MGTpiA protein exhibited catalytic activity that was similar to TPI from rabbit muscle, reducing NAD(+) to NADH. The MGTpiA was also found to be a surface-exposed protein by western blotting and immunofluorescence assays. In addition, cytadherence inhibition assays confirmed that the cytadherence of M. gallisepticum to the DF-1 cells was significantly inhibited by the anti-MGTpiA serum. The results of the study suggested that MGTpiA plays an important role in the metabolism and closely related to the M. gallisepticum pathogenicity.

Disruption of the membrane nuclease gene (MBOVPG45_0215) of Mycoplasma bovis greatly reduces cellular nuclease activity

Although the complete genome sequences of three strains of Mycoplasma bovis are available, few studies have examined gene function in this important pathogen. Mycoplasmas lack the biosynthetic machinery for the de novo synthesis of nucleic acid precursors, so nucleases are likely to be essential for them to acquire nucleotide precursors. Three putative membrane nucleases have been annotated in the genome of M. bovis strain PG45, MBOVPG45_0089 and MBOVPG45_0310, both of which have the thermonuclease (TNASE_3) functional domain, and MBOVPG45_0215 (mnuA), which has an exonuclease/endonuclease/phosphatase domain. While previous studies have demonstrated the function of TNASE_3 domain nucleases in several mycoplasmas, quantitative comparisons of the contributions of different nucleases to cellular nuclease activity have been lacking. Mapping of a library of 319 transposon mutants of M. bovis PG45 by direct genome sequencing identified mutants with insertions in MBOVPG45_0310 (the Δ0310 mutant) and MBOVPG45_0215 (the Δ0215 mutant). In this study, the detection of the product of MBOVPG45_0215 in the Triton X-114 fraction of M. bovis cell lysates, its cell surface exposure, and its predicted signal peptide suggested that it is a surface-exposed lipoprotein nuclease. Comparison of a ΔmnuA mutant with wild-type M. bovis on native and denatured DNA gels and in digestion assays using double-stranded phage λ DNA and closed circular plasmid DNA demonstrated that inactivation of this gene abolishes most of the cellular exonuclease and endonuclease activity of M. bovis. This activity could be fully restored by complementation with the wild-type mnuA gene, demonstrating that MnuA is the major cellular nuclease of M. bovis.

IMPORTANCE

Nucleases are thought to be important contributors to virulence and crucial for the maintenance of a nutritional supply of nucleotides in mycoplasmas that are pathogenic in animals. This study demonstrates for the first time that of the three annotated cell surface nuclease genes in an important pathogenic mycoplasma, the homologue of the thermostable nuclease identified in Gram-positive bacteria is responsible for the majority of the nuclease activity detectable in vitro.

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.

Flux balance analysis (FBA) is a powerful approach for genome-scale metabolic modeling. It provides metabolic engineers with a tool for manipulating, predicting, and optimizing metabolism for biotechnological and biomedical purposes. However, we posit that it can also be used as tool for fundamental research in understanding and curating metabolic networks. Specifically, by using a genetic algorithm integrated with FBA, we developed a curation approach to identify missing reactions, incomplete reactions, and erroneous reactions. Additionally, it was possible to take advantage of the ensemble information from the genetic algorithm to identify the most critical reactions for curation. We tested our strategy using Mycoplasma gallisepticum as our model organism. Using the genome annotation as the basis, the preliminary genome-scale metabolic model consisted of 446 metabolites involved in 380 reactions. Carrying out our analysis, we found over 80 incorrect reactions and 16 missing reactions. Based upon the guidance of the algorithm, we were able to curate and resolve all discrepancies. The model predicted an average bacterial growth rate of 0.358±0.12 h⁻¹ compared to the experimentally observed 0.244±0.03 h⁻¹. Thus, our approach facilitated the curation of a genome-scale metabolic network and generated a high quality metabolic model.

A whole-cell computational model predicts phenotype from genotype

Understanding how complex phenotypes arise from individual molecules and their interactions is a primary challenge in biology that computational approaches are poised to tackle. We report a whole-cell computational model of the life cycle of the human pathogen Mycoplasma genitalium that includes all of its molecular components and their interactions. An integrative approach to modeling that combines diverse mathematics enabled the simultaneous inclusion of fundamentally different cellular processes and experimental measurements. Our whole-cell model accounts for all annotated gene functions and was validated against a broad range of data. The model provides insights into many previously unobserved cellular behaviors, including in vivo rates of protein-DNA association and an inverse relationship between the durations of DNA replication initiation and replication. In addition, experimental analysis directed by model predictions identified previously undetected kinetic parameters and biological functions. We conclude that comprehensive whole-cell models can be used to facilitate biological discovery.

'Ca. Liberibacter asiaticus' proteins orthologous with pSymA-encoded proteins of Sinorhizobium meliloti: hypothetical roles in plant host interaction

Sinorhizobium meliloti strain 1021, a nitrogen-fixing, root-nodulating bacterial microsymbiont of alfalfa, has a 3.5 Mbp circular chromosome and two megaplasmids including 1.3 Mbp pSymA carrying nonessential 'accessory' genes for nitrogen fixation (nif), nodulation and host specificity (nod). A related bacterium, psyllid-vectored 'Ca. Liberibacter asiaticus,' is an obligate phytopathogen with a reduced genome that was previously analyzed for genes orthologous to genes on the S. meliloti circular chromosome. In general, proteins encoded by pSymA genes are more similar in sequence alignment to those encoded by S. meliloti chromosomal orthologs than to orthologous proteins encoded by genes carried on the 'Ca. Liberibacter asiaticus' genome. Only two 'Ca. Liberibacter asiaticus' proteins were identified as having orthologous proteins encoded on pSymA but not also encoded on the chromosome of S. meliloti. These two orthologous gene pairs encode a Na(+)/K+ antiporter (shared with intracellular pathogens of the family Bartonellacea) and a Co++, Zn++ and Cd++ cation efflux protein that is shared with the phytopathogen Agrobacterium. Another shared protein, a redox-regulated K+ efflux pump may regulate cytoplasmic pH and homeostasis. The pSymA and 'Ca. Liberibacter asiaticus' orthologs of the latter protein are more highly similar in amino acid alignment compared with the alignment of the pSymA-encoded protein with its S. meliloti chromosomal homolog. About 182 pSymA encoded proteins have sequence similarity (≤ E-10) with 'Ca. Liberibacter asiaticus' proteins, often present as multiple orthologs of single 'Ca. Liberibacter asiaticus' proteins. These proteins are involved with amino acid uptake, cell surface structure, chaperonins, electron transport, export of bioactive molecules, cellular homeostasis, regulation of gene expression, signal transduction and synthesis of amino acids and metabolic cofactors. The presence of multiple orthologs defies mutational analysis and is consistent with the hypothesis that these proteins may be of particular importance in host/microbe interaction and their duplication likely facilitates their ongoing evolution.

Comparison of the 'Ca. Liberibacter asiaticus' genome adapted for an intracellular lifestyle with other members of the Rhizobiales

An intracellular plant pathogen 'Candidatus Liberibacter asiaticus,' a member of the Rhizobiales, is related to Sinorhizobium meliloti, Bradyrhizobium japonicum, nitrogen fixing endosymbionts, Agrobacterium tumefaciens, a plant pathogen, and Bartonella henselae, an intracellular mammalian pathogen. Whole chromosome comparisons identified at least 50 clusters of conserved orthologous genes found on the chromosomes of all five metabolically diverse species. The intracellular pathogens 'Ca. Liberibacter asiaticus' and Bartonella henselae have genomes drastically reduced in gene content and size as well as a relatively low content of guanine and cytosine. Codon and amino acid preferences that emphasize low guanosine and cytosine usage are globally employed in these genomes, including within regions of microsynteny and within signature sequences of orthologous proteins. The length of orthologous proteins is generally conserved, but not their isoelectric points, consistent with extensive amino acid substitutions to accommodate selection for low GC content. The 'Ca. Liberibacter asiaticus' genome apparently has all of the genes required for DNA replication present in Sinorhizobium meliloti except it has only two, rather than three RNaseH genes. The gene set required for DNA repair has only one rather than ten DNA ligases found in Sinorhizobium meliloti, and the DNA PolI of 'Ca. Liberibacter asiaticus' lacks domains needed for excision repair. Thus the ability of 'Ca. Liberibacter asiaticus' to repair mutations in its genome may be impaired. Both 'Ca. Liberibacter asiaticus and Bartonella henselae lack enzymes needed for the metabolism of purines and pyrimidines, which must therefore be obtained from the host. The 'Ca. Liberibacter asiaticus' genome also has a greatly reduced set of sigma factors used to control transcription, and lacks sigma factors 24, 28 and 38. The 'Ca. Liberibacter asiaticus' genome has all of the hallmarks of a reduced genome of a pathogen adapted to an intracellular lifestyle.

Cloning and expression of mitochondrial malate dehydrogenase of Clonorchis sinensis

The NAD-dependent mitochondrial malate dehydrogenase (mMDH, EC1.1.1.37) plays pivotal roles in tricarboxylic acid and is crucial for the survival and pathogenecity of parasites. A cDNA, which was identified by high throughput sequencing from the cDNA library constructed from adult Clonorchis sinensis, encoded a putative peptide of 341 amino acids with more than 50% identity with mMDHs from other organisms. The mMDH was expressed in Escherichia coli as the recombinant protein with a GST tag and purified by glutathione-Sepharose 4B column. The recombinant mMDH showed MDH activity of 63.6 U/mg, without lactate dehydrogenase activity and NADPH selectivity. The kinetic constants of recombinant mMDH were determined.

Regulation of Carbon Metabolism in the Mollicutes and Its Relation to Virulence

The mollicutes are cell wall-less bacteria that live in close association with their eukaryotic hosts. Their genomes are strongly reduced and so are their metabolic capabilities. A survey of the available genome sequences reveals that the mollicutes are capable of utilizing sugars as source of carbon and energy via glycolysis. The pentose phosphate pathway is incomplete in these bacteria, and genes encoding enzymes of the tricarboxylic acid cycle are absent from the genomes. Sugars are transported by the phosphotransferase system. As in related bacteria, the phosphotransferase system does also seem to play a regulatory role in the mollicutes as can be concluded from the functionality of the regulatory HPr kinase/phosphorylase. In Mycoplasma pneumoniae, the activity of HPr kinase is triggered in the presence of glycerol. This carbon source may be important for the mollicutes since it is available in epithelial tissues and its metabolism results in the formation of hydrogen peroxide, the major virulence factor of several mollicutes. In plant-pathogenic mollicutes such as Spiroplasma citri, the regulation of carbon metabolism is crucial in the adaptation to life in plant tissues or the insect vectors. Thus, carbon metabolism seems to be intimately linked to pathogenicity in the mollicutes.

Enzymatic and physico-chemical characteristics of recombinant cMDH and mMDH of Clonorchis sinensis

The cytosol and mitochondrial malate dehydrogenases (MDHs, EC 1.1.1.37) of Clonorchis sinensis were expressed in Escherichia coli as a fusion protein with a 6xHis and GST tag, respectively. The cytosol MDH of Clonorchis sinensis (Cs-cMDH) has higher resistibility to acid than mitochondrial MDH (Cs-mMDH). The Cs-cMDH also has higher heat resistibility and thermal stability than Cs-mMDH. Although there is only 22.8% identity between the amino acid sequences of Cs-cMDH and Cs-mMDH, they share several conserved residues. There are some differences between the circular dichroism spectra of Cs-cMDH and Cs-mMDH, but they have approximate percentages of helix. 4,4'-Bisdimethylamino diphenylcarbinol can decrease the Cs-mMDH activity but not the Cs-cMDH activity. Paraziquantel, metronidazole and albendazole did not inhibit the enzymes' activity, but adenosine 5'-monophosphate showed competitive inhibition to enzyme, with the Ki for Cs-cMDH and Cs-mMDH being 2.81 and 0.49 mM, respectively.

Essential genes of a minimal bacterium

Mycoplasma genitalium has the smallest genome of any organism that can be grown in pure culture. It has a minimal metabolism and little genomic redundancy. Consequently, its genome is expected to be a close approximation to the minimal set of genes needed to sustain bacterial life. Using global transposon mutagenesis, we isolated and characterized gene disruption mutants for 100 different nonessential protein-coding genes. None of the 43 RNA-coding genes were disrupted. Herein, we identify 382 of the 482 M. genitalium protein-coding genes as essential, plus five sets of disrupted genes that encode proteins with potentially redundant essential functions, such as phosphate transport. Genes encoding proteins of unknown function constitute 28% of the essential protein-coding genes set. Disruption of some genes accelerated M. genitalium growth.

Suspected utility of enzymes with multiple activities in the small genome Mycoplasma species: the replacement of the missing "household" nucleoside diphosphate kinase gene and activity by glycolytic kinases

The small genome Mollicutes whose DNAs are completely sequenced (Mycoplasma genitalium, Mycoplasma pneumoniae, Mycoplasma pulmonis, and Ureaplasma urealyticum [parvum]) lack a gene (ndk) for the presumably essential nucleoside diphosphate kinase (NDPK). We hypothesized that other activities might replace NDPK activity. We found in M. genitalium G37(T), Mycoplasma pneumoniae FH(T), Mycoplasma fermentans PG18(T), and Mycoplasma capricolum subsp. capricolum Kid(T) that their 6-phosphofructokinases (6-PFKs), phosphoglycerate kinases (PGKs), pyruvate kinases (PKs), and acetate kinases (AKs), besides reactant ADP/ATP, could use other ribo- and deoxyribo-purine and pyrimidine NDPs and NTPs. These activities could compensate for the absence of an orthologous ndk gene in the Mycoplasmataceae. They suggest a metabolically varied and consequential role for unrelated and perhaps unsuspected "replacement" or compensatory enzymes that may confound metabolic prediction. We partially purified and biochemically characterized the PKs, 6-PFKs, PGKs, and AKs from M. capricolum subsp. capricolum Kid(T) and M. fermentans PG18(T).

Probing the specificity of a trypanosomal aromatic alpha-hydroxy acid dehydrogenase by site-directed mutagenesis

The aromatic l-alpha-hydroxy acid dehydrogenase (AHDAH) from Trypanosoma cruzi has over 50% sequence identity with cytosolic malate dehydrogenases (cMDHs), yet it is unable to reduce oxaloacetate. Molecular modeling of the three-dimensional structure of AHADH using the pig cMDH as template directed the construction of several mutants. AHADH shares with MDHs the essential catalytic residues H195 and R171 (using Eventoff's numbering). The AHADH A102R mutant became able to reduce oxaloacetate, while remaining fully active towards aromatic alpha-oxoacids. The Y237G mutant diminished its affinity for all of the natural substrates, whereas the double mutant A102R/Y237G was more active than Y237G and had similar activity with oxaloacetate and with aromatic substrates. The present results reinforce our proposal that AHADH arose by a moderate number of point mutations from a cMDH no longer present in the parasite.

Ureaplasma urealyticum: an opportunity for combinatorial genomics

Examination of genomic or enzymatic activity data alone neither provides a complete picture of metabolic function or potential nor confidently reveals sites amenable to inhibition. Furthermore, in some cases, gene annotation and in aqua assays disagree by describing gene annotation without enzyme activity and enzyme activity without homologous annotation. The newly sequenced genome of Ureaplasma urealyticum (parvum) is another prokaryote example of the class Mollicutes where such confounding differences are observed. The little-considered role of some proteins as multifunctional enzymes - substitutes for 'missing' genes - could both partially explain the apparent anomalies and relate to any inaccurate deductions of inhibitor function. A combinatorial analysis involving available evidence of genomic sequence, transcription, translational phenomena, structure and enzymatic activity gives the best picture of the organism's vital metabolic alternatives.

The proteome of Mycoplasma genitalium. Chaps-soluble component

Mycoplasma genitalium is the smallest member of the class Mollicutes, with a genome size of 580 kb. It has the potential to express 480 gene products, and is therefore considered to be an excellent model to assess: (a) the minimum metabolism required by a free living cell; and (b) proteomic technologies and the information obtained by proteome analysis. Here, we report on the most complete proteome observed at 73% (expected proteome), and analysed at 33% (reported proteome). The use of four overlapping pH windows in conjunction with SDS/PAGE has allowed 427 distinct proteins to be resolved in association with the exponential growth of M. genitalium. Proof of expression for 201 proteins of sufficient abundance on silver stained two-dimensional gels was obtained using peptide mass fingerprinting (PMF) of which 158 were identified. The potential for gene product modification in even the simplest known self-replicating organism was quantified at a ratio of 1.22 : 1, more proteins than genes. A reduction in protein expression of 42% was observed for post-exponentially-grown cells. DnaK, GroEL, DNA gyrase, and a cytadherence accessory protein were significantly elevated, while some ribosomal proteins were reduced in relative abundance. The strengths and weaknesses of techniques employed were assessed with respect to the observed and predicted proteome derived from DNA sequence information. Proteomics was shown to provide a perspective into the biochemical and metabolic activities of this organism, beyond that obtainable by sequencing alone.

Convergent evolution of Trichomonas vaginalis lactate dehydrogenase from malate dehydrogenase

Lactate dehydrogenase (LDH) is present in the amitochondriate parasitic protist Trichomonas vaginalis and some but not all other trichomonad species. The derived amino acid sequence of T. vaginalis LDH (TvLDH) was found to be more closely related to the cytosolic malate dehydrogenase (MDH) of the same species than to any other LDH. A key difference between the two T. vaginalis sequences was that Arg91 of MDH, known to be important in coordinating the C-4 carboxyl of oxalacetate/malate, was replaced by Leu91 in LDH. The change Leu91Arg by site-directed mutagenesis converted TvLDH into an MDH. The reverse single amino acid change Arg91Leu in TvMDH, however, gave a product with no measurable LDH activity. Phylogenetic reconstructions indicate that TvLDH arose from an MDH relatively recently.

Molecular biology and pathogenicity of mycoplasmas

The recent sequencing of the entire genomes of Mycoplasma genitalium and M. pneumoniae has attracted considerable attention to the molecular biology of mycoplasmas, the smallest self-replicating organisms. It appears that we are now much closer to the goal of defining, in molecular terms, the entire machinery of a self-replicating cell. Comparative genomics based on comparison of the genomic makeup of mycoplasmal genomes with those of other bacteria, has opened new ways of looking at the evolutionary history of the mycoplasmas. There is now solid genetic support for the hypothesis that mycoplasmas have evolved as a branch of gram-positive bacteria by a process of reductive evolution. During this process, the mycoplasmas lost considerable portions of their ancestors' chromosomes but retained the genes essential for life. Thus, the mycoplasmal genomes carry a high percentage of conserved genes, greatly facilitating gene annotation. The significant genome compaction that occurred in mycoplasmas was made possible by adopting a parasitic mode of life. The supply of nutrients from their hosts apparently enabled mycoplasmas to lose, during evolution, the genes for many assimilative processes. During their evolution and adaptation to a parasitic mode of life, the mycoplasmas have developed various genetic systems providing a highly plastic set of variable surface proteins to evade the host immune system. The uniqueness of the mycoplasmal systems is manifested by the presence of highly mutable modules combined with an ability to expand the antigenic repertoire by generating structural alternatives, all compressed into limited genomic sequences. In the absence of a cell wall and a periplasmic space, the majority of surface variable antigens in mycoplasmas are lipoproteins. Apart from providing specific antimycoplasmal defense, the host immune system is also involved in the development of pathogenic lesions and exacerbation of mycoplasma induced diseases. Mycoplasmas are able to stimulate as well as suppress lymphocytes in a nonspecific, polyclonal manner, both in vitro and in vivo. As well as to affecting various subsets of lymphocytes, mycoplasmas and mycoplasma-derived cell components modulate the activities of monocytes/macrophages and NK cells and trigger the production of a wide variety of up-regulating and down-regulating cytokines and chemokines. Mycoplasma-mediated secretion of proinflammatory cytokines, such as tumor necrosis factor alpha, interleukin-1 (IL-1), and IL-6, by macrophages and of up-regulating cytokines by mitogenically stimulated lymphocytes plays a major role in mycoplasma-induced immune system modulation and inflammatory responses.

Purification, cloning, and preliminary characterization of a Spiroplasma citri ribosomal protein with DNA binding capacity

The rpsB-tsf-x operon of Spiroplasma citri encodes ribosomal protein S2 and elongation factor Ts, two components of the translational apparatus, and an unidentified X protein. A potential DNA-binding site (a 20-base pair (bp) inverted repeat sequence) is located at the 3' end of rpsB. Southwestern analysis of S. citri proteins, with a 30-bp double-stranded oligonucleotide probe (IRS), containing the 20-bp inverted repeat sequence and the genomic flanking sequences, detected an IRS-binding protein of 46 kDa (P46). P46 protein, which displays preferential affinity for the IRS, was purified from S. citri by a combination of affinity and gel filtration chromatographies. The native form of P46 seems to be homomultimeric as estimated by SDS-polyacrylamide gel electrophoresis analysis and gel filtration. A 3.5-kilobase pair S. citri DNA fragment comprising the P46 gene and flanking sequences was cloned and sequenced. Sequence analysis of this DNA fragment indicated that the P46 gene is located within the S10-spc operon of S. citri at the position of the gene coding for ribosomal protein L29 in the known S10-spc operons. The similarity between the N-terminal domain of P46 and the L29 ribosomal protein family and the presence of a 46-kDa IRS-binding protein in S. citri ribosomes indicated that P46 is the L29 ribosomal protein of S. citri. We suggest that P46 is a bifunctional protein with an L29 N-terminal domain and a C-terminal domain involved in IRS binding.

The comparative metabolism of the mollicutes (Mycoplasmas): the utility for taxonomic classification and the relationship of putative gene annotation and phylogeny to enzymatic function in the smallest free-living cells

Mollicutes or mycoplasmas are a class of wall-less bacteria descended from low G + C% Gram-positive bacteria. Some are exceedingly small, about 0.2 micron in diameter, and are examples of the smallest free-living cells known. Their genomes are equally small; the smallest in Mycoplasma genitalium is sequenced and is 0.58 mb with 475 ORFs, compared with 4.639 mb and 4288 ORFs for Escherichia coli. Because of their size and apparently limited metabolic potential, Mollicutes are models for describing the minimal metabolism necessary to sustain independent life. Mollicutes have no cytochromes or the TCA cycle except for malate dehydrogenase activity. Some uniquely require cholesterol for growth, some require urea and some are anaerobic. They fix CO2 in anaplerotic or replenishing reactions. Some require pyrophosphate not ATP as an energy source for reactions, including the rate-limiting step of glycolysis: 6-phosphofructokinase. They scavenge for nucleic acid precursors and apparently do not synthesize pyrimidines or purines de novo. Some genera uniquely lack dUTPase activity and some species also lack uracil-DNA glycosylase. The absence of the latter two reactions that limit the incorporation of uracil or remove it from DNA may be related to the marked mutability of the Mollicutes and their tachytelic or rapid evolution. Approximately 150 cytoplasmic activities have been identified in these organisms, 225 to 250 are presumed to be present. About 100 of the core reactions are graphically linked in a metabolic map, including glycolysis, pentose phosphate pathway, arginine dihydrolase pathway, transamination, and purine, pyrimidine, and lipid metabolism. Reaction sequences or loci of particular importance are also described: phosphofructokinases, NADH oxidase, thioredoxin complex, deoxyribose-5-phosphate aldolase, and lactate, malate, and glutamate dehydrogenases. Enzymatic activities of the Mollicutes are grouped according to metabolic similarities that are taxonomically discriminating. The arrangements attempt to follow phylogenetic relationships. The relationships of putative gene assignments and enzymatic function in My. genitalium, My. pneumoniae, and My. capricolum subsp. capricolum are specially analyzed. The data are arranged in four tables. One associates gene annotations with congruent reports of the enzymatic activity in these same Mollicutes, and hence confirms the annotations. Another associates putative annotations with reports of the enzyme activity but from different Mollicutes. A third identifies the discrepancies represented by those enzymatic activities found in Mollicutes with sequenced genomes but without any similarly annotated ORF. This suggests that the gene sequence is significantly different from those already deposited in the databanks and putatively annotated with the same function. Another comparison lists those enzymatic activities that are both undetected in Mollicutes and not associated with any ORF. Evidence is presented supporting the theory that there are relatively small gene sequences that code for functional centers of multiple enzymatic activity. This property is seemingly advantageous for an organism with a small genome and perhaps under some coding restraint. The data suggest that a concept of "remnant" or "useless genes" or "useless enzymes" should be considered when examining the relationship of gene annotation and enzymatic function. It also suggests that genes in addition to representing what cells are doing or what they may do, may also identify what they once might have done and may never do again.

Proteome analysis: genomics via the output rather than the input code

A knowledge of the 'proteome,' total protein output encoded by a genome, provides information on (1) if and when predicted gene products are translated, (2) the relative concentrations of gene products, and (3) the extent of posttranslational modification, none of which can be accurately predicted from the nucleic acid sequence alone. The current status of proteome analysis is reviewed with respect to some of the techniques employed, automation, relevance to genomic studies, mass spectrometry and bioinformatics, limitations, and recent improvements in resolution and sensitivity for the detection of protein expression in whole cells, tissues, or organisms. The concept of 'proteomic contigs' is introduced for the first time. Traditional approaches to genomic analysis call upon a number of strategies to produce contiguous DNA sequence information, while 'proteomic contigs' are derived from multiple molecular mass and isoelectric point windows in order to construct a picture of the total protein expression within living cells. In higher eukaryotes, the latter may require several dozen image subsets of protein spots to be stitched together using advanced image analysis. The utility of both experimental and theoretical peptide-mass fingerprinting (PMF) and associated bioinformatics is outlined. A previously unknown motif within the peptide sequence of Elongation Factor Tu from Thermus aquaticus was discovered using PMF. This motif was shown to possess potential significance in maintaining structural integrity of the entire molecule.