Ab Initio QM/MM Modeling of the Rate-Limiting Proton Transfer Step in the Deamination of Tryptamine by Aromatic Amine Dehydrogenase

On the Mechanisms of Hypohalous Acid Formation and Electrophilic Halogenation by Non-Native Halogenases

Enzymatic electrophilic halogenation is a mild tool for functionalization of diverse organic compounds. Only a few groups of native halogenases are capable of catalyzing such a reaction. In this study, we used a mechanism-guided strategy to discover the electrophilic halogenation activity catalyzed by non-native halogenases. As the ability to form a hypohalous acid (HOX) is key for halogenation, flavin-dependent monooxygenases/oxidases capable of forming C4a-hydroperoxyflavin (FlC4a-OOH), such as dehalogenase, hydroxylases, luciferase and pyranose-2-oxidase (P2O), and flavin reductase capable of forming H2O2 were explored for their abilities to generate HOX in situ. Transient kinetic analyses using stopped-flow spectrophotometry/fluorometry and product analysis indicate that FlC4a-OOH in dehalogenases, selected hydroxylases and luciferases, but not in P2O can form HOX; however, the HOX generated from FlC4a-OOH cannot halogenate their substrates. Remarkably, in situ H2O2 generated by P2O can form HOI and also iodinate various compounds. Because not all enzymes capable of forming FlC4a-OOH can react with halides to form HOX, QM/MM calculations, site-directed mutagenesis and structural analysis were carried out to elucidate the mechanism underlying HOX formation and characterize the active site environment. Our findings shed light on identifying new halogenase scaffolds besides the currently known enzymes and have invoked a new mode of chemoenzymatic halogenation.

Screening and evaluating microorganisms with broad-spectrum biogenic amine-degrading ability from naturally fermented dry sausage collected from Northeast China

This study aimed to screen autochthonous strains with broad-spectrum biogenic amine (BA) degradation ability from traditional dry sausages and to evaluate their BA-degrading ability in dry sausages. A total of 120 strains were isolated from dry sausages collected from various regions in Northeast China, and 35 of 120 isolates were identified as non-BA producing strains by the in vitro agar method. The random amplified polymorphic DNA polymerase chain reaction technique genotyped these 35 isolates into 18 biotypes. Moreover, high performance liquid chromatography (HPLC) quantification showed that six strains (Latilactobacillus sakei MDJ6; Lactiplantibacillus plantarum SH7; Weissella hellenica DQ9; Staphylococcus saprophyticus JX18 and SYS8; and Macrococcus caseolyticus SYS11) of the 18 biotypes exhibited broad-spectrum BA-degrading ability, all of which had various levels of amine oxidase activity with monoamine oxidase and diamine oxidase activities ranged of 6.60-619.04 and 26.32-352.81 U/mg protein, respectively. These six strains were subsequently inoculated into dry sausages and the results showed that they exhibited varying degrees of BA-degrading ability, of which strain Lat. sakei MDJ6 allowed to have less BA production on dry sausage with a final concentration of 61.33 mg/kg.

MoBioTools: A toolkit to setup quantum mechanics/molecular mechanics calculations

We present a toolkit that allows for the preparation of QM/MM input files from a conformational ensemble of molecular geometries. The package is currently compatible with trajectory and topology files in Amber, CHARMM, GROMACS and NAMD formats, and has the possibility to generate QM/MM input files for Gaussian (09 and 16), Orca (≥4.0), NWChem and (Open)Molcas. The toolkit can be used in command line, so that no programming experience is required, although it presents some features that can also be employed as a python application programming interface. We apply the toolkit in four situations in which different electronic-structure properties of organic molecules in the presence of a solvent or a complex biological environment are computed: the reduction potential of the nucleobases in acetonitrile, an energy decomposition analysis of tyrosine interacting with water, the absorption spectrum of an azobenzene derivative integrated into a voltage-gated ion channel, and the absorption and emission spectra of the luciferine/luciferase complex. These examples show that the toolkit can be employed in a manifold of situations for both the electronic ground state and electronically excited states. It also allows for the automatic correction of the active space in the case of CASSCF calculations on an ensemble of geometries, as it is shown for the azobenzene derivative photoswitch case.

A QM/MM Evaluation of the Missing Step in the Reduction Mechanism of HMG-CoA by Human HMG-CoA Reductase

Statins are important drugs in the regulation of cholesterol levels in the human body that have as a primary target the enzyme β-hydroxy-β-methylglutaryl-CoA reductase (HMGR). This enzyme plays a crucial role in the mevalonate pathway, catalyzing the four-electron reduction of HMG-CoA to mevalonate. A second reduction step of this reaction mechanism has been the subject of much speculation in the literature, with different conflicting theories persisting to the present day. In this study, the different mechanistic hypotheses were evaluated with atomic-level detail through a combination of molecular dynamics simulations (MD) and quantum mechanics/molecular mechanics (QM/MM) calculations. The obtained Gibbs free activation and Gibbs free reaction energy (15 kcal mol−1 and −40 kcal mol−1) show that this hydride step takes place with the involvement of a cationic His405 and Lys639, and a neutral Glu98, while Asp715 remains in an anionic state. The results provide an atomic-level portrait of this step, clearly demonstrating the nature and protonation state of the amino acid residues involved, the energetics associated, and the structure and charge of the key participating atoms in the several intermediate states, finally elucidating this missing step.

Dissecting the low catalytic capability of flavin-dependent halogenases

Although flavin-dependent halogenases (FDHs) are attractive biocatalysts, their practical applications are limited because of their low catalytic efficiency. Here, we investigated the reaction mechanisms and structures of tryptophan 6-halogenase (Thal) from Streptomyces albogriseolus using stopped-flow, rapid-quench flow, quantum/mechanics molecular mechanics calculations, crystallography, and detection of intermediate (hypohalous acid [HOX]) liberation. We found that the key flavin intermediate, C4a-hydroperoxyflavin (C4aOOH-FAD), formed by Thal and other FDHs (tryptophan 7-halogenase [PrnA] and tryptophan 5-halogenase [PyrH]), can react with I-, Br-, and Cl- but not F- to form C4a-hydroxyflavin and HOX. Our experiments revealed that I- reacts with C4aOOH-FAD the fastest with the lowest energy barrier and have shown for the first time that a significant amount of the HOX formed leaks out as free HOX. This leakage is probably a major cause of low product coupling ratios in all FDHs. Site-saturation mutagenesis of Lys79 showed that changing Lys79 to any other amino acid resulted in an inactive enzyme. However, the levels of liberated HOX of these variants are all similar, implying that Lys79 probably does not form a chloramine or bromamine intermediate as previously proposed. Computational calculations revealed that Lys79 has an abnormally lower pKa compared with other Lys residues, implying that the catalytic Lys may act as a proton donor in catalysis. Analysis of new X-ray structures of Thal also explains why premixing of FDHs with reduced flavin adenine dinucleotide generally results in abolishment of C4aOOH-FAD formation. These findings reveal the hidden factors restricting FDHs capability which should be useful for future development of FDHs applications.

Exploring Free Energy Profiles of Enantioselective Organocatalytic Aldol Reactions under Full Solvent Influence

We present a computational study on the enantioselectivity of organocatalytic proline-catalyzed aldol reactions between aldehydes in dimethylformamide (DMF). To explore the free energy surface of the reaction, we apply two-dimensional metadynamics on top of ab initio molecular dynamics (AIMD) simulations with explicit solvent description on the DFT level of theory. We avoid unwanted side reactions by utilizing our newly developed hybrid AIMD (HyAIMD) simulation scheme, which adds a simple force field to the AIMD simulation to prevent unwanted bond breaking and formation. Our condensed phase simulation results are able to nicely reproduce the experimental findings, including the main stereoisomer that is formed, and give a correct qualitative prediction of the change in syn:anti product ratio with different substituents. Furthermore, we give a microscopic explanation for the selectivity. We show that both the explicit description of the solvent and the inclusion of entropic effects are vital to a good outcome-metadynamics simulations in vacuum and static nudged elastic band (NEB) calculations yield significantly worse predictions when compared to the experiment. The approach described here can be applied to a plethora of other enantioselective or organocatalytic reactions, enabling us to tune the catalyst or determine the solvent with the highest stereoselectivity.

Electronic structure benchmark calculations of CO2 fixing elementary chemical steps in RuBisCO using the projector-based embedding approach

Ribulose 1,5-bisphosphate carboxylase-oxygenase (RuBisCO) is the main enzyme involved in atmospheric carbon dioxide (CO₂ ) fixation in the biosphere. This enzyme catalyzes a set of five chemical steps that take place in the same active-site within magnesium (II) coordination sphere. Here, a set of electronic structure benchmark calculations have been carried out on a reaction path proposed by Gready et al. by means of the projector-based embedding approach. Activation and reaction energies for all main steps catalyzed by RuBisCO have been calculated at the MP2, SCS-MP2, CCSD, and CCSD(T)/aug-cc-pVDZ and cc-pVDZ levels of theory. The treatment of the magnesium cation with post-HF methods is explored to determine the nature of its involvement in the mechanism. With the high-level ab initio values as a reference, we tested the performance of a set of density functional theory (DFT) exchange-correlation (xc) functionals in reproducing the reaction energetics of RuBisCO carboxylase activity on a set of model fragments. Different DFT xc-functionals show large variation in activation and reaction energies. Activation and reaction energies computed at the B3LYP level are close to the reference SCS-MP2 results for carboxylation, hydration and protonation reactions. However, for the carbon-carbon bond dissociation reaction, B3LYP and other functionals give results that differ significantly from the ab initio reference values. The results show the applicability of the projector-based embedding approach to metalloenzymes. This technique removes the uncertainty associated with the selection of different DFT xc-functionals and so can overcome some of inherent limitations of DFT calculations, complementing, and potentially adding to modeling of enzyme reaction mechanisms with DFT methods.

Embedded Mean-Field Theory for Solution-Phase Transition-Metal Polyolefin Catalysis

Decreasing the wall-clock time of quantum mechanics/molecular mechanics (QM/MM) calculations without sacrificing accuracy is a crucial prerequisite for widespread simulation of solution-phase dynamical processes. In this work, we demonstrate the use of embedded mean-field theory (EMFT) as the QM engine in QM/MM molecular dynamics (MD) simulations to examine polyolefin catalysts in solution. We show that employing EMFT in this mode preserves the accuracy of hybrid-functional DFT in the QM region, while providing up to 20-fold reductions in the cost per SCF cycle, thereby increasing the accessible simulation time-scales. We find that EMFT reproduces DFT-computed binding energies and optimized bond lengths to within chemical accuracy, as well as consistently ranking conformer stability. Furthermore, solution-phase EMFT/MM simulations provide insight into the interaction strength of strongly coordinating and bulky counterions.

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.

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.

Quantum Mechanics/Molecular Mechanics (QM/MM) Calculations Support a Concerted Reaction Mechanism for the Zika Virus NS2B/NS3 Serine Protease with Its Substrate

Zika virus (ZIKV) is mainly transmitted to humans by Aedes species mosquitoes and is associated with serious pathological disorders including microcephaly in newborns and Guillain-Barré syndrome in adults. Currently, there is no vaccine or anti-ZIKV drug available for preventing or controlling ZIKV infection. An attractive drug target for ZIKV treatment is a two-compartment (NS2B/NS3) serine protease that processes viral polyprotein during infection. Here, conventional molecular dynamics simulations of the ZIKV protease in complex with peptide substrate (TGKRS) sequence at the C-terminus of NS2B show that the substrate is in the active conformation for the cleavage reaction by ZIKV protease. Hybrid quantum mechanics/molecular mechanics (QM/MM) umbrella sampling simulations (PM6/ff14SB) of acylation results reveal that proton transfer from S135 to H51 and nucleophilic attack on the substrate by S135 are concerted. The rate-limiting step involves the formation of a tetrahedral intermediate. In addition, the single-point energy QM/MM calculations, precisely at the level of coupled cluster theory (LCCSD(T)/(aug)-cc-pVTZ), were performed to correct the potential energy profiles for the first step of the acylation process. The average computed activation barrier at this level of theory is 16.3 kcal mol^-1. Therefore, the computational approaches presented here are helpful for further designing of NS2B/NS3 inhibitors based on transition-state analogues.

QM/MM molecular modelling on mutation effect of chorismate synthase enzyme catalysis

Chorismate synthase (CS) catalyzes the conversion of 5-enolpyruvylshikimate-3-phosphate (EPSP) to chorismate which is a key intermediate in the biosynthesis of aromatic amino acids. CS enzyme is a new target for antibacterial drugs. Even though several reaction mechanisms have been proposed, the catalytic mechanism is still unclear. QM/MM adiabatic mapping calculations were performed in order to investigate roles of this enzyme. High-accuracy SCS-MP2/aVDZ/CHARMM27 calculations indicated that the reaction pathway has three steps; (i) proton transfer from reduced flavin mononucleotide (FMNH₂) to D339, (ii) proton transfer from EPSP to FMNH^- and (iii) phosphate elimination. Adiabatic mapping calculations indicated that H110 and R48 residues play essential catalyst roles for CS enzyme catalysis by transition state (TS) and product stabilizations via charge polarization and hydrogen bonding to EPSP and/or FMNH₂. A high accuracy calculation - SCS-MP2/aVDZ/CHARMM27 method was employed to obtain the accurate reaction mechanism pathway and to evaluate the effect of amino acid residues in the active site on the enzyme catalysis. The potential energy barriers of the reactions of H110A and R48A were found to increase. The CS catalysis was consequently slowed down due to missing the TS and product stabilizations.