How To Produce Methane Precursor in the Upper Ocean by An Untypical Non‐Heme Fe‐Dependent Methylphosphonate Synthase?

Tutorial Review on the Set-Up and Running of Quantum Mechanical Cluster Models for Enzymatic Reaction Mechanisms

Enzymes turnover substrates into products with amazing efficiency and selectivity and as such have great potential for use in biotechnology and pharmaceutical applications. However, details of their catalytic cycles and the origins surrounding the regio- and chemoselectivity of enzymatic reaction processes remain unknown, which makes the engineering of enzymes and their use in biotechnology challenging. Computational modelling can assist experimental work in the field and establish the factors that influence the reaction rates and the product distributions. A popular approach in modelling is the use of quantum mechanical cluster models of enzymes that take the first- and second coordination sphere of the enzyme active site into consideration. These QM cluster models are widely applied but often the results obtained are dependent on model choice and model selection. Herein, we show that QM cluster models can give highly accurate results that reproduce experimental product distributions and free energies of activation within several kcal mol-1, regarded that large cluster models with >300 atoms are used that include key hydrogen bonding interactions and charged residues. In this tutorial review, we give general guidelines on the set-up and applications of the QM cluster method and discuss its accuracy and reproducibility. Finally, several representative QM cluster model examples on metal-containing enzymes are presented, which highlight the strength of the approach.

Mechanistic Insights into the N-Hydroxylations Catalyzed by the Binuclear Iron Domain of SznF Enzyme: Key Piece in the Synthesis of Streptozotocin

SznF, a member of the emerging family of heme-oxygenase-like (HO-like) di-iron oxidases and oxygenases, employs two distinct domains to catalyze the conversion of Nω-methyl-L-arginine (L-NMA) into N-nitroso-containing product, which can subsequently be transformed into streptozotocin. Using unrestricted density functional theory (UDFT) with the hybrid functional B3LYP, we have mechanistically investigated the two sequential hydroxylations of L-NMA catalyzed by SznF's binuclear iron central domain. Mechanism B primarily involves the O-O bond dissociation, forming Fe(IV)=O, induced by the H+/e- introduction to the FeA side of μ-1,2-peroxo-Fe2(III/III), the substrate hydrogen abstraction by Fe(IV)=O, and the hydroxyl rebound to the substrate N radical. The stochastic addition of H+/e- to the FeB side (mechanism C) can transition to mechanism B, thereby preventing enzyme deactivation. Two other competing mechanisms, involving the direct O-O bond dissociation (mechanism A) and the addition of H2O as a co-substrate (mechanism D), have been ruled out.

Reaction mechanism of the PuDddK dimethylsulfoniopropionate lyase and cofactor effects of various transition metal ions

The microbial cleavage of dimethylsulfoniopropionate (DMSP) produces volatile dimethyl sulfide (DMS) via the lyase pathway, playing a crucial role in the global sulfur cycle. Herein, the DMSP decomposition catalyzed by PuDddK (a DMSP lyase) devised with various transition metal ion cofactors are investigated using density functional calculations. The PuDddK reaction has been demonstrated to employ a concerted β-elimination mechanism, where the substrate α-proton abstraction by the deprotonated Tyr64 occurs simultaneously with the C_β-S bond cleavage and C_α = C_β double bond formation. The PuDddK enzymes with diverse divalent metal ions (Ni²⁺, Mn²⁺, Fe²⁺, Co²⁺, Zn²⁺, and Cu²⁺) incorporated prefer DMSP as a monodentate ligand. The cases of Ni²⁺, Mn²⁺, Fe²⁺, Co²⁺, and Zn²⁺ with the same 3His-1Glu ligands have close reaction energy barriers, indicating that the lyase activity may be hardly affected by the divalent transition metal type with the same ligand type and number. The coordination loss of one histidine in Cu²⁺, forming a 2His-1Glu architecture, leads to a lower activity, revealing that the 3His-1Glu ligand set used by DddK appears to be a scaffold capable of more efficiently catalyzing the DMSP decomposition. Further analysis reveals that the inactivation of Fe³⁺-dependent PuDddK is derived from an electron transfer from the Tyr64 phenolate to Fe³⁺, with the implication that the PuDddK activity may be primarily affected by the redox effects induced by a strongly oxidizing transition metal ion (like Fe³⁺).

Mechanism for the Halogenation and Azidation of Lysine Catalyzed by Non-heme Iron BesD Enzyme

Selective halogenation is important in synthetic chemistry. BesD, a new member of the non-heme Fe(II)/α-ketoglutarate (αKG)-dependent halogenase family, can activate the sp3 C-H bond and halogenate lysine, in particular without a carrier protein. Using the density functional calculations, a chlorination mechanism in BesD has been proposed, mainly including the formation of Cl-Fe(IV)=O through the αKG decarboxylation, the isomerization of Cl-Fe(IV)=O, the substrate hydrogen abstraction by Fe(IV)=O, and the rebound of chloro to the substrate carbon radical. The hydrogen abstraction is rate-limiting. The isomerization of Cl-Fe(IV)=O is essential for the hydrogen abstraction and the chiral selectivity. The BesD-catalyzed bromination and azidation of lysine adopt the same mechanism as the chlorination. The hardly-changed overall barriers indicate that the introduced ligands (X) do not affect the reaction rate significantly, implying that the X-introduced reactions catalyzed by BesD may be extended to other X anions.

The Se-S Bond Formation in the Covalent Inhibition Mechanism of SARS-CoV-2 Main Protease by Ebselen-like Inhibitors: A Computational Study

The inhibition mechanism of the main protease (Mpro) of SARS-CoV-2 by ebselen (EBS) and its analog with a hydroxyl group at position 2 of the benzisoselenazol-3(2H)-one ring (EBS-OH) was studied by using a density functional level of theory. Preliminary molecular dynamics simulations on the apo form of Mpro were performed taking into account both the hydrogen donor and acceptor natures of the Nδ and Nε of His41, a member of the catalytic dyad. The potential energy surfaces for the formation of the Se-S covalent bond mediated by EBS and EBS-OH on Mpro are discussed in detail. The EBS-OH shows a distinctive behavior with respect to EBS in the formation of the noncovalent complex. Due to the presence of canonical H-bonds and noncanonical ones involving less electronegative atoms, such as sulfur and selenium, the influence on the energy barriers and reaction energy of the Minnesota hybrid meta-GGA functionals M06, M06-2X and M08HX, and the more recent range-separated hybrid functional wB97X were also considered. The knowledge of the inhibition mechanism of Mpro by the small protease inhibitors EBS or EBS-OH can enlarge the possibilities for designing more potent and selective inhibitor-based drugs to be used in combination with other antiviral therapies.

The platination mechanism of RNase A by arsenoplatin: insight from the theoretical study

A detailed metalation process of the bovine pancreatic ribonuclease (RNase A) by a novel multitarget anti-cancer agent arsenoplatin-1, ([Pt(μ-NHC(CH₃)O)₂ClAs(OH)₂]), performed at DFT level and using different models size is provided.

Mechanistic Insights into the Oxidative Rearrangement Catalyzed by the Unprecedented Dioxygenase ChaP Involved in Chartreusin Biosynthesis

ChaP is a non-heme iron-dependent dioxygenase belonging to the vicinal oxygen chelate (VOC) enzyme superfamily that catalyzes the final α-pyrone ring formation in the biosynthesis of chartreusin. In contrast to other common dioxygenases, for example, 2,3-catechol dioxygenase which uses the dioxygen molecule as the oxidant, ChaP requires the flavin-activated oxygen (O₂^2-) as the equivalent. Previous experiments showed that the ChaP-catalyzed ring rearrangement contains two successive C-C bond cleavages and one lactonization; however, the detailed reaction mechanism is unknown. In this work, on the basis of the recently obtained crystal structure of ChaP, the computational model was constructed and the catalytic mechanism of ChaP was explored by performing quantum mechanical/molecular mechanical (QM/MM) calculations. Our calculation results reveal that ChaP uses the proximal oxygen in iron-coordinated HOO^- to attack the carbonyl carbon of the substrate, whereas the previous proposal that Asp49 acts as a base to deprotonate the iron-coordinated HOO^- to generate O₂^2- is unlikely. In the first stage reaction, owing to the coordination of the substrate with iron, the substrate is activated by accepting an electron from iron and the resulting oxy-radical intermediate formed by O-O cleavage can easily trigger the ring rearrangement. In the final decarboxylation, the phenolic anion of the substrate cooperatively accepts the proton of iron-coordinated HOO^- to facilitate the attack of the distal oxygen, and the proton-coupled electron transfer (PCET) from the substrate to the Fe^IV═O plays a key role for the decarboxylation. These findings may provide useful information for understanding the ChaP-catalyzed oxidative rearrangement and other flavin-dependent non-heme dioxygenases.