BhrPETase catalyzed polyethylene terephthalate depolymerization: A quantum mechanics/molecular mechanics approach

Polyethylene terephthalate (PET) is a widely used material in our daily life, particularly in areas such as packaging, fibers, and engineering plastics. However, PET waste can accumulate in the environment and pose a great threat to our ecosystem. Recently enzymatic conversion has emerged as an efficient and green strategy to address the PET crisis. Here, using a theoretical approach combining molecular dynamics simulation and quantum mechanics/molecular mechanics calculations, the depolymerization mechanism of the thermophilic cutinase BhrPETase was fully deciphered. Surprisingly, unlike the previously studied cutinase LCC^ICCG, our results indicate that the first step, catalytic triad assisted nucleophilic attack, is the rate-determining step. The corresponding Boltzmann weighted average energy barrier is 18.2 kcal/mol. Through extensive comparison between BhrPETase and LCC^ICCG, we evidence that key features like charge C_His@N1 and angle A_{PET@C1-Ser@O1-His@H1} significantly impact the depolymerization efficiency of BhrPETase. Non-covalent bond interaction and distortion/interaction analysis inform new insights on enzyme engineer and may aid the recycling of enzymatic PET waste. This study will aid the advancement of the plastic bio-recycling economy and promote resource conservation and reuse.