The database of epoxide hydrolases and haloalkane dehalogenases: one structure, many functions

Metagenomic Mining for Esterases in the Microbial Community of Los Rueldos Acid Mine Drainage Formation

Acid mine drainage (AMD) systems are extremely acidic and are metal-rich formations inhabited by relatively low-complexity communities of acidophiles whose enzymes remain mostly uncharacterized. Indeed, enzymes from only a few AMD sites have been studied. The low number of available cultured representatives and genome sequences of acidophiles inhabiting AMDs makes it difficult to assess the potential of these environments for enzyme bioprospecting. In this study, using naïve and in silico metagenomic approaches, we retrieved 16 esterases from the α/β-hydrolase fold superfamily with the closest match from uncultured acidophilic Acidobacteria, Actinobacteria (Acidithrix, Acidimicrobium, and Ferrimicrobium), Acidiphilium, and other Proteobacteria inhabiting the Los Rueldos site, which is a unique AMD formation in northwestern Spain with a pH of ∼2. Within this set, only two polypeptides showed high homology (99.4%), while for the rest, the pairwise identities ranged between 4 and 44.9%, suggesting that the diversity of active polypeptides was dominated not by a particular type of protein or highly similar clusters of proteins, but by diverse non-redundant sequences. The enzymes exhibited amino acid sequence identities ranging from 39 to 99% relative to homologous proteins in public databases, including those from other AMDs, thus indicating the potential novelty of proteins associated with a specialized acidophilic community. Ten of the 16 hydrolases were successfully expressed in Escherichia coli. The pH for optimal activity ranged from 7.0 to 9.0, with the enzymes retaining 33–68% of their activities at pH 5.5, which was consistent with the relative frequencies of acid residues (from 54 to 67%). The enzymes were the most active at 30–65°C, retaining 20–61% of their activity under the thermal conditions characterizing Los Rueldos (13.8 ± 0.6°C). The analysis of the substrate specificity revealed the capacity of six hydrolases to efficiently degrade (up to 1,652 ± 75 U/g at pH 8.0 and 30°C) acrylic- and terephthalic-like [including bis(2-hydroxyethyl)-terephthalate, BHET] esters, and these enzymes could potentially be of use for developing plastic degradation strategies yet to be explored. Our assessment uncovers the novelty and potential biotechnological interest of enzymes present in the microbial populations that inhibit the Los Rueldos AMD system.

The modular structure of α/β-hydrolases

The α/β-hydrolase fold family is highly diverse in sequence, structure and biochemical function. To investigate the sequence-structure-function relationships, the Lipase Engineering Database (https://led.biocatnet.de) was updated. Overall, 280 638 protein sequences and 1557 protein structures were analysed. All α/β-hydrolases consist of the catalytically active core domain, but they might also contain additional structural modules, resulting in 12 different architectures: core domain only, additional lids at three different positions, three different caps, additional N- or C-terminal domains and combinations of N- and C-terminal domains with caps and lids respectively. In addition, the α/β-hydrolases were distinguished by their oxyanion hole signature (GX-, GGGX- and Y-types). The N-terminal domains show two different folds, the Rossmann fold or the β-propeller fold. The C-terminal domains show a β-sandwich fold. The N-terminal β-propeller domain and the C-terminal β-sandwich domain are structurally similar to carbohydrate-binding proteins such as lectins. The classification was applied to the newly discovered polyethylene terephthalate (PET)-degrading PETases and MHETases, which are core domain α/β-hydrolases of the GX- and the GGGX-type respectively. To investigate evolutionary relationships, sequence networks were analysed. The degree distribution followed a power law with a scaling exponent γ = 1.4, indicating a highly inhomogeneous network which consists of a few hubs and a large number of less connected sequences. The hub sequences have many functional neighbours and therefore are expected to be robust toward possible deleterious effects of mutations. The cluster size distribution followed a power law with an extrapolated scaling exponent τ = 2.6, which strongly supports the connectedness of the sequence space of α/β-hydrolases. DATABASE: Supporting data about domains from other proteins with structural similarity to the N- or C-terminal domains of α/β-hydrolases are available in Data Repository of the University of Stuttgart (DaRUS) under doi: https://doi.org/10.18419/darus-458.

Epoxide hydrolysis as a model system for understanding flux through a branched reaction scheme

The epoxide hydrolase StEH1 catalyzes the hydrolysis of trans-methylstyrene oxide to 1-phenyl-propane-1,2-diol. The (S,S)-epoxide is exclusively transformed into the (1R,2S)-diol, while hydrolysis of the (R,R)-epoxide results in a mixture of product enantiomers. In order to understand the differences in the stereoconfigurations of the products, the reactions were studied kinetically during both the pre-steady-state and steady-state phases. A number of closely related StEH1 variants were analyzed in parallel, and the results were rationalized by structure-activity analysis using the available crystal structures of all tested enzyme variants. Finally, empirical valence-bond simulations were performed in order to provide additional insight into the observed kinetic behaviour and ratios of the diol product enantiomers. These combined data allow us to present a model for the flux through the catalyzed reactions. With the (R,R)-epoxide, ring opening may occur at either C atom and with similar energy barriers for hydrolysis, resulting in a mixture of diol enantiomer products. However, with the (S,S)-epoxide, although either epoxide C atom may react to form the covalent enzyme intermediate, only the pro-(R,S) alkylenzyme is amenable to subsequent hydrolysis. Previously contradictory observations from kinetics experiments as well as product ratios can therefore now be explained for this biocatalytically relevant enzyme.

Comparative Ecology of Capsular Exophiala Species Causing Disseminated Infection in Humans

Exophiala spinifera and Exophiala dermatitidis (Fungi: Chaetothyriales) are black yeast agents potentially causing disseminated infection in apparently healthy humans. They are the only Exophiala species producing extracellular polysaccharides around yeast cells. In order to gain understanding of eventual differences in intrinsic virulence of the species, their clinical profiles were compared and found to be different, suggesting pathogenic strategies rather than coincidental opportunism. Ecologically relevant factors were compared in a model set of strains of both species, and significant differences were found in clinical and environmental preferences, but virulence, tested in Galleria mellonella larvae, yielded nearly identical results. Virulence factors, i.e., melanin, capsule and muriform cells responded in opposite direction under hydrogen peroxide and temperature stress and thus were inconsistent with their hypothesized role in survival of phagocytosis. On the basis of physiological profiles, possible natural habitats of both species were extrapolated, which proved to be environmental rather than animal-associated. Using comparative genomic analyses we found differences in gene content related to lipid metabolism, cell wall modification and polysaccharide capsule production. Despite the fact that both species cause disseminated infections in apparently healthy humans, it is concluded that they are opportunists rather than pathogens.

N-terminal truncation contributed to increasing thermal stability of mannanase Man1312 without activity loss

BACKGROUND

The disordered residues on distal loops affect the molecular structural stability and on some occasions have regulatory roles in catalytic reaction. To increase understanding of the influence of distal residue mutation, this study explored the thermostability and enzymatic activity of mannanase Man1312 deletion mutants. The focus was on residues located on the N-terminal region because they are more disordered and changeable. The effects of N-terminal truncation on enzymatic activity and thermal dynamics were investigated by spectrophotometry, circular dichroism and differential scanning calorimetry assays.

RESULTS

The deletion mutants on V3, N7 and Q11 showed a marked increase in stability, while the enzymatic activity was significantly improved when triplet deletion was carried out. Triplet deletion MandVNQ showed around double the stability of its corresponding single-site and double-site deletion mutants. The Tm value of MandVNP was about 8 °C higher than that of Man1312. MandVNP had improved characteristics of Topt by 10 °C, t1/2 by 10 min and catalytic activity by 11% in comparison with Man1312. Analysis of spectra and modeling showed that MandVNQ had increased helix and strand contents.

CONCLUSION

N-terminal truncation had positive effects on the thermostability and activity of mannanase.

Estimating the success of enzyme bioprospecting through metagenomics: current status and future trends

Recent reports have suggested that the establishment of industrially relevant enzyme collections from environmental genomes has become a routine procedure. Across the studies assessed, a mean number of approximately 44 active clones were obtained in an average size of approximately 53,000 clones tested using naïve screening protocols. This number could be significantly increased in shorter times when novel metagenome enzyme sequences obtained by direct sequencing are selected and subjected to high-throughput expression for subsequent production and characterization. The pre-screening of clone libraries by naïve screens followed by the pyrosequencing of the inserts allowed for a 106-fold increase in the success rate of identifying genes encoding enzymes of interest. However, a much longer time, usually on the order of years, is needed from the time of enzyme identification to the establishment of an industrial process. If the hit frequency for the identification of enzymes performing at high turnover rates under real application conditions could be increased while still covering a high natural diversity, the very expensive and time-consuming enzyme optimization phase would likely be significantly shortened. At this point, it is important to review the current knowledge about the success of fine-tuned naïve- and sequence-based screening protocols for enzyme selection and to describe the environments worldwide that have already been subjected to enzyme screen programmes through metagenomic tools. Here, we provide such estimations and suggest the current challenges and future actions needed before environmental enzymes can be successfully introduced into the market.

Optimization of translation profiles enhances protein expression and solubility

mRNA is translated with a non-uniform speed that actively coordinates co-translational folding of protein domains. Using structure-based homology we identified the structural domains in epoxide hydrolases (EHs) and introduced slow-translating codons to delineate the translation of single domains. These changes in translation speed dramatically improved the solubility of two EHs of metagenomic origin in Escherichia coli. Conversely, the importance of transient attenuation for the folding, and consequently solubility, of EH was evidenced with a member of the EH family from Agrobacterium radiobacter, which partitions in the soluble fraction when expressed in E. coli. Synonymous substitutions of codons shaping the slow-transiting regions to fast-translating codons render this protein insoluble. Furthermore, we show that low protein yield can be enhanced by decreasing the free folding energy of the initial 5’-coding region, which can disrupt mRNA secondary structure and enhance ribosomal loading. This study provides direct experimental evidence that mRNA is not a mere messenger for translation of codons into amino acids but bears an additional layer of information for folding, solubility and expression level of the encoded protein. Furthermore, it provides a general frame on how to modulate and fine-tune gene expression of a target protein.

AromaDeg, a novel database for phylogenomics of aerobic bacterial degradation of aromatics

Understanding prokaryotic transformation of recalcitrant pollutants and the in-situ metabolic nets require the integration of massive amounts of biological data. Decades of biochemical studies together with novel next-generation sequencing data have exponentially increased information on aerobic aromatic degradation pathways. However, the majority of protein sequences in public databases have not been experimentally characterized and homology-based methods are still the most routinely used approach to assign protein function, allowing the propagation of misannotations. AromaDeg is a web-based resource targeting aerobic degradation of aromatics that comprises recently updated (September 2013) and manually curated databases constructed based on a phylogenomic approach. Grounded in phylogenetic analyses of protein sequences of key catabolic protein families and of proteins of documented function, AromaDeg allows query and data mining of novel genomic, metagenomic or metatranscriptomic data sets. Essentially, each query sequence that match a given protein family of AromaDeg is associated to a specific cluster of a given phylogenetic tree and further function annotation and/or substrate specificity may be inferred from the neighboring cluster members with experimentally validated function. This allows a detailed characterization of individual protein superfamilies as well as high-throughput functional classifications. Thus, AromaDeg addresses the deficiencies of homology-based protein function prediction, combining phylogenetic tree construction and integration of experimental data to obtain more accurate annotations of new biological data related to aerobic aromatic biodegradation pathways. We pursue in future the expansion of AromaDeg to other enzyme families involved in aromatic degradation and its regular update.

Database URL:http://aromadeg.siona.helmholtz-hzi.de

Formation of epichlorohydrin, a known rodent carcinogen, following oral administration of 1,3-dichloro-2-propanol in rats

The observed toxicity and carcinogenicity of 1,3-dichloro-2-propanol (DCP) in rodents is thought to be due to the formation of reactive metabolites, epichlorohydrin (ECH) and dichloroacetone (DCA). However, there is no direct evidence for the formation of these metabolites from exposure to DCP in rodents due to the challenges of measuring these reactive intermediates directly in vivo. The objective of this work was to investigate the metabolism of DCP to ECH and DCA in vivo by first developing a sensitive analytical method in a suitable biological matrix and analyzing samples from rats administered DCP. DCA reacted rapidly in vitro in rat blood, plasma, and liver homogenate, precluding its detection. Because ECH rapidly disappeared in liver homogenate, but was relatively long-lived in plasma and blood in vitro, blood was selected for analysis of this metabolite. Following a single oral dose of 50 mg/kg DCP in male or female Harlan Sprague-Dawley rats, ECH was detected in blood with a maximum concentration reached at ≤13.7 min. ECH was cleared rapidly with a half-life of ca. 33 and 48 min in males and females, respectively. Following a single oral dose of 25 mg/kg ECH in male and female rats, the elimination half-life of ECH was ca. 34 and 20 min, respectively; the oral bioavailability of ECH was low (males, 5.2%; females, 2.1%), suggesting extensive first pass metabolism of ECH following oral administration. The area under the concentration vs time curve for ECH following oral administration of DCP and intravenous administration of ECH was used to estimate the percent of the DCP dose converted to ECH in rats. On the basis of this analysis, we concluded that in male and female rats following oral administration of 50 mg/kg DCP, ≥1.26% or ≥1.78% of the administered dose was metabolized to ECH, respectively.

Compositional profile of α / β-hydrolase fold proteins in mangrove soil metagenomes: prevalence of epoxide hydrolases and haloalkane dehalogenases in oil-contaminated sites

The occurrence of genes encoding biotechnologically relevant α/β-hydrolases in mangrove soil microbial communities was assessed using data obtained by whole-metagenome sequencing of four mangroves areas, denoted BrMgv01 to BrMgv04, in São Paulo, Brazil. The sequences (215 Mb in total) were filtered based on local amino acid alignments against the Lipase Engineering Database. In total, 5923 unassembled sequences were affiliated with 30 different α/β-hydrolase fold superfamilies. The most abundant predicted proteins encompassed cytosolic hydrolases (abH08; ∼ 23%), microsomal hydrolases (abH09; ∼ 12%) and Moraxella lipase-like proteins (abH04 and abH01; < 5%). Detailed analysis of the genes predicted to encode proteins of the abH08 superfamily revealed a high proportion related to epoxide hydrolases and haloalkane dehalogenases in polluted mangroves BrMgv01-02-03. This suggested selection and putative involvement in local degradation/detoxification of the pollutants. Seven sequences that were annotated as genes for putative epoxide hydrolases and five for putative haloalkane dehalogenases were found in a fosmid library generated from BrMgv02 DNA. The latter enzymes were predicted to belong to Actinobacteria, Deinococcus-Thermus, Planctomycetes and Proteobacteria. Our integrated approach thus identified 12 genes (complete and/or partial) that may encode hitherto undescribed enzymes. The low amino acid identity (< 60%) with already-described genes opens perspectives for both production in an expression host and genetic screening of metagenomes.

New generation of biocatalysts for organic synthesis

The use of enzymes as catalysts for the preparation of novel compounds has received steadily increasing attention over the past few years. High demands are placed on the identification of new biocatalysts for organic synthesis. The catalysis of more ambitious reactions reflects the high expectations of this field of research. Enzymes play an increasingly important role as biocatalysts in the synthesis of key intermediates for the pharmaceutical and chemical industry, and new enzymatic technologies and processes have been established. Enzymes are an important part of the spectrum of catalysts available for synthetic chemistry. The advantages and applications of the most recent and attractive biocatalysts--reductases, transaminases, ammonia lyases, epoxide hydrolases, and dehalogenases--will be discussed herein and exemplified by the syntheses of interesting compounds.

The lid domain of the MCP hydrolase DxnB2 contributes to the reactivity toward recalcitrant PCB metabolites

DxnB2 and BphD are meta-cleavage product (MCP) hydrolases that catalyze C-C bond hydrolysis of the biphenyl metabolite 2-hydroxy-6-oxo-6-phenylhexa-2,4-dienoic acid (HOPDA). BphD is a bottleneck in the bacterial degradation of polychlorinated biphenyls (PCBs) by the Bph catabolic pathway due in part to inhibition by 3-Cl HOPDAs. By contrast, DxnB2 from Sphingomonas wittichii RW1 catalyzes the hydrolysis of 3-Cl HOPDAs more efficiently. X-ray crystallographic studies of the catalytically inactive S105A variant of DxnB2 complexed with 3-Cl HOPDA revealed a binding mode in which C1 through C6 of the dienoate are coplanar. The chlorine substituent is accommodated by a hydrophobic pocket that is larger than the homologous site in BphDLB400 from Burkholderia xenovorans LB400. The planar binding mode observed in the crystalline complex was consistent with the hyper- and hypsochromically shifted absorption spectra of 3-Cl and 3,9,11-triCl HOPDA, respectively, bound to S105A in solution. Moreover, ES(red), an intermediate possessing a bathochromically shifted spectrum observed in the turnover of HOPDA, was not detected, suggesting that substrate destabilization was rate-limiting in the turnover of these PCB metabolites. Interestingly, electron density for the first α-helix of the lid domain was poorly defined in the dimeric DxnB2 structures, unlike in the tetrameric BphDLB400. Structural comparison of MCP hydrolases identified the NC-loop, connecting the lid to the α/β-hydrolase core domain, as a determinant in the oligomeric state and suggests its involvement in catalysis. Finally, an increased mobility of the DxnB2 lid may contribute to the enzyme's ability to hydrolyze PCB metabolites, highlighting how lid architecture contributes to substrate specificity in α/β-hydrolases.

Computer-Aided Protein Directed Evolution: a Review of Web Servers, Databases and other Computational Tools for Protein Engineering

The combination of computational and directed evolution methods has proven a winning strategy for protein engineering. We refer to this approach as computer-aided protein directed evolution (CAPDE) and the review summarizes the recent developments in this rapidly growing field. We will restrict ourselves to overview the availability, usability and limitations of web servers, databases and other computational tools proposed in the last five years. The goal of this review is to provide concise information about currently available computational resources to assist the design of directed evolution based protein engineering experiment.

Structural classification by the Lipase Engineering Database: a case study of Candida antarctica lipase A

Background

The Lipase Engineering Database (LED) integrates information on sequence, structure and function of lipases, esterases and related proteins with the α/β hydrolase fold. A new superfamily for Candida antarctica lipase A (CALA) was introduced including the recently published crystal structure of CALA. Since CALA has a highly divergent sequence in comparison to other α/β hydrolases, the Lipase Engineering Database was used to classify CALA in the frame of the already established classification system. This involved the comparison of CALA to similar structures as well as sequence-based comparisons against the content of the LED.

Results

The new release 3.0 (December 2009) of the Lipase Engineering Database contains 24783 sequence entries for 18585 proteins as well as 656 experimentally determined protein structures, including the structure of CALA. In comparison to the previous release [1] with 4322 protein and 167 structure entries this update represents a significant increase in data volume. By comparing CALA to representative structures from all superfamilies, a structure from the deacetylase superfamily was found to be most similar to the structure of CALA. While the α/β hydrolase fold is conserved in both proteins, the major difference is found in the cap region. Sequence alignments between both proteins show a sequence similarity of only 15%. A multisequence alignment of both protein families was used to create hidden Markov models for the cap region of CALA and showed that the cap region of CALA is unique among all other proteins of the α/β hydrolase fold. By specifically comparing the substrate binding pocket of CALA to other binding pockets of α/β hydrolases, the binding pocket of Candida rugosa lipase was identified as being highly similar. This similarity also applied to the lid of Candida rugosa lipase in comparison to the potential lid of CALA.

Conclusion

The LED serves as a valuable tool for the systematic analysis of single proteins or protein families. The updated release 3.0 was used for the evaluation of α/β hydrolases. The HTML version of the database with new features is available at http://www.led.uni-stuttgart.de and provides sequences, structures and a set of analysis tools including phylogenetic trees and HMM profiles

Characterization of the SgcF epoxide hydrolase supporting an (R)-vicinal diol intermediate for enediyne antitumor antibiotic C-1027 biosynthesis

C-1027 is a chromoprotein antitumor antibiotic consisting of an apoprotein and the C-1027 chromophore. The C-1027 chromophore possesses four distinct structural moieties-an enediyne core, a deoxy aminosugar, a benzoxazolinate, and an (S)-3-chloro-5-hydroxy-beta-tyrosine-the latter two of which are proposed to be appended to the enediyne core via a convergent biosynthetic strategy. Here we report the in vitro characterization of SgcF, an epoxide hydrolase from the C-1027 biosynthetic gene cluster that catalyzes regio- and stereospecific hydrolysis of styrene oxide, serving as an enediyne core epoxide intermediate mimic, to form a vicinal diol. Abolishment of C-1027 production in the DeltasgcF mutant strain Streptomyces globisporus SB1010 unambiguously establishes that sgcF plays an indispensable role in C-1027 biosynthesis. SgcF efficiently hydrolyzes (S)-styrene oxide, displaying an apparent K(m) of 0.6 +/- 0.1 mM and k(cat) of 48 +/- 1 min(-1), via attack at the alpha-position to exclusively generate the (R)-phenyl vicinal diol, consistent with the stereochemistry of the C-1027 chromophore. These findings support the role of SgcF in the proposed convergent pathway for C-1027 biosynthesis, unveiling an (R)-vicinal diol as a key intermediate. Interestingly, SgcF can also hydrolyze (R)-styrene oxide to afford preferentially the (R)-phenyl vicinal diol via attack at the beta-position, albeit with significantly reduced efficiency (apparent K(m) of 2.0 +/- 0.4 mM and k(cat) = 4.3 +/- 0.3 min(-1)). Although the latter activity unlikely contributes to C-1027 biosynthesis in vivo, such enantioconvergence arising from complementary regioselective hydrolysis of a racemic substrate could be exploited to engineer epoxide hydrolases with improved regio- and/or enantiospecificity.

Effect of an exponential feeding regime on the production of Rhodotorula araucariae epoxide hydrolase in Yarrowia lipolytica

AIMS

To evaluate the effect of and exponential feeding regime on the production of epoxide hydrolase (EH) enzyme in recombinant Yarrowia lipolytica in comparison to a constant feed strategy.

METHODS AND RESULTS

An exponential feed model was developed and fermentations were fed at six different exponential rates. A twofold increase in EH productivity and a 15% increase in volumetric EH activity was obtained by applying exponential glucose feed rates in fed-batch cultivation. These responses were modelled to obtain a theoretical optimum feed rate that was validated in duplicate fermentations. The model optimum of 0.06 h(-1) resulted in a volumetric EH activity of c. 5500 U l(-1) h(-1) and a maximum activity of 206,000 U l(-1). This correlated well with model predictions, with a variance of <10%.

CONCLUSIONS

The use of an exponential feed strategy at a rate of 0.06 h(-1) yielded best results for all key responses which show a clear improvement over a constant feed strategy.

SIGNIFICANCE AND IMPACT OF THE STUDY

The study was the first evaluation of an exponential feed strategy on recombinant Y. lipolytica for the production of EH enzyme. The results suggest a strategy for the commercial production of a valuable pharmaceutical enzyme.

Molecular engineering of epoxide hydrolase and its application to asymmetric and enantioconvergent hydrolysis

Safety and regulatory issues favor increasing use of enantiopure compounds in pharmaceuticals. Enantiopure epoxides and diols are valuable intermediates in organic synthesis for the production of optically active pharmaceuticals. Enantiopure epoxide can be prepared using epoxide hydrolase (EH)-catalyzed asymmetric hydrolysis of its racemate. Enantioconvergent hydrolysis of racemic epoxides by EHs possessing complementary enantioselectivity and regioselectivity can lead to the formation of enantiopure vicinal diols with high yield. EHs are cofactor-independent and easy-to-use catalysts. EHs will attract much attention as commercial biocatalysts for the preparation of enantiopure epoxides and diols. In this paper, recent progress in molecular engineering of EHs is reviewed. Some examples and prospects of asymmetric and enantioconvergent hydrolysis reactions are discussed as supplements to molecular engineering to improve EH performance.

The chemical defensome: environmental sensing and response genes in the Strongylocentrotus purpuratus genome

Metazoan genomes contain large numbers of genes that participate in responses to environmental stressors. We surveyed the sea urchin Strongylocentrotus purpuratus genome for homologs of gene families thought to protect against chemical stressors; these genes collectively comprise the 'chemical defensome.' Chemical defense genes include cytochromes P450 and other oxidases, various conjugating enzymes, ATP-dependent efflux transporters, oxidative detoxification proteins, and transcription factors that regulate these genes. Together such genes account for more than 400 genes in the sea urchin genome. The transcription factors include homologs of the aryl hydrocarbon receptor, hypoxia-inducible factor, nuclear factor erythroid-derived 2, heat shock factor, and nuclear hormone receptors, which regulate stress-response genes in vertebrates. Some defense gene families, including the ABCC, the UGT, and the CYP families, have undergone expansion in the urchin relative to other deuterostome genomes, whereas the stress sensor gene families do not show such expansion. More than half of the defense genes are expressed during embryonic or larval life stages, indicating their importance during development. This genome-wide survey of chemical defense genes in the sea urchin reveals evolutionary conservation of this network combined with lineage-specific diversification that together suggest the importance of these chemical stress sensing and response mechanisms in early deuterostomes. These results should facilitate future studies on the evolution of chemical defense gene networks and the role of these networks in protecting embryos from chemical stress during development.

Diversity and biocatalytic potential of epoxide hydrolases identified by genome analysis

Epoxide hydrolases play an important role in the biodegradation of organic compounds and are potentially useful in enantioselective biocatalysis. An analysis of various genomic databases revealed that about 20% of sequenced organisms contain one or more putative epoxide hydrolase genes. They were found in all domains of life, and many fungi and actinobacteria contain several putative epoxide hydrolase-encoding genes. Multiple sequence alignments of epoxide hydrolases with other known and putative alpha/beta-hydrolase fold enzymes that possess a nucleophilic aspartate revealed that these enzymes can be classified into eight phylogenetic groups that all contain putative epoxide hydrolases. To determine their catalytic activities, 10 putative bacterial epoxide hydrolase genes and 2 known bacterial epoxide hydrolase genes were cloned and overexpressed in Escherichia coli. The production of active enzyme was strongly improved by fusion to the maltose binding protein (MalE), which prevented inclusion body formation and facilitated protein purification. Eight of the 12 fusion proteins were active toward one or more of the 21 epoxides that were tested, and they converted both terminal and nonterminal epoxides. Four of the new epoxide hydrolases showed an uncommon enantiopreference for meso-epoxides and/or terminal aromatic epoxides, which made them suitable for the production of enantiopure (S,S)-diols and (R)-epoxides. The results show that the expression of epoxide hydrolase genes that are detected by analyses of genomic databases is a useful strategy for obtaining new biocatalysts.

EHPred: an SVM-based method for epoxide hydrolases recognition and classification

A two-layer method based on support vector machines (SVMs) has been developed to distinguish epoxide hydrolases (EHs) from other enzymes and to classify its subfamilies using its primary protein sequences. SVM classifiers were built using three different feature vectors extracted from the primary sequence of EHs: the amino acid composition (AAC), the dipeptide composition (DPC), and the pseudo-amino acid composition (PAAC). Validated by 5-fold cross tests, the first layer SVM classifier can differentiate EHs and non-EHs with an accuracy of 94.2% and has a Matthew's correlation coefficient (MCC) of 0.84. Using 2-fold cross validation, PAAC-based second layer SVM can further classify EH subfamilies with an overall accuracy of 90.7% and MCC of 0.87 as compared to AAC (80.0%) and DPC (84.9%). A program called EHPred has also been developed to assist readers to recognize EHs and to classify their subfamilies using primary protein sequences with greater accuracy.

Cloning, biochemical properties, and distribution of mycobacterial haloalkane dehalogenases

Haloalkane dehalogenases are enzymes that catalyze the cleavage of the carbon-halogen bond by a hydrolytic mechanism. Genomes of Mycobacterium tuberculosis and M. bovis contain at least two open reading frames coding for the polypeptides showing a high sequence similarity with biochemically characterized haloalkane dehalogenases. We describe here the cloning of the haloalkane dehalogenase genes dmbA and dmbB from M. bovis 5033/66 and demonstrate the dehalogenase activity of their translation products. Both of these genes are widely distributed among species of the M. tuberculosis complex, including M. bovis, M. bovis BCG, M. africanum, M. caprae, M. microti, and M. pinnipedii, as shown by the PCR screening of 48 isolates from various hosts. DmbA and DmbB proteins were heterologously expressed in Escherichia coli and purified to homogeneity. The DmbB protein had to be expressed in a fusion with thioredoxin to obtain a soluble protein sample. The temperature optimum of DmbA and DmbB proteins determined with 1,2-dibromoethane is 45 degrees C. The melting temperature assessed by circular dichroism spectroscopy of DmbA is 47 degrees C and DmbB is 57 degrees C. The pH optimum of DmbA depends on composition of a buffer with maximal activity at 9.0. DmbB had a single pH optimum at pH 6.5. Mycobacteria are currently the only genus known to carry more than one haloalkane dehalogenase gene, although putative haloalkane dehalogenases can be inferred in more then 20 different bacterial species by comparative genomics. The evolution and distribution of haloalkane dehalogenases among mycobacteria is discussed.

How to find soluble proteins: a comprehensive analysis of alpha/beta hydrolases for recombinant expression in E. coli

Background

In screening of libraries derived by expression cloning, expression of active proteins in E. coli can be limited by formation of inclusion bodies. In these cases it would be desirable to enrich gene libraries for coding sequences with soluble gene products in E. coli and thus to improve the efficiency of screening. Previously Wilkinson and Harrison showed that solubility can be predicted from amino acid composition (Biotechnology 1991, 9(5):443–448). We have applied this analysis to members of the alpha/beta hydrolase fold family to predict their solubility in E. coli. alpha/beta hydrolases are a highly diverse family with more than 1800 proteins which have been grouped into homologous families and superfamilies.

Results

The predicted solubility in E. coli depends on hydrolase size, phylogenetic origin of the host organism, the homologous family and the superfamily, to which the hydrolase belongs. In general small hydrolases are predicted to be more soluble than large hydrolases, and eukaryotic hydrolases are predicted to be less soluble in E. coli than prokaryotic ones. However, combining phylogenetic origin and size leads to more complex conclusions. Hydrolases from prokaryotic, fungal and metazoan origin are predicted to be most soluble if they are of small, medium and large size, respectively. We observed large variations of predicted solubility between hydrolases from different homologous families and from different taxa.

Conclusion

A comprehensive analysis of all alpha/beta hydrolase sequences allows more efficient screenings for new soluble alpha/beta hydrolases by the use of libraries which contain more soluble gene products. Screening of hydrolases from families whose members are hard to express as soluble proteins in E. coli should first be done in coding sequences of organisms from phylogenetic groups with the highest average of predicted solubility for proteins of this family. The tools developed here can be used to identify attractive target genes for expression using protein sequences published in databases. This analysis also directs the design of degenerate, family- specific primers to amplify new members from homologous families or superfamilies with a high probability of soluble alpha/beta hydrolases.

Cloning, recombinant expression and biochemical characterisation of novel esterases from Bacillus sp. associated with the marine sponge Aplysina aerophoba

Two novel esterases (EstB1 and EstB2) were isolated from a genomic library of Bacillus sp. associated with the marine sponge Aplysina aerophoba. EstB1 shows low identity (26-44%) with the published hydrolases of the genus Bacillus, whereas EstB2 shows high identity (73-74%) with the carboxylesterases from B. cereus and B. anthracis. Both esterases were efficiently expressed in Escherichia coli under the control of T7 promoter using the vector pET-22b(+). Recombinant EstB1 was purified in a single step to electrophoretic homogeneity by IMAC. A method for the refolding of inclusion bodies formed by the recombinant EstB2 was established to obtain active enzyme. Substrate specificity of the two enzymes towards p-nitrophenyl and methyl esters and the respective kinetic parameters K(m) and V(max) were determined. The temperature optima of EstB1 and EstB2 were determined to be in the range of 30-50 degrees C and 20-35 degrees C, respectively. The pH optima were found to be in the range of 6.5-7.5 and 6.5-8.0, respectively. Both enzymes showed the highest stability in up to 50% (v/v) DMSO followed by methanol, ethanol and 2-propanol. The influence of high NaCl and KCl concentrations was tested. The inhibition effect of 10-50 mM Zn(2+) and 50 mM Mg(2+) and Ca(2+) ions was observed for both esterases. One to five millimolar PMSF deactivated the enzymes, whereas beta-mercaptoethanol, DTT and EDTA had no effect on the enzymes activity.