QM/MM Study of a VIM-1 Metallo-β-Lactamase Enzyme: The Catalytic Reaction Mechanism

Mechanism of Nucleic Acid Phosphodiester Bond Cleavage by Human Endonuclease V: MD and QM/MM Calculations Reveal a Versatile Metal Dependence

Human endonuclease V (EndoV) catalytically removes deaminated nucleobases by cleaving the phosphodiester bond as part of RNA metabolism. Despite being implicated in several diseases (cancers, cardiovascular diseases, and neurological disorders) and potentially being a useful tool in biotechnology, details of the human EndoV catalytic pathway remain unclear due to limited experimental information beyond a crystal structure of the apoenzyme and select mutational data. Since a mechanistic understanding is critical for further deciphering the central roles and expanding applications of human EndoV in medicine and biotechnology, molecular dynamics (MD) simulations and quantum mechanics/molecular mechanics (QM/MM) calculations were used to unveil the atomistic details of the catalytic pathway. Due to controversies surrounding the number of metals required for nuclease activity, enzyme-substrate models with different numbers of active site metals and various metal-substrate binding configurations were built based on structural data for other nucleases. Subsequent MD simulations revealed the structure and stability of the human EndoV-substrate complex for a range of active site metal binding architectures. Four unique pathways were then characterized using QM/MM that vary in metal number (one versus two) and modes of substrate coordination [direct versus indirect (water-mediated)], with several mechanisms being fully consistent with experimental structural, kinetic, and mutational data for related nucleases, including members of the EndoV family. Beyond uncovering key roles for several active site amino acids (D240 and K155), our calculations highlight that while one metal is essential for human EndoV activity, the enzyme can benefit from using two metals due to the presence of two suitable metal binding sites. By directly comparing one- versus two-metal-mediated P-O bond cleavage reactions within the confines of the same active site, our work brings a fresh perspective to the "number of metals" controversy.

State-of-the-art local correlation methods enable affordable gold standard quantum chemistry for up to hundreds of atoms

In this feature, we review the current capabilities of local electron correlation methods up to the coupled cluster model with single, double, and perturbative triple excitations [CCSD(T)], which is a gold standard in quantum chemistry. The main computational aspects of the local method types are assessed from the perspective of applications, but the focus is kept on how to achieve chemical accuracy (i.e., <1 kcal mol-1 uncertainty), as well as on the broad scope of chemical problems made accessible. The performance of state-of-the-art methods is also compared, including the most employed DLPNO and, in particular, our local natural orbital (LNO) CCSD(T) approach. The high accuracy and efficiency of the LNO method makes chemically accurate CCSD(T) computations accessible for molecules of hundreds of atoms with resources affordable to a broad computational community (days on a single CPU and 10-100 GB of memory). Recent developments in LNO-CCSD(T) enable systematic convergence and robust error estimates even for systems of complicated electronic structure or larger size (up to 1000 atoms). The predictive power of current local CCSD(T) methods, usually at about 1-2 order of magnitude higher cost than hybrid density functional theory (DFT), has become outstanding on the palette of computational chemistry applicable for molecules of practical interest. We also review more than 50 LNO-based and other advanced local-CCSD(T) applications for realistic, large systems across molecular interactions as well as main group, transition metal, bio-, and surface chemistry. The examples show that properly executed local-CCSD(T) can contribute to binding, reaction equilibrium, rate constants, etc. which are able to match measurements within the error estimates. These applications demonstrate that modern, open-access, and broadly affordable local methods, such as LNO-CCSD(T), already enable predictive computations and atomistic insight for complicated, real-life molecular processes in realistic environments.

Probing the interaction between the metallo-β-lactamase SMB-1 and ampicillin by multispectral approaches combined with molecular dynamics

Metallo-β-lactamases (MβLs) hydrolyze and inactivate β-lactam antibiotics, are a pivotal mechanism conferring resistance against bacterial infections. SMB-1, a novel B3 subclass of MβLs from Serratia marcescens could deactivate almost all β-lactam antibiotics including ampicillin (AMP), which has posed a serious threat to public health. To illuminate the mechanism of recognition and interaction between SMB-1 and AMP, various fluorescence spectroscopy techniques and molecular dynamics simulation were employed. The results of quenching spectroscopy unraveled that AMP could make SMB-1 fluorescence quenching that mechanism was the static quenching; the synchronous and three-dimensional fluorescence spectra validated that the microenvironment and conformation of SMB-1 were altered after interaction with AMP. The molecular dynamics results demonstrated that the whole AMP enters the binding pocket of SMB-1, even though with a relatively bulky R1 side chain. Loop1 and loop2 in SMB-1 undergo significant fluctuations, and α2 (71-73) and local α5 (186-188) were turned into random coils, promoting zinc ion exposure consistent with circular dichroism spectroscopy results. The binding between them was driven by a combination of enthalpy and entropy changes, which was dominated by electrostatic force in agreement with the fluorescence observations. The present study brings structural insights and solid foundations for the design of new substrates for β-lactamases and the development of effective antibiotics that are resistant to superbugs.

Engineering Glucosamine-6-Phosphate Synthase to Achieve Efficient One-Step Biosynthesis of Glucosamine

As an important functional monosaccharide, glucosamine (GlcN) is widely used in fields such as medicine, food nutrition, and health care. Here, we report a distinct GlcN biosynthesis method that utilizes engineered Bacillus subtilis glucosamine-6-phosphate synthase (BsGlmS) to convert D-fructose to directly generate GlcN. The best variant obtained by using a combinatorial active-site saturation test/iterative saturation mutagenesis (CAST/ISM) strategy was a quadruple mutant S596D/V597G/S347H/G299Q (BsGlmS-BK19), which has a catalytic activity 1736-fold that of the wild type toward D-fructose. Upon using mutant BK19 as a whole-cell catalyst, D-fructose was converted into GlcN with 65.32% conversion in 6 h, whereas the wild type only attained a conversion rate of 0.31% under the same conditions. Molecular docking and molecular dynamics simulations were implemented to provide insights into the mechanism underlying the enhanced activity of BK19. Importantly, the BsGlmS-BK19 variant specifically catalyzes D-fructose without the need for phosphorylated substrates, representing a significant advancement in GlcN biosynthesis.

Cysteine Radical and Glutamate Collaboratively Enable C-H Bond Activation and C-N Bond Cleavage in a Glycyl Radical Enzyme HplG

Hydroxyprolines are abundant in nature and widely utilized by many living organisms. Isomerization of trans-4-hydroxy-d-proline (t4D-HP) to generate 2-amino-4-ketopentanoate has been found to need a glycyl radical enzyme HplG, which catalyzes the cleavage of the C-N bond, while dehydration of trans-4-hydroxy-l-proline involves a homologous enzyme of HplG. Herein, molecular dynamics simulations and quantum mechanics/molecular mechanics (QM/MM) calculations are employed to understand the reaction mechanism of HplG. Two possible reaction pathways of HplG have been explored to decipher the origin of its chemoselectivity. The QM/MM calculations reveal that the isomerization proceeds via an initial hydrogen shift from the Cγ site of t4D-HP to a catalytic cysteine radical, followed by cleavage of the Cδ-N bond in t4D-HP to form a radical intermediate that captures a hydrogen atom from the cysteine. Activation of the Cδ-H bond in t4D-HP to bring about dehydration of t4D-HP possesses an extremely high energy barrier, thus rendering the dehydration pathway implausible in HplG. On the basis of the current calculations, conserved residue Glu429 plays a pivotal role in the isomerization pathway: the hydrogen bonding between it and t4D-HP weakens the hydroxyalkyl Cγ-Hγ bond, and it acts as a proton acceptor to trigger the cleavage of the C-N bond in t4D-HP. Our current QM/MM calculations rationalize the origin of the experimentally observed chemoselectivity of HplG and propose an H-bond-assisted bond activation strategy in radical-containing enzymes. These findings have general implications on radical-mediated enzymatic catalysis and expand our understanding of how nature wisely and selectively activates the C-H bond to modulate catalytic selectivity.

Effect of lactic acid bacteria fermentation on cow milk allergenicity and antigenicity: A review

Cow milk is a major allergenic food. The potential prevention and treatment effects of lactic acid bacteria (LAB)-fermented dairy products on allergic symptoms have garnered considerable attention. Cow milk allergy (CMA) is mainly attributed to extracellular and/or cell envelope proteolytic enzymes with hydrolysis specificity. Numerous studies have demonstrated that LAB prevents the risk of allergies by modulating the development and regulation of the host immune system. Specifically, LAB and its effectors can enhance intestinal barrier function and affect immune cells by interfering with humoral and cellular immunity. Fermentation hydrolysis of allergenic epitopes is considered the main mechanism of reducing CMA. This article reviews the linear epitopes of allergens in cow milk and the effect of LAB on these allergens and provides insight into the means of predicting allergenic epitopes by conventional laboratory analysis methods combined with molecular simulation. Although LAB can reduce CMA in several ways, the mechanism of action remains partially clarified. Therefore, this review additionally attempts to summarize the main mechanism of LAB fermentation to provide guidance for establishing an effective preventive and treatment method for CMA and serve as a reference for the screening, research, and application of LAB-based intervention.

Impacts of QM region sizes and conformation numbers on modelling enzyme reactions: a case study of polyethylene terephthalate hydrolase

A quantum mechanics/molecular mechanics (QM/MM) approach is a broadly used tool in computational enzymology. Treating the QM region with a high-level DFT method is one of the important branches. Here, taking leaf-branch compost cutinase-catalyzed polyethylene terephthalate depolymerization as an example, the convergence behavior of energy barriers as well as key structural and charge features with respect to the size of the QM region (up to 1000 atoms) is systematically investigated. BP86/6-31G(d)//CHARMM and M06-2X/6-311G(d,p)//CHARMM level of theories were applied for geometry optimizations and single-point energy calculations, respectively. Six independent enzyme conformations for all the four catalytic steps (steps (i)-(iv)) were considered. Most of the twenty-four cases show that at least 500 QM atoms are needed while only two rare cases show that ∼100 QM atoms are sufficient for convergence when only a single conformation was considered. This explains why most previous studies showed that 500 or more QM atoms are required while a few others showed that ∼100 QM atoms are sufficient for DFT/MM calculations. More importantly, average energy barriers and key structural/charge features from six conformations show an accelerated convergence than that in a single conformation. For instance, to reach energy barrier convergence (within 2.0 kcal mol-1) for step (ii), only ∼100 QM atoms are required if six conformations are considered while 500 or more QM atoms are needed with a single conformation. The convergence is accelerated to be more rapid if hundreds and thousands of conformations were considered, which aligns with previous findings that only several dozens of QM atoms are required for convergence with semi-empirical QM/MM MD simulations.

Unveiling the Catalytic Mechanism of a Processive Metalloaminopeptidase

Intracellular leucine aminopeptidases (PepA) are metalloproteases from the family M17. These enzymes catalyze peptide bond cleavage, removing N-terminal residues from peptide and protein substrates, with consequences for protein homeostasis and quality control. While general mechanistic studies using model substrates have been conducted on PepA enzymes from various organisms, specific information about their substrate preferences and promiscuity, choice of metal, activation mechanisms, and the steps that limit steady-state turnover remain unexplored. Here, we dissected the catalytic and chemical mechanisms of PaPepA: a leucine aminopeptidase from Pseudomonas aeruginosa. Cleavage assays using peptides and small-molecule substrate mimics allowed us to propose a mechanism for catalysis. Steady-state and pre-steady-state kinetics, pH rate profiles, solvent kinetic isotope effects, and biophysical techniques were used to evaluate metal binding and activation. This revealed that metal binding to a tight affinity site is insufficient for enzyme activity; binding to a weaker affinity site is essential for catalysis. Progress curves for peptide hydrolysis and crystal structures of free and inhibitor-bound PaPepA revealed that PaPepA cleaves peptide substrates in a processive manner. We propose three distinct modes for activity regulation: tight packing of PaPepA in a hexameric assembly controls substrate length and reaction processivity; the product leucine acts as an inhibitor, and the high concentration of metal ions required for activation limits catalytic turnover. Our work uncovers catalysis by a metalloaminopeptidase, revealing the intricacies of metal activation and substrate selection. This will pave the way for a deeper understanding of metalloenzymes and processive peptidases/proteases.

Deciphering the Catalytic Mechanism of Virginiamycin B Lyase with Multiscale Methods and Molecular Dynamics Simulations

Due to the emergence of antibiotic resistance, the need to explore novel antibiotics and/or novel strategies to counter antibiotic resistance is of utmost importance. In this work, we explored the molecular and mechanistic details of the degradation of a streptogramin B antibiotic by virginiamycin B (Vgb) lyase of Staphylococcus aureus using classical molecular dynamics simulations and multiscale quantum mechanics/molecular mechanics methods. Our results were in line with available experimental kinetic information. Although we were able to identify a stepwise mechanism, in the wild-type enzyme, the intermediate is short-lived, showing a small barrier to decay to the product state. The impact of point mutations on the reaction was also assessed, showing not only the importance of active site residues to the reaction catalyzed by Vgb lyase but also of near positive and negative residues surrounding the active site. Using molecular dynamics simulations, we also predicted the most likely protonation state of the 3-hydroxypicolinic moiety of the antibiotic and the impact of mutants on antibiotic binding. All this information will expand our understanding of linearization reactions of cyclic antibiotics, which are crucial for the development of novel strategies that aim to tackle antibiotic resistance.

Catalytic Reaction Mechanism of Glyoxalase II: A Quantum Mechanics/Molecular Mechanics Study

Methylglyoxal (MG) is a reactive and toxic compound produced in carbohydrate, lipid, and amino acid metabolism. The glyoxalase system is the main detoxifying route for MG and consists of two enzymes, glyoxalase I (GlxI) and glyoxalase II (GlxII). GlxI catalyzes the formation of S-d-lactoylglutathione from hemithioacetal, and GlxII converts this intermediate to d-lactate. A relationship between the glyoxalase system and some diseases like diabetes has been shown, and inhibiting enzymes of this system may be an effective means of controlling certain diseases. A detailed understanding of the reaction mechanism of an enzyme is essential to the rational design of competitive inhibitors. In this work, we use quantum mechanics/molecular mechanics (QM/MM) calculations and energy refinement utilizing the big-QM and QM/MM thermodynamic cycle perturbation methods to propose a mechanism for the GlxII reaction that starts with a nucleophilic attack of the bridging OH^- group on the substrate. The coordination of the substrate to the Zn ions places its electrophilic center close to the hydroxide group, enabling the reaction to proceed. Our estimated reaction energies are in excellent agreement with experimental data, thus demonstrating the reliability of our approach and the proposed mechanism. Additionally, we examined alternative protonation states of Asp-29, Asp-58, Asp-134, and the bridging hydroxide ion in the catalytic process. However, these give less favorable reactions, a poorer reproduction of the crystal structure geometry of the active site, and higher root-mean-squared deviations of the active site residues in molecular dynamics simulations.

Interplay between the Enamine and Imine Forms of the Hydrolyzed Imipenem in the Active Sites of Metallo-β-lactamases and in Water Solution

Deactivation of the β-lactam antibiotics in the active sites of the β-lactamases is among the main mechanisms of bacterial antibiotic resistance. As drugs of last resort, carbapenems are efficiently hydrolyzed by metallo-β-lactamases, presenting a serious threat to human health. Our study reveals mechanistic aspects of the imipenem hydrolysis by bizinc metallo-β-lactamases, NDM-1 and L1, belonging to the B1 and the B3 subclasses, respectively. The results of QM(PBE0-D3/6-31G**)/MM simulations show that the enamine product with the protonated nitrogen atom is formed as the major product in NDM-1 and as the only product in the L1 active site. In NDM-1, there is also another reaction pathway that leads to the formation of the (S)-enantiomer of the imine form of the hydrolyzed imipenem; this process occurs with the higher energy barriers. The absence of the second pathway in L1 is due to the different amino acid composition of the active site loop. In L1, the hydrophobic Pro226 residue is located above the pyrroline ring of imipenem that blocks protonation of the carbon atom. Electron density analysis is performed at the stationary points to compare reaction pathways in L1 and NDM-1. Tautomerization from the enamine to the imine form likely happens in solution after the dissociation of the hydrolyzed imipenem from the active site of the enzyme. Classical molecular dynamics simulations of the hydrolyzed imipenem in solution, both with the neutral enamine and the negatively charged N-C₂-C₃ fragment, demonstrate a huge diversity of conformations. The vast majority of conformations blocks the C₃-atom from the side required for the (S)-imine formation upon tautomerization. Thus, according to our calculations, formation of the (R)-imine is more likely. QM(PBE0-D3/6-31G**)/MM molecular dynamics simulations of the hydrolyzed imipenem with the negatively charged N-C₂-C₃ fragment followed by the Laplacian bond order analysis demonstrate that the N═C₂-C₃^- resonance structure is the most pronounced that facilitates formation of the imine form. The proposed mechanism of the enzymatic enamine formation and its subsequent tautomerization to the imine form in solution is in agreement with the recent spectroscopic and NMR studies.

The Impact of DFT Functional, Cluster Model Size, and Implicit Solvation on the Structural Description of Single-Metal-Mediated DNA Phosphodiester Bond Cleavage: The Case Study of APE1

Phosphodiester bond hydrolysis in nucleic acids is a ubiquitous reaction that can be facilitated by enzymes called nucleases, which often use metal ions to achieve catalytic function. While a two-metal-mediated pathway has been well established for many enzymes, there is growing support that some enzymes require only one metal for the catalytic step. Using human apurinic/apyrimidinic endonuclease (APE1) as a prototypical example and cluster models, this study clarifies the impact of DFT functional, cluster model size, and implicit solvation on single-metal-mediated phosphodiester bond cleavage and provides insight into how to efficiently model this chemistry. Initially, a model containing 69 atoms built from a high-resolution X-ray crystal structure is used to explore the reaction pathway mapped by a range of DFT functionals and basis sets, which provides support for the use of standard functionals (M06-2X and B3LYP-D3) to study this reaction. Subsequently, systematically increasing the model size to 185 atoms by including additional amino acids and altering residue truncation points highlights that small models containing only a few amino acids or β carbon truncation points introduce model strains and lead to incorrect metal coordination. Indeed, a model that contains all key residues (general base and acid, residues that stabilize the substrate, and amino acids that maintain the metal coordination) is required for an accurate structural depiction of the one-metal-mediated phosphodiester bond hydrolysis by APE1, which results in 185 atoms. The additional inclusion of the broader enzyme environment through continuum solvation models has negligible effects. The insights gained in the present work can be used to direct future computational studies of other one-metal-dependent nucleases to provide a greater understanding of how nature achieves this difficult chemistry.

A Molecular Mechanics Energy Partitioning Software for Biomolecular Systems

The partitioning of the molecular mechanics (MM) energy in calculations involving biomolecular systems is important to identify the source of major stabilizing interactions, e.g., in ligand–protein interactions, or to identify residues with considerable contributions in hybrid multiscale calculations, i.e., quantum mechanics/molecular mechanics (QM/MM). Here, we describe Energy Split, a software program to calculate MM energy partitioning considering the AMBER Hamiltonian and parameters. Energy Split includes a graphical interface plugin for VMD to facilitate the selection of atoms and molecules belonging to each part of the system. Energy Split is freely available at or can be easily installed through the VMD Store.

Quantum Mechanics/Molecular Mechanics Studies on the Catalytic Mechanism of a Novel Esterase (FmtA) of Staphylococcus aureus

FmtA is a novel esterase that shares the penicillin-binding protein (PBP) core structural folding but found to hydrolyze the removal of d-Ala from teichoic acids. Molecular docking, dynamics, and MM-GBSA of FmtA and its variants S127A, K130A, Y211A, D213A, and K130AY211A, in the presence or absence of wall teichoic acid (WTA), suggest that active site residues S127, K130, Y211, D213, N343, and G344 play a role in substrate binding. Quantum mechanics (QM)/molecular mechanics (MM) calculations reveal that during WTA catalysis, K130 deprotonates S127, and the nucleophilic S127 attacks the carbonyl carbon of d-Ala bound to WTA. The tetrahedral intermediate (TI) complex is stabilized by hydrogen bonding to the oxyanion holes. The TI complex displays a high energy gap and collapses to an energetically favorable acyl-enzyme complex.