Acylation and deacylation mechanism and kinetics of penicillin G reaction with StreptomycesR61 DD‐peptidase

Abstract

Two quantum mechanical (QM)-cluster models are built for studying the acylation and deacylation mechanism and kinetics of Streptomyces R61 DD-peptidase with the penicillin G at atomic level detail. DD-peptidases are bacterial enzymes involved in the cross-linking of peptidoglycan to form the cell wall, necessary for bacterial survival. The cross-linking can be inhibited by antibiotic beta-lactam derivatives through acylation, preventing the acyl-enzyme complex from undergoing further deacylation.The deacylation step was predicted to be rate-limiting. Transition state and intermediate structures are found using density functional theory in this study, and thermodynamic and kinetic properties of the proposed mechanism are evaluated. The acylenzyme complex is found lying in a deep thermodynamic sink, and deacylation is indeed the severely rate-limiting step, leading to suicide inhibition of the peptidoglycan cross-linking. The usage of QM-cluster models is a promising technique to understand, improve, and design antibiotics to disrupt function of the Streptomyces R61

QM-Cluster Model Study of the Guaiacol Hydrogen Atom Transfer and Oxygen Rebound with Cytochrome P450 Enzyme GcoA

“…Our group has employed Residue Interaction Networks (RINs) , to algorithmically generate QM-cluster models. RINs have been used in the context of QM-cluster modeling to study various enzymatic reactions. − The in-house software toolkit Residue Interaction Network-based ResidUe Selector (RINRUS) is under development to simplify and automate computational enzymology workflows. The size and shape of the QM regions are very important to consider when constructing models in computational enzymology , but general strategies have not been well characterized.…”

The Glycine N-Methyltransferase Case Study: Another Challenge for QM-Cluster Models?

“…This dielectric constant value has been previously determined as appropriate for simulating the less-polarized environment within an enzyme active site. , Implicit solvation was incorporated into all geometry optimizations and harmonic frequency calculations. The computational methodology has been successfully employed in several enzyme studies which have reliable comparisons to experimental results. ,− Unscaled harmonic vibrational frequency calculations were used to identify all stationary points as either minima (no imaginary frequency) or transition states (only one imaginary frequency). Reactants and products corresponding to the methyl transfer TS were located by following the intrinsic reaction coordinate (IRC). , The same “freeze code” scheme of Gaussian 16 used in other enzyme studies by our group is used in this study.…”

The influence of model building schemes and molecular dynamics sampling on QM-cluster models: the chorismate mutase case study

“…We will use the F-SAPT interaction energies between chorismate and surrounding residue fragments to rank incremental QM-cluster model building. This work uses the zeroth-order formulation of F-SAPT, F-SAPT0, described by the equation: (1) elec + E (1) exch + [E (2) ind + E (2) exch-ind + dE (2) HF ] ind…”

“…Our group has developed the Residue Interaction Network ResidUe Selector (RINRUS) software toolkit to facilitate studying the reaction mechanisms of enzymes with quantum chemistry. [2][3][4] Instead of relying on chemical intuition or distance-based criteria to prioritize the critical fragments within the enzyme active site, RINRUS algorithmically constructs enzyme models based on several possible qualitative and quantitative criteria. RINRUS infrastructure was first developed to build QMcluster models of enzymes, but adapting the code to also build QM/MM enzyme models is in progress.…”