Characterization and engineering of a two-enzyme system for plastics depolymerization
Edited by Alexis T. Bell, University of California, Berkeley, CA, and approved August 26, 2020 (received for review April 11, 2020)
Significance
Deconstruction of recalcitrant polymers, such as cellulose or chitin, is accomplished in nature by synergistic enzyme cocktails that evolved over millions of years. In these systems, soluble dimeric or oligomeric intermediates are typically released via interfacial biocatalysis, and additional enzymes often process the soluble intermediates into monomers for microbial uptake. The recent discovery of a two-enzyme system for polyethylene terephthalate (PET) deconstruction, which employs one enzyme to convert the polymer into soluble intermediates and another enzyme to produce the constituent PET monomers (MHETase), suggests that nature may be evolving similar deconstruction strategies for synthetic plastics. This study on the characterization of the MHETase enzyme and synergy of the two-enzyme PET depolymerization system may inform enzyme cocktail-based strategies for plastics upcycling.
Abstract
Plastics pollution represents a global environmental crisis. In response, microbes are evolving the capacity to utilize synthetic polymers as carbon and energy sources. Recently, Ideonella sakaiensis was reported to secrete a two-enzyme system to deconstruct polyethylene terephthalate (PET) to its constituent monomers. Specifically, the I. sakaiensis PETase depolymerizes PET, liberating soluble products, including mono(2-hydroxyethyl) terephthalate (MHET), which is cleaved to terephthalic acid and ethylene glycol by MHETase. Here, we report a 1.6 Å resolution MHETase structure, illustrating that the MHETase core domain is similar to PETase, capped by a lid domain. Simulations of the catalytic itinerary predict that MHETase follows the canonical two-step serine hydrolase mechanism. Bioinformatics analysis suggests that MHETase evolved from ferulic acid esterases, and two homologous enzymes are shown to exhibit MHET turnover. Analysis of the two homologous enzymes and the MHETase S131G mutant demonstrates the importance of this residue for accommodation of MHET in the active site. We also demonstrate that the MHETase lid is crucial for hydrolysis of MHET and, furthermore, that MHETase does not turnover mono(2-hydroxyethyl)-furanoate or mono(2-hydroxyethyl)-isophthalate. A highly synergistic relationship between PETase and MHETase was observed for the conversion of amorphous PET film to monomers across all nonzero MHETase concentrations tested. Finally, we compare the performance of MHETase:PETase chimeric proteins of varying linker lengths, which all exhibit improved PET and MHET turnover relative to the free enzymes. Together, these results offer insights into the two-enzyme PET depolymerization system and will inform future efforts in the biological deconstruction and upcycling of mixed plastics.
Sign up for PNAS alerts.
Get alerts for new articles, or get an alert when an article is cited.
Synthetic polymers pervade all aspects of modern life, due to their low cost, high durability, and impressive range of tunability. Originally developed to avoid the use of animal-based products, plastics have now become so widespread that their leakage into the biosphere and accumulation in landfills is creating a global-scale environmental crisis. Indeed, plastics have been found widespread in the world’s oceans (1–4), in the soil (5), and more recently, microplastics have been observed entrained in the air (6). The leakage of plastics into the environment on a planetary scale has led to the subsequent discovery of multiple biological systems able to convert man-made polymers for use as a carbon and energy source (7–12). These plastic-degrading systems offer a starting point for biotechnology applications toward a circular materials economy (12–16).
Among synthetic polymers manufactured today, polyethylene terephthalate (PET) is the most abundant polyester, which is made from petroleum-derived terephthalic acid (TPA) and ethylene glycol (EG). Given the prevalence of esterase enzymes in nature, PET biodegradation has been studied for nearly two decades, with multiple cutinase enzymes reported to perform depolymerization (17–26). In 2016, Yoshida et al. (10) reported the discovery and characterization of the soil bacterium, Ideonella sakaiensis 201-F6, which employs a two-enzyme system to depolymerize PET to TPA and EG, which are further catabolized as a carbon and energy source. Characterization of I. sakaiensis revealed the PETase enzyme, which is a cutinase-like serine hydrolase that attacks the PET polymer, liberating bis-(hydroxyethyl) terephthalate (BHET), mono(2-hydroxyethyl) terephthalate (MHET), and TPA. PETase cleaves BHET to MHET and EG, and the soluble MHET product is further hydrolyzed by MHETase to produce TPA and EG. Multiple crystal structures and biochemical studies of I. sakaiensis PETase (27–33) revealed an open active-site architecture that is able to bind to PET oligomers. The PETase enzyme likely follows the canonical serine hydrolase catalytic mechanism (34), but open questions remain regarding the mobility of certain residues during the catalytic cycle (27).
Conversely, the structure and function of the MHETase enzyme is far less characterized, with only two published studies focused on MHETase structure and engineering to date (35, 36). These studies report structures at 2.1 to 2.2 Å resolution, wherein the similarity to ferulic acid esterase (FAE) is noted (37, 38). Informed by these structures, engineering efforts aimed to improve turnover of BHET by MHETase, which is a nonnative substrate of the wild-type. Beyond these studies and the original report of MHETase from Yoshida et al. (10), several questions remain regarding the MHETase mechanism and PETase-MHETase synergy. To that end, here we combine structural, computational, biochemical, and bioinformatics approaches to reveal molecular insights into the MHETase structure, mechanism of hydrolysis, the evolution of MHETase activity from FAEs, and engineering of the two-enzyme system for PET depolymerization.
Results
Structural Characterization of MHETase Reveals a Core Domain Similar to That of PETase.
Four crystal structures of MHETase were obtained with the highest-resolution data (6QZ3) extending to 1.6 Å with a benzoate molecule in the active site (Fig. 1 and SI Appendix, Table S1). These data reveal a catalytic domain that adopts the α/β-hydrolase fold typical of a serine hydrolase (Fig. 1A, light gray), with an extensive lid domain (Fig. 1A, dark gray) that partially covers the active site and hosts a well-coordinated Ca2+ cation. A similar Ca2+ binding site was characterized for Aspergillus oryzae FaeB (AoFaeB), wherein it was hypothesized to have a role in stabilizing the lid domain (38). Overall, the structure of MHETase is most similar to that of FAEs, as discussed previously (35, 36). The structural conservation between the hydrolase domains of MHETase and PETase is striking (Fig. 1D and SI Appendix, Fig. S1), and despite the large insertion of ∼240 residues representing the lid domain, residues Ser225, Asp492, and His528 effectively reconstitute the catalytic triad (Fig. 1B). In fact, the terminal residues of the lid domain converge to within hydrogen-bonding distance of each other (Tyr252-Ala469, 2.9 Å), creating a compact linkage to the hydrolase domain. The lid domain of MHETase is exceptionally large, as average lid domains in α/β-hydrolases tend to be ∼100 residues (34), and is more typical of a lid from tannase family members (vide infra). The equivalent connection site in PETase is occupied by a seven-residue loop.
Fig. 1.
In addition, we determined two apo structures with alternative packing (6QZ2 and 6QZ4), one structure with a fully occupied benzoic acid ligand (6QZ3), and one with partially occupied benzoic acid (6QZ1). We observed that residue Phe415 adopts a “closed” orientation on substrate binding consistent with prior substrate bound structures (PDB ID codes 6QGA and 6QGB) (35), and the partially occupied site results in an intermediate dual “open/closed” conformation (SI Appendix, Fig. S2 A–C) (35). The only other amino acid with side-chain positioning correlated with ligand binding is Gln410. When Phe415 is in the open position, the side chain of Gln410 pivots toward the active site to a position wherein the heavy atoms would be as close as 1.8 Å to those of Phe415 if it were in the closed conformation (SI Appendix, Fig. S2D).
Given the difference in overall isoelectric point (pI) between PETase (9.65) and MHETase (5.11), we generated electrostatic surface profiles for comparison (Fig. 1E). As previously reported, PETase has a highly polarized surface charge (29), whereas MHETase exhibits a more heterogeneous and acidic surface charge distribution. MHETase contains five disulfide bonds (Fig. 1A). One of the MHETase disulfides is located at the active site, connecting cysteines (Cys224 and Cys529) adjacent to the catalytic residues (Ser225 and His528, respectively), which is conserved in tannase family members (38). PETase lacks a structurally equivalent disulfide, and the aligning residues in PETase (Trp159 and Ser238) are the same two residues mutated by Austin et al. (29) to yield a PETase substrate binding groove similar to that of cutinases, resulting in improved activity on a crystalline PET substrate.
Molecular Simulations of the MHETase Reaction Predict Deacylation Is Rate-Limiting.
The MHETase structure suggests a serine hydrolase mechanism for MHET hydrolysis (34). To elucidate the detailed reaction mechanism, we first constructed a Michaelis complex in silico utilizing the CHARMM (39) molecular simulation package (details in SI Appendix, Supplementary Materials and Methods). To examine MHETase dynamics and ligand stability, classic molecular dynamics (MD) simulations were conducted with NAMD (40) (all simulations totaling 2.25 µs) utilizing the CHARMM forcefield (41). Given the observed dual occupancy for Phe415 positioning in the crystal structures, we simulated in triplicate (each simulation 150 ns in length) the four combinations of Phe415 position (open and closed) and active site occupancy (empty active site and MHET-bound). In each case wherein Phe415 begins in the closed state (starting configurations from 6QZ3 structure, with coordinates for residues 56 to 61 from 6QZ4), Phe415 opens in the first 10 ns and rarely returns to the closed state; simulations that begin with Phe415 open (built from 6QZ4 structure) all remain open. To examine the effect of calcium binding, a fifth scenario absent of either MHET or Ca2+ was modeled in triplicate 150-ns simulations (the prior four scenarios each include bound Ca2+). These trajectories show evidence for lid stabilization upon Ca2+ binding mainly in the immediate vicinity of the calcium binding site (SI Appendix, Fig. S3). When bound at the active site, the carboxylate motif of MHET exhibits stable hydrogen bonds with Arg411 and Ser416, while the carbonyl is stabilized via hydrogen bonds to the oxyanion hole residues, Glu226 and Gly132 (SI Appendix, Figs. S4 and S5). In all three simulations with MHET bound and Phe415 open, MHET maintains these interactions and remains bound at the active site throughout the duration of the 150-ns simulation and primed for hydrolysis (hydrogen bond distance between Ser225 and His528 = 2.0 ± 0.2 Å; nucleophilic attack distance between Ser225 and MHET = 3.1 ± 0.3 Å; hydrogen bond distance between Asp492 and His528 = 1.8 ± 0.1 Å). Further analysis of the MD simulations, including time traces for these important interactions, is available in the SI Appendix, Supplementary Materials and Methods.
Serine hydrolases catalyze a two-step reaction involving formation of an acyl-enzyme intermediate (acylation) that is released hydrolytically in the second step (deacylation) (34). We utilized the Amber software package (42) to perform hybrid quantum mechanics/molecular mechanics (QM/MM) 2D umbrella sampling with semiempirical force field SCC-DFTB (43) to study the catalytic steps. Judicious selection of a reaction coordinate is critical for kinetically meaningful barrier calculations. We chose the forming and breaking C-O bonds to map the free-energy landscape for both reaction steps, informed by transition path sampling studies of other hydrolase enzymes (44, 45).
In acylation, the catalytic serine (Ser225) is deprotonated by His528, activating it for nucleophilic attack upon the carbonyl C of MHET, liberating EG and forming the acyl-enzyme intermediate (AEI) (Fig. 2 A–C and Movie S1). The minimum free-energy path (MFEP) computed from QM/MM 2D umbrella sampling (along the forming C-O bond between the MHET carbonyl C and Ser225 and the breaking MHET C-O ester bond) predicts an acylation free-energy barrier (∆G‡) of 13.9 ± 0.17 kcal/mol with an overall reaction free energy (∆Greaction) of −5.2 ± 0.04 kcal/mol (Fig. 2G and SI Appendix, Fig. S6A). Although serine hydrolases have at times been considered to proceed through metastable tetrahedral intermediates along the reaction pathway for acylation and deacylation (46–48), the acylation MFEP calculated from 2D umbrella sampling does not indicate intermediate configurations with metastability.
Movie S1.
Fig. 2.
Following cleavage of EG from MHET, classic MD simulation reveals that EG leaves the active site in the presence of the AEI (Fig. 2H and Movie S2). In one simulation, EG initially maintains a hydrogen bond with His528 for ∼100 ps, then dislodges from the active site, and is free in solution within 1 ns. Three identical simulations were initiated, and EG exits the active site within 4 ns in each. An important implication of this observation is that the deacylation reaction proceeds without EG in the active site. This allows greater access for water molecules to approach the charged nitrogen of His528 for deacylation (SI Appendix, Fig. S7).
Movie S2.
Deacylation involves nucleophilic attack by a water molecule on the AEI, liberating TPA (Fig. 2 D–F and Movie S3). His528 plays a similar role as in acylation, deprotonating the catalytic water and transferring this proton to the catalytic serine, regenerating Ser225 for another catalytic cycle. The MFEP computed from QM/MM 2D umbrella sampling (along the forming C-O bond between MHET and water and the breaking AEI C-O bond) reveals a deacylation free-energy barrier (∆G‡) of 19.8 ± 0.10 kcal/mol and an overall reaction free energy (∆Greaction) of 2.6 ± 0.07 kcal/mol (Fig. 2I and SI Appendix, Fig. S6B). Together, the two catalytic steps are exergonic by −2.6 ± 0.08 kcal/mol. Deacylation is predicted to be the rate-limiting step, with a rate of 7.1 ± 1.1 × 10−2 s−1 (from transition state theory, at 30 °C, and assuming a transmission coefficient of 1), more than four orders-of-magnitude slower than acylation (1.02 ± 0.28 × 103 s−1). As in acylation, metastable configurations along the MFEP are not observed.
Movie S3.
Bioinformatics Analysis Suggests That MHETase Evolved from an FAE.
Beyond structural and mechanistic investigations, we were also interested in understanding potential MHETase evolutionary ancestry and identifying other MHET-active enzymes from natural diversity. MHETase belongs to the tannase family (PFAM ID: PF07519), which consists of fungal and bacterial FAEs, fungal and bacterial tannases, and several bacterial homologs of unknown function (49). To elucidate sequence relationships between MHETase and tannase family enzymes, we performed bioinformatic analyses of 6,671 tannase family sequences retrieved from the National Center for Biotechnology Information (NCBI) via PSI-BLAST (50). MHETase shares low sequence similarity (<53%) with most sequences in the family, with the exception of homologs from Comamonas thiooxydans strains DS1, DF1, and DF2 (strain: NCBI:txid363952, protein:GenBank WP_080747404.1) (51) and Hydrogenophaga sp. PML113 (strain: NCBI:txid1899350, protein:GenBank WP_083293388.1), which exhibit 81% and 73% identity to MHETase, respectively (Fig. 3A). Since initial identification of the homologous C. thiooxydans sequence (WP_080747404.1), this entry was removed from GenBank, as discussed in SI Appendix, Supplementary Materials and Methods.
Fig. 3.
Using the multiple sequence alignment of 6,671 tannase family sequences, we performed conservation analysis with MHETase sequence positions as a reference (SI Appendix, Figs. S8 and S9), which shows that most positions in the active site are highly conserved. Notable exceptions are positions 257, 411, 415, and 416, which exhibit low conservation scores and are less conserved than 80% of MHETase positions overall (SI Appendix, Fig. S8 B and C). It is noteworthy (vide infra) that position 131 is a well-conserved glycine in 91% of tannase family sequences but serine appears at position 131 in MHETase (Fig. 3B). Furthermore, the 10 cysteine positions in MHETase that form 5 disulfide bonds are highly conserved in the tannase family (SI Appendix, Fig. S10A). Although a sixth disulfide bond exists in AoFaeB (38), less than 8% of tannase family sequences exhibit this sixth disulfide bond, and the sixth disulfide bond positions are variable among this set (SI Appendix, Fig. S10B). One cysteine of the sixth disulfide bond in AoFaeB is a single residue variation found in MHETase (38), whereas the other sits on a loop where a 15-residue deletion is found in MHETase.
With this large dataset, we further conducted phylogenetic analysis of 120 sequences selected from tannase family sequences that were clearly annotated as tannases or FAEs in GenBank, including MHETase (SI Appendix, Table S2). In the phylogenetic tree (SI Appendix, Fig. S11), bacterial and fungal enzymes form paraphyletic groups, and within these groups, there are separate FAE and tannase subgroups. MHETase and the C. thiooxydans and Hydrogenophaga sp. PML113 homologs are found within a group of proteobacterial FAEs (bootstrap value > 95%). In addition, when separate profile hidden Markov models (pHMM) are constructed with the annotated tannase family FAE and tannase sequences (52), and then aligned with MHETase, a higher alignment score is achieved with the FAE pHMM than with the tannase pHMM (456.8 vs. 396.8), suggesting that MHETase is more similar to FAEs than tannases.
Biochemical Characterization of Active-Site MHETase Mutants and Homologs Reveals Important Residues for MHET Hydrolytic Activity.
From the bioinformatics analyses, we selected the MHETase homologs from C. thiooxydans and Hydrogenophaga sp. PML113 (Fig. 3A) to test for MHET hydrolysis activity, which, along with the I. sakaiensis MHETase, were produced in Escherichia coli and purified. Activity assays were performed for each enzyme to determine MHET turnover rates (Fig. 3E). The turnover rate (kcat) for MHETase is 27.6 ± 2.6 s−1, as compared to 9.5 ± 0.8 s−1 and 3.8 ± 2.5 s−1 for the C. thiooxydans and Hydrogenophaga sp. PML113 enzymes, respectively. The enzymes were also evaluated over a range of substrate concentrations to determine the Michaelis–Menten kinetic parameters (Fig. 3 G–J and SI Appendix, Table S3). FAEs have been shown to exhibit concentration-dependent substrate inhibition (53) in addition to the likely product inhibition of the enzyme (35). MHETase and both homologs also display this behavior. Using a substrate inhibition model (details in SI Appendix, Supplementary Materials and Methods), evaluation of the substrate-dependent reaction kinetics shows that MHETase more efficiently accepts MHET as a substrate than the C. thiooxydans and Hydrogenophaga sp. PML113 homologs, demonstrated by a Km value of 23.17 ± 1.65 µM as compared to values of 174.70 ± 4.75 µM and 41.09 ± 3.38 µM, respectively (Fig. 3 G–I and SI Appendix, Table S3). However, MHETase is also the most susceptible to substrate inhibition with a Kk value of 307.30 ± 20.65 µM. Despite the difference in affinity for MHET, MHETase and the C. thiooxydans enzyme exhibit similar maximal reaction rates, while the enzyme from Hydrogenophaga sp. PML113 is slower. The MHETase reaction efficiency, reported as kcat/Km, is ∼10-fold higher than for the C. thiooxydans enzyme and ∼20-fold higher than the Hydrogenophaga sp. PML113 enzyme.
Homology models of both the C. thiooxydans and Hydrogenophaga sp. PML113 enzymes were constructed with SWISS-MODEL (54) using the MHETase structure as a template (PDB ID code 6QZ3), and the active site aligned with a modeled MHET-bound MHETase structure (Fig. 3D). As noted, position 131 is a serine in MHETase, but a glycine in the two homologs (C. thiooxydans, Gly179 and Hydrogenophga sp. PML113, Gly121) (Fig. 3B). The C. thiooxydans enzyme is otherwise identical within 6 Å of the docked MHET ligand, whereas the Hydrogenophaga sp. PML113 enzyme also exhibits a serine in the equivalent position to the MHETase residue Phe415 (Fig. 3C). An S131G mutant of MHETase was constructed to examine the role of this residue in MHET hydrolytic activity, and steady-state enzyme kinetics were evaluated. The MHETase S131G mutant does not demonstrate concentration-dependent substrate inhibition, as is observed for the wild-type enzyme, which is likely due to the poor affinity for MHET. The S131G mutant has a Km value approximately eightfold higher than that of wild-type MHETase and the reaction efficiency is reduced to ∼3% that of the wild-type (Fig. 3J and SI Appendix, Table S3), illustrating the importance of this residue in MHET turnover.
Focusing on residues within the coordination sphere of the docked MHET ligand, amino acids and their frequencies across the tannase family were compared to MHETase (SI Appendix, Figs. S8 and S9). In position 495, MHETase features a phenylalanine, while isoleucine is also a common residue in this position across the tannase family. Palm et al. (35) demonstrated that Phe495 has a modest effect on activity by mutation to alanine. We constructed and evaluated an MHETase mutant with isoleucine in this position (F495I), which dramatically impairs activity, lowering the turnover rate from 27.6 ± 2.6 s−1 to 1.3 ± 0.7 s−1 (Fig. 3E). In position 226, which is part of the conserved “lipase box” motif in serine hydrolases (55), MHETase exhibits a glutamate, while threonine and asparagine are more common among tannase family members. Mutation of this lipase box residue to threonine (E226T) yielded a ∼50% reduction in MHET activity relative to the wild-type MHETase (Fig. 3E). Mutation of the catalytic serine (S225A), as expected, produced an inactive enzyme.
Unique Structural Features between MHETase and PETase Determine Substrate Specificity and Stability.
Given the structural similarities of the MHETase and PETase core domains (Fig. 1 C and D), we were interested in understanding the role of unique MHETase features, namely the lid domain and the active site disulfide bond between Cys224 and Cys529, on substrate specificity and MHET hydrolytic activity. Accordingly, the lid was both added to PETase (“lidded PETase”) and removed from MHETase (“lidless MHETase”). Given the natural substrate specificities of wild-type PETase and MHETase, we hypothesized that the former could confer MHET activity, but abolish PET hydrolytic potential, whereas the latter was expected to have the opposite effect. The lidded PETase was created by replacing the seven-residue loop of PETase (Trp185:Phe191, PETase numbering) with Gly251:Thr472 from MHETase. In control experiments, wild-type PETase exhibited no detectable activity on MHET, and the lidded PETase is not able to degrade amorphous PET film. However, meager activity of lidded PETase was detected on MHET (kcat = 0.11 ± 0.02 s−1) (Fig. 3F). The lidless MHETase was created by replacing the MHETase lid (Gly251:Thr472) with the seven-residue loop of PETase (Trp185:Phe191, PETase numbering). This construction results in an exposed MHETase active site, possibly allowing for acquired PET hydrolytic activity. The resulting enzyme has a kcat value on MHET of 0.05 ± 0.03 s−1, 1,000-fold lower than the rate for wild-type MHETase, demonstrating that the lid domain is crucial for MHET hydrolytic activity (Fig. 3F). The lidless MHETase enzyme was also unable to degrade amorphous PET film over 96 h, despite the more accessible active site.
Similarly, variants of lidless MHETase were generated to remove the active site disulfide and replace the two sites with tryptophan and serine (lidless MHETase C224W/C529S) (SI Appendix, Fig. S1B) to reconstitute the PETase active site, or with histidine and phenylalanine (lidless MHETase C224H/C529F), matching the active site of the double-mutant PETase variant previously shown to exhibit improved PET hydrolytic activity on crystalline PET (29). The lidless MHETase C224W/C529S mutant, which reconstitutes the wild-type active site motif of PETase, displays the same turnover rate (within error) as the lidded PETase mutant on MHET (kcat = 0.10 ± 0.06 s−1), while the lidless MHETase C224H/C529F mutant is even less active on MHET (kcat = 0.06 ± 0.03 s−1) (Fig. 3F). We also generated a PETase variant to recreate the active site disulfide found in MHETase (PETase W159C/S238C). The PETase mutant exhibited very low MHET hydrolytic activity (kcat = 0.03 ± 0.3 s−1), and similarly had no activity on BHET or amorphous PET film.
To delineate the effects of engineering the lid and removing the active site disulfide bond, we also generated three MHETase mutants altering only the active site disulfide motif (C224A/C529A, C224W/C529S, and C224H/C529F). We hypothesized that removal of this disulfide bond may diminish the thermal stability of MHETase. However, each of these variants either expressed in inclusion bodies or did not express at all. Attempts were also made to introduce disulfide motifs into MHETase that are found in PETase (G489C/S530C) or in AoFaeB. To recapitulate the AoFaeB disulfide, the mutations include both a point mutation (S136C) as well as the insertion of a 15-residue loop from AoFaeB that harbors the partnering cysteine residue. As with the active site disulfide mutants, these mutants either expressed in inclusion bodies or did not express at all. A variant was also created that included both the PETase-like disulfide (G489C/S530C) and the AoFaeB modification (S136C with 15-residue loop from AoFaeB). This last variant, with seven total disulfides, was successfully expressed and had very low activity on MHET (kcat = 0.16 ± 0.14 s−1) (Fig. 3F). A full list of clones, mutants, and primers can be found in Dataset S1.
MHETase Is Catalytically Inactive on MHE-Isophthalate and MHE-Furanoate.
We evaluated the substrate specificity of MHETase using the monohydroxyethyl monomer unit of two additional compounds. Specifically, assays were performed with mono(2-hydroxyethyl)-isophthalate (MHEI) (SI Appendix, Fig. S12) and mono(2-hydroxyethyl)-furanoate (MHEF) (SI Appendix, Fig. S13). Isophthalate is a common comonomer in industrial PET formulations used to modify crystallinity, such that MHEI could be released from polyester depolymerization. PETase has been demonstrated to deconstruct other aromatic polyesters (29), including polyethylene furanoate (PEF), yielding MHEF as a product of the enzymatic hydrolysis reaction. Over the course of 24 h at 30 °C, no MHETase activity was detected for either substrate using substrate concentrations from 25 to 250 μM, in contrast with complete hydrolysis of MHET (SI Appendix, Fig. S14) in the same time using identical reaction conditions.
To explain the inability of MHETase to act on MHEI and MHEF, we conducted flexible ligand/flexible receptor docking simulations and predicted 10 binding orientations for each molecule (MHET, MHEI, and MHEF) in MHETase. These docking simulations indicate that MHET binds to MHETase with a binding free energy of −7.13 kcal/mol and in a catalytically primed configuration. This binding mode features the carbonyl C of MHET within 3.2 Å of Ser225-O, which itself is within 2.90 Å of His528-N(e) and His528-N(d) is 3.93 Å from Asp492-O (SI Appendix, Fig. S15). For MHEI and MHEF, no binding modes were predicted that exhibit similarly favorable binding free energies, feature the MHET carbonyl C within range for attack by Ser225, and stabilize the carbonyl of the ester in the oxyanion hole, suggesting that MHETase will not readily act on these molecules.
PETase and MHETase Act Synergistically during PET Depolymerization.
While MHET is susceptible to hydrolysis by a number of PET-degrading cutinases, I. sakaiensis requires the action of two enzymes for PET degradation to liberate TPA and EG (10). Given the turnover rates for MHETase reported here, depolymerization by PETase is likely the rate-limiting step when the enzymes are employed together. To investigate the action of the two-enzyme system, we thus measured the extent of hydrolysis of a commercial amorphous PET substrate over 96 h at 30 °C using PETase and MHETase at varying concentrations (Fig. 4A and SI Appendix, Table S4). As expected, MHETase alone has no activity on PET film. Over the range of enzyme loadings tested (0 to 2.0 mg enzyme/g PET), degradation by PETase alone, as determined by concentration of product released (the sum of BHET, MHET, and TPA), scales with enzyme loading. Upon addition of MHETase in the reaction, at any loading tested (0.1 to 1.0 mg MHETase/g PET), product release still scales with PETase loading, but at a markedly higher level than with PETase alone (Fig. 4A). The overall trend of degradation within the range of enzyme loadings tested, which shows increasing levels of constituent monomers released as concentration of both enzymes increases, is indicative that these reactions are enzyme-limited under these conditions, rather than substrate-limited. The synergy study does not strongly indicate that any particular ratio of PETase to MHETase results in optimal degradation over the enzyme loadings tested, but rather that degradation scales with PETase loading and the presence of MHETase, even at low concentrations relative to PETase, improves total degradation.
Fig. 4.
Chimeric Proteins of MHETase and PETase Improves PET Degradation and MHET Hydrolysis Rates.
In light of the highly synergistic relationship between PETase and MHETase on amorphous PET, where increasing loading of each enzyme results in more constituent monomer release, we next examined how proximity of the two enzymes influences hydrolytic activity. Chimeric proteins covalently linking the C terminus of MHETase to the N terminus of PETase, using flexible glycine-serine linkers of 8, 12, and 20 total glycine and serine residues, were generated and assayed for degradation of amorphous PET (Fig. 4B). Varying linker lengths were explored to understand the effect of increased mobility between the two domains (56). Furthermore, two enzyme loadings were compared: The lower loading corresponding to ∼0.08 mg PETase/g PET and 0.16 mg MHETase/g PET and the higher enzyme loading corresponding to 0.25 mg PETase/g PET and 0.5 mg MHETase/g PET (Fig. 4 C and D). At both loadings, when comparing the extent of degradation achieved by PETase alone, MHETase alone, and an equimolar mix of PETase and MHETase, the chimeric proteins outperform PETase, as well as the mixed reaction containing both PETase and MHETase unlinked in solution. Furthermore, the chimeras demonstrate a higher catalytic activity on MHET (Fig. 4E). Chimeric constructs linking the C terminus of PETase to the N terminus of MHETase did not successfully express protein (Fig. 4B). SEM analysis of digested amorphous PET film confirms degradation (SI Appendix, Fig. S16).
Discussion
The ability to degrade polymers to their monomeric units is important for subsequent reuse in new products, which is a critical technical advance needed to enable a global circular materials economy. In biological systems, complete depolymerization to monomers can be necessary for microbial uptake and growth, as in I. sakaiensis wherein MHETase is the enzymatic partner to PETase, together allowing for the complete degradation of PET to TPA and EG for catabolism (10). Prior studies presenting MHETase crystal structures focused upon understanding and tuning substrate specificity, particularly the rational engineering of MHETase to impart BHET hydrolysis activity (35, 36). Drawing inspiration from our structural analyses, this complementary study offers further insights into the two-enzyme PETase/MHETase system.
The recent structural report from Palm et al. (35) highlighted several important amino acid contributions to substrate specificity in MHETase, specifically focusing on active site residues. Of note, they pointed out the importance of Phe415 for substrate binding via an “induced fit” mechanism and highlighted Arg411 with respect to hydrogen bonding of the MHET carboxylate group, both of which are proposed to be drivers of substrate specificity. In addition, beyond engineering a starting point for BHET activity in MHETase for further optimization, the potential for MHEF turnover was suggested, given the proposed utility of PEF as a bio-based PET replacement (57). In our previous work (29), we demonstrated that PETase effectively depolymerizes PEF, but the results here do not indicate the same for MHETase on MHEF, and docking simulations agree with the observed patterns in MHETase selectivity. Despite success with predicting a low-energy catalytically competent binding mode for MHET to MHETase, we were only able to predict one binding mode of MHEF to MHETase with the MHEF nucleophilic carbonyl in the oxyanion hole, but in this pose the MHEF carboxylate moiety is not in range to interact with Arg411, suggesting that further active site engineering will be necessary to enable MHEF turnover. Similarly, only one binding mode for MHEI was predicted wherein the catalytic triad was oriented for catalysis but, akin to MHEF, the nonlinearity of the molecule prevents simultaneous interaction with the oxyanion hole and R411.
The enzyme kinetics studies presented here reveal a substantial reduction in activity for the S131G, E226T, and F495I MHETase mutants, indicating that these positions play important roles in substrate specificity and catalytic efficiency. A previous study also demonstrated greatly reduced hydrolytic activity by a F415S variant (36). Additionally, two homologs identified via bioinformatics analysis from C. thiooxydans and Hydrogenophaga sp. PML113 exhibit extremely similar active site environments (Fig. 3D), with the only exceptions being variations at positions 131 and 415 (MHETase numbering), and these homologs display reaction efficiencies (kcat/Km) reduced by an order-of-magnitude (SI Appendix, Table S3). Furthermore, as the amino acids at these positions in wild-type MHETase are less common in tannase family sequences (SI Appendix, Fig. S9), and mutation to the more common amino acids led to a reduction in activity, this suggests that these two sequence positions were specifically evolved in MHETase to accommodate MHET.
Two-enzyme systems for complete PET degradation have been examined previously, either derived from a single microorganism (e.g., Thermobifida fusca) (58) or screened from multiple sources for optimal activity (25, 59). The enzyme synergy results for the I. sakaiensis PETase/MHETase system on amorphous PET display a clear performance improvement when MHETase is included in the reaction. Namely, overall degradation scales with PETase loading within the tested range (0 to 2.0 mg PETase/g PET), but the inclusion of MHETase in the degradation reaction markedly improves depolymerization and this synergistic enhancement also scales with MHETase loading.
The presence of confirmed MHETase homologs in C. thiooxydans and Hydrogenophaga sp. PML113 suggests that these bacteria may harbor abilities for TPA catabolism (SI Appendix, Fig. S17). Bioinformatics analysis was thus conducted to query the genomes of the strains compared to known TPA catabolic genes from I. sakaiensis (10), Comamonas sp. E6 (60, 61), Delftia tsuruhatensis (62), Paraburkholderia xenovorans (63), Rhodococcus jostii RHA1 (64), and Rhodococcus sp. DK17 (65), including putative PETases, terephthalate transporter genes, two-component terephthalate dioxygenases, the 1,2-dihydroxy-3,5-cyclohexadiene-1,4-dicarboxylate dehydrogenase, and the three types of protocatechuate (PCA) dioxygenases (PCA-2,3, PCA-3,4, and PCA-4,5-dioxygenases) (SI Appendix, Table S5). This analysis revealed that neither C. thiooxydans nor Hydrogenophaga sp. PML113 harbor putative PETase genes. Interestingly, both strains exhibit genes encoding for TPA catabolic enzymes and transporters highly homologous to those of I. sakaiensis, Comamonas sp. E6, and D. tsuruhatensis (in all cases above 60% identity) (SI Appendix, Figs. S18 and S19), suggesting that they are highly likely able to turnover TPA to PCA, a common central intermediate in aerobic aromatic catabolic pathways (66). Each strain also contains annotated PCA-4,5-dioxygenases (SI Appendix, Table S5). Further experimental work will be required to understand if either of these bacteria exhibit the ability to depolymerize PET, perhaps through another type of mechanism than via ester hydrolases, or perhaps like Comamonas sp. E6, they are primarily able to consume soluble, xenobiotic intermediates. Perhaps these strains could serve as useful sources of TPA catabolic genes for synthetic biology efforts associated with biological plastics recycling and upcycling (67).
The enzymatic deconstruction of recalcitrant natural polymers, such as cellulose, hemicellulose, and chitin, is accomplished in nature by the action of mixtures of synergistic enzymes secreted from microbes (68, 69). For example, as observed in fungal cellulase systems for cellulose depolymerization, these mixtures typically contain a subset of enzymes to act directly on solid polymeric substrates via interfacial enzyme mechanisms, and complementary enzymes (e.g., β-glucosidases) that further process solubilized intermediates to monomeric constituents (e.g., cellobiose hydrolysis to glucose). Given that natural microbial systems evolved over millions of years to optimally degrade recalcitrant polymers, perhaps it is thus not surprising, in hindsight, that a soil bacterium such as I. sakaiensis evolved the ability to utilize a crystalline polyester substrate with, to our collective knowledge, a two-enzyme system (10, 16). Extending the analogy of cellulase enzymes and plant cell wall deconstruction for breaking down diverse polysaccharides simultaneously, the concept of deconstructing synthetic polymers in the form of mixed plastics waste with advanced enzyme mixtures is an exciting research direction beyond PET to other polyesters, natural fibers (e.g., cellulose from cotton, proteins from wool) (70), polyamides, polyurethanes (71), and other polymers susceptible to enzymatic depolymerization. Going forward, the design of multienzyme systems for depolymerization of mixed polymer wastes is a promising and fruitful area for continued investigation.
Methods
Plasmid Construction.
pET21b(+)-based expression plasmids for I. sakaiensis genes, homologous genes, and mutants were generated as further described in Dataset S1.
Protein Expression and Purification.
E. coli-based protein expression and chromatographic purification is described in SI Appendix, Supplementary Materials and Methods.
Crystallization and Structure Determination.
MHETase was crystallized in four conditions, including a seleno-methionine–labeled version for single-wavelength anomalous diffraction phasing. All X-ray data collections were performed at Beamline I03 at the Diamond Light Source. Detailed methods and statistics are provided in SI Appendix, Supplementary Materials and Methods and Table S1.
Molecular Simulations.
MD simulations were performed for solvated MHETase both in the free state and with MHET bound at the active site. All systems were built in CHARMM (39) and simulations utilized the CHARMM forcefield (41). Classic MD simulations were run with NAMD (40); QM/MM simulations, including 2D umbrella sampling free-energy calculations, were run in Amber (42, 72). Additional simulation details are in SI Appendix, Supplementary Materials and Methods.
Bioinformatics.
A total of 6,671 tannase family sequences were retrieved via PSI-BLAST against the NCBI nonredundant database (50). Phylogenetic analyses were conducted with MEGA7 (73). Additional details are in SI Appendix, Supplementary Materials and Methods.
MHETase Kinetics and Turnover Experiments.
MHETase and mutant enzymes were incubated with MHET, MHEI, or MHEF and reactions quenched with methanol and a heat treatment at 85 °C for 10 min. Hydrolysis extent was measured by HPLC as described in SI Appendix, Supplementary Materials and Methods and Table S6.
Molecular Docking.
MHET, MHEI, and MHEF docking into MHETase were modeled and prepared using tools in Schrödinger. Substrate docking simulations were conducted using Induced Fit Docking simulations in Schrödinger as described in SI Appendix, Supplementary Materials and Methods.
Ligand Synthesis.
MHET, MHEI, and MHEF were prepared via the coupling and subsequent deprotection of a monotBoc-protected EG with the respective acyl chlorides as further described in SI Appendix, Supplementary Materials and Methods and Figs. S12–S14.
MHETase Synergy with PETase.
The effect of MHETase loading and PETase loading on amorphous PET film after 96 h was measured as total product release (MHET, BHET, and TPA) via HPLC, as described in SI Appendix, Supplementary Materials and Methods.
MHETase-PETase Chimeras.
Chimeric constructs covalently linking MHETase to PETase were generated and incubated with either MHET or amorphous PET film as described in SI Appendix, Supplementary Materials and Methods.
Data Availability
The atomic coordinates and structure factors have been deposited in the Protein Data Bank, https://www.wwpdb.org/ (PDB ID codes 6QZ1, 6QZ2, 6QZ3, and 6QZ4).
Acknowledgments
We thank staff at the Diamond Light Source for beamtime and support, and McKenzie Hawkins for technical support with HPLC analysis. This work was authored in part by Alliance for Sustainable Energy, the manager and operator of the National Renewable Energy Laboratory (NREL) for the US Department of Energy (DOE) under Contract DE-AC36-08GO28308. Funding was provided by the DOE, Office of Energy Efficiency and Renewable Energy, Advanced Manufacturing Office and Bioenergy Technologies Office. This work was performed as part of the Bio-Optimized Technologies to keep Thermoplastics out of Landfills and the Environment Consortium and was supported by the Advanced Manufacturing Office and Bioenergy Technologies Office under contract DE-AC36-08GO28308 with the NREL, operated by Alliance for Sustainable Energy, LLC. G.T.B., B.C.K., E.E., B.S.D., N.A.R., G.D., and C.W.J. also acknowledge the NREL Laboratory Directed Research and Development program for funding. J.E.M. acknowledges Research England for E3 funding and the Biotechnology and Biological Sciences Research Council for Grant BB/P011918/1. R.G. and H.P.A. were funded through an NREL subcontract and University of Portsmouth Faculty of Science bursary. F.L.K. acknowledges support from the NSF Graduate Research Fellowship Program (Grant 3900101301). This work is also based on work supported by (while C.M.P. is serving at) the NSF. Any opinions, findings and conclusions or recommendations expressed in this material are those of the authors and do not necessarily reflect the views of the NSF. H.L.W. acknowledges the DOE for Grant DE-SC0011297TDD, the NIH for Grant 1R01GM129519-01, and the NSF for Grant CHE-1464946. J.E.G. acknowledges support from NSF Grant CBET-1552355 (to C.M.P.). V.C. and E.T. acknowledge the NSF for Grant MCB-1714556. Computer time was provided by Extreme Science and Engineering Discovery Environment allocation MCB-090159 at the San Diego Supercomputing Center and the Texas Advanced Computing Center, and by the National Renewable Energy Laboratory Computational Sciences Center supported by the DOE Office of Energy Efficience and Renewable Energy under Contract DE-AC36-08GO28308.
Supporting Information
Appendix (PDF)
- Download
- 3.08 MB
Dataset_S01 (XLSX)
- Download
- 19.35 KB
Movie S1.
MHETase acylation reaction. Aimless Shooting trajectory demonstrating the basic elements of the acylation reaction of MHET by I. sakaiensis MHETase. In acylation, the acyl-enzyme intermediate (AEI) is formed, and ethylene glycol is released. MHET substrate is shown in purple-colored sticks, and the catalytic residues (His528, Asp492, and nucleophile Ser225) are shown in white sticks. MHETase backbone is shown in transparent white cartoon. Total simulation time shown is 2 ps. These path sampling simulations constitute completely unrestrained instances of the chemical reaction and were undertaken purely for illustrative purposes. Most hydrogens not directly involved in the reaction are not shown for clarity; the one exception is the hydrogen on His528 that hydrogen bonds with Asp492 throughout acylation. Movie was created in VMD (57).
- Download
- 5.48 MB
Movie S2.
Release of acylation product ethylene glycol. Molecular dynamics (MD) trajectory of the acyl-enzyme intermediate (AEI) immediately following the acylation reaction. The AEI features a covalent bond between the nucleophile Ser225 and the terephthalic acid moiety of MHET. Upon leaving the enzyme active site, EG is fully immersed in water, which is present in the simulation but not shown for clarity. Color schemes and representations are the same as in Supplementary Movie S1, with the exception that the protein backbone shown in cartoon is now opaque. The simulation length is 1 ns and is performed with a purely classical forcefield (no QM region). 500 frames from this simulation are utilized in the movie; coordinates for the cartoon protein backbone are smoothed over 6 frames, while coordinates for the stick representations are smoothed over 4 frames. Movie was created in VMD (57).
- Download
- 3.58 MB
Movie S3.
MHETase deacylation reaction. Aimless Shooting trajectory demonstrating the basic elements of the deacylation reaction. The acyl-enzyme intermediate (AEI) formed in the acylation reaction is broken and terephthalic acid (TPA) is released. The catalytic residues Ser225 and His528 are regenerated in this step and TPA released. Color schemes and representations are the same as in Supplementary Movie S1. Total simulation time shown is 2 ps. Most hydrogens not directly involved in the reaction are not shown for clarity; the one exception is the hydrogen on His528 that hydrogen bonds with Asp492 throughout deacylation. Movie was created in VMD (57).
- Download
- 5.30 MB
References
1
K. L. Law et al., Plastic accumulation in the North Atlantic subtropical gyre. Science 329, 1185–1188 (2010).
2
A. Cózar et al., Plastic debris in the open ocean. Proc. Natl. Acad. Sci. U.S.A. 111, 10239–10244 (2014).
3
J. R. Jambeck et al., Marine pollution. Plastic waste inputs from land into the ocean. Science 347, 768–771 (2015).
4
J. B. Lamb et al., Plastic waste associated with disease on coral reefs. Science 359, 460–462 (2018).
5
A. A. de Souza Machado, W. Kloas, C. Zarfl, S. Hempel, M. C. Rillig, Microplastics as an emerging threat to terrestrial ecosystems. Glob. Change Biol. 24, 1405–1416 (2018).
6
S. Allen et al., Atmospheric transport and deposition of microplastics in a remote mountain catchment. Nat. Geosci. 12, 679 (2019).
7
J.-M. Restrepo-Flórez, A. Bassi, M. R. Thompson, Microbial degradation and deterioration of polyethylene—A review. Intl. Biodeter. Biodegrad. 88, 83–90 (2014).
8
J. Yang, Y. Yang, W.-M. Wu, J. Zhao, L. Jiang, Evidence of polyethylene biodegradation by bacterial strains from the guts of plastic-eating waxworms. Environ. Sci. Technol. 48, 13776–13784 (2014).
9
Y. Yang et al., Biodegradation and mineralization of polystyrene by plastic-eating mealworms: Part 2. Role of gut microorganisms. Environ. Sci. Technol. 49, 12087–12093 (2015).
10
S. Yoshida et al., A bacterium that degrades and assimilates poly(ethylene terephthalate). Science 351, 1196–1199 (2016).
11
P. Bombelli, C. J. Howe, F. Bertocchini, Polyethylene bio-degradation by caterpillars of the wax moth Galleria mellonella. Curr. Biol. 27, R292–R293 (2017).
12
P. Dvořák, P. I. Nikel, J. Damborský, V. de Lorenzo, Bioremediation 3.0: Engineering pollutant-removing bacteria in the times of systemic biology. Biotechnol. Adv. 35, 845–866 (2017).
13
N. Wierckx et al., Plastic waste as a novel substrate for industrial biotechnology. Microb. Biotechnol. 8, 900–903 (2015).
14
R. Wei, W. Zimmermann, Biocatalysis as a green route for recycling the recalcitrant plastic polyethylene terephthalate. Microb. Biotechnol. 10, 1302–1307 (2017).
15
R. Wei, W. Zimmermann, Microbial enzymes for the recycling of recalcitrant petroleum-based plastics: How far are we? Microb. Biotechnol. 10, 1308–1322 (2017).
16
I. Taniguchi et al., Biodegradation of PET: Current status and application aspects. ACS Catal. 9, 4089–4105 (2019).
17
R. Araújo et al., Tailoring cutinase activity towards polyethylene terephthalate and polyamide 6,6 fibers. J. Biotechnol. 128, 849–857 (2007).
18
Å. M. Ronkvist, W. Xie, W. Lu, R. A. Gross, Cutinase-catalyzed hydrolysis of poly (ethylene terephthalate). Macromolecules 42, 5128–5138 (2009).
19
E. Herrero Acero et al., Enzymatic surface hydrolysis of PET: Effect of structural diversity on kinetic properties of cutinases from Thermobifida. Macromolecules 44, 4632–4640 (2011).
20
D. Ribitsch et al., Characterization of a new cutinase from Thermobifida alba for PET-surface hydrolysis. Biocat. Biotrans. 30, 2–9 (2012).
21
S. Sulaiman et al., Isolation of a novel cutinase homolog with polyethylene terephthalate-degrading activity from leaf-branch compost by using a metagenomic approach. Appl. Environ. Microbiol. 78, 1556–1562 (2012).
22
D. Ribitsch et al., Fusion of binding domains to Thermobifida cellulosilytica cutinase to tune sorption characteristics and enhancing PET hydrolysis. Biomacromolecules 14, 1769–1776 (2013).
23
C. Roth et al., Structural and functional studies on a thermostable polyethylene terephthalate degrading hydrolase from Thermobifida fusca. Appl. Microbiol. Biotechnol. 98, 7815–7823 (2014).
24
D. Ribitsch et al., Enhanced cutinase-catalyzed hydrolysis of polyethylene terephthalate by covalent fusion to hydrophobins. Appl. Environ. Microbiol. 81, 3586–3592 (2015).
25
A. M. de Castro, A. Carniel, J. Nicomedes Junior, A. da Conceição Gomes, É. Valoni, Screening of commercial enzymes for poly(ethylene terephthalate) (PET) hydrolysis and synergy studies on different substrate sources. J. Ind. Microbiol. Biotechnol. 44, 835–844 (2017).
26
D. Danso et al., New insights into the function and global distribution of polyethylene terephthalate (PET)-degrading bacteria and enzymes in marine and terrestrial metagenomes. Appl. Environ. Microbiol. 84, e02773-17 (2018).
27
X. Han et al., Structural insight into catalytic mechanism of PET hydrolase. Nat. Commun. 8, 2106 (2017).
28
S. Joo et al., Structural insight into molecular mechanism of poly(ethylene terephthalate) degradation. Nat. Commun. 9, 382 (2018).
29
H. P. Austin et al., Characterization and engineering of a plastic-degrading aromatic polyesterase. Proc. Natl. Acad. Sci. U.S.A. 115, E4350–E4357 (2018).
30
T. Fecker et al., Active site flexibility as a hallmark for efficient PET degradation by I. sakaiensis PETase. Biophys. J. 114, 1302–1312 (2018).
31
M. Furukawa, N. Kawakami, K. Oda, K. Miyamoto, Acceleration of enzymatic degradation of poly (ethylene terephthalate) by surface coating with anionic surfactants. ChemSusChem 11, 4018–4025 (2018).
32
Y. Ma et al., Enhanced poly(ethylene terephthalate) hydrolase activity by protein engineering. Engineering 4, 888–893 (2018).
33
H. F. Son et al., Rational protein engineering of thermo-stable PETase from Ideonella sakaiensis for highly efficient PET degradation. ACS Catal. 9, 3519–3526 (2019).
34
A. Rauwerdink, R. J. Kazlauskas, How the same core catalytic machinery catalyzes 17 different reactions: The serine-histidine-aspartate catalytic triad of α/β-hydrolase fold enzymes. ACS Catal. 5, 6153–6176 (2015).
35
G. J. Palm et al., Structure of the plastic-degrading Ideonella sakaiensis MHETase bound to a substrate. Nat. Commun. 10, 1717 (2019).
36
H.-Y. Sagong et al., Decomposition of PET film by MHETase using Exo-PETase function. ACS Catal. (2020) in press.
37
J. A. Hermoso et al., The crystal structure of feruloyl esterase A from Aspergillus niger suggests evolutive functional convergence in feruloyl esterase family. J. Mol. Biol. 338, 495–506 (2004).
38
K. Suzuki et al., Crystal structure of a feruloyl esterase belonging to the tannase family: A disulfide bond near a catalytic triad. Proteins 82, 2857–2867 (2014).
39
B. R. Brooks et al., CHARMM: The biomolecular simulation program. J. Comput. Chem. 30, 1545–1614 (2009).
40
J. C. Phillips et al., Scalable molecular dynamics with NAMD. J. Comput. Chem. 26, 1781–1802 (2005).
41
R. B. Best et al., Optimization of the additive CHARMM all-atom protein force field targeting improved sampling of the backbone φ, ψ and side-chain χ(1) and χ(2) dihedral angles. J. Chem. Theory Comput. 8, 3257–3273 (2012).
42
D. A. Case et al., The Amber biomolecular simulation programs. J. Comput. Chem. 26, 1668–1688 (2005).
43
M. Gaus, A. Goez, M. Elstner, Parametrization and benchmark of DFTB3 for organic molecules. J. Chem. Theory Comput. 9, 338–354 (2013).
44
B. C. Knott et al., The mechanism of cellulose hydrolysis by a two-step, retaining cellobiohydrolase elucidated by structural and transition path sampling studies. J. Am. Chem. Soc. 136, 321–329 (2014).
45
H. B. Mayes et al., Who’s on base? Revealing the catalytic mechanism of inverting family 6 glycoside hydrolases. Chem. Sci. (Camb.) 7, 5955–5968 (2016).
46
J. Liu, A. Hamza, C.-G. Zhan, Fundamental reaction mechanism and free energy profile for (-)-cocaine hydrolysis catalyzed by cocaine esterase. J. Am. Chem. Soc. 131, 11964–11975 (2009).
47
A. J. T. Smith et al., Structural reorganization and preorganization in enzyme active sites: Comparisons of experimental and theoretically ideal active site geometries in the multistep serine esterase reaction cycle. J. Am. Chem. Soc. 130, 15361–15373 (2008).
48
Y. Zhou, S. Wang, Y. Zhang, Catalytic reaction mechanism of acetylcholinesterase determined by Born-Oppenheimer ab initio QM/MM molecular dynamics simulations. J. Phys. Chem. B 114, 8817–8825 (2010).
49
A. L. Mitchell et al., InterPro in 2019: Improving coverage, classification and access to protein sequence annotations. Nucleic Acids Res. 47, D351–D360 (2019).
50
S. F. Altschul et al., Gapped BLAST and PSI-BLAST: A new generation of protein database search programs. Nucleic Acids Res. 25, 3389–3402 (1997).
51
K. D. Narayan, S. K. Pandey, S. K. Das, Characterization of Comamonas thiooxidans sp. nov., and comparison of thiosulfate oxidation with Comamonas testosteroni and Comamonas composti. Curr. Microbiol. 61, 248–253 (2010).
52
S. R. Eddy, Profile hidden Markov models. Bioinformatics 14, 755–763 (1998).
53
V. F. Crepin, C. B. Faulds, I. F. Connerton, A non-modular type B feruloyl esterase from Neurospora crassa exhibits concentration-dependent substrate inhibition. Biochem. J. 370, 417–427 (2003).
54
A. Waterhouse et al., SWISS-MODEL: Homology modelling of protein structures and complexes. Nucleic Acids Res. 46, W296–W303 (2018).
55
K.-E. Jaeger et al., Bacterial lipases. FEMS Microbiol. Rev. 15, 29–63 (1994).
56
X. Chen, J. L. Zaro, W.-C. Shen, Fusion protein linkers: Property, design and functionality. Adv. Drug Deliv. Rev. 65, 1357–1369 (2013).
57
A. Gandini, A. J. Silvestre, C. P. Neto, A. F. Sousa, M. Gomes, The furan counterpart of poly (ethylene terephthalate): An alternative material based on renewable resources. J. Polym. Sci. A Polym. Chem. 47, 295–298 (2009).
58
M. Barth et al., A dual enzyme system composed of a polyester hydrolase and a carboxylesterase enhances the biocatalytic degradation of polyethylene terephthalate films. Biotechnol. J. 11, 1082–1087 (2016).
59
A. Carniel, É. Valoni, J. Nicomedes, A. C. Gomes, A. M. d. Castro, Lipase from Candida antarctica (CALB) and cutinase from Humicola insolens act synergistically for PET hydrolysis to terephthalic acid. Process Biochem. 59, 84–90 (2017).
60
M. Sasoh et al., Characterization of the terephthalate degradation genes of Comamonas sp. strain E6. Appl. Environ. Microbiol. 72, 1825–1832 (2006).
61
M. Hosaka et al., Novel tripartite aromatic acid transporter essential for terephthalate uptake in Comamonas sp. strain E6. Appl. Environ. Microbiol. 79, 6148–6155 (2013).
62
T. Shigematsu, K. Yumihara, Y. Ueda, S. Morimura, K. Kida, Purification and gene cloning of the oxygenase component of the terephthalate 1,2-dioxygenase system from Delftia tsuruhatensis strain T7. FEMS Microbiol. Lett. 220, 255–260 (2003).
63
P. S. G. Chain et al., Burkholderia xenovorans LB400 harbors a multi-replicon, 9.73-Mbp genome shaped for versatility. Proc. Natl. Acad. Sci. U.S.A. 103, 15280–15287 (2006).
64
H. Hara, L. D. Eltis, J. E. Davies, W. W. Mohn, Transcriptomic analysis reveals a bifurcated terephthalate degradation pathway in Rhodococcus sp. strain RHA1. J. Bacteriol. 189, 1641–1647 (2007).
65
K. Y. Choi et al., Molecular and biochemical analysis of phthalate and terephthalate degradation by Rhodococcus sp. strain DK17. FEMS Microbiol. Lett. 252, 207–213 (2005).
66
G. Fuchs, M. Boll, J. Heider, Microbial degradation of aromatic compounds—From one strategy to four. Nat. Rev. Microbiol. 9, 803–816 (2011).
67
T. Tiso et al., Bio-upcycling of polyethylene terephthalate. bioRxiv:10.1101/2020.03.16.993592 (18 March 2020).
68
C. M. Payne et al., Fungal cellulases. Chem. Rev. 115, 1308–1448 (2015).
69
V. G. Eijsink, G. Vaaje-Kolstad, K. M. Vårum, S. J. Horn, Towards new enzymes for biofuels: Lessons from chitinase research. Trends Biotechnol. 26, 228–235 (2008).
70
F. Quartinello et al., Highly selective enzymatic recovery of building blocks from wool-cotton-polyester textile waste blends. Polymers (Basel) 10, 1107 (2018).
71
A. Magnin, E. Pollet, V. Phalip, L. Avérous, Evaluation of biological degradation of polyurethanes. Biotechnol. Adv. 39, 107457 (2020).
72
R. C. Walker, M. F. Crowley, D. A. Case, The implementation of a fast and accurate QM/MM potential method in Amber. J. Comput. Chem. 29, 1019–1031 (2008).
73
S. Kumar, G. Stecher, K. Tamura, MEGA7: Molecular evolutionary genetics analysis version 7.0 for bigger datasets. Mol. Biol. Evol. 33, 1870–1874 (2016).
Information & Authors
Information
Published in
Classifications
Copyright
Copyright © 2020 the Author(s). Published by PNAS. This open access article is distributed under Creative Commons Attribution License 4.0 (CC BY).
Data Availability
The atomic coordinates and structure factors have been deposited in the Protein Data Bank, https://www.wwpdb.org/ (PDB ID codes 6QZ1, 6QZ2, 6QZ3, and 6QZ4).
Submission history
Published online: September 28, 2020
Published in issue: October 13, 2020
Keywords
Acknowledgments
We thank staff at the Diamond Light Source for beamtime and support, and McKenzie Hawkins for technical support with HPLC analysis. This work was authored in part by Alliance for Sustainable Energy, the manager and operator of the National Renewable Energy Laboratory (NREL) for the US Department of Energy (DOE) under Contract DE-AC36-08GO28308. Funding was provided by the DOE, Office of Energy Efficiency and Renewable Energy, Advanced Manufacturing Office and Bioenergy Technologies Office. This work was performed as part of the Bio-Optimized Technologies to keep Thermoplastics out of Landfills and the Environment Consortium and was supported by the Advanced Manufacturing Office and Bioenergy Technologies Office under contract DE-AC36-08GO28308 with the NREL, operated by Alliance for Sustainable Energy, LLC. G.T.B., B.C.K., E.E., B.S.D., N.A.R., G.D., and C.W.J. also acknowledge the NREL Laboratory Directed Research and Development program for funding. J.E.M. acknowledges Research England for E3 funding and the Biotechnology and Biological Sciences Research Council for Grant BB/P011918/1. R.G. and H.P.A. were funded through an NREL subcontract and University of Portsmouth Faculty of Science bursary. F.L.K. acknowledges support from the NSF Graduate Research Fellowship Program (Grant 3900101301). This work is also based on work supported by (while C.M.P. is serving at) the NSF. Any opinions, findings and conclusions or recommendations expressed in this material are those of the authors and do not necessarily reflect the views of the NSF. H.L.W. acknowledges the DOE for Grant DE-SC0011297TDD, the NIH for Grant 1R01GM129519-01, and the NSF for Grant CHE-1464946. J.E.G. acknowledges support from NSF Grant CBET-1552355 (to C.M.P.). V.C. and E.T. acknowledge the NSF for Grant MCB-1714556. Computer time was provided by Extreme Science and Engineering Discovery Environment allocation MCB-090159 at the San Diego Supercomputing Center and the Texas Advanced Computing Center, and by the National Renewable Energy Laboratory Computational Sciences Center supported by the DOE Office of Energy Efficience and Renewable Energy under Contract DE-AC36-08GO28308.
Notes
This article is a PNAS Direct Submission.
Authors
Competing Interests
Competing interest statement: A patent application was filed on this work.
Metrics & Citations
Metrics
Citation statements
Altmetrics
Citations
Cite this article
117 (41) 25476-25485,
Export the article citation data by selecting a format from the list below and clicking Export.
Cited by
Loading...
View Options
View options
PDF format
Download this article as a PDF file
DOWNLOAD PDFLogin options
Check if you have access through your login credentials or your institution to get full access on this article.
Personal login Institutional LoginRecommend to a librarian
Recommend PNAS to a LibrarianPurchase options
Purchase this article to access the full text.