Include Dispersion in Quantum Chemical Modeling of Enzymatic Reactions: The Case of Isoaspartyl Dipeptidase

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.

Benchmarking quantum mechanical methods for calculating reaction energies of reactions catalyzed by enzymes

To assess the accuracy of different quantum mechanical methods for biochemical modeling, the reaction energies of 20 small model reactions (chosen to represent chemical steps catalyzed by commonly studied enzymes) were calculated. The methods tested included several popular Density Functional Theory (DFT) functionals, second-order Møller Plesset perturbation theory (MP2) and its spin-component scaled variant (SCS-MP2), and coupled cluster singles and doubles and perturbative triples (CCSD(T)). Different basis sets were tested. CCSD(T)/aug-cc-pVTZ results for all 20 reactions were used to benchmark the other methods. It was found that MP2 and SCS-MP2 reaction energy calculation results are similar in quality to CCSD(T) (mean absolute error (MAE) of 1.2 and 1.3 kcal mol−1, respectively). MP2 calculations gave a large error in one case, and are more subject to basis set effects, so in general SCS-MP2 calculations are a good choice when CCSD(T) calculations are not feasible. Results with different DFT functionals were of reasonably good quality (MAEs of 2.5–5.1 kcal mol−1), whereas popular semi-empirical methods (AM1, PM3, SCC-DFTB) gave much larger errors (MAEs of 11.6–14.6 kcal mol−1). These results should be useful in guiding methodological choices and assessing the accuracy of QM/MM calculations on enzyme-catalyzed reactions.

Adsorption Properties of Hydrated Cr3+ Ions on Schiff-base Covalent Organic Frameworks: A DFT Study

Considering the superior physiochemical property, increasing efforts have been devoted to exploiting the covalent organic frameworks (COFs) materials on the environmental remediation of heavy metal ions. Water pollution caused by Cr³⁺ metal ions is of special concern for scientists and engineers. Notwithstanding all the former efforts made, it is surprising that very little is known about the interaction mechanisms between the hydrated Cr³⁺ metal ions and COF materials. In present context, density functional theory (DFT) method is used to elucidate geometric and electronic properties with the purpose of putting into theoretical perspective the application values and interaction mechanisms for COF materials on Cr³⁺ capture. The results showed that all the five selected Schiff-base COFs materials displayed good adsorption performance on Cr³⁺ removal while the phenazine-linked and imine-COFs possessed the most favorable adsorption capacity due to the optimal chemical units and frameworks. The hydration effect was found to play a two-side role in the adsorption process and interaction mechanisms, involving coordination, hydrogen bonds, as well as weak non-covalent interactions, have been illuminated to explain the observed different adsorption behaviors. This study provides a general guidance for the design and selection of efficient COF materials as high-capacity Cr³⁺ adsorbents.

Contributions of IQA electron correlation in understanding the chemical bond and non-covalent interactions

The quantum topological energy partitioning method Interacting Quantum Atoms (IQA) has been applied for over a decade resulting in an enlightening analysis of a variety of systems. In the last three years we have enriched this analysis by incorporating into IQA the two-particle density matrix obtained from Møller–Plesset (MP) perturbation theory. This work led to a new computational and interpretational tool to generate atomistic electron correlation and thus topologically based dispersion energies. Such an analysis determines the effects of electron correlation within atoms and between atoms, which covers both bonded and non-bonded “through -space” atom–atom interactions within a molecule or molecular complex. A series of papers published by us and other groups shows that the behavior of electron correlation is deeply ingrained in structural chemistry. Some concepts that were shown to be connected to bond correlation are bond order, multiplicity, aromaticity, and hydrogen bonding. Moreover, the concepts of covalency and ionicity were shown not to be mutually excluding but to both contribute to the stability of polar bonds. The correlation energy is considerably easier to predict by machine learning (kriging) than other IQA terms. Regarding the nature of the hydrogen bond, correlation energy presents itself in an almost contradicting way: there is much localized correlation energy in a hydrogen bond system, but its overall effect is null due to internal cancelation. Furthermore, the QTAIM delocalization index has a connection with correlation energy. We also explore the role of electron correlation in protobranching, which provides an explanation for the extra stabilization present in branched alkanes compared to their linear counterparts. We hope to show the importance of understanding the true nature of the correlation energy as the foundation of a modern representation of dispersion forces for ab initio, DFT, and force field calculations.

Mechanism and Inhibitor Exploration with Binuclear Mg Ketol-Acid Reductoisomerase: Targeting the Biosynthetic Pathway of Branched-Chain Amino Acids

Binuclear Mg ketol-acid reductoisomerase (KARI), which converts (S)-2-acetolactate into (R)-2,3-dihydroxyisovalerate, is responsible for the second step of the biosynthesis of branched-chain amino acids in plants and microorganisms and thus serves as a key inhibition target potentially without effects on mammals. Here, through the use of density functional calculations and a chemical model, the KARI-catalyzed reaction has been demonstrated to include the initial deprotonation of the substrate C2 hydroxy group, bridged by the two Mg ions, alkyl migration from the C2-alkoxide carbon atom to the C3-carbonyl carbon atom, and hydride transfer from a nicotinamide adenine dinucleotide phosphate [NAD(P)H] cofactor to C2. A dead-end mechanism with a hydride transferred to the C3 carbonyl group has been ruled out. The nucleophilicity (migratory aptitude) of the migrating carbon atom and the provision of additional negative charge to the di-Mg coordination sphere have significant effects on the steps of alkyl migration and hydride transfer, respectively. Other important mechanistic characteristics are also revealed. Inspired by the mechanism, an inhibitor (2-carboxylate-lactic acid) was designed and predicted by barrier analysis to be effective in inactivating KARI, hence probably enriching the antifungal and antibacterial library. Two types of slow substrate analogues (2-trihalomethyl acetolactic acids and 2-glutaryl lactic acid) were also found.

How does Mo-dependent perchlorate reductase work in the decomposition of oxyanions?

The mechanisms of Mo-dependent perchlorate reductase (PcrAB)-catalyzed decomposition of perchlorate, bromate, iodate, and nitrate were revealed by density functional calculations.

Revisiting the catalytic mechanism of Mo-Cu carbon monoxide dehydrogenase using QM/MM and DFT calculations

Previous density functional theory (DFT) studies have shown that the release of the produced carbon dioxide (CO2) from an active-site cluster is a thermodynamically or kinetically difficult step in the enzymatic carbon monoxide (CO) oxidation catalyzed by Mo-Cu carbon monoxide dehydrogenase (Mo-Cu CODH). To better understand the effect of the protein environment on this difficult CO2 release step as well as other reaction steps, we applied hybrid quantum mechanics and molecular mechanics (QM/MM) calculations to the Mo-Cu CODH enzyme. The results show that in the first step, the equatorial Mo[double bond, length as m-dash]O group in the active-site cluster attacks the nearby CO molecule bound to the Cu site. Afterward, a stable thiocarbonate intermediate is formed in which the CO2 molecule is embedded and the copper-S(μ-sulfido) bond is broken. A free CO2 molecule, i.e., the final product, is then released from the active-site cluster, not directly from the thiocarbonate intermediate but via a previously formed intermediate that also contains CO2 but retains the Cu-S(μ-sulfido) bond. In contrast to the previous DFT results, the calculated barrier for this process was low in our QM/MM calculations. An additional QM/MM analysis of the barrier height showed that the effect of the protein environment on this barrier lowering is not very large. We found that the reason for the low barrier obtained by QM/MM is that the barrier for CO2 release is already not high at the DFT level. These results allow us to conclude that the CO oxidation reaction passes through the formation of a thiocarbonate intermediate, and that the subsequent CO2 release is kinetically not difficult. Nevertheless, the protein environment has an important role to play in making the latter process thermodynamically favored. No low-barrier pathway for the product release could be obtained for the reaction of n-butylisocyanide, which is consistent with the experimental fact that n-butylisocyanide inhibits Mo-Cu CODH.

How does binuclear zinc amidohydrolase FwdA work in the initial step of methanogenesis: From formate to formyl-methanofuran

The initial step of methanogenesis is the fixation of CO₂ to formyl-methanofuran (formyl-MFR) catalyzed by formyl-MFR dehydrogenase, which can be divided into two half reactions. Herein, the second half reaction catalyzed by FwdA (formyl-methanofuran dehydrogenase subunit A), i.e., from formate to formyl-methanofuran, has been investigated using density functional theory and a chemical model based on the X-ray crystal structure. The calculations indicate that, compared with other well-known di-zinc hydrolases, the FwdA reaction employs a reverse mechanism, including the nucleophilic attack of MFR amine on formate carbon leading to a tetrahedral gem-diolate intermediate, two steps of proton transfer from amine to formate moieties assisted by the Asp385, and the CO bond dissociation to form the formyl-MFR product. The second step of proton transfer from the amine moiety to the Asp385 is rate-limiting with an overall barrier of 21.2 kcal/mol. The two zinc ions play an important role in stabilizing the transition states and intermediates, in particular the negative charge at the formate moiety originated from the nucleophilic attack of the MFR amine. The work here appends a crucial piece in the methanogenic mechanistics and advances the understanding of the global carbon cycle.

Asymmetric abstraction of two chemically-equivalent methylene hydrogens: significant enantioselectivity of endoperoxide presented by fumitremorgin B endoperoxidase

The unique enantioselectivity for an R-chiral endoperoxy ring in verruculogen biosynthesis originates from asymmetric abstraction of two chemically-equivalent methylene hydrogens.

Machine Learning of Dynamic Electron Correlation Energies from Topological Atoms

We present an innovative method for predicting the dynamic electron correlation energy of an atom or a bond in a molecule utilizing topological atoms. Our approach uses the machine learning method Kriging (Gaussian Process Regression with a non-zero mean function) to predict these dynamic electron correlation energy contributions. The true energy values are calculated by partitioning the MP2 two-particle density-matrix via the Interacting Quantum Atoms (IQA) procedure. To our knowledge, this is the first time such energies have been predicted by a machine learning technique. We present here three important proof-of-concept cases: the water monomer, the water dimer, and the van der Waals complex H₂···He. These cases represent the final step toward the design of a full IQA potential for molecular simulation. This final piece will enable us to consider situations in which dispersion is the dominant intermolecular interaction. The results from these examples suggest a new method by which dispersion potentials for molecular simulation can be generated.

Recent Trends in Quantum Chemical Modeling of Enzymatic Reactions

The quantum chemical cluster approach is a powerful method for investigating enzymatic reactions. Over the past two decades, a large number of highly diverse systems have been studied and a great wealth of mechanistic insight has been developed using this technique. This Perspective reviews the current status of the methodology. The latest technical developments are highlighted, and challenges are discussed. Some recent applications are presented to illustrate the capabilities and progress of this approach, and likely future directions are outlined.

From NAD+ to Nickel Pincer Complex: A Significant Cofactor Evolution Presented by Lactate Racemase

Lactate racemase (LarA), a new nickel enzyme discovered recently, catalyzes the racemization between d- and l-lactates with a novel nickel pincer cofactor (Ni-PTTMN) derived from nicotinic acid. In this study, by using DFT and a 200-atom active-site model, LarA is revealed to employ a modified proton-coupled hydride-transfer mechanism in which a hydride is transferred to a cofactor pyridine carbon from the substrate α-carbon along with proton transfer from the substrate hydroxy group to a histidine, and then moved back from the opposite side. Tyr294 and Lys298 provide significant acceleration effects by orientating substrates and stabilizing the negative charge developing at the substrate hydroxy oxygen. The barrier was determined to be 12.0 kcal mol^-1 , which reveals enhanced racemase activity relative to the LarA reaction using NAD⁺ -like cofactors. Compared with NAD⁺ , Ni-PTTMN has a stronger hydride-addition reactivity in moderate and high environmental polarity and may fit perfectly the moderately polar active site of LarA.

Effects of metal-ion replacement on pyrazinamidase activity: A quantum mechanical study

Pyrazinamidase (PZase), a metalloenzyme, is responsible for acidic modification of pyrazinamide (PZA), a drug used in tuberculosis treatment. The metal coordination site of the enzyme is able to coordinate various divalent metal cofactors. Previous experimental studies have demonstrated that metal ions, such as Co²⁺, Mn²⁺, and Zn²⁺, are able to reactivate metal-depleted PZase, while others including Cu²⁺, Fe²⁺, and Mg²⁺, cannot restore activity. In this study, we investigated binding of various metal ions to the metal coordination site (MCS) of the enzyme using quantum mechanical calculations. We calculated the metal-ligand (residue) binding energy and the atomic partial charges in the presence of various ions. The results indicated that the tendency of alkaline earth metals to bind to PZase MCS is very low and not suitable for enzyme structural and catalytic function. In contrast, Co²⁺ and Ni²⁺ ions have very high binding affinity and are favorable to the structural and functional properties of the enzyme. Furthermore, we observed that the rate at which Ni²⁺, Co²⁺ and Fe²⁺ ions in PZase MCS polarize the OH bond of coordinated water molecules is much higher than the polarization rate created by other ions. This finding suggests that the coordination of Ni²⁺, Co²⁺, or Fe²⁺ to PZase facilitates the deprotonation of coordinated water molecules to generate a nucleophile that catalyzes the enzymatic reaction.

The Decarboxylation of α,β-Unsaturated Acid Catalyzed by Prenylated FMN-Dependent Ferulic Acid Decarboxylase and the Enzyme Inhibition

Ferulic acid decarboxylase (Fdc1) is able to catalyze the decarboxylation of α,β-unsaturated acids using a novel cofactor, prenylated flavin mononucleotide (PrFMN). Using density functional theory calculations, we here have investigated the Fdc1 reaction mechanism with the substrate of α-methylcinnamic acid. It is demonstrated that Fdc1 employs a 1,3-dipolar cycloaddition mechanism involving four concerted steps, where the Glu282 acts as a crucial proton donor to protonate the α carbon (C_α). The last step, the decomposition of a pyrrolidine species, is rate-limiting with an overall barrier of 18.9 kcal mol^-1. Furthermore, when α-hydroxycinnamic acid is used, the Glu282 is found to have another face to transport the hydroxyl proton to the C_β atom to promote the tautomerization from enol intermediate to ketone species leading to the inhibition of the Fdc1 enzyme. The PrFMN roles are also discussed in detail.

A Projector-Embedding Approach for Multiscale Coupled-Cluster Calculations Applied to Citrate Synthase

Projector-based embedding has recently emerged as a robust multiscale method for the calculation of various electronic molecular properties. We present the coupling of projector embedding with quantum mechanics/molecular mechanics modeling and apply it for the first time to an enzyme-catalyzed reaction. Using projector-based embedding, we combine coupled-cluster theory, density-functional theory (DFT), and molecular mechanics to compute energies for the proton abstraction from acetyl-coenzyme A by citrate synthase. By embedding correlated ab initio methods in DFT we eliminate functional sensitivity and obtain high-accuracy profiles in a procedure that is straightforward to apply.

μ3-Oxo stabilized by three metal cations is a sufficient nucleophile for enzymatic hydrolysis of phosphate monoesters

Diverse species have previously been proposed to be effective nucleophiles in the enzymatic hydrolysis of phosphate esters. A novel penta-metal cluster (two Fe(3+) and three Ca(2+)) was recently discovered in the active site of PhoX alkaline phosphatase, with the revelation of the architecture of μ3-oxo bridging one Ca(2+) and two antiferromagnetically coupled Fe(3+). In this work, using density functional theory calculations, the μ3-oxo stabilized by three cations has been demonstrated to be a new type of effective nucleophile. The calculations give strong support to the "ping-pong" mechanism involving the nucleophilic attack of the μ3-oxo on the substrate phosphor and the subsequent hydrolysis of the covalent phospho-enzyme intermediate. A base mechanism with the μ3-oxo acting as a general base to activate an additional water molecule has further been demonstrated to be inaccessible. The results advance the understanding of the enzymatic hydrolysis of phosphate esters and may give important inspiration for the exploration of multinuclear biomimetic catalysts.