Inter‐conversion of catalytic abilities in a bifunctional carboxyl/feruloyl‐esterase from earthworm gut metagenome

Summary Carboxyl esterases (CE) exhibit various reaction specificities despite of their overall structural similarity. In present study we have exploited functional metagenomics, saturation mutagenesis and experimental protein evolution to explore residues that have a significant role in substrate discrimination. We used an enzyme, designated 3A6, derived from the earthworm gut metagenome that exhibits CE and feruloyl esterase (FAE) activities with p‐nitrophenyl and cinnamate esters, respectively, with a [(kcat/Km)]CE/[(kcat/Km)]FAE factor of 17. Modelling‐guided saturation mutagenesis at specific hotspots (Lys281, Asp282, Asn316 and Lys317) situated close to the catalytic core (Ser143/Asp273/His305) and a deletion of a 34‐AA–long peptide fragment yielded mutants with the highest CE activity, while cinnamate ester bond hydrolysis was effectively abolished. Although, single to triple mutants with both improved activities (up to 180‐fold in kcat/Km values) and enzymes with inverted specificity ((kcat/Km)CE/(kcat/Km)FAE ratio of ∼0.4) were identified, no CE inactive variant was found. Screening of a large error‐prone PCR‐generated library yielded by far less mutants for substrate discrimination. We also found that no significant changes in CE activation energy occurs after any mutation (7.3 to −5.6 J mol−1), whereas a direct correlation between loss/gain of FAE function and activation energies (from 33.05 to −13.7 J mol−1) was found. Results suggest that the FAE activity in 3A6 may have evolved via introduction of a limited number of ‘hot spot’ mutations in a common CE ancestor, which may retain the original hydrolytic activity due to lower restrictive energy barriers but conveys a dynamic energetically favourable switch of a second hydrolytic reaction.


Introduction
Carboxyl-esterases (CEs; EC 3.1.1.X) are ubiquitous a/bhydrolases with representatives in all three domains of life, Eukarya, Bacteria and Archaea (Bornscheuer, 2002). These proteins catalyse the cleavage and formation of ester bonds whose components differ in chain nature and length. Within the a/b-hydrolase superfamily, CEs belong to the best structurally and functionally characterized, with more than 6016 members divided in 90 subfamilies after a comprehensive non-redundant search in the nucleotide collections, reference genomic sequences, whole genome shotgun reads and environmental samples databases. Site-directed mutagenesis (Manco et al., 2001;Reyes-Duarte et al., 2005), saturation mutagenesis (Mee-Hie Cho et al., 2006;Wang et al., 2006;Nakagawa et al., 2007) and a number of directed evolution methods (Funke et al., 2005;Ivancic et al., 2007) have successfully been applied to identify the molecular determinants of substrate specificity and enantio-selectivity and to gain some knowledge on molecular evolution of this super-family of proteins.
We focused in this study on esterases able to attack the ester bond between hydroxyl-cinnamic acids, such as ferulic, p-coumaric and sinnapic acid, and sugars present in intricate structure of the plant cell wall (Fazary and Ju, 2007). The hydrolysis of these compounds is catalysed by feruloyl esterases (FAEs; EC 3.1.1.73). About 30 FAEs have been described and based on structural and functional factors a division into four types (A-D) has been established (Crepin et al., 2004;Wong, 2006;Benoit et al., 2008). They have been purified and characterized from bacteria (Wang et al., 2004;Aurilia et al., 2007; and fungi (Shin and Chen, 2007;Kanauchi et al., 2008;Moukouli et al., 2008;Koseki et al., 2009), and showed significant variations in activities. Despite all of them having the typical a/b fold and the common CE-like catalytic triad, Ser-Asp/Glu-His (Fazary and Ju, 2007), little is known about the factors determining substrate specificity of these enzymes compared with 'common' CEs, which has inter alia motivated present work. Only few studies have demonstrated that few mutations at the ligand-binding cavity (Tarbouriech et al., 2005) or in insertions such as lids (Hermoso et al., 2004;Faulds et al., 2005;Koseki et al., 2005) have important effects in controlling the binding and hydrolysis of esters with m-methoxy groups. Following on from this, Levasseur's group studied the evolutionary relationships between fungal lipases and FAEs through the study of their phylogenies (Levasseur et al., 2006). Authors concluded that lipases appear to be primitive functions from which FAE functionality could have been derived bought duplication events probably induced by environmental stresses, i.e. those caused by new substrates. Although, authors suggest that modifications within the active site might thus be positively selected, causing the shift in functionality at the organismal level (Hermoso et al., 2004;Levasseur et al., 2006), no experimental evidence for such hypothesis has been provided.
In this work, the first systematic study was performed to investigate the structural and energetic determinants for the transition between CEs and FAEs. Here, we used a novel bi-functional CE/FAE protein (3A6) that catalyses hydrolytic cleavage of 'common' and 'cinnamate' esters, isolated from the microbial community of the gut of the earthworm Aporrectodea caliginosa. Invertebrate guts are certainly one of the most diverse environments yet studied and, in particular, the earthworm gut associated-organisms require an efficient enzymatic machinery to cope with a continuous flux of various polymeric materials (including cinnamate esters) from the soil and from plant litter passing through intestinal tract, which make it attractive for isolating hydrolytic activities. We show here by examining saturation mutagenesis and error-prone PCR libraries, that a reduced number of 'hot-spot' mutations could alter the overall substrate hydrolysis and, more importantly, evidence is presented for the first time that suggests that the acquisition of FAE phenotype, from an ancestor CE activity, is energetically favourable. Moreover, to the best of our knowledge, this work is the first example for a successful conversion of an FAE into a common CE.

General characteristics of the 3A6 protein
The fosmid carrying the 3A6 gene was isolated by screening a metagenomic library from the cellulose-enriched microbial community from the gut of the earthworm A. caliginosa, based on its ability to hydrolyse methyl ferulate (MF). The analysis of 3081-bp-long subcloned DNA fragment (GC-content 42.2%; Fig. S1) revealed the presence of two oppositely oriented open reading frames (ORFs), the first of which (1026 bp, positions 1144 to 119) encodes a putative 341-amino-acid protein with a predicted molecular mass of 37 502 Da that exhibited high homology (54% identity and 71% of similarity) to proteins of the esterase/lipase superfamily (pfam00135) ( Table S1). The 3A6 gene was subcloned into a plasmid vector and after that, inserted into the pET41 Ek/LIC vector for expression in Escherichia coli. The protein, purified to homogeneity as a single monomer of 38 000 Ϯ 3200 Da, was tested for its catalytic ability to hydrolyse a series of commercially available substrates (Table 1). We show that besides   hydrolysing p-nitrophenyl (pNP) esters, 3A6 efficiently catalysed the hydrolysis of cinnamates with a (kcat/Km) factor of~17:1 with pNP acetate (pNPC2) and MF as substrates. Overall, the enzyme was able to cleavage most efficiently methyl sinapinate, followed by MF (2.5-fold), and methyl-p-coumarate (16-fold), but not methyl caffeate. As judged by the kcat/Km values, the enzyme functions with short chain substrates from 2-to 1500-fold better than the longer substrates. The enzyme was also able to hydrolyse, to some extent, p-nitrophenyl 5-O-trans-feruloyl-a-Larabinofuranoside (kcat of 1 s -1 ) and (5-O-(trans-feruloyl)a-L-arabinofuranosyl)-(1,3)-b-D-xylopyranosyl-(1,4)-Dxylopyranose (kcat of 10 s -1 ). According to its substrate specificity, the enzyme was classified as an FAE type A (Crepin et al., 2004). The protein showed maximal activity at 45-50°C and pH 7.5-8.5 (Fig. S2). The substrate specificity showed strikingly different results with those reported in literature: nearly all characterized a/b hydrolases that share sequence similarity with 3A6 cleaved p-nitrophenyl esters but not cinnamates (Table S1). Following on from this, one of the main objectives of this study was to complement the kinetic data with mechanistic and engineering protocols so as to facilitate examination the essential residues that might hinder the recognition of cinnamate esters by steric or catalytic effects in the 3A6 protein.

Saturation mutagenesis strategy for altering the enzyme's specificity
A model of the enzyme was obtained by homology using the crystal structure of carboxyl-esterase Este1 of metagenomic origin (PDB 2C7B) as a template (Fig. S3), given the high level of identity (38% of the overall sequence, close to 70% in the active-site region). The quality of the obtained model was judged to be adequate for our purposes as judged by the Ramachandran plots (not shown), and a putative active-site pocket was identified (Ser 143 /Asp 273 /His 305 ) and further confirmed by mutagenic analysis: mutations of those residues by Gly, Gln and Asn, in the same order, results in a dramatic effect on the kcat/Km for pNPC2 and MF cleavage (> 10 000-fold). To further identify key residues that might define additional stabilization, subtraction or orientation effects on substrate preference, a detailed multiple (structural) alignment of 3A6 with homologous proteins was performed (Figs S4 and S5). The attention was focused on identifying regions that might have a fundamental role in substrate specificity and preference and that are not found in homologue proteins. The regions that were selected for mutagenesis were: Asn 109 -Leu 111 , Lys 281 -Asp 282 and Tyr 315 -Lys 317 (Fig. S3).
Saturation mutagenesis libraries (Asn 109 X, Glu 110 X, Leu 111 X, Lys 281 X, Asp 282 X, Tyr 315 X, Asn 316 X and Lys 317 X) were initially constructed, and crude cell extracts from 100 colonies (found to be sufficient to ensure that the 20 AA were represented) from each library were screened for activity with both pNPC 2 and MF. The results corrected for cell growth level, are shown in Fig. 1 (mutations at positions 109-111 did not affect hydrolytic proficiency and are not shown). In each case, a minor proportion of the library members were inactive with both substrates, showing that only few amino acid substitutions were not tolerated for activity. For each library some of the variants showed altered improved activity but no change in substrate preference. However, a number of particularly interesting variants were identified in the K 281 X, D 282 X, N 316 X and K 317 X libraries, suggesting that they may play important roles in controlling substrate recognition. It is particularly interesting to note that some of the subpopulations of variants in the libraries appeared to have a significant preference for the pNP substrate, with some of N 316X and K 317 X variants having no measurable activity against MF, whereas some (i.e. K 281 X and D 282 X) showed higher overall activity for the MF screening substrate. Interestingly, pNPC2 inactive variants were not detected, suggesting that CE activity is more robust to amino acid substitutions. Detailed kinetic characterization of fourteen purified enzymes selected on the basis of their substrate preferences was carried out in order to accurately determine the importance of each of the identified amino acid substitutions in the discrimination between the pNPC2 and MF screening substrates ( Table 2).

Construction of a CE-selective enzyme.
Screening of the saturation mutagenesis libraries ( Fig. 1) identified three substitutions (populations B, D and E) that yield pNPC2 selective catalysts. In the initial screen, members of these populations showed similar or higher activity with the pNPC2 substrate than the wild-type enzyme, but most importantly for our selection criteria, those appeared to have little or no measurable activity with MF. Sequence analysis revealed that three simple amino acid substitutions (D 282 L, N 316 STOP and K 317 H) are responsible for the excellent substrate discrimination: from (kcat/Km)pNPC2/(kcat/Km)MF of~17 for the wild-type to 69 (for D 282 L) to 3735 (for N 316 STOP) and 2 ¥ 10 6 (for K 317 H). This result can be explained by the up to 10-fold greater Km values coupled with a significant (from 1.5-to 29 000-fold) reduction in the kcat values for MF. This result is entirely consistent with the apparent~3 to > 700-fold discrimination displayed by populations B, D and E in the initial screen ( Fig. 1). Our data suggest that neither of the Y 315 X variants were individually responsible for discrimination of esterase and feruloyl esterase activities, although mutations at this residue equally increased both activities.   Construction of FAE-selective enzymes. As shown in Fig. 1, populations A, C and F were identified as leading to the most MF-selective variants. Sequencing identified 10 new mutations. The K 281 I mutant displayed the highest (34-fold increase) hydrolytic efficiency for MF, mainly due to a 14-fold reduction in binding capacity (Table 2). Overall, K 281 I substitution produced the only enzyme variant more selective towards MF with a (kcat/Km)pNPC2/(kcat/Km)MF ratio of 0.7 (24 times lower than wild-type), which is in a good agreement with the apparent 16-fold discrimination displayed by population A in the screen (Fig. 1). For the other variants (K 281 N, K 281 T, K 281 S, D 282 E, N 316 L, K 317 N, K 317 G, K 317 L and K 317 D) there is a large combined effect of lower kcat/Km for pNP esters and high kcat/Km for MF that led to (kcat/Km)pNPC2/ (kcat/Km)MF ratios ranging from 4-to 14-fold lower than wild-type. Overall, the cleavage of MF was mostly positively affected by the mutations.

Gain-of-function by site-directed mutagenesis
In an attempt to improve the activity and selectivity levels a number of double and triple mutants were created by Quick-Change site-directed mutagenesis and expressed recombinantly in E. coli. As shown in the

Maximizing FAE activity of 3A6 by deletion
Although single mutations in key regions were sufficient enough to significantly alter hydrolytic activity and control substrate specificity, we attempted to delete the Phe 11 -Lys 29 and Gly 178 -Gly 211 regions, to check for their catalytic influence. Both regions, not found in homologous proteins, were selected based on structural alignments (Figs S3-S5). The first one is located on the N-terminus of the protein and its deletion caused a complete loss of enzyme activity and will not be discussed further. The second insert is located on a long loop with extremely low sequence conservation. Circular dichroism (CD) analysis of pure protein, named 3A6-I, revealed that the deletion variant exhibited a slightly modified CD spectrum com-pared with the 3A6 wild-type protein (Fig. S6) and moreover, it showed an excellent substrate discrimination (Table 1). Interestingly, compared with parental and saturation mutants, the 3A6-I protein was the only enzyme variant, which appeared to demonstrate a significant, 13-fold preference for longer substrates, with pNP butyrate being the preferred substrate over pNPC 2 (Table 1). This variant showed an unexpected residual activity towards the cinnamates, with a (kcat/Km)pNP/ (kcat/Km)MF ratio higher than 10 7 . It is therefore a CE-selective enzyme. Even though we do not know the exact role in vivo of this extension, it is clear that it has a substrate recognition role.

Substrate discrimination by random mutagenesis
The mutational strategy has therefore produced a number of enzymes with high activity levels and We first observed that the Km value for pNPC2 was 12-fold higher than that of the native enzyme, while maintaining that for MF (Table 1). Furthermore, mutant showed a 39-fold increase of (kcat/Km)pNP/(kcat/Km)MF ratio (from 17 to 1200). This result is explained by the~14-fold lower and~70fold higher kcat values for MF and pNPC2 respectively. These data suggest that there is a strong bias towards CE phenotype in the A8P4 mutant.
Above results unambiguously confirmed that mutations at specific hot spots of the polypeptide (i.e. K 281 , D 282 , N 316 and K 317 ) may play a more effective supportive role in substrate discrimination in the parental 3A6 protein com-pared with random mutations within the whole protein length, since only one mutant with substrate discrimination ability was identified by epPCR. Figure 2 summarized the effect of single to triple mutations and insert deletion on substrate preference.

Substrate preference and its association with activation energies
The experimental work on the in vitro evolution of CE and FAE activities was followed by calculating the differential binding (DDGB) and activation (DDGA) free energy profiles for the cleavage of pNPC2 and MF, respectively, by the 'parental' and 'mutant' 3A6 variants. Difference energy diagrams (energy of the variant minus energy of the parent 3A6) are shown in Fig. 3A and B, and were used to understand the thermodynamics of the changes in enzyme reaction specificity, since the changes in substrate hydrolysis and preference are brought about by the relative changes in the free energies of the transition-state barriers. First, we see that the changes in DDGB and DDGA from the native to double or triple mutants are not the sum of free energies of the single saturation mutants, either for binding or catalysis. Thus, the contribution of K 281 I, D 282 E, D 282 L, N 316 STOP, K 317 H and K 317 D is synergic rather than mutually independent for the cleavage of both substrates, each of them positively or negatively affecting free energy values.
The effects on mutants on DDGB and DDGA are rather complex. As shown, major differences are observed for the activation free energies of MF compared with that shown for pNPC2 ( Fig. 3A and B). This suggests that FAE phenotype is more mutation affected by introducing single point mutations in specific 'hot spot' residues (see also Table 2). This was further confirmed by plotting the DDGA for pNPC2 and MF against the (kcat/Km)pNPC2/(kcat/Km)MF ratio for each of the variants. Surprisingly, whereas the DDGA for pNPC2 do not varied significantly (-5.6 to 7.3 J mol -1 ), that of MF do with a sharp increase (from -13.7 to 33.1 J mol -1 ) at increasing the CE phenotype (Fig. 3C). Similar situation was found when analysing the DDGB, although at lower level (Fig. 3D). Data suggest that mutations producing FAE-selective enzymes, i.e. those with (kcat/Km)pNP/(kcat/Km)MF ratio lower than 17, may be more favourable energetically (decrease in DDG). The low dependence of pNPC2-associated free energy with cata-

Fig. 2.
Schematic representation of the complementary substrate preference of 3A6-like variants. The preferred substrate specificity of wild-type and variants are shown as the logarithm of the ratio of catalytic efficiencies (kcat/Km) towards pNPC2 and MF substrates. A value of zero thus represents an enzyme, which is non-discriminatory between the two substrates.
lytic efficiency is consistent with the fact that any of the mutations abolished the CE activity during the screening library tests (see Fig. 1 and Table 2).

Discussion
A detailed analysis of a newly identified FAE enzyme of the type A from metagenomic library is presented here. The recombinant 3A6 enzyme appeared to be a monomer of~37 kDa. Besides hydrolysing ester bonds, 3A6 efficiently catalysed the hydrolysis of pNP and cinnamate esters with a kcat/Km factor of~17. The ability to cleave both common and feruloylated esters is significantly different to that shown by homologous proteins known in bibliography and databases: those enzymes behave as common esterases and the capacity to hydrolyse cinnamates has not been reported in any of them. Here, creation of a number of complementary CE and FAEselective variants with improved activity phenotypes was accomplished using a combination of modelling-guided saturation mutagenesis and site-directed mutagenesis. Specifically, modelling and structural alignment were used to identify residues responsible for controlling discrimination of substrate targets with similar type of bonds.  lower than 1.0). Therefore, only four amino acids are required to produce both the most CE-and FAE-selective enzymes. This is a remarkable result which highlights the relative ease by which specificities of esterases can be interchanged, in a good agreement with previous studies (Levasseur et al., 2006). Since K 281 , D 282 , N 316 and K 317 are situated at specific loci within the catalytic core, their effect in promoting substrate promiscuity can be analysed in combination with the 3D model. K 281 I exerts a perturbation in the final access tunnel to the catalytic centre, generating a high (34-fold) or mild (1.4-fold) improvement in catalytic efficiency for FAE-and CE-like substrates respectively ( Fig. 4A and B). Saturation mutagenesis at the position 281 confirmed Ile as the best possible amino acid substitution for improving FAE phenotype. The Ile is a smaller residue making more accessible substrate channel that may correlate with higher values of kcat, whereas the substitution of a hydrophilic (Lys) by a hydrophobic (Ile) residue may correlate with enhanced binding of the substrates (lower Km values). This observation can also be extrapolated to the variant D 282 L (hydrophilic→hydrophobic); however, here, the slightly larger size of Asp than that of Leu appears to have more impact towards smaller substrates such as pNPC2 ( Fig. 4C and D). Lys 317 His appear to have a significant negative impact on the properties of the enzyme since the mutation was associated with a complete loss of FAE activity. The effect on pNP esters was mostly associated with a twofold decrease in substrate affinity while slightly affecting the reaction rates (kcat). Since it was suggested that the presence of at least one m/pmethoxy group in the cinnamate-like substrates is required for binding (Tarbouriech et al., 2005), one could argue that the mutation at His 317 may cause a partial unwinding of the loop and shortening the distance between binding residues and cinnamate-like substrates while maintaining the activity with common substrates (Fig. 4E and F). The mutation Asn 316 STOP produced a deletion of the C-terminal tail (25 amino acids) ( Fig. 4G and H). This fragment constitutes an a-helix located back to the catalytic cavity. Although, the deletion does not seem to be directly relevant to the catalytic core, we suggest that it may confer a higher flexibility, making catalysis more efficient, in particular, with shorter substrates. To further assess whether the higher flexibility is produced at expenses of lower structural stability, as suggested by the model, we studied the biophysical parameters of this protein variant after pre-incubation with different concentrations of guanidinium chloride (GdmCl) and at different temperatures. We showed that this variant tends to misfold to a large degree after pre-incubation with Ն 0.64 M GdmCl whereas the wildtype protein was mostly chemically stable (misfolding occurred Ն 2.1 M GdmCl) (Fig. S7).
Our results also show that the hybrid 3A6 protein can also be converted into a common CE enzyme through the deletion of 34-AA-long loop. The functional cost of the Gly 178 -Gly 211 insert was offset by a large change in the substrate channel that enables higher efficiency of exclusion of the shorter substrates (Fig. 4I). This configuration apparently provides an effective solution to re-channelling longer substrates to the active site. It is also likely to hamper the diffusion of shorter-chain pNP esters to and from the active site, explaining the observed difference in hydrolytic rates and binding efficiency and the inactivity with the cinnamate substrates. We assume that this additional insertion may have been acquired as a rudimentary cap in providing considerable structural and functional substrate variability, which may also corroborate their role in the phylogenetic separation of the 3A6-like proteins (Fig. S1). The present study suggests that mutations at specific amino acids were most effective for promoting substrate specificity compared with the random in vitro evolution. Whatever the case, independently of the mutagenesis method applied, the CE phenotype was present in all variants, in contrast to the gain or loss of FAE activity. This view is supported by the fact that after point mutations the activity could be fully abolished towards cinnamates, but not towards pNPC 2. Moreover, substrate discrimination appears to be associated to energy barriers. Free energy contribution analysis revealed that the energy barriers to affect CE-like activity were mutationindependent, whereas it is energetically more favourable (negative DDG) to create variants with improved FAE phenotype (Fig. 3C). In this context, it is possible to assume that the 3A6 enzyme may have a CE-like origin and that the newly evolved enzyme accumulated mutations in the vicinity of the original active site leading to an increase of its activity towards cinnamates (energetically favourable - Fig. 3C) without compromising its original esterase activity (energetically independent - Fig. 3C). We hypothesize that this scenario might be perhaps a consequence of the environmental conditions in vivo: earthworm gut microorganisms require an efficient enzymatic machinery to cope with a continuous flux of complex polymeric materials from the soil and from plant litter passing through the intestinal tract, thus demanding an instant response to adapt their enzymatic machinery to the changing environmental conditions. In this context, a substrate selection pressure applied may allow the creation/evolution of enzyme variants with broader substrate range but not necessarily high specificity (Levasseur et al., 2006). Our saturation mutagenesis experiments on 3A6 may have thus mimicked a possible scenario of natural evolution towards creating cinammate-selective enzymes from a 'common' carboxyl esterase and directly evidenced the close evolutionary relationship between both enzymes. Finally, the results obtained here are intriguing from both an academic and biotechnological point of view and also highlights the importance of uncultured microbial resources (Vieites et al., 2009) to access novel functionalities.

Experimental procedures
A full description of experimental procedures is available in Supporting information.

Protein samples
All hydrolases used in the present study (wild-type and variants) were cloned into pET-41 Ek/LIC vector, expressed with an N-terminal fusion to 6xHis tag and purified as described in Supporting information.

Enzyme characterization
Unlike indicated otherwise, hydrolytic activities were routinely measured and kinetic parameters determined as described by López-Cortés and colleagues (2007). Reaction conditions were: [E]o = 0-12 nM, [substrate] ranging from 0 to 50 mM, 100 mM Tris-HCl, pH 8.5, T = 40°C. Molecular mass of the 3A6 protein of 37.502 Da was considered for kinetic parameter calculations. Kinetic parameters were calculated from the Hanes-Woolf plot. Results shown are the average of three independent assays Ϯ the standard deviation.

Library construction and FAE screen
Earthworms were collected at the surface level (0-20 cm) from the Ecological Soil Station of the Lomonosov Moscow State University (Solnechnogorskiy District, Moscow Region, Russia). The worms were maintained in a terrarium at 12-15°C during 3 months and were fed with sterilized birch litter. The gut contents were collected after transferring worms onto a wet filter paper and keeping them at 4°C for 1 h. Approximately 0.5 g of the gut content was inoculated into the 1 l Erlenmeyer flask containing 500 ml of Getchinson medium (K2HPO3, 1.3 g; MgSO4¥7H2O, 0.3 g; CaCl2¥6H2O, 0.1 g; FeCl2¥6H2O, 0.01 g; NaNO3, 2.5 g; water 1 l; pH 7.2-7.4). The only carbon and energy source was the cellulose supplemented as a disk of filter paper of approximately 20 cm in diameter (Whatman 3MM), which was submerged in the medium. The enrichment was performed without shaking at 12°C until the filter paper was totally degraded (10 days). The DNA was extracted from 50 ml of this enrichment using G'NOME DNA Extraction Kit (Qbiogene). A fosmid library using CopyControl Genomic Library Production Kit was established in pCCFOS vector and E. coli EPI300-T1 R according to the protocols of the supplier (Epicentre). Fosmid clones, in total about 30 000, were picked with Qpix2 (Genetix, UK) and deposited in 384-microtiter plates containing Luria-Bertani (LB) medium with chloramphenicol (12.5 mg ml -1 ) and 15% (v/v) glycerol. To screen for FAE activity the clones were replicated onto large (22.5 ¥ 22.5 cm) square-shaped LB agar plates containing chloramphenicol (12.5 mg ml -1 ) (Qtray, Genetix, UK). In total, each plate contained 2304 clones. After overnight incubation at 37°C, the plates were overlaid with 20 ml of 5 mM EPPS buffer containing 0.4% (w/v) agarose, 0.456 mM phenol red and 320 ml of MF (120 mg ml -1 in dimethylformamide). A hydrolase-positive colony exhibiting a strong yellow halo after 2 min was picked, and the insert containing the hydrolase was sequenced after subcloning in pUC19 vector.

Rationale for mutagenesis
To understand structural backgrounds of substrate specificity we used four different mutagenesis approaches.
For site-directed mutagenesis, mutations were introduced into the gene 3A6 using the QuikChange mutagenesis kit (Stratagene) and the plasmid pUC19-3A6 as template with appropriate oligonucleotide pairs (s. details in Table S2). The resulting plasmids variants were then transferred into E. coli XL10 Gold and selected on the LB agar supplemented with 100 mg ml -1 ampicillin (Amp).
Fragment deletion in the 3A6 gene was carried out using a reverse PCR deletion, with pGEM-3A6 as template, PfuTurbo DNA polymerase and the pair of primers DGly 178 -Gly 211 Fwd and DGly 178 -Gly 211 Rev (Table S2). Conditions were as follows: 95°C -120 s, 30¥[95°C -45 s, 50°C -60 s, 72°C -120 s], 72°C -500 s. The PCR product was purified, diluted to a concentration of 5 ng ml -1 , further self-ligated with T4 DNA ligase in a total volume of 100 ml at room temperature, and transformed into E. coli TOP10. The resulting transformants were plated onto fresh LB agar plate containing 100 mg ml -1 Amp.
Saturation mutagenesis libraries were constructed using the QuikChange site-directed mutagenesis protocol (Stratagene), the pGEM-3A6 as emplate and appropriate pair of degenerated primers listed in Table S2. For each position subjected to the mutagenesis, the mutated plasmids were transformed in E. coli TOP10 (Invitrogen) and the resulting transformants plated onto fresh LB agar plate containing 100 mg ml -1 Amp.
Error-prone PCR mutagenesis (epPCR) was carried out using Taq polymerase (Sigma Chemical Co.) and pUC19-3A6 as template. The reaction was performed in 50 ml volume and contained 3% dimethyl sulfoxide (DMSO), 5 mM MnCl2, 1.5 mM MgCl2, 0.3 mM dNTPs, 2.5 U Taq polymerase, 5 ng of template and 4.5 pmol of oligonucleotides pet41Fwd and pet41Rev (Table S2). The concentration of template and MnCl2 was adjusted to achieve a mutation rate from 1 to 3 mutations per kb. The amplification program was as follows: 2 min at 95°C, 27 s at 94°C, 27 s at 53°C, followed by 28 cycles of 3 min at 74°C, and 10 min at 74°C. The amplified PCR products were purified from a 0.75% agarose gel using QIAEX II gel extraction kit from Qiagen, cloned into the plasmid pGEM (Promega) and transformed into E. coli TOP10 (Invitrogen) as recommended by the supplier, and the resulting transformants were plated onto LB agar plate containing 100 mg ml -1 Amp.
Clones from saturation mutagenesis and epPCR libraries were picked and growth in 96-well microtiter plates containing LB-100 mg ml -1 Amp. Plates were covered with gaspermeable films (Abgene) and incubated at 37°C with shaking at 300 r.p.m. for 16-20 h. The cultures were lysed in the plate wells using 50 ml of FastBreak reagent (Novagen) and 0.1 unit ml -1 of DN aseI, at room temperature for 30 min with shaking at 300 r.p.m. These lysates were used to determine CE-and FAE enzymatic activities in a BioTek Synergy HT spectrophotometer, with pNPC2 and MF as substrates respectively (López-Cortés et al., 2007).

Free energy calculations
The quantitative effects of mutations on kinetics are expressed as changes in free energy compared with that of the wild-type enzymes and were calculated as described elsewhere (Numata and Kimura, 2001;Bethel et al., 2006;Williams et al., 2006). Differences in the free energy change (DDG) caused by mutations were calculated using the equations: where R is the gas constant (8.31 J kmol -1 ), T is the absolute temperature (317 K) and DDGA and DDGB are the differential free activation and binding energies respectively. In order to analyse the mutual interplay of amino acids involved in pNPC2 and MF binding and cleavage, the free energy barriers to substrate specificity introduced by single, double and triple mutations were analysed. When the sum of the DDG values for single mutants is equal to that of the double/triple mutant [DDG(X + Y) = DDG(X) + DDG(Y)] the sites function independently.

Fig. S1
. Phylogenetic affiliation of 3A6 protein with related putative and experimentally analysed hydrolytic enzymes. Fig. S2. Temperature (A) and pH (B) optima of the wild-type 3A6 protein. Fig. S3. Overall three-dimensional structure of the 3A6 protein, as obtained by homology modelling. Fig. S4. Alignment of 3A6 with protein PDB 2C7B. Fig. S5. Multiple alignment of the C-terminal part of 3A6 esterase and related enzymes. Fig. S6. Far-UV CD spectra of 3A6 wild-type protein and 3A6-3I variant. Fig. S7. Stability effect by a single mutation Asn 316 STOP. Table S1. Annotation table of proteins homologous to the 3A6 protein sequence (BLAST search). Table S2. List of primers used in this study.
Please note: Wiley-Blackwell are not responsible for the content or functionality of any supporting materials supplied by the authors. Any queries (other than missing material) should be directed to the corresponding author for the article.