RESUMO
Docking studies play a critical role in the current workflow of drug discovery. However, limitations may often arise through factors including inadequate ligand sampling, a lack of protein flexibility, scoring function inadequacies (e.g., due to metals, co-factors, etc.), and difficulty in retaining explicit water molecules. Herein, we present a novel CHARMM-based induced fit docking (CIFDock) workflow that can circumvent these limitations by employing all-atom force fields coupled to enhanced sampling molecular dynamics procedures. Self-guided Langevin dynamics simulations are used to effectively sample relevant ligand conformations, side chain orientations, crystal water positions, and active site residue motion. Protein flexibility is further enhanced by dynamic sampling of side chain orientations using an expandable rotamer library. Steps in the procedure consisting of fixing individual components (e.g., the ligand) while sampling the other components (e.g., the residues in the active site of the protein) allow for the complex to adapt to conformational changes. Ultimately, all components of the complex-the protein, ligand, and waters-are sampled simultaneously and unrestrained with SGLD to capture any induced fit effects. This modular flexible docking procedure is automated using CHARMM scripting, interfaced with SLURM array processing, and parallelized to use the desired number of processors. We validated the CIFDock procedure by performing cross-docking studies using a data set comprised of 21 pharmaceutically relevant proteins. Five variants of the CHARMM-based SWISSDOCK scoring functions were created to quantify the results of the final generated poses. Results obtained were comparable to, or in some cases improved upon, commercial docking program data.
Assuntos
Simulação de Acoplamento Molecular , Proteínas/química , Ligantes , Termodinâmica , Água/químicaRESUMO
Accurately predicting free energy differences is essential in realizing the full potential of rational drug design. Unfortunately, high levels of accuracy often require computationally expensive QM/MM Hamiltonians. Fortuitously, the cost of employing QM/MM approaches in rigorous free energy simulation can be reduced through the use of the so-called "indirect" approach to QM/MM free energies, in which the need for QM/MM simulations is avoided via a QM/MM "correction" at the classical endpoints of interest. Herein, we focus on the computation of QM/MM binding free energies in the context of the SAMPL8 Drugs of Abuse host-guest challenge. Of the 5 QM/MM correction coupled with force-matching submissions, PM6-D3H4/MM ranked submission proved the best overall QM/MM entry, with an RMSE from experimental results of 2.43 kcal/mol (best in ranked submissions), a Pearson's correlation of 0.78 (second-best in ranked submissions), and a Kendall [Formula: see text] correlation of 0.52 (best in ranked submissions).
Assuntos
Simulação de Dinâmica Molecular , Proteínas , Ligantes , Ligação Proteica , Teoria Quântica , TermodinâmicaRESUMO
Poly(ethylene terephthalate) (PET) is one of the most abundantly produced synthetic polymers and is accumulating in the environment at a staggering rate as discarded packaging and textiles. The properties that make PET so useful also endow it with an alarming resistance to biodegradation, likely lasting centuries in the environment. Our collective reliance on PET and other plastics means that this buildup will continue unless solutions are found. Recently, a newly discovered bacterium, Ideonella sakaiensis 201-F6, was shown to exhibit the rare ability to grow on PET as a major carbon and energy source. Central to its PET biodegradation capability is a secreted PETase (PET-digesting enzyme). Here, we present a 0.92 Å resolution X-ray crystal structure of PETase, which reveals features common to both cutinases and lipases. PETase retains the ancestral α/ß-hydrolase fold but exhibits a more open active-site cleft than homologous cutinases. By narrowing the binding cleft via mutation of two active-site residues to conserved amino acids in cutinases, we surprisingly observe improved PET degradation, suggesting that PETase is not fully optimized for crystalline PET degradation, despite presumably evolving in a PET-rich environment. Additionally, we show that PETase degrades another semiaromatic polyester, polyethylene-2,5-furandicarboxylate (PEF), which is an emerging, bioderived PET replacement with improved barrier properties. In contrast, PETase does not degrade aliphatic polyesters, suggesting that it is generally an aromatic polyesterase. These findings suggest that additional protein engineering to increase PETase performance is realistic and highlight the need for further developments of structure/activity relationships for biodegradation of synthetic polyesters.
Assuntos
Proteínas de Bactérias/química , Burkholderiales/enzimologia , Esterases/química , Polietilenotereftalatos/química , Proteínas de Bactérias/genética , Burkholderiales/genética , Cristalografia por Raios X , Esterases/genética , Engenharia de Proteínas , Especificidade por SubstratoRESUMO
Polyethylene terephthalate (PET), the most abundantly produced polyester plastic, can be depolymerized by the Ideonella sakaiensis PETase enzyme. Based on multiple PETase crystal structures, the reaction has been proposed to proceed via a two-step serine hydrolase mechanism mediated by a serine-histidine-aspartate catalytic triad. To elucidate the multi-step PETase catalytic mechanism, we use transition path sampling and likelihood maximization to identify optimal reaction coordinates for the PETase enzyme. We predict that deacylation is likely rate-limiting, and the reaction coordinates for both steps include elements describing nucleophilic attack, ester bond cleavage, and the "moving-histidine" mechanism. We find that the flexibility of Trp185 promotes the reaction, providing an explanation for decreased activity observed in mutations that restrict Trp185 motion. Overall, this study uses unbiased computational approaches to reveal the detailed reaction mechanism necessary for further engineering of an important class of enzymes for plastics bioconversion.