Discovery and rational engineering of PET hydrolase with both mesophilic and thermophilic PET hydrolase properties

Jul 09, 2023

Nature Communications volume 14, Article number: 4556 (2023) Cite this article

1678 Accesses

2 Altmetric

Metrics details

Excessive polyethylene terephthalate (PET) waste causes a variety of problems. Extensive research focused on the development of superior PET hydrolases for PET biorecycling has been conducted. However, template enzymes employed in enzyme engineering mainly focused on IsPETase and leaf-branch compost cutinase, which exhibit mesophilic and thermophilic hydrolytic properties, respectively. Herein, we report a PET hydrolase from Cryptosporangium aurantiacum (CaPETase) that exhibits high thermostability and remarkable PET degradation activity at ambient temperatures. We uncover the crystal structure of CaPETase, which displays a distinct backbone conformation at the active site and residues forming the substrate binding cleft, compared with other PET hydrolases. We further develop a CaPETaseM9 variant that exhibits robust thermostability with a Tm of 83.2 °C and 41.7-fold enhanced PET hydrolytic activity at 60 °C compared with CaPETaseWT. CaPETaseM9 almost completely decompose both transparent and colored post-consumer PET powder at 55 °C within half a day in a pH-stat bioreactor.

Plastics are synthetic polymers, with advantages of chemical resistance, light weight, and low cost. They have been extensively used worldwide since the 1950s and have become essential in our daily lives1,2. However, because plastics do not naturally decompose, billions of tons of plastic waste in landfills, floating in waste islands in the ocean, and dispersed as microplastics have led to severe global pollution3,4,5,6,7. Governments and plastic manufacturers around the world are aware of these issues, and research and development focused on chemical and biological recycling strategies for plastic waste have accelerated in recent years8,9,10.

Polyethylene terephthalate (PET), a polyester composed of units of terephthalic acid (TPA) and ethylene glycol, is a widely used plastic packaging material and is relatively easy to recycle. Studies on the biological degradation of PET have been actively conducted in the last two decades8,11,12,13. PET biorecycling technology includes hydrolysis of PET by microorganisms or enzymes and its conversion into high-value-added chemicals, which ultimately establishes a circular economy for PET11,14,15,16.

To date, various enzymes capable of hydrolyzing PET have been discovered and characterized biochemically and structurally. In particular, PETase from Ideonella sakaiensis 201-F6 (IsPETase) exhibits the highest PET hydrolytic activity at ambient temperatures and is considered a promising enzyme for PET waste treatment17. Accordingly, various enzyme engineering studies are currently underway to improve the PET hydrolysis efficiency and thermal stability of IsPETase, and variants with improved performance have been reported18,19,20,21,22,23,24. Recently, metagenome-derived leaf-branch compost cutinase (LCC), which exhibits thermophilic properties, has been discovered25. Its variant has shown to exhibit high PET depolymerization activity in a bioreactor system at 72 °C, rendering it a candidate for biological PET recycling26. Moreover, efforts for discovering new PET-degrading enzymes are being made, and enzymes from various microorganisms and environmental metagenomes, such as Thermobifida fusca cutinases (TfCut1,2)27, Thermomonospora curvata DSM43183 cutinases (Tcur1278,0390)28, Bacterium HR29 (BhrPETase)29, Fusarium solani pisi (FsCut)30, Humicola insolens (HiC)31,32, PET233, PET533, PET633, and PHL734 have been reported. And recently, thermotolerant PET hydrolases were discovered by genome mining in bioinformatics35.

For feasible application of enzymatic PET hydrolysis, the inherently high catalytic activity of PET hydrolases is crucial factor as a starting point. Also, exploiting the thermophilic PET hydrolase properties is also an efficient strategy for development of superior PET hydrolases, because high-temperature operation near the glass transition temperature of PET material is favorable to the PET degradation performance36. Thus, the discovery of promising PET-degrading enzymes that have both high catalytic activity and high thermostability properties is strongly needed to provide insight into the catalysis reaction and expand the scope of efficient PET biocatalysts.

There is a rapid progress in research on developing efficient PET-degrading enzymes, and numerous variants have emerged. Nonetheless, recently developed PET hydrolases do not outperform the previously reported enzymes in terms of enzymatic catalysis and thermostability. Thus, template enzymes used for enzyme engineering are mainly focused on IsPETase and LCC, which have mesophilic and thermophilic hydrolytic properties, respectively. In the present study, we report a PET hydrolase from Cryptosporangium aurantiacum (CaPETase) that exhibits remarkable PET degradation activity at ambient temperatures and high thermostability. The crystal structure of CaPETase indicates that the enzyme has a unique active site compared with other known PET hydrolases. Through structure-guided rational protein engineering, we developed a robust CaPETase variant that exhibits markedly enhanced enzyme activity and thermostability. The applicability of this variant in the recycling industry was further demonstrated by evaluating its activity with a pH-stat bioreactor.

To discover a PET hydrolase, we performed a sequence homology analysis using the National Center for Biotechnology Information (NCBI) database and selected 10 PETase candidates (see “Methods” for details). A phylogenetic tree was constructed for the 10 selected PETase candidates and 17 reported PET hydrolases. The 27 enzymes were divided into two groups: one group contained mesophilic enzymes, such as IsPETase, and the other group contained thermophilic enzymes, such as TfCut2 and LCC (Fig. 1a, Supplementary Fig. 1 and Supplementary Table 1). The phylogenetic tree was further separated into 10 subgroups, and 3 selected PETase candidates, namely, RZL00883.1, KOX11336.1, and SHM40309.1, formed discrete lineages with low phylogenetic relationships to the reported PET hydrolases (Fig. 1a). To characterize the 10 selected PETase candidates (PCs), we first attempted to produce them in a signal peptide-truncated form. Eight of the 10 PETase candidates were successfully produced, except for KOX11336.1 and MAM88718.1. We also measured the melting temperature (Tm) of the eight PETase candidates to determine the thermostability of these enzymes. The candidates exhibited a range of Tm values from 38.6 °C to 70.5 °C (Fig. 1b and Supplementary Fig. 2). We then checked PET hydrolytic activity of the eight PETase candidates at a wide range of temperatures from 30 to 60 °C using several PET samples, such as post-consumer transparent PET powder (PC-PETTransparent), semi-crystalline PET powder (Cry-PET, Goodfellow Cambridge Ltd) (Cat. No. ES306000), and amorphous PET film (AF-PET, Goodfellow Cambridge Ltd) (Cat. No. ES301445). Among these candidates, PC3, PC7, PC8 and PC10 showed relatively low or undetectable levels of PET hydrolytic activity across the tested conditions (Fig. 1c and Supplementary Fig. 3). PC6, which have the highest Tm value of 70.5 °C, produced only negligible amounts of PET hydrolysis products at 30 °C, but showed the optimal PET hydrolytic activity at 50 °C (Fig. 1b, c and Supplementary Figs. 2, 3). PC4 and PC5 showed relatively high PET hydrolytic activity at 50 °C and 40 °C, respectively (Fig. 1c and Supplementary Fig. 3). Surprisingly, compared to other PETase candidates, PC2 exhibited significantly high PET hydrolytic activity across a broad range of reaction conditions. Notably, PC2 exhibited superior activity at 30 °C and produced the highest amount of PET hydrolysis products across all three PET substrates conditions compared to other candidates. (Fig. 1c and Supplementary Fig. 3). In addition, the pH profile results for the eight PETase candidates also showed that PC2 showed the highest level of PET hydrolytic activity (Supplementary Fig. 4). It is noteworthy that PC2 showed remarkable PET hydrolytic activity compared with the other enzymes, and also exhibited high thermostability with a Tm value of 66.8 °C (Fig. 1b, c and Supplementary Fig. 2). Moreover, PC2 had the highest soluble expression level compared with the other enzymes (Fig. 1b). These results indicate that PC2 has excellent properties for efficient PET degradation, including enzyme activity, thermostability, and protein expression levels. Thus, we selected PC2 (accession code: SHM40309.1, PETase from Cryptosporangium aurantiacum, CaPETase) as the most robust PET hydrolase among the eight PETase candidates tested. The measurements of changes of the Tm value and activity by addition of metal ions showed that PC2 is not a metal ion-dependent enzyme (Supplementary Fig. 5).

a Maximum likelihood phylogenetic tree and percentage identity matrix of the 10 selected PETase candidate (PC1–PC10) and 17 reported PET hydrolase sequences. Bootstrap values for 1000 replications are shown at the branching edges. The colored bar represents the level of the percent identity of the enzymes, and detailed percent identity values are listed in Supplementary Fig. 31. b Protein yield and the Tm values of the 8 PCs. c PET hydrolytic activity of the eight PCs. The reaction was performed with post-consumer transparent PET powder (PC-PETTransparent, 15 mg mL−1 with 500 nM enzyme), semi-crystalline PET powder (Cry-PET, 15 mg mL−1 with 2 µM enzyme), and amorphous PET film (AF-PET, 15 mg mL−1 with 2 µM enzyme) in 50 mM Glycine-NaOH pH 9.0 buffer at various temperatures (30 °C, 40 °C, 50 °C, 60 °C) for 3 days. Reactions were performed in triplicate; Data are presented as mean values  ±  SD. d Comparison of the PET hydrolytic activity of CaPETase, IsPETase, LCC, and TfCut2. The reaction was conducted with Cry-PET (15 mg mL−1 with 2 µM enzyme) in 50 mM Glycine-NaOH (pH 9.0) at various temperatures (30 °C, 40 °C, 50 °C, 60 °C) for 12 h. Reactions were performed in triplicate; Data are presented as mean values  ±  SD.

Next, we compared the PET hydrolytic activity of CaPETase with that of well-known PET hydrolases, such as IsPETase, TfCut2, and LCC, over a broad temperature range from 30 to 60 °C using Cry-PET as a substrate. In reactions at 30 °C, CaPETase showed significantly higher PET hydrolytic activity than LCC and TfCut2 (Fig. 1d and Supplementary Fig. 3). Moreover, CaPETase exhibited 1.4-fold higher activity than IsPETase, which is known to have the highest PET hydrolytic activity at ambient temperature among the reported PET hydrolases (Fig. 1d)17. In particular, the PET hydrolytic activity of CaPETase was 3.1-fold higher than that of IsPETase at 40 °C (Fig. 1d), likely because CaPETase has much higher thermostability than IsPETase. However, the PET hydrolytic activity of CaPETase at 60 °C dramatically decreased and reversed compared with that of LCC at temperatures of 50 °C and 60 °C (Fig. 1d). These results indicate that CaPETase is a promising PET hydrolase that exhibits high PET decomposition ability and thermostability. Considering that it is important to make improvements without the loss of enzyme activity and thermostability in the development of superior PET-degrading enzymes37, we propose that CaPETase represents a more efficient template enzyme for enzyme engineering than other enzymes with extreme mesophilic and thermophilic properties, such as IsPETase and LCC, respectively.

To provide a structural basis for high PET hydrolytic activity of CaPETase, we determined its crystal structure at a resolution of 1.36 Å (Supplementary Table 2). CaPETase shows an α/β hydrolase fold and a nine-stranded β-sheet at the center surrounded by six α-helices and two 310-helices (Fig. 2a and Supplementary Fig. 6). Sequence-independent pairwise superposition of CaPETase with three distinctive PET hydrolases, namely, IsPETase, LCC, and TfCut2, generated global root mean square deviation values of 0.69, 0.65, and 0.53 Å, respectively. Formation of one conserved disulfide bond (DS, C279/C297) and lack of an extended loop in the substrate binding site of CaPETase suggest that the enzyme originated from an ancestor of TfCut2 and LCC rather than IsPETase (Supplementary Fig. 7). Interestingly, CaPETase exhibits a somewhat different backbone structure at the active site compared with other PET hydrolases (Supplementary Fig. 8). Because the structural comparison using a Cartesian coordinate system is known to be subjective for distinguishing the detailed structural differences of the main chains38, we further analyzed the φ-ψ torsion angles of the main chains of these four PETases (Supplementary Fig. 9) and found that there were local differences in backbone torsion angles between CaPETase and other PET hydrolases (Fig. 2a and Supplementary Fig. 10). Interestingly, CaPETase also showed significant differences in the backbone torsion angles at the five connecting loops (β3–α1, β4–α2, β6–β7, β7–α4, and β8–α5) that form an active site, whereas comparisons of the corresponding loops between the other three PET hydrolases exhibited less differences (Fig. 2a and Supplementary Fig. 11), suggesting that CaPETase has a unique active site conformation. There were some differences in the network of residues extending from the active site to the nearby spatial environment compared with that of the other PET hydrolases. Among them, we observed unique differences affecting the backbone torsion angles of these loops. Near the β3–α1 loop, distinct residues positioned in the β4–α2 loop and a W105–L108–G124 network force, which form a unique side-chain internal network, appear to influence the conformation of the β3–α1 loop (Supplementary Fig. 12). In fact, the β3–α1 loop has high root mean square fluctuation values near the active sites of other PET hydrolases in molecular dynamic simulations39,40. A unique A192–G212–F248 network is formed under the β8–α5 loop, where catalytic H246 is located (Supplementary Fig. 13). At the corresponding F248 position in CaPETase, LCC and TfCut2 have an alanine residue, whereas IsPETase has a cysteine residue that forms a second disulfide bond. Therefore, the positioning of a bulky F248 might cause significant torsional differences in the β8–α5 loop and β7–α4 loop of CaPETase (Supplementary Fig. 13). Finally, an R176–W200–F209 network appears to trigger conformational differences in the α3-helix and β6–β7 loop (Supplementary Fig. 14). Importantly, the α3-helix contains the catalytic S169, and the β6–β7 loop was previously annotated as a wobbling tryptophan-containing loop in PETase from Rhizobacter gummiphilus (Supplementary Fig. 14)41. We further analyzed backbone fluctuations of these four PET hydrolases using molecular dynamic simulations and CaPETase exhibits quite unique backbone fluctuation profile (Fig. 2b and Supplementary Fig. 15). CaPETase has more stable β6–β7 and β7–α4 loops than mesophilic IsPETase, and particularly, the enzyme shows high stability at the front region of the β8–α5 catalytic loop where H246 is located (Fig. 2b). However, the end region of β8–α5 which corresponds to the extended loop of IsPETase, and the front region of β3–α1 loop showed the highest and lowest flexibility among the four homologs, respectively (Fig. 2b). To our interest, the differences of the backbone fluctuation profile was localized exactly to the unique internal network affecting the backbone torsion angles of these loops. Thus, we believe that the unique backbone conformation of CaPETase allows the enzyme to maintain high activity while stabilizing several flexible loops of the mesophilic PET hydrolase.

a Comparison of the backbone torsion angle differences between CaPETase and IsPETase, LCC, and TfCut2. The structure of five connecting loops forming the active site of CaPETase is displayed as a putty tube representation of the same diameter in PyMoL. The structure is colored according to the Euclidean distance values between the two Ramachandran points of the aligned residues. Colors of white to red designate low to high Euclidean distance values, respectively. The catalytic triad of CaPETase is shown as a stick model with a cyan-color circle. b MD simulations show unique backbone fluctuation profile of CaPETase. Cα atom root-mean-square fluctuations (RMSF, Å) of the CaPETase, IsPETase, TfCut2, LCC during MD simulations. c Comparison of the residues forming the substrate binding cleft of CaPETase, LCC, and TfCut2. The highlighted residues are shown as a stick model. d Distinct residues in the substrate binding site of CaPETase. Distinct and conserved residues are presented in magenta and light blue, respectively. e PET hydrolytic activity of the variants. PC-PETTransparent (15 mg mL−1) were incubated with 500 nM enzymes at 40 °C for 24 h in 50 mM Glycine-NaOH buffer pH 9.0. Total amount of released products and the Tm value of the variants are shown as bars and red-colored dots, respectively. Reactions were performed in triplicate; Data are presented as mean values  ±  SD.

In addition to the unique backbone conformation at the active site, residues forming the substrate binding cleft of CaPETase showed significant differences compared with other thermophilic PET hydrolases (Fig. 2c, d). In the vicinity of the wobbling W194, CaPETase possesses unique G196 and L133 residues, where highly conserved residues are located in other PET hydrolases (Fig. 2c, d and Supplementary Fig. 7). Mutating these residues to the corresponding residues in other PET hydrolases, such as G196L, L133Y, and L133Q, had a negative effect on enzyme activity and/or thermostability (Fig. 2e). However, the G196T mutation exhibited enhanced thermostability (Fig. 2e), which may result from the formation of a hydrogen bond between G196T and N195. CaPETase also contains a unique I102 residue in the β3–α1 loop showing the largest torsion differences, whereas other PET hydrolases contain a highly conserved threonine residue at the corresponding position, which probably enables CaPETase to form a relatively wider substrate binding cleft (Fig. 2c, d and Supplementary Fig. 16). Replacement of I102 with threonine resulted in decreased enzyme activity, confirming that the residue contributes to high enzyme activity (Fig. 2e). Furthermore, CaPETase has unique residues, such as Q107, W168, and T250, at the regions of the β3–α1 loop, β8–α5 loop, and α3, whereas most of the corresponding residues are highly conserved in other thermophilic PET hydrolases (Fig. 2c, d and Supplementary Fig. 7). Mutating these residues to the conserved residues in other thermophilic PET hydrolases decreased enzymatic activity and/or stability, indicating that the combined positioning of these residues is necessary to create an optimal substrate binding site for CaPETase with a unique shape and polarity (Fig. 2e). One exception was the Q107S mutation, which resulted in no noticeable differences in enzyme activity or thermostability (Fig. 2e). Taken together, we suggest that along with unique backbone torsion angles, the positioning of distinct residues at the substrate binding site enable CaPETase to form an optimal substrate binding site for high PET hydrolytic activity.

Although CaPETase has high PET hydrolytic activity and thermostability, its performance is still insufficient for industrial applications. We conducted rational protein engineering of CaPETase to further enhance the PET hydrolytic activity and thermostability of the enzyme using various strategies, such as introducing disulfide bonds and hydrogen bonds and modifying the protein surface charge (Supplementary Fig. 17). The thermostability of the variants was monitored by measuring Tm values, and the PET hydrolytic activity of the variants was measured using post-consumer transparent PET powder (PC-PETTransparent) at ambient temperature (40 °C). We introduced four disulfide bonds, namely, G76C/A143C (DS1), L180C/A202C (DS2), T204C/R233C (DS3), and R242C/S291C (DS4). The DS2 and DS4 mutations increased the Tm value by approximately 3 °C, whereas the DS1 and DS3 mutations decreased the Tm value compared with CaPETaseWT (Fig. 3a and Supplementary Fig. 18). Moreover, the introduction of the DS2 and DS4 mutations increased PET hydrolytic activity by more than 20% compared with CaPETaseWT (Fig. 3a and Supplementary Fig. 18). These results indicate that the DS2 and DS4 mutations were successfully formed in CaPETaseWT and exerted positive effects on enzyme activity and thermostability. We also attempted to improve the thermostability of CaPETase by introducing noncovalent bonds, such as hydrogen bonds and salt bridges, and designed seven mutations, namely, V129T (NC1), P136S (NC2), A192T (NC3), R198K (NC4), V203T (NC5), A252N (NC6), and A257S(NC7). Of these, the NC1 and NC4 mutations increased Tm values by approximately 2 °C and enhanced PET hydrolytic activity by 30% compared with CaPETaseWT (Fig. 3a and Supplementary Fig. 18). Finally, in an attempt to improve the protein adsorption ability to the PET surface by modifying the protein surface charge, we designed five mutations to render the protein surface hydrophobic, i.e., N109A (HP1), R151A (HP2), R157A (HP3), R160A (HP4), and R233A (HP5), and four mutations to render the protein surface positive, i.e., T86R (SC1), A155R (SC2), T275R (SC3), and M294R (SC4). Unfortunately, most mutations did not show significant changes or even negative effects on thermostability or enzyme activity; however, the HP1 mutation increased the Tm value by 3.2 °C, and the SC2 mutation enhanced PET hydrolytic activity by 20% compared with CaPETaseWT (Fig. 3a and Supplementary Fig. 18). Taken together, we introduced eight point-mutations that resulted in improved thermostability and PET hydrolytic activity, i.e., DS2, DS4, NC1, NC4, HP1, and SC2, among the 20 rationally designed mutations tested (Fig. 3a and Supplementary Fig. 18). There were also ambiguous mutations that only improved enzyme activity or thermostability, such as NC2, HP2, and HP5. We excluded these mutations from the final selection to develop a much superior variant without compromising enzymatic activity or thermostability (Fig. 3a and Supplementary Fig. 18)37.

a Single-point mutations of CaPETase. Released PET hydrolysis products and the Tm values of the single-point mutations are presented. The reactions were performed using 500 nM enzymes with post-consumer transparent PET powder (PC-PETTransparent, 15 mg mL−1) as the substrate in 50 mM Glycine-NaOH buffer (pH 9.0) for 24 h at 40 °C. Reaction was carried out in triplicate; error bars represent the s.d. of the replicate measurement. b Combinatorial mutations of CaPETase. PET hydrolytic activity of the variants generated using the combinatorial strategy. Released PET hydrolysis products per hour and the Tm values of the combinatorial variants are presented. The reactions were performed using 500 nM enzymes with PC-PET (15 mg mL−1) as the substrate in 100 mM Glycine-NaOH buffer (pH 9.0) at 40 °C for 24 h and 60 °C for 6 h, respectively. Reaction was carried out in triplicate; error bars represent the s.d. of the replicate measurement. c Comparison of PET hydrolysis activity between CaPETaseM9 and LCCICCG at various temperatures. The reactions were carried out at different temperatures using PC-PETTransparent (12.5 mg mL−1) with 1 µM enzyme and Cry-PET (12.5 mg mL−1) with 4 µM enzyme under the 200 mM Glycine-NaOH buffer pH 9.0. Reactions were performed in triplicate; Data are presented as mean values  ±  SD.

We sequentially integrated the six mutations described above to develop a superior CaPETase variant with higher thermostability and PET hydrolytic activity. First, we combined the DS2 and DS4 mutations, and the resulting CaPETaseDS2/DS4 variant showed a synergistic effect on thermostability with a Tm value of 74.3 °C (ΔTm = 7.4 °C) (Fig. 3b and Supplementary Fig. 19). Moreover, the variant enhanced PET hydrolytic activity by 1.35- and 4.45-fold at 40 °C and 60 °C, respectively, compared with CaPETaseWT (Fig. 3b and Supplementary Fig. 19). Next, we set up CaPETaseDS2/DS4 as a scaffold for the next combination. We integrated the NC1/NC4, HP1, and SC2 mutations individually into CaPETaseDS2/DS4 using our engineering strategy. The addition of the NC1/NC4 mutation increased the Tm value significantly by 3.9 °C and increased PET hydrolytic activity at both 40 °C and 60 °C (Fig. 3b and Supplementary Fig. 19). It showed 3.8- and 17-fold enhanced activity at 60 °C compared with CaPETaseDS2/DS4 and CaPETaseWT, respectively (Fig. 3b and Supplementary Fig. 19). The addition of the HP1 mutation increased the Tm value by 2.7 °C and enhanced PET hydrolytic activity by 1.7-fold at 60 °C compared with CaPETaseDS2/DS4 (Fig. 3b and Supplementary Fig. 19). When the SC2 mutation was integrated into the CaPETaseDS2/DS4 variant, we observed no noticeable improvements in activity and thermostability; however, there was a slight increase in PET hydrolytic activity at 60 °C (Fig. 3b and Supplementary Fig. 19).

These results suggest that all four mutations (NC1, NC4, HP1, and SC2) had a positive effect on the thermostability and activity of CaPETaseDS2/DS4; thus, we combined the four mutations into CaPETaseDS2/DS4 to generate CaPETaseDS2/DS4/NC1/NC4/HP1/SC2 (CaPETaseM8). Surprisingly, when all four mutations were added to CaPETaseDS2/DS4, a synergistic effect on thermostability and enzyme activity was observed, and CaPETaseM8 exhibited significantly enhanced thermostability with a Tm value of 80.7 °C and 1.5- and 25.8-fold enhanced PET hydrolytic activity at 40 °C and 60 °C, respectively, compared with CaPETaseWT (Fig. 3b and Supplementary Fig. 19).

As mentioned above, the G196T mutation resulted in positive effects on both thermostability and enzyme activity (Fig. 2e); thus, we finally generated CaPETaseDS2/DS4/NC1/NC4/HP1/SC2/G196T (CaPETaseM9) by integrating the G196T mutation into CaPETaseM8. CaPETaseM9 exhibited a Tm value of 83.2 °C, which corresponds to a 16.7 °C increase in Tm compared with CaPETaseWT. Moreover, the PET hydrolytic activity of CaPETaseM9 increased by 1.7- and 31.2-fold at 40 °C and 60 °C, respectively, compared with CaPETaseWT (Fig. 3b and Supplementary Fig. 19). These results indicate a positive effect of G196T on CaPETaseWT was applied similarly to CaPETaseM8.

CaPETaseM9 showed much higher activity at all temperature conditions from 30 to 70 °C, and particularly, showed 41.7-fold higher specific activity at 60 °C than CaPETaseWT (Supplementary Fig. 20). The result was also reproduced in a scale-up system of 50-mL shaking flasks (Supplementary Fig. 21). These results demonstrated the improved enzyme activity and reinforced thermostability of CaPETaseM9. The improved thermostability of the variant was further verified through heat inactivation experiments, where CaPETaseM9 maintained its activity even after incubation at 60 °C for 12 h, whereas CaPETaseWT showed complete loss of activity within an hour (Supplementary Fig. 22).

We then compared the PET hydrolytic activity of CaPETaseM9 with LCCICCG towards PC-PET and Cry-PET at temperatures ranging from 30 to 60 °C. CaPETaseM9 showed significantly higher PET hydrolytic activity compared to LCCICCG, at 30 °C and 40 °C (Fig. 3c). At 50 °C and 60 °C, CaPETaseM9 showed quite similar activity compared with LCCICCG (Fig. 3c).

To provide structural insights into the enhanced PET-degrading capacity of CaPETaseM9, we determined its crystal structure at a resolution of 1.53 Å (Fig. 4 and Supplementary Table 2). The formation of the introduced DS2 and DS4 disulfide bonds was clearly observed in CaPETaseM9, and the S–S interatomic length of both disulfide bonds was within the optimal disulfide bond length range (Fig. 4 and Supplementary Fig. 23). Interestingly, DS4 was located in the vicinity of one of the calcium binding sites of Cut190 and the mutation point of IsPETase R280A42,43, and the formation of DS4 also caused significant changes in the surface electrostatic potential and neighboring region conformation (Fig. 4 and Supplementary Figs. 23 and 24). The side chain of the mutated V129T was flipped to form a hydrogen bond with the adjacent T131 and D132 residues, thereby further stabilizing the β4–α3 connecting loop (Fig. 4 and Supplementary Figs. 23 and 25). With respect to R198K, the mutated lysine residue moved inward to form hydrogen bonds with the main chains of N195 and D222, thereby stabilizing the wobbly tryptophan-containing loop (Fig. 4 and Supplementary Fig. 23). The mutated G196T formed a water-mediated hydrogen bond with the adjacent N195 residue, resulting in further stabilization of the wobbly tryptophan-containing loop (Fig. 4 and Supplementary Fig. 23). The A155R mutation changed the hydrophobic surface to a positive charge, which seems to increase the attachment of the enzyme to the PET surface, as suggested by a previous report (Fig. 4 and Supplementary Fig. 23)44. Finally, the N109A mutation appeared to strengthen internal hydrophobic interactions (Fig. 4 and Supplementary Figs. 23 and 26).

The crystal structure of CaPETaseM9 is shown as a cartoon diagram, and the mutated residues are shown as a stick or a surface electrostatic potential model.

To evaluate the industrial applicability of CaPETaseM9, we conducted a PET decomposition experiment in a pH-stat bioreactor using PC-PETTransparent as a substrate (Supplementary Fig. 27). The bioreactor was operated at 55 °C using 2.70 mgenzyme ∙ gPET−1, and the pH was continuously titrated at 8.0 by adding NaOH. The decomposition rate was measured by monitoring released amounts of MHET and TPA. After a short lag phase of an hour, which was required for initial hydrophilization, the PET degradation rate increased exponentially, and 50% of PC-PETTransparent was depolymerized within 4 h (Fig. 5a). In the second half of the reaction, the degradation rate slightly decreased because of a decrease in the amount of substrate; however, a final degradation rate of 94.1% was achieved after 12 h (Fig. 5a). This was a significant result in terms of showing that a high depolymerization rate of 90% or more could be achieved even at 55 °C, which is a temperature condition relatively lower than the Tg temperature. This result also indicated that CaPETaseM9 has significant PET hydrolytic activity and thermostability comparable to other benchmark biocatalysts. We also performed decomposition of a post-consumer colored PET powder (PC-PETColored), which is known to be relatively difficult to recycle because of the presence of colors, additives, multilayer structure, labels and other complexities45. Interestingly, the depolymerization rate of PC-PETColored was almost identical to that of PC-PETTransparent, showing 50% depolymerization within 4 h (Fig. 5b); however, the final depolymerization rate of PC-PETColored was 89.2% at 12 h, which was slightly lower than that of PC-PETTransparent (Fig. 5b). This is probably due to the impurities present in PC-PETColored. These results demonstrate that unlike other recycling methods, biorecycling of PET plastic can be achieved regardless of the color of PET plastic.

Decomposition of post-consumer transparent PET powder (PC-PETTransparent) (a) and post-consumer colored PET powder (PC-PETColored) (b) in a pH-stat bioreactor using CaPETaseM9. Reactions were performed in triplicate independently; Data are presented as mean values  ±  SD. c Complete degradation of a post-consumer PET container using CaPETaseM9 at 60 °C. Reactions were performed in triplicate; Data are presented as mean values  ±  SD.

Finally, we determined whether untreated post-consumer PET containers can be depolymerized by CaPETaseM9. As depolymerization proceeded, the PET film became opaque and thin, and the PET film disappeared completely in 3 days (Fig. 5c). These results suggest that CaPETaseM9 can be utilized for decomposing PET plastics with various physical properties.

As plastic waste is emerging as an environmental issue, research on enzymatic PET depolymerization for plastic waste is expanding as a sustainable alternative for plastic waste treatment and recycling. Therefore, research focused on the discovery of promising PET-degrading enzymes has been actively conducted, and to date, more than 20 PET-degrading enzymes have been biochemically identified and characterized46,47. Moreover, studies on developing robust PET-degrading enzymes through various engineering strategies, such as rational designing, directed evolution, and computation-guided approaches, have been conducted. Although reports on the development of efficient PET-degrading enzyme variants are rapidly increasing, the target enzymes for the development of superior variants are focused on a few PET hydrolases, such as IsPETase and LCC, which have excellent catalytic and thermostability properties, respectively. In this work, we reported a PET-degrading enzyme, CaPETase, with both mesophilic and thermophilic PET hydrolase properties; It exhibits higher activity than IsPETase under the ambient conditions and has a Tm value of 66.8 °C. A positioning of serine for tryptophan wobbling and an extended loop for enzyme-substrate access were proposed as key structural features of IsPETase for its high activity at ambient temperature21,22,48. Although CaPETase shows high activity at ambient temperature, the enzyme does not share these structural features of mesophilic PET hydrolase (Supplementary Fig. 7). Moreover, its backbone fluctuation profile showed that overall flexibility at the active site is lower than that of IsPETase (Fig. 2b), suggesting that it has high activity at ambient conditions through a structural feature different from that of IsPETase. Structural analysis of CaPETase with three representative PET-degrading enzymes, such as IsPETase, TfCut2, and LCC, showed that CaPETase has unique backbone torsions of the active site loops caused by three major internal networking differences. One of the distinct internal networks is located under the β8–α5 loop (Supplementary Fig. 13). It is noteworthy to note that the unique disulfide bond of IsPETase corresponds to this region and that replacements of the two cysteine residues with alanine severely decrease its BHET and PET hydrolytic activity22. Compared with LCC and TfCut2 having a A-A-T network, CaPETase has a unique G212-F248-A192 network having the bulky hydrophobic side chain ‘anchoring’ to the core region (Supplementary Fig. 13). Presumably, the anchoring phenylalanine in the network allowed CaPETase to maintain lower perturbation of the catalytic H246 position in the front region of β8–α5 loop, although the enzyme shows a relatively higher backbone fluctuation at the middle region of the β8–α5 loop (Fig. 2b, c). The middle region of the β8–α5 loop in CaPETase contains a unique T250 residue, which interacts with W168 and the β3–α1 loop (Fig. 2b, c). Formation of this unique environment extends to the unique W105-L108-G124 network at the interface between β3–α1 and β4–α2 loops (Supplementary Fig. 12). The R176-W200-F209 network of β6–β7 loop is also noteworthy, as this region has been reported to be important for the activity of mesophilic enzymes41,42,48. Taken together, the unique structural features and backbone dynamics profile of CaPETase seem to contribute to its stability and high activity at the same time.

Through the rational design of CaPETase, combinatorial scaffold engineering was performed to increase thermal stability and catalytic efficiency of enzyme, and as a result, considerable improvement was carried out at the Tm value (ΔTm = 16.7 °C) and catalytic efficacy (increased by 41.7-folds). As part of further efforts to assess the capacity of CaPETaseM9, comparisons with benchmark enzyme and pH-stat bioreactor operation elucidated that CaPETaseM9 has sufficient PET decomposition performance for industrial applications. Thus, through structure-guided rational protein engineering of CaPETase, we developed CaPETaseM9 with dramatically enhanced PET hydrolytic activity and thermostability and demonstrated that this variant has sufficient PET depolymerization ability for industrial applications. Our engineering strategy could provide insights for the improvement of catalytic activity and thermostability for other enzymes along with the rational design used in this study. The emergence of a CaPETase having both mesophilic and thermophilic PET hydrolase properties will broaden the spectrum of enzyme for enzyme engineering and accelerate the implementation of biological recycling of PET.

Post-consumer transparent PET powder (PC-PETTransparent) and post-consumer colored PET powder (PC-PETColored) used in this study was produced by grinding PET flakes purchased from a recycling company (Yoochang R&C, Republic of Korea). Prior to milling to produce PET powder, the PET flakes were washed with water and dried at a temperature of 160 °C for 6 h. The extensively dried PET flakes were fed into the twin-screw extruder with set temperatures to 260 °C in the extruder zones, 275 °C in the discharge zones at 275 °C, and 275 °C in the die plate, and an extrusion speed of 50 rpm. The filament was then sequentially cooled into two water baths at temperatures of 75 °C and 60 °C, and excess water was removed using air wipers. Then, a single hole die with a diameter of 3.0 mm was employed for pelletizing, resulting in pellets of approximately 3 mm in length. The obtained PET pellets were micronized using a cryogenic Micro Ball Mill. Particle size distribution analysis of the PET samples was conducted using a HORIBA LA 960S2 particle size analyzer in accordance with ISO 13320. The particle size distributions, crystallinities, and molecular weights (Mn, Mw, Mz) of PC-PETTransparent and PC-PETColored are presented in Supplementary Figs. 28, 29 and Supplementary Table. 3.

Sequences were retrieved from the NCBI nonredundant (nr) protein sequence database via PSI-BLAST (parameters: E-value of 1.00E-58, % identity >40%) using I. sakaiensis PETase (Accession code: GAP38373.1) to query the model. Partial genes were removed from the retrieved sequences, and then sequences of the PETase candidates were randomly selected from the dataset. Multiple sequence alignment was performed using Clustal omega49. Phylogenetic reconstruction was performed via maximum likelihood (ML) using MEGA11 with the Dayhoff w/freq. model50,51. The tree with the highest log likelihood (−10401.10) is shown. Initial trees for heuristic search were obtained automatically by applying Neighbor-Join and BioNJ algorithms to a matrix of pairwise distances estimated using the JTT model, and then the topology with a superior log likelihood value was selected. Discrete gamma distribution was used to model evolutionary rate differences among sites (5 categories [+G, parameter = 1.6961]). The rate variation model allowed for some sites to be evolutionarily invariable ([+I], 4.28% sites). This analysis involved 27 amino acid sequences, and there was a total of 444 positions in the final dataset. Bootstrap values were obtained from 1000 replicates52. Sequence alignment of enzymes used in phylogenetic tree analysis was depicted as a graphical illustration using ESPript3 (Supplementary Fig. 1)53. Pairwise identity and similarity of the proteins were calculated using the Clustal omega method, and a sequence identity matrix was created using excel software.

All proteins used in this study were codons optimized for E. coli expression. The signal peptides sequences were removed, and NdeI and XhoI restriction sites were added to the 5′ and 3′ ends of the target gene sequences, respectively. These DNA sequences were synthesized by Twist Bioscience without signal peptides sequences. The signal peptide cleavage sites for each sequence were predicted using SignalP 3.0 server. Ten PETase candidate genes as well as TfCut2 and LCC were subcloned into a pET21b expression vector, and IsPETase was subcloned into a pET15b expression vector. Site-directed mutagenesis experiments were performed via polymerase chain reaction (PCR) with the oligonucleotides listed in Supplementary Table 4. All variants were verified by sequencing.

The constructed vectors were transformed into E. coli strain Rosetta gami-B (DE3), which was grown in lysogeny broth containing 100 mg L−1 ampicillin at 37 °C until the optical density of the culture broth reached a value of 0.6 at 600 nm. After induction by the addition of 0.5 mM isopropyl β-D-1-thiogalactopyranoside, the culture was further incubated at 18 °C for 20 h. Cells were harvested via centrifugation at 4000 × g for 15 min at 20 °C. The cell pellet of the 10 PETase candidates, LCC, and TfCut2 were resuspended in lysis buffer A (40 mM Tris-HCl, pH 8.0), and the cell pellet of IsPETase was resuspended in lysis buffer B (50 mM Na2HPO4-HCl, pH 7.0). Cell lysis was performed via ultrasonication, and the cell lysates were centrifuged at 13,500 × g for 30 min. The cell-free crude enzymes were purified using a Ni-NTA agarose column (Qiagen). After washing with buffers A and B containing 30 mM imidazole, the target proteins were eluted using 300 mM imidazole in buffers A and B. All protein purification experiments were performed at 4 °C. The purified proteins were verified via sodium dodecyl sulfate polyacrylamide gel electrophoresis and quantified using a BioTekTM Epoch microplate spectrophotometer and Gen5TM microplate data analysis software. The purified enzymes were concentrated using an Amicon Ultra Centrifugal filter unit (molecular weight cutoff value of 10,000, Millipore, Billerica, MA, USA). The expression and purification of the variants of CaPETase were performed under the same conditions as the wild-type.

The thermal stability of the enzymes was determined by establishing melting curves with a protein thermal shift dye (Applied Biosystems) using a StepOnePlus Real-time PCR instrument (Thermo Fisher Scientific). The reaction mixture contained 5 µg of proteins, 100 mM Na2HPO4-HCl (pH 7.0), and 20 µL of 1× protein thermal shift dye. Signal changes representing protein denaturation were monitored by increasing the temperature from 25 °C to 99 °C. Melting temperatures were estimated from the first derivative curve. The first derivative curves of RFU for CaPETaseWT and its variants are presented in Supplementary Fig. 30.

Crystallization of CaPETaseWT and CaPETaseM9 was initially performed with commercially available spares-matrix screens Wizard I and II (Rigaku), PEG/Ion, Index (Hampton Research), and Wizard CRYO I and II (Rigaku) using the sitting-drop vapor diffusion method at 20 °C. For each experiment, 1.0 µL of the protein solution (CaPETaseWT: 18.4 mg mL−1 in 40 mM Tris, pH 8.0, and CaPETaseM9: 22.3 mg mL−1 in 40 mM Tris, pH 8.0) was mixed with 1.0 µL of reservoir solution and the drop was equilibrated against 50 µL of reservoir solution. The crystal of CaPETaseWT appeared under conditions of 12% (v/v) PEG 3350 and 0.1 M sodium malonate (pH 4.0). The crystal of CaPETaseM9 was formed under conditions of 10% (v/v) polyethylene glycol monomethyl ether 5,000, 0.1 M HEPES, pH 7.0, and 5% tacsimate, pH 7.0. The highest quality crystals of CaPETaseWT and CaPETaseM9 were transferred into a cryoprotectant solution composed of the corresponding conditions described above with 25% (v/v) glycerol. The crystals were fished out with a loop and flash-frozen with a spout of liquid nitrogen. Data were collected using a 7 A beamline of the Pohang Accelerator Laboratory (Republic of Korea) under 100 K54. All diffraction data were indexed, integrated, and scaled using HKL2000 software55. The CaPETaseWT crystals belonged to a space group P21212 with a unit cell parameter: a = 82.31 Å, b = 82.45 Å, c = 87.39 Å, and α = β = γ = 90°. The two CaPETaseWT molecules were in an asymmetric unit, and the crystal volume per unit of protein was 2.68 Å3 Da−1, corresponding to a solvent content of 52.25%56. The structure of CaPETaseWT was determined via molecular replacement with the CCP4 version of MOLREP using the structure of cutinase 1 from Thermobifida cellulosilytica (PDB code: 5LUI) as a search model57,58. The structure model was built using the WinCoot program, and structure refinement was performed using the CCP4 REFMAC5 program59,60. The crystal of CaPETaseM9 belonged to the space group P21212 with a unit cell parameter: a = 41.02 Å, b = 112.31 Å, c = 112.54 Å, and α = β = γ = 90°. The asymmetric unit contained two molecules of CaPETaseM9, and the crystal volume per unit of protein was 2.3468 Å3 Da−1, corresponding to a solvent content of 45.37%56. The structure of CaPETaseM9 was determined via molecular replacement with the CCP4 version of MOLREP using the structure of CaPETaseWT as a search model57. Model building and refinement were performed as described above. The statistics are summarized in Supplementary Table 2. The structures of CaPETaseWT and CaPETaseM9 were deposited in the Protein Data Bank, with PDB codes 7YM9 and 7YME, respectively.

The deposited structures of 5XJH (G30-S290), 4EB0 (S36-Q293), 4CG1 (A1-F261), and 7YM9 (D43-F299) were used for IsPETase, LCC, TfCut2, and CaPETase residues and inner cavity water coordinates. All protein models were prepared with protonation states at pH 9, as calculated by propka361. The systems were prepared with ff19SB and OPC parameters for protein and explicit water and solvated in a cubic box of 70 Å length with 0.1 M Na+ and Cl- ions. We minimized the systems using gromacs with conjugate gradient algorithm, and NVT-equilibrated for 1 ns at 323.15 K, and then NPT-equilibrated for 1 ns with a force constraint of 1000 kJ∙mol−1∙nm−2. Production simulations were performed for 200 ns at 323.15 K, and rms fluctuations were calculated with 10 ps interval trajectories. All simulations were performed using Gromacs 2021.3 version package62.

In the preliminary PET degradation assay of eight PETase candidates, PC-PETTransparent (15 mg mL−1) was soaked in 50 mM Glycine-NaOH pH 9.0 buffer with 500 nM enzyme. Cry-PET (15 mg mL−1) and hole-punched amorphous PET film (⌀6 mm, roughly 10 mg) were soaked in 50 mM Glycine-NaOH pH 9.0 buffer with 2 µM enzyme. The reaction mixtures (1 mL) were then incubated at temperatures of 30 °C, 40 °C, 50 °C, and 60 °C for 3 days. The reaction was terminated by heating the mixture at 95 °C for 10 min. The reaction mixture with PET powder was filtered using a PVDF syringe filter (0.22 µm) to eliminate the residual substrate. The samples were then analyzed using HPLC. All experiments were performed in triplicate. The pH profiles on eight PETase candidates were performed using PC-PETTransparent (15 mg mL−1) in 50 mM sodium phosphate buffer (pH 6-7), 50 mM Tris-HCl buffer (pH 7-9), and 50 mM Glycine-NaOH buffer (pH 9-10) at 30 °C for 3 days. All experiments were performed in triplicate. To compare the PET hydrolytic activity of CaPETase with LCC, IsPETase, and TfCut2, 15 mg of Cry-PET was prepared and submerged in 1 mL of 50 mM glycine-NaOH pH 9.0 buffer with 2 µM enzyme. The reaction mixture was incubated at 30 °C, 40 °C, 50 °C, and 60 °C for 12 h. Then, the reaction mixture was treated as described above and analyzed using HPLC. All experiments were performed in triplicate.

The PET hydrolytic activity of CaPETaseWT and its variants was determined using PC-PETTransparent. Enzyme reactions were performed with 500 nM enzymes in 1 mL of 50 mM glycine-NaOH (pH 9.0) buffer at 40 °C for 24 h. The reaction mixture was treated as described above and analyzed using HPLC. All experiments were performed in triplicate. To evaluate the PET-hydrolyzing activity of the combinatorial variants of CaPETase, 500 nM enzyme was incubated with 100 mM glycine-NaOH (pH 9.0) and 15 mg of PC-PETTransparent at 40 °C for 24 h and 60 °C for 6 h, respectively. The reaction mixture was treated as described above and twofold diluted samples were analyzed using HPLC. All experiments were performed in triplicate.

The comparison of PET hydrolytic activities of CaPETaseWT and CaPETaseM9 were performed using PC-PETTransparent (15 mg mL−1) at 30, 40, 50, 60, 70 °C with 1 µM enzymes in 1 mL of 50 mM glycine-NaOH (pH 9.0) buffer. To further evaluate CaPETaseM9 for PET depolymerization, the enzyme (1 µM) was added to 50 mL of 200 mM Glycine-NaOH pH 9.0, containing 100 mg of PC-PETTransparent in a 250-mL flask. The flask was incubated in a shaking incubator at 40 °C, 50 °C, 60 °C and 180 rpm (40 °C for 72 h; 50 °C for 48 h; and 60 °C for 12 h). The reaction mixtures were sampled at different time points to quantify the total amount of released PET monomers. The samples were filtered as described above and 4-fold diluted samples were analyzed using HPLC. All experiments were performed in triplicate.

To compare the PET hydrolytic activity of CaPETaseM9 with LCCICCG, a membrane enclosed enzymatic catalysis (MEEC) system were used. In this system, the PET degradation reaction was induced on a cellulose dialysis membrane (Spectra/Por 2 Trial Kit, 12–14 kD, flat width 25 mm; Repligen, USA), and PC-PETTransparent and Cry-PET were used as substrate. 50 ( ± 2) mg of PC-PETTransparent and Cry-PET substrate were immersed in a 4 mL of 200 mM Glycine-NaOH pH 9.0 with 1 µM and 4 µM enzyme, respectively. Each reaction was packaged in a 4-mL dialysis bag was immersed in 200 mL of 200 mM Glycine-NaOH pH 9.0. The reaction was performed at 30 °C, 40 °C, 50 °C, and 60 °C with shaking at 80 rpm in a small shaking incubator (JEIO TECH, Republic of Korea). Sample was taken from the external equilibrium solution, diluted 5-fold, and then analyzed by HPLC. All experiments were performed in triplicate.

Protein heat inactivation experiments were performed using CaPETaseWT and CaPETaseM9. Enzyme mixtures (500 nM enzyme in 1 mL of 50 mM glycine-NaOH, pH 9.0) were incubated at 60 °C and 70 °C for 10 min to 12 h and 5 min to 12 h, respectively. Then, 15 mg of PC-PETTransparent was added to the heat-inactivated mixture, and the reaction mixture was incubated at 40 °C for 1 day. The mixture was filtered using a PVDF syringe filter (0.22 µm) and the samples were analyzed using HPLC.

A customized stirred bioreactor (MARADO-PDA, BIOCNS Co. Ltd, Korea) equipped with two standard Rushton impellers was designed for PET depolymerization. It was interfaced with a personal computer for data recording. The temperature was regulated using a refrigeration and heating bath circulator (RW3-1025, JEIO TECH, Daejeon, Republic of Korea). CaPETaseM9 (10.2 mg) was added to a 150-mL solution containing 10 mM Na2HPO4 (pH 8.0) and 3.75 g PC-PET and incubated at 55 °C and 200 rpm. The pH was regulated at 8.0 by the adding 0.3575 M NaOH (40%; extra pure grade, Duksan, Republic of Korea) solution using a PID controller. The reactor solution was sampled during the reaction and filtered through a PVDF syringe filter (0.22 µm) to remove residual PET powder. The resulting samples were diluted with 50 mM glycine-NaOH (pH 9.0) and analyzed using HPLC to quantify the amount of released PET monomer. The PET depolymerization rate was calculated by considering the sum of the hydrolysis products (TPA and MHET). All experiments were performed in triplicate.

Fragments of a post-consumer PET container (salad container) were prepared in a circular form (⌀14 mm, 50 mg). The reaction was performed using a membrane-enclosed enzymatic catalysis method. The dialysis membrane bag (Carl Roth, Product number 0653.1, MWCO 14000) was loaded with 1 µM of CaPETaseM9 dissolved in 4 mL of 50 mM glycine-NaOH (pH 9.0) buffer and post-consumer PET fragment. The filled membrane bag was immersed into a container containing 196 mL of 50 mM glycine-NaOH (pH 9.0) buffer and incubated on a shaker for 84 h at 60 °C and 80 rpm. Samples were removed from the solution on the outside at different time points (0, 12, 24, 30, 36, 48, 60, 72, and 84 h) to quantify the total amount of released PET monomers. The samples were analyzed using HPLC. All experiments were performed in triplicate.

Samples were analyzed using a CBM-20A (Shimadzu, Kyoto, Japan) connected to a UV/Vis detector (SPD-20A) and a C18 column (Shimadzu Shim-pack GIS ODS-I C18, 5 µm, 4.6 × 150 mm). Mobile phases A (0.1% formic acid) and B (20% acetonitrile) were used at a flow rate of 1 mL min−1. The elution program was as follows: 0–2 min, 60%–80% buffer B, linear gradient; 2–8 min, 80%–100% buffer B, linear gradient; 8–16 min, 100% buffer B. The hydrolysis products (TPA, MHET, and BHET) were separated at 40 °C and detected at 260 nm. Amounts of TPA, MHET, and BHET were measured using standard curves prepared from commercial TPA (Sigma-Aldrich, MO, USA, Cat. No.: 185361), BHET (Sigma-Aldrich, MO, USA, Cat. No.: 465151), and MHET (Ambeed, IL, USA, Cat. No.: A875019).

The thermal crystallinity of PET samples was estimated using a differential scanning calorimetry (DSC) instrument (Q2000, TA instrument). For the DSC thermographs, 5 mg of samples were equilibrated at 30 °C and heated to 300 °C at a heating rate of 10 °C min−1. The sample was quenched to 30 °C at 10 °C min−1. Heat of fusion (ΔHm) and cold crystallization (ΔHc) were determined by integrating the areas (J g-1) under the peaks. The enthalpy of melting for a fully crystalline PET (∆Hf) was 140.1 J/g. The percent of crystallinity was calculated based on the following equation:

DSC curves for the crystallinity of the PET samples are presented in Supplementary Fig. 29.

The molecular weight of the PET samples was measured using gel-permeation chromatography (GPC) using an EcoSEC HLC-8320 GPC instrument (Tosoh Bioscience, Tokyo, Japan). The TSKgel SuperAWM-H column was used with a mobile phase of 0.01 N sodium trifluoroacetate (NaTFA) in hexafluoroisopropanol (HFIP) at a flow rate 0.3 mL min−1 at 40 °C, and 30 min run time per sample. The injected sample volume was 20 μL. The molecular weights (Mn, Mw, Mp) were determined relative to polymethylmethacrylate (PMMA) standards. The PET samples and PMMA were dissolved in 6% v/v hexafluoroisopropanol (HFIP) in chloroform to a concentration of 3 mg mL−1, and the solution was filtered through a polytetrafluoroethylene (PTFE) membrane with a pore size of 0.45 μm.

Further information on research design is available in the Nature Portfolio Reporting Summary linked to this article.

The refined models of CaPETaseWT, CaPETaseM9 generated in this study have been deposited in Protein Data Bank with PDB codes 7YM9 and 7YME. We have used the following published structure, PDB code 5LUI, as a search model for the initial molecular replacement. The raw data for figures generated in this study are provided in the Source Data file. The percent identity values for sequences used in phylogenetic analysis were provided in the Supplementary Fig. 31. The data that support the findings of this study are available from the corresponding author upon request. Source data are provided with this paper.

Geyer, R., Jambeck, J. R. & Law, K. L. Production, use, and fate of all plastics ever made. Sci. Adv. 3, e1700782 (2017).

Article ADS PubMed PubMed Central Google Scholar

Gibb, B. C. Plastics are forever, Nature Publishing Group, (2019).

Hopewell, J., Dvorak, R. & Kosior, E. Plastics recycling: challenges and opportunities. Philos. Trans. R. Soc. B: Biol. Sci. 364, 2115–2126 (2009).

Article CAS Google Scholar

Jambeck, J. R. et al. Marine pollution. Plastic waste inputs from land into the ocean. Science 347, 768–771 (2015).

Article ADS CAS PubMed Google Scholar

Eriksen, M. et al. Plastic pollution in the world’s oceans: more than 5 trillion plastic pieces weighing over 250,000 tons afloat at sea. PloS one 9, e111913 (2014).

Article ADS PubMed PubMed Central Google Scholar

Qi, Y. et al. Macro-and micro-plastics in soil-plant system: effects of plastic mulch film residues on wheat (Triticum aestivum) growth. Sci. Total Environ. 645, 1048–1056 (2018).

Article ADS CAS PubMed Google Scholar

Borrelle, S. B. et al. Predicted growth in plastic waste exceeds efforts to mitigate plastic pollution. Science 369, 1515–1518 (2020).

Article ADS CAS PubMed Google Scholar

Ellis, L. D. et al. Chemical and biological catalysis for plastics recycling and upcycling. Nat. Catal. 4, 539–556 (2021).

Article CAS Google Scholar

Kakadellis, S. & Rosetto, G. Achieving a circular bioeconomy for plastics. Science 373, 49–50 (2021).

Article ADS CAS PubMed Google Scholar

Coates, G. W. & Getzler, Y. D. Chemical recycling to monomer for an ideal, circular polymer economy. Nat. Rev. Mater. 5, 501–516 (2020).

Article ADS CAS Google Scholar

Webb, H. K., Arnott, J., Crawford, R. J. & Ivanova, E. P. Plastic degradation and its environmental implications with special reference to poly (ethylene terephthalate). Polymers 5, 1–18 (2012).

Article Google Scholar

Kuczenski, B. & Geyer, R. Material flow analysis of polyethylene terephthalate in the US, 1996–2007. Resour., Conserv. Recycling 54, 1161–1169 (2010).

Article Google Scholar

Soong, Y.-H. V., Sobkowicz, M. J. & Xie, D. Recent advances in biological recycling of polyethylene terephthalate (PET) plastic wastes. Bioengineering 9, 98 (2022).

Article CAS PubMed PubMed Central Google Scholar

Zimmermann, W. Biocatalytic recycling of polyethylene terephthalate plastic. Philos. Trans. R. Soc. A 378, 20190273 (2020).

Article ADS CAS Google Scholar

Tiso, T. et al. Towards bio-upcycling of polyethylene terephthalate. Metab. Eng. 66, 167–178 (2021).

Article CAS PubMed Google Scholar

Kawai, F., Kawabata, T. & Oda, M. Current knowledge on enzymatic PET degradation and its possible application to waste stream management and other fields. Appl. Microbiol. Biotechnol. 103, 4253–4268 (2019).

Article CAS PubMed PubMed Central Google Scholar

Yoshida, S. et al. A bacterium that degrades and assimilates poly (ethylene terephthalate). Science 351, 1196–1199 (2016).

Article ADS CAS PubMed Google Scholar

Ma, Y. et al. Enhanced poly (ethylene terephthalate) hydrolase activity by protein engineering. Engineering 4, 888–893 (2018).

Article CAS Google Scholar

Liu, B. et al. Protein crystallography and site‐direct mutagenesis analysis of the poly (ethylene terephthalate) hydrolase PETase from Ideonella sakaiensis. ChemBioChem 19, 1471–1475 (2018).

Article CAS PubMed Google Scholar

Austin, H. P. et al. Characterization and engineering of a plastic-degrading aromatic polyesterase. Proc. Natl. Acad. Sci. 115, E4350–E4357 (2018).

Article CAS PubMed PubMed Central Google Scholar

Han, X. et al. Structural insight into catalytic mechanism of PET hydrolase. Nat. Commun. 8, 1–6 (2017).

Article ADS MathSciNet Google Scholar

Joo, S. et al. Structural insight into molecular mechanism of poly (ethylene terephthalate) degradation. Nat. Commun. 9, 1–12 (2018).

Article ADS Google Scholar

Cui, Y. et al. Computational redesign of a PETase for plastic biodegradation under ambient condition by the GRAPE strategy. ACS Catal. 11, 1340–1350 (2021).

Article CAS Google Scholar

Lu, H. et al. Machine learning-aided engineering of hydrolases for PET depolymerization. Nature 604, 662–667 (2022).

Article ADS CAS PubMed Google Scholar

Sulaiman, S. 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 (2012).

Article ADS CAS PubMed PubMed Central Google Scholar

Tournier, V. et al. An engineered PET depolymerase to break down and recycle plastic bottles. Nature 580, 216–219 (2020).

Article ADS CAS PubMed Google Scholar

Then, J. et al. Ca2+ and Mg2+ binding site engineering increases the degradation of polyethylene terephthalate films by polyester hydrolases from Thermobifida fusca. Biotechnol. J. 10, 592–598 (2015).

Article CAS PubMed Google Scholar

Wei, R. et al. Functional characterization and structural modeling of synthetic polyester-degrading hydrolases from Thermomonospora curvata. AMB express 4, 1–10 (2014).

Article MathSciNet CAS Google Scholar

Xi, X. et al. Secretory expression in bacillus subtilis and biochemical characterization of a highly thermostable polyethylene terephthalate hydrolase from bacterium HR29. Enzym. Microb. Technol. 143, 109715 (2021).

Article CAS Google Scholar

Heredia, A. Biophysical and biochemical characteristics of cutin, a plant barrier biopolymer. Biochimica et. Biophys. Acta (BBA)-Gen. Subj. 1620, 1–7 (2003).

Article CAS Google Scholar

Ronkvist, Å. M., Xie, W., Lu, W. & Gross, R. A. Cutinase-catalyzed hydrolysis of poly(ethylene terephthalate). Macromolecules 42, 5128 (2009).

Article ADS CAS Google Scholar

Ribitsch, D. et al. Characterization of a new cutinase from Thermobifida alba for PET-surface hydrolysis. Biocatal. Biotransform. 30, 2 (2012).

Article CAS Google Scholar

Danso, D. 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–02717 (2018).

Article CAS PubMed PubMed Central Google Scholar

Sonnendecker, C. et al. Low carbon footprint recycling of post‐consumer PET plastic with a metagenomic polyester hydrolase. ChemSusChem 15, e202101062 (2022).

CAS PubMed PubMed Central Google Scholar

Erickson, E. et al. Sourcing thermotolerant poly (ethylene terephthalate) hydrolase scaffolds from natural diversity. Nat. Commun. 13, 7850 (2022).

Article ADS CAS PubMed PubMed Central Google Scholar

Thomsen, T. B., Hunt, C. J. & Meyer, A. S. Influence of substrate crystallinity and glass transition temperature on enzymatic degradation of polyethylene terephthalate (PET). N. Biotechnol. 69, 28–35 (2022).

Article CAS PubMed Google Scholar

Lee, S. H. et al. Balance-directed protein engineering of IsPETase enhances both PET hydrolysis activity and thermostability. Preprint at https://doi.org/10.1101/2022.11.14.516528 (2022).

Kufareva, I. & Abagyan, R. Methods of protein structure comparison. Methods Mol. Biol. 857, 231–257 (2012).

Charupanit, K., Tipmanee, V., Sutthibutpong, T. & Limsakul, P. In silico identification of potential sites for a plastic-degrading enzyme by a reverse screening through the protein sequence space and molecular dynamics simulations. Molecules 27, 3353 (2022).

Article CAS PubMed PubMed Central Google Scholar

Fecker, T. et al. Active site flexibility as a hallmark for efficient PET degradation by I. Sakaiensis Petase. Biophys. J. 114, 1302 (2018).

Article ADS CAS PubMed PubMed Central Google Scholar

Sagong, H. Y. et al. Implications for the PET decomposition mechanism through similarity and dissimilarity between PETases from Rhizobacter gummiphilus and Ideonella sakaiensis. J. Hazard Mater. 416, https://doi.org/10.1016/j.jhazmat.2021.126075 (2021).

Son, H. F. et al. Rational protein engineering of thermo-stable PETase from Ideonella sakaiensis for highly efficient PET degradation. ACS Catal. 9, 3519–3526 (2019).

Article CAS Google Scholar

Oda, M. et al. Enzymatic hydrolysis of PET: functional roles of three Ca2+ ions bound to a cutinase-like enzyme, Cut190*, and its engineering for improved activity. Appl. Microbiol. Biotechnol. 102, 10067–10077 (2018).

Article CAS PubMed Google Scholar

Nakamura, A., Kobayashi, N., Koga, N. & Iino, R. Positive Charge Introduction on the Surface of Thermostabilized PET Hydrolase Facilitates PET Binding and Degradation. ACS Catal. 11, 8550–8564 (2021).

Article CAS Google Scholar

Crippa, M. et al. A circular economy for plastics: Insights from research and innovation to inform policy and funding decisions. (European Comission Publications, 2019).

Carr, C. M., Clarke, D. J. & Dobson, A. D. Microbial polyethylene terephthalate hydrolases: current and future perspectives. Front. Microbiol. 11, 571265 (2020).

Article PubMed PubMed Central Google Scholar

Magalhães, R. P., Cunha, J. M. & Sousa, S. F. Perspectives on the Role of Enzymatic Biocatalysis for the Degradation of Plastic PET. Int. J. Mol. Sci. 22, 11257 (2021).

Article PubMed PubMed Central Google Scholar

Chen, C. C. et al. General Features to Enhance Enzymatic Activity of Poly(Ethylene Terephthalate) Hydrolysis. Nat. Catal. 4, 425 (2021).

Article CAS Google Scholar

Sievers, F. et al. Fast, scalable generation of high‐quality protein multiple sequence alignments using Clustal Omega. Mol. Syst. Biol. 7, 539 (2011).

Article PubMed PubMed Central Google Scholar

Schwartz, R. Matrices for detecting distant relationships. Atlas of protein sequence and structure 5, 353–359 (1978).

Tamura, K., Stecher, G. & Kumar, S. MEGA11: molecular evolutionary genetics analysis version 11. Mol. Biol. evolution 38, 3022–3027 (2021).

Article CAS Google Scholar

Felsenstein, J. Confidence limits on phylogenies: an approach using the bootstrap. evolution 39, 783–791 (1985).

Article PubMed Google Scholar

Robert, X. & Gouet, P. Deciphering key features in protein structures with the new ENDscript server. Nucleic Acids Res. 42, W320–W324 (2014).

Article CAS PubMed PubMed Central Google Scholar

Park, S.-Y., Ha, S.-C. & Kim, Y.-G. The protein crystallography beamlines at the pohang light source II. Biodesign 5, 30–34 (2017).

Google Scholar

Otwinowski, Z. & Minor, W. [20] Processing of X-ray diffraction data collected in oscillation mode. Methods Enzymol. 276, 307–326 (1997).

Article CAS PubMed Google Scholar

Matthews, B. W. Solvent content of protein crystals. J. Mol. Biol. 33, 491–497 (1968).

Article CAS PubMed Google Scholar

Vagin, A. & Teplyakov, A. Molecular replacement with MOLREP. Acta Crystallogr D. Biol. Crystallogr 66, 22–25 (2010).

Article CAS PubMed Google Scholar

Ribitsch, D. et al. Small cause, large effect: structural characterization of cutinases from Thermobifida cellulosilytica. Biotechnol. Bioeng. 114, 2481–2488 (2017).

Article CAS PubMed Google Scholar

Emsley, P. & Cowtan, K. Coot: model-building tools for molecular graphics. Acta Crystallogr D. Biol. Crystallogr 60, 2126–2132 (2004).

Article PubMed Google Scholar

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).

Article CAS PubMed Google Scholar

Olsson, M. H., Søndergaard, C. R., Rostkowski, M. & Jensen, J. H. PROPKA3: consistent treatment of internal and surface residues in empirical p K a predictions. J. Chem. theory Comput. 7, 525–537 (2011).

Article CAS PubMed Google Scholar

Páll, S. et al. Heterogeneous parallelization and acceleration of molecular dynamics simulations in GROMACS. J. Chem. Phys. 153, 134110 (2020).

Article ADS PubMed Google Scholar

Download references

This research was supported by the Bio & Medical Technology Development Program of the National Research Foundation (NRF) funded by the Ministry of Science & ICT (NRF-2020M3A9I5037635) and the Cooperative Research Program for Agricultural Science & Technology Development (Project No. PJ01492602), Rural Development Administration, Republic of Korea. This research was also supported by the Technology Innovation Program (20018132, Development of the biodegradable polybutylene plastics, PBAT and PBS, from biomass) funded By the Ministry of Trade, Industry & Energy (MOTIE, Korea).

These authors contributed equally: Hwaseok Hong, Dongwoo Ki.

School of Life Sciences, BK21 FOUR KNU Creative BioResearch Group, KNU Institute for Microorganisms, Kyungpook National University, Daegu, 41566, Republic of Korea

Hwaseok Hong, Dongwoo Ki, Hogyun Seo, Jiyoung Park & Kyung-Jin Kim

Institute of Biotechnology, CJ CheilJedang Co., Suwon-si, Gyeonggi-do, 16495, Republic of Korea

Jaewon Jang

Zyen Co, Daegu, 41566, Republic of Korea

Kyung-Jin Kim

You can also search for this author in PubMed Google Scholar

You can also search for this author in PubMed Google Scholar

You can also search for this author in PubMed Google Scholar

You can also search for this author in PubMed Google Scholar

You can also search for this author in PubMed Google Scholar

You can also search for this author in PubMed Google Scholar

K.-J.K. conceived the project and supervised the research. H.H. and D.K. performed the crystallographic experiment and S.H. contributed to the computational analysis. H.H. and D.K. designed the enzyme variants and performed the PET degradation experiments. H.H., D.K., and J.P. analyzed the data. H.H. and K.-J.K. participated in the manuscript writing. All author discussed results and S.H. and J.J. provided guidance of this manuscript.

Correspondence to Kyung-Jin Kim.

The authors declare no competing interest.

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 Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article’s Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article’s Creative Commons licence and your intended use is not permitted 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 licence, visit http://creativecommons.org/licenses/by/4.0/.

Reprints and Permissions

Hong, H., Ki, D., Seo, H. et al. Discovery and rational engineering of PET hydrolase with both mesophilic and thermophilic PET hydrolase properties. Nat Commun 14, 4556 (2023). https://doi.org/10.1038/s41467-023-40233-w

Download citation

Received: 04 January 2023

Accepted: 13 July 2023

Published: 28 July 2023

DOI: https://doi.org/10.1038/s41467-023-40233-w

Anyone you share the following link with will be able to read this content:

Sorry, a shareable link is not currently available for this article.

Provided by the Springer Nature SharedIt content-sharing initiative

By submitting a comment you agree to abide by our Terms and Community Guidelines. If you find something abusive or that does not comply with our terms or guidelines please flag it as inappropriate.