Previous Article | Next Article ![]()
Journal of Bacteriology, September 2006, p. 6135-6142, Vol. 188, No. 17
0021-9193/06/$08.00+0 doi:10.1128/JB.00207-06
Copyright © 2006, American Society for Microbiology. All Rights Reserved.
Fundación Instituto Leloir, Buenos Aires, Argentina,1 Lehrstuhl für Organische Chemie und Biochemie Technische Universität München, Lichtenbergstraße 4, D-85748 Garching, Germany2
Received 7 February 2006/ Accepted 7 June 2006
|
|
|---|
|
|
|---|
![]() View larger version (13K): [in a new window] |
FIG. 1. Biosynthesis of riboflavin in eubacteria. I, GTP cyclohydrolase II; II, 2,5-diamino-6-ribosylamino-4(3H)-pyrimidinone 5'-phosphate deaminase; III, 5-amino-6-ribosylamino-2,4(1H,3H)-pyrimidinedione 5'-phosphate reductase; IV, hypothetical phosphatase; V, 3,4-dihydroxy-2-butanone 4-phosphate synthase; VI, 6,7-dimethyl-8-ribityllumazine synthase; VII, riboflavin synthase; 1, GTP; 2, 5-amino-6-ribitylamino-2,4(1H,3H)-pyrimidinedione; 3, ribulose 5'-phosphate; 4, 3,4-dihydroxy-2-butanone 4-phosphate; 5, 6,7-dimethyl-8-ribityllumazine; 6, riboflavin.
|
Lumazine synthases from a variety of eubacteria (including Escherichia coli, Bacillus subtilis, Mycobacterium tuberculosis, and the hyperthermophile Aquifex aeolicus), archaea (Methanococcus jannaschii), fungi (Saccharomyces cerevisiae, Schizosaccharomyces pombe, and Magnaporthe grisea), and a plant (spinach) have been studied in some detail (11, 16, 22, 32, 34, 35, 40-42, 44, 56, 57). The enzymes from fungi and from M. tuberculosis are C5-symmetric homopentamers, and the lumazine synthases of plants, most eubacteria, and archaea are 532-symmetric, hollow capsids, which are best described as dodecamers of pentamers. The subunit folds of these enzymes and the topology of the pentamer moieties are closely similar. The topologically equivalent active sites (5 in the case of the pentameric enzymes and 60 in the case of the icosahedral enzymes) are invariably located at interfaces between adjacent subunits in the pentamer moieties. Recently, it was found that pentameric riboflavin synthases of archaea are closely related to 6,7-dimethyl-8-ribityllumazine synthases (15, 25, 45).
Brucellosis is a disease of humans and livestock that is caused by closely related Brucella species adapted to intracellular life within the cells of a variety of mammals; the main pathogenic species for domestic animals are Brucella abortus, Brucella melitensis, and Brucella suis. Goldbaum and coworkers have shown that an 18-kDa B. abortus antigen with sequence similarity to lumazine synthases is a serological marker of active disease in human brucellosis patients (17, 19). Immunization with the protein has been shown to induce both cellular and humoral immune responses in mice. Moreover, the generation of protective immunity has also been observed in this model (2). Hence, the 18-kDa protein is of considerable immunologic interest and has been suggested to be a general carrier for the engineering of subunit vaccines (36). In preliminary enzymatic studies, this protein was shown to catalyze the formation of 6,7-dimethyl-8-ribityllumazine (compound 5), albeit at a low rate. Recently, the catalytic properties of this protein have been analyzed in closer detail (31). The three-dimensional structure of the 18-kDa antigen has been studied in considerable detail by Goldbaum and coworkers, and the protein has been reported to be a homodecamer which is best described as a D5-symmetric dimer of pentamers (31, 58). The 18-kDa antigen of B. abortus is the first lumazine synthase that has been reported to exhibit this quaternary arrangement.
Brucella spp. are known to share a genome topology characterized by two circular chromosomes with approximate sizes of 2.1 MDa (chromosome I) and 1.2 MDa (chromosome II). The recently published B. abortus genome sequence indicates that the 18-kDa antigen is specified by the ribH2 gene, located on chromosome II of B. abortus. Moreover, a second gene (designated ribH1) with sequence similarity to lumazine synthase was found on chromosome I.
This paper describes the identification, cloning, and expression of the ribH1 gene of B. abortus and provides a biochemical characterization of the encoded type I lumazine synthase. Based on this study, we suggest that the decameric protein encoded by the ribH2 gene be designated type II lumazine synthase. We also describe the genomic organization of the ribH1 and ribH2 genes and provide a phylogenetic analysis of lumazine synthases and related pentameric riboflavin synthases derived from different organisms.
|
|
|---|
Cloning and expression of the ribH1 gene of B. abortus. The ribH1 gene was amplified by PCR using B. abortus chromosomal DNA as the template and the oligonucleotides RibH1-Ndel-Vo 5'-ATAATAATACATATGGAGTTTCTCATGTCCAAGCAC-3' and RibH1-HindIII-Hi 5'-TATTATTATAAGCTTAGGCTCCGAATTTTTTGCGCAGGC-3' as primers. The amplicon was digested with NdeI and HindIII. The resulting fragment was purified with an agarose gel extraction kit (QIAGEN, Hilden, Germany) and was then ligated to the plasmid pT7-7 (51), which had been treated with the same restriction enzymes. The ligation mixture was transformed into E. coli strain XL1-Blue (3), resulting in the recombinant strain XL-1-pT7-7-BARibH1. The reisolated plasmid was sequenced by the method of Sanger (47) (GATC Biotech, Konstanz, Germany) and was transformed into E. coli BL21(DE3) competent cells (50) (Stratagene, La Jolla, CA), resulting in the recombinant E. coli strain BL21(DE3)-pT7-7-BARibH1.
Bacterial culture. E. coli BL21(DE3)-pT7-7-BARibH1 was grown to an optical density (at 600 nm) of 1.0 in LB medium containing 100 µg of ampicillin per ml at 37°C with shaking (150 rpm). An aliquot (5 ml) of this culture was diluted into 500 ml of medium, and incubation with shaking was continued to an optical density of 1.0. Isopropyl-ß-thiogalactoside was added to a final concentration of 1 mM, and the suspension was incubated for 4 h at 37°C with shaking (150 rpm).
Purification of type I lumazine synthase of B. abortus. Frozen bacterial cell mass was thawed in 50 mM potassium phosphate, pH 7.0. The suspension was ultrasonically treated and centrifuged. The supernatant was loaded on top of a Q-Sepharose column (1.8 by 25 cm) that had been equilibrated with 50 mM potassium phosphate, pH 7.0. The column was developed with a linear gradient of 0.05 to 1.0 M potassium phosphate, pH 7.0. Fractions were combined and concentrated by ultrafiltration. The solution was loaded on top of a Superdex-200 column (2.6 by 60 cm), which was developed with 100 mM potassium phosphate, pH 7.0. Fractions were combined and concentrated by ultrafiltration.
Assay of lumazine synthase activity. Steady-state kinetic experiments were performed as described previously (12). Corrected initial rates were fitted to the Hill equation (equation 1) or equation 2 (see Results), which represents a classical substrate inhibition model, with the program package Origin or the program DynaFit, version 3.28.024 (33).
Analytical ultracentrifugation. Experiments were performed with an analytical ultracentrifuge, Optima XL-A, from Beckman Instruments (Palo Alto, CA), equipped with absorbance and interference optics. Aluminum double-sector cells equipped with quartz windows were used throughout. The protein concentration was monitored photometrically at 280 nm. For boundary sedimentation experiments, a solution containing 100 mM potassium phosphate, pH 5.0, 300 mM sodium chloride, and 1.2 mg of protein per ml was centrifuged at 59,000 rpm and 20°C. Sedimentation equilibrium experiments were performed with a solution containing 100 mM potassium phosphate, pH 5.0, 300 mM sodium chloride, and 0.4 mg of protein per ml. Samples were centrifuged at 10,000 rpm and 4°C for 72 h. The partial specific volume was estimated from the amino acid composition, resulting in a value of 0.733 ml g1 (37).
Phylogenetic analysis. Sequences used for the phylogenetic analysis of lumazine synthase proteins were obtained from the HMM library, the genome assignment server Superfamily 1.69, and the NCBI database. The amino acid sequences were aligned with the ClustalX (version 1.81) program for multiple sequence alignment (53). Phylogenetic analysis of this alignment was inferred with the maximum-likelihood heuristic algorithm implemented by PHYML, version 2.4.4 (21), under the JTT substitution model (26). The reliability of tree nodes was analyzed by generating 1,000 bootstrap trees (8).
|
|
|---|
![]() View larger version (43K): [in a new window] |
FIG. 2. Topology of riboflavin biosynthesis operons. Sequences were derived from fully sequenced genomes available from GenBank. Orientations of the genes involved in riboflavin biosynthesis or adjacent to those genes are given as arrows (not drawn to scale). Organisms (accession numbers) are as follows: Brucella meliloti 16 M chromosome I (NC 003317), Brucella suis 1330 chromosome I (NC 004310), Brucella abortus biovar 1 strain 9-941 chromosome I (NC 006932), Sinorhizobium meliloti 1021 (NC_003047), Mesorhizobium loti MAFF303099 (NC 002678), Bartonella henselae strain Houston-1 (NC 005956), Rhodopseudomonas palustris CGA009 (NC 005296), Bradyrhizobium japonicum USDA 110 (NC 004463), Caulobacter crescentus CB15 (NC_002696), Xanthomonas axonopodis pv. citri strain 306 (NC 003919), Xanthomonas oryzae KACC10331 (NC 006834), Xylella fastidiosa Temecula1 (NC_004556), Pseudomonas syringae pv. tomato strain DC3000 (NC_004578), Pseudomonas putida KT2440 (NC 002947), Methylococcus capsulatus strain Bath (NC 002977), Acinetobacter sp. ADP1 (NC 005966), Escherichia coli K-12 (NC_000913), Shigella sonnei Ss046 (NC 007384), Shigella flexneri 2a strain 301 (NC_004337), Salmonella enterica subsp. enterica serovar Choleraesuis strain SC-B67 (NC_006905), Yersinia pestis KIM (NC_004088), Photorhabdus luminescens subsp. laumondii TTO1 (NC 005126), Haemophilus influenzae 86-028NP (NC_007146), Geobacter sulfurreducens PCA (NC 002939), Aquifex aeolicus VF5 (NC 000918), Helicobacter pylori J99 (NC_000921), Candidatus (Blochmannia) floridanus (NC_005061), Buchnera aphidicola strain Sg (Schizaphis graminum) (NC_004061), Caulobacter crescentus CB15 (NC_002696), Ralstonia solanacearum GMI1000 (NC 003295), Chromobacterium violaceum ATCC 12472 (NC 005085), Azoarcus sp. EbN1 (NC 006513), Nitrosomonas europaea ATCC 19718 (NC 004757), Bordetella bronchiseptica RB50 (NC 002927), B. meliloti 16 M chromosome II (NC 003318), Brucella suis 1330 chromosome II (NC 004311), B. abortus biovar 1 strain 9-941 chromosome II (NC 006933), M. loti MAFF303099 (NC 002678), S. meliloti 1021 (NC_003047), R. palustris CGA009 (NC 005296), B. japonicum USDA 110 (NC 004463), P. syringae pv. tomato strain DC3000 (NC_004578).
|
Figure 2 shows the genomic localizations of ribH genes from a variety of organisms. Besides Brucella, other organisms (Mesorhizobium loti, Sinorhizobium meliloti, Rhodopseudomonas palustris, Bradyrhizobium japonicum, and Pseudomonas syringae) also have two different ORFs with sequence similarity to lumazine synthase (designated ribH1 and ribH2). This analysis could not detect any organism presenting only the ribH2 gene; moreover, the ribH2 gene was never found within an operon, with the exception of in Caulobacter crescentus. In the C. crescentus genome, we found a rib operon with a lumazine synthase resembling the ribH2 gene, whereas a gene with similarity to the ribH1 genes is located at a different site on the same chromosome; this putative ribH1 gene is probably not part of a larger transcription unit. Nine organisms, including Brucella, have shortened rib operons without ribA genes (bifunctional GTP cyclohydrolase II/3,4-dihydroxy-2-butanone 4-phosphate synthase); isolated ribA genes are located elsewhere on chromosome I. Interestingly, Pseudomonas syringae harbors a ribA gene within a rib operon and another gene coding for an additional monofunctional GTP cyclohydrolase II. In Sinorhizobium meliloti, the genes ribD/ribB and ribH/nusB are separated by three ORFs with unknown functions and opposite directions spanning 2,555 bp. Buchnera aphidicola is the only organism in this study that has separate genes for deaminase and reductase (ribD1 and ribD2), though arranged in tandem.
In many cases, the genes in the respective operons overlap by a few base pairs. Notably, in almost all operons (32 of 33 analyzed genomes), the nusB gene appears downstream of the rib genes, in most cases followed by thiL, which specifies thiamine monophosphate kinase (Fig. 2). The rib genes are invariably oriented in the same direction as nusB and thiL, although the functional relevance of this association is still unknown.
The amino acid sequence identity between the type I and type II lumazine synthases of B. abortus is 21% (Fig. 3). Interestingly, both lumazine synthases have higher sequence similarities with other lumazine synthases, e.g., B. subtilis (type I, 33%; type II, 24%) and E. coli (type I, 37%; type II, 29%) than with each other.
![]() View larger version (62K): [in a new window] |
FIG. 3. Sequence alignment of lumazine synthases from eubacteria and yeast and an archaeal riboflavin synthase. BabRibH2, type II lumazine synthase of B. abortus (P61711); BabRibH1, type I lumazine synthase of B. abortus (Q57DY1); BsuLS, lumazine synthase of Bacillus subtilis (P11998); AaeLS, lumazine synthase of Aquifex aeolicus (O66529); SpoLS, lumazine synthase of Schizosaccharomyces pombe (Q9UUB1); MjaRS, riboflavin synthase of Methanocaldococcus jannaschii (Q58584). Identical residues are printed in white on black, and similar residues are shaded in gray.
|
The type I lumazine synthase of B. abortus sediments at an apparent velocity of 5.9S at 20°C (Fig. 4A). In comparison, it should be noted that pentameric lumazine synthases from the yeasts S. cerevisiae and S. pombe have similar apparent sedimentation coefficients, of 5.5S and 5.0S, respectively (11, 42). Sedimentation equilibrium experiments indicated a molecular mass of 88 kDa, with an ideal monodisperse model used for calculation (Fig. 4B). The calculated subunit molecular mass of 17,599 Da indicates a pentameric mass of 88 kDa, in excellent agreement with the experimental data.
![]() View larger version (26K): [in a new window] |
FIG. 4. Analytical ultracentrifugation analysis of the type I lumazine synthase of B. abortus. (A) For boundary sedimentation centrifugation analysis, a solution containing 100 mM potassium phosphate, pH 5.0, 300 mM sodium chloride, and 1.2 mg of protein per ml was centrifuged at 59,000 rpm and 20°C for 3 h. The protein concentration was monitored photometrically (280 nm) at intervals of 5 min. (B) For sedimentation equilibrium centrifugation analysis, a solution containing 100 mM potassium phosphate, pH 5.0, 300 mM sodium chloride, and 0.4 mg of protein per ml was centrifuged at 10,000 rpm and 4°C for 72 h.
|
![]() | (1) |
|
View this table: [in a new window] |
TABLE 1. Properties of lumazine synthases
|
![]() | (2) |
As shown in Table 1, the affinity for 3,4-dihydroxy-2-butanone 4-phosphate (substrate 4) is in the same range as the values reported for other lumazine synthases. In contrast, the Km value for substrate 2 is about 10 times higher than that observed for other orthologs.
Table 2 summarizes the genomic localization and the structural and biochemical characterization of B. abortus type I and type II lumazine synthases.
|
View this table: [in a new window] |
TABLE 2. Summary of genomic, structural, and biochemical knowledge about B. abortus lumazine synthases
|
![]() View larger version (41K): [in a new window] |
FIG. 5. Phylogenetic analysis of the quaternary arrangements of lumazine synthases or related pentameric riboflavin synthases derived from different organisms, showing a distribution into four clearly defined branches: (i) type II lumazine synthases (decameric), (ii) archaeal pentameric riboflavin synthases, (iii) archaeal lumazine synthases, and (iv) type I eubacterial, fungal, and plant lumazine synthases. The corresponding organisms are indicated as follows (references describing the structures of the enzymes involved in this study are given in brackets): Aae, A. aeolicus (lumazine synthase [56, 57]); Afu, Archaeoglobus fulgidus; Ape, Aeropyrum pernix; Bab, B. abortus; Basu, B. subtilis (lumazine synthase [34, 35]); Bja, Bradyrhizobium japonicum; Bme, B. melitensis; Bsu, B. suis (RibH2 protein [58]); Ccr, Caulobacter crescentus; Eco, E. coli (lumazine synthase [42]); Fac, Ferroplasma acidarmanus; Mac, Methanosarcina acetivorans; Mgr, M. grisea (lumazine synthase [44]); Mja, Methanococcus jannaschii (riboflavin synthase [15, 45], lumazine synthase [22]); Mka, Methanopyrus kandleri; Mlo, M. loti; Mma, Methanosarcina mazei; Mth, Methanobacterium thermoautotrophicum; Mtu, M. tuberculosis (lumazine synthase [41]); Nfa, Nocardia farcinica; Pae, Pyrobaculum aerophilum; Pfl, Pseudomonas fluorescens; Pfu, Pyrococcus furiosus; Psy, P. syringae; Rpa, R. palustris; Sce, S. cerevisiae (lumazine synthase [40]); Sipo, Silicibacter pomeroyi; Sme, S. meliloti; Sol, S. oleracea (lumazine synthase [44]); Spo, S. pombe (lumazine synthase [16, 32]); Sso, Sulfolobus solfataricus; Sto, Sulfolobus tokodaii.
|
-Proteobacteria (Brucella, Rhizobium, and Rhodobacter) diverge from pentameric fungal and yeast lumazine synthases (S. cerevisiae, S. pombe, and M. grisea). |
|
|---|
-proteobacterial speciation by lateral gene transfer. At this time, there is no experimental evidence that would permit a decision. Recent studies have shown the existence of a family of pentameric lumazine synthase-like riboflavin synthases (Fig. 5 and 6). These enzymes have been found exclusively in archaea, which have also been shown to be devoid of riboflavin synthases of the trimeric eubacterium/yeast/plant type. The cavity harboring the active site of the pentameric riboflavin synthase of M. jannaschii is similar to that of the lumazine synthases. The binding mode of the acceptor lumazine molecule in pentameric riboflavin synthase closely resembles the pyrimidinedione-binding site of the lumazine synthases (Fig. 6). Notably, these archaeal riboflavin synthases are evolutionarily old and have no detectable lumazine synthase activity (13, 15, 45). Completely sequenced archaeal genomes typically comprise sets of two similar genes, coding for a lumazine synthase-like riboflavin synthase and a regular lumazine synthase. Interestingly, present-day lumazine synthases have retained the capacity to bind riboflavin. For example, the lumazine synthase of the yeast S. pombe is yellow colored due to the presence of a tightly bound riboflavin (11).
![]() View larger version (39K): [in a new window] |
FIG. 6. Structural comparison of type I and type II lumazine synthases and archaeal riboflavin synthase. (Top) Structural superposition of riboflavin synthase from M. jannaschii (MjaRS; PDB entry code 2B99 [45]), type I lumazine synthase from S. pombe (SpoLS; PDB entry code 1KYZ [16]), and type II lumazine synthase from B. abortus (BabRibH2; PDB entry code 1T13 [31]). The active sites are formed by two adjacent monomers of SpoLS with bound riboflavin (green; residues E17 to D112 and S113 to L158), BabRibH2 with bound 5-nitro-6-(D-ribitylamino)-2,4(1H,3H)-pyrimidinedione (gray; residues S12 to E106 and T107 to L156), and MjaRS with bound 6,7-dioxo-8-ribityllumazine, resembling the acceptor lumazine molecule (red; residues T2 to M90 and T91 to Y135). Secondary structure element labeling refers to SpoLS. The five topologically equivalent active sites of pentameric lumazine synthases are located at the interfaces between adjacent monomers of the pentamer; two of them (A and B) are shown. Ligands are drawn in the respective colors. (Bottom) Enlarged ligand-binding sites of all three enzymes, with Trp27 from SpoLS, Trp22 from BabRibH2, and Phe12 from MjaRS.
|
It should be noted that the eubacterial type I lumazine synthases that have been reported in the literature have low catalytic activities, in the range of about 200 to 300 nmol mg1 min1 when assayed near the optimum growth temperature of the cognate species (Table 1). By comparison with these data, the activity of the B. abortus type I lumazine synthase is more than 10-fold lower.
The generally low catalytic activity of lumazine synthases does not result from a particularly large free energy barrier of the reaction catalyzed. Quite to the contrary, the condensation of 3,4-dihydroxy-2-butanone 4-phosphate (substrate 4) and 5-amino-6-ribitylamino-2,4(1H,3H)-pyrimidinedione (substrate 2) is characterized by such a low energy barrier that it can proceed at an appreciable rate at room temperature in aqueous solution and at neutral pH (30). Thus, it is mandatory that blank assays without enzyme be run in order to correct for the contribution of the uncatalyzed reaction in lumazine synthase activity measurements (30). In the case of the lumazine synthase of B. abortus, with its particularly low catalytic rate, correction for the contribution of the uncatalyzed reaction becomes most important.
The low catalytic activities of the lumazine synthases are in no way unique among the enzymes of the riboflavin pathway. All activity values for the entire riboflavin pathways in E. coli, B. subtilis, and yeast are in the range of nanomoles per milligram per minute (5-7, 11, 12, 14, 20, 24, 28, 38, 46, 48). Since, at least in the case of lumazine synthase and riboflavin synthase, the inherent free energy barriers of the catalyzed reactions cannot be the reason for these low rates, we must assume, for lack of other arguments, that the selective pressure controlling this pathway favors the evolution of catalysts with low reaction rates. In fact, riboflavin is required in only small amounts, and excess production would unnecessarily deplete the precursor pools. Still, it remains an open question why the catalytic activity of the B. abortus type I enzyme is at the lower end of all documented lumazine synthases.
An earlier study showed that the decameric arrangement of type II B. abortus lumazine synthase is related to a very high Km for 3,4-dihydroxy-2-butanone 4-phosphate (substrate 4) (31). The in vivo concentration of 3,4-dihydroxy-2-butanone 4-phosphate is unknown. However, unless we assume that it is in the same numerical range as the Km of the type II enzyme, we must assume that the bulk of the 6,7-dimethyl-8-ribityllumazine (substrate 6) would be generated by the enzyme with the lower Km value, i.e., the type I lumazine synthase, whereas the type II enzyme could at best supply a minor amount of the overall riboflavin production.
Thus, the question of which selective pressures could have prevented the loss of the ribH2 gene (in case both lumazine synthase genes are evolutionarily old in Brucella-related organisms) or could have favored its more recent acquisition by horizontal gene transfer is still open. Interestingly, the type II lumazine synthase is an immunodominant antigen of B. abortus, and there is unpublished evidence that links this protein to Brucella virulence, suggesting that the type II lumazine synthase has evolved for a new, as-yet-unknown function.
We acknowledge R. Ugalde, D. Comerci, and Ines Marchesini for genomic B. abortus DNA and for early genomic analysis and Diana Posadas for helping us in the phylogenetic analysis.
|
|
|---|
| |||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
Copyright © 2009 by the American Society for Microbiology. For an alternate route to Journals.ASM.org, visit: http://intl-journals.asm.org | More Info»