Biochemical characterization of an enantioselective esterase from Brevundimonas sp. LY-2
© The Author(s) 2017
Received: 18 January 2017
Accepted: 14 June 2017
Published: 19 June 2017
Lactofen, a member of the diphenylether herbicides, has high activity and is commonly used to control broadleaf weeds. As a post-emergent herbicide, it is directly released to the environment, and easily caused the pollution. This herbicide is degraded in soil mainly by microbial activity, but the functional enzyme involved in the biodegradation of lactofen is still not clear now.
A novel esterase gene lacH, involved in the degradation of lactofen, was cloned from the strain Brevundimonas sp. LY-2. The gene contained an open reading frame of 921 bp, and a putative signal peptide at the N-terminal was identified with the most likely cleavage site between Ala 28 and Ala 29. The encoded protein, LacH, could catalyze the hydrolysis of lactofen to form acifluorfen. Phylogenetic analysis showed that LacH belong to family V of bacterial lipolytic enzymes. Biochemical characterization analysis showed that LacH was a neutral esterase with an optimal pH of 7.0 and an optimal temperature of 40 °C toward lactofen. Besides, the activity of LacH was strongly inhibited by Hg2+ and Zn2+. LacH preferred short chain p-nitrophenyl esters (C2–C6), exhibited maximum activity toward p-nitrophenyl acetate. Furthermore, the enantioselectivity of LacH during lactofen hydrolysis was also studied, and the results show that R-(−)-lactofen was degraded faster than S-(+)-lactofen, indicating the occurrence of enantioselectivity in the enzymatic reaction.
Our studies characterized a novel esterase involved in the biodegradation of diphenylether herbicide lactofen. The esterase showed enantioselectivity during lactofen degradation, which revealed the occurrence of enzyme-mediated enantioselective degradation of chiral herbicides.
KeywordsEsterase Lactofen Brevundimonas sp. Enantioselective degradation
Many pesticides in current use contain chiral structures, which consist of enantiomers . Enantiomers of a chiral compound exhibit identical physical and chemical properties . However, their toxicities, biological activities and environmental fates vary because biological processes usually show high enantioselectivity [3, 4]. Enantioselective analysis is required for a full understanding of the biological behavior of chiral compounds. Information on the enantioselective degradation dynamics of chiral pesticides will help us to evaluate the impacts of such pesticides to environment. Enantioselective degradation in soils has been observed for various chiral pesticides, such as mecoprop, dichlorprop, metalaxyl and malathion [5–7]. When one enantiomer is preferentially degraded, the enantiomeric ratio [ER, the ratio of the concentration of (+)-enantiomers to (−)-enantiomers] deviates from the original value [2, 8]. The occurrence of such selective degradation involves the mediation of bacteria, enzymes, or other biological entities.
Diphenyl ether herbicides are widely used to control broadleaf weeds in cereal crops, soybeans, peanuts, and potatoes . This class of herbicides has been proved to inhibit protoporphyrinogen oxidase, thereby resulting in the accumulation of protoporphyrin and blockage of chlorophyll formation [10, 11]. Their frequent occurrence in natural waters and soils indicate that they may be important environmental contaminants [12–14]. Lactofen is an important member of the diphenyl ether family. It has one chiral center in the alkyl moiety and consists of two enantiomers. This herbicide is applied as a racemic mixture, although its herbicidal activity is almost entirely attributed to the S-enantiomer . The enantioselective degradation of lactofen in soil and sediment has been studied by Diao et al. [15, 16]. The S-enantiomer was preferentially degraded either in soil or sediment, resulting in a relative enrichment of the R-enantiomer. However, to our knowledge, no studies have focused on the enantioselectivity of the degrading-enzyme involved in lactofen metabolism so far.
The lactofen-degrading strain Brevundimonas sp. LY-2 was isolated from enrichment cultures inoculated with lactofen-contaminated soil sample in our lab. This strain could degrade about 80% of 50 mg L−1 lactofen in 5 days of incubation in flasks. The metabolic behaviors of the herbicide in the media were described . In this study, we have found that the degradation process of lactofen by LY-2 is enantioselective, with R-(−)-lactofen being degraded faster than the S-(+)-enantiomer, which implies that the enzyme involved in lactofen degradation probably has the enantioselectivity between the different enantiomers. In this study, a novel esterase gene lacH involved in lactofen degradation was cloned from this strain, and the biochemical properties of the purified enzyme LacH were determined. The enantioselectivity of LacH during lactofen hydrolysis was also investigated.
Cloning and sequence analysis of the esterase gene
The total DNA was extracted from strain Brevundimonas sp. LY-2, and the genomic library was constructed. A positive clone that produced a transparent halo around the colony was screened from approximately 3000 transformants. The sequencing results show that the inserted fragment in the transformant was 3453 bp containing three complete ORFs. These ORFs were each subcloned into the pMD18-T vector and then transformed into competent Escherichia coli DH5α cells. One ORF, designated as lacH, was confirmed to be the target gene encoding the lactofen-hydrolyzing enzyme. This ORF consists of 921 bp that encode a protein with 306 amino acids. A putative signal peptide at the N-terminal was identified by the SignalP 4.1 server, with the most likely cleavage site situated between amino acids Ala 28 and Ala 29. Thus, the encoded protein, LacH, is a secretory protein. The BLAST results showed that LacH shared moderate identities with some characterized α/β-hydrolase fold proteins, and showed the highest identity (71%) with esterase sys410 (AFE88176), a thermostable pyrethroid-hydrolyzing enzyme isolated from uncultured bacterium through the metagenomic approach . Among the other characterized proteins, LacH shared 48% identity with an α/β-hydrolase fold protein from Caulobacter crescentus CB15 (AAK24201) , 40% identity with a triacylglycerol lipase from Moraxella sp. TA144 (CAA37863) , and only 27% identity with a triacylglycerol lipase from Psychrobacter immobilis B10 (CAA47949) .
Expression and purification of the recombinant LacH
Characteristics of LacH
Substrate specificity and activity of LacH
Substrate specificity and kinetic constants of LacH
Specific activity (μmol min−1 mg−1)
k cat (s−1)
K m (mM)
k cat/K m (s−1 mM−1)
2.97 ± 0.61
0.87 ± 0.09
0.067 ± 0.002
118 ± 5.3
8.39 ± 0.28
0.147 ± 0.008
88.1 ± 3.2
7.48 ± 0.45
0.185 ± 0.009
15.8 ± 1.2
3.47 ± 0.13
0.208 ± 0.012
4.90 ± 0.18
0.767 ± 0.06
0.335 ± 0.016
Degradation efficiency by LacH with lactofen as the substrate was tested by high-performance liquid chromatography (HPLC) analysis (Additional file 1: Figure S1). The hydrolysis rate of lactofen was 59.3% under assay conditions of pH 7.0 and 40 °C for 15 min (Additional file 1: Table S1). The metabolite of lactofen hydrolysis by LacH was confirmed to be acifluorfen (Additional file 1: Figure S1).
Enantioselectivity of LacH
Esterases (EC 184.108.40.206) and lipases (EC 220.127.116.11) are widely distributed in animals, plants and microorganisms. They have played an important role in the hydrolysis of a wide range of xenobiotics containing ester bonds. Esterases and lipases found in bacteria have been divided into eight different families based on their amino acid sequences and biochemical properties . The esterase LacH that was cloned from strain LY-2 was responsible for converting lactofen to acifluorfen. Cheng et al.  studied the acute toxicity of lactofen and its metabolite acifluorfen in the aquatic algae Scenedesmus obliquus, and the results indicated that the toxicity of the acifluorfen was lower than that of the lactofen. Wang et al.  evaluated the toxic effects of the herbicide lactofen and acifluorfen to the aquatic plant Lemna minor, and they found L. minor growth was inhibited in the order of lactofen > acifluorfen. Both of these studies have found that acifluorfen had a lower toxicity to aquatic organisms than lactofen. To the best of our knowledge, LacH is the first identified esterase for lactofen hydrolysis in microorganisms. LacH belonged to the α/β-hydrolase fold protein family and contained both the typical Ser-His-Asp/Glu (Ser128, Asp233 and His286) catalytic triad and the conserved pentapeptide sequence GXSXG (residues 126–130) of α/β-hydrolase fold proteins [22, 23]. Amino acid sequence alignment and phylogenetic analysis suggested that LacH belonged to family V of lipolytic enzymes.
The substrate spectrum of LacH was found to be broad. It was able to hydrolyze various p-nitrophenyl esters of short-medium chain fatty acids (C2–C8). The catalytic efficiency value (k cat/K m) demonstrated that p-nitrophenyl acetate (C2) was the most efficient catalytic substrate. The specificity profile indicated that LacH should be classified as an esterase because maximal activity was observed on the short-chain fatty acid esters, while the activity of long-chain fatty acid esters was low. LacH could also hydrolyze other pesticides with similar carboxyl ester, such as fenoxaprop-P-ethyl, cyhalofop-butyl, fluazifop-P-butyl, quizalofop-P-ethyl and fluoroglycofen (Additional file 1: Table S1).
The degradation process of lactofen by Brevundimonas sp. LY-2 is enantioselective, with R-(−)-lactofen being degraded faster than the S-(+)-enantiomer (Additional file 1: Table S2), which implies that the enzyme involved in lactofen conversion can differentiate between the enantiomers. Previous studies have shown that the herbicidally active S-(+)-lactofen was preferentially degraded either in soil or sediment, resulting in relative enrichment of the R-(−)-form [15, 16]. However, studies on enantioselectivity of biotransformation enzyme involved in lactofen metabolism by microorganisms are limited. In this study, we demonstrated the occurrence of enantioselectivity in biotransformation of lactofen by LacH and the rising of ER indicated the R-(−)-lactofen was preferentially degraded over the S-(+)-form.
Up to now, little is known about the functional enzyme involved in the degradation of lactofen by pure culture. In this paper, a novel esterase LacH was functionally expressed and characterized from Brevundimonas sp. LY-2. This is the first report of a functional enzyme involved in the degradation of diphenylether herbicide. In this paper, the enantioselectivity of LacH during lactofen hydrolysis was studied, and the results show that R-(−)-lactofen was degraded faster than S-(+)-lactofen, indicating the occurrence of enantioselectivity in the enzymatic reaction, which may help in understanding the mechanism of enantioselective degradation of chiral pesticides. In the future, for a better understanding of enantioselective degradation, detailed studies on recognition mechanism of enzymes to differentiate between enantiomers will be needed.
Lactofen (99% purity) was purchased from Sigma-Aldrich Chemical Co. (Shanghai, China). All the p-nitrophenyl esters were purchased from Sigma. Methanol, n-hexane, and isopropanol were of pure chromatographic grade. All other chemicals used were of analytical grade.
Bacterial strains and plasmids
Brevundimonas sp. LY-2 was deposited in our laboratory. E. coli DH5α and E. coli BL21 (DE3) were used as the host for gene cloning and protein expression, respectively. The pUC118 and pET-29a (+) were used to construct the genomic library and express the target protein, respectively.
Screening of genomic library
DNA manipulation was performed as described by Sambrook and Russell . To construct the genomic library, genomic DNA of Brevundimonas sp. LY-2 was subjected to partial digestion with Sau3AI. Fractions containing approximately 2–4 kb DNA fragments were pooled, ligated into the BamHI site of the plasmid pUC118, and transformed into E. coli DH5α. The library was plated onto LB agar containing 100 mg L−1 ampicillin and 100 mg L−1 lactofen, and incubated at 37 °C for approximately 24 h. The transformants harboring the plasmid containing the functional gene and producing transparent halos were screened and further tested by HPLC analysis to determine their lactofen-degrading ability. Positive recombinant plasmids were extracted and submitted for sequencing.
Nucleotide and deduced amino acid sequence analyses were performed using Omiga Software 2.0. BlastN and BlastP were used for the nucleotide sequence and deduced amino acid identity searches (http://www.ncbi.nlm.nih.gov/Blast), respectively. The signal peptide of the amino acid sequence was predicted using SignalP 4.1 (http://www.cbs.dtu.dk/services/SignalP/). Multiple sequence alignment was performed with Clustal W program . Phylogenetic analysis was carried out via the neighbor-joining method using MEGA 6.0. Bootstrapping of 1000 replicates was performed to estimate the confidence levels of phylogenetic reconstructions .
Gene expression and purification of the recombinant LacH
The open reading frame (ORF) of lacH without the signal peptide was amplified by PCR using primers F1 (5′-GGAATTCCATATGGCGACCGCCGAACCGG-3′; NdeI site was underlined) and R1 (5′-CCGCTCGAGGTCCGCCAGGAACCGGTC-3′; XhoI site was underlined). The PCR product, digested by NdeI and XhoI, was inserted into pET-29a (+) and the recombinant plasmid was then transformed into E. coli BL21 (DE3). The fusion protein was obtained when cells in mid-log phase (optical density at 600 nm = 0.6) were induced with 1.0 mM isopropyl-β-d-thiogalactopyranoside at 30 °C for 3 h. The harvested cells were washed and disrupted by sonication. Cell debris was removed by centrifugation. The supernatant was loaded onto a His-Bind resin (Novagen). The target protein was eluted with 100 mM imidazole after elution of non-target proteins with 25 and 50 mM imidazole. The enzyme was dialysis against 50 mM sodium phosphate buffer (pH 7.0) for 24 h, and concentrated using an Amicon ultrafiltration tube. The protein concentration was quantified by the Bradford method using bovine serum albumin (BSA) as the standard. The purified protein was analyzed by 12% SDS-PAGE gel electrophoresis. The molecular mass markers used were rabbit muscle phosphorylase b (97.2 kDa), BSA (66.4 kDa), hen egg white ovalbumin (44.3 kDa), bovine carbonic anhydrase (29.0 kDa), soybean trypsin inhibitor (20.1 kDa), and hen egg white lysozyme (14.3 kDa). The molecular mass of the native protein was determined by gel filtration on a Superdex 200 column. The molecular mass markers used were thyroglobulin (669 kDa), ferritin (440 kDa), bovine gamma globulin (158 kDa), chicken ovalbumin (44 kDa), carbonic anhydrase (29 kDa), and ribonuclease A (13.7 kDa).
The optimal temperature for lactofen degradation by pure enzyme was determined in 50 mM sodium phosphate buffer (pH 7.0) at different temperatures from 15 to 70 °C. To determine thermostability, the enzyme was pre-incubated at various temperatures ranging from 15 to 80 °C for 30 min, and the residual activity was determined. The optimal pH was measured at 40 °C with pH ranging from 4.0 to 10.0. The buffers used were 50 mM citric acid–sodium citrate (pH 4.0–pH 6.0), sodium phosphate (pH 6.0–pH 8.0), Tris–HCl (pH 8.0–pH 9.0), and glycine–NaOH (pH 9.0–pH 10.0). To investigate the effects of metal ions on enzyme activity, the enzyme was pre-incubated with various metal salts at 40 °C for 30 min. Metal salts (FeSO4, CdCl2, CaCl2, MnCl2, MgSO4, CuSO4, ZnSO4, HgCl2 and NiSO4) were added at a final concentration of 1 mM. The residual activity was assayed and expressed as a percentage of the activity obtained in the absence of added metal ions.
Substrate specificity and kinetics study
Substrate specificity of LacH was determined using lactofen, p-nitrophenyl acetate (C2), p-nitrophenyl butyrate (C4), p-nitrophenyl caproate (C6), p-nitrophenyl caprylate (C8), and p-nitrophenyl laurate (C12) as substrates. Hydrolytic activity toward p-nitrophenyl esters was assayed in accordance with the method described by Gao et al. . Hydrolytic activity toward lactofen was performed as described below. For kinetic studies, stock solutions of each substrate were appropriately diluted into five different concentrations around the K m values. Kinetic values were obtained from Lineweaver–Burk plots against various substrate concentrations.
HPLC–MS/MS analysis of lactofen degradation by LacH
Hydrolysis of lactofen was assayed in 50 mM sodium phosphate buffer (pH 7.0) at 40 °C with a final concentration of 50 mg L−1 lactofen. Briefly, 1.5 μL of lactofen stock solution (100,000 mg L−1, dissolved in methanol) was added to 3 mL sodium phosphate buffer (50 mM, pH 7.0). The reaction was initiated by the addition of 50 μL of the purified LacH (57.5 μg mL−1) and incubation for 15 min. The substrate residues were measured by HPLC. Control samples containing boiled enzyme solution were treated and analyzed in the same way. The metabolite of lactofen degradation by LacH was further confirmed by MS/MS .
For lactofen determination, the solution mixture was extracted with an equal volume of dichloromethane. The organic layer was dried over anhydrous Na2SO4, and dichloromethane was removed using a stream of nitrogen at room temperature. The residues were re-dissolved in methanol. All samples were analyzed by HPLC equipped with a Zorbax C-18 ODS Spherex column (250 mm × 4.6 mm). The mobile phase was pure methanol and the flow rate was 1 mL min−1. The UV wavelength for detection was 230 nm, and the injection volume was 20 μL.
The lactofen enantiomers were separated and quantified using HPLC with a Chiralcel OD-H column (250 mm × 4.6 mm). The mobile phase was n-hexane–isopropanol (97:3, vol/vol), and the flow rate was 1 mL min−1. The UV wavelength for detection was 230 nm, and the injection volume was 20 μL.
Nucleotide sequence accession number
The nucleotide sequence of the lacH gene was deposited in the GenBank database under accession number KF286657.
Brevundimonas sp. LY-2 was deposited to the China General Microbiological Culture Collection Center (CGMCC3651).
Statistical analysis and reproducibility
All experiments were performed in triplicate. The values were expressed as mean ± SD in the figures.
XH and JZ designed the experiments. JZ, MZ, DY, JY and HZ performed the experiments and analyzed the data. XH and JZ wrote the paper. All authors read and approved the final manuscript.
The authors would like to thank Pro Jianhe Xu at East China University of Science and Technology for his valuable advices about this paper.
The authors declare that they have no competing interests.
Availability of data and materials
The datasets supporting the conclusions of this article are included within the article and its additional files.
This work was funded by The National Natural Science Foundation of China (41671317), the Fundamental Research Funds for the Central Universities (KYZ201635).
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Open AccessThis article is distributed under the terms of the Creative Commons Attribution 4.0 International License (http://creativecommons.org/licenses/by/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated.
- Sekhon BS. Chiral pesticides. J Pestic Sci. 2009;34:1–12.View ArticleGoogle Scholar
- Liu WP, Gan JY, Schlenk D, Jury WA. Enantioselectivity in environmental safety of current chiral insecticides. Proc Natl Acad Sci. 2005;102:701–6.View ArticleGoogle Scholar
- Garrison AW. Probing the enantioselectivity of chiral pesticides. Environ Sci Technol. 2006;40:16–23.View ArticleGoogle Scholar
- Müller TA, Kohler HP. Chirality of pollutants-effects on metabolism and fate. Appl Microbiol Biotechnol. 2004;64:300–16.View ArticleGoogle Scholar
- Romero E, Matallo MB, Pena A, Sánchez-Rasero F, Schmitt-Kopplin P, Dios G. Dissipation of racemic mecoprop and dichlorprop and their pure R -enantiomers in three calcareous soils with and without peat addition. Environ Pollut. 2001;111:209–15.View ArticleGoogle Scholar
- Buerge IJ, Poiger T, Müller MD, Buser H-R. Enantioselective degradation of metalaxyl in soils: chiral preference changes with soil pH. Environ Sci Technol. 2003;37:2668–74.View ArticleGoogle Scholar
- Sun MJ, Liu DH, Zhou GX, Li JD, Qiu XX, Zhou ZQ, Wang P. Enantioselective degradation and chiral stability of malathion in environmental samples. J Agric Food Chem. 2011;60:372–9.View ArticleGoogle Scholar
- Hegeman WJM, Laane RWPM. Enantiomeric enrichment of chiral pesticides in the environment. Rev Environ Contam Toxicol. 2002;173:85–116.Google Scholar
- Ahrens WH. Herbicide handbook. Champaign: Weed Science Society of America; 1994.Google Scholar
- Geoffroy L, Teisseire H, Couderchet M, Vernet G. Effect of oxyfluorfen and diuron alone and in mixture on antioxidative enzymes of Scenedesmus obliquus. Pestic Biochem Physiol. 2002;72:178–85.View ArticleGoogle Scholar
- Jinno H, Hatakeyama N, Hanioka N, Yoda R, Nishimura T, Ando M. Cytotoxic and porphyrinogenic effects of diphenyl ethers in cultured rat hepatocytes: chlornitrofen (CNP), CNP-amino, chlomethoxyfen and bifenox. Food Chem Toxicol. 1999;37:69–74.View ArticleGoogle Scholar
- Laganà A, Fago G, Fasciani L, Marino A, Mosso M. Determination of diphenyl-ether herbicides and metabolites in natural waters using high-performance liquid chromatography with diode array tandem mass spectrometric detection. Anal Chim Acta. 2000;414:79–94.View ArticleGoogle Scholar
- Sheu HL, Sung YH, Melwanki MB, Huang SD. Determination of diphenylether herbicides in water samples by solid-phase microextraction coupled to liquid chromatography. J Sep Sci. 2006;29:2647–52.View ArticleGoogle Scholar
- Sondhia S. Persistence of oxyfluorfen in soil and detection of its residues in rice crop. Toxicol Environ Chem. 2009;91:425–33.View ArticleGoogle Scholar
- Diao JL, Xu P, Wang P, Lu DH, Lu YL, Zhou ZQ. Enantioselective degradation in sediment and aquatic toxicity to daphnia magna of the herbicide lactofen enantiomers. J Agric Food Chem. 2010;58:2439–45.View ArticleGoogle Scholar
- Diao JL, Lv CG, Wang XQ, Dang ZH, Zhu WT, Zhou ZQ. Influence of soil properties on the enantioselective degradation of the herbicide lactofen in soils. J Agric Food Chem. 2009;57:5865–71.View ArticleGoogle Scholar
- Liang B, Zhao YK, Lu P, Li SP, Huang X. Biotransformation of the diphenyl ether herbicide lactofen and purification of a lactofen esterase from Brevundimonas sp. LY-2. J Agric Food Chem. 2010;58:9711–5.View ArticleGoogle Scholar
- Fan XJ, Liu XL, Huang R, Liu YH. Identification and characterization of a novel thermostable pyrethroid-hydrolyzing enzyme isolated through metagenomic approach. Microb Cell Fact. 2002;11:33.View ArticleGoogle Scholar
- Iniesta AA, McGrath PT, Reisenauer A, McAdams HH, Shapiro LA. Phospho-signaling pathway controls the localization and activity of a protease complex critical for bacterial cell cycle progression. Proc Natl Acad Sci. 2006;103:10935–40.View ArticleGoogle Scholar
- Feller G, Thiry M, Gerday C. Nucleotide sequence of the lipase Gene lip3 from the antarctic psychotroph Moraxella TA144. Biochim Biophys Acta. 1991;1088:323–4.View ArticleGoogle Scholar
- Arpigny JL, Feller G, Gerday C. Cloning, sequence and structural features of a lipase from the antarctic facultative psychrophile Psychrobacter immobilis B10. Biochim Biophys Acta. 1993;1171:331–3.View ArticleGoogle Scholar
- Ollis DL, Cheah E, Cygler M, Dijkstra B, Frolow F. The α/β hydrolase fold. Protein Eng. 1992;5:197–211.View ArticleGoogle Scholar
- Bornscheuer UT. Microbial carboxyl esterases: classification, properties and application in biocatalysis. FEMS Microbiol Rev. 2002;26:73–81.View ArticleGoogle Scholar
- Arpigny JL, Jaeger KE. Bacterial lipolytic enzymes: classification and properties. Biochem J. 1999;343:177–83.View ArticleGoogle Scholar
- Cheng C, Huang L, Ma R, Zhou ZQ, Diao JL. Enantioselective toxicity of lactofen and its metabolites in Scenedesmus obliquus. Algal Res. 2015;10:72–9.View ArticleGoogle Scholar
- Wang F, Liu DH, Qu H, Chen L, Zhou ZQ, Wang P. A full evaluation for the enantiomeric impacts of lactofen and its metabolites on aquatic macrophyte Lemna minor. Water Res. 2016;15:55–63.View ArticleGoogle Scholar
- Sambrook J, Russell DW. Molecular cloning: a laboratory manual. 3rd ed. New York: Cold Spring Harbor; 2001.Google Scholar
- Thompson JD, Higgins DG, Gibson TJ. CLUSTL 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:4673–80.View ArticleGoogle Scholar
- Tamura K, Stecher G, Peterson D, Filipski A, Kumar S. MEGA6: molecular evolutionary genetics analysis version 6.0. Mol Biol Evol. 2013;30:2725–9.View ArticleGoogle Scholar
- Gao RJ, Feng Y, Ishikawa K, Ishida H, Ando S, Kosugib Y, Cao SG. Cloning, purification and properties of a hyperthermophilic esterase from archaeon Aeropyrum pernix K1. J Mol Catal B Enzym. 2003;24–25:1–8.View ArticleGoogle Scholar