Previous Article | Next Article ![]()
Applied and Environmental Microbiology, August 2007, p. 4959-4965, Vol. 73, No. 15
0099-2240/07/$08.00+0 doi:10.1128/AEM.02993-06
Copyright © 2007, American Society for Microbiology. All Rights Reserved.

Chao Yang,1,2,
Hong Qu,3
Zheng Liu,1,2
Q. S. Fu,4 and
Chuanling Qiao1*
State Key Laboratory of Integrated Management of Pest Insects & Rodents, Institute of Zoology, Chinese Academy of Sciences, Beijing 100101, China,1 Graduate School of the Chinese Academy of Sciences, Beijing 100049, China,2 College of Life Sciences, Peking University, Beijing 100871, China,3 Department of Civil and Environmental Engineering, Stanford University, Stanford, California 94305-40204
Received 26 December 2006/ Accepted 31 May 2007
|
|
|---|
|
|
|---|
Various attempts have been made to isolate malathion-degrading microorganisms. Matsumura et al. isolated a malathion-degrading fungus, Trichoderma viride, and a bacterium, Pseudomonas sp., from contaminated soil, which rapidly degrades malathion through carboxylesteratic hydrolysis as well as desmethylation processes (22). Lewis et al. isolated a fungus, Aspergillus oryzae, from a freshwater pond, which transforms malathion to ß-malathion monoacid and malathion dicarboxylic acid, indicating the presence of carboxylesteratic hydrolysis activity (19). Degradation of malathion by other bacterial species was also reported (11). Although many studies have been carried out on the microbial degradation of malathion, there is a lack of information on the genetic and enzymatic aspects in the degradation of malathion.
The aldo-keto reductases (AKRs) are a growing superfamily of approximately 120 enzymes, currently composed of 15 families from a wide variety of plants, animals, and prokaryotes (17). They catalyze the reversible NAD(P)H-dependent reduction of a wide range of ketones and/or aldehydes containing substrates including aliphatic and polycyclic aldehydes, aldoses, and lipid-derived aldehydes to the corresponding alcohols, which are endogenously produced intermediates in many metabolic pathway (20) and involved in the detoxification of carbonyl group-containing xenobiotics (2); methylglyoxal (39); antiemetic, antitumor drugs, and the carcinogen 4-methylnitrosamino-1-(3-pyridyl)-1-butanone NNK (1, 21); polycyclic aromatic hydrocarbons (5); and naloxone, naltrexone, and dihydromorphinone (40). The AKR superfamily is also found in a wide range of microorganisms, and the well-reported superfamily members include the xylose reductases, 2,5-diketo-D-gluconic acid reductases (2,5-DKGR A), and ß-keto ester reductases, etc. These enzymes can metabolize many endogenous intermediate products, such as aldehydes and ketones, and exogenous compounds, including plant phytoallexins, toxins, and anthropogenic chemicals in the environment, as well as aldehydic products of lipid peroxidation (7). However, members of the AKR superfamily capable of metabolizing exogenous OP compounds have not been reported to date.
In the present study, a gene homologous to the AKR superfamily and responsible for the reduction of malathion as well as other thion OP compounds was cloned from a soil bacterium and expressed in Escherichia coli. Kinetic constants of the recombinant enzyme were also determined against various OP compounds.
|
|
|---|
|
View this table: [in a new window] |
TABLE 1. Bacterial strains and plasmids used in this study
|
Cloning an akr5f1 gene and data analysis.
Routine DNA manipulations were carried out as described by Sambrook et al. (26). A size-fractionated genomic library was constructed from 4- to 10-kb DNA fragments (obtained from a partial Sau3AI digest) cloned into the BamHI site of plasmid pBluescript KS II (pBS) and transformed into E. coli DH5
(16, 43). Screenings for positive clones were performed spectrophotometrically by monitoring the absorbance at 410 nm resulting from the free thiol group produced from the P
S reduction of malathion (16, 18), and the inserted fragment from a positive clone was sequenced (Invitrogen). BlastN was used for the nucleotide sequence identity search, and BLASTP was used for the deduced amino acid identity search (www.ncbi.nlm.nih.gov/BLAST). Open reading frames (ORFs) were identified using the NCBI ORF finder tool. Promoter prediction was done online (www.fruitfly.org/seq_tools/promoter.html). The 5' region of the deduced AKR gene was checked for the presence of ribosome binding site as described by Stormo et al. (32). Multiple alignments of amino acid sequences, construction of a neighbor-joining phylogenetic tree with a homology rate model, and a bootstrap analysis for evaluation of the phylogenetic topology were accomplished by using the Clustal X program (33) and MEGA 3 software (www.megasoftware.net).
Plasmid construction for overexpression in E. coli.
Based on the two deduced ORFs, ra and r, in the sequenced fragment, oligonucleotide primers were synthesized to amplify the intact coding regions. The restriction sites for EcoRI and HindIII (underlined) were incorporated into the forward and the reverse primer sequences, respectively. The following two pairs of primers were used: pair I (forward, 5' GGCGAATTCATGCAAACTGTAAAAC 3'; reverse, 5' CCGAAGCTTTTGCGCCTTAAACATC 3') for amplification of ORF ra and pair II (forward, 5' GAATTCATGCAAAAACGTTATCTG 3'; reverse, 5' AAGCTTTTAACGGCCTACGCGAGC3') for amplification of ORF r. Two PCR fragments were subcloned into EcoRI-HindIII sites of pET28a(+) (Novagen) to generate the recombinant plasmids pRA and pR, and the plasmids were subsequently transformed into E. coli BL21(DE3). The gene homologous to ORF ra was further proved to possess OP-metabolizing activity by enzyme assay (16, 18).
Overexpression and purification of fused AKR5F1.
AKR5F1 was overexpressed in E. coli BL21(DE3) by using the His-Bind protein fusion and purification system, which results in expression of a His-binding fusion protein. Optimal production of the fusion proteins was obtained when mid-log-phase cells (optical density at 600 nm of 0.5) were induced with 0.1 mM isopropyl-ß-D-thiogalactopyranoside (IPTG) for 16 h at 30°C. Harvested cells were washed and disrupted by sonication (12), and the soluble fraction was loaded onto a His-Bind resin (Novagen). Fusion proteins were eluted with 300 mM imidazole in 50 mM Tris-HCl (pH 6.5).
Enzyme assays.
Cell extracts were prepared as previously described (12). Protein concentrations were determined according to the Bradford method (Bio-Rad) using bovine serum albumin as the standard. AKR5F1 activities were measured spectrophotometrically by monitoring the absorbance at 410 nm due to the formation of a free thiol group by using Ellman's reagent (16, 18). The reactions were carried out in 50 mM Tris-HCl (pH 6.5) containing 0.3 mM malathion (or other OPs), 0.3 mM NAD(P)H, 60 µg of purified AKR5F1, and 1 mM dithionitrobenzoic acid at 28°C with ß-mercaptoethanol as a standard. The reaction systems without NAD(P)H were used as the controls.
Kinetic measurement.
All kinetic measurements were carried out in 50 mM Tris-HCl buffer (pH 6.5). Besides malathion, five other thion OPs, including chlorpyrifos, parathion, methyl parathion, dimethoate, and methidathion, were tested. The Michaelis-Menten kinetic constants for each substrate were obtained by measuring the initial hydrolysis rate at different substrate concentrations (0.08 to 0.8 mM) with a constant enzyme concentration of 60 µg/ml. The change in absorbance was measured with a VERSAmax microplate reader for 10 min at 30°C. All assays were performed in triplicate. The kinetic constants (Km and kcat) were calculated by the Lineweaver-Burk equation by using the Microcal origin software program.
Bioinformatic analysis.
Sequence alignments were obtained with Clustal X (33) program. The putative structure for AKR5F1 was predicted as a function of the sequence homology between this protein and the data obtained from the structures of AKR2B5 previously crystallized. (AKR2B5 is a complex of xylose reductase from Candida tenuis [Protein Data Bank accession no. 1MI3].) The predicted structure has been generated using SwissModel (http://swissmodel.expasy.org) (10, 25, 28) and was refined with the Discover program in Insight II (2000) software package (Accelrys, Inc.). The interactions of the AKR5F1-OP-NADH complex was modeled by Docking program.
Chemicals.
Malathion [O,O-dimethyl S-(1,2-dicarbethoxyethyl) phosphorodithioate] with 95% purity, chlorpyrifos (O,O-diethyl O-3,5,6-trichloro-2-pyridyl phosphorothioate) with 95% purity, parathion (O,O-diethyl p-nitrophenyl phosphorothioate) with 92% purity, and methyl parathion (O,O-dimethyl O-p-nitrophenyl phosphorothioate) with 86% purity were obtained from Institute for the Control of Agrochemicals, Ministry of Agriculture, China. Dimethoate (O,O-dimethyl S-methyl-carbamoylmethyl phosphorodithioate) with 94% purity and methidathion [O,O-dimethyl-S-(2-methoxy-1,3,4-thiadiazol-5 (4H)-onyl-(4)-methyl)-dithiophosphate] with 95% purity were a gift from China Agricultural University.
Nucleotide sequence accession numbers.
The nucleotide sequences of the 16S rRNA and gene akr5f1 of strain F51-1-2 have been deposited in the GenBank database under accession no. DQ277701 and DQ268871.
|
|
|---|
The substrate utilization of strain F51-1-2 was compared with that of referred strains in the BIOLOG-GN database, and strain F51-1-2 had the greatest similarity index of 0.68 with a Klebsiella sp.
Cloning and expression of the akr5f1 gene.
As an enzyme assay can be used to determine the metabolism of OPs using cell extracts of strain F51-1-2, it is practical to screen the fractionated genomic library. A positive clone showing enzymatic activity was screened from the genomic library and sequenced. The sequence of the inserted fragment was 4,918 bp and contained two complete ORFs, designated as r and ra, respectively. After induction with IPTG, E. coli BL21(DE3) harboring pRA with ORF ra showed the reductase activity against the tested OPs but not the same cells harboring pR with ORF r, which indicated that ORF ra encoded the true reductase for OP metabolism. Sequence BLAST analysis indicated ORF ra was approximately 87% identical at the nucleotide level and 95% identical at the amino acid level to a gene encoding 2,5-DKGR A (AE014075) from E. coli CFT073 (36). Based on sequence alignments, ORF ra was putatively identified as a member of the AKR superfamily and named akr5f1 by the AKR superfamily website (http://www.med.upenn.edu/akr/) (17). A phylogenetic tree of strain F51-1-2 based on multiple alignments of amino acid sequences is shown in Fig. 1.
![]() View larger version (36K): [in a new window] |
FIG. 1. Phylogenetic relationship of AKRs from bacteria. Bootstrap values obtained with 1,000 repetitions are indicated as percentages at all branches. GenBank accession numbers are given in brackets.
|
The high expression of His-tagged AKR5F1 fusion protein was achieved with 1 mM IPTG induction, and the fusion protein was purified with the His-Bind column. A clear band corresponding to a molecular mass of 32 kDa was observed by sodium dodecyl sulfate-polyacyrlamide gel electrophoresis (Fig. 2). Approximately 19 mg of purified protein was obtained from 500 ml of culture.
![]() View larger version (105K): [in a new window] |
FIG. 2. Sodium dodecyl sulfate-polyacrylamide gel electrophoresis analysis of purified AKR5F1. Lane a, protein markers (kDa); lane b, purified AKR5F1; lane c, total proteins from E. coli BL21 cells harboring pRA with IPTG induction; lane d, total proteins from the same cells without plasmid. Purified AKR5F1 is indicated by an arrow.
|
|
View this table: [in a new window] |
TABLE 2. Kinetic constants of His-tagged AKR5F1 toward various OPs with NADH as a cofactora
|
The kinetic constants are not significantly different among the tested OPs. Km is relatively higher and kcat is relatively lower than those of the corresponding substrates hydrolyzed by OPH. OPH shared kinetic constants (the Km and kcat values of OPH are 50.6 ± 12.2 µM and 23.5 ± 0.2 min–1, respectively, for parathion and 32.9 ± 1.7 µM and 5.46 ± 0.05 min–1, respectively, for methyl parathion) (16). The OP substrates with a P—S bond such as malathion, dimethoate, and methidathion, which cannot efficiently be hydrolyzed by OPH/OPDA, can be reduced by AKR5F1. In addition, OP substrates with either a P—O bond such as parathion and methyl parathion or a P—S bond such as malathion can be reduced by AKR5F1 in moderate efficiency, suggesting that the reduction of OPs by AKR5F1 may share similar mechanism.
Structures of OPH have been well documented (34), and three pockets known as the small, large, and leaving group subsites to accommodate the phosphate substituents have been identified (4). The substrate binding pockets observed in the OPDA structure are similar to those found in OPH, except that the dimension of the large subsite in the OPDA is smaller than that of OPH (42). Due to the reduction of this subsite, OPDA processes methyl substrate more efficiently than the diethyl equivalent, while the two compounds are processed in a similar manner by OPH (3, 14, 16). Unlike OPH/OPDA, AKR5F1 has only one large pocket, which holds the O,O-dimethyl/diethyl phosphorodithioate moiety of OPs as well as NADH. The dimension of the pocket binds dimethyl-substituted substrates slightly better than diethyl-substituted substrates (Table 3). The inferences derived from structural modeling were supported by the kinetic properties of purified AKR5F1. The enzyme exhibited higher kcat values for dimethyl-substituted substrates such as malathion, methyl parathion, dimethoate, and methidathion than diethyl-substituted substrates such as parathion and chlorpyrifos (Table 2).
|
View this table: [in a new window] |
TABLE 3. Distance parameters for the interactions of AKR5F1 and OPs with NADH cofactor
|
/ß-barrel consisting of eight
-helices and eight ß-strands with one active pocket inside. This type of structural fold has been demonstrated in many other enzymes (8) belonging to the AKR family (15, 37, 38). In the malathion-bound structure, the O,O-dimethyl phosphorodithioate moiety of malathion and NADH interact with the active site of the enzyme, while the 1,2-dicarbethoxyethyl moiety is exposed to the solvent (Fig. 3). The thiophosphoryl sulfur makes hydrogen bond interactions with Tyr49 and His107. The residues involved in forming the catalytic pocket in AKR5F1 are almost conserved and correspond to those in AKR5C1, with the exception of Glu188, which presents a carboxylic group instead of an amidic group of Gln192 in AKR5C1, suggesting that AKR5F1 might exhibit a manner of stereo-specific reduction similar to that of AKR5C1, which proceeds via a two-step reaction mechanism involving transferring a hydride ion (H–) from NADH to the substrate thiophosphoryl phosphorus, followed by transferring a proton (H+) from Tyr49 or His107 to the substrate thiophosphoryl sulfur. The data of distance parameters (Table 3) suggest that other OP compounds might make similar contact with AKR5F1 and NADH as described above.
![]() View larger version (37K): [in a new window] |
FIG. 3. Modeled three-dimensional structure of AKR5F1 from Klebsiella sp. strain F51-1-2 with the malathion-NADH complex. The backbone of AKR5F1 is represented in blue. The residues involved in the substrate binding pocket are highlighted by blue sticks. NADH and malathion are denoted by colored sticks. The hydrogen bonds formed between AKR5F1 and substrate are indicated by green broken lines.
|
![]() View larger version (56K): [in a new window] |
FIG. 4. Sequence alignment of amino acid sequences of AKR5F1 from Klebsiella sp. strain F51-1-2 (NCBI accession no. DQ268871), AKR2B5 from Corynebacterium (NCBI accession no. AF074484), and AKR5C1 from Candida tenuis (NCBI accession no. M12799). The identity or similarity of residues is represented by symbols beneath the sequences, with asterisks representing the highest level of similarity, followed by colons and then periods. The residues involved in formation of the binding pocket are highlighted in gray.
|
In most of the studies on microbial degradation of organophosphorus pesticides, the first reaction was hydrolysis of the phosphotriester bond (23). Conversion of thion (P
S) to oxon (P
O) OPs has also been reported (24). Recently, we found a Bacillus sp. strain that transforms parathion and methyl parathion by reducing the nitro group and a nitroreductase involved in the transformation (41). In the present study, we report a transformation pathway that has not been previously described: the reduction of the P
S bond of the thion OPs by Klebsiella sp. Moreover, we present the identification of the gene/enzyme system involved in the transformation. These results suggest that several types of reductases may play important roles in microbial degradation of organophosphates. Here, for the first time we report a novel AKR AKR5F1 that can reduce a kind of the thion OPs. The enzyme could be potentially utilized for reducing a double bond linking phosphate and sulfur in the OP molecules, which will drive the fields of both organophosphate degradation and AKR enzymology forward.
Our findings demonstrated that the novel AKR5F1 functions in OP metabolism in the presence of NADH. Moreover, the novel AKR gene involved in OP metabolism provided valuable information for exploring the course of evolution of OP metabolic genes in bacteria. Elucidation of the reductive pathway and further characterization of the enzymes involved in degradation of OPs in strain F51-1-2 represent areas for further investigation.
We thank Wensheng Lan and Yulan Yuan for their technical assistance.
Published ahead of print on 15 June 2007. ![]()
H.J. and C.Y. contributed equally to this work. ![]()
|
|
|---|
/ß pairs. Photosynth. Res. 49:223-235.[CrossRef]
| |||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
Copyright © 2009 by the American Society for Microbiology. For an alternate route to Journals.ASM.org, visit: http://intl-journals.asm.org | More Info»