Influence of Wobbling Tryptophan and Mutations on PET Degradation Explored by QM/MM Free Energy Calculations.

Plastic-degrading enzymes, particularly poly(ethylene terephthalate) (PET) hydrolases, have garnered significant attention in recent years as potential eco-friendly solutions for recycling plastic waste. However, understanding of their PET-degrading activity and influencing factors remains incomplete, impeding the development of uniform approaches for enhancing PET hydrolases for industrial applications. A key aspect of PET hydrolase engineering is optimizing the PET-hydrolysis reaction by lowering the associated free energy barrier. However, inconsistent findings have complicated these efforts. Therefore, our goal is to elucidate various aspects of enzymatic PET degradation by means of quantum mechanics/molecular mechanics (QM/MM) reaction simulations and analysis, focusing on the initial reaction step, acylation, in two thermophilic PET hydrolases, LCC and PES-H1, along with their highly active variants, LCCIG and PES-H1FY. Our findings highlight the impact of semiempirical QM methods on proton transfer energies, affecting the distinction between a two-step reaction involving a metastable tetrahedral intermediate and a one-step reaction. Moreover, we uncovered a concerted conformational change involving the orientation of the PET benzene ring, altering its interaction with the side-chain of the "wobbling" tryptophan from T-stacking to parallel π-π interactions, a phenomenon overlooked in prior research. Our study thus enhances the understanding of the acylation mechanism of PET hydrolases, in particular by characterizing it for the first time for the promising PES-H1FY using QM/MM simulations. It also provides insights into selecting a suitable QM method and a reaction coordinate, valuable for future studies on PET degradation processes.

Multiple Substrate Binding Mode-Guided Engineering of a Thermophilic PET Hydrolase.

Thermophilic polyester hydrolases (PES-H) have recently enabled biocatalytic recycling of the mass-produced synthetic polyester polyethylene terephthalate (PET), which has found widespread use in the packaging and textile industries. The growing demand for efficient PET hydrolases prompted us to solve high-resolution crystal structures of two metagenome-derived enzymes (PES-H1 and PES-H2) and notably also in complex with various PET substrate analogues. Structural analyses and computational modeling using molecular dynamics simulations provided an understanding of how product inhibition and multiple substrate binding modes influence key mechanistic steps of enzymatic PET hydrolysis. Key residues involved in substrate-binding and those identified previously as mutational hotspots in homologous enzymes were subjected to mutagenesis. At 72 °C, the L92F/Q94Y variant of PES-H1 exhibited 2.3-fold and 3.4-fold improved hydrolytic activity against amorphous PET films and pretreated real-world PET waste, respectively. The R204C/S250C variant of PES-H1 had a 6.4 °C higher melting temperature than the wild-type enzyme but retained similar hydrolytic activity. Under optimal reaction conditions, the L92F/Q94Y variant of PES-H1 hydrolyzed low-crystallinity PET materials 2.2-fold more efficiently than LCC ICCG, which was previously the most active PET hydrolase reported in the literature. This property makes the L92F/Q94Y variant of PES-H1 a good candidate for future applications in industrial plastic recycling processes.

Dataset of AMBER force field parameters of drugs, natural products and steroids for simulations using GROMACS.

We provide general AMBER force field (GAFF) parameters for 160 organic molecules including drugs, natural products, and steroids, which can be employed without further processing in molecular dynamics (MD) simulations using GROMACS. We determined these parameters based on quantum mechanical (QM) calculations involving geometry optimization at the HF6-31G* level of theory. For each molecule we provide a coordinate file of the three-dimensional molecular structure, the topology and the parameter file. The applicability of these parameters was demonstrated by MD simulations of these molecules bound to the active site of the main protease of the coronavirus SARS-CoV-2, 3CLpro, which is a main player during viral replication causing COVID-19.

Novel inhibitors of the main protease enzyme of SARS-CoV-2 identified via molecular dynamics simulation-guided in vitro assay.

For the COVID-19 pandemic caused by SARS-CoV-2, there are currently no effective drugs or vaccines to treat this coronavirus infection. In this study, we focus on the main protease enzyme of SARS-CoV-2, 3CLpro, which is critical for viral replication. We employ explicit solvent molecular dynamics simulations of about 150 compounds docked into 3CLpro's binding site and that had emerged as good main protease ligands from our previous in silico screening of over 1.2 million compounds. By incoporating protein dynamics and applying a range of structural descriptors, such as the ability to form specific contacts with the catalytic dyad residues of 3CLpro and the structural fluctuations of the ligands in the binding site, we are able to further refine our compound selection. Fourteen compounds including estradiol shown to be the most promising based on our calculations were procured and screened against recombinant 3CLpro in a fluorescence assay. Eight of these compounds have significant activity in inhibiting the SARS-CoV-2 main protease. Among these are corilagin, a gallotannin, and lurasidone, an antipsychotic drug, which emerged as the most promising natural product and drug, respectively, and might thus be candidates for drug repurposing for the treatment of COVID-19. In addition, we also tested the inhibitory activity of testosterone, and our results reveal testosterone as possessing moderate inhibitory potency against the 3CLpro enzyme, which may thus provide an explanation why older men are more severely affected by COVID-19.