Get 20M+ Full-Text Papers For Less Than $1.50/day. Start a 14-Day Trial for You or Your Team.

Learn More →

Rational engineering of a native hyperthermostable lactonase into a broad spectrum phosphotriesterase

Rational engineering of a native hyperthermostable lactonase into a broad spectrum... www.nature.com/scientificreports OPEN Rational engineering of a native hyperthermostable lactonase into a broad spectrum phosphotriesterase Received: 29 August 2017 1 1,4 2 1,3 1 Pauline Jacquet , Julien Hiblot , David Daudé , Céline Bergonzi , Guillaume Gotthard , 1 1 3 Accepted: 16 November 2017 Nicholas Armstrong , Eric Chabrière & Mikael Elias Published: xx xx xxxx The redesign of enzyme active sites to alter their function or specificity is a difficult yet appealing challenge. Here we used a structure-based design approach to engineer the lactonase SsoPox from Sulfolobus solfataricus into a phosphotriesterase. The five best variants were characterized and their structure was solved. The most active variant, αsD6 (V27A-Y97W-L228M-W263M) demonstrates a large increase in catalytic efficiencies over the wild-type enzyme, with increases of 2,210-fold, 163-fold, 58- fold, 16-fold against methyl-parathion, malathion, ethyl-paraoxon, and methyl-paraoxon, respectively. Interestingly, the best mutants are also capable of degrading fensulfothion, which is reported to be an inhibitor for the wild-type enzyme, as well as others that are not substrates of the starting template or previously reported W263 mutants. The broad specificity of these engineered variants makes them promising candidates for the bioremediation of organophosphorus compounds. Analysis of their structures reveals that the increase in activity mainly occurs through the destabilization of the active site loop involved in substrate binding, and it has been observed that the level of disorder correlates with the width of the enzyme specificity spectrum. This finding supports the idea that active site conformational flexibility is essential to the acquisition of broader substrate specificity. Enzymes are extraordinary molecules that can accelerate chemical reactions, including those that spontaneously occur extremely slowly, and can make them happen many times per second . Some of these reactions are of high biotechnological or industrial interest and enzymes are, therefore, the focus of numerous studies . However, enzymes are often not directly compatible with industrial processes for various reasons ranging from their 3,4 5 sub-optimality or lack of stability under certain conditions . The properties of enzymes can be tuned by molec- 5–10 ular engineering to fit these new requirements using high throughput methods such as directed evolution and 11–13 /or low throughput methods, reducing the size of the library by focusing on positions that can be saturated , or 14,15 16 17 by using phylogenetically-inferred or structure-based mutations to create combinatorial , reduced libraries . Here, we focus on the molecular engineering of an enzymatic bio-decontaminant of organophosphates com- pounds (OPs), including nerve agents and insecticides. These compounds are neurotoxic because of their ability to irreversibly inhibit the acetylcholine esterase, a key enzyme in the nervous system . Massively used in agri- culture as pesticides , they were modified before World War II to increase their toxicity and develop chemical warfare agents . OPs still represent major environmental, health and security concerns, as they cause significant 21,22 19,20 pollution , and intoxication, and have been used in terrorist attacks . However, since current methods for removing them are cost prohibitive and lead to ecological concerns, enzymatic remediation is a promising way of 23–25 developing efficient bio-decontamination methods . Phosphotriesterase-like Lactonases (PLLs) are appealing candidates, being natural lactonases endowed with promiscuous phosphotriesterase activity . They have long been related to the bacterial phosphotriesterases (PTEs), enzymes which are capable of hydrolyzing insecticides such as parathion with high catalytic efficiency . 28,29 PTEs and PLLs belong to the amidohydrolase superfamily , and share the same (β/α) barrel topology. Some CNRS UMR 7278, IRD198, INSERM U1095, APHM, Institut Hospitalier Universitaire Méditerranée-Infection, Aix- Marseille Université, 19-21 Bd Jean Moulin, 13005, Marseille, France. Gene&GreenTK, IHU Méditerranée Infection, 19-21 Bd Jean Moulin, 13005, Marseille, France. University of Minnesota, Department of Biochemistry, Molecular Biology and Biophysics & Biotechnology Institute, St. Paul, MN, 55108, USA. Present address: MPI for Medical Research, Chemical Biology department (EPFL), Heidelberg, Germany. Pauline Jacquet, Julien Hiblot and David Daudé contributed equally to this work. Correspondence and requests for materials should be addressed to E.C. (email: eric.chabriere@univ-amu.fr) or M.E. (email: mhelias@umn.edu) Scientific REPO R TS | 7: 16745 | DOI:10.1038/s41598-017-16841-0 1 www.nature.com/scientificreports/ Data collections asA6 open asA6 closed asB5 asC6 asD6 asA1 PDB ID 5VRK 5VRI 5W3U 5W3Z 5W3W 5VSA Wavelength (Å) 0.97934 0.93340 0.97625 1.00882 1.00882 0.99987 Resolution (Å) 1.4 2.15 2.5 2.55 2.95 2.0 Space group P2 P2 2 2 P2 2 2 P2 2 2 P2 2 2 P2 2 2 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 Unit cell dimensions a (Å) 49.28 87.17 82.93 81.65 82.34 86.6 b (Å) 137.20 103.61 105.95 105.80 105.01 103.100 c (Å) 49.36 152.61 152.83 154.53 153.15 151.8 α (°) 90.0 β (°) 98.7 γ (°) 90.0 No. observed reflections 419571 (75591) 408826 (24712) 206443 (23317) 194083 (10849) 106247 (5379) 703275 (96197) No. unique reflections 121915 (22376) 74124 (4744) 47092 (5182) 43204 (2416) 28163 (1372) 91103 (12398) Completeness (%) 96.1 (94.3) 97.8 (95.5) 99.5 (99.8) 97.3 (98.0) 98.3 (99.9) 98.6% (99.7%) Rmerge (%) 8.2 (33.4) 9.8 (44.6) 6.1 (48.7) 5.55 (42.8) 6.6 (46.0) 9.0% (48.8%) Rmeasure (%) 9.7 (39.5) 10.8 (49.3) 6.9 (55.4) 6.2 (48.5) 7.7 (52.9) 9.6% (52.6%) I/σ(I) 9.07 (3.39) 14.31 (3.47) 12.84 (2.50) 15.68 (2.67) 12.89 (2.75) 15.59 (4.23) Last resolution shell 1.5–1.4 2.2–2.15 2.6–2.5 2.6–2.55 3.0–2.95 2.1–2.0 Redundancy 3.44 (3.38) 5.51 (5.21) 4.38 (4.50) 4.49 (4.49) 3.77 (3.92) 7.72 (7.76) Refinement statistics Resolution range (Å) 48.7128–1.40 49.0556–2.1500 49.6449–2.500 46.313–2.55 49.667–2.95 48.812–2.0 No. Reflections 115818 70416 44736 41042 26754 86547 Rfree /Rwork 16.84 / 13.36 21.35/17.23 27.17/21.22 24.15/19.97 28.31/24.68 20.55/16.32 No. protein atoms 5368 10122 9310 9657 9683 9987 No. water molecules 634 659 76 154 23 733 Average B factor (Å2) 19.095 24.984 86.6 80.375 95.992 29.885 rmsd from ideal Bond lengths (Å) 0.0088 0.0028 0.0050 0.0110 0.0070 0.0211 Bond angles (°) 1.0703 0.7309 0.8490 1.1730 0.9630 0.6950 Table 1. Data collection and refinement statistics of SsoPox variants structures. 30–38 PLL representatives are extremely thermostable and are particularly promising candidates for biotechno- 39,40 41 logical use by virtue of their high stability and resistance towards denaturing agents . PLLs, however, exhibit much lower phosphotriesterase activity than PTEs. Engineering experiments with the aim of increasing PLLs’ phosphotriesterase activity were performed. In particular, because PTEs and PLLs’ active sites mainly differ in 16,26,42 their loop length (7 and 8), sequences and conformations , loop grafting experiments were performed and 16,34 proved to be difficult. This approach resulted in aggregation-prone or inactive enzymes , but was successful when using a stepwise approach . Recently, more classical engineering techniques were successful in increasing the phosphotriesterase activity of the PLL DrOPH . 32,45 In this study, we worked on the PLL SsoPox from the hyperthermophilic archaea Sulfolobus solfataricus . SsoPox exhibits high acyl-homoserine lactones and oxo-lactones hydrolysis abilities , and promiscuous, low phosphotriesterase activity . Previous work showed that engineering can successfully improve SsoPox catalytic 11,47 efficiency against phosphotriesters such as ethyl-paraoxon . These previous studies highlighted the role of a key active site loop 8 residue, W263, which modulates the active site loop conformational space . In this study, we aimed to further improve the phosphotriesterase activity of SsoPox, including against phosphothionoesters. SsoPox is a very appealing candidate for the bioremediation of phosphotriesters because of its unique thermal stability, and its ability to resist aging, solvent and protease treatments . On the other hand, SsoPox is a poor phosphothionoesterase, exhibiting a strong preference (100-fold) for methyl-paraoxon over methyl-parathion, a feature referred to as the thiono-effect . We used a rational strategy: taking advantage of the structural similarity between SsoPox and the bacterial PTE from the mesophilic Brevundimonas diminuta (BdPTE) , we designed and produced a structure-based combinatorial library, and screened this library for improved phosphotriesterase activity with the aim of obtaining heat stable, highly active variants with broad specificity. We obtained sev - eral improved variants, including against methyl-parathion, with a 2,210-fold improvement in activity. All these variants but variant αsB5 contained a mutation at position 263. Interestingly, the addition of the W263 muta- tion using site-saturation mutagenesis on an αsB5 background revealed an incompatibility of these mutations, instead of the expected increase in phosphotrieseterase activity. Of these improved variants, we obtained a mutant 4 −1 −1 that lost the thiono-effect against methyl-parathion, exhibiting a > 2,000-fold (k /K ~10 M .s ) as referred cat M to the wild-type enzyme, and a methyl-paraoxon/methyl-parathion activity ratio of ~1. The X-ray structures of several improved variants reveal that the selected mutations increase the active site loop conformational flexibility and reshape the active site. The best variants were extensively characterized against eight commercially available Scientific REPO R TS | 7: 16745 | DOI:10.1038/s41598-017-16841-0 2 www.nature.com/scientificreports/ Figure 1. Structure-based design approach. (A) Pipeline used for the selection of improved mutants. (B) Structural alignments of SsoPox (red) and BdPTE structures (yellow) guided the mutational design. (C) Active site superposition of BdPTE (green) and SsoPox (cyan) active sites. Residues with structural equivalent (orange sticks) in both structures were used to design mutations of SsoPox’s residues into the BdPTE corresponding residues. Positions with no structural equivalent (red sticks) were designed using modelling tools (purple sticks) to mimic the active site cavity size, shape and chemical properties of the BdPTE crystal structure model (Figure S1). (D) The mutations data base consists of mutations to the BdPTE sequence when there is a structural equivalence (orange), and of mutations designed when there is no structural equivalence (red). The third set of mutations (black) adds more diversity at two selected positions, 258 and 263. (E) Once the library had been validated, primers carrying mutations were used to shuffle them and generate a gene library with random combinations of selected mutations. (F) The library was screened for paraoxonase activity to identify enzymes with improved proficiency. insecticides, along with previously reported SsoPox monovariants harboring mutation at residue 263 and their specificity spectrum was determined. Results and Discussion e l Th actonase SsoPox was engineered for higher phosphotriesterase activity using structure-based combinatorial libraries. By comparing structures of enzymes with similar topology, it has been possible to redesign, using mod- elling tools, the active site cavity of SsoPox to mimic as closely as possible that of BdPTE (Fig. 1; Table S1). This method consisted of two main steps: (i) the identification of mutations corresponding to equivalent positions at the enzyme active sites using structural alignment and (ii) the rational selection of mutations, in a position with no equivalent residue, which mimicked the shape and chemical nature of the target enzyme cavity. A mutation dataset was thus obtained and used to develop a combinatorial library consisting of random combinations of our pre-selected mutations. Scientific REPO R TS | 7: 16745 | DOI:10.1038/s41598-017-16841-0 3 www.nature.com/scientificreports/ Figure 2. Chemical structure of substrates used in this study. The chemical structure of ethyl-paraoxon (I), ethyl-parathion (II), methyl-paraoxon (III), methyl-parathion (IV), malathion (V), chlorpyrifos (VI), diazinon (VII), fenitrothion (VIII), fensulfothion (IX), coumaphos (X), 3-oxo-C10 AHL (l) (XI), undecanoic-γ-lactone (XII) and undecanoic-δ-lactone (XIII) generated using ChemDraw software. Structure-based design produced mutants with improved phosphotriesterase activity. The structure-based design considerably reshaped the active site cavity of SsoPox, theoretically bringing it closer to that of BdPTE, both in shape and nature (Figure S1). A total of 14 key positions within the active site were iden- tified, and were mutated to specific residues, or degenerated into several possible amino acids (see Methods). The paraoxonase activity of 184 randomly selected clones was screened and the 14 best were sequenced. Most of these 14 were found to possess the same mutations (e.g. W263F/L/M/A; 73% of cases) and (Figure S2A; Table S3). Of them, the most active variants, αsA1 (C258L-I261F-W263A), αsA6 (F46L-C258A-W263M-I280T), αsB5 (V27A-I76T-Y97W-Y99F-L130P-L226V), αsC6 (L72I-Y99F-I122L-L228M-F229S-W263L), and αsD6 (V27A-Y97W-L228M-W263M) were subject to kinetic and structural characterization. It should be noted that αsA1 was obtained in a previous engineering effort , and that four of the five selected variants contained a sub- stitution of the key residue W263. We previously highlighted the role of this residue in substrate binding and modulation of the active site loop 8 conformation, resulting in an increase in activity against promiscuous sub- strates . Variant αsB5 does not harbor any substitution at position W263. We therefore decided to test the W263 substitution on an αsB5 background, using saturation mutagenesis strategy. The five variants (α sA1, αsA6, αsB5, αsC6 & αsD6) were characterized against several phosphotriesters, including ethyl-paraoxon (I), ethyl-parathion (II), methyl-paraoxon (III), methyl-parathion (IV) and mala- thion (V) (Fig. 2). Overall, all variant exhibits improved catalytic efficiencies against all tested substrates, with the exception of αsA1 with malathion (~2-fold decrease) (Fig. 3A; Table 2). e Th best improved variant, αsD6, exhib- its a 2,210-fold increase in catalytic efficiency against methyl-parathion (Table  2). The largest improvements were observed for ethyl/methyl-parathion, two bad substrates for wild-type SsoPox. Interestingly, while the wt enzyme shows a clear preference for small substituents , most selected variants lost this preference, possibly indicating an enlargement of the active site cavity. These improvements were in the range of what was previously obtained but with higher intensive mutation protocols and for only one OP substrate . Finally, it was noted that two variants (αsA6 & αsC6) presented a substrate inhibition for some thiono-phosphotriesters. The five selected variants were further evaluated for their lactonase activity (Table  2). Interestingly, variants αsA6 and αsC6 exhibit enhanced lactonase activity against γ- and δ-lactones, while αsA1, αsB5 and αsD6 show reduced lactonase catalytic efficiencies. This emphasizes the fact that the improvement of the phosphotriesterase activity does not necessarily compromise the cognate, lactonase activity of the enzyme. These mutations there- fore dramatically increased a new activity (phosphotriesterase) without the complete loss of the native/original function (lactonase) . These mutants are die ff rent from those generated on another PLLs, through loop inser - tion. While the stepwise insertion of residues in loop 7 can significantly increase the ability of PLLs to degrade phosphotriesters, it also drastically decreases their lactonase activity . Insertion in loop 7 has been previously described as a key evolutionary event in the transition from lactonase to phosphotriesterase . Scientific REPO R TS | 7: 16745 | DOI:10.1038/s41598-017-16841-0 4 www.nature.com/scientificreports/ Figure 3. Activity characterization of the improved variants and bioremediation potential. (A) Phosphotriesterase catalytic efficiency comparison between wt-SsoPox and the best selected variants. Data for wt-SsoPox are from Hiblot et al. . (B) Determination of the reaction time for hydrolyzing 95% using the αsD6 −2 −3 variant of a solution of pesticide at 250 µM. Two enzyme-to-substrate ratios ( [E]/[S]), 10 and 10 , were considered. Structural analysis of improved variants reveals increased mobility of the active site loop 8. e cr Th ystal structures of variants α sA1, αsA6, αsB5, αsC6 and αsD6 were solved (Table 1). While most structures were solved at high or medium resolutions (1.4–2.55 Å), mutant asD6 could only be solved at low resolution (2.95 Å). Given this resolution, we limited our interpretations and retained this data for the consistency of the study as it exhibits structural features (loop 8 disorder) that are consistent with all the other mutants. Overall, the mutants’ structures are similar to that of wt-SsoPox. However, the active site cavities are larger for the mutants: most selected mutations replace residues with smaller ones: W263 is mutated into M/L/I, V27 into A, I76 into T, L130 into P, Y99 into F. Only two selected substitutions relate to bulkier side chains: Y97W and L228M. This enlargement is difficult to quantify even with the structure of these mutants because of the extremely mobile nature of loop 8 in the variants. Enlarging the active site was a main objective of the active site redesign of SsoPox, to enable the bulkier phosphotriesters to bind within the natural lactonase enzyme. Additionally, the loop 8 conformations differ in mutants, as compared to the wt -enzyme. Because part of loop 8, including position 263, is located at the enzyme dimer interface, altered loop 8 conformations modulate the relative orientation of both monomers. In the case of the mutants characterized in this study, the dimer reorientation yields significant dis- placement up to 5.2 (e Å .g., α sA1, between equivalent carbon α positions), as compared to wt-SsoPox (Figure S5). 11,42 Similar reorientations were previously observed upon substrate binding or mutations of W263 . Although the active sites of the improved variant structures superpose well onto the wt-structure, all selected mutants display altered loop 8 conformations. Changes are subtle for some of the variants (e.g. αsA1), but are large for others (e.g. αsC6) (Fig. 4). Notably, mutant αsA6 was crystallized in two different conformations: one in which loop 8 adopts a wt-like conformation, referred to as closed conformation (CC), and another one where loop 8 is unfolded, referred to as open conformation (OC) (Figure S6). In the two other mutants, αsD6 and αsB5, loop 8 could not be modeled due to the lack of electronic density, which was likely to be due to the high level of motion of this enzyme region. Analysis of the normalized thermal motion B-factor supports this hypothesis. It also confirms that loop 8 is highly mobile in all selected mutants, as compared to wt-SsoPox (Fig. 5). The higher mobility of loop 8 may also partly explain the improved ability of the variants to hydrolyze phosphotriesters with large substituents. Improved mutant αsD6 lost the thiono-effect. SsoPox presents a marked preference for oxono-OP substrates as compared to thiono-ones (>100-fold in catalytic efficiency) . Conversely, PTEs do not exhibit such 27,51 a drastic preference, paraoxon being a slightly better substrate than parathion for the enzyme . Interestingly, Scientific REPO R TS | 7: 16745 | DOI:10.1038/s41598-017-16841-0 5 www.nature.com/scientificreports/ Enhancement/ −1 −1 −1 Substrate SsoPox k (s ) K (µM) K (µM) k /K (M .s ) wt cat M I cat M 1 4 2 wt° 1.26 ± 0.13 × 10 2.43 ± 0.37 × 10 — 5.12 ± 1.31 × 10 1 1 2 4 αsA1 2.70 ± 0.23 × 10 8.00 ± 1.55 × 10 — 3.37 ± 0.71 × 10 65 3 3 αsA6 5.87 ± 0.47 1.63 ± 0.63 × 10 — 3.61 ± 1.43 × 10 7 1 2 4 αsB5 2.06 ± 0.14 × 10 6.01 ± 0.96 × 10 — 3.43 ± 1.46 × 10 67 1 3 4 αsC6 4.27 ± 0.25 × 10 1.50 ± 0.16 × 10 — 2.86 ± 0.35 × 10 55 1 3 4 Ethyl-paraoxon (I) αsD6 3.22 ± 0.18 × 10 1.09 ± 0.12 × 10 — 2.95 ± 1.55 × 10 58 αsB5 W263M * * — 5.62 ± 0.11 × 10 11 αsB5 W263L * * — 8.07 ± 0.28 × 10 16 1 3 3 αsB5 W263I 4.09 ± 0.89 × 10 6.64 ± 1.75 × 10 — 6.16 ± 5.07 × 10 12 + 3 W263L * * — 2.37 ± 0.33 × 10 5 + 2 3 W263M 6.82 ± 0.57 9.31 ± 1.63 × 10 — 7.33 ± 1.42 × 10 14 wt° ND ND ND ND ND αsA1 ND ND ND ND ND −2 2 2 αsA6 3.20 ± 0.20 × 10 1.34 ± 0.18 × 10 2.39 ± 0.35 × 10 ND 250 ± 1020 −1 2 2 αsB5 2.69 ± 0.15 × 10 2.72 ± 0.49 × 10 — 9.86 ± 3.14 × 10 ND −3 1 2 αsC6 8.53 ± 0.30 × 10 7.70 ± 1.10 × 10 — 1.10 ± 0.16 × 10 ND Ethyl-parathion (II) −1 2 3 αsD6 7.35 ± 0.36 × 10 5.34 ± 0.62 × 10 — 1.38 ± 0.57 × 10 ND −2 2 2 αsB5 W263M 5.24 ± 0.23 × 10 2.02 ± 0.31 × 10 — 2.59 ± 0.74 × 10 ND −2 2 2 αsB5 W263L 4.35 ± 0.16 × 10 3.46 ± 0.37 × 10 — 1.26 ± 0.44 × 10 ND −2 2 2 αsB5 W263I 4.74 ± 0.21 × 10 2.61 ± 0.37 × 10 — 1.82 ± 0.57 × 10 ND W263L ND ND — ND ND W263M ND ND — ND ND 3 3 wt° 2.71 ± 0.64 2.14 ± 0.68 × 10 — 1.27 ± 0.70 × 10 1 1 3 4 αsA1 4.35 ± 0.87 × 10 1.90 ± 0.52 × 10 — 2.29 ± 0.78 × 10 18 2 4 αsA6 5.89 ± 0.55 5.46 ± 1.02 × 10 — 1.08 ± 0.23 × 10 8 Methyl-paraoxon (III) αsB5 * * — 4.31 ± 0.14 × 10 3 1 2 4 αsC6 2.36 ± 0.37 × 10 7.59 ± 2.04 × 10 — 3.11 ± 0.97 × 10 25 1 3 4 αsD6 4.25 ± 0.52 × 10 2.08 ± 0.34 × 10 — 2.04 ± 0.42 × 10 16 −3 2 wt° 1.10 ± 0.02 × 10 1.21 ± 0.10 × 10 — 9.09 ± 0.90 1 −2 2 1 αsA1 1.30 ± 0.50 × 10 3.53 ± 0.34 × 10 — 3.68 ± 1.46 × 10 4 −2 2 1 αsA6 1.55 ± 0.15 × 10 2.54 ± 0.39 × 10 1 520 ± 384 6.10 ± 1.11 × 10 7 Methyl-parathion (IV) −1 2 2 αsB5 1.50 ± 0.04 × 10 1.58 ± 0.15 × 10 — 9.49 ± 0.94 × 10 104 −3 2 1 αsC6 3.00 ± 0.10 × 10 1.21 ± 0.15 × 10 — 2.48 ± 0.32 × 10 3 2 4 αsD6 6.89 ± 0.35 3.43 ± 0.44 × 10 — 2.01 ± 0.28 × 10 2210 −4 2 wt° 8.90 ± 0.40 × 10 1.60 ± 0.29 × 10 — 5.56 ± 1.26 1 −4 2 αsA1 7.10 ± 0.40 × 10 2.22 ± 0.35 × 10 — 3.20 ± 0.54 0.5 αsA6 ND ND ND ND ND −2 2 1 αsB5 1.84 ± 0.16 × 10 9.13 ± 1.61 × 10 — 2.02 ± 1.01 × 10 4 −2 1 2 αsC6 2.47 ± 0.44 × 10 3.20 ± 1.90 × 10 820 ± 343 7.72 ± 4.78 × 10 139 −1 2 2 Malathion (V) αsD6 1.00 ± 0.02 × 10 1.13 ± 0.10 × 10 — 9.08 ± 2.00 × 10 163 −3 2 1 αsB5 W263M 8.92 ± 0.67 × 10 8.61 ± 1.32 × 10 — 1.04 ± 0.51 × 10 2 −2 2 1 αsB5 W263L 1.11 ± 0.08 × 10 9.28 ± 1.32 × 10 — 1.20 ± 0.60 × 10 2 −2 2 1 αsB5 W263I 1.18 ± 0.08 × 10 7.80 ± 1.08 × 10 — 1.51 ± 0.71 × 10 3 W263L ND ND — ND ND W263M ND ND — ND ND wt ND ND — ND ND αsD6 * * — 3.12 ± 0.09 × 10 ND αsB5 * * — 2.07 ± 0.05 × 10 ND αsB5 W263M * * — 1.64 ± 0.03 × 10 ND Chlorpyrifos (VI) αsB5 W263L * * — 1.46 ± 0.03 × 10 ND αsB5 W263I * * — 1.42 ± 0.03 × 10 ND W263L ND ND — ND ND W263M ND ND — ND ND Continued Scientific REPO R TS | 7: 16745 | DOI:10.1038/s41598-017-16841-0 6 www.nature.com/scientificreports/ Enhancement/ −1 −1 −1 Substrate SsoPox k (s ) K (µM) K (µM) k /K (M .s ) wt cat M I cat M wt ND ND — ND ND 3 3 αsD6 8.93 ± 0.79 2.30 ± 0.31 × 10 — 3.89 ± 2.56 × 10 ND 2 3 αsB5 1.26 ± 0.06 4.10 ± 0.65 × 10 — 3.03 ± 0.98 × 10 ND −1 2 2 αsB5 W263M 3.89 ± 0.38 × 10 3.91 ± 1.04 × 10 — 9.92 ± 3.65 × 10 ND Diazinon (VII) −1 2 3 αsB5 W263L 5.61 ± 0.23 × 10 4.74 ± 0.53 × 10 — 1.18 ± 0.44 × 10 ND −1 2 3 αsB5 W263I 4.01 ± 0.30 × 10 3.24 ± 0.72 × 10 — 1.24 ± 0.42 × 10 ND W263L ND ND — ND ND W263M ND ND — ND ND wt ND ND — ND ND −1 2 3 αsD6 7.94 ± 0.50 × 10 5.28 ± 0.80 × 10 — 1.50 ± 0.62 × 10 ND −2 2 1 αsB5 1.40 ± 0.20 × 10 3.14 ± 1.43 × 10 — 9.11 ± 1.40 × 10 ND −2 2 1 αsB5 W263M 2.22 ± 0.12 × 10 3.31 ± 0.53 × 10 — 6.69 ± 2.31 × 10 ND Fenitrothion (VIII) −2 2 1 αsB5 W263L 2.63 ± 0.13 × 10 6.54 ± 0.75 × 10 — 4.01 ± 1.79 × 10 ND −2 2 1 αsB5 W263I 2.01 ± 0.08 × 10 3.97 ± 0.46 × 10 — 5.05 ± 1.81 × 10 ND −3 2 W263L 3.66 ± 0.15 × 10 4.15 ± 0.44 × 10 — 8.82 ± 3.32 ND −4 2 W263M 9.49 ± 0.43 × 10 1.57 ± 0.28 × 10 — 6.06 ± 1.55 ND wt ND ND — ND ND −2 2 1 αsB5 4.28 ± 0.39 × 10 9.41 ± 1.72 × 10 — 4.55 ± 2.30 × 10 ND −1 3 1 αsD6 4.06 ± 1.60 × 10 8.60 ± 3.97 × 10 — 4.72 ± 4.04 × 10 ND αsB5 W263M * * — 1.12 ± 0.46 × 10 ND Fensulfothion (IX) −2 3 1 αsB5 W263L 4.94 ± 1.35 × 10 2.57 ± 1.04 × 10 — 1.92 ± 1.30 × 10 ND αsB5 W263I * * — 6.50 ± 0.07 ND W263L ND ND — ND ND W263M ND ND — ND ND wt ND ND — ND ND αsD6 * * — 1.64 ± 0.02 × 10 ND αsB5 * * — 7.82 ± 0.14 × 10 ND αsB5 W263M * * — 1.15 ± 0.02 × 10 ND Coumaphos (X) αsB5 W263L * * — 1.16 ± 0.05 × 10 ND αsB5 W263I * * — 1.11 ± 0.04 × 10 ND W263L * * — 6.01 ± 0.18 ND W263M * * — 4.80 ± 0.32 ND wt° * * — 3.16 × 10 1 αsA1 * * — 1.41 × 10 0.004 αsA6 * * — 1.78 × 10 0.06 3-oxo-C10-AHL (l) (XI) αsB5 ND ND — ND ND αsC6 * * — 5.84 × 10 0.02 αsD6 ND ND — ND ND wt° * * — 2.36 × 10 1 αsA1 * * — 1.63 × 10 0.7 αsA6 * * — 4.40 × 10 186 Undecanoic-γ-lactone (XII) αsB5 * * — 9.05 × 10 0.4 αsC6 * * — 5.73 × 10 243 αsD6 * * — 1.60 × 10 7 wt° * * — 7.86 × 10 1 αsA1 * * — 2.55 × 10 0.03 αsA6 * * — 9.65 × 10 123 Undecanoic-δ-lactone (XIII) αsB5 * * — 5.33 × 10 0.08 αsC6 * * — 3.60 × 10 5 αsD6 * * — 7.50 × 10 0.9 Table 2. Kinetic characterization of SsoPox variants. *Corresponds to undetermined values because enzyme saturation could not be reached. ND corresponds to undetectable values because substrates were not hydrolyzed 46 + 11 or too poorly hydrolyzed to be observable. °Data taken from Hiblot et al. . Data taken from Hiblot et al. . some selected variants retained a strong preference for oxono-phosphotriesters (~10-fold), whereas others such as αsD6 lost the thiono-effect (Figure S7). The wild-type enzyme exhibits a low K value with methyl-parathion −3 −1 (121 µ M), but also a very low k value (1.10 × 10 s ) (Table 2). This may indicate a stronger interaction cat Scientific REPO R TS | 7: 16745 | DOI:10.1038/s41598-017-16841-0 7 www.nature.com/scientificreports/ between the thiono moiety and the metal cations (low K value), which would subsequently impair catalysis (low k values) . Interestingly, the improved variants do not exhibit significant difference in K values for cat M thiono-phosphotriesters, but show a dramatic increase in k values, as compared to wild-type enzymes (e.g. cat −1 −3 −1 k value of 6.89 s and 1.10 × 10 s , as compared to methyl-parathion for the αsD6 variant and wild-type cat enzymes, respectively). The two selected variants, namely αsB5 and αsD6, which lost the thiono-effect have only two mutations in common: V27A and Y97W. Both residues are located relatively close to the bi-metallic active site, yet do not interact directly with the metals (V27-Fe: 6.7 Å, Y97-Co: 4.3 Å; Fig.  6). Because other PLLs have 33,53 been shown to exhibit high charge coupling between the β metal and the conserved tyrosine residue , Y97 is a likely candidate for the thiono-effect, and will be examined in future studies. Mutations of W263 and αsB5 are incompatible. αsB5 (V27A-I76T-Y97W-Y99F-L130P-L226V) is the only selected mutant that does not contain a substitution of position W263. Because W263 substitutions were previ- ously shown to be key for increasing SsoPox’s phosphotriesterase activity, we used saturation mutagenesis to introduce W263 mutations to the background of αsB5 . The best selected mutants against paraoxon (Figure S2C) were found to be αsB5-W263I/L/M. Surprisingly, none of these mutants demonstrated an increase in phosphotriesterase activity in comparison with αsB5, but rather a decrease in activity (Table 2). This is intriguing because, taken individually, these mutations have been shown to greatly improve phosphotriesterase activity. For example, with ethyl-paraoxon as a substrate, W263M and αsB5 produce a catalytic efficiency improvement of 14 and 67-fold, respectively, while the combination αsB5-W263M increases activity only 11-fold (Table 2). This example of negative epistasis could be due to the mode of action of these mutations. Structures reveal that αsB5 mutations have a destabilizing effect on loop 8, as illustrated by the B factor analysis. The same was observed for W263M in a previous study . Combining these two dest- abilizing sets of mutations may have harmed the active site integrity, as previously described with another enzyme as “conformational active site disorder” , and in particular the necessary alignment and pre-organization of the catalytic residues. The destabilizing effect of these mutations can also be observed on the overall enzyme stability, with α sB5, W263M, αsB5-W263M exhibiting T values of 70.4 °C, 85.3 °C and 69.1 °C, respectively, as compared to 106 °C for the wt enzyme (Table S4). Here, the effects of the combined 5 mutations harbored by variant α sB5 together with the previ- ously identified positive substitutions at residue W263 were not additive underlining that single beneficial mutations do not necessarily induce positive epistatic effects. Synergistic effects of mutations are oen ft complex to predict and most engineering strategies usually lead to an optimization plateau difficult to overcome . Improvement of engineered variants over W263 mutation. In a previous report, we noted that muta- tion of W263 into all the other 19 residues lead to an increase of the promiscuous paraoxonase activity . Here we note that most of the improved variants (i.e. αsA1, αsA6, αsB5, αsC6 & αsD6) (i) exhibit higher paraoxonase activity than W263L/M, and (ii) foremost are capable of hydrolyzing phosphotriesters that are not substrates for the wild-type enzyme or W263 mutants, such as ethyl-parathion, chlorpyrifos or fensulfothion (Table 2). The activity spectrum of these engineered variants against the tested phosphotriesters is therefore much wider than the wild-type enzyme and W263 mutants. Relation between active site loop disorder and broad enzymatic specificity. The organophospho - rus compound bioremediation potential of several mutants was investigated. In particular, the best variants, αsD6 and αsB5, alongside other improved variants αsB5-W263I/L/M, W263L/M and the wt-enzyme, were character- ized for their phosphotriesterase activity with ethyl-paraoxon (I), ethyl-parathion (II), malathion (V), chlorpy- rifos (VI), diazinon (VII), fenitrothion (VIII), fensulfothion (IX) and coumaphos (X) (Fig. 2). Interestingly, we found that most mutants were able to hydrolyze fensulfothion, previously reported as an inhibitor of the wild-type enzyme . Fensulfothion was indeed previously shown to bind head-to-tail into the wild-type enzyme active side , a non-productive binding mode. Changes in the active site loop 8 conformational ensemble of engineered mutants might have allowed for the productive binding of fensulfothion. All mutants selected and used in this study are destabilized as compared to the wt-enzyme (Fig. 7; Table S4). Structures reveal that these mutations, located on loop 8, contribute to increasing the mobility of the loop. Interestingly, the ability of the tested mutants to hydrolyze a wide range of phosphotriesters, including com- pounds that were inhibitors for the wt-enzyme, correlates with a lower T (Fig. 7). This correlation suggest that the active site loop 8 degree of disorder may modulate the enzyme’s ability to bind and hydrolyze a variety of phosphotriesters, a wider conformational ensemble of loop 8 being correlated with a broader enzyme specificity. This observation is consistent with previous studies highlighting the importance of the destabilization of active 54,56 site loops to evolve new activities and new substrate specificity . Such flexible active site loops involved in enzymatic specificity supports the notion of fold polarity: a portion of the active site (e.g. the loop) is weakly connected to the enzyme scao ff ld and thereby makes it possible for new functions to evolve with few changes . Mutant αsD6 is an active, broad spectrum phosphotriester biodecontaminant. Mutant αsD6 was further investigated for bioremediation considerations. Organophosphate-based pesticides are usually spread with concentrations in the millimolar range leading to micromolar-ranged contaminations. We thus investigated the potential of αsD6 for decontaminating pesticide solutions at 250 µM to simulate a contamination of ground- −2 −3 water or runoff waters. We investigated two [E]/[S] ratios, 10 and 10 respectively, for decontaminating OP solutions and determined the time required to hydrolyze 95% of pesticide preparations (Fig. 3B; Figure S4). When −3 considering the [E]/[S] ratio of 10 , 95% degradation was achieved for four substrates (ethyl-paraoxon, diazinon, fenitrothion and coumaphos) within 10 minutes. Two substrates were decontaminated within an hour (ethyl-parathion and chlorpyrifos) and fensulfothion was degraded in three hours. Scientific REPO R TS | 7: 16745 | DOI:10.1038/s41598-017-16841-0 8 www.nature.com/scientificreports/ Figure 4. Loop 8 conformations in improved mutants. Active site loop 8 region of mutant structures superimposed with the wt-SsoPox structure. The active site location is indicated by the presence of the bi- metallic catalytic center shown as spheres. The wt-SsoPox structure loop 8 conformation is shown as dark green cartoon. (A) αsA1′s loop 8 is shown in orange cartoon, (B) αsD6′s loop 8 is shown in green cartoon, and is not entirely visible in the electronic density maps, (C) αsC6′s loop 8 is shown in blue cartoon, (D) αsB5′s loop 8 is shown in yellow cartoon, and is not entirely visible in the electronic density maps, (E) αsA6-CC’s loop 8 is shown in magenta cartoon, (F) αsA6-OC’s loop 8 is shown in cyan cartoon. −2 When increasing enzyme concentration 10-fold ([E]/[S] ratio of 10 ), five substrates (ethyl-paraoxon, ethyl-parathion, diazinon, fenitrothion and coumaphos) were decontaminated within two minutes, while fensul- fothion needed ~20 minutes and the two other tested (chlorpyrifos and malathion) required 40 minutes. Given Scientific REPO R TS | 7: 16745 | DOI:10.1038/s41598-017-16841-0 9 www.nature.com/scientificreports/ Figure 5. Mutants’ active site loop 8 exhibit higher conformational flexibility. (A) Cartoon representation of the superposition of the structures (monomers) of wt-Ssopox (dark green), αsA1 (orange), αsB5 (yellow), αsD6 (green), αsC6 (blue), αsA6-CC (magenta) and αsA6-OC (cyan). (B) Smoothed ribbon representation of the superposition of the structures (monomers) of of wt-Ssopox (green), αsA1 (orange), αsB5 (yellow), αsD6 (green), αsC6 (blue), αsA6-CC (magenta) and αsA6-OC (cyan). (C) Positional distributions of normalized B-factor values (the x-axis represents residue number) for wt-SsoPox (black line), αsA1 (orange), αsB5 (yellow), αsD6 (green), αsC6 (blue), αsA6-OC (blue line) and (magenta line). A zoomed inset highlights the loop 7 and 8 sequence region. both its catalytic proficiency and high thermal stability (T = 82.5 °C), αsD6 is a promising candidate for organ- ophosphorus compound bioremediation and, in particular, the decontamination of water runoffs, soils, food products and materials. Moreover, SsoPox variants were previously shown to be compatible with immobilization 48,58 steps including alginate beads and polyurethane-based coatings , and could be of prime interest for the devel- opment of filtration devices for water treatment purposes. Methods Structure-based identification of mutations. Sharing the same (α/β) barrel topology, SsoPox and BdPTE structures (Pdb ID 2vc5 and 1dpm, respectively) were superimposed (Fig. 1) using PyMol , making it possible to identify structurally equivalent residues at their respective active sites. In an effort to reshape the SsoPox active site, side chains that were well superimposed between the two enzymes were mutated using Coot into the residue present in BdPTE , with the exception of V27 which was mutated in Ala (and not Gly) to avoid increased entropy (Fig. 1; Table S1). However, due to the major differences in loop 8 and 7 between PLLs and PTEs, numerous residues were not superimposable. These active site residues (Y99, L228, F229 and W263) were therefore mutated in an effort to mimic the Bd PTE active site cavity in terms of shape and chemical nature, as Scientific REPO R TS | 7: 16745 | DOI:10.1038/s41598-017-16841-0 10 www.nature.com/scientificreports/ Figure 6. Localization of mutations harbored by SsoPox monovariants, αsA1, αsA6, αsB5, αsC6 and αsD6. Mutations are highlighted in red spheres using the SsoPox wild-type (PDB 3UF9) as a template. Figure 7. Overview of SsoPox wild-type and variants properties. Catalytic efficiencies towards organophosphate-based pesticides are presented using a color code ranging from white to dark blue. Melting temperatures are presented using a color code ranging from red to green. illustrated by the chimeric reconstruction of SsoPox carrying all these mutations (Figure S1). Others mutations were also implemented in the dataset, e.g. C258A and C258L and different substitutions of the key position W263 (F/L/M/A) that have been shown to improve the phosphodiesterase activity . Synthesis of the combinatorial library. The SsoPox coding gene was amplified from the previously described pET22b-SsoPox plasmid . For each of the 14 mutations, a mutagenesis primer of 30–33 bp (~10–15 bp from each mutation side) was synthesized (Table S2). For close mutations, primers can share several mutations and different primers were synthesized and mixed to keep an equivalent statistical probability for each of the 14 mutations. Typically, 2 pmol of primer mixture and 100 ng of DNaseI (TaKaRa)-generated fragments of the SsoPox gene were assembled as previously described , followed by nested PCR with external cloning primers (SsoPox-lib-pET-5′/3′) (Table S2). The PCR product was then cloned into a customized version of the pET - 32b(-ΔTrx) plasmid and then electroporated into E. cloni cells (Lucigen, USA). Aer a ft gar-plate growth, plasmid extraction was performed to create the plasmid bank. Variant genes were PCR amplified using T7-prom and Scientific REPO R TS | 7: 16745 | DOI:10.1038/s41598-017-16841-0 11 www.nature.com/scientificreports/ pET-RP primers (Table S2) and sequenced. The combinatorial library shows an average of 5.3 ± 3.3 library muta- tions by sequencing 10 randomly picked colonies, and 0.43 ± 0.53 random mutations per gene. Screening of the library. The plasmid library was used to transform the Escherichia coli strain BL21(DE )-pGro7/GroEL (TaKaRa) to obtain colonies expressing a library of SsoPox variants. Randomly picked clones (184), representing a coverage of 3.8% of the library, were grown on a 96-well plate in 500 µ L of ZYP medium as previously described . Production of proteins and chaperones was induced after five hours of cul- ture at 37 °C by reducing the temperature to 25 °C, adding CoCl (0.2 mM) and adding arabinose (0.2%, w/v). After overnight growth, cell lysate was used to perform activity screening with 100 µ M of paraoxon substrate 35,49 (Figure S2B) after partial purification of the protein with a heating step of 15 minutes incubation at 70 °C . The screening was performed in 50 mM HEPES pH 8, 150 mM NaCl, 0.2 mM CoCl . Kinetics of paraoxon/ methyl-parathion hydrolysis were monitored by following absorbance at 405 nm for 10 minutes using a microplate reader (Synergy HT, BioTek, USA) and the Gen5.1 software. Saturation library. Saturation mutagenesis directed on residue 263 of the αsB5 variant encoding gene, was per- formed using NNS degenerated primers (W263NNS direct: 5′-TGCACCATTGATNNSGGCACCGCAAAACCG -3′ and W263NNS reverse: 5′-CGGTTTTGCGGTGCCSNNATCAATGGTGCA-3′) and pET32b-ΔTrx carry- ing αsB5 coding gene as template. PCR amplification was performed with 2.5 U of PfuTurbo DNA polymerase (Agilent) according to the manufacturer’s recommendations [95 °C, 5 minutes; 20x (95 °C, 30 seconds; 55 °C, one minute; 68 °C, 15 minutes); 68 °C, 25 minutes]. DNA was digested with DpnI to remove the methylated parental template. E. coli BL21(DE )-pGro7/GroEL (TaKaRa) competent cells were transformed with the plasmid pool and −1 −1 plated on LB-agar supplemented with 100 µ g.mL ampicillin and 34 µ g.mL chloramphenicol. Over a theoret- −1 ical diversity of 20 sequences, 181 variants were collected in a microplate containing LB (100 µg.mL ampicillin −1 and 34 µ g.mL chloramphenicol) for characterization. This oversampling by a factor of nearly 10-fold ensure a >99% probability to get the 20 possible substitutions. Screening was performed as described above (Figure S2C). Clones were produced in liquid ZYP medium and their activity was screened with 100 µM of paraoxon substrate. The plasmids corresponding to the most interest- ing variants were extracted and the SsoPox variant encoding genes were sequenced. Production-purification of SsoPox wild-type and variants. The genes coding for the SsoPox wild-type and its variant were cloned in a pET32b-Δtrx plasmid. Productions were performed using the E. coli BL21(DE )-pGro7/GroEL (TaKaRa) chaperone expressing strain. Starter cultures were produced in an −1 −1 auto-inducible ZYP medium (supplemented with 100 µ g.mL ampicillin and 34 µ g.mL chloramphenicol). When OD reached a value of 0.8–1, CoCl was added (final concentration 0.2 mM) as well as L-arabinose 600nm 2 −1 (final concentration 2 g.L ) to induce the production of chaperones GroEL/ES and the temperature was decreased to 23 °C for 16–20 hours. Cells were harvested by centrifugation (4400 g, 4 °C, 20 minutes) and resus- −1 pended in lysis buffer (50 mM HEPES pH 8.0, 150 mM NaCl, 0.2 mM CoCl , 0.25 mg.mL lysozyme, 0.1 mM −1 PMSF and 10 µ g.mL DNAseI) and were stored at -80 °C. Cells were thawed and lysed in three steps of 30 sec- onds of sonication (Qsonica, Q700; Amplitude 45). Cell debris was removed by centrifugation (20000 g, 4 °C, 15 minutes). As SsoPox and its variants are hyperthermostable, a pre-purification step was performed by heat- ing the lysate for 30 minutes at 80 °C. Precipitated host proteins were removed by centrifugation (20000 g, 4 °C, 15 minutes). SsoPox and its variants were collected by ammonium sulfate precipitation (75%) and resuspended in 8 ml of activity bue ff r (50 mM HEPES pH 8.0, 150 mM NaCl, 0.2 mM CoCl ). The remaining ammonium sulfate was eliminated by injection on a desalting column (HiPrep 26/10 desalting, GE Healthcare; ÄKTA Avant) and TM concentrated to 2 mL for separation on exclusion size chromatography (HiLoad 16/600 Superdex 75 pg, GE Healthcare; ÄKTA Avant). Final purity was monitored by SDS-PAGE and protein concentration was measured with a NanoDrop 2000 spectrophotometer (Thermo Scientific). Kinetic assays. Generalities. Catalytic parameters were evaluated in triplicate at 25 °C, and recorded using a microplate reader (Synergy HT, BioTek, USA) and the Gen5.1 software, in a 6.2 mm path length cell for 200 µ L reaction in 96-well plate, as previously explained . Catalytic parameters were obtained by fitting the data to the Michaelis-Menten (MM) equation using Graph-Pad Prism 6 software. When V could not be reached in the max experiments, catalytic efficiency was obtained by fitting the linear part of MM plot to a linear regression using Graph-Pad Prism 6 software. For some Sso Pox variants, the MM plot was fitted to the substrate inhibition equa - tion using Graph-Pad Prism 6 software enabling us to determine a K for several substrates. Consequently, the calculated catalytic efficiencies in these conditions are true only at low substrate concentrations. In some other cases, saturation could not be reached, therefore k /K values were determined using linear regression. For catM Coumaphos, data were fitted to one-phase decay non-linear regression. Phosphotriesterase activity characterization. e k Th inetic assays were carried out in activity buer ff . The conditions used for pesticides were as follows: ethyl-paraoxon (I), ethyl-parathion (II), methyl-paraoxon (III) and −1 −1 methyl-parathion (IV) hydrolysis was recorded at 405 nm (ε = 17000 M .cm ). The concentration interval was between 0.05 and 2 mM from 100 mM stock solution in ethanol . Malathion (V) hydrolysis was monitored at −1 −1 46 412 nm, (ε = 13400 M .cm ) . The concentration interval was between 0.05 and 2 mM from 100 mM stock solution in DMSO. DTNB was added to the activity buffer. Chlorpyrifos (VI) hydrolysis was recorded at 310 nm −1 −1 51 (ε = 5562 M .cm ) . The concentration interval was between 0.05 and 2 mM from 100 mM stock solution in −1 −1 27 1-propanol. Diazinon (VII) hydrolysis was recorded at 228 nm (ε = 3300 M .cm ) . The concentration interval was between 0.05 and 2 mM from 200 mM stock solution in ethanol. Fenitrothion (VIII) hydrolysis was recorded −1 −1 52 at 358 nm (ε = 18700 M .cm ) . The concentration interval was between 0.05 and 2 mM from 200 mM stock Scientific REPO R TS | 7: 16745 | DOI:10.1038/s41598-017-16841-0 12 www.nature.com/scientificreports/ −1 −1 27 solution in ethanol. Fensulfothion (IX) hydrolysis was recorded at 284 nm (ε = 8000 M .cm ) . The concentra- tion interval was between 0.05 and 2 mM from 200 mM stock solution in methanol. Coumaphos (X) hydrolysis was monitored by fluorescence of released chlorferon (excitation 360/40 and emission 460/40) . A linear corre- lation between fluorescence and chlorferon was found (Figure S3). Measurements were taken at 5 µ M from a 200 mM stock solution in DMSO. We assumed that , thus enabling us to estimate . KS [] k /K M catM Lactonase activity characterization. Lactonase kinetics were performed using a previously described protocol . e Th time course hydrolysis of lactones were performed in lac buffer (2.5 mM Bicine pH 8.3, 150 mM NaCl, 0.2 mM CoCl , 0.25 mM Cresol purple and 0.5% DMSO) over a concentration range 0–2 mM for 3-oxo-C10 AHLs (XI) or 0–5 mM for undecanoic- γ/δ -lactones (XII, XIII). Cresol purple (pK 8.3 at 25 °C) is a pH indicator used to fol- low lactone ring hydrolysis by acidification of the medium. Molar coefficient extinction at 577 nm was evaluated recording absorbance of the buffer over an acetic acid range of concentration 0–0.35 mM. Bioremediation of pesticide solutions at 250 µM. For all OP substrates but malathion, experiments were performed in triplicate at 25 °C using a microplate reader (Synergy HT, BioTek, USA). Wavelengths were chosen as described in the kinetic assays section. Degradation of malathion was followed by GC/MS in triplicate. OP concentration was fixed at 250 µ M. All the substrates were soluble at this concentration. The best variants of SsoPox for each substrate were used, namely αsD6 for all the OPs. Two enzyme-to-substrate ratios were [E]/[S] −2 −3 −1 −1 used, 10 and 10 , corresponding to enzyme concentrations of 90 µg.ml and 9 µg.ml respectively. The reac- tions were monitored until the plateau was reached. Experimental measures were obtained using Gen5.1 sowa ft re, then analyzed with GraphPad Prism 6 software. Curves were then fitted using One-Phase Decay non-linear regression with the equation (1): () −kt YY =− (0 Plateaue ) ×+ Plateau (1) where Y0 = 0% and Plateau = 100%. The rate constant k was determined and the time required to observe a 95% decontamination was calculated accordingly. The curves and results of the fits are shown in Figure S4. −3 Degradation at 95% was confirmed using a final point by GC/MS for all pesticides with a 10 enzyme-to-substrate ratio. 100 µ L of activity buer ff solution was first extracted with 100 µ L chloroform. Organic extracts were analyzed by using a Clarus 500 gas chromatograph equipped with a SQ8S MS detector (Perkin −1 Elmer, Courtaboeuf, France). 1 µ L of organic extract was volatilized at 220 °C (split 15 mL.min ) in a deactivated FocusLiner with quartz wool (SGE, Ringwood, Australia) and compounds separated on an Elite-5MS column (30 m, −1 0.25 mm i.d., 0.25 mm film thickness) for 12 minutes using a temperature gradient (80–280 °C at 30 °C.min , five −1 minutes’ hold). Helium flowing at 2 mL.min was used as the carrier gas. The MS inlet line was set at 280 °C and the electron ionization source at 280 °C and 70 eV. Full scan monitoring was performed from 40 to 400 m/z in order to identify chemicals by spectral database search using MS Search 2.0 operated with the Standard Reference Database 1 A (National Institute of Standards and Technology, Gaithersburg, MD, USA). m/z is the mass-to-charge ratio of the base peak fragment detected for each molecule. In the case of weak pesticide signals, extracted ion chromato- grams were generated with base peak ions to confirm the presence of chemicals (Paraoxon (I ) and Parathion ethyl (II) m/z 109; Malathion (V) and Fenitrothion (VIII) m/z 125; Diazinon (VII) m/z 137; Coumaphos (X) m/z 97; Fensulfothion (IX) m/z 293; and Chlorpyrifos (VI) m/z 197). Selected Ion Recording using base peaks ions was applied in order to specifically monitor pesticides and collect peak areas for kinetics. All samples were analyzed over short periods of time to avoid signal drift. All data were processed using Turbomass 6.1 (Perkin Elmer). 42,46 Crystallization. Crystallization assays were performed as previously described . Crystallization was performed using the hanging drop vapor diffusion method in 96-well plates (Greiner Microplate, 96 well, PS, F-bottom) on ViewDrop II seals (TPP Labtech). Equal volumes (0.5 µ L) of protein and reservoir solutions were mixed, and the resulting drops were equilibrated against a 150 µ L reservoir solution containing 20–30% (w/v) PEG 8000 and 50 mM Tris-HCl buffer (pH 8). Crystals appeared aer f ft ew days at 4 °C. Data collection and structure refinement. Crystals were first transferred to a cryoprotectant solution con- sisting of the reservoir solution and 20% (v/v) glycerol. Crystals were then flash-cooled in liquid nitrogen. X-ray diffraction data were collected at 100 K using synchrotron radiation at the ID23-1 beam line (ESRF, Grenoble, France) and using an ADSC Q315r detector. X-ray diffraction data were integrated and scaled with the XDS package (Table 1) . The phases were obtained using the native structure of SsoPox (PDB code 2vc5) as a starting model, 62,63 performing a molecular replacement with MOLREP or PHASER . The models were built with Coot and refined 60,64 using REFMAC . We note that three structures presented in this work exhibit one disordered monomer (and low corresponding electronic density): αsB5 (monomer D, total of 4 monomers), αsD6 (monomer C, out of a total of 4 monomers) and αsC6 (monomer C, total of 4 monomers). Structure illustrations were performed using PyMOL . Relative B-factor analysis. e o Th ccupancies of all residues in all tested structures were set to 1 for this anal- ysis. For residues with alternate conformations, the sums of occupancies were set to 1. Structures were re-refined 60,64 with REFMAC . The relative B-factor values were obtained by normalizing the B-factor values of each residue 11,54,57 by the average B-factor of the whole structure as previously described . References 1. Wolfenden, R. Benchmark Reaction Rates, the Stability of Biological Molecules in Water, and the Evolution of Catalytic Power in Enzymes. Annu. Rev. Biochem. 80, 645–667 (2011). Scientific REPO R TS | 7: 16745 | DOI:10.1038/s41598-017-16841-0 13 www.nature.com/scientificreports/ 2. Singh, R. K., Tiwari, M. K., Singh, R. & Lee, J.-K. From Protein Engineering to Immobilization: Promising Strategies for the Upgrade of Industrial Enzymes. Int. J. Mol. Sci. 14, 1232–1277 (2013). 3. Bar-Even, A. et al. The Moderately Efficient Enzyme: Evolutionary and Physicochemical Trends Shaping Enzyme Parameters. Biochemistry (Mosc.) 50, 4402–4410 (2011). 4. Tokuriki, N. et al. Diminishing returns and tradeoffs constrain the laboratory optimization of an enzyme. Nat. Commun. 3, 1257 (2012). 5. Burton, S. G., Cowan, D. A. & Woodley, J. M. The search for the ideal biocatalyst. Nat. Biotechnol. 20, 37–45 (2002). 6. Gupta, R. D. & Tawfik, D. S. Directed enzyme evolution via small and effective neutral drift libraries. Nat. Methods 5, 939–942 (2008). 7. Bommarius, A. S. Biocatalysis: A Status Report. Annu. Rev. Chem. Biomol. Eng. 6, 319–345 (2015). 8. Currin, A., Swainston, N., Day, P. J. & Kell, D. B. Synthetic biology for the directed evolution of protein biocatalysts: navigating sequence space intelligently. Chem Soc Rev 44, 1172–1239 (2015). 9. Reetz, M. T. Directed Evolution of Selective Enzymes: Catalysts for Organic Chemistry and Biotechnology. (Wiley-VCH Verlag GmbH & Co. KGaA, 2016), https://doi.org/10.1002/9783527655465. 10. Directed Enzyme Evolution: Advances and Applications. (Springer International Publishing, 2017). https://doi.org/10.1007/978-3- 319-50413-1. 11. Hiblot, J., Gotthard, G., Elias, M. & Chabriere, E. Differential Active Site Loop Conformations Mediate Promiscuous Activities in the Lactonase SsoPox. PLoS ONE 8, e75272 (2013). 12. Reetz, M. T. Laboratory Evolution of Stereoselective Enzymes: A Prolific Source of Catalysts for Asymmetric Reactions. Angew. Chem. Int. Ed. 50, 138–174 (2011). 13. Daudé, D., Topham, C. M., Remaud‐Siméon, M. & André, I. Probing impact of active site residue mutations on stability and activity of Neisseria polysaccharea amylosucrase. Protein Sci. 22, 1754–1765 (2013). 14. Alcolombri, U., Elias, M. & Tawfik, D. S. Directed Evolution of Sulfotransferases and Paraoxonases by Ancestral Libraries. J. Mol. Biol. 411, 837–853 (2011). 15. Gonzalez, D. et al. Ancestral mutations as a tool for solubilizing proteins: The case of a hydrophobic phosphate-binding protein. FEBS Open Bio 4, 121–127 (2014). 16. Hawwa, R., Larsen, S. D., Ratia, K. & Mesecar, A. D. Structure-Based and Random Mutagenesis Approaches Increase the Organophosphate-Degrading Activity of a Phosphotriesterase Homologue from Deinococcus radiodurans. J. Mol. Biol. 393, 36–57 (2009). 17. Herman, A. & Tawfik, D. S. Incorporating Synthetic Oligonucleotides via Gene Reassembly (ISOR): a versatile tool for generating targeted libraries. Protein Eng. Des. Sel. 20, 219–226 (2007). 18. Rochu, D., Chabrière, E. & Masson, P. Human paraoxonase: A promising approach for pre-treatment and therapy of organophosphorus poisoning. Toxicology 233, 47–59 (2007). 19. Krieger, R. Handbook of Pesticide Toxicology, Two-Volume Set: Principles and Agents. (Academic Press, 2001). 20. Gupta, R. C. Handbook of Toxicology of Chemical Warfare Agents. (Academic Press, 2009). 21. Sapozhnikova, Y. et al. Evaluation of pesticides and metals in fish of the Dniester River, Moldova. Chemosphere 60, 196–205 (2005). 22. Sapozhnikova, Y., Bawardi, O. & Schlenk, D. Pesticides and PCBs in sediments and fish from the Salton Sea, California, USA. Chemosphere 55, 797–809 (2004). 23. LeJeune, K. E., Wild, J. R. & Russell, A. J. Nerve agents degraded by enzymatic foams. Nature 395, 27–28 (1998). 24. Jacquet, P. et al. Current and emerging strategies for organophosphate decontamination: special focus on hyperstable enzymes. Environ. Sci. Pollut. Res. 1–19 https://doi.org/10.1007/s11356-016-6143-1 (2016). 25. Poirier, L., Jacquet, P., Elias, M., Daudé, D. & Chabrière, E. Decontamination of organophosphorus compounds: Towards new alternatives. Ann. Pharm. Fr. https://doi.org/10.1016/j.pharma.2017.01.004 (2017). 26. Afriat, L., Roodveldt, C., Manco, G. & Tawfik, D. S. The Latent Promiscuity of Newly Identified Microbial Lactonases Is Linked to a Recently Diverged Phosphotriesterase†. Biochemistry (Mosc.) 45, 13677–13686 (2006). 27. Dumas, D. P., Caldwell, S. R., Wild, J. R. & Raushel, F. M. Purification and properties of the phosphotriesterase from Pseudomonas diminuta. J. Biol. Chem. 264, 19659–19665 (1989). 28. Seibert, C. M. & Raushel, F. M. Structural and Catalytic Diversity within the Amidohydrolase Superfamily†. Biochemistry (Mosc.) 44, 6383–6391 (2005). 29. Roodveldt, C. & Tawfik, D. S. Shared Promiscuous Activities and Evolutionary Features in Various Members of the Amidohydrolase Superfamily. Biochemistry (Mosc.) 44, 12728–12736 (2005). 30. Hiblot, J., Gotthard, G., Chabriere, E. & Elias, M. Structural and Enzymatic characterization of the lactonase SisLac from Sulfolobus islandicus. PLoS ONE 7, e47028 (2012). 31. Porzio, E., Merone, L., Mandrich, L., Rossi, M. & Manco, G. A new phosphotriesterase from Sulfolobus acidocaldarius and its comparison with the homologue from Sulfolobus solfataricus. Biochimie 89, 625–636 (2007). 32. Merone, L., Mandrich, L., Rossi, M. & Manco, G. A thermostable phosphotriesterase from the archaeon Sulfolobus solfataricus: cloning, overexpression and properties. Extremophiles 9, 297–305 (2005). 33. Hawwa, R., Aikens, J., Turner, R. J., Santarsiero, B. D. & Mesecar, A. D. Structural basis for thermostability revealed through the identification and characterization of a highly thermostable phosphotriesterase-like lactonase from Geobacillus stearothermophilus. Arch. Biochem. Biophys. 488, 109–120 (2009). 34. Xiang, D. F. et al. Functional Identification of Incorrectly Annotated Prolidases from the Amidohydrolase Superfamily of Enzymes. Biochemistry (Mosc.) 48, 3730–3742 (2009). 35. Chow, J. Y. et al. Directed Evolution of a Thermostable Quorum-quenching Lactonase from the Amidohydrolase Superfamily. J. Biol. Chem. 285, 40911–40920 (2010). 36. Hiblot, J., Gotthard, G., Champion, C., Chabriere, E. & Elias, M. Crystallization and preliminary X-ray diffraction analysis of the lactonase VmoLac from Vulcanisaeta moutnovskia. Acta Crystallograph. Sect. F Struct. Biol. Cryst. Commun. 69, 1235–1238 (2013). 37. Hiblot, J., Bzdrenga, J., Champion, C., Chabriere, E. & Elias, M. Crystal structure of VmoLac, a tentative quorum quenching lactonase from the extremophilic crenarchaeon Vulcanisaeta moutnovskia. Sci. Rep. 5, 8372 (2015). 38. Bzdrenga, J. et al. SacPox from the thermoacidophilic crenarchaeon Sulfolobus acidocaldarius is a proficient lactonase. BMC Res. Notes 7, 333 (2014). 39. Bzdrenga, J. et al. Biotechnological applications of quorum quenching enzymes. Chem. Biol. Interact. https://doi.org/10.1016/j. cbi.2016.05.028 (2016). 40. Rémy, B., Plener, L., Elias, M., Daudé, D. & Chabrière, E. Des enzymes pour bloquer la communication bactérienne, une alternative aux antibiotiques? Ann. Pharm. Fr. https://doi.org/10.1016/j.pharma.2016.06.005 (2016). 41. Vieille, C. & Zeikus, G. J. Hyperthermophilic Enzymes: Sources, Uses, and Molecular Mechanisms for Thermostability. Microbiol. Mol. Biol. Rev. 65, 1–43 (2001). 42. Elias, M. et al. Structural Basis for Natural Lactonase and Promiscuous Phosphotriesterase Activities. J. Mol. Biol. 379, 1017–1028 (2008). 43. Hoque, M. A. et al. Stepwise Loop Insertion Strategy for Active Site Remodeling to Generate Novel Enzyme Functions. ACS Chem. Biol. https://doi.org/10.1021/acschembio.7b00018 (2017). Scientific REPO R TS | 7: 16745 | DOI:10.1038/s41598-017-16841-0 14 www.nature.com/scientificreports/ 44. Meier, M. M. et al. Molecular Engineering of Organophosphate Hydrolysis Activity from a Weak Promiscuous Lactonase Template. J. Am. Chem. Soc. 135, 11670–11677 (2013). 45. Elias, M. et al. Crystallization and preliminary X-ray diffraction analysis of the hyperthermophilic Sulfolobus solfataricus phosphotriesterase. Acta Crystallograph. Sect. F Struct. Biol. Cryst. Commun. 63, 553–555 (2007). 46. Hiblot, J., Gotthard, G., Chabriere, E. & Elias, M. Characterisation of the organophosphate hydrolase catalytic activity of SsoPox. Sci. Rep. 2 (2012). 47. Del Giudice, I. et al. An efficient thermostable organophosphate hydrolase and its application in pesticide decontamination. Biotechnol. Bioeng. 113, 724–734 (2016). 48. Rémy, B. et al. Harnessing hyperthermostable lactonase from Sulfolobus solfataricus for biotechnological applications. Sci. Rep. 6 (2016). 49. Vecchio, P. D. et al. Structural determinants of the high thermal stability of SsoPox from the hyperthermophilic archaeon Sulfolobus solfataricus. Extremophiles 13, 461–470 (2009). 50. Afriat-Jurnou, L., Jackson, C. J. & Tawfik, D. S. Reconstructing a Missing Link in the Evolution of a Recently Diverged Phosphotriesterase by Active-Site Loop Remodeling. Biochemistry (Mosc.) 51, 6047–6055 (2012). 51. Jackson, C. J. et al. Anomalous scattering analysis of Agrobacterium radiobacter phosphotriesterase: the prominent role of iron in the heterobinuclear active site. Biochem. J. 397, 501 (2006). 52. Tawfik, D. S. Accuracy-rate tradeoffs: how do enzymes meet demands of selectivity and catalytic efficiency? Curr. Opin. Chem. Biol. 21, 73–80 (2014). 53. Xiang, D. F. et al. Functional Annotation and Three-Dimensional Structure of Dr0930 from Deinococcus radiodurans: A Close Relative of Phosphotriesterase in the Amidohydrolase Superfamily. Biochemistry (Mosc.) 48, 2237–2247 (2009). 54. Dellus-Gur, E. et al. Negative Epistasis and Evolvability in TEM-1 β-Lactamase—The Thin Line between an Enzyme’s Conformational Freedom and Disorder. J. Mol. Biol. 427, 2396–2409 (2015). 55. Reetz, M. T. The Importance of Additive and Non-Additive Mutational Effects in Protein Engineering. Angew. Chem. Int. Ed. 52, 2658–2666 (2013). 56. Campbell, E. et al. The role of protein dynamics in the evolution of new enzyme function. Nat. Chem. Biol. 12, 944–950 (2016). 57. Dellus-Gur, E., Toth-Petroczy, A., Elias, M. & Tawfik, D. S. What Makes a Protein Fold Amenable to Functional Innovation? Fold Polarity and Stability Trade-offs. J. Mol. Biol. 425, 2609–2621 (2013). 58. Guendouze, A. et al. Effect of quorum quenching lactonase in clinical isolates of Pseudomonas aeruginosa and comparison with quorum sensing inhibitors. Front. Microbiol. 8 (2017). 59. DeLano, W. L. The PyMOL molecular graphics system (2002). 60. Emsley, P. & Cowtan, K. Coot: model-building tools for molecular graphics. Acta Crystallogr. D Biol. Crystallogr. 60, 2126–2132 (2004). 61. Kabsch, W. Automatic processing of rotation diffraction data from crystals of initially unknown symmetry and cell constants. J. Appl. Crystallogr. 26, 795–800 (1993). 62. McCoy, A. J. et al. Phaser crystallographic software. J. Appl. Crystallogr. 40, 658–674 (2007). 63. Vagin, A. & Teplyakov, A. An approach to multi-copy search in molecular replacement. Acta Crystallogr. D Biol. Crystallogr. 56, 1622–1624 (2000). 64. Murshudov, G. N., Vagin, A. A. & Dodson, E. J. Refinement of Macromolecular Structures by the Maximum-Likelihood Method. Acta Crystallogr. D Biol. Crystallogr. 53, 240–255 (1997). Acknowledgements This work was granted by DGA, France (REI. 2009 34 0045). G.G. and P.J. are PhD students granted by DGA. Financial support was also obtained from the MnDrive Initiative (M.E.) to support this work. This work was supported by “Investissements d’avenir” program (Méditerranée Infection 10-IAHU-03) of the French Agence Nationale de la Recherche (ANR). Author Contributions P.J., J.H., D.D., C.B., N.A., G.G. performed the experiments, P.J., J.H., D.D., C.B., N.A., G.G., E.C., M.E. analyzed the data, P.J., J.H., D.D., C.B., M.E. wrote the manuscript and prepared the figures, E.C., D.D. and M.E. designed the experiments. Additional Information Supplementary information accompanies this paper at https://doi.org/10.1038/s41598-017-16841-0. Competing Interests: M.E. and E.C. have a patent WO2014167140 A1 licensed to Gene&GreenTK. D.D. and E.C. report personal fees from Gene&GreenTK during the conduct of the study. Publisher's note: Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations. Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Cre- ative Commons license, and indicate if changes were made. The images or other third party material in this article are included in the article’s Creative Commons license, unless indicated otherwise in a credit line to the material. If material is not included in the article’s Creative Commons license and your intended use is not per- mitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this license, visit http://creativecommons.org/licenses/by/4.0/. © The Author(s) 2017 Scientific REPO R TS | 7: 16745 | DOI:10.1038/s41598-017-16841-0 15 http://www.deepdyve.com/assets/images/DeepDyve-Logo-lg.png Scientific Reports Springer Journals

Rational engineering of a native hyperthermostable lactonase into a broad spectrum phosphotriesterase

Download PDF
15 pages

Loading next page...
1
 
/lp/springer_journal/rational-engineering-of-a-native-hyperthermostable-lactonase-into-a-9GkDpujNzZ

References (54)

Publisher
Springer Journals
Copyright
Copyright © 2017 by The Author(s)
Subject
Science, Humanities and Social Sciences, multidisciplinary; Science, Humanities and Social Sciences, multidisciplinary; Science, multidisciplinary
eISSN
2045-2322
DOI
10.1038/s41598-017-16841-0
Publisher site
See Article on Publisher Site

Abstract

www.nature.com/scientificreports OPEN Rational engineering of a native hyperthermostable lactonase into a broad spectrum phosphotriesterase Received: 29 August 2017 1 1,4 2 1,3 1 Pauline Jacquet , Julien Hiblot , David Daudé , Céline Bergonzi , Guillaume Gotthard , 1 1 3 Accepted: 16 November 2017 Nicholas Armstrong , Eric Chabrière & Mikael Elias Published: xx xx xxxx The redesign of enzyme active sites to alter their function or specificity is a difficult yet appealing challenge. Here we used a structure-based design approach to engineer the lactonase SsoPox from Sulfolobus solfataricus into a phosphotriesterase. The five best variants were characterized and their structure was solved. The most active variant, αsD6 (V27A-Y97W-L228M-W263M) demonstrates a large increase in catalytic efficiencies over the wild-type enzyme, with increases of 2,210-fold, 163-fold, 58- fold, 16-fold against methyl-parathion, malathion, ethyl-paraoxon, and methyl-paraoxon, respectively. Interestingly, the best mutants are also capable of degrading fensulfothion, which is reported to be an inhibitor for the wild-type enzyme, as well as others that are not substrates of the starting template or previously reported W263 mutants. The broad specificity of these engineered variants makes them promising candidates for the bioremediation of organophosphorus compounds. Analysis of their structures reveals that the increase in activity mainly occurs through the destabilization of the active site loop involved in substrate binding, and it has been observed that the level of disorder correlates with the width of the enzyme specificity spectrum. This finding supports the idea that active site conformational flexibility is essential to the acquisition of broader substrate specificity. Enzymes are extraordinary molecules that can accelerate chemical reactions, including those that spontaneously occur extremely slowly, and can make them happen many times per second . Some of these reactions are of high biotechnological or industrial interest and enzymes are, therefore, the focus of numerous studies . However, enzymes are often not directly compatible with industrial processes for various reasons ranging from their 3,4 5 sub-optimality or lack of stability under certain conditions . The properties of enzymes can be tuned by molec- 5–10 ular engineering to fit these new requirements using high throughput methods such as directed evolution and 11–13 /or low throughput methods, reducing the size of the library by focusing on positions that can be saturated , or 14,15 16 17 by using phylogenetically-inferred or structure-based mutations to create combinatorial , reduced libraries . Here, we focus on the molecular engineering of an enzymatic bio-decontaminant of organophosphates com- pounds (OPs), including nerve agents and insecticides. These compounds are neurotoxic because of their ability to irreversibly inhibit the acetylcholine esterase, a key enzyme in the nervous system . Massively used in agri- culture as pesticides , they were modified before World War II to increase their toxicity and develop chemical warfare agents . OPs still represent major environmental, health and security concerns, as they cause significant 21,22 19,20 pollution , and intoxication, and have been used in terrorist attacks . However, since current methods for removing them are cost prohibitive and lead to ecological concerns, enzymatic remediation is a promising way of 23–25 developing efficient bio-decontamination methods . Phosphotriesterase-like Lactonases (PLLs) are appealing candidates, being natural lactonases endowed with promiscuous phosphotriesterase activity . They have long been related to the bacterial phosphotriesterases (PTEs), enzymes which are capable of hydrolyzing insecticides such as parathion with high catalytic efficiency . 28,29 PTEs and PLLs belong to the amidohydrolase superfamily , and share the same (β/α) barrel topology. Some CNRS UMR 7278, IRD198, INSERM U1095, APHM, Institut Hospitalier Universitaire Méditerranée-Infection, Aix- Marseille Université, 19-21 Bd Jean Moulin, 13005, Marseille, France. Gene&GreenTK, IHU Méditerranée Infection, 19-21 Bd Jean Moulin, 13005, Marseille, France. University of Minnesota, Department of Biochemistry, Molecular Biology and Biophysics & Biotechnology Institute, St. Paul, MN, 55108, USA. Present address: MPI for Medical Research, Chemical Biology department (EPFL), Heidelberg, Germany. Pauline Jacquet, Julien Hiblot and David Daudé contributed equally to this work. Correspondence and requests for materials should be addressed to E.C. (email: eric.chabriere@univ-amu.fr) or M.E. (email: mhelias@umn.edu) Scientific REPO R TS | 7: 16745 | DOI:10.1038/s41598-017-16841-0 1 www.nature.com/scientificreports/ Data collections asA6 open asA6 closed asB5 asC6 asD6 asA1 PDB ID 5VRK 5VRI 5W3U 5W3Z 5W3W 5VSA Wavelength (Å) 0.97934 0.93340 0.97625 1.00882 1.00882 0.99987 Resolution (Å) 1.4 2.15 2.5 2.55 2.95 2.0 Space group P2 P2 2 2 P2 2 2 P2 2 2 P2 2 2 P2 2 2 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 Unit cell dimensions a (Å) 49.28 87.17 82.93 81.65 82.34 86.6 b (Å) 137.20 103.61 105.95 105.80 105.01 103.100 c (Å) 49.36 152.61 152.83 154.53 153.15 151.8 α (°) 90.0 β (°) 98.7 γ (°) 90.0 No. observed reflections 419571 (75591) 408826 (24712) 206443 (23317) 194083 (10849) 106247 (5379) 703275 (96197) No. unique reflections 121915 (22376) 74124 (4744) 47092 (5182) 43204 (2416) 28163 (1372) 91103 (12398) Completeness (%) 96.1 (94.3) 97.8 (95.5) 99.5 (99.8) 97.3 (98.0) 98.3 (99.9) 98.6% (99.7%) Rmerge (%) 8.2 (33.4) 9.8 (44.6) 6.1 (48.7) 5.55 (42.8) 6.6 (46.0) 9.0% (48.8%) Rmeasure (%) 9.7 (39.5) 10.8 (49.3) 6.9 (55.4) 6.2 (48.5) 7.7 (52.9) 9.6% (52.6%) I/σ(I) 9.07 (3.39) 14.31 (3.47) 12.84 (2.50) 15.68 (2.67) 12.89 (2.75) 15.59 (4.23) Last resolution shell 1.5–1.4 2.2–2.15 2.6–2.5 2.6–2.55 3.0–2.95 2.1–2.0 Redundancy 3.44 (3.38) 5.51 (5.21) 4.38 (4.50) 4.49 (4.49) 3.77 (3.92) 7.72 (7.76) Refinement statistics Resolution range (Å) 48.7128–1.40 49.0556–2.1500 49.6449–2.500 46.313–2.55 49.667–2.95 48.812–2.0 No. Reflections 115818 70416 44736 41042 26754 86547 Rfree /Rwork 16.84 / 13.36 21.35/17.23 27.17/21.22 24.15/19.97 28.31/24.68 20.55/16.32 No. protein atoms 5368 10122 9310 9657 9683 9987 No. water molecules 634 659 76 154 23 733 Average B factor (Å2) 19.095 24.984 86.6 80.375 95.992 29.885 rmsd from ideal Bond lengths (Å) 0.0088 0.0028 0.0050 0.0110 0.0070 0.0211 Bond angles (°) 1.0703 0.7309 0.8490 1.1730 0.9630 0.6950 Table 1. Data collection and refinement statistics of SsoPox variants structures. 30–38 PLL representatives are extremely thermostable and are particularly promising candidates for biotechno- 39,40 41 logical use by virtue of their high stability and resistance towards denaturing agents . PLLs, however, exhibit much lower phosphotriesterase activity than PTEs. Engineering experiments with the aim of increasing PLLs’ phosphotriesterase activity were performed. In particular, because PTEs and PLLs’ active sites mainly differ in 16,26,42 their loop length (7 and 8), sequences and conformations , loop grafting experiments were performed and 16,34 proved to be difficult. This approach resulted in aggregation-prone or inactive enzymes , but was successful when using a stepwise approach . Recently, more classical engineering techniques were successful in increasing the phosphotriesterase activity of the PLL DrOPH . 32,45 In this study, we worked on the PLL SsoPox from the hyperthermophilic archaea Sulfolobus solfataricus . SsoPox exhibits high acyl-homoserine lactones and oxo-lactones hydrolysis abilities , and promiscuous, low phosphotriesterase activity . Previous work showed that engineering can successfully improve SsoPox catalytic 11,47 efficiency against phosphotriesters such as ethyl-paraoxon . These previous studies highlighted the role of a key active site loop 8 residue, W263, which modulates the active site loop conformational space . In this study, we aimed to further improve the phosphotriesterase activity of SsoPox, including against phosphothionoesters. SsoPox is a very appealing candidate for the bioremediation of phosphotriesters because of its unique thermal stability, and its ability to resist aging, solvent and protease treatments . On the other hand, SsoPox is a poor phosphothionoesterase, exhibiting a strong preference (100-fold) for methyl-paraoxon over methyl-parathion, a feature referred to as the thiono-effect . We used a rational strategy: taking advantage of the structural similarity between SsoPox and the bacterial PTE from the mesophilic Brevundimonas diminuta (BdPTE) , we designed and produced a structure-based combinatorial library, and screened this library for improved phosphotriesterase activity with the aim of obtaining heat stable, highly active variants with broad specificity. We obtained sev - eral improved variants, including against methyl-parathion, with a 2,210-fold improvement in activity. All these variants but variant αsB5 contained a mutation at position 263. Interestingly, the addition of the W263 muta- tion using site-saturation mutagenesis on an αsB5 background revealed an incompatibility of these mutations, instead of the expected increase in phosphotrieseterase activity. Of these improved variants, we obtained a mutant 4 −1 −1 that lost the thiono-effect against methyl-parathion, exhibiting a > 2,000-fold (k /K ~10 M .s ) as referred cat M to the wild-type enzyme, and a methyl-paraoxon/methyl-parathion activity ratio of ~1. The X-ray structures of several improved variants reveal that the selected mutations increase the active site loop conformational flexibility and reshape the active site. The best variants were extensively characterized against eight commercially available Scientific REPO R TS | 7: 16745 | DOI:10.1038/s41598-017-16841-0 2 www.nature.com/scientificreports/ Figure 1. Structure-based design approach. (A) Pipeline used for the selection of improved mutants. (B) Structural alignments of SsoPox (red) and BdPTE structures (yellow) guided the mutational design. (C) Active site superposition of BdPTE (green) and SsoPox (cyan) active sites. Residues with structural equivalent (orange sticks) in both structures were used to design mutations of SsoPox’s residues into the BdPTE corresponding residues. Positions with no structural equivalent (red sticks) were designed using modelling tools (purple sticks) to mimic the active site cavity size, shape and chemical properties of the BdPTE crystal structure model (Figure S1). (D) The mutations data base consists of mutations to the BdPTE sequence when there is a structural equivalence (orange), and of mutations designed when there is no structural equivalence (red). The third set of mutations (black) adds more diversity at two selected positions, 258 and 263. (E) Once the library had been validated, primers carrying mutations were used to shuffle them and generate a gene library with random combinations of selected mutations. (F) The library was screened for paraoxonase activity to identify enzymes with improved proficiency. insecticides, along with previously reported SsoPox monovariants harboring mutation at residue 263 and their specificity spectrum was determined. Results and Discussion e l Th actonase SsoPox was engineered for higher phosphotriesterase activity using structure-based combinatorial libraries. By comparing structures of enzymes with similar topology, it has been possible to redesign, using mod- elling tools, the active site cavity of SsoPox to mimic as closely as possible that of BdPTE (Fig. 1; Table S1). This method consisted of two main steps: (i) the identification of mutations corresponding to equivalent positions at the enzyme active sites using structural alignment and (ii) the rational selection of mutations, in a position with no equivalent residue, which mimicked the shape and chemical nature of the target enzyme cavity. A mutation dataset was thus obtained and used to develop a combinatorial library consisting of random combinations of our pre-selected mutations. Scientific REPO R TS | 7: 16745 | DOI:10.1038/s41598-017-16841-0 3 www.nature.com/scientificreports/ Figure 2. Chemical structure of substrates used in this study. The chemical structure of ethyl-paraoxon (I), ethyl-parathion (II), methyl-paraoxon (III), methyl-parathion (IV), malathion (V), chlorpyrifos (VI), diazinon (VII), fenitrothion (VIII), fensulfothion (IX), coumaphos (X), 3-oxo-C10 AHL (l) (XI), undecanoic-γ-lactone (XII) and undecanoic-δ-lactone (XIII) generated using ChemDraw software. Structure-based design produced mutants with improved phosphotriesterase activity. The structure-based design considerably reshaped the active site cavity of SsoPox, theoretically bringing it closer to that of BdPTE, both in shape and nature (Figure S1). A total of 14 key positions within the active site were iden- tified, and were mutated to specific residues, or degenerated into several possible amino acids (see Methods). The paraoxonase activity of 184 randomly selected clones was screened and the 14 best were sequenced. Most of these 14 were found to possess the same mutations (e.g. W263F/L/M/A; 73% of cases) and (Figure S2A; Table S3). Of them, the most active variants, αsA1 (C258L-I261F-W263A), αsA6 (F46L-C258A-W263M-I280T), αsB5 (V27A-I76T-Y97W-Y99F-L130P-L226V), αsC6 (L72I-Y99F-I122L-L228M-F229S-W263L), and αsD6 (V27A-Y97W-L228M-W263M) were subject to kinetic and structural characterization. It should be noted that αsA1 was obtained in a previous engineering effort , and that four of the five selected variants contained a sub- stitution of the key residue W263. We previously highlighted the role of this residue in substrate binding and modulation of the active site loop 8 conformation, resulting in an increase in activity against promiscuous sub- strates . Variant αsB5 does not harbor any substitution at position W263. We therefore decided to test the W263 substitution on an αsB5 background, using saturation mutagenesis strategy. The five variants (α sA1, αsA6, αsB5, αsC6 & αsD6) were characterized against several phosphotriesters, including ethyl-paraoxon (I), ethyl-parathion (II), methyl-paraoxon (III), methyl-parathion (IV) and mala- thion (V) (Fig. 2). Overall, all variant exhibits improved catalytic efficiencies against all tested substrates, with the exception of αsA1 with malathion (~2-fold decrease) (Fig. 3A; Table 2). e Th best improved variant, αsD6, exhib- its a 2,210-fold increase in catalytic efficiency against methyl-parathion (Table  2). The largest improvements were observed for ethyl/methyl-parathion, two bad substrates for wild-type SsoPox. Interestingly, while the wt enzyme shows a clear preference for small substituents , most selected variants lost this preference, possibly indicating an enlargement of the active site cavity. These improvements were in the range of what was previously obtained but with higher intensive mutation protocols and for only one OP substrate . Finally, it was noted that two variants (αsA6 & αsC6) presented a substrate inhibition for some thiono-phosphotriesters. The five selected variants were further evaluated for their lactonase activity (Table  2). Interestingly, variants αsA6 and αsC6 exhibit enhanced lactonase activity against γ- and δ-lactones, while αsA1, αsB5 and αsD6 show reduced lactonase catalytic efficiencies. This emphasizes the fact that the improvement of the phosphotriesterase activity does not necessarily compromise the cognate, lactonase activity of the enzyme. These mutations there- fore dramatically increased a new activity (phosphotriesterase) without the complete loss of the native/original function (lactonase) . These mutants are die ff rent from those generated on another PLLs, through loop inser - tion. While the stepwise insertion of residues in loop 7 can significantly increase the ability of PLLs to degrade phosphotriesters, it also drastically decreases their lactonase activity . Insertion in loop 7 has been previously described as a key evolutionary event in the transition from lactonase to phosphotriesterase . Scientific REPO R TS | 7: 16745 | DOI:10.1038/s41598-017-16841-0 4 www.nature.com/scientificreports/ Figure 3. Activity characterization of the improved variants and bioremediation potential. (A) Phosphotriesterase catalytic efficiency comparison between wt-SsoPox and the best selected variants. Data for wt-SsoPox are from Hiblot et al. . (B) Determination of the reaction time for hydrolyzing 95% using the αsD6 −2 −3 variant of a solution of pesticide at 250 µM. Two enzyme-to-substrate ratios ( [E]/[S]), 10 and 10 , were considered. Structural analysis of improved variants reveals increased mobility of the active site loop 8. e cr Th ystal structures of variants α sA1, αsA6, αsB5, αsC6 and αsD6 were solved (Table 1). While most structures were solved at high or medium resolutions (1.4–2.55 Å), mutant asD6 could only be solved at low resolution (2.95 Å). Given this resolution, we limited our interpretations and retained this data for the consistency of the study as it exhibits structural features (loop 8 disorder) that are consistent with all the other mutants. Overall, the mutants’ structures are similar to that of wt-SsoPox. However, the active site cavities are larger for the mutants: most selected mutations replace residues with smaller ones: W263 is mutated into M/L/I, V27 into A, I76 into T, L130 into P, Y99 into F. Only two selected substitutions relate to bulkier side chains: Y97W and L228M. This enlargement is difficult to quantify even with the structure of these mutants because of the extremely mobile nature of loop 8 in the variants. Enlarging the active site was a main objective of the active site redesign of SsoPox, to enable the bulkier phosphotriesters to bind within the natural lactonase enzyme. Additionally, the loop 8 conformations differ in mutants, as compared to the wt -enzyme. Because part of loop 8, including position 263, is located at the enzyme dimer interface, altered loop 8 conformations modulate the relative orientation of both monomers. In the case of the mutants characterized in this study, the dimer reorientation yields significant dis- placement up to 5.2 (e Å .g., α sA1, between equivalent carbon α positions), as compared to wt-SsoPox (Figure S5). 11,42 Similar reorientations were previously observed upon substrate binding or mutations of W263 . Although the active sites of the improved variant structures superpose well onto the wt-structure, all selected mutants display altered loop 8 conformations. Changes are subtle for some of the variants (e.g. αsA1), but are large for others (e.g. αsC6) (Fig. 4). Notably, mutant αsA6 was crystallized in two different conformations: one in which loop 8 adopts a wt-like conformation, referred to as closed conformation (CC), and another one where loop 8 is unfolded, referred to as open conformation (OC) (Figure S6). In the two other mutants, αsD6 and αsB5, loop 8 could not be modeled due to the lack of electronic density, which was likely to be due to the high level of motion of this enzyme region. Analysis of the normalized thermal motion B-factor supports this hypothesis. It also confirms that loop 8 is highly mobile in all selected mutants, as compared to wt-SsoPox (Fig. 5). The higher mobility of loop 8 may also partly explain the improved ability of the variants to hydrolyze phosphotriesters with large substituents. Improved mutant αsD6 lost the thiono-effect. SsoPox presents a marked preference for oxono-OP substrates as compared to thiono-ones (>100-fold in catalytic efficiency) . Conversely, PTEs do not exhibit such 27,51 a drastic preference, paraoxon being a slightly better substrate than parathion for the enzyme . Interestingly, Scientific REPO R TS | 7: 16745 | DOI:10.1038/s41598-017-16841-0 5 www.nature.com/scientificreports/ Enhancement/ −1 −1 −1 Substrate SsoPox k (s ) K (µM) K (µM) k /K (M .s ) wt cat M I cat M 1 4 2 wt° 1.26 ± 0.13 × 10 2.43 ± 0.37 × 10 — 5.12 ± 1.31 × 10 1 1 2 4 αsA1 2.70 ± 0.23 × 10 8.00 ± 1.55 × 10 — 3.37 ± 0.71 × 10 65 3 3 αsA6 5.87 ± 0.47 1.63 ± 0.63 × 10 — 3.61 ± 1.43 × 10 7 1 2 4 αsB5 2.06 ± 0.14 × 10 6.01 ± 0.96 × 10 — 3.43 ± 1.46 × 10 67 1 3 4 αsC6 4.27 ± 0.25 × 10 1.50 ± 0.16 × 10 — 2.86 ± 0.35 × 10 55 1 3 4 Ethyl-paraoxon (I) αsD6 3.22 ± 0.18 × 10 1.09 ± 0.12 × 10 — 2.95 ± 1.55 × 10 58 αsB5 W263M * * — 5.62 ± 0.11 × 10 11 αsB5 W263L * * — 8.07 ± 0.28 × 10 16 1 3 3 αsB5 W263I 4.09 ± 0.89 × 10 6.64 ± 1.75 × 10 — 6.16 ± 5.07 × 10 12 + 3 W263L * * — 2.37 ± 0.33 × 10 5 + 2 3 W263M 6.82 ± 0.57 9.31 ± 1.63 × 10 — 7.33 ± 1.42 × 10 14 wt° ND ND ND ND ND αsA1 ND ND ND ND ND −2 2 2 αsA6 3.20 ± 0.20 × 10 1.34 ± 0.18 × 10 2.39 ± 0.35 × 10 ND 250 ± 1020 −1 2 2 αsB5 2.69 ± 0.15 × 10 2.72 ± 0.49 × 10 — 9.86 ± 3.14 × 10 ND −3 1 2 αsC6 8.53 ± 0.30 × 10 7.70 ± 1.10 × 10 — 1.10 ± 0.16 × 10 ND Ethyl-parathion (II) −1 2 3 αsD6 7.35 ± 0.36 × 10 5.34 ± 0.62 × 10 — 1.38 ± 0.57 × 10 ND −2 2 2 αsB5 W263M 5.24 ± 0.23 × 10 2.02 ± 0.31 × 10 — 2.59 ± 0.74 × 10 ND −2 2 2 αsB5 W263L 4.35 ± 0.16 × 10 3.46 ± 0.37 × 10 — 1.26 ± 0.44 × 10 ND −2 2 2 αsB5 W263I 4.74 ± 0.21 × 10 2.61 ± 0.37 × 10 — 1.82 ± 0.57 × 10 ND W263L ND ND — ND ND W263M ND ND — ND ND 3 3 wt° 2.71 ± 0.64 2.14 ± 0.68 × 10 — 1.27 ± 0.70 × 10 1 1 3 4 αsA1 4.35 ± 0.87 × 10 1.90 ± 0.52 × 10 — 2.29 ± 0.78 × 10 18 2 4 αsA6 5.89 ± 0.55 5.46 ± 1.02 × 10 — 1.08 ± 0.23 × 10 8 Methyl-paraoxon (III) αsB5 * * — 4.31 ± 0.14 × 10 3 1 2 4 αsC6 2.36 ± 0.37 × 10 7.59 ± 2.04 × 10 — 3.11 ± 0.97 × 10 25 1 3 4 αsD6 4.25 ± 0.52 × 10 2.08 ± 0.34 × 10 — 2.04 ± 0.42 × 10 16 −3 2 wt° 1.10 ± 0.02 × 10 1.21 ± 0.10 × 10 — 9.09 ± 0.90 1 −2 2 1 αsA1 1.30 ± 0.50 × 10 3.53 ± 0.34 × 10 — 3.68 ± 1.46 × 10 4 −2 2 1 αsA6 1.55 ± 0.15 × 10 2.54 ± 0.39 × 10 1 520 ± 384 6.10 ± 1.11 × 10 7 Methyl-parathion (IV) −1 2 2 αsB5 1.50 ± 0.04 × 10 1.58 ± 0.15 × 10 — 9.49 ± 0.94 × 10 104 −3 2 1 αsC6 3.00 ± 0.10 × 10 1.21 ± 0.15 × 10 — 2.48 ± 0.32 × 10 3 2 4 αsD6 6.89 ± 0.35 3.43 ± 0.44 × 10 — 2.01 ± 0.28 × 10 2210 −4 2 wt° 8.90 ± 0.40 × 10 1.60 ± 0.29 × 10 — 5.56 ± 1.26 1 −4 2 αsA1 7.10 ± 0.40 × 10 2.22 ± 0.35 × 10 — 3.20 ± 0.54 0.5 αsA6 ND ND ND ND ND −2 2 1 αsB5 1.84 ± 0.16 × 10 9.13 ± 1.61 × 10 — 2.02 ± 1.01 × 10 4 −2 1 2 αsC6 2.47 ± 0.44 × 10 3.20 ± 1.90 × 10 820 ± 343 7.72 ± 4.78 × 10 139 −1 2 2 Malathion (V) αsD6 1.00 ± 0.02 × 10 1.13 ± 0.10 × 10 — 9.08 ± 2.00 × 10 163 −3 2 1 αsB5 W263M 8.92 ± 0.67 × 10 8.61 ± 1.32 × 10 — 1.04 ± 0.51 × 10 2 −2 2 1 αsB5 W263L 1.11 ± 0.08 × 10 9.28 ± 1.32 × 10 — 1.20 ± 0.60 × 10 2 −2 2 1 αsB5 W263I 1.18 ± 0.08 × 10 7.80 ± 1.08 × 10 — 1.51 ± 0.71 × 10 3 W263L ND ND — ND ND W263M ND ND — ND ND wt ND ND — ND ND αsD6 * * — 3.12 ± 0.09 × 10 ND αsB5 * * — 2.07 ± 0.05 × 10 ND αsB5 W263M * * — 1.64 ± 0.03 × 10 ND Chlorpyrifos (VI) αsB5 W263L * * — 1.46 ± 0.03 × 10 ND αsB5 W263I * * — 1.42 ± 0.03 × 10 ND W263L ND ND — ND ND W263M ND ND — ND ND Continued Scientific REPO R TS | 7: 16745 | DOI:10.1038/s41598-017-16841-0 6 www.nature.com/scientificreports/ Enhancement/ −1 −1 −1 Substrate SsoPox k (s ) K (µM) K (µM) k /K (M .s ) wt cat M I cat M wt ND ND — ND ND 3 3 αsD6 8.93 ± 0.79 2.30 ± 0.31 × 10 — 3.89 ± 2.56 × 10 ND 2 3 αsB5 1.26 ± 0.06 4.10 ± 0.65 × 10 — 3.03 ± 0.98 × 10 ND −1 2 2 αsB5 W263M 3.89 ± 0.38 × 10 3.91 ± 1.04 × 10 — 9.92 ± 3.65 × 10 ND Diazinon (VII) −1 2 3 αsB5 W263L 5.61 ± 0.23 × 10 4.74 ± 0.53 × 10 — 1.18 ± 0.44 × 10 ND −1 2 3 αsB5 W263I 4.01 ± 0.30 × 10 3.24 ± 0.72 × 10 — 1.24 ± 0.42 × 10 ND W263L ND ND — ND ND W263M ND ND — ND ND wt ND ND — ND ND −1 2 3 αsD6 7.94 ± 0.50 × 10 5.28 ± 0.80 × 10 — 1.50 ± 0.62 × 10 ND −2 2 1 αsB5 1.40 ± 0.20 × 10 3.14 ± 1.43 × 10 — 9.11 ± 1.40 × 10 ND −2 2 1 αsB5 W263M 2.22 ± 0.12 × 10 3.31 ± 0.53 × 10 — 6.69 ± 2.31 × 10 ND Fenitrothion (VIII) −2 2 1 αsB5 W263L 2.63 ± 0.13 × 10 6.54 ± 0.75 × 10 — 4.01 ± 1.79 × 10 ND −2 2 1 αsB5 W263I 2.01 ± 0.08 × 10 3.97 ± 0.46 × 10 — 5.05 ± 1.81 × 10 ND −3 2 W263L 3.66 ± 0.15 × 10 4.15 ± 0.44 × 10 — 8.82 ± 3.32 ND −4 2 W263M 9.49 ± 0.43 × 10 1.57 ± 0.28 × 10 — 6.06 ± 1.55 ND wt ND ND — ND ND −2 2 1 αsB5 4.28 ± 0.39 × 10 9.41 ± 1.72 × 10 — 4.55 ± 2.30 × 10 ND −1 3 1 αsD6 4.06 ± 1.60 × 10 8.60 ± 3.97 × 10 — 4.72 ± 4.04 × 10 ND αsB5 W263M * * — 1.12 ± 0.46 × 10 ND Fensulfothion (IX) −2 3 1 αsB5 W263L 4.94 ± 1.35 × 10 2.57 ± 1.04 × 10 — 1.92 ± 1.30 × 10 ND αsB5 W263I * * — 6.50 ± 0.07 ND W263L ND ND — ND ND W263M ND ND — ND ND wt ND ND — ND ND αsD6 * * — 1.64 ± 0.02 × 10 ND αsB5 * * — 7.82 ± 0.14 × 10 ND αsB5 W263M * * — 1.15 ± 0.02 × 10 ND Coumaphos (X) αsB5 W263L * * — 1.16 ± 0.05 × 10 ND αsB5 W263I * * — 1.11 ± 0.04 × 10 ND W263L * * — 6.01 ± 0.18 ND W263M * * — 4.80 ± 0.32 ND wt° * * — 3.16 × 10 1 αsA1 * * — 1.41 × 10 0.004 αsA6 * * — 1.78 × 10 0.06 3-oxo-C10-AHL (l) (XI) αsB5 ND ND — ND ND αsC6 * * — 5.84 × 10 0.02 αsD6 ND ND — ND ND wt° * * — 2.36 × 10 1 αsA1 * * — 1.63 × 10 0.7 αsA6 * * — 4.40 × 10 186 Undecanoic-γ-lactone (XII) αsB5 * * — 9.05 × 10 0.4 αsC6 * * — 5.73 × 10 243 αsD6 * * — 1.60 × 10 7 wt° * * — 7.86 × 10 1 αsA1 * * — 2.55 × 10 0.03 αsA6 * * — 9.65 × 10 123 Undecanoic-δ-lactone (XIII) αsB5 * * — 5.33 × 10 0.08 αsC6 * * — 3.60 × 10 5 αsD6 * * — 7.50 × 10 0.9 Table 2. Kinetic characterization of SsoPox variants. *Corresponds to undetermined values because enzyme saturation could not be reached. ND corresponds to undetectable values because substrates were not hydrolyzed 46 + 11 or too poorly hydrolyzed to be observable. °Data taken from Hiblot et al. . Data taken from Hiblot et al. . some selected variants retained a strong preference for oxono-phosphotriesters (~10-fold), whereas others such as αsD6 lost the thiono-effect (Figure S7). The wild-type enzyme exhibits a low K value with methyl-parathion −3 −1 (121 µ M), but also a very low k value (1.10 × 10 s ) (Table 2). This may indicate a stronger interaction cat Scientific REPO R TS | 7: 16745 | DOI:10.1038/s41598-017-16841-0 7 www.nature.com/scientificreports/ between the thiono moiety and the metal cations (low K value), which would subsequently impair catalysis (low k values) . Interestingly, the improved variants do not exhibit significant difference in K values for cat M thiono-phosphotriesters, but show a dramatic increase in k values, as compared to wild-type enzymes (e.g. cat −1 −3 −1 k value of 6.89 s and 1.10 × 10 s , as compared to methyl-parathion for the αsD6 variant and wild-type cat enzymes, respectively). The two selected variants, namely αsB5 and αsD6, which lost the thiono-effect have only two mutations in common: V27A and Y97W. Both residues are located relatively close to the bi-metallic active site, yet do not interact directly with the metals (V27-Fe: 6.7 Å, Y97-Co: 4.3 Å; Fig.  6). Because other PLLs have 33,53 been shown to exhibit high charge coupling between the β metal and the conserved tyrosine residue , Y97 is a likely candidate for the thiono-effect, and will be examined in future studies. Mutations of W263 and αsB5 are incompatible. αsB5 (V27A-I76T-Y97W-Y99F-L130P-L226V) is the only selected mutant that does not contain a substitution of position W263. Because W263 substitutions were previ- ously shown to be key for increasing SsoPox’s phosphotriesterase activity, we used saturation mutagenesis to introduce W263 mutations to the background of αsB5 . The best selected mutants against paraoxon (Figure S2C) were found to be αsB5-W263I/L/M. Surprisingly, none of these mutants demonstrated an increase in phosphotriesterase activity in comparison with αsB5, but rather a decrease in activity (Table 2). This is intriguing because, taken individually, these mutations have been shown to greatly improve phosphotriesterase activity. For example, with ethyl-paraoxon as a substrate, W263M and αsB5 produce a catalytic efficiency improvement of 14 and 67-fold, respectively, while the combination αsB5-W263M increases activity only 11-fold (Table 2). This example of negative epistasis could be due to the mode of action of these mutations. Structures reveal that αsB5 mutations have a destabilizing effect on loop 8, as illustrated by the B factor analysis. The same was observed for W263M in a previous study . Combining these two dest- abilizing sets of mutations may have harmed the active site integrity, as previously described with another enzyme as “conformational active site disorder” , and in particular the necessary alignment and pre-organization of the catalytic residues. The destabilizing effect of these mutations can also be observed on the overall enzyme stability, with α sB5, W263M, αsB5-W263M exhibiting T values of 70.4 °C, 85.3 °C and 69.1 °C, respectively, as compared to 106 °C for the wt enzyme (Table S4). Here, the effects of the combined 5 mutations harbored by variant α sB5 together with the previ- ously identified positive substitutions at residue W263 were not additive underlining that single beneficial mutations do not necessarily induce positive epistatic effects. Synergistic effects of mutations are oen ft complex to predict and most engineering strategies usually lead to an optimization plateau difficult to overcome . Improvement of engineered variants over W263 mutation. In a previous report, we noted that muta- tion of W263 into all the other 19 residues lead to an increase of the promiscuous paraoxonase activity . Here we note that most of the improved variants (i.e. αsA1, αsA6, αsB5, αsC6 & αsD6) (i) exhibit higher paraoxonase activity than W263L/M, and (ii) foremost are capable of hydrolyzing phosphotriesters that are not substrates for the wild-type enzyme or W263 mutants, such as ethyl-parathion, chlorpyrifos or fensulfothion (Table 2). The activity spectrum of these engineered variants against the tested phosphotriesters is therefore much wider than the wild-type enzyme and W263 mutants. Relation between active site loop disorder and broad enzymatic specificity. The organophospho - rus compound bioremediation potential of several mutants was investigated. In particular, the best variants, αsD6 and αsB5, alongside other improved variants αsB5-W263I/L/M, W263L/M and the wt-enzyme, were character- ized for their phosphotriesterase activity with ethyl-paraoxon (I), ethyl-parathion (II), malathion (V), chlorpy- rifos (VI), diazinon (VII), fenitrothion (VIII), fensulfothion (IX) and coumaphos (X) (Fig. 2). Interestingly, we found that most mutants were able to hydrolyze fensulfothion, previously reported as an inhibitor of the wild-type enzyme . Fensulfothion was indeed previously shown to bind head-to-tail into the wild-type enzyme active side , a non-productive binding mode. Changes in the active site loop 8 conformational ensemble of engineered mutants might have allowed for the productive binding of fensulfothion. All mutants selected and used in this study are destabilized as compared to the wt-enzyme (Fig. 7; Table S4). Structures reveal that these mutations, located on loop 8, contribute to increasing the mobility of the loop. Interestingly, the ability of the tested mutants to hydrolyze a wide range of phosphotriesters, including com- pounds that were inhibitors for the wt-enzyme, correlates with a lower T (Fig. 7). This correlation suggest that the active site loop 8 degree of disorder may modulate the enzyme’s ability to bind and hydrolyze a variety of phosphotriesters, a wider conformational ensemble of loop 8 being correlated with a broader enzyme specificity. This observation is consistent with previous studies highlighting the importance of the destabilization of active 54,56 site loops to evolve new activities and new substrate specificity . Such flexible active site loops involved in enzymatic specificity supports the notion of fold polarity: a portion of the active site (e.g. the loop) is weakly connected to the enzyme scao ff ld and thereby makes it possible for new functions to evolve with few changes . Mutant αsD6 is an active, broad spectrum phosphotriester biodecontaminant. Mutant αsD6 was further investigated for bioremediation considerations. Organophosphate-based pesticides are usually spread with concentrations in the millimolar range leading to micromolar-ranged contaminations. We thus investigated the potential of αsD6 for decontaminating pesticide solutions at 250 µM to simulate a contamination of ground- −2 −3 water or runoff waters. We investigated two [E]/[S] ratios, 10 and 10 respectively, for decontaminating OP solutions and determined the time required to hydrolyze 95% of pesticide preparations (Fig. 3B; Figure S4). When −3 considering the [E]/[S] ratio of 10 , 95% degradation was achieved for four substrates (ethyl-paraoxon, diazinon, fenitrothion and coumaphos) within 10 minutes. Two substrates were decontaminated within an hour (ethyl-parathion and chlorpyrifos) and fensulfothion was degraded in three hours. Scientific REPO R TS | 7: 16745 | DOI:10.1038/s41598-017-16841-0 8 www.nature.com/scientificreports/ Figure 4. Loop 8 conformations in improved mutants. Active site loop 8 region of mutant structures superimposed with the wt-SsoPox structure. The active site location is indicated by the presence of the bi- metallic catalytic center shown as spheres. The wt-SsoPox structure loop 8 conformation is shown as dark green cartoon. (A) αsA1′s loop 8 is shown in orange cartoon, (B) αsD6′s loop 8 is shown in green cartoon, and is not entirely visible in the electronic density maps, (C) αsC6′s loop 8 is shown in blue cartoon, (D) αsB5′s loop 8 is shown in yellow cartoon, and is not entirely visible in the electronic density maps, (E) αsA6-CC’s loop 8 is shown in magenta cartoon, (F) αsA6-OC’s loop 8 is shown in cyan cartoon. −2 When increasing enzyme concentration 10-fold ([E]/[S] ratio of 10 ), five substrates (ethyl-paraoxon, ethyl-parathion, diazinon, fenitrothion and coumaphos) were decontaminated within two minutes, while fensul- fothion needed ~20 minutes and the two other tested (chlorpyrifos and malathion) required 40 minutes. Given Scientific REPO R TS | 7: 16745 | DOI:10.1038/s41598-017-16841-0 9 www.nature.com/scientificreports/ Figure 5. Mutants’ active site loop 8 exhibit higher conformational flexibility. (A) Cartoon representation of the superposition of the structures (monomers) of wt-Ssopox (dark green), αsA1 (orange), αsB5 (yellow), αsD6 (green), αsC6 (blue), αsA6-CC (magenta) and αsA6-OC (cyan). (B) Smoothed ribbon representation of the superposition of the structures (monomers) of of wt-Ssopox (green), αsA1 (orange), αsB5 (yellow), αsD6 (green), αsC6 (blue), αsA6-CC (magenta) and αsA6-OC (cyan). (C) Positional distributions of normalized B-factor values (the x-axis represents residue number) for wt-SsoPox (black line), αsA1 (orange), αsB5 (yellow), αsD6 (green), αsC6 (blue), αsA6-OC (blue line) and (magenta line). A zoomed inset highlights the loop 7 and 8 sequence region. both its catalytic proficiency and high thermal stability (T = 82.5 °C), αsD6 is a promising candidate for organ- ophosphorus compound bioremediation and, in particular, the decontamination of water runoffs, soils, food products and materials. Moreover, SsoPox variants were previously shown to be compatible with immobilization 48,58 steps including alginate beads and polyurethane-based coatings , and could be of prime interest for the devel- opment of filtration devices for water treatment purposes. Methods Structure-based identification of mutations. Sharing the same (α/β) barrel topology, SsoPox and BdPTE structures (Pdb ID 2vc5 and 1dpm, respectively) were superimposed (Fig. 1) using PyMol , making it possible to identify structurally equivalent residues at their respective active sites. In an effort to reshape the SsoPox active site, side chains that were well superimposed between the two enzymes were mutated using Coot into the residue present in BdPTE , with the exception of V27 which was mutated in Ala (and not Gly) to avoid increased entropy (Fig. 1; Table S1). However, due to the major differences in loop 8 and 7 between PLLs and PTEs, numerous residues were not superimposable. These active site residues (Y99, L228, F229 and W263) were therefore mutated in an effort to mimic the Bd PTE active site cavity in terms of shape and chemical nature, as Scientific REPO R TS | 7: 16745 | DOI:10.1038/s41598-017-16841-0 10 www.nature.com/scientificreports/ Figure 6. Localization of mutations harbored by SsoPox monovariants, αsA1, αsA6, αsB5, αsC6 and αsD6. Mutations are highlighted in red spheres using the SsoPox wild-type (PDB 3UF9) as a template. Figure 7. Overview of SsoPox wild-type and variants properties. Catalytic efficiencies towards organophosphate-based pesticides are presented using a color code ranging from white to dark blue. Melting temperatures are presented using a color code ranging from red to green. illustrated by the chimeric reconstruction of SsoPox carrying all these mutations (Figure S1). Others mutations were also implemented in the dataset, e.g. C258A and C258L and different substitutions of the key position W263 (F/L/M/A) that have been shown to improve the phosphodiesterase activity . Synthesis of the combinatorial library. The SsoPox coding gene was amplified from the previously described pET22b-SsoPox plasmid . For each of the 14 mutations, a mutagenesis primer of 30–33 bp (~10–15 bp from each mutation side) was synthesized (Table S2). For close mutations, primers can share several mutations and different primers were synthesized and mixed to keep an equivalent statistical probability for each of the 14 mutations. Typically, 2 pmol of primer mixture and 100 ng of DNaseI (TaKaRa)-generated fragments of the SsoPox gene were assembled as previously described , followed by nested PCR with external cloning primers (SsoPox-lib-pET-5′/3′) (Table S2). The PCR product was then cloned into a customized version of the pET - 32b(-ΔTrx) plasmid and then electroporated into E. cloni cells (Lucigen, USA). Aer a ft gar-plate growth, plasmid extraction was performed to create the plasmid bank. Variant genes were PCR amplified using T7-prom and Scientific REPO R TS | 7: 16745 | DOI:10.1038/s41598-017-16841-0 11 www.nature.com/scientificreports/ pET-RP primers (Table S2) and sequenced. The combinatorial library shows an average of 5.3 ± 3.3 library muta- tions by sequencing 10 randomly picked colonies, and 0.43 ± 0.53 random mutations per gene. Screening of the library. The plasmid library was used to transform the Escherichia coli strain BL21(DE )-pGro7/GroEL (TaKaRa) to obtain colonies expressing a library of SsoPox variants. Randomly picked clones (184), representing a coverage of 3.8% of the library, were grown on a 96-well plate in 500 µ L of ZYP medium as previously described . Production of proteins and chaperones was induced after five hours of cul- ture at 37 °C by reducing the temperature to 25 °C, adding CoCl (0.2 mM) and adding arabinose (0.2%, w/v). After overnight growth, cell lysate was used to perform activity screening with 100 µ M of paraoxon substrate 35,49 (Figure S2B) after partial purification of the protein with a heating step of 15 minutes incubation at 70 °C . The screening was performed in 50 mM HEPES pH 8, 150 mM NaCl, 0.2 mM CoCl . Kinetics of paraoxon/ methyl-parathion hydrolysis were monitored by following absorbance at 405 nm for 10 minutes using a microplate reader (Synergy HT, BioTek, USA) and the Gen5.1 software. Saturation library. Saturation mutagenesis directed on residue 263 of the αsB5 variant encoding gene, was per- formed using NNS degenerated primers (W263NNS direct: 5′-TGCACCATTGATNNSGGCACCGCAAAACCG -3′ and W263NNS reverse: 5′-CGGTTTTGCGGTGCCSNNATCAATGGTGCA-3′) and pET32b-ΔTrx carry- ing αsB5 coding gene as template. PCR amplification was performed with 2.5 U of PfuTurbo DNA polymerase (Agilent) according to the manufacturer’s recommendations [95 °C, 5 minutes; 20x (95 °C, 30 seconds; 55 °C, one minute; 68 °C, 15 minutes); 68 °C, 25 minutes]. DNA was digested with DpnI to remove the methylated parental template. E. coli BL21(DE )-pGro7/GroEL (TaKaRa) competent cells were transformed with the plasmid pool and −1 −1 plated on LB-agar supplemented with 100 µ g.mL ampicillin and 34 µ g.mL chloramphenicol. Over a theoret- −1 ical diversity of 20 sequences, 181 variants were collected in a microplate containing LB (100 µg.mL ampicillin −1 and 34 µ g.mL chloramphenicol) for characterization. This oversampling by a factor of nearly 10-fold ensure a >99% probability to get the 20 possible substitutions. Screening was performed as described above (Figure S2C). Clones were produced in liquid ZYP medium and their activity was screened with 100 µM of paraoxon substrate. The plasmids corresponding to the most interest- ing variants were extracted and the SsoPox variant encoding genes were sequenced. Production-purification of SsoPox wild-type and variants. The genes coding for the SsoPox wild-type and its variant were cloned in a pET32b-Δtrx plasmid. Productions were performed using the E. coli BL21(DE )-pGro7/GroEL (TaKaRa) chaperone expressing strain. Starter cultures were produced in an −1 −1 auto-inducible ZYP medium (supplemented with 100 µ g.mL ampicillin and 34 µ g.mL chloramphenicol). When OD reached a value of 0.8–1, CoCl was added (final concentration 0.2 mM) as well as L-arabinose 600nm 2 −1 (final concentration 2 g.L ) to induce the production of chaperones GroEL/ES and the temperature was decreased to 23 °C for 16–20 hours. Cells were harvested by centrifugation (4400 g, 4 °C, 20 minutes) and resus- −1 pended in lysis buffer (50 mM HEPES pH 8.0, 150 mM NaCl, 0.2 mM CoCl , 0.25 mg.mL lysozyme, 0.1 mM −1 PMSF and 10 µ g.mL DNAseI) and were stored at -80 °C. Cells were thawed and lysed in three steps of 30 sec- onds of sonication (Qsonica, Q700; Amplitude 45). Cell debris was removed by centrifugation (20000 g, 4 °C, 15 minutes). As SsoPox and its variants are hyperthermostable, a pre-purification step was performed by heat- ing the lysate for 30 minutes at 80 °C. Precipitated host proteins were removed by centrifugation (20000 g, 4 °C, 15 minutes). SsoPox and its variants were collected by ammonium sulfate precipitation (75%) and resuspended in 8 ml of activity bue ff r (50 mM HEPES pH 8.0, 150 mM NaCl, 0.2 mM CoCl ). The remaining ammonium sulfate was eliminated by injection on a desalting column (HiPrep 26/10 desalting, GE Healthcare; ÄKTA Avant) and TM concentrated to 2 mL for separation on exclusion size chromatography (HiLoad 16/600 Superdex 75 pg, GE Healthcare; ÄKTA Avant). Final purity was monitored by SDS-PAGE and protein concentration was measured with a NanoDrop 2000 spectrophotometer (Thermo Scientific). Kinetic assays. Generalities. Catalytic parameters were evaluated in triplicate at 25 °C, and recorded using a microplate reader (Synergy HT, BioTek, USA) and the Gen5.1 software, in a 6.2 mm path length cell for 200 µ L reaction in 96-well plate, as previously explained . Catalytic parameters were obtained by fitting the data to the Michaelis-Menten (MM) equation using Graph-Pad Prism 6 software. When V could not be reached in the max experiments, catalytic efficiency was obtained by fitting the linear part of MM plot to a linear regression using Graph-Pad Prism 6 software. For some Sso Pox variants, the MM plot was fitted to the substrate inhibition equa - tion using Graph-Pad Prism 6 software enabling us to determine a K for several substrates. Consequently, the calculated catalytic efficiencies in these conditions are true only at low substrate concentrations. In some other cases, saturation could not be reached, therefore k /K values were determined using linear regression. For catM Coumaphos, data were fitted to one-phase decay non-linear regression. Phosphotriesterase activity characterization. e k Th inetic assays were carried out in activity buer ff . The conditions used for pesticides were as follows: ethyl-paraoxon (I), ethyl-parathion (II), methyl-paraoxon (III) and −1 −1 methyl-parathion (IV) hydrolysis was recorded at 405 nm (ε = 17000 M .cm ). The concentration interval was between 0.05 and 2 mM from 100 mM stock solution in ethanol . Malathion (V) hydrolysis was monitored at −1 −1 46 412 nm, (ε = 13400 M .cm ) . The concentration interval was between 0.05 and 2 mM from 100 mM stock solution in DMSO. DTNB was added to the activity buffer. Chlorpyrifos (VI) hydrolysis was recorded at 310 nm −1 −1 51 (ε = 5562 M .cm ) . The concentration interval was between 0.05 and 2 mM from 100 mM stock solution in −1 −1 27 1-propanol. Diazinon (VII) hydrolysis was recorded at 228 nm (ε = 3300 M .cm ) . The concentration interval was between 0.05 and 2 mM from 200 mM stock solution in ethanol. Fenitrothion (VIII) hydrolysis was recorded −1 −1 52 at 358 nm (ε = 18700 M .cm ) . The concentration interval was between 0.05 and 2 mM from 200 mM stock Scientific REPO R TS | 7: 16745 | DOI:10.1038/s41598-017-16841-0 12 www.nature.com/scientificreports/ −1 −1 27 solution in ethanol. Fensulfothion (IX) hydrolysis was recorded at 284 nm (ε = 8000 M .cm ) . The concentra- tion interval was between 0.05 and 2 mM from 200 mM stock solution in methanol. Coumaphos (X) hydrolysis was monitored by fluorescence of released chlorferon (excitation 360/40 and emission 460/40) . A linear corre- lation between fluorescence and chlorferon was found (Figure S3). Measurements were taken at 5 µ M from a 200 mM stock solution in DMSO. We assumed that , thus enabling us to estimate . KS [] k /K M catM Lactonase activity characterization. Lactonase kinetics were performed using a previously described protocol . e Th time course hydrolysis of lactones were performed in lac buffer (2.5 mM Bicine pH 8.3, 150 mM NaCl, 0.2 mM CoCl , 0.25 mM Cresol purple and 0.5% DMSO) over a concentration range 0–2 mM for 3-oxo-C10 AHLs (XI) or 0–5 mM for undecanoic- γ/δ -lactones (XII, XIII). Cresol purple (pK 8.3 at 25 °C) is a pH indicator used to fol- low lactone ring hydrolysis by acidification of the medium. Molar coefficient extinction at 577 nm was evaluated recording absorbance of the buffer over an acetic acid range of concentration 0–0.35 mM. Bioremediation of pesticide solutions at 250 µM. For all OP substrates but malathion, experiments were performed in triplicate at 25 °C using a microplate reader (Synergy HT, BioTek, USA). Wavelengths were chosen as described in the kinetic assays section. Degradation of malathion was followed by GC/MS in triplicate. OP concentration was fixed at 250 µ M. All the substrates were soluble at this concentration. The best variants of SsoPox for each substrate were used, namely αsD6 for all the OPs. Two enzyme-to-substrate ratios were [E]/[S] −2 −3 −1 −1 used, 10 and 10 , corresponding to enzyme concentrations of 90 µg.ml and 9 µg.ml respectively. The reac- tions were monitored until the plateau was reached. Experimental measures were obtained using Gen5.1 sowa ft re, then analyzed with GraphPad Prism 6 software. Curves were then fitted using One-Phase Decay non-linear regression with the equation (1): () −kt YY =− (0 Plateaue ) ×+ Plateau (1) where Y0 = 0% and Plateau = 100%. The rate constant k was determined and the time required to observe a 95% decontamination was calculated accordingly. The curves and results of the fits are shown in Figure S4. −3 Degradation at 95% was confirmed using a final point by GC/MS for all pesticides with a 10 enzyme-to-substrate ratio. 100 µ L of activity buer ff solution was first extracted with 100 µ L chloroform. Organic extracts were analyzed by using a Clarus 500 gas chromatograph equipped with a SQ8S MS detector (Perkin −1 Elmer, Courtaboeuf, France). 1 µ L of organic extract was volatilized at 220 °C (split 15 mL.min ) in a deactivated FocusLiner with quartz wool (SGE, Ringwood, Australia) and compounds separated on an Elite-5MS column (30 m, −1 0.25 mm i.d., 0.25 mm film thickness) for 12 minutes using a temperature gradient (80–280 °C at 30 °C.min , five −1 minutes’ hold). Helium flowing at 2 mL.min was used as the carrier gas. The MS inlet line was set at 280 °C and the electron ionization source at 280 °C and 70 eV. Full scan monitoring was performed from 40 to 400 m/z in order to identify chemicals by spectral database search using MS Search 2.0 operated with the Standard Reference Database 1 A (National Institute of Standards and Technology, Gaithersburg, MD, USA). m/z is the mass-to-charge ratio of the base peak fragment detected for each molecule. In the case of weak pesticide signals, extracted ion chromato- grams were generated with base peak ions to confirm the presence of chemicals (Paraoxon (I ) and Parathion ethyl (II) m/z 109; Malathion (V) and Fenitrothion (VIII) m/z 125; Diazinon (VII) m/z 137; Coumaphos (X) m/z 97; Fensulfothion (IX) m/z 293; and Chlorpyrifos (VI) m/z 197). Selected Ion Recording using base peaks ions was applied in order to specifically monitor pesticides and collect peak areas for kinetics. All samples were analyzed over short periods of time to avoid signal drift. All data were processed using Turbomass 6.1 (Perkin Elmer). 42,46 Crystallization. Crystallization assays were performed as previously described . Crystallization was performed using the hanging drop vapor diffusion method in 96-well plates (Greiner Microplate, 96 well, PS, F-bottom) on ViewDrop II seals (TPP Labtech). Equal volumes (0.5 µ L) of protein and reservoir solutions were mixed, and the resulting drops were equilibrated against a 150 µ L reservoir solution containing 20–30% (w/v) PEG 8000 and 50 mM Tris-HCl buffer (pH 8). Crystals appeared aer f ft ew days at 4 °C. Data collection and structure refinement. Crystals were first transferred to a cryoprotectant solution con- sisting of the reservoir solution and 20% (v/v) glycerol. Crystals were then flash-cooled in liquid nitrogen. X-ray diffraction data were collected at 100 K using synchrotron radiation at the ID23-1 beam line (ESRF, Grenoble, France) and using an ADSC Q315r detector. X-ray diffraction data were integrated and scaled with the XDS package (Table 1) . The phases were obtained using the native structure of SsoPox (PDB code 2vc5) as a starting model, 62,63 performing a molecular replacement with MOLREP or PHASER . The models were built with Coot and refined 60,64 using REFMAC . We note that three structures presented in this work exhibit one disordered monomer (and low corresponding electronic density): αsB5 (monomer D, total of 4 monomers), αsD6 (monomer C, out of a total of 4 monomers) and αsC6 (monomer C, total of 4 monomers). Structure illustrations were performed using PyMOL . Relative B-factor analysis. e o Th ccupancies of all residues in all tested structures were set to 1 for this anal- ysis. For residues with alternate conformations, the sums of occupancies were set to 1. Structures were re-refined 60,64 with REFMAC . The relative B-factor values were obtained by normalizing the B-factor values of each residue 11,54,57 by the average B-factor of the whole structure as previously described . References 1. Wolfenden, R. Benchmark Reaction Rates, the Stability of Biological Molecules in Water, and the Evolution of Catalytic Power in Enzymes. Annu. Rev. Biochem. 80, 645–667 (2011). Scientific REPO R TS | 7: 16745 | DOI:10.1038/s41598-017-16841-0 13 www.nature.com/scientificreports/ 2. Singh, R. K., Tiwari, M. K., Singh, R. & Lee, J.-K. From Protein Engineering to Immobilization: Promising Strategies for the Upgrade of Industrial Enzymes. Int. J. Mol. Sci. 14, 1232–1277 (2013). 3. Bar-Even, A. et al. The Moderately Efficient Enzyme: Evolutionary and Physicochemical Trends Shaping Enzyme Parameters. Biochemistry (Mosc.) 50, 4402–4410 (2011). 4. Tokuriki, N. et al. Diminishing returns and tradeoffs constrain the laboratory optimization of an enzyme. Nat. Commun. 3, 1257 (2012). 5. Burton, S. G., Cowan, D. A. & Woodley, J. M. The search for the ideal biocatalyst. Nat. Biotechnol. 20, 37–45 (2002). 6. Gupta, R. D. & Tawfik, D. S. Directed enzyme evolution via small and effective neutral drift libraries. Nat. Methods 5, 939–942 (2008). 7. Bommarius, A. S. Biocatalysis: A Status Report. Annu. Rev. Chem. Biomol. Eng. 6, 319–345 (2015). 8. Currin, A., Swainston, N., Day, P. J. & Kell, D. B. Synthetic biology for the directed evolution of protein biocatalysts: navigating sequence space intelligently. Chem Soc Rev 44, 1172–1239 (2015). 9. Reetz, M. T. Directed Evolution of Selective Enzymes: Catalysts for Organic Chemistry and Biotechnology. (Wiley-VCH Verlag GmbH & Co. KGaA, 2016), https://doi.org/10.1002/9783527655465. 10. Directed Enzyme Evolution: Advances and Applications. (Springer International Publishing, 2017). https://doi.org/10.1007/978-3- 319-50413-1. 11. Hiblot, J., Gotthard, G., Elias, M. & Chabriere, E. Differential Active Site Loop Conformations Mediate Promiscuous Activities in the Lactonase SsoPox. PLoS ONE 8, e75272 (2013). 12. Reetz, M. T. Laboratory Evolution of Stereoselective Enzymes: A Prolific Source of Catalysts for Asymmetric Reactions. Angew. Chem. Int. Ed. 50, 138–174 (2011). 13. Daudé, D., Topham, C. M., Remaud‐Siméon, M. & André, I. Probing impact of active site residue mutations on stability and activity of Neisseria polysaccharea amylosucrase. Protein Sci. 22, 1754–1765 (2013). 14. Alcolombri, U., Elias, M. & Tawfik, D. S. Directed Evolution of Sulfotransferases and Paraoxonases by Ancestral Libraries. J. Mol. Biol. 411, 837–853 (2011). 15. Gonzalez, D. et al. Ancestral mutations as a tool for solubilizing proteins: The case of a hydrophobic phosphate-binding protein. FEBS Open Bio 4, 121–127 (2014). 16. Hawwa, R., Larsen, S. D., Ratia, K. & Mesecar, A. D. Structure-Based and Random Mutagenesis Approaches Increase the Organophosphate-Degrading Activity of a Phosphotriesterase Homologue from Deinococcus radiodurans. J. Mol. Biol. 393, 36–57 (2009). 17. Herman, A. & Tawfik, D. S. Incorporating Synthetic Oligonucleotides via Gene Reassembly (ISOR): a versatile tool for generating targeted libraries. Protein Eng. Des. Sel. 20, 219–226 (2007). 18. Rochu, D., Chabrière, E. & Masson, P. Human paraoxonase: A promising approach for pre-treatment and therapy of organophosphorus poisoning. Toxicology 233, 47–59 (2007). 19. Krieger, R. Handbook of Pesticide Toxicology, Two-Volume Set: Principles and Agents. (Academic Press, 2001). 20. Gupta, R. C. Handbook of Toxicology of Chemical Warfare Agents. (Academic Press, 2009). 21. Sapozhnikova, Y. et al. Evaluation of pesticides and metals in fish of the Dniester River, Moldova. Chemosphere 60, 196–205 (2005). 22. Sapozhnikova, Y., Bawardi, O. & Schlenk, D. Pesticides and PCBs in sediments and fish from the Salton Sea, California, USA. Chemosphere 55, 797–809 (2004). 23. LeJeune, K. E., Wild, J. R. & Russell, A. J. Nerve agents degraded by enzymatic foams. Nature 395, 27–28 (1998). 24. Jacquet, P. et al. Current and emerging strategies for organophosphate decontamination: special focus on hyperstable enzymes. Environ. Sci. Pollut. Res. 1–19 https://doi.org/10.1007/s11356-016-6143-1 (2016). 25. Poirier, L., Jacquet, P., Elias, M., Daudé, D. & Chabrière, E. Decontamination of organophosphorus compounds: Towards new alternatives. Ann. Pharm. Fr. https://doi.org/10.1016/j.pharma.2017.01.004 (2017). 26. Afriat, L., Roodveldt, C., Manco, G. & Tawfik, D. S. The Latent Promiscuity of Newly Identified Microbial Lactonases Is Linked to a Recently Diverged Phosphotriesterase†. Biochemistry (Mosc.) 45, 13677–13686 (2006). 27. Dumas, D. P., Caldwell, S. R., Wild, J. R. & Raushel, F. M. Purification and properties of the phosphotriesterase from Pseudomonas diminuta. J. Biol. Chem. 264, 19659–19665 (1989). 28. Seibert, C. M. & Raushel, F. M. Structural and Catalytic Diversity within the Amidohydrolase Superfamily†. Biochemistry (Mosc.) 44, 6383–6391 (2005). 29. Roodveldt, C. & Tawfik, D. S. Shared Promiscuous Activities and Evolutionary Features in Various Members of the Amidohydrolase Superfamily. Biochemistry (Mosc.) 44, 12728–12736 (2005). 30. Hiblot, J., Gotthard, G., Chabriere, E. & Elias, M. Structural and Enzymatic characterization of the lactonase SisLac from Sulfolobus islandicus. PLoS ONE 7, e47028 (2012). 31. Porzio, E., Merone, L., Mandrich, L., Rossi, M. & Manco, G. A new phosphotriesterase from Sulfolobus acidocaldarius and its comparison with the homologue from Sulfolobus solfataricus. Biochimie 89, 625–636 (2007). 32. Merone, L., Mandrich, L., Rossi, M. & Manco, G. A thermostable phosphotriesterase from the archaeon Sulfolobus solfataricus: cloning, overexpression and properties. Extremophiles 9, 297–305 (2005). 33. Hawwa, R., Aikens, J., Turner, R. J., Santarsiero, B. D. & Mesecar, A. D. Structural basis for thermostability revealed through the identification and characterization of a highly thermostable phosphotriesterase-like lactonase from Geobacillus stearothermophilus. Arch. Biochem. Biophys. 488, 109–120 (2009). 34. Xiang, D. F. et al. Functional Identification of Incorrectly Annotated Prolidases from the Amidohydrolase Superfamily of Enzymes. Biochemistry (Mosc.) 48, 3730–3742 (2009). 35. Chow, J. Y. et al. Directed Evolution of a Thermostable Quorum-quenching Lactonase from the Amidohydrolase Superfamily. J. Biol. Chem. 285, 40911–40920 (2010). 36. Hiblot, J., Gotthard, G., Champion, C., Chabriere, E. & Elias, M. Crystallization and preliminary X-ray diffraction analysis of the lactonase VmoLac from Vulcanisaeta moutnovskia. Acta Crystallograph. Sect. F Struct. Biol. Cryst. Commun. 69, 1235–1238 (2013). 37. Hiblot, J., Bzdrenga, J., Champion, C., Chabriere, E. & Elias, M. Crystal structure of VmoLac, a tentative quorum quenching lactonase from the extremophilic crenarchaeon Vulcanisaeta moutnovskia. Sci. Rep. 5, 8372 (2015). 38. Bzdrenga, J. et al. SacPox from the thermoacidophilic crenarchaeon Sulfolobus acidocaldarius is a proficient lactonase. BMC Res. Notes 7, 333 (2014). 39. Bzdrenga, J. et al. Biotechnological applications of quorum quenching enzymes. Chem. Biol. Interact. https://doi.org/10.1016/j. cbi.2016.05.028 (2016). 40. Rémy, B., Plener, L., Elias, M., Daudé, D. & Chabrière, E. Des enzymes pour bloquer la communication bactérienne, une alternative aux antibiotiques? Ann. Pharm. Fr. https://doi.org/10.1016/j.pharma.2016.06.005 (2016). 41. Vieille, C. & Zeikus, G. J. Hyperthermophilic Enzymes: Sources, Uses, and Molecular Mechanisms for Thermostability. Microbiol. Mol. Biol. Rev. 65, 1–43 (2001). 42. Elias, M. et al. Structural Basis for Natural Lactonase and Promiscuous Phosphotriesterase Activities. J. Mol. Biol. 379, 1017–1028 (2008). 43. Hoque, M. A. et al. Stepwise Loop Insertion Strategy for Active Site Remodeling to Generate Novel Enzyme Functions. ACS Chem. Biol. https://doi.org/10.1021/acschembio.7b00018 (2017). Scientific REPO R TS | 7: 16745 | DOI:10.1038/s41598-017-16841-0 14 www.nature.com/scientificreports/ 44. Meier, M. M. et al. Molecular Engineering of Organophosphate Hydrolysis Activity from a Weak Promiscuous Lactonase Template. J. Am. Chem. Soc. 135, 11670–11677 (2013). 45. Elias, M. et al. Crystallization and preliminary X-ray diffraction analysis of the hyperthermophilic Sulfolobus solfataricus phosphotriesterase. Acta Crystallograph. Sect. F Struct. Biol. Cryst. Commun. 63, 553–555 (2007). 46. Hiblot, J., Gotthard, G., Chabriere, E. & Elias, M. Characterisation of the organophosphate hydrolase catalytic activity of SsoPox. Sci. Rep. 2 (2012). 47. Del Giudice, I. et al. An efficient thermostable organophosphate hydrolase and its application in pesticide decontamination. Biotechnol. Bioeng. 113, 724–734 (2016). 48. Rémy, B. et al. Harnessing hyperthermostable lactonase from Sulfolobus solfataricus for biotechnological applications. Sci. Rep. 6 (2016). 49. Vecchio, P. D. et al. Structural determinants of the high thermal stability of SsoPox from the hyperthermophilic archaeon Sulfolobus solfataricus. Extremophiles 13, 461–470 (2009). 50. Afriat-Jurnou, L., Jackson, C. J. & Tawfik, D. S. Reconstructing a Missing Link in the Evolution of a Recently Diverged Phosphotriesterase by Active-Site Loop Remodeling. Biochemistry (Mosc.) 51, 6047–6055 (2012). 51. Jackson, C. J. et al. Anomalous scattering analysis of Agrobacterium radiobacter phosphotriesterase: the prominent role of iron in the heterobinuclear active site. Biochem. J. 397, 501 (2006). 52. Tawfik, D. S. Accuracy-rate tradeoffs: how do enzymes meet demands of selectivity and catalytic efficiency? Curr. Opin. Chem. Biol. 21, 73–80 (2014). 53. Xiang, D. F. et al. Functional Annotation and Three-Dimensional Structure of Dr0930 from Deinococcus radiodurans: A Close Relative of Phosphotriesterase in the Amidohydrolase Superfamily. Biochemistry (Mosc.) 48, 2237–2247 (2009). 54. Dellus-Gur, E. et al. Negative Epistasis and Evolvability in TEM-1 β-Lactamase—The Thin Line between an Enzyme’s Conformational Freedom and Disorder. J. Mol. Biol. 427, 2396–2409 (2015). 55. Reetz, M. T. The Importance of Additive and Non-Additive Mutational Effects in Protein Engineering. Angew. Chem. Int. Ed. 52, 2658–2666 (2013). 56. Campbell, E. et al. The role of protein dynamics in the evolution of new enzyme function. Nat. Chem. Biol. 12, 944–950 (2016). 57. Dellus-Gur, E., Toth-Petroczy, A., Elias, M. & Tawfik, D. S. What Makes a Protein Fold Amenable to Functional Innovation? Fold Polarity and Stability Trade-offs. J. Mol. Biol. 425, 2609–2621 (2013). 58. Guendouze, A. et al. Effect of quorum quenching lactonase in clinical isolates of Pseudomonas aeruginosa and comparison with quorum sensing inhibitors. Front. Microbiol. 8 (2017). 59. DeLano, W. L. The PyMOL molecular graphics system (2002). 60. Emsley, P. & Cowtan, K. Coot: model-building tools for molecular graphics. Acta Crystallogr. D Biol. Crystallogr. 60, 2126–2132 (2004). 61. Kabsch, W. Automatic processing of rotation diffraction data from crystals of initially unknown symmetry and cell constants. J. Appl. Crystallogr. 26, 795–800 (1993). 62. McCoy, A. J. et al. Phaser crystallographic software. J. Appl. Crystallogr. 40, 658–674 (2007). 63. Vagin, A. & Teplyakov, A. An approach to multi-copy search in molecular replacement. Acta Crystallogr. D Biol. Crystallogr. 56, 1622–1624 (2000). 64. Murshudov, G. N., Vagin, A. A. & Dodson, E. J. Refinement of Macromolecular Structures by the Maximum-Likelihood Method. Acta Crystallogr. D Biol. Crystallogr. 53, 240–255 (1997). Acknowledgements This work was granted by DGA, France (REI. 2009 34 0045). G.G. and P.J. are PhD students granted by DGA. Financial support was also obtained from the MnDrive Initiative (M.E.) to support this work. This work was supported by “Investissements d’avenir” program (Méditerranée Infection 10-IAHU-03) of the French Agence Nationale de la Recherche (ANR). Author Contributions P.J., J.H., D.D., C.B., N.A., G.G. performed the experiments, P.J., J.H., D.D., C.B., N.A., G.G., E.C., M.E. analyzed the data, P.J., J.H., D.D., C.B., M.E. wrote the manuscript and prepared the figures, E.C., D.D. and M.E. designed the experiments. Additional Information Supplementary information accompanies this paper at https://doi.org/10.1038/s41598-017-16841-0. Competing Interests: M.E. and E.C. have a patent WO2014167140 A1 licensed to Gene&GreenTK. D.D. and E.C. report personal fees from Gene&GreenTK during the conduct of the study. Publisher's note: Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations. Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Cre- ative Commons license, and indicate if changes were made. The images or other third party material in this article are included in the article’s Creative Commons license, unless indicated otherwise in a credit line to the material. If material is not included in the article’s Creative Commons license and your intended use is not per- mitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this license, visit http://creativecommons.org/licenses/by/4.0/. © The Author(s) 2017 Scientific REPO R TS | 7: 16745 | DOI:10.1038/s41598-017-16841-0 15

Journal

Scientific Reports – Springer Journals

Published: Dec 1, 2017

There are no references for this article.