|
|
||||||||
Division of Microbiology, National Center for Toxicological Research, US Food and Drug Administration, Jefferson, AR 72079-9502, USA
Correspondence
Huizhong Chen
huizhong.chen{at}fda.hhs.gov
| ABSTRACT |
|---|
|
|
|---|
| INTRODUCTION |
|---|
|
|
|---|
Azoreductase from Ent. faecalis (AzoA) is one of the most active oxygen-tolerant FMN-dependent azoreductases characterized to date, showing a broad spectrum of substrate specificity. Recently, we reported the 3D structure of AzoA determined at 2.07 Å (0.207 nm) resolution (Liu et al., 2007
). That study provided direct evidence for a 46 kDa asymmetrical homodimer with one FMN molecule bound to each monomer. Dimeric AzoA has two separate active sites that are located at the interfaces between the two monomers, and the binding site of FMN lies inside each active site, indicating that the dimeric arrangement is required for active holoenzyme. In the active site, the FMN sits inside a positively charged cavity formed by several stretches of residues; the si face of the isoalloxazine ring of FMN is solvent-accessible, whereas its re face is buried in the protein. The structural environment of the active site could be correlated with the mechanical view adopted previously for azoreductase (Ito et al., 2006
; Liu et al., 2007
; Wang et al., 2007
).
AzoA absolutely requires reduced FMN for enzyme activity. Despite the large body of biochemical and 3D structural information available on bacterial FMN-dependent azoreductases (Chen et al., 2004
, 2005
; Ito et al., 2006
; Liu et al., 2007
; Wang et al., 2007
), little is known about the role of specific amino acid residues involved in the flavin binding and catalytic mechanism. In this study, a computational protein–FMN interaction analysis was conducted, which identified Trp-105 on loop
9 as a potential determinant for FMN stabilization. We utilized site-directed mutagenesis with the equivalents of Trp-105 from related proteins of known 3D structure and then performed kinetic analysis to evaluate the relative contributions of the side-chains on position 105 to the FMN binding and catalysis.
| METHODS |
|---|
|
|
|---|
Site-directed mutagenesis.
Site-directed mutagenesis was performed using the QuickChange II XL site-directed mutagenesis kit (Stratagene) according to the manufacturer's instructions. The mutagenic oligonucleotide primers are listed in Table 1
. The template was pAZOA, which consists of the azoA gene cloned into pET-11a (Chen et al., 2004
). Complete sequences were determined on a model 3100 Genetic Analyzer (Applied Biosystems) using T7 promoter and terminator primers to confirm that only the intended mutations had been introduced. pAZOA and its mutant plasmids were used for transformation of E. coli BL21-Gold (DE3)pLysS.
|
SDS-PAGE analysis.
SDS-PAGE was carried out in Laemmli's buffer (Laemmli, 1970
) with 12.5 % polyacrylamide. Perfect protein markers (Novagen) were used. Electrophoresis was performed in a Hoefer SE 260 Mighty Small II Mini Vertical unit (Amersham Biosciences). Gels were stained for proteins with Coomassie brilliant blue R-250 (Bio-Rad).
Enzyme assay of the mutant proteins.
Azoreductase activity was determined by measuring the decrease in A430 of MR with a Hewlett Packard 8453 UV–visible spectrophotometer at room temperature (23 °C). A typical reaction mixture (2.0 ml) contained 25 mM potassium phosphate buffer (pH 7.1), 25 µM MR and 0.1 mM NADH. The reaction was initiated by addition of the enzyme. Initial velocity was determined by monitoring the change in the amount of substrate in the first 2 min in a glass cuvette of 1.0 cm light path. The molar absorption coefficient for MR at 430 nm (
=23.4 mM–1 cm–1) was used. Proteins were quantified using the bicinchoninic acid assay (Pierce) with BSA as the standard (Chen et al., 2004
).
Kinetic constants of the mutant proteins.
Initial velocities of the enzymic reaction were determined by varying the concentrations of one substrate, MR (from 5 to 40 µM) or NADH (from 50 to 400 µM), while the concentration of the other substrate was kept constant (NADH 400 µM or MR 40 µM). Km and Vmax values were obtained from Lineweaver–Burk plots. Each data point was analysed three times.
Computational analysis
Homology modelling.
The homology models for mutants were generated using SWISS-MODEL (http://swissmodel.expasy.org/workspace/index.php) with the wild-type protein structure (2HPV) as a template. The structures were minimized by the GROMOS 43B1 force field; algorithm for minimization, steepest descent; maximum force to stop (kJ mol–1), 1.0; number of steps, 200 (Schwede et al., 2003
). Basic structure validation was checked using PROCHECK (European Bioinformatics Institute, Cambridge, UK).
In silico analysis.
To align the AzoA (2HPV) sequence with related proteins of known 3D structure in the PDB (Protein Data Bank), SAS (http://www.ebi.ac.uk/thornton-srv/databases/sas/) was used (Milburn et al., 1998
). The protein–FMN interactions were verified and plotted using the LPC (Ligand Protein Contacts) online server (http://bip.weizmann.ac.il/oca-bin/lpccsu) (Sobolev et al., 1999
) and LIGPLOT v.4.4.2 (Wallace et al., 1995
), respectively. SuperPose (version 1.0) was used for structural superposition and RMSD (root-mean-square deviation) calculation of the wild-type AzoA and the modelled mutant proteins (Maiti et al., 2004
). Two scoring functions, DrugScoreCSD (http://pc1664.pharmazie.uni-marburg.de/drugscore/) and the stand-alone X-Score (Velec et al., 2005
), were adopted to score FMN binding affinity in the protein–ligand complexes. PyMOL (0.99RC6) (http://www.pymol.org) was used to visualize the per-atom score contributions and all other structure figures. LPC and HBPLUS (McDonald & Thornton, 1994
) were used to make careful predictions of hydrogen bonds.
| RESULTS |
|---|
|
|
|---|
|
|
To apply structural information derived from the related proteins of known 3D structure in the PDB to the site-directed mutagenesis of Trp-105, a FASTA alignment from a search of the sequence of AzoA (2HPV) was generated by SAS analysis (Fig. 2
). The retrieved and aligned sequences from the PDB can be clustered by the ligand, FMN or FAD, into two groups. The proteins 2d5i (PDB) and 1t5b (PDB) each have a FMN ligand, and show a favourable sequence identity of 37.4 % of the entire protein sequence and appreciable alignment of secondary structure with AzoA. Note that those proteins with a FAD ligand, despite a relatively low sequence identity (<27.1 %) to AzoA, have appreciable overall alignments, based on the secondary structure of the flavin binding region. Human quinone oxidoreductase 2 (NQO2) has been known to reduce MR, indicating that it can also act as a four-electron transfer oxidoreductase like AzoA. NQO2 is expressed selectively in organs such as liver and kidney. It may function physiologically as one of several vitamin K reductases (Wu et al., 1997
). Moreover, although there is no significant conservation of residues in the flavin binding region, particularly in positions 103–107 of AzoA, the corresponding residues are in contact with the flavin in the structure of the proteins. In AzoA, Trp-105 is on the turning point of loop
9, made up of residues 104–110, which connects β-sheet β3 and
-sheet
3 (Fig. 2
). The equivalents of Trp-105 of AzoA in the related proteins of known 3D structure are either tyrosine or glutamine. Tyrosine is the most conserved among the FMN binding proteins in Fig. 2
.
|
2), suggesting that 2 FMN bind to 1 mutant AzoA (homodimer) (Table 3
|
|
|
|
carbons of the mutant proteins were less than 0.06 Å.
The FMN (oxidized form) binding affinity of each mutant was calculated using two different scoring functions, the knowledge-based DrugScore and the empirical-based X-score. The scores of the FMN binding affinity and the per-atom score contributions, together with the LPC analysis, can be used as a measure of the relative stability of the bound FMN in the structures. Specifically, the visualization of the per-atom score contributions gives us further insight into the overall structural impact of the mutations. The computational scores show a quantitative correlation with the experimentally determined FMN binding affinities. The analysis shows that the two substitutions with the smaller side-chains, W105A and W105G, cause the largest decrease in FMN binding affinity relative to that of the wild-type. In the case of mutant W105A, the score was decreased from 5.30 to 5.24 (wild-type vs W105A) in DrugScore. Strong unfavourable interactions are observed at position Met-104, where the N3-C2=O2 locus of FMN comes too close to the residue (Fig. 4
). The distance between M104-CE and FMN-N3 is 2.11 Å (5.20 Å in the wild-type). In LPC analysis, Gly-105 cannot contribute interactions for the stabilization of FMN (Table 3
). This mutation results in an inactive enzyme that was unable to bind to FMN, together with the W105A mutant.
The modelled structures also suggested that changing the side-chain size at position 105 could significantly affect the solvent exposure of the bound FMN cofactor. On LPC analysis, the exposed surface of the FMN systematically increased from 134.4 to 156 Å2 (wild-type vs W105G) as the size of the side-chain at position 105 was decreased, with the change in exposed area being localized mainly to the dimethylbenzene portion of the isoalloxazine ring.
| DISCUSSION |
|---|
|
|
|---|
Although azoreductases share a set of common features with flavodoxins with respect to the flavin binding, a significant difference is in the flavin isoalloxazine ring binding site. In flavodoxin, the isoalloxazine rings are sandwiched between aromatic amino acids, tryptophan and tyrosine, flanking the re face and si face, respectively (Lostao et al., 2003
; Murray & Swenson, 2003
). AzoA, however, does not have the combined net effect of Trp and Tyr on FMN binding. No aromatic residues seem to be in the si face directly affecting flavin ring binding (Fig. 1
). Tyr-129 is in a position at a distance of 5.34 Å from the OH of tyrosine to N5 of FMN and at a 6 ° angle to the si face of FMN, so that this aromatic residue cannot contribute to the flavin ring binding. Therefore, if AzoA Trp-105 is a homologue of the Trp involved in FMN stabilization of flavodoxin (Lostao et al., 1997
), the position of AzoA Trp-105 is intriguing. This aromatic amino acid residue appears to interact with all moieties of FMN, with its indole functional group being involved in the interactions. This contribution to FMN binding is probably due to the combination of the residue position and its side-chain specificity in AzoA. The presence of an aromatic side-chain at the flavin binding site is a common feature of many flavin binding proteins and is considered essential for the flavin cofactor binding and hydrophobic active site environment (Ito et al., 2006
; Liu et al., 2007
; Wang et al., 2007
). It is known to promote a stable and highly negative redox potential for the electron-transporting activity of the enzyme (Li & Meighen, 1995
).
The FMN binding affinity and enzyme activity of the mutants were revealed to be correlated closely with the side-chain length of the introduced amino acids. AzoA did not appear to be tolerant of the substitution of Trp-105 by short side-chain amino acids. Complete loss of flavin of the mutants W105A and W105G was probably due to a decrease in the number of FMN–protein interactions caused by the introduced short side-chain amino acids, resulting in no enzyme activity. On the other hand, substitution of Trp-105 by the aromatic amino acids Tyr and Phe did not affect the bound FMN content when compared to the wild-type AzoA (Table 3
). In the mutant W105Y, the free phenolic group of the Tyr-105 side-chain interacts with all the moieties of FMN, similar to Trp-105 of the wild-type AzoA. However, in spite of binding an approximately equal amount of FMN as compared to the wild-type, the mutants with bulky side-chains showed decreased catalytic activity. The results may be related to structural perturbations induced by the bulky side-chain substitutions. At least in part, local adjustments within the FMN binding site are likely to be responsible for the lower enzyme activities observed in the mutants. In fact, conformational changes could affect not only the binding of FMN, NADH and/or substrates like MR, but also the interactions among them. However, since structural changes were localized within the FMN binding area, it is most likely that the mutations caused changes in the position of the bound FMN, which negatively affected the interactions between FMN and NADH and/or substrate, resulting in the decrease in enzyme activity. Although none of the substitutions introduced significant changes to the mutants, relatively small changes in the protein structure within the immediate vicinity of the substitution were observed, leaving FMN possibly more accessible to solvent. The exposed surface of the FMN was gradually increased as the size of the side-chain at position 105 was decreased, with the change in exposed area being localized mainly in the dimethylbenzene portion of the isoalloxazine ring.
In conclusion, we performed this study to elucidate the structural and functional role of the residue Trp-105 in the AzoA reductase. The results demonstrated that Trp-105 plays important roles in both the FMN binding and catalytic ability of AzoA. This study provides a step forward in determining the catalytic properties of AzoA in FMN stabilization and enzyme activity. Further mutagenesis studies are needed to elucidate the role of other potential amino acid residues in the enzyme.
| ACKNOWLEDGEMENTS |
|---|
Edited by: H. L. Drake
| REFERENCES |
|---|
|
|
|---|
Chen, H., Wang, R. F. & Cerniglia, C. E. (2004). Molecular cloning, overexpression, purification, and characterization of an aerobic FMN-dependent azoreductase from Enterococcus faecalis. Protein Expr Purif 34, 302–310.[CrossRef][Medline]
Chen, H., Hopper, S. L. & Cerniglia, C. E. (2005). Biochemical and molecular characterization of an azoreductase from Staphylococcus aureus, a tetrameric NADPH-dependent flavoprotein. Microbiology 151, 1433–1441.
Chung, K. T. (1983). The significance of azo-reduction in the mutagenesis and carcinogenesis of azo dyes. Mutat Res 114, 269–281.[CrossRef][Medline]
Deller, S., Sollner, S., Trenker-El-Toukhy, R., Jelesarov, I., Gubitz, G. M. & Macheroux, P. (2006). Characterization of a thermostable NADPH : FMN oxidoreductase from the mesophilic bacterium Bacillus subtilis. Biochemistry 45, 7083–7091.[CrossRef][Medline]
Eckburg, P. B., Bik, E. M., Bernstein, C. N., Purdom, E., Dethlefsen, L., Sargent, M., Gill, S. R., Nelson, K. E. & Relman, D. A. (2005). Diversity of the human intestinal microbial flora. Science 308, 1635–1638.
Hanauer, S. B. (1996). Inflammatory bowel disease. N Engl J Med 334, 841–848.
Ito, K., Nakanishi, M., Lee, W. C., Sasaki, H., Zenno, S., Saigo, K., Kitade, Y. & Tanokura, M. (2006). Three-dimensional structure of AzoR from Escherichia coli. An oxidoreductase conserved in microorganisms. J Biol Chem 281, 20567–20576.
Laemmli, U. K. (1970). Cleavage of structural proteins during the assembly of the head of bacteriophage T4. Nature 227, 680–685.[CrossRef][Medline]
Li, Z. & Meighen, E. A. (1995). Tryptophan 250 on the alpha subunit plays an important role in flavin and aldehyde binding to bacterial luciferase. Effects of W
Y mutations on catalytic function. Biochemistry 34, 15084–15090.[CrossRef][Medline]
Liu, Z. J., Chen, H., Shaw, N., Hopper, S. L., Chen, L., Chen, S., Cerniglia, C. E. & Wang, B. C. (2007). Crystal structure of an aerobic FMN-dependent azoreductase (AzoA) from Enterococcus faecalis. Arch Biochem Biophys 463, 68–77.[CrossRef][Medline]
Lostao, A., Gomez-Moreno, C., Mayhew, S. G. & Sancho, J. (1997). Differential stabilization of the three FMN redox forms by tyrosine 94 and tryptophan 57 in flavodoxin from Anabaena and its influence on the redox potentials. Biochemistry 36, 14334–14344.[CrossRef][Medline]
Lostao, A., Daoudi, F., Irun, M. P., Ramon, A., Fernandez-Cabrera, C., Romero, A. & Sancho, J. (2003). How FMN binds to Anabaena apoflavodoxin: a hydrophobic encounter at an open binding site. J Biol Chem 278, 24053–24061.
Maiti, R., Van Domselaar, G. H., Zhang, H. & Wishart, D. S. (2004). SuperPose: a simple server for sophisticated structural superposition. Nucleic Acids Res 32, W590–W594.
McDonald, I. K. & Thornton, J. M. (1994). Satisfying hydrogen bonding potential in proteins. J Mol Biol 238, 777–793.[CrossRef][Medline]
Milburn, D., Laskowski, R. A. & Thornton, J. M. (1998). Sequences annotated by structure: a tool to facilitate the use of structural information in sequence analysis. Protein Eng 11, 855–859.
Murray, T. A. & Swenson, R. P. (2003). Mechanism of flavin mononucleotide cofactor binding to the Desulfovibrio vulgaris flavodoxin. 1. Kinetic evidence for cooperative effects associated with the binding of inorganic phosphate and the 5'-phosphate moiety of the cofactor. Biochemistry 42, 2307–2316.[CrossRef][Medline]
Platzek, T., Lang, C., Grohmann, G., Gi, U. S. & Baltes, W. (1999). Formation of a carcinogenic aromatic amine from an azo dye by human skin bacteria in vitro. Hum Exp Toxicol 18, 552–559.
Schwede, T., Kopp, J., Guex, N. & Peitsch, M. C. (2003). SWISS-MODEL: an automated protein homology-modeling server. Nucleic Acids Res 31, 3381–3385.
Sobolev, V., Sorokine, A., Prilusky, J., Abola, E. E. & Edelman, M. (1999). Automated analysis of interatomic contacts in proteins. Bioinformatics 15, 327–332.
Stolz, A. (2001). Basic and applied aspects in the microbial degradation of azo dyes. Appl Microbiol Biotechnol 56, 69–80.[CrossRef][Medline]
Velec, H. F., Gohlke, H. & Klebe, G. (2005). DrugScore (CSD)-knowledge-based scoring function derived from small molecule crystal data with superior recognition rate of near-native ligand poses and better affinity prediction. J Med Chem 48, 6296–6303.[CrossRef][Medline]
Wallace, A. C., Laskowski, R. A. & Thornton, J. M. (1995). LIGPLOT: a program to generate schematic diagrams of protein-ligand interactions. Protein Eng 8, 127–134.
Wang, C. J., Hagemeier, C., Rahman, N., Lowe, E., Noble, M., Coughtrie, M., Sim, E. & Westwood, I. (2007). Molecular cloning, characterisation and ligand-bound structure of an azoreductase from Pseudomonas aeruginosa. J Mol Biol 373, 1213–1228.[CrossRef][Medline]
Wu, K., Knox, R., Sun, X. Z., Joseph, P., Jaiswal, A. K., Zhang, D., Deng, P. S. & Chen, S. (1997). Catalytic properties of NAD(P)H : quinone oxidoreductase-2 (NQO2), a dihydronicotinamide riboside dependent oxidoreductase. Arch Biochem Biophys 347, 221–228.[CrossRef][Medline]
Received 22 April 2008;
revised 29 May 2008;
accepted 9 June 2008.
| ||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
| HOME | HELP | FEEDBACK | SUBSCRIPTIONS | ARCHIVE | SEARCH | TABLE OF CONTENTS |
| INT J SYST EVOL MICROBIOL | MICROBIOLOGY | J GEN VIROL |
| J MED MICROBIOL | ALL SGM JOURNALS | |