Advertisement
JBC

HOME HELP FEEDBACK SUBSCRIPTIONS ARCHIVE SEARCH TABLE OF CONTENTS
 QUICK SEARCH:   [advanced]


     


Originally published In Press as doi:10.1074/jbc.M510254200 on November 24, 2005

J. Biol. Chem., Vol. 281, Issue 11, 7605-7613, March 17, 2006
This Article
Right arrow Abstract Freely available
Right arrow Full Text (PDF)
Right arrow All Versions of this Article:
281/11/7605    most recent
M510254200v1
Right arrow Submit a Letter to Editor
Right arrow Alert me when this article is cited
Right arrow Alert me when eLetters are posted
Right arrow Alert me if a correction is posted
Right arrow Citation Map
Services
Right arrow Email this article to a friend
Right arrow Similar articles in this journal
Right arrow Similar articles in PubMed
Right arrow Alert me to new issues of the journal
Right arrow Download to citation manager
Right arrowRequest Permissions
Citing Articles
Right arrow Citing Articles via HighWire
Right arrow Citing Articles via Google Scholar
Google Scholar
Right arrow Articles by Chen, S.-C.
Right arrow Articles by Liaw, S.-H.
Right arrow Search for Related Content
PubMed
Right arrow PubMed Citation
Right arrow Articles by Chen, S.-C.
Right arrow Articles by Liaw, S.-H.
Social Bookmarking
 Add to CiteULike   Add to Complore   Add to Connotea   Add to Del.icio.us   Add to Digg   Add to Reddit   Add to Technorati  
What's this?

Crystal Structure of a Bifunctional Deaminase and Reductase from Bacillus subtilis Involved in Riboflavin Biosynthesis*

Sheng-Chia Chen{ddagger}§1, Yuan-Chih Chang{ddagger}§1, Chao-Hsiung Lin{ddagger}, Chun-Hung Lin||, and Shwu-Huey Liaw{ddagger}**2

From the {ddagger}Structural Biology Program, §Institute of Biochemistry, and ||Faculty of Life Science, National Yang-Ming University, Taipei 11221, Taiwan, Institute of Biological Chemistry, Academia Sinica, Taipei 11529, Taiwan, and **Department of Medical Research and Education, Taipei Veterans General Hospital, Taipei 11217, Taiwan

Received for publication, September 19, 2005 , and in revised form, November 3, 2005.


    ABSTRACT
 TOP
 ABSTRACT
 INTRODUCTION
 EXPERIMENTAL PROCEDURES
 RESULTS AND DISCUSSION
 REFERENCES
 
Bacterial RibG is an attractive candidate for development of antimicrobial drugs because of its involvement in the riboflavin biosynthesis. The crystal structure of Bacillus subtilis RibG at 2.41-Å resolution displayed a tetrameric ring-like structure with an extensive interface of ~2400 Å2/monomer. The N-terminal deaminase domain belongs to the cytidine deaminase superfamily. A structure-based sequence alignment of a variety of nucleotide deaminases reveals not only the unique signatures in each family member for gene annotation but also putative substrate-interacting residues for RNA-editing deaminases. The strong structural conservation between the C-terminal reductase domain and the pharmaceutically important dihydrofolate reductase suggests that the two reductases involved in the riboflavin and folate biosyntheses evolved from a single ancestral gene. Together with the binding of the essential cofactors, zinc ion and NADPH, the structural comparison assists substrate modeling into the active-site cavities allowing identification of specific substrate recognition. Finally, the present structure reveals that the deaminase and the reductase are separate functional domains and that domain fusion is crucial for the enzyme activities through formation of a stable tetrameric structure.


    INTRODUCTION
 TOP
 ABSTRACT
 INTRODUCTION
 EXPERIMENTAL PROCEDURES
 RESULTS AND DISCUSSION
 REFERENCES
 
Flavin coenzymes are ubiquitous in all organisms because of their involvements in central metabolic pathways. Plants and many microorganisms obtain the precursor riboflavin by biosynthesis, whereas animals depend on nutritional sources. Numerous pathogenic microorganisms are unable to take up flavins from the environment and hence are absolutely dependent on their endogenous production. Therefore, the enzymes involved in riboflavin biosynthesis have the potential to become attractive candidates for the design of new defenses against antibiotic-resistant pathogens.

During riboflavin biosynthesis (1), GTP cyclohydrolase II first catalyzes the hydrolytic C-8 release of GTP to yield formate and pyrophosphate as side products. The product, 2,5-diamino-6-ribosylamino-4(3H)-pyrimidinone 5'-phosphate (compound 1) is converted into 5-amino-6-ribitylamino-2,4(1H,3H)-pyrimidinedione 5'-phosphate (compound 4) by deamination of the pyrimidine ring and NAD(P)H-dependent reduction of the ribose (Fig. 1). The deamination and reduction steps have been shown to proceed in the opposite order in yeast and Escherichia coli (2, 3). Most eubacteria contain a bifunctional protein; for instance Bacillus subtilis RibG (BsRibG)3 is composed of an N-terminal deaminase domain (D domain) and a C-terminal reductase domain (R domain) (4). In contrast, in fungi, plants and most archaea, these two enzymes are separate (5-7). In yeast, Rib7 is the corresponding reductase, whereas the C-terminal domain of Rib2 is responsible for deamination, in addition to N-terminal pseudouridine synthase activity.

The D domain of bacterial RibG was expected to belong to the cytidine deaminase (CDA) superfamily because of the two conserved signatures, H(C)XE and PCX8-9C. The CDA superfamily consists of the mononucleotide deaminases involved in nucleotide metabolism, and the RNA/DNA-editing deaminases involved in gene diversity and in antivirus defense (8, 9). The RNA/DNA-editing deaminases include A- to -I tRNA-specific adenosine deaminases (TADs), A- to -I adenosine deaminases acting on RNA, and C- to -U cytidine deaminases acting on RNA/DNA. These deaminases catalyze the hydrolytic deamination of cytosine, guanine, and adenine moieties and several of their therapeutically useful analogues. The available member structures reveal a virtually identical zinc-assisted deamination mechanism with the consensus histidine and cysteines acting as the zinc ligands, whereas the glutamate serves as a proton shuttle (10-14). However, the question of how nature has evolved the CDA fold into these various deaminases requires additional study.

Furthermore, the structural fold of the R domain was predicted by 3D-PSSM (15) to be similar to dihydrofolate reductase (DHFR). DHFR catalyzes the NADPH-utilizing reduction of dihydrofolate to tetrahydrofolate. Many DHFR inhibitors, such as methotrexate, pyrimethamine, and trimethoprim, have long been used clinically in the treatment of cancer, rheumatoid arthritis, malaria, and bacterial and fungal infection (16). Therefore, the R domain may become an important target for new drug design. To gain structural insights into the inhibitor design, substrate specificity, and evolution, we have solved the BsRibG structure at 2.41-Å resolution.


    EXPERIMENTAL PROCEDURES
 TOP
 ABSTRACT
 INTRODUCTION
 EXPERIMENTAL PROCEDURES
 RESULTS AND DISCUSSION
 REFERENCES
 
Protein Preparation—The His-tagged BsRibG was expressed using the pQE30 vector (Qiagen) in E. coli BL21 pLysS. The recombinant protein contains 12 additional vector residues (MRGSH6GS) at the N terminus. One-liter Luria broth (LB) cultures were grown at 37 °C to an A600 of 0.8 and induced with the addition of 1 mM isopropyl beta-D-thiogalactopyranoside. Cells were grown for another 5 h at 37°C before harvest. Cell pellets were resuspended in 40 ml of cold buffer A (20 mM Tris-HCl, 200 mM NaCl, pH 7.5) and lysed by French press. After removal of cellular debris by centrifugation at 39,000 x g at 4 °C for 30 min, the supernatant was applied to a 5-ml nickel-nitrilotriacetic acid column (Ni-NTA, Qiagen). The resin was washed with 40 mM imidazole in buffer A, and the protein was eluted with buffer A containing 500 mM imidazole. The protein fractions were collected, dialyzed against buffer B (20 mM Tris-HCl, pH 8.5), and loaded onto a SepharoseTM Q column (Amersham Biosciences). After washing with 230 mM NaCl in buffer B, the BsRibG was eluted with 260 mM NaCl and dialyzed against 20 mM HEPES (pH 7.5) and 5 mM dithiothreitol.


Figure 1
View larger version (7K):
[in this window]
[in a new window]
 
FIGURE 1.
The deamination and reduction steps in the riboflavin biosynthesis.

 
Protein Characterization—The enzyme activity assay was carried out as described previously (4). The substrate was prepared by GTP hydrolysis using the recombinant E. coli GTP cyclohydrolase II. The molecular mass in solution was estimated by a Beckman-Coulter XL-A analytical Ultracentrifuge with an An60Ti rotor. Sedimentation velocity was performed at 20 °C and 40,000 rpm with standard double sector centerpieces. The UV absorption of the cells was scanned every 5 min for 2 h, and the data were analyzed using the SedFit program (17).

Protein Crystallization—The initial crystallization screening was performed with screening kits using the hanging-drop vapor diffusion method at 22 °C. The hanging drops were mixtures of 2 µl of reservoir solution and 2 µl of protein solution. Protein crystals could be obtained under several reservoir solutions, and the best crystals were grown in 26.6% polyethylene glycol 400, 190 mM MgCl2, 5% glycerol, and 95 mM HEPES (pH 7.5) with a protein solution of 20-25 mg/ml. Crystals appeared and reached their final dimensions in 1 week at 15 °C. The NADPH derivative was prepared by soaking crystals for 4 days in reservoir solution containing 10 mM NADPH. X-ray diffraction data were collected and processed at beamlines BL12B2 at SPring-8 (Harima, Japan) and NW12 at the Photon Factory (Tsukuba, Japan). The crystals belong to the P212121 space group with 1 tetramer/asymmetric unit.

Structure Determination—The BsRibG structure was determined using single-wavelength anomalous dispersion (SAD) of the endogenous zinc ion. Four Zn2+ positions were identified and refined with SOLVE (18). The initial phase was improved by direct-method phasing refinement using OASIS (19) before density modification. The initial electron density map showed an interpretable density of ~70% of the tetramer. About 60% of the initial tetrameric model was automatically built using RESOLVE (18). Because the structural fold was predicted to be similar to yeast cytosine deaminase (yCD) and Thermotoga maritima DHFR (TmDHFR) (10, 20), both structures were used to assist in the manual building of the atomic model using TURBO-FRODO (21), and refinement was carried out against the native data to 2.41-Å resolution using crystallography NMR software (22). The statistics for data collection and refinement are summarized in Table 1. About 88% of the residues are in the most favored regions of the Ramachandran plot, with the remaining in the additional allowed regions.


View this table:
[in this window]
[in a new window]
 
TABLE 1
Statistics for data collection and structural refinement

Values in parentheses are for the highest resolution shell. SAD, single-wavelength anomalous dispersion.

 


Figure 2
View larger version (40K):
[in this window]
[in a new window]
 
FIGURE 2.
Sedimentation velocity analysis of BsRibG. The enzyme concentration was 0.3 mg/ml in 20 mM HEPES (pH 7.5). The three panels represent the optical traces (top), the residues of the model fitting (middle), and the sedimentation coefficient distribution (bottom), respectively. The ultracentrifugation analysis demonstrated that BsRibG exists in solution as a tetramer with a protein concentration ranging from 80 to 300 µg/ml.

 
Sequence Alignment—To identify the unique signatures for each family member for gene annotation, sequence similarity searches were conducted by PSI-BLAST (23), and multiple sequence alignment of the homologous sequences was performed by ClustalW (24). Because of a conservative hydrophobic core and the consensus HXE and PCX2-9C, a structural-based sequence alignment between the CDA members was feasible and was carried out by manual editing according to the available structures. The conserved residues in each family member were mapped onto the known structures to reveal their potential involvement in structural integrity or the enzyme catalysis. Figs. 3, 4, A and B, and 6A were generated by MolScript (25) and Raster3D (26), Fig. 5A by BobScript (25), and Fig. 5B by LigPlot (27).


    RESULTS AND DISCUSSION
 TOP
 ABSTRACT
 INTRODUCTION
 EXPERIMENTAL PROCEDURES
 RESULTS AND DISCUSSION
 REFERENCES
 
Overall Structure—Analytical ultracentrifugation experiments clearly demonstrated that BsRibG exists as a tetramer in solution as well as in crystal form, where the enzyme forms a tetrameric ring-like structure (Figs. 2 and 3a). The current tetrameric model contains residues 1-359 for each subunit. Each monomer is composed of two separate functional domains (Fig. 3b). The N-terminal D domain (residues 1-143) consists of a central five-stranded beta-sheet (beta1-beta5) with beta1 running antiparallel to the others. The beta-sheet is sandwiched by two helices ({alpha}A and {alpha}E) on one side and by three helices ({alpha}B, {alpha}C, and {alpha}D) on the other side. The C-terminal helix {alpha}F extends away from the D domain and connects to the R domain. The R domain (residues 146-359) is composed of a large nine-stranded beta-sheet (betaA-betaH and betaD') with the C-terminal strand betaH running antiparallel to the others. The beta-sheet is flanked by five {alpha}-helices ({alpha}B, {alpha}C, {alpha}D', {alpha}E, and {alpha}F). The secondary structure elements of the R domain are numbered as for DHFRs (28).

The four subunits in the crystal asymmetric unit did not show significant differences between each individual domain except for several loops (root-mean-square deviations (r.m.s.d.) of 0.38-0.46 and 0.59-0.73 Å for the backbones of residues 2-139 and residues 146-359, respectively). However, the relative orientations between the D and R domains are slightly different, resulting in a weak noncrystallographic symmetry. Molecules A and B interact with each other through their D domains with a buried surface area of ~650 Å2 per D domain (the D interface) (Fig. 4A), whereas molecule A makes extensive contacts with molecule C through their respective R domains with a buried area of ~1750 Å2 per R domain (the R interface) (Fig. 4B). There are no contacts between molecules A and D in the tetramer.

The D interface is made up mainly of the N-terminal two helices {alpha}A and {alpha}B, the beta2 strand, and the connected loops (residues 4-19 and 35-61). There are 14 direct hydrogen bonds between the protein atoms across the interface and hydrophobic patches formed by Leu8, Leu12, Ile36, Met39, Leu43, and Met57. The anti-parallel disposition of the two pseudodyad-related {alpha}A helices also contributes to the interface by dipole-dipole interactions. Interestingly, the side chains of His56 from the two D domains stack very well, with a distance of 3.3-3.4 Å between the aromatic rings. The R interface mainly is made up of: the two large loops, LbetaA-{alpha}B (residues 156-169) and LbetaF-{alpha}G (residues 314-338); and the C-terminal residues 339-358. There are 18 direct hydrogen bonds between the protein atoms across the interface, and extensive hydrophobic patches are by Try169 formed, Pro314, Lys315, Leu316, Ile317, Leu325, Phe331, Met334, Val337, Leu339, Leu340, Phe342, Ile345, Ile352, and Leu354.


Figure 3
View larger version (63K):
[in this window]
[in a new window]
 
FIGURE 3.
Structure of BsRibG. Ribbon views of the tetramer (a) and monomer (b). The tightly bound zinc ion in D domain is shown as a sphere (magenta) with the cofactor NADPH (black) and the modeled substrates (cyan) as ball- and -stick representations. Molecules A-D in the tetramer, as referenced in the text, are colored in red, blue, green, and black, respectively. The D and R domains are made up of a central beta-sheet flanked by {alpha}-helices. These two domains have few contacts with each other and hence can fold independently.

 
Structural Conservation in the CDA Superfamily—As expected, the D domain displays high structural homology to the available structures of the CDA members including yCD, B. subtilis guanine deaminse (BsGD), CDAs, T4 bacteriophage dCMP deaminse (T4dCMPD), and Aquifex aeolicus TADA (AaTADA), and subdomains 2 and 4 of the chicken AICAR transformylase domain (10-14, 29). The TADAs in prokaryotes and the TAD2/TAD3 heterodimers in eukaryotes are responsible for the I34 alteration at the wobble position of the tRNA anticodon, whereas TAD1 creates the unique 1-methyl-I37 in eukaryotic tRNAAla (8). Detailed structural comparisons reveal a common three-layer {alpha}/beta/{alpha} structure, in which five beta-strands (beta1-beta5) and three helices ({alpha}A-{alpha}C) correspond closely, whereas the remainder varies across the different deaminases (Fig. 5A). The main chain atoms of the 65-70 structurally equivalent residues are overlaid with an r.m.s.d. of 0.75 to 1.35 Å and with 8-24% sequence identity. These structural elements in the CDA fold form a conservative hydrophobic core, which is also preserved in the AICAR transformylase domain, implying that the hydrophobic core has been highly conserved throughout evolution.

The active site of the D domain contains one tightly bound endogenous zinc ion, for which the anomalous data provided sufficient phase information for structure determination. The zinc ion is tetrahedrally coordinated by His49 N{epsilon}1 (2.0 Å), Cys74 S{gamma} (2.4 Å), Cys83 S{gamma} (2.3 Å), and a water molecule (2.0 Å).The zinc-bound water molecule interacts with Glu51 O{epsilon}2 (2.5 Å). The active-site architecture resembles those of the CDA members, which share a similar zinc-assisted deamination mechanism with a virtually identical interaction network between the common moiety of the pyrimidine ring of the substrate, the zinc ion, the zinc-bound water molecule, the zinc ligands, and the base glutamate (Fig. 5B). In addition, the active-site cavities of these deaminases are mainly made up of the C-terminal tail and the loops connecting the {alpha}A-beta1, beta2-{alpha}B, beta3-{alpha}C, and beta4-{alpha}D (Fig. 5C).

Based on the structural comparison of the active-site cavities, the substrate of the D domain was modeled into the active site through superposition of the nucleobase rings because of the highly conserved interaction networks surrounding the target amino group (Fig. 5B). The model was then subjected to energy minimization with crystallography NMR software. Simulation of the complex structure suggested that the nucleophilic OH-2 group of the pyrimidine ring coordinates to the catalytic zinc ion and interacts with Glu51O{epsilon}2 and Cys74 N. The NH-3 group hydrogen bonds with Glu51 O{epsilon}1, the O-4 atom with Ala50 N and His42 N{delta}, and the NH2-5 with Asn23 O{delta}1 and His42 N{delta}1. The two hydroxyl groups of the ribose have close contacts with the side chains of Asp101 and Asn103. The phosphate moiety forms salt bridges with His76 and Lys79, located in the unique insertion between the two zinc ligand cysteines, and these interactions are essential for the deamination activity because the enzyme cannot utilize the dephosphorylated form as substrate (3). These predicted substrate-binding residues are all highly conserved in the eubacterial RibGs. However, the fungal deaminases such as yeast Rib2 (yRib2) apparently contain different substrate-interacting residues because of their distinct substrate, which has an open ribityl group instead of a cyclic ribose (Figs. 1 and 5C).


Figure 4
View larger version (51K):
[in this window]
[in a new window]
 
FIGURE 4.
Subunit interfaces. Stereo views of the D interface (A) and the R interface (B). Two distinct subunit interfaces are formed by the D and R domains with total buried areas of ~1300 and ~3500 Å2, respectively. The D interface is made up mainly of the beta2 strand and helices {alpha}A and {alpha}B, whereas the R interface is made up of the two large loops, LbetaA-{alpha}B and LbetaF-betaG, and the C-terminal tail.

 
Structural Divergence in the CDA Superfamily—The CDA members exist as an oligomer. All of the available member structures except for BsRibG utilize helices {alpha}B-{alpha}D and surrounding loops for oligomerization. BsGD, yCD, and AaTADA display similar dimeric structures (10, 11, 14). The swapping of the C-terminal segment in BsGD causes additional dimeric contacts including the beta5 strand and helices {alpha}A and {alpha}E. CDAs form tetramers in which one subunit interacts with the other three subunits (12). The hexameric T4dCMPD contains two types of intersubunit interfaces (13). In contrast, the D interface of BsRibG is very distinct from the others and possesses the fewest contacts, with a total buried area of ~1300 Å2 (Fig. 4A). In addition to the loops and helices {alpha}A and {alpha}B, RibG includes the beta2 strand in the D interface. The distinct intersubunit orientation in the D interface also separates the active sites away from each other with an inter-zinc distance of 30 Å. Notably, the shortest inter-zinc distance in the other five deaminases is about 14-15 Å.


Figure 5
View larger version (57K):
[in this window]
[in a new window]
 
FIGURE 5.
Structural conservation and divergence of the CDA superfamily. A, stereo view of structural superposition of BsRibG (red), yCD (blue) (10), T4dCMPD (green) (13), and AaTADA (yellow) (14). The zinc ion is displayed as a sphere (magenta) with the yCD inhibitor 3,4-dihydrouracil (DHU)(cyan) shown as ball- and -stick representations. These deaminases share the conserved beta-sheet and helices {alpha}A-{alpha}C and even a part of the {alpha}E helix. T4dCMPD contains a ~60-residue insertion, which folds into two helices ({alpha}B' and {alpha}B'') and flexibleloops. The C-terminal tail of yCD folds backward to limit the pocket size, whereas those of the remaining members swing away to enlarge the active-site cavity. B, superposition of the active sites of BsRibG (magenta), yCD (cyan), and T4dCMPD (green). The residue numbering is labeled in the same color for each protein. The deaminases display highly conserved interaction networks surrounding the target amino group of the nucleobase ring. In contrast, each member contains its own unique substrate recognition residues. C, multiple sequence alignment of some CDA members. The Arabidopsis thaliana deaminase (At363) (6) involved in riboflavin biosynthesis is also included. The GenBankTM accession codes are listed in the right column. Secondary structure elements for BsRibG are labeled (s.s). The number of residues in gaps is indicated in parentheses, and the protein length is in brackets. The superfamily signatures, H(C)XE and PCX2-9C, are shaded in cyan, whereas the residues for the conserved hydrophobic core are in yellow. In addition, the substrate-binding residues in the known complex structures are shaded in red, and those that are predicted are in blue. The unique signatures for each member are highlighted in italics, and the residues involved in loss-of-function point mutants of AID are shaded in magenta.

 
The C-terminal segment beyond the beta4 strand is quite diverse and may make a major contribution to the structural plasticity and functional diversity among the CDA members (Fig. 5A). For instance, the mononucleobase deaminase yCD and the tRNA-editing deaminase AaTADA unexpectedly superimposed very well, with an r.m.s.d. of 1.08 Å for 115 C{alpha} atoms with 22% sequence identity. The AaTADA structure was solved by molecular replacement using yCD as a search model (14). The major structural difference around the active-site cavity is the C-terminal helix. In yCD as well as BsGD, genetic changes to alter the substrate specificity are through an introduction of substrate recognition residues at the C-terminal tail (Asp155 in yCD and Tyr156 in BsGD), which then forms a "flap," capping and hence narrowing the opening of the active-site cavity upon substrate binding to limit the pocket size for the nucleobase (10, 11). In contrast, the C-terminal tail in AaTADA, as well as in CDAs, dCMPDs, and the D domain of BsRibG, swings away to enlarge the active-site cavity for their larger substrates (12-14). In addition to the diverse C-terminal segment, each family member has some unique substrate recognition residues; these are discussed below in detail.

Structure-based Sequence Alignment of the CDA Members—In combination with a sequence-structure analysis, a structure-based sequence alignment of the CDA members was constructed (Fig. 5C). Multiple sequence alignments reveal that there are unique member signatures that are useful for gene annotation. The unique substrate recognition residues for the cytosine deaminases are WXXDI at the C terminus, and those for the guanosine deaminases are FDD between helices {alpha}D1 and {alpha}D2 and a tyrosine at the C terminus. The CDAs also have unique substrate-binding residues, namely FXV at the N terminus of the beta1 strand and NXE at the C terminus of the beta2 strand. The absence of the {alpha}D helix between the beta4 and beta5 strands causes both strands to become antiparallel and the C terminus to be in the opposite direction. dCMPDs possess GYNG on the beta2 strand and a diverse extra 20-60-residue insertion between the beta2 strand and the {alpha}B helix. Interestingly, most phosphate-interacting residues are located in this insertion and are not conserved across dCMPDs. RibGs have NPXVG at the N terminus of the beta1 strand and an extra 8-9-residue loop between the two zinc-ligating cysteines (PCX8-9C). The extra loop is highly conserved in RibGs and is predicted to interact with the phosphate moiety of the substrate.

Even though the CDA members display substantial sequence diversity (15-25% sequence identity), our structure-sequence analysis suggests that comparative modeling is feasible. For instance, bacterial TADAs and eukaryotic TAD2s share a unique signature, EXPVG, at the N terminus of the beta1 strand, in which the glutamate is also conserved in TAD1. Our comparative modeling suggests that this glutamate residue may serve as an adenine recognition residue for these A- to -I tRNA deaminases.


Figure 6
View larger version (28K):
[in this window]
[in a new window]
 
FIGURE 6.
The NADPH-binding site in the R domain. A, the 2Fo-Fc electron density map for NADPH contoured at 1.5 {sigma} level and is shown in cyan. B, schematic diagram of BsRibG interactions with the NADPH cofactor. Hydrogen bonds are presented as dashed lines; the interatomic distances are given in angstroms. "Radiating" spheres indicate hydrophobic contacts between the cofactor and the surrounding residues.

 
Possible Effects of the AID Point Mutants in Hyper-IgM Patients—To date, apolipoprotein B mRNA-editing catalytic subunit 1 (APOBEC1) and its sequence homologues, activation-induced deaminase (AID), APOBEC2, and the tandem repeats APOBEC3A to -3H, are the only identified C- to -U RNA/DNA-editing deaminases in humans (9). AID is an essential B cell-specific factor required for antibody maturation, whereas several APOBECs are involved in defense against a broad range of retroviruses. The AID and APOBEC1 structures have been modeled with an N-terminal catalytic domain, a linker, and a C-terminal pseudocatalytic domain based on the CDA structures (30, 31). However, both proteins seem to contain the {alpha}D helix (Fig. 5C), and hence the beta4 and beta5 strands are parallel and the direction of the C terminus is opposite from that of CDAs. In addition, both AID and APOBEC1 are too short to accommodate two CDA folds in one protein molecule. We have modeled residues 1-160 of human AID into a CDA fold consisting of the five beta-strands surrounded by helices {alpha}A-{alpha}E. APOBECs contain more conserved aromatic residues than other members with a unique signature TWY(F)XSWSPCX2C around the beta3 strand.


Figure 7
View larger version (56K):
[in this window]
[in a new window]
 
FIGURE 7.
Structural conservation and divergence between RibG and DHFR. A, stereo view of structural superposition of the R domain of BsRibG (red) and TmDHFR (green) (20). The R domain contains an extra ~25-residue insertion, which folds into the {alpha}D' helix and the betaD' strand. B, multiple sequence alignment. Three other reductases involved in the riboflavin biosynthesis from A. thaliana (At599), yeast (5), and the Archaea Methanocaldococcus jannaschii (7) and two DHFRs from hyperthermophilic T. maritima and mesophilic E. coli are also included. The SwissProt accession codes are listed in the right column. Residues involved in NADPH binding are shaded in magenta, whereas the residues for the conserved hydrophobic core are in yellow. The substrate-binding residues in the known complex structures are shaded in red, and those that are predicted are in blue. DHFRs and RibG share four conserved regions for NADPH binding (motifs A-D). Their residues involved in substrate recognition are located in similar spatial positions in the betaA strand, the {alpha}B helix, and the loop between the {alpha}C helix and the betaC strand.

 
Most residues involved in loss-of-function point mutants of AID in patients with hyper-IgM syndrome type 2 (32) are highly conserved in APOBECs (Fig. 5C). His56 and Cys87 ligate the zinc ion. Arg24, prior the beta1 strand, and Arg112 to at the C terminus of the beta4 strand may interact with the phosphate group of the edited cytidylate, with respect to corresponding Ser21 and Tyr153 the in T4dCMPD. Phe151 in the {alpha}E helix may form close contacts with the substrate. The conserved FFX3R motif in the {alpha}E helix in bacterial TADAs has been shown to play a critical role in A34- to -I34 deamination (33). In contrast, Trp80, Leu106, and Met139, located in the beta3, beta4, and beta5 strands, respectively, may be responsible for construction of the conserved hydrophobic core. The truncated C147X mutant could not be expressed in 293T cells (32), implying that this mutant might not form a stable CDA fold. The 181X and 190X mutants retain deamination activity but are defective for class switch recombination and normal nuclear distribution, indicating the functional roles of the extra C-terminal segment in AID. Similarly, Ser3 and Lys10 at the N-terminal {alpha}A helix are expected to be exposed at the protein surface and might be important for nuclear location signaling or for interaction with associate proteins.

NADPH Binding in the R Domain—NADPH was identified by its strong electron density and is firmly embedded at the active site of the R domain (Fig. 6A). The enzyme did not show significant structural changes upon NADPH binding, with an r.m.s.d. of 0.4 Å for all of the backbone atoms. The cofactor is bound in an extended conformation with extensive interactions with the protein. The adenine N-3 and N-6 atoms interact with Thr221 and Ser298, respectively (Fig. 6B). The phosphate group at the ribose O-3' makes close contacts with Thr221. These interactions seem not to contribute significantly to the cofactor binding, because NADH as well as NADPH can serve as the coenzyme (4). The pyrophosphate moiety interacts with a constellation of backbone NH groups from Gly194, Thr195, Gly292, Ser293, Ala294, and Val295 and the side chain of Thr195. These residues are located at the N termini of helices {alpha}C and {alpha}E, and the helix dipoles may partially neutralize the negative charges of the phosphate groups. The nicotiamide amide group complementarily interacts with the amide backbone of Ala153, whereas the nicotinamide ribose OH-2' and OH-3' contact with Asp199 and Gly165, respectively.

Structural Conservation and Divergence between BsRibG and DHFRs—A structural homology search by DALI (34) revealed that the R domain displays significant structural similarity to DHFRs (Fig. 7A). Despite low sequence identity, which is as little as 20%, all the DHFR structures exhibit a virtually identical fold consisting of a central eight-stranded beta-sheet (designated betaA-betaH) and four flanking helices ({alpha}B, {alpha}C, {alpha}E, and {alpha}F) (28). The R domain of BsRibG contains an extra 25-residue insertion between the betaD strand and the {alpha}E helix, forming the {alpha}D' helix and the betaD' strand. This insertion is quite a distance away from the active site and is diverse across the bacterial and fungal reductases.

TmDHFR is the most stable DHFR isolated thus far. The most prominent feature of TmDHFR is its highly stable dimer, which has been shown to make a major contribution to the high intrinsic stability (20). The R-R structure of BsRibG is similar to the dimeric TmDHFR structure. Both TmDHFR and BsRibG utilize similar regions for the subunit association. The contact areas of the R interface are larger than those of TmDHFR, ~3500 Å2 versus ~3050 Å2, perhaps due to a larger LbetaF-betaG loop in BsRibG. Structural superposition of RibG and TmDHFR reveals an r.m.s.d. of 1.5 Å for 120 C{alpha} atoms with 25% sequence identity. The strong conservation of the tertiary structures suggests that the two reductases involved in the riboflavin and folate biosyntheses are descended from a single ancestral gene and there by define a new superfamily.

DHFRs and BsRibG share four conserved regions for NADPH binding, motifs A-D (Fig. 7B). In motif A, the two hydrogen bonds between the amide groups of the nicotinamide ring and an alanine residue at the C terminus of the betaA strand (Ala153 in BsRibG) are strictly conserved. Both proteins use motifs B and D for the pyrophosphate binding, and motif C for the phosphate group at the adenosine O-3'. The interactions between DHFRs and the 3'-phosphate group are so extensive via 3-4 residues that the enzyme has a strict NADPH dependence. Both DHFR and RibG utilize many main chain atoms for NADPH binding, and hence the corresponding sequences of the NADPH-interacting residues have diverged during evolution, but their spatial positions have remained convergent in the four regions. A conserved cis peptide bond occurs at the two consecutive glycine residues in motif D. In contrast, both enzymes form diverse interactions with the remaining areas including the adenine ring and the riboses.

The enzymatic mechanism of DHFR has been chemically and structurally studied in detail (28, 35). A complete catalytic scheme through five kinetically observable intermediates has been proposed. Similarly, RibG is expected to catalyze the reduction of a cyclic ribosyl into an open ribityl group by hydride transfer from the C-4 atom of the nicotinamide ring of NAD(P)H to the C-1' of the ribose with concomitant protonation of O-5'. The binding architecture of the nicotinamide ring is virtually identical in DHFR and RibG. Thus, based on the structural comparison, the substrate of the R domain was modeled into the active site with subsequent energy minimization. Simulation of the complex structure suggested that the O-2 atom of the pyrimidine ring interacts with Lys151 N{zeta} and the NH-3 and O-4 with Thr172 O{gamma}1. Thr195 and Asp199 make close contacts with the two hydroxyl groups of the ribose. Arg176, Arg183, and Arg206 form salt bridges with the phosphate group, interactions that are essential for the reduction activity because the enzyme cannot utilize the dephosphorylated form as substrate (3). These substrate-binding residues are highly conserved in the reductases and are located in spatial positions similar to those in DHFRs (Fig. 7B). Lys151 is highly conserved in most eubacteria and plants but not in fungi and archaea, implying that this residue distinguishes the O-2 atom from the NH2-2 group of the pyrimidine moiety, which is the only difference in the substrates (Fig. 1). A highly conserved glutamate in motif D (Glu290 in BsRibG) may assist in proton transfer during the hydride transfer.

Fusion of Two Enzyme Domains in RibG—There is no evidence for any dependence of the active sites within the RibG tetramer. The active sites of the D and R domains within the same molecule are ~40 Å apart if measured from the zinc ion to the nicotiamide C-4 atom. The openings of the D and R active sites point in the opposite directions. Unlike the bifunctional enzyme of DHFR and thymidylate synthase (36), the present BsRibG structure does not show a channel for substrate transport from the D domain to the R domain for the sequential reaction. In addition, the active sites of the dimeric TmDHFR have been shown to be independent (20). Unlike the monomeric DHFRs, TmDHFR does not show significant structural changes on binding of NADPH and the inhibitor methotrexate, because the active-site loops (LbetaA-{alpha}B and LbetaF-betaG) in the monomeric enzymes are rearranged on the dimeric interfaces. Similarly, NADPH binding does not induce any significant conformational changes in BsRibG, and the corresponding loops are also involved in the R interface.

Previous deletion mutants have demonstrated that the N-terminal 147 residues and the C-terminal 248 residues of BsRibG were sufficient for their respective enzyme activities (4). However, these truncated proteins could not be isolated because of poor stability. The present structure reveals that the D and R domains make only a few contacts, including four hydrogen bonds by Glu21, Gln135, His139 Gln187, Lys212, Gln306 and Lys357 and small hydrophobic contacts by Lys136, Phe140, Tyr147, and Met285.If the {alpha}F helix (residues 135-143) is considered as an interdomain linker, the D and R domain is even less interactive. Therefore, even though the two enzyme domains can fold independently, the domain fusion is crucial for the enzyme activities through formation of a stable tetrameric structure.


    FOOTNOTES
 
The atomic coordinates and structure factors (codes 2B3Z and 2D5N) have been deposited in the Protein Data Bank, Research Collaboratory for Structural Bioinformatics, Rutgers University, New Brunswick, NJ (http://www.rcsb.org/).

* This study was supported by National Science Council Grant NSC94-2311-B010-017. The costs of publication of this article were defrayed in part by the payment of page charges. This article must therefore be hereby marked "advertisement" in accordance with 18 U.S.C. Section 1734 solely to indicate this fact. Back

1 These authors contributed equally to this work. Back

2 To whom correspondence should be addressed: Faculty of Life Science, National Yang-Ming University, Taipei, Taiwan 11221. Tel.: 886-2-2826-7278; Fax: 886-2-2820-2449; E-mail: shliaw{at}ym.edu.tw.

3 The abbreviations used are: BsRibG, B. subtilis RibG; BsGD, B. subtilis guanine deaminase; CDA, cytidine deaminase; yCD, yeast cytosine deaminase; T4dCMPD, T4 bacteriophage dCMP deaminase; TAD, tRNA-specific adenosine deaminase; AaTADA, A. aeolicus tRNA-specific adenosine deaminase; DHFR, dihydrofolate reductase; TmDHFR, T. maritima dihydrofolate reductase; AID, activation-induced deaminase; APOBEC, apolipoprotein B mRNA-editing catalytic subunit; r.m.s.d., root-mean-square deviation; AICAR, aminoimidazole-4-carboxamide ribonucleotide. Back


    ACKNOWLEDGMENTS
 
We thank Dr. R. Kirby for critical reading of the manuscript and Dr. H. Chang for performing the ultracentrifuge experiments.



    REFERENCES
 TOP
 ABSTRACT
 INTRODUCTION
 EXPERIMENTAL PROCEDURES
 RESULTS AND DISCUSSION
 REFERENCES
 

  1. Bacher, A., Eberhardt, S., Fischer, M., Kis, K., and Richter, G. (2000) Annu. Rev. Nutr. 20, 153-167[CrossRef][Medline] [Order article via Infotrieve]
  2. Oltmanns, O., and Bacher, A. (1972) J. Bacteriol. 110, 818-822[Abstract/Free Full Text]
  3. Burrows, R. B., and Brown, G. M. (1978) J. Bacteriol. 136, 657-667[Abstract/Free Full Text]
  4. Richter, G., Fischer, M., Krieger, C., Eberhardt, S., Lüttgen, H., Gerstenschläger, I., and Bacher, A. (1997) J. Bacteriol. 179, 2022-2028[Abstract/Free Full Text]
  5. Baur, A., Schaaff-Gerstenschlager, I., Boles, E., Miosga, T., Rose, M., and Zimmermann, F. K. (1993) Yeast 9, 289-293[CrossRef][Medline] [Order article via Infotrieve]
  6. Fisher, M., Romisch, W., Saller, S., Illarionov, B., Richter, G., Rohdich, F., Eisenreich, W., and Bacher, A. (2004) J. Biol. Chem. 279, 36299-36308[Abstract/Free Full Text]
  7. Graupner, M., Xu, H., and White, R. H. (2002) J. Bacteriol. 184, 1952-1957[Abstract/Free Full Text]
  8. Keegan, L. P., Leroy, A., Sproul, D., and O'Connell, M. A. (2004) Genome Biol. 5, 209-218[CrossRef][Medline] [Order article via Infotrieve]
  9. Harris, R. S., and Liddament, M. T. (2004) Nat. Rev. Immunol. 4, 868-877[CrossRef][Medline] [Order article via Infotrieve]
  10. Ko, T. P., Lin, J. J., Hu, C. Y., Hsu, Y. H., Wang, A. H. J., and Liaw, S. H. (2003) J. Biol. Chem. 278, 19111-19117[Abstract/Free Full Text]
  11. Liaw, S. H., Chang, Y. J., Lai, C. T., Chang, H. C., and Chang, G. G. (2004) J. Biol. Chem. 279, 35479-35485[Abstract/Free Full Text]
  12. Johansson, E., Mejlhede, N., Neuhard, J., and Larsen, S. (2002) Biochemistry 41, 2563-2570[CrossRef][Medline] [Order article via Infotrieve]
  13. Almog, R., Maley, F., Maley, G. F., Maccoll, R., and Van Roey, P. (2004) Biochemistry 43, 13715-13723[CrossRef][Medline] [Order article via Infotrieve]
  14. Kuratani, M., Ishii, R., Bessho, Y., Fukunaga, R., Sengoku, T., Shirouzu, M., Sekine, S., and Yokoyama, S. (2005) J. Biol. Chem. 280, 16002-16008[Abstract/Free Full Text]
  15. Kelley, L. A., MacCallum, R. M., and Sternberg, M. J. (2000) J. Mol. Biol. 299, 499-520[Medline] [Order article via Infotrieve]
  16. Kompis, I. M., Islam, K, and Then, R. L. (2005) Chem. Rev. 105, 593-620[CrossRef][Medline] [Order article via Infotrieve]
  17. Schuck, P., Perugini, M. A., Gonzales, N. R., Howlett, G. J., and Schubert, D. (2002) Biophys. J. 82, 1096-1111[Medline] [Order article via Infotrieve]
  18. Terwilliger, T. C. (2003) Methods Enzymol. 374, 22-37[Medline] [Order article via Infotrieve]
  19. Wang, J.W., Chen, J. R., Gu, Y. X., Zheng, C. D., Jiang, F., Fan, H. F., Terwilliger, T. C., and Hao, Q. (2004) Acta Crystallogr. Sect. D 60, 1244-1253[CrossRef][Medline] [Order article via Infotrieve]
  20. Dams, T., Auerbach, G., Bader, G., Jacob, U., Ploom, T., Huber, R., and Jaenicke, R. (2000) J. Mol. Biol. 297, 659-672[CrossRef][Medline] [Order article via Infotrieve]
  21. Roussel, A., and Cambillau, C. (1991) in Silicon Graphics Geometry Partners Directory 81, Silicon Graphics Corp., Mountain View, CA
  22. Brünger, A. T., Adams, P. D., Clore, G. M., DeLano, W. L., Gros, P., Grosse-Kunstleve, R. W., Jiang, J. S., Kuszewski, J., Nilges, M., Pannu, N. S., Read, R. J., Rice, L. M., Simonson, T., and Warren, G. L. (1998) Acta Crystallogr. Sect. D 54, 905-921[CrossRef][Medline] [Order article via Infotrieve]
  23. Schäffer, A. A., Aravind, L., Madden, T. L., Shavirin, S., Spouge, J. L., Wolf, Y. I., Koonin, E. V., and Altschul, S. F. (2001) Nucleic Acids Res. 29, 2995-3005
  24. Chenna, R., Sugawara, H., Koike, T., Lopez, R., Gibson, T. J., Higgins, D. G., and Thompson, J. D. (2003) Nucleic Acids Res. 31, 3497-3500[Abstract/Free Full Text]
  25. Esnouf, R.M. (1999) Acta Crystallogr. Sect. D 55, 938-940[CrossRef][Medline] [Order article via Infotrieve]
  26. Merritt, E. A., and Bacon, D. J. (1997) Methods Enzymol. 277, 505-524[Medline] [Order article via Infotrieve]
  27. Wallace, A. C., Laskowski, R. A., and Thornton, J. M. (1995) Protein Eng. 8, 127-134[Abstract/Free Full Text]
  28. Schnell, J. R., Dyson, H. J., and Wright, P. E. (2004) Annu. Rev. Biophys. Biomol. Struct. 33, 119-140[CrossRef][Medline] [Order article via Infotrieve]
  29. Wolan, D. W., Greasly, S. E., Beardsley, G. P., and Wilson, I. A. (2002) Biochemistry 41, 15505-15513[CrossRef][Medline] [Order article via Infotrieve]
  30. Navaratnam, N., Fujino, T., Bayliss, J., Jarmuz, A., How, A., Richardson, N., Somasekaram, A., Bhattacharya, S., Carter, C., and Scott, J. (1998) J. Mol. Biol. 275, 695-714[CrossRef][Medline] [Order article via Infotrieve]
  31. Xu, H., Svarovskaia, E. S., Barr, R., Zhang, Y., Khan, M. A., Strebel, K., and Pathak, V. K. (2004) Proc. Natl. Acad. Sci. U. S. A. 101, 5652-5657[Abstract/Free Full Text]
  32. Ta, V. T., Nagaoka, H., Catalan, N., Durandy, A., Fischer, A., Imai, K., Nonoyama, S., Tashiro, J., Ikegawa, M., Ito, S., Kinoshita, K., Muramatsu, M., and Honjo, T. (2003) Nat. Immunol. 4, 843-848[CrossRef][Medline] [Order article via Infotrieve]
  33. Elias, Y., and Huang, R. H. (2005) Biochemistry 44, 12057-12065[CrossRef][Medline] [Order article via Infotrieve]
  34. Dietmann, S., Park, J., Notredame, C., Heger, A., Lappe, M., and Holm, L. (2001) Nucleic Acids Res. 29, 55-57[Abstract/Free Full Text]
  35. Fierke, C. A., Johnson, K. A., and Benkovic, S. J. (1987) Biochemistry 26, 4085-4092[CrossRef][Medline] [Order article via Infotrieve]
  36. Yuvaniyama, J., Chitnumsub, P., Kamchonwongpaisan, S., Vanichtanankul, J., Sirawaraporn, W., Taylor, P., Walkinshaw, M. D., and Yuthavong, Y. (2003) Nat. Struct. Biol. 10, 357-365[CrossRef][Medline] [Order article via Infotrieve]

Add to CiteULike CiteULike   Add to Complore Complore   Add to Connotea Connotea   Add to Del.icio.us Del.icio.us   Add to Digg Digg   Add to Reddit Reddit   Add to Technorati Technorati    What's this?


This article has been cited by other articles:


Home page
J. Biol. Chem.Home page
S.-C. Chen, Y.-H. Lin, H.-C. Yu, and S.-H. Liaw
Complex Structure of Bacillus subtilis RibG: THE REDUCTION MECHANISM DURING RIBOFLAVIN BIOSYNTHESIS
J. Biol. Chem., January 16, 2009; 284(3): 1725 - 1731.
[Abstract] [Full Text] [PDF]


Home page
Proc. Natl. Acad. Sci. USAHome page
T. Toyomasu, M. Tsukahara, A. Kaneko, R. Niida, W. Mitsuhashi, T. Dairi, N. Kato, and T. Sassa
Fusicoccins are biosynthesized by an unusual chimera diterpene synthase in fungi
PNAS, February 27, 2007; 104(9): 3084 - 3088.
[Abstract] [Full Text] [PDF]


This Article
Right arrow Abstract Freely available
Right arrow Full Text (PDF)
Right arrow All Versions of this Article:
281/11/7605    most recent
M510254200v1
Right arrow Submit a Letter to Editor
Right arrow Alert me when this article is cited
Right arrow Alert me when eLetters are posted
Right arrow Alert me if a correction is posted
Right arrow Citation Map
Services
Right arrow Email this article to a friend
Right arrow Similar articles in this journal
Right arrow Similar articles in PubMed
Right arrow Alert me to new issues of the journal
Right arrow Download to citation manager
Right arrowRequest Permissions
Citing Articles
Right arrow Citing Articles via HighWire
Right arrow Citing Articles via Google Scholar
Google Scholar
Right arrow Articles by Chen, S.-C.
Right arrow Articles by Liaw, S.-H.
Right arrow Search for Related Content
PubMed
Right arrow PubMed Citation
Right arrow Articles by Chen, S.-C.
Right arrow Articles by Liaw, S.-H.
Social Bookmarking
 Add to CiteULike   Add to Complore   Add to Connotea   Add to Del.icio.us   Add to Digg   Add to Reddit   Add to Technorati  
What's this?


HOME HELP FEEDBACK SUBSCRIPTIONS ARCHIVE SEARCH TABLE OF CONTENTS
 All ASBMB Journals   Molecular and Cellular Proteomics 
 Journal of Lipid Research   ASBMB Today 
Copyright © 2006 by the American Society for Biochemistry and Molecular Biology.
Advertisement
spacer
Advertisement
Advertisement