Understanding Complex Mechanisms of Enzyme Reactivity: The Case of Limonene-1,2-Epoxide Hydrolases

Reaction Mechanisms and Kinetics of Nanozymes: Insights from Theory and Computation

"Nanozymes" usually refers to inorganic nanomaterials with enzyme-like catalytic activities. The research into nanozymes is one of the hot topics on the horizon of interdisciplinary science involving materials, chemistry, and biology. Although great progress has been made in the design, synthesis, characterization, and application of nanozymes, the study of the underlying microscopic mechanisms and kinetics is still not straightforward. Density functional theory (DFT) calculations compute the potential energy surfaces along the reaction coordinates for chemical reactions, which can give atomistic-level insights into the micro-mechanisms and kinetics for nanozymes. Therefore, DFT calculations have been playing an increasingly important role in exploring the mechanisms and kinetics for nanozymes in the past years. The calculations either predict the microscopic details for the catalytic processes to complement the experiments or further develop theoretical models to depict the physicochemical rules. In this review, the corresponding research progress is summarized. Particularly, the review focuses on the computational studies that closely interplay with the experiments. The relevant experimental results without DFT calculations will be also briefly discussed to offer a historic overview of how the computations promote the understanding of the microscopic mechanisms and kinetics of nanozymes.

Efficient decolorization of melanoidin in raw molasses wastewater by thermophilic esterase in actual extreme conditions

This work was developed to explore the versatility of thermophilic esterase for decolorizing raw molasses wastewater at high temperature and acidic pH. Combining covalent crosslinking method with deep eutectic solvent, a thermophilic esterase from Pyrobaculum calidifontis was immobilized on chitosan/macroporous resin composite carrier. The application of this immobilized thermophilic esterase eliminated 92.35% of colorants in raw molasses wastewater, achieving maximal decolorization efficiency across all the enzymes tested. Strikingly, this immobilized thermophilic esterase was capable of engaging in continuous activity for a 5-day period while removing 76.23% of pigments from samples. It effectively and continuously eliminated BOD₅ and COD, effectively and directly facilitating raw molasses wastewater decolorization under extreme conditions more readily than control group. In addition, this thermophilic esterase was believed to achieve decolorization through an addition reaction that disrupted conjugated system of melanoidins. Together, these results highlight an efficient and practical means of achieving enzyme-based molasses wastewater decolorization.

Exploring the Catalytic Mechanism of the RNA Cap Modification by nsp16-nsp10 Complex of SARS-CoV-2 through a QM/MM Approach

The inhibition of key enzymes that may contain the viral replication of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) have assumed central importance in drug discovery projects. Nonstructural proteins (nsps) are essential for RNA capping and coronavirus replication since it protects the virus from host innate immune restriction. In particular, nonstructural protein 16 (nsp16) in complex with nsp10 is a Cap-0 binding enzyme. The heterodimer formed by nsp16-nsp10 methylates the 5′-end of virally encoded mRNAs to mimic cellular mRNAs and thus it is one of the enzymes that is a potential target for antiviral therapy. In this study, we have evaluated the mechanism of the 2′-O methylation of the viral mRNA cap using hybrid quantum mechanics/molecular mechanics (QM/MM) approach. It was found that the calculated free energy barriers obtained at M062X/6-31+G(d,p) is in agreement with experimental observations. Overall, we provide a detailed molecular analysis of the catalytic mechanism involving the 2′-O methylation of the viral mRNA cap and, as expected, the results demonstrate that the TS stabilization is critical for the catalysis.

Exploration of the structural requirements of Aurora Kinase B inhibitors by a combined QSAR, modelling and molecular simulation approach

Aurora kinase B plays an important role in the cell cycle to orchestrate the mitotic process. The amplification and overexpression of this kinase have been implicated in several human malignancies. Therefore, Aurora kinase B is a potential drug target for anticancer therapies. Here, we combine atom-based 3D-QSAR analysis and pharmacophore model generation to identify the principal structural features of acylureidoindolin derivatives that could potentially be responsible for the inhibition of Aurora kinase B. The selected CoMFA and CoMSIA model showed significant results with cross-validation values (q²) of 0.68, 0.641 and linear regression values (r²) of 0.971, 0.933 respectively. These values support the statistical reliability of our model. A pharmacophore model was also generated, incorporating features of reported crystal complex structures of Aurora kinase B. The pharmacophore model was used to screen commercial databases to retrieve potential lead candidates. The resulting hits were analyzed at each stage for diversity based on the pharmacophore model, followed by molecular docking and filtering based on their interaction with active site residues and 3D-QSAR predictions. Subsequently, MD simulations and binding free energy calculations were performed to test the predictions and to characterize interactions at the molecular level. The results suggested that the identified compounds retained the interactions with binding residues. Binding energy decomposition identified residues Glu155, Trp156 and Ala157 of site B and Leu83 and Leu207 of site C as major contributors to binding affinity, complementary to 3D-QSAR results. To best of our knowledge, this is the first comparison of WaterSwap field and 3D-QSAR maps. Overall, this integrated strategy provides a basis for the development of new and potential AK-B inhibitors and is applicable to other protein targets.

Computer-aided understanding and engineering of enzymatic selectivity

Enzymes offering chemo-, regio-, and stereoselectivity enable the asymmetric synthesis of high-value chiral molecules. Unfortunately, the drawback that naturally occurring enzymes are often inefficient or have undesired selectivity toward non-native substrates hinders the broadening of biocatalytic applications. To match the demands of specific selectivity in asymmetric synthesis, biochemists have implemented various computer-aided strategies in understanding and engineering enzymatic selectivity, diversifying the available repository of artificial enzymes. Here, given that the entire asymmetric catalytic cycle, involving precise interactions within the active pocket and substrate transport in the enzyme channel, could affect the enzymatic efficiency and selectivity, we presented a comprehensive overview of the computer-aided workflow for enzymatic selectivity. This review includes a mechanistic understanding of enzymatic selectivity based on quantum mechanical calculations, rational design of enzymatic selectivity guided by enzyme-substrate interactions, and enzymatic selectivity regulation via enzyme channel engineering. Finally, we discussed the computational paradigm for designing enzyme selectivity in silico to facilitate the advancement of asymmetric biosynthesis.

Tuning Proton Transfer Thermodynamics in SARS-CoV-2 Main Protease: Implications for Catalysis and Inhibitor Design

The catalytic reaction in SARS-CoV-2 main protease is activated by a proton transfer (PT) from Cys145 to His41. The same PT is likely also required for the covalent binding of some inhibitors. Here we use a multiscale computational approach to investigate the PT thermodynamics in the apo enzyme and in complex with two potent inhibitors, N3 and the α-ketoamide 13b. We show that with the inhibitors the free energy cost to reach the charge-separated state of the active-site dyad is lower, with N3 inducing the most significant reduction. We also show that a few key sites (including specific water molecules) significantly enhance or reduce the thermodynamic feasibility of the PT reaction, with selective desolvation of the active site playing a crucial role. The approach presented is a cost-effective procedure to identify the enzyme regions that control the activation of the catalytic reaction and is thus also useful to guide the design of inhibitors.

Decolorization of molasses alcohol wastewater by thermophilic hydrolase with practical application value

The aim of this work was to explore the ability of cutinase in the decolorization of molasses wastewater. Thermophilic cutinase from Thermobifida alba eliminated 76.1-78.2% of colorants and exhibited the highest decolorization efficiency amongst all of the cutinases tested. Cutinase from Thermobifida alba was immobilized on an affordable and efficient modified chitosan carrier and achieved a decolorization yield of 79.3-81.2%. This cutinase removed 66.3-71.1% of pigments and lasted continuously for 5 days. Importantly, it was also shown to continuously and effectively remove COD and BOD₅. Compared to other enzymes, the immobilized cutinase from Thermobifida alba had the advantage of being low-cost and had a high expression level and activity. The results confirmed the decolorization occurred by destroying the conjugated system of melanoidins via an addition reaction by cutinase from Thermobifida alba. Thus, this study contributes a more practical and efficient approach to enzymatic decolorization of molasses wastewater.

Structure-function investigation of 3-methylaspartate ammonia lyase reveals substrate molecular determinants for the deamination reaction

The enzymatic reactions leading to the deamination of β-lysine, lysine, or 2-aminoadipic acid are of great interest for the metabolic conversion of lysine to adipic acid. Enzymes able to carry out these reactions are not known, however ammonia lyases (EC 4.3.1.-) perform deamination on a wide range of substrates. We have studied 3-methylaspartate ammonia lyase (MAL, EC 4.3.1.2) as a potential candidate for protein engineering to enable deamination towards β-lysine, that we have shown to be a competitive inhibitor of MAL. We have characterized MAL activity, binding and inhibition properties on six different compounds that would allow to define the molecular determinants necessary for MAL to deaminate our substrate of interest. Docking calculations showed that β-lysine as well as the other compounds investigated could fit spatially into MAL catalytic pocket, although they probably are weak or very transient binders and we identified molecular determinants involved in the binding of the substrate. The hydrophobic interactions formed by the methyl group of 3-methylaspartic acid, together with the presence of the amino group on carbon 2, play an essential role in the appropriate binding of the substrate. The results showed that β-lysine is able to fit and bind in MAL catalytic pocket and can be potentially converted from inhibitor to substrate of MAL upon enzyme engineering. The characterization of the binding and inhibition properties of the substrates tested here provide the foundation for future and more extensive studies on engineering MAL that could lead to a MAL variant able to catalyse this challenging deamination reaction.

Assessment of the Cruzain Cysteine Protease Reversible and Irreversible Covalent Inhibition Mechanism

Reversible and irreversible covalent ligands are advanced cysteine protease inhibitors in the drug development pipeline. K777 is an irreversible inhibitor of cruzain, a necessary enzyme for the survival of the Trypanosoma cruzi (T. cruzi) parasite, the causative agent of Chagas disease. Despite their importance, irreversible covalent inhibitors are still often avoided due to the risk of adverse effects. Herein, we replaced the K777 vinyl sulfone group with a nitrile moiety to obtain a reversible covalent inhibitor (Neq0682) of cysteine protease. Then, we used advanced experimental and computational techniques to explore details of the inhibition mechanism of cruzain by reversible and irreversible inhibitors. The isothermal titration calorimetry (ITC) analysis shows that inhibition of cruzain by an irreversible inhibitor is thermodynamically more favorable than by a reversible one. The hybrid Quantum Mechanics/Molecular Mechanics (QM/MM) and Molecular Dynamics (MD) simulations were used to explore the mechanism of the reaction inhibition of cruzain by K777 and Neq0682. The calculated free energy profiles show that the Cys25 nucleophilic attack and His162 proton transfer occur in a single step for a reversible inhibitor and two steps for an irreversible covalent inhibitor. The hybrid QM/MM calculated free energies for the inhibition reaction correspond to -26.7 and -5.9 kcal mol^-1 for K777 and Neq0682 at the MP2/MM level, respectively. These results indicate that the ΔG of the reaction is very negative for the process involving K777, consequently, the covalent adduct cannot revert to a noncovalent protein-ligand complex, and its binding tends to be irreversible. Overall, the present study provides insights into a covalent inhibition mechanism of cysteine proteases.

Computational Design of Enantiocomplementary Epoxide Hydrolases for Asymmetric Synthesis of Aliphatic and Aromatic Diols

The use of enzymes in preparative biocatalysis often requires tailoring enzyme selectivity by protein engineering. Herein we explore the use of computational library design and molecular dynamics simulations to create variants of limonene epoxide hydrolase that produce enantiomeric diols from meso‐epoxides. Three substrates of different sizes were targeted: cis‐2,3‐butene oxide, cyclopentene oxide, and cis‐stilbene oxide. Most of the 28 designs tested were active and showed the predicted enantioselectivity. Excellent enantioselectivities were obtained for the bulky substrate cis‐stilbene oxide, and enantiocomplementary mutants produced (S,S)‐ and (R,R)‐stilbene diol with >97 % enantiomeric excess. An (R,R)‐selective mutant was used to prepare (R,R)‐stilbene diol with high enantiopurity (98 % conversion into diol, >99 % ee). Some variants displayed higher catalytic rates (k_cat) than the original enzyme, but in most cases K_M values increased as well. The results demonstrate the feasibility of computational design and screening to engineer enantioselective epoxide hydrolase variants with very limited laboratory screening.

Molecular Dynamics Simulation Framework to Probe the Binding Hypothesis of CYP3A4 Inhibitors

The Cytochrome P450 family of heme-containing proteins plays a major role in catalyzing phase I metabolic reactions, and the CYP3A4 subtype is responsible for the metabolism of many currently marketed drugs. Additionally, CYP3A4 has an inherent affinity for a broad spectrum of structurally diverse chemical entities, often leading to drug-drug interactions mediated by the inhibition or induction of the metabolic enzyme. The current study explores the binding of selected highly efficient CYP3A4 inhibitors by docking and molecular dynamics (MD) simulation protocols and their binding free energy calculated using the WaterSwap method. The results indicate the importance of binding pocket residues including Phe57, Arg105, Arg106, Ser119, Arg212, Phe213, Thr309, Ser312, Ala370, Arg372, Glu374, Gly481 and Leu483 for interaction with CYP3A4 inhibitors. The residue-wise decomposition of the binding free energy from the WaterSwap method revealed the importance of binding site residues Arg106 and Arg372 in the stabilization of all the selected CYP3A4-inhibitor complexes. The WaterSwap binding energies were further complemented with the MM(GB/PB)SA results and it was observed that the binding energies calculated by both methods do not differ significantly. Overall, our results could guide towards the use of multiple computational approaches to achieve a better understanding of CYP3A4 inhibition, subsequently leading to the design of highly specific and efficient new chemical entities with suitable ADMETox properties and reduced side effects.

The reaction mechanism of Zika virus NS2B/NS3 serine protease inhibition by dipeptidyl aldehyde: a QM/MM study

Zika virus (ZIKV) infection has become a global public health problem, associated with microcephaly in newborns and Guillain-Barré syndrome in adults. Currently, there are no commercially available anti-ZIKV drugs. The viral protease NS2B/NS3, which is involved in viral replication and maturation, is a potential drug target. Peptidomimetic aldehyde inhibitors bind covalently to the catalytic S135 of the NS3 protease. Here, we apply hybrid quantum mechanics/molecular mechanics (QM/MM) free-energy simulations at the PDDG-PM3/ff14SB level to investigate the inhibition mechanism of the ZIKV protease by a dipeptidyl aldehyde inhibitor (acyl-KR-aldehyde). The results show that proton transfer from the catalytic S135 to H51 occurs in concert with nucleophilic addition on the aldehyde warhead by S135. The anionic covalent complex between the dipeptidyl aldehyde and the ZIKV protease is analogous to the tetrahedral intermediate for substrate hydrolysis. Spontaneous protonation by H51 forms the hemiacetal. In addition, we use correlated ab initio QM/MM potential energy path calculations at levels up to LCCSD(T)/(aug)-cc-pVTZ to obtain accurate potential energy profiles of the reaction, which also support a concerted mechanism. These results provide detailed insight into the mechanism of ZIKV protease inhibition by a peptidyl aldehyde inhibitor, which will guide in the design of inhibitors.

Projector-Based Embedding Eliminates Density Functional Dependence for QM/MM Calculations of Reactions in Enzymes and Solution

Combined quantum mechanics/molecular mechanics (QM/MM) methods are increasingly widely utilized in studies of reactions in enzymes and other large systems. Here, we apply a range of QM/MM methods to investigate the Claisen rearrangement of chorismate to prephenate, in solution, and in the enzyme chorismate mutase. Using projector-based embedding in a QM/MM framework, we apply treatments up to the CCSD(T) level. We test a range of density functional QM/MM methods and QM region sizes. The results show that the calculated reaction energetics are significantly more sensitive to the choice of density functional than they are to the size of the QM region in these systems. Projector-based embedding of a wave function method in DFT reduced the 13 kcal/mol spread in barrier heights calculated at the DFT/MM level to a spread of just 0.3 kcal/mol, essentially eliminating dependence on the functional. Projector-based embedding of correlated ab initio methods provides a practical method for achieving high accuracy for energy profiles derived from DFT and DFT/MM calculations for reactions in condensed phases.

Visualizing protein-ligand binding with chemical energy-wise decomposition (CHEWD): application to ligand binding in the kallikrein-8 S1 Site

Kallikrein-8, a serine protease, is a target for structure-based drug design due to its therapeutic potential in treating Alzheimer's disease and is also useful as a biomarker in ovarian cancer. We present a binding assessment of ligands to kallikrein-8 using a residue-wise decomposition of the binding energy. Binding of four putative inhibitors of kallikrein-8 is investigated through molecular dynamics simulation and ligand binding energy evaluation with two methods (MM/PBSA and WaterSwap). For visualization of the residue-wise decomposition of binding energies, chemical energy-wise decomposition or CHEWD is introduced as a plugin to UCSF Chimera and Pymol. CHEWD allows easy comparison between ligands using individual residue contributions to the binding energy. Molecular dynamics simulations indicate one ligand binds stably to the kallikrein-8 S1 binding site. Comparison with other members of the kallikrein family shows that residues responsible for binding are specific to kallikrein-8. Thus, ZINC02927490 is a promising lead for development of novel kallikrein-8 inhibitors.