Computational biotransformation of polyethylene terephthalate by depolymerase: A QM/MM approach

Despite increasing environmental concerns on ever-lasting Polyethylene Terephthalate (PET), its global production is continuously growing. Effective strategies that can completely remove PET from environment are urgently desired. Here biotransformation processes of PET by one of the most effective enzymes, leaf-branch compost cutinase (LCC), were systematically explored with Molecular Dynamics and Quantum Mechanics/Molecular Mechanics approaches. We found that four concerted steps are required to complete the whole catalytic cycle. The last concerted step, deacylation, was determined as the rate-determining step with Boltzmann-weighted average barrier of 13.6 kcal/mol and arithmetic average of 16.1 ± 2.9 kcal/mol. Interestingly, unprecedented fluctuations of hydrogen bond length during LCC catalyzed transformation process toward PET were found. This fluctuation was also observed in enzyme IsPETase, indicating that it may widely exist in other catalytic triad (Ser-His-Asp) containing enzymes as well. In addition, possible features (bond, angle, dihedral angle and charge) that influence the catalytic reaction were identified and correlations between activation energies and key features were established. Our results present new insights into catalytic mechanism of hydrolases and shed light on the efficient recycling of the ever-lasting PET.