Structural analysis of FAD synthetase from Corynebacterium ammoniagenes
© Frago et al; licensee BioMed Central Ltd. 2008
Received: 07 May 2008
Accepted: 23 September 2008
Published: 23 September 2008
The prokaryotic FAD synthetase family – a group of bifunctional enzymes that catalyse riboflavin phosphorylation and FMN adenylylation within a single polypeptide chain- was analysed in terms of sequence and structure.
Sequences of nearly 800 prokaryotic species were aligned. Those related with bifunctional FAD synthetase activities showed conservation of several consensus regions and highly conserved residues. A 3D model for the FAD synthetase from Corynebacterium ammoniagenes (CaFADS) was generated. This model confirms that the N-terminal and C-terminal domains are related to nucleotydyltransferases and riboflavin kinases, respectively. Models for the interaction of CaFADS with its substrates were also produced, allowing location of all the protein substrates in their putative binding pockets. These include two independent flavin binding sites for each CaFADS activity.
For the first time, the putative presence of a flavin binding site for the adenylylation activity, independent from that related with the phosphorylation activity, is shown. Additionally, these models suggest the functional relevance of some residues putatively involved in the catalytic processes. Their relevant roles were analysed by site-directed mutagenesis. A role was confirmed for H28, H31, S164 and T165 in the stabilisation of the P groups and the adenine moiety of ATP and, the P of FMN for the adenylylation. Similarly, T208, N210 and E268 appear critical for accommodation of the P groups of ATP and the ribityl end of RF in the active site for the phosphorylation process. Finally, the C-terminal domain was shown to catalyse the phosphorylation process on its own, but no reaction at all was observed with the individually expressed N-terminal domain.
The activities of riboflavin (RF) phosphorylation and FMN adenylylation are present among all kingdoms of living organisms. However, whereas eukaryotes use two different enzymes for FMN (flavokinase) and FAD (adenylyltransferase) production [1–6], prokaryotic organisms depend on an enzyme of ~38 kDa, FAD synthetase (Flavin adenine dinucleotide synthetase, FADS), that exhibits both activities . Thus, FADS constitutes a bifunctional prokaryotic protein family of enzymes with both, ATP:RF 5'-phosphotransferase [Riboflavin kinase (RFK) (EC 184.108.40.206) (Flavokinase)] and ATP:FMN-adenylyltransferase [FMN adenylyltransferase (EC 220.127.116.11)] activities [7, 8] that catalyses the 5'-phosphorylation of RF to FMN and, subsequently, the adenylylation of FMN to FAD.
RF phosphorylation RF + ATP → FMN + ADP
FMN adenylylation FMN + ATP → FAD + PPi
The FMN adenylylation reaction of FADS is a reversible process, while phosphorylation of RF appears to be irreversible . Additionally, FADS has been shown to function with a broad variety of RF isoesters, and it is, therefore, widely used in the preparation of FMN and FAD analogues [8, 10]. A kinetic mechanism for the steady-state reaction of the FADS from Corynebacterium ammoniagenes (CaFADS) has been proposed , suggesting an order in substrate binding and product release. The FMN intermediate produced by the first reaction appears to be released by the enzyme to later rebind the enzyme as a substrate for FAD production. Both enzyme activities present important differences in their substrate requirements, including concentration dependence, specificity for divalent cations and optimal pH or temperature [7, 11, 12].
The C-terminal region of FADS (~150 residues) shares considerably sequence similarity with eukaryotic monofunctional RFKs. Likewise, the N-terminal region has been proposed to present remote similarity to nucleotydyltransferases (NTs), suggesting its involvement in the adenylylation process . The structure of bifunctional FADS has only been reported for the Thermotoga maritima enzyme (TmFADS), both free and in complex with some substrates [14, 15]. This structure shows that the enzyme is folded in two domains and it also confirms the presence of one ATP-binding site in each of the domains and a single flavin-binding site [14, 15]. So far, neither functional studies of the TmFADS nor FADS structures from any other species have been reported.
A sequence analysis of the bifunctional FADS family and an in silico structural model for the CaFADS are here presented. Based on these data, a preliminary functional analysis of the CaFADS has been performed. Several residues are shown to be exclusively involved in either the first or second catalytic event. Additionally, the C- and N-terminal domains have been independently produced and assayed for their activities. The structural information here provided might help to envisage the rational design of selective antimicrobial drugs with the function of inhibiting FMN or FAD production and therefore, the availability in the cell of enzymes depending on these cofactors.
Sequence searches and alignments
Prediction of a three-dimensional model for the structure of FAD synthetase from C. ammoniagenes
The structural model for CaFADS was constructed using homology modelling procedures based on multiple structure-base sequence alignments (including all the proteins in the PDB) as implemented in the Geno3D Web server . This program uses distance geometry, simulated annealing and energy minimisation algorithms to build the protein models. The structural quality of the models was checked using the PROCHECK validation program . The structures of TmFADS (1mrz, 1s4m, 1t6x, 1t6y) [14, 15] and those of RFKs from H. sapiens (1q9s, 1p4m) [21, 22] and S. pombe (1n06, 1n08)  were suitable as templates. Six different combinations of templates were used and three models were generated by each combination.
Prediction of substrate binding sites
Structures used to model CaFADS and the interaction with its ligands.
(%) CONSERVATIVE SUBSTITUTIONS
T. maritima complexed
FMN, ADP, AMP
GCT B. subtilis
NMNAT M. jannaschii
NMNAT H. sapiens
NMNAT M. thermoautotropicum
PPAT E. coli
PPAT T. thermophilus
ADP, Mg2+, FMN
ADP, Mg2+, FMN
Cloning and production of C. ammoniagenes wild-type and mutated FAD synthetase forms in E. coli
Following the Ethical Rules established at the Universidad de Zaragoza, strain DSM 20305 (ATCC 6872) from C. ammoniagenes was grown at 30°C in a culture medium containing 0.5% glucose, 0.5% NaCl, 1% tryptone and 0.5% yeast extract, pH 7.2–7.4. Genomic DNA was isolated from the cells as previously described . The CaFADS encoding sequence was then amplified by the polymerase chain reaction using the forward 5'-GTAAGCCATGGATATTTGGTACGGAACAG-3' and reverse 5'-CCATCGATAGCGGATCCGGCATATAC-3' primers, synthesised on the basis of the published nucleotide sequence of the enzyme . The primers incorporate the NcoI and BamHI sites respectively, and have the transcription start codon GTG replaced with ATG. The amplified 1237 bp DNA fragment was then cloned in the NcoI/BamHI sites of the pET-28a(+) expression vector. The resulting recombinant plasmid, PET28a-FADS, was used to transform E. coli BL21(DE3) competent cells. The transformants were selected on solid LB/agar medium containing 30 μg/ml kanamicin. The fragment was sequenced to confirm the accuracy of the amplification. Mutations were introduced into the cloned structural gene encoding the WT CaFADS by site-directed mutagenesis. The H28A, H28D, H31A, H31D, R161A, R161D, S164A, S164D, T165A, T165D, T208A, T208D, N210A, N210D, E268A and E268D FADS mutants were generated using the QuikChange mutagenesis kit (Stratagene) in combination with the adequate synthetic oligonucleotides. The gene codifying for the individual N-terminal domain (PET28a-Δ(183–338)FADS) was obtained replacing the codon codifying for Leu183 in PET28a-FADS with a stop codon, following the methodology above mentioned for production of site-directed mutants. The C-terminal domain was individually cloned (PET28a-Δ(1–182)FADS) using the protocol described for WT CaFADS with the forward primer 5'-CCAACTGGGCCATGGGGCGGCAC-3' and the reverse primer above used for the full length protein cloning. The forward primer used to clone the C-terminal domain incorporates the start codon ATG in a NcoI site substituting for Leu183. Mutations were verified by DNA sequence analysis. WT and mutant FADS proteins were over-expressed in E. coli.
Purification of recombinant FADS from C. ammoniagenes
E. coli cells containing the recombinant pET28a-FADS plasmid (or its mutants) were grown in LB. During the exponential E. coli growing phase, FADS was induced by overnight incubation at 37°C with IPTG (1 mM). Cells were harvested by centrifugation and stored at -20°C. In a typical purification, around 25 g of cells (over-expressing recombinant FADS) were thawed, resuspended in 125 ml of cell disruption buffer (50 mM Tris/HCl, pH 8.0, 1 mM EDTA, 12 mM β-mercaptoethanol, 1 μM PMSF) and broken by ultrasonic treatment at 4°C using 16 cycles of 1 min in a DRH UP200 DR Hielsher sonicator. The cell debris was removed by centrifugation at 40000 g during 45 min. The resulting yellow FADS-containing supernatant was fractioned with 45% ammonium sulphate. The mixture was again centrifuged for 45 min at 4000 g and the supernatant was loaded onto a Phenyl-Sepharose High performance (Amersham Biosciences, GE Healthcare) column equilibrated with 50 mM Tris/HCl, pH 8.0, 45% ammonium sulphate. The column was washed with the same buffer containing only 17% ammonium sulphate until most of the yellow colour washed out of the column. The enzyme was eluted using a 17→0% ammonium sulphate reversed-gradient in the same buffer. Fractions of 6–8 ml were collected while the absorbance at 280 nm was recorded. Fractions presenting absorbance at 280 nm were analysed by SDS-PAGE and, those containing the over-expressed FADS were pooled and dialysed against 50 mM Tris/HCl, pH 8.0. After dialysis the protein was loaded onto a DEAE-cellulose (DE52, Whatman, England) column equilibrated with 50 mM Tris/HCl, pH 8.0. The column was washed with the equilibration buffer containing 0.1 M NaCl and subsequently eluted using a 0.1→0.5 M NaCl gradient. Fractions of 5 ml were collected while recording the absorbance at 280 nm and, subsequently analysed by SDS-PAGE. The enzyme was considered pure when a single band of the expected size was observed in SDS-PAGE. Pure FADS fractions were pooled together, dialysed against 50 mM Tris/HCl, pH 8.0 and stored at -20°C.
Measurement of riboflavin kinase and adenylyltransferase activities of FAD synthetase
Conversion of RF into FMN and, of FMN into FAD were qualitatively assayed by addition of the different FADS variants (final enzyme concentration ~1 μM, calculated using a theoretical ε280 nm = 27.8 mM-1cm-1) to a solution (final volume, 150 μl) containing 50 μM flavin (either RF or FMN), 0.2 mM ATP and either 0.8 mM or 10 mM MgCl2, in 50 mM Tris/HCl, pH 8.0. All these chemicals were obtained from SIGMA and were from the highest purity available. After 30 minutes of incubation at 37°C, the reaction was stopped by boiling the preparations for 5 minutes. Transformation of RF into FMN and FAD, or FMN into FAD, was visualised by resolving the products of the reaction at room temperature and in the dark by TLC on Silica Gel SIL-G-25 (20 cm × 20 cm, thickness 0.25 mm) plates . The moving phase was a solution of butanol:acetic acid:water (12:3:5). Flavin TLC spots were visually examined and scanned by determining their fluorescence under an ultraviolet light. The different flavins move up the plate at different rates due to differences in their partioning behaviour between the mobile liquid phase and the stationary phase. The weaker intensity of the FAD fluorescence with regard to equivalent amounts of RF or FMN relates with its flavin fluorescence quenching upon stacking of the adenine portion . An additional band, running slightly faster than FMN, is observed in some experiments. The nature of this band is not known, but it is only present in the experiments containing simultaneously FMN and FAD.
Results and discussion
Sequence analysis of the FAD synthetase family
No matches for bifunctional FADS sequences were obtained when performing the SIB-BLAST search of the CaFADS sequence on the UniProt Knowledgebase eukaryota subsection. When searching the bacteria+archaea database subsection (2,929,087 sequences searched), the 777 sequences displaying E-values ≤ 0.009 were chosen for further analysis. Redundant sequences from the same species presenting sequence identity over 80% were eliminated. The remaining sequences included some proteins noticeably shorter as well as proteins lacking the typical consensus sequences found in FADSs. These sequences (less than 5%), most of them displaying the lowest E-values of the selected sequences, were removed and separately analysed. They could be sorted out into several groups: i) proteins that present only the N-terminal region of FADS, ii) proteins that present only the C-terminal region of FADS, iii) proteins that lack the first conserved motif in the N-terminal region, iv) proteins with a conserved N-terminal region and a C-terminal region that shares no similarity with the bifunctional FADS family, and v) proteins with a conserved C-terminal region and a N-terminal region that shares no similarity with the bifunctional FADS family (in Rhodococcus sp. RHA1, the protein N-terminal region is 40% identical to 3,4-dihydroxy-2-butanone-4-phosphate synthase, an enzyme involved in RF biosynthesis). Some of these divergent proteins belong to organisms that also possess a typical prokaryotic FADS sequence.
After the manual refinement, a multiple sequence alignment over the remaining 500 FADS-like sequences, displaying E-values ≤ 2xe-6 with CaFADS, was carried out. Only FADSs from the genus Corynebacterium displayed sequence identities over 50% with the CaFADS (up to 60% identity in the case of the FADS from C. efficiens). However, the multiple alignment pointed out several highly conserved residues and motifs (Figure 1). Motifs 21-Φ-G-Ω-F-D-G-Φ-H-Ω-G-H-Ω-32 and 162-Φ-S-S-[TS]-Ω-[IV]-R-Ω-Ω-Φ-Ω-Ψ-G-174 (Φ, Ψ and Ω, denoting hydrophobic, polar and any residue, respectively) (N1 and N5 in Figure 1) corresponded to two consensus sequences in the FADS family that are also conserved in NTs remotely related to FADS (see below). Two additional consensus sequences in between the above mentioned motifs were also observed. They included residues 50-Ω-Φ-Φ-Ω-F-Ψ-P-[HQ]-P-Ω-59 and 120-Φ-Φ-Φ-G-Ω-[DN]-[FYH]-Ω-[FY]-G-Ω-129 (N2 and N4 in Figure 1) and, apparently, only the motif 123-G-Ω-Ψ-125 presents an equivalent in NT sequences (see below). Multiple sequence alignment also showed a number of residues highly conserved, or conservatively substituted, between residues 1 and 186: I71, R77, G85, I86, D87, F94, F98, Y106, V107, A179, L183 and G184 (CaFADS numbering). These residues are mainly included in region N3 and at the end of N5 in Figure 1 and, can be considered as forming additional consensus sequences. Noticeably, only R77 and A179 are conserved in NTs.
The C-terminal region of the bifunctional FADS sequences (residues 187–338) also included several consensus sequences, 191-G-Ω-V-Φ-Ψ-G-Ω-Ω-Ω-G-Ω-Ψ-Φ-203, 205-G-[FY]-P-T-[ALIV]-N-210, 223-P-Ω-Ω-G-[VI]-[YF]-225, 253-Ψ-Φ-G-Ω-Ψ-P-T-[FVI]-260, 267-[LVI]-E-Ω-[HFYN]-Φ-[FL]-[DN]-[FWY]-Ψ-Φ-[DNE]-[LIAV]-Y-[GDN]-280 and 291-[LI]-R-Ω-[EQMN]-Ω-[KRT]-[FY]-296, and the highly conserved residues L303 and D310 (see C6–C10 in Figure 1), being all of them also conserved in the sequences of RFKs (see below).
Structural model of the FAD synthetase from C. ammoniagenes
R.m.s.d. (Å) between the three predicted models for CaFADS and the structures used for their production.
T. maritima FADS (1mrz)
T. maritima FADS, RF (1s4m)
H. sapiens RFK, ADP, Mg2+, FMN (1p4m)
S. pombe RFK, ADP, Zn2+ (1n08)
According to the model, the protein folds in two domains. The N-terminal domain (residues 18–186) lacks the first 18 residues of the sequence (not present in the TmFADS template), while the C-terminal domain (residues 187–318) lacks the last 20 residues (templates were shorter in sequence or presented very low sequence similarity in those residues) (Figures 2 and 3). The main structural arrangements predicted for the CaFADS are present in the structure reported for the TmFADS [14, 15], but our model also showed the disposition of some loops and a 310-helix at the C-terminal domain that are missed in the TmFADS structure. Additionally, the C-terminal domain of CaFADS showed an insertion around residue 230 with regard to TmFADS and many members of the FADS family.
The kinase domain: the C-terminal domain of C. ammoniagenesFAD synthetase
The models for CaFADS suggest flexibility in the loops connecting β6–β7 and β9–β10 (Figure 2B). Both of these regions are not solved in the structures reported for the TmFADS [14, 15]. Flexibility is also observed in regions involved in substrate binding (Figure 4B), in agreement with conformational changes observed in RFKs upon substrate interaction . These conformational changes provide additional interactions, for both the ribityl tail and the isoalloxazine ring and, help to anchor the reactive substrates and the catalytic residues in the most adequate disposition for catalysis [21–23]. This analysis also suggests that despite RFKs and FADSs surely share a similar mechanism in the RF phosphorylation reaction, variations in the structural details for substrate binding and catalysis can be expected.
The adenylyltransferase domain: the N-terminal domain of C. ammoniagenesFAD synthetase
Flexible regions can be predicted by comparison of the CaFADS models (Figure 2B, Table 2). These, including 55–67, 85–97, 111–119 and 127–135, correspond to low similarity regions between NT and FADS families (Figures 2B, 3A and 5A). Structures of NMNAT in complex with NAD+ , and PPAT in complex with dPCoA  or Ppant , indicate that the binding of the adenylylation substrate, NMN+ or Ppant moieties, is allocated in the pocket formed by the above mentioned regions. Superposition of these complexes on our CaFADS model, followed by fitting and sculpting the FMN molecule along the nicotinamide nucleotide moiety of NAD+ in the NMNAT structure, allowed production of putative models for the interaction of the FADS N-terminal domain with its FMN substrate (Figure 6A). In this model, the FMN isoalloxazine ring would be placed in a cavity surrounded by the above mentioned flexible regions. The entrance of this cavity would be flanked by the F54, F62 and D94 hydrophobic side-chains, the consensus region 54-FPHP-58 would fix the position of the ribityl of FMN, while the flavin ring environment would be contributed by hydrophobic interactions with V59, L98, Y106 and F128 (Figure 6A). All these positions are highly conserved in the FADS family, but do not have equivalence in NTs (Figures 1 and 3A). The model suggests that the hydrophobic isoalloxazine portion of the flavin ring might be embedded inside the cavity, while the hydrophilic edge would be partially exposed to the solvent. The size of the cavity makes it difficult to predict the precise location of the isoalloxazine ring. Various possibilities are allowed (Figure 6A) and, final disposition of loops and side-chains would depend on the flavin binding event and conformational changes not predicted in our model. This flavin binding mode could represent a novel flavin binding motif so far not described. Additionally, while sequence motifs N1, N4 and N5 are also conserved in the sequences of eukaryotic FMN adenylyltransferases, residues and regions proposed for flavin binding in the FADS family do not present homologous in the sequences of the monofunctional enzymes (data not shown).
This is the first structural model where a putative binding site for the flavin is suggested at the adenylylation site of FADS, indicating the presence of two independent flavin binding sites in the FADS family. An early hypothetical model for the active site of FADS proposed a single flavin binding site for both activities and two independent ATP binding sites . No structural information of this family of enzymes was reported at that time, but that model accounted for the kinetic mechanism experimentally derived . The 3D structures and models today available for FADS indicate that the phosphorylation and adenylylation sites are far away (Figure 6B). This fact is in agreement with experimental observations indicating that the FMN produced in the phosphorylation process was released from the protein, to subsequently rebind to initiate the adenylylation process . Superposition of our CaFADS model in complex with all its substrates with the dimer reported in the TmFADS crystallographic asymmetric unit also indicates that phosphorylation and adenylylation active sites from different monomers are far away to share the FMN molecule. However, this dimeric model suggests that the 229–243 loop and the terminal α7 of the C-terminal domain might contribute to the formation of the ATP and the flavin binding pockets in the adenylylation domain (Figure 6B).
Production of C. ammoniagenesFAD synthetase forms
E. coli BL21(DE3) cells transformed with the recombinant PET28a-FADS vector produced a high level expression of active CaFADS. This CaFADS was purified in complex with one of the products of its reaction, FAD, as confirmed by HPLC (data not shown). To remove the FAD a hydrophobic chromatography was required. The sequence and structural analysis suggested the presence of several residues and motifs in CaFADS highly conserved in this family of enzymes and potentially involved in the Catalytic sites of the RFK and NT activities. Some of these residues were mutated, namely H28, H31, R161, S164, T165, T208, N210 and E268 to Ala and Asp. These CaFADS variants were expressed in E. coli, with yields similar to the WT FADS, and were purified to homogeneity following the same protocol. Additionally, plasmids individually expressing the N-terminal and C-terminal CaFADS domains were also produced, PET28a-Δ(183–338)FADS and PET28a-Δ(1–182)FADS, respectively, and transformed in E. coli, showing a level of expression similar to the full length protein.
C. ammoniagenesFAD synthetase residues involved in FMN and FAD production
Replacement of H28 with either Ala or Asp and H31 with Ala prevented conversion of FMN into FAD, confirming their implication in the CaFADS NT activity. Both of these His residues belong to a consensus region 27-Φ-H-Ω-G-H-Ω-32 found both in FADS and NT families (Figure 1 and 3A). Therefore, H28 and H31 will surely be major determinants in the stabilisation of the P groups and the adenine moiety of ATP in the reaction transition state in CaFADS, in agreement with previous observations in NTs [31, 32, 34]. Replacement of H28 with Ala also prevented the conversion of RF into FMN, pointing to a possible influence of the N-terminal domain in the RFK activity. His28 appears situated at the N-terminal of α1, in contact with α5 and α6, which, according to our model, will interact with the C-terminal domain from the same polypeptide chain, but also might influence the C-terminal of the second monomer in a putative dimeric structure (Figure 6B). However, it is difficult to envisage how the H28A mutation might affect the phosphorylation reaction with the available information. This will need further structural and biochemical characterisation. Replacement of H31 with Asp produced an enzyme still able to produce FMN and FAD, suggesting that an Asp at this position somewhat replaces the function of the His. Replacements were also introduced in the second major consensus sequence shared by the FADS and the NT families, 162-ΦSST-165. Replacement of S164 or T165 with Ala produced enzyme forms able to transform RF into FMN and FMN into FAD, while the introduction of a Asp residue at such positions clearly prevented FAD formation. This latter observation was consistent with these residues being involved in the stabilisation at the NT catalytic site of the negatively charged β-P and γ-P of ATP or the ribityl end, as suggested by the above structural analysis (Figures 3A and 5B). Finally, the putative role of the structural disposition of the positive charge of R161 in the stabilisation of the adenine and α-P and β-P of ATP in the active site was studied by replacement of R161 with either Ala or Asp. The produced FADS mutants were still able to produce FMN and FAD under the assayed conditions, suggesting the positive charge at position of R161 is not playing a critical role in CaFADS catalysis.
Taking into account the fact that T208, N210, and E268 are conserved in sequence and structural position in the close environment of the ATP, not only in the FADS family but also in RFKs, and that the terminal oxygen of the ribityl of RF binding sites (Figures 1, 3B and 4), a catalytic role for them in the RF phosphorylation process can be envisaged. These CaFADS residues were replaced with Ala or Asp to prove such a role. All these mutations prevented FMN production, but allowed transformation of FMN into FAD, confirming their implication in the RFK activity of FADS. These data are, therefore, in agreement with the proposed structural-based role for T208 and N210 in accommodating the active portion of the ATP molecule, the P groups, and the metal in the catalytic centre (Figure 4B). The lack of FMN and FAD production when replacing E268 with Ala confirms it is a key residue in the CaFADS RFK activity. Analysis of the interaction of the CaFADS model with FMN (or RF) also indicated that in CaFADS E268 side-chain might H-bond to the terminal OH of the ribityl chain of RF (Figure 4B). This Glu is invariant in the FADSs and RFKs families (Figures 1 and 3B). In this later family this residue has been proposed to act as a catalytic base [21–23]. Moreover, the fact that incubation of E268D with RF did not produce any accumulation of FMN, while traces of FAD were detected, is consistent with the RFK activity taking place at a very slow rate in this mutant, thus avoiding FMN accumulation. This suggests that although the introduced mutation is conservative and the Asp side-chain might still act as a base, its shorter side-chain prevents the optimal organisation of the transition state in the active centre and, therefore, reduces the catalytic efficiency of the RFK activity. Again, this is consistent with this residue not presenting substitutions, even conservative ones, in the analysed proteins. Nevertheless, other explanations might be also possible, as higher catalytic efficiency of NT activity than that of RFK or sensitivity of the TLC analysis. Further work should be done by the analysis of the kinetic parameters for both activities, as well as of the interaction parameters for the different substrates, for all these CaFADS mutants to provide additional insights in their specific role in the enzyme catalytic mechanisms.
Activity of the independently produced C. ammoniagenesFAD synthetase N- and C- terminal domains
Crude extracts of the E. coli cells expressing the individually cloned N-terminal and C-terminal domains of CaFADS, Δ(183–338)FADS and Δ(1–182)FADS, respectively, were prepared and assayed for NT and RFK activities (Figures 7B, 7C). Analysis of the products of the reaction by TLC concluded that the C-terminal domain, Δ(1–182)FADS, was able to qualitatively catalyse the phosphorylation process. However, crude extracts over-expressing the N-terminal domain, Δ(183–338)FADS, were not able to catalyse the conversion of FMN into FAD. Extracts containing both domains were mixed and assayed for both reactions. Conversion of RF into FMN was produced in these samples, but they were not able to mimic CaFADS in the adenylylation process. This is in agreement with a possible C-terminal domain role in contributing to the stabilisation of the flavin and ATP substrates in the corresponding N-terminal domain pockets (Figure 6B).
A sequence analysis has been carried out in the bifunctional FADS family. An important number of consensus regions and conserved residues have been determined, despite the low sequence identity among members of this family. Sequence divergences found among pathogenic organisms, such as Mycobacterium or Corynebacterium, and monofunctional eukaryotic enzymes will probably be reflected in structural divergences that might affect the catalytic mechanism, making possible the design of specific antimicrobial drugs. Models have been proposed, not only for the CaFADS structure, but also for the complexes formed with all its substrates. This includes de interaction with the FMN substrate at the adenylylation site, in a region not previously reported as involved in the enzyme activity and, in a pocket that might constitute a novel flavin binding motif. As expected, these models supported that the N-terminal presents similar folding to NTs and is involved in the adenylylation of FMN, whereas the C-terminal is similar to RFKs and responsible for the phosphorylation. Models of CaFADS allowed identification of several residues putatively involved in the catalysis. Site-directed mutagenesis on some of these residues confirmed their participation in substrate stabilisation and/or the catalytic mechanism. The particular behaviour of some of these mutants and of the independently expressed NT and RFK domains suggests one of the domains might influence the activity in the other one. This study sheds important and new information on the structural mechanism of substrate recognition and catalysis while waiting for the experimental structure of the CaFADS and provides a platform for future investigations into the mechanism of enzymes of the bifunctional FADS family.
This work has been supported by the Spanish Ministry of Sciences (Grant BIO2007-65890-C02-01 to M.M.).
- Barile M, Brizio C, Valenti D, De Virgilio C, Passarella S: The riboflavin/FAD cycle in rat liver mitochondria. Eur J Biochem. 2000, 267 (15): 4888-4900.View ArticlePubMedGoogle Scholar
- McCormick DB, Oka M, Bowers-Komro DM, Yamada Y, Hartman HA: Purification and properties of FAD synthetase from liver. Methods Enzymol. 1997, 280: 407-413.View ArticlePubMedGoogle Scholar
- Oka M, McCormick DB: Complete purification and general characterization of FAD synthetase from rat liver. J Biol Chem. 1987, 262 (15): 7418-7422.PubMedGoogle Scholar
- Merrill AH, McCormick DB: Affinity chromatographic purification and properties of flavokinase (ATP:riboflavin 5'-phosphotransferase) from rat liver. J Biol Chem. 1980, 255 (4): 1335-1338.PubMedGoogle Scholar
- Brizio C, Galluccio M, Wait R, Torchetti EM, Bafunno V, Accardi R, Gianazza E, Indiveri C, Barile M: Over-expression in Escherichia coli and characterization of two recombinant isoforms of human FAD synthetase. Biochem Biophys Res Commun. 2006, 344 (3): 1008-1016.View ArticlePubMedGoogle Scholar
- Galluccio M, Brizio C, Torchetti EM, Ferranti P, Gianazza E, Indiveri C, Barile M: Over-expression in Escherichia coli, purification and characterization of isoform 2 of human FAD synthetase. Protein Expr Purif. 2007, 52 (1): 175-181.View ArticlePubMedGoogle Scholar
- Manstein DJ, Pai EF: Purification and characterization of FAD synthetase from Brevibacterium ammoniagenes. J Biol Chem. 1986, 261 (34): 16169-16173.PubMedGoogle Scholar
- Grill S, Busenbender S, Pfeiffer M, Kohler U, Mack M: The bifunctional flavokinase/flavin adenine dinucleotide synthetase from Streptomyces davawensis produces inactive flavin cofactors and is not involved in resistance to the antibiotic roseoflavin. J Bacteriol. 2008, 190 (5): 1546-1553.PubMed CentralView ArticlePubMedGoogle Scholar
- Efimov I, Kuusk V, Zhang X, McIntire WS: Proposed steady-state kinetic mechanism for Corynebacterium ammoniagenes FAD synthetase produced by Escherichia coli. Biochemistry. 1998, 37 (27): 9716-9723.View ArticlePubMedGoogle Scholar
- Murthy YV, Massey V: Syntheses and applications of flavin analogs as active site probes for flavoproteins. Methods Enzymol. 1997, 280: 436-460.View ArticlePubMedGoogle Scholar
- Hagihara T, Fujio T, Aisaka K: Cloning of FAD synthetase gene from Corynebacterium ammoniagenes and its application to FAD and FMN production. Appl Microbiol Biotechnol. 1995, 42 (5): 724-729.View ArticlePubMedGoogle Scholar
- Nakagawa S, Igarashi A, Ohta T, Hagihara T, Fujio T, Aisaka K: Nucleotide sequence of the FAD synthetase gene from Corynebacterium ammoniagenes and its expression in Escherichia coli. Biosci Biotechnol Biochem. 1995, 59 (4): 694-702.View ArticlePubMedGoogle Scholar
- Krupa A, Sandhya K, Srinivasan N, Jonnalagadda S: A conserved domain in prokaryotic bifunctional FAD synthetases can potentially catalyze nucleotide transfer. Trends Biochem Sci. 2003, 28 (1): 9-12.View ArticlePubMedGoogle Scholar
- Wang W, Kim R, Jancarik J, Yokota H, Kim SH: Crystal structure of a flavin-binding protein from Thermotoga maritima. Proteins. 2003, 52 (4): 633-635.View ArticlePubMedGoogle Scholar
- Wang W, Kim R, Yokota H, Kim SH: Crystal structure of flavin binding to FAD synthetase of Thermotoga maritima. Proteins. 2005, 58 (1): 246-248.View ArticlePubMedGoogle Scholar
- Altschul SF, Madden TL, Schaffer AA, Zhang J, Zhang Z, Miller W, Lipman DJ: Gapped BLAST and PSI-BLAST: a new generation of protein database search programs. Nucleic Acids Res. 1997, 25 (17): 3389-3402.PubMed CentralView ArticlePubMedGoogle Scholar
- Thompson JD, Higgins DG, Gibson TJ: CLUSTAL W: improving the sensitivity of progressive multiple sequence alignment through sequence weighting, position-specific gap penalties and weight matrix choice. Nucleic Acids Res. 1994, 22 (22): 4673-4680.PubMed CentralView ArticlePubMedGoogle Scholar
- Crooks GE, Hon G, Chandonia JM, Brenner SE: WebLogo: a sequence logo generator. Genome Res. 2004, 14 (6): 1188-1190.PubMed CentralView ArticlePubMedGoogle Scholar
- Combet C, Jambon M, Deleage G, Geourjon C: Geno3D: automatic comparative molecular modelling of protein. Bioinformatics. 2002, 18 (1): 213-214.View ArticlePubMedGoogle Scholar
- Laskowski RA, MacArthur MW, Moss DS, Thornton JM: PROCHECK: a program to check the stereochemical quality of protein structures. J Appl Cryst. 1993, 26: 283-291.View ArticleGoogle Scholar
- Karthikeyan S, Zhou Q, Mseeh F, Grishin NV, Osterman AL, Zhang H: Crystal structure of human riboflavin kinase reveals a beta barrel fold and a novel active site arch. Structure. 2003, 11 (3): 265-273.View ArticlePubMedGoogle Scholar
- Karthikeyan S, Zhou Q, Osterman AL, Zhang H: Ligand binding-induced conformational changes in riboflavin kinase: structural basis for the ordered mechanism. Biochemistry. 2003, 42 (43): 12532-12538.View ArticlePubMedGoogle Scholar
- Bauer S, Kemter K, Bacher A, Huber R, Fischer M, Steinbacher S: Crystal structure of Schizosaccharomyces pombe riboflavin kinase reveals a novel ATP and riboflavin-binding fold. J Mol Biol. 2003, 326 (5): 1463-1473.View ArticlePubMedGoogle Scholar
- Jones TA, Zou JY, Cowan SW, Kjeldgaard M: Improved methods for building protein models in electron density maps and the location of errors in these models. Acta Crystallogr A. 1991, 47 (Pt 2): 110-119.View ArticlePubMedGoogle Scholar
- Guex N, Peitsch MC: SWISS-MODEL and the Swiss-PdbViewer: an environment for comparative protein modeling. Electrophoresis. 1997, 18 (15): 2714-2723.View ArticlePubMedGoogle Scholar
- Delano WL: They PyMOL molecular graphics system. 2002, DeLano Scientific, San Carlos, CA, USA, [http://www.pymol.org]Google Scholar
- Pospiech A, Neumann B: A versatile quick-prep of genomic DNA from Gram-positive bacteria. Trends in Genetics. 1995, 11 (6): 217-218.View ArticlePubMedGoogle Scholar
- Bauer K, Gros L, Sauer W: Thin Layer Chromatography: An Introduction. 1991, Germany: Hüthig VerlagGoogle Scholar
- Berg van den PAW, Feenstra KA, Mark AE, Berendsen HJC, Visser AJWG: Dynamic conformations of flavin adenine dinucleotide: Simulated molecular dynamics of the flavin cofactor related to the time-resolved fluorescence characteristics. J Phys Chem B. 2002, 106 (34): 8858-8869.View ArticleGoogle Scholar
- Weber CH, Park YS, Sanker S, Kent C, Ludwig ML: A prototypical cytidylyltransferase: CTP:glycerol-3-phosphate cytidylyltransferase from Bacillus subtilis. Structure. 1999, 7 (9): 1113-1124.View ArticlePubMedGoogle Scholar
- Garavaglia S, D'Angelo I, Emanuelli M, Carnevali F, Pierella F, Magni G, Rizzi M: Structure of human NMN adenylyltransferase. A key nuclear enzyme for NAD+ homeostasis. J Biol Chem. 2002, 277 (10): 8524-8530.View ArticlePubMedGoogle Scholar
- D'Angelo I, Raffaelli N, Dabusti V, Lorenzi T, Magni G, Rizzi M: Structure of nicotinamide mononucleotide adenylyltransferase: a key enzyme in NAD+ biosynthesis. Structure. 2000, 8 (9): 993-1004.View ArticlePubMedGoogle Scholar
- Saridakis V, Christendat D, Kimber MS, Dharamsi A, Edwards AM, Pai EF: Insights into ligand binding and catalysis of a central step in NAD+ synthesis: structures of Methanobacterium thermoautotrophicum NMN adenylyltransferase complexes. J Biol Chem. 2001, 276 (10): 7225-7232.View ArticlePubMedGoogle Scholar
- Izard T, Geerlof A: The crystal structure of a novel bacterial adenylyltransferase reveals half of sites reactivity. Embo J. 1999, 18 (8): 2021-2030.PubMed CentralView ArticlePubMedGoogle Scholar
- Takahashi H, Inagaki E, Fujimoto Y, Kuroishi C, Nodake Y, Nakamura Y, Arisaka F, Yutani K, Kuramitsu S, Yokoyama S, et al: Structure and implications for the thermal stability of phosphopantetheine adenylyltransferase from Thermus thermophilus. Acta Crystallogr D Biol Crystallogr. 2004, 60 (Pt 1): 97-104.View ArticlePubMedGoogle Scholar
- Mizuguchi KDC, Blundell TL, Johnson MS, Overington JP: JOY: protein sequence-structure representation and analysis. Bioinformatics. 1998, 14: 617-623.View ArticlePubMedGoogle Scholar
- Izard T: The crystal structures of phosphopantetheine adenylyltransferase with bound substrates reveal the enzyme's catalytic mechanism. J Mol Biol. 2002, 315 (4): 487-495.View ArticlePubMedGoogle Scholar
This article is published under license to BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.