The fragmentation-recombination mechanism of the enzyme glutamate mutase studied by QM/MM simulations

Hydrolysis Mechanism of Carbamate Methomyl by a Novel Esterase PestE: A QM/MM Approach

Methomyl is one of the most important carbamates that has caused potential hazardous effects on both human beings and the environment. Here, we systematically investigated the hydrolysis mechanism of methomyl catalyzed by esterase PestE using molecular dynamics simulations (MD) and quantum mechanics/molecular mechanics (QM/MM) calculations. The hydrolysis mechanism involves two elementary steps: (Ⅰ) serine-initiated nucleophilic attack and (Ⅱ) C-O bond cleavage. Our work elicits the atomic level details of the hydrolysis mechanism and free energy profiles along the reaction pathway. The Boltzmann-weighted average potential barriers are 19.1 kcal/mol and 7.5 kcal/mol for steps Ⅰ and Ⅱ, respectively. We identified serine-initiated nucleophilic attack as the rate determining-step. The deep learning-based k_cat prediction model indicated that the barrier of the rate-determining step is 15.4 kcal/mol, which is in good agreement with the calculated results using Boltzmann-weighted average method. We have elucidated the importance of the protein-substrate interactions and the roles of the key active site residues during the hydrolysis process through noncovalent interactions analysis and electrostatic potential (ESP) analysis. The results provide practical value for achieving efficient degradation of carbamates by hydrolases.

C-N Bond Forming Radical Rebound Is the Enantioselectivity-Determining Step in P411-Catalyzed Enantioselective C(sp3)-H Amination: A Combined Computational and Experimental Investigation

Engineered metalloenzymes represent promising catalysts for stereoselective C-H functionalization reactions. Recently, P450 enzymes have been evolved to allow for new-to-nature intramolecular C(sp3)-H amination reactions via a nitrene transfer mechanism, giving rise to diamine derivatives with excellent enantiocontrol. To shed light on the origin of enantioselectivity, a combined computational and experimental study was carried out. Hybrid quantum mechanics/molecular mechanics calculations were performed to investigate the activation energies and enantioselectivities of both the hydrogen atom transfer (HAT) and the subsequent C-N bond forming radical rebound steps. Contrary to previously hypothesized enantioinduction mechanisms, our calculations show that the radical rebound step is enantioselectivity-determining, whereas the preceding HAT step is only moderately stereoselective. Furthermore, the selectivity in the initial HAT is ablated by rapid conformational change of the radical intermediate prior to C-N bond formation. This finding is corroborated by our experimental study using a set of enantiomerically pure, monodeuterated substrates. Furthermore, classical and ab initio molecular dynamics simulations were carried out to investigate the conformational flexibility of the carbon-centered radical intermediate. This key radical species undergoes a facile conformational change in the enzyme active site from the pro-(R) to the pro-(S) configuration, whereas the radical rebound is slower due to the spin-state change and ring strain of the cyclization process, thereby allowing stereoablative C-N bond formation. Together, these studies revealed an underappreciated enantioinduction mechanism in biocatalytic C(sp3)-H functionalizations involving radical intermediates, opening up new avenues for the development of other challenging asymmetric C(sp3)-H functionalizations.

Revealing the Mechanism of Isethionate Sulfite-Lyase by QM/MM Calculations

Isethionate sulfite-lyase (IseG) is a recently characterized glycyl radical enzyme (GRE) that catalyzes radical-mediated C-S bond cleavage of isethionate to produce acetaldehyde and sulfite. Herein, we use quantum mechanical/molecular mechanical (QM/MM) calculations to investigate the detailed catalytic reaction mechanism of IseG. Our calculations indicate that a previously proposed direct 1,2-elimination mechanism is disfavored. Instead, we suggest a new 1,2-migration mechanism for this enzymatic reaction: a key stepwise 1,2-SO₃^- radical migration occurs after the catalytically active cysteinyl radical grabs a hydrogen atom from isethionate, followed by hydrogen atom transfer from cysteine to a 1-hydroxylethane-1-sulfonate radical intermediate. Finally, the elimination of sulfite from 1-hydroxylethane-1-sulfonate to result in the final product is likely to occur outside the enzyme. Glu468 in the active site is found to help orient the substrate rather than grabbing a proton from the hydroxyl group of the substrate. Our findings help reveal the mechanisms of radical-mediated C-S bond cleavage of organosulfonates catalyzed by GREs and expand the understanding of radical-based enzymatic catalysis.

Computational study on the detoxifying mechanism of DDT metabolized by cytochrome P450 enzymes

Predicting the detoxifying mechanism and potential toxic derivatives of xenobiotic substances is significant for risk assessment. The present study delineated the detoxifying mechanism of 1-chloro-4-[2,2,2-trichloro-1-(4-chlorophenyl)ethyl]benzene (DDT) metabolized by human P450 enzymes using a combination of molecular dynamic (MD), quantum mechanics/molecular mechanics (QM/MM) and density functional theory (DFT). This study highlights that DDT can be metabolized by P450 enzymes through the hydrogen abstraction and electrophilic addition mechanism, and the main derivatives are epoxides (2,3-oxide-DDT and 3,4-oxide-DDT), DDE and dicofol. The epoxides are unstable and the C-O bond cleavage easily occurs by the reaction with hydronium ion or hydroxyl radicals, yielding endocrine disruptor hydroxylated DDT. The eco-toxicity evaluation indicates that the derivatives of DDT are less toxic than DDT, and the solubility increase of the derivatives can accelerate their excretion from the body. The study can provide an understanding of the biotransformation of DDT by the P450 enzymes in human livers.

A QM/MM study of the catalytic mechanism of SAM methyltransferase RlmN from Escherichia coli

RlmN is a radical S-adenosylmethionine (SAM) enzyme that catalyzes the C2 methylation of adenosine 2503 (A2503) in 23S rRNA and adenosine 37 (A37) in several Escherichia coli transfer RNAs (tRNA). The catalytic reaction of RlmN is distinctly different from that of typical SAM-dependent methyltransferases that employs an S_N 2 mechanism, but follows a ping-pong mechanism which involves the intermediate methylation of a conserved cysteine residue. Recently, the x-ray structure of a key intermediate in the RlmN reaction has been reported, allowing us to perform combined quantum mechanics and molecular mechanics (QM/MM) calculations to delineate the reaction details of RlmN at atomic level. Starting from the Cross-Linked RlmN C118A-tRNA complex, the possible mechanisms for both the formation and the resolution of the cross-linked species (IM2) have been illuminated. On the basis of our calculations, IM2 is formed by the attack of the C355-based methylene radical on the sp² -hybridized C2 of the adenosine ring, corresponding to energy barrier of 14.4 kcal/mol, and the resolution of IM2 is confirmed to follow a radical fragmentation mechanism. The cleavage of C'-S' bond of mC355-A37 cross-link is in concert with the deprotonation of C2 by C118 residue, which is the rate-limiting step with an energy barrier of 17.4 kcal/mol. Moreover, the cleavage of C'-S' bond of IM2 can occur independently, that is, it does not require the loss of an electron of IM2 and the formation of disulfide bond between C355 and C118 as precondition. These findings would deepen the understanding of the catalysis of RlmN.

Tryptophan lyase (NosL): mechanistic insights into amine dehydrogenation and carboxyl fragment migration by QM/MM calculations

QM/MM calculations suggest two feasible pathways for the breaking of the C–C bond of the substrate. The breaking of the Cα–Cβ bond leads to the final product, whereas the cleavage of the Cα–C bond will terminate in the EPR-trapped radical intermediate.

How Many Conformations of Enzymes Should Be Sampled for DFT/MM Calculations? A Case Study of Fluoroacetate Dehalogenase

The quantum mechanics/molecular mechanics (QM/MM) method (e.g., density functional theory (DFT)/MM) is important in elucidating enzymatic mechanisms. It is indispensable to study "multiple" conformations of enzymes to get unbiased energetic and structural results. One challenging problem, however, is to determine the minimum number of conformations for DFT/MM calculations. Here, we propose two convergence criteria, namely the Boltzmann-weighted average barrier and the disproportionate effect, to tentatively address this issue. The criteria were tested by defluorination reaction catalyzed by fluoroacetate dehalogenase. The results suggest that at least 20 conformations of enzymatic residues are required for convergence using DFT/MM calculations. We also tested the correlation of energy barriers between small QM regions and big QM regions. A roughly positive correlation was found. This kind of correlation has not been reported in the literature. The correlation inspires us to propose a protocol for more efficient sampling. This saves 50% of the computational cost in our current case.

Glutamate 338 is an electrostatic facilitator of C-Co bond breakage in a dynamic/electrostatic model of catalysis by ornithine aminomutase

How cobalamin-dependent enzymes promote C–Co homolysis to initiate radical catalysis has been debated extensively. For the pyridoxal 5′-phosphate and cobalamin-dependent enzymes lysine 5,6-aminomutase and ornithine 4,5-aminomutase (OAM), large-scale re-orientation of the cobalamin-binding domain linked to C–Co bond breakage has been proposed. In these models, substrate binding triggers dynamic sampling of the B₁₂-binding Rossmann domain to achieve a catalytically competent ‘closed’ conformational state. In ‘closed’ conformations of OAM, Glu338 is thought to facilitate C–Co bond breakage by close association with the cobalamin adenosyl group. We investigated this using stopped-flow continuous-wave photolysis, viscosity dependence kinetic measurements, and electron paramagnetic resonance spectroscopy of a series of Glu338 variants. We found that substrate-induced C–Co bond homolysis is compromised in Glu388 variant forms of OAM, although photolysis of the C–Co bond is not affected by the identity of residue 338. Electrostatic interactions of Glu338 with the 5′-deoxyadenosyl group of B₁₂ potentiate C–Co bond homolysis in ‘closed’ conformations only; these conformations are unlocked by substrate binding. Our studies extend earlier models that identified a requirement for large-scale motion of the cobalamin domain. Our findings indicate that large-scale motion is required to pre-organize the active site by enabling transient formation of ‘closed’ conformations of OAM. In ‘closed’ conformations, Glu338 interacts with the 5′-deoxyadenosyl group of cobalamin. This interaction is required to potentiate C–Co homolysis, and is a crucial component of the approximately 10¹² rate enhancement achieved by cobalamin-dependent enzymes for C–Co bond homolysis.

Exploring the substrate-assisted acetylation mechanism by UDP-linked sugar N-acetyltransferase from QM/MM calculations: the role of residue Asn84 and the effects of starting geometries

Based on the QM/MM calculation, we revised the proposed mechanism ofN-acetyltransferase and explore the role of Asn84 and the effects of starting geometries.

With or without light: comparing the reaction mechanism of dark-operative protochlorophyllide oxidoreductase with the energetic requirements of the light-dependent protochlorophyllide oxidoreductase

The addition of two electrons and two protons to the C₁₇=C₁₈ bond in protochlorophyllide is catalyzed by a light-dependent enzyme relying on NADPH as electron donor, and by a light-independent enzyme bearing a (Cys)₃Asp-ligated [4Fe–4S] cluster which is reduced by cytoplasmic electron donors in an ATP-dependent manner and then functions as electron donor to protochlorophyllide. The precise sequence of events occurring at the C₁₇=C₁₈ bond has not, however, been determined experimentally in the dark-operating enzyme. In this paper, we present the computational investigation of the reaction mechanism of this enzyme at the B3LYP/6-311+G(d,p)//B3LYP/6-31G(d) level of theory. The reaction mechanism begins with single-electron reduction of the substrate by the (Cys)₃Asp-ligated [4Fe–4S], yielding a negatively-charged intermediate. Depending on the rate of Fe–S cluster re-reduction, the reaction either proceeds through double protonation of the single-electron-reduced substrate, or by alternating proton/electron transfer. The computed reaction barriers suggest that Fe–S cluster re-reduction should be the rate-limiting stage of the process. Poisson–Boltzmann computations on the full enzyme–substrate complex, followed by Monte Carlo simulations of redox and protonation titrations revealed a hitherto unsuspected pH-dependence of the reaction potential of the Fe–S cluster. Furthermore, the computed distributions of protonation states of the His, Asp and Glu residues were used in conjuntion with single-point ONIOM computations to obtain, for the first time, the influence of all protonation states of an enzyme on the reaction it catalyzes. Despite exaggerating the ease of reduction of the substrate, these computations confirmed the broad features of the reaction mechanism obtained with the medium-sized models, and afforded valuable insights on the influence of the titratable amino acids on each reaction step. Additional comparisons of the energetic features of the reaction intermediates with those of common biochemical redox intermediates suggest a surprisingly simple explanation for the mechanistic differences between the dark-catalyzed and light-dependent enzyme reaction mechanisms.

Computational evidence for the detoxifying mechanism of epsilon class glutathione transferase toward the insecticide DDT

A combined quantum mechanics/molecular mechanics (QM/MM) computation of the detoxifying mechanism of an epsilon class glutathione transferases (GSTs) toward organochlorine insecticide DDT, 1,1,1-trichloro-2,2-bis(p-chlorophenyl)ethane, has been carried out. The exponential average barrier of the proton transfer mechanism is 15.2 kcal/mol, which is 27.6 kcal/mol lower than that of the GS-DDT conjugant mechanism. It suggests that the detoxifying reaction proceeds via a proton transfer mechanism where GSH acts as a cofactor rather than a conjugate. The study reveals that the protein environment has a strong effect on the reaction barrier. The experimentally proposed residues Arg112, Glu116 and Phe120 were found to have a strong influence on the detoxifying reaction. The influence of residues Pro13, Cys15, His53, Ile55, Glu67, Ser68, Phe115, and Leu119 was detected as well. It is worth noticing that Ile55 facilitates the detoxifying reaction most. On the basis of the structure of DDT, structure 2, (BrC6H4)2CHCCl3, is the best candidate among all the tested structures in resisting the detoxification of enzyme agGSTe2.

Reaction mechanism of the bicopper enzyme peptidylglycine α-hydroxylating monooxygenase

Peptidylglycine α-hydroxylating monooxygenase is a noninteracting bicopper enzyme that stereospecifically hydroxylates the terminal glycine of small peptides for its later amidation. Neuroendocrine messengers, such as oxytocin, rely on the biological activity of this enzyme. Each catalytic turnover requires one oxygen molecule, two protons from the solvent, and two electrons. Despite this enzyme having been widely studied, a consensus on the reaction mechanism has not yet been found. Experiments and theoretical studies favor a pro-S abstraction of a hydrogen atom followed by the rebinding of an OH group. However, several hydrogen-abstracting species have been postulated; because two protons are consumed during the reaction, several protonation states are available. An electron transfer between the copper atoms could play a crucial role for the catalysis as well. This leads to six possible abstracting species. In this study, we compare them on equal footing. We perform quantum mechanics/molecular mechanics calculations, considering the glycine hydrogen abstraction. Our results suggest that the most likely mechanism is a protonation of the abstracting species before the hydrogen abstraction and another protonation as well as a reduction before OH rebinding.

Role of active site residues in promoting cobalt-carbon bond homolysis in adenosylcobalamin-dependent mutases revealed through experiment and computation

Adenosylcobalamin (AdoCbl) serves as a source of reactive free radicals that are generated by homolytic scission of the coenzyme's cobalt-carbon bond. AdoCbl-dependent enzymes accelerate AdoCbl homolysis by ∼10(12)-fold, but the mechanism by which this is accomplished remains unclear. We have combined experimental and computational approaches to gain molecular-level insight into this process for glutamate mutase. Two residues, glutamate 330 and lysine 326, form hydrogen bonds with the adenosyl group of the coenzyme. A series of mutations that impair the enzyme's ability to catalyze coenzyme homolysis and tritium exchange with the substrate by 2-4 orders of magnitude were introduced at these positions. These mutations, together with the wild-type enzyme, were also characterized in silico by molecular dynamics simulations of the enzyme-AdoCbl-substrate complex with AdoCbl modeled in the associated (Co-C bond formed) or dissociated [adenosyl radical with cob(II)alamin] state. The simulations reveal that the number of hydrogen bonds between the adenosyl group and the protein side chains increases in the homolytically dissociated state, with respect to the associated state, for both the wild-type and mutant enzymes. The mutations also cause a progressive increase in the mean distance between the 5'-carbon of the adenosyl radical and the abstractable hydrogen of the substrate. Interestingly, the distance between the 5'-carbon and substrate hydrogen, determined computationally, was found to inversely correlate with the log k for tritium exchange (r = 0.93) determined experimentally. Taken together, these results point to a dual role for these residues: they both stabilize the homolytic state through electrostatic interactions between the protein and the dissociated coenzyme and correctly position the adenosyl radical to facilitate the abstraction of hydrogen from the substrate.

Radical use of Rossmann and TIM barrel architectures for controlling coenzyme B12 chemistry

The ability of enzymes to harness free-radical chemistry allows for some of the most amazing transformations in nature, including reduction of ribonucleotides and carbon skeleton rearrangements. Enzyme cofactors involved in this chemistry can be large and complex, such as adenosylcobalamin (coenzyme B(12)), simpler, such as S-adenosylmethionine and an iron-sulfur cluster (i.e., poor man's B(12)), or very small, such as one nonheme iron atom coordinated by protein ligands. Although the chemistry catalyzed by these enzyme-bound cofactors is unparalleled, it does come at a price. The enzyme must be able to control these radical reactions, preventing unwanted chemistry and protecting the enzyme active site from damage. Here, we consider a set of radical folds: the (β/α)(8) or TIM barrel, combined with a Rossmann domain for coenzyme B(12)-dependent chemistry. Using specific enzyme examples, we consider how nature employs the common TIM barrel fold and its Rossmann domain partner for radical-based chemistry.

Role of tunneling in the enzyme glutamate mutase

The role of quantum mechanical atom tunneling during the conversion of glutamate to methylaspartate catalyzed by glutamate mutase is investigated by quantum mechanical/molecular mechanical (QM/MM) simulations based on coupled cluster and density functional calculations. The use of instanton theory allows us to calculate the tunneling contributions of up to 78 atoms in the active site. We calculate kinetic isotope effects (KIEs) and compare them to experimental data. The simulations lead to deuterium KIEs of 10 for the hydrogen abstraction from glutamate substrate and 16 for the hydrogen abstraction from methylaspartate substrate, which are consistent with the experimental results. The hydrogen abstraction from methylaspartate has higher primary deuterium and tritium (46.1) KIEs than the abstraction from glutamate. The tunneling effect increases the reaction rate by a factor of 12.3 for the hydrogen abstraction from methylaspartate at 0. Tunneling is supported by the environment by preparing the enzyme through classical motions. Consideraton of the tunneling contributions of more and more atoms around the active center shows that the motions at the ribose ring play a central role during the tunneling enhancement of the hydrogen transfers. Our simulations give new insight into the catalytic process in glutamate mutase and the way enzymes use tunneling effects for a successful catalysis.

Experimental study of hydrogen bonding potentially stabilizing the 5'-deoxyadenosyl radical from coenzyme B12

Coenzyme B(12) can assist radical enzymes that accomplish the vicinal interchange of a hydrogen atom with a functional group. It has been proposed that the Co-C bond homolysis of coenzyme B(12) to cob(II)alamin and the 5'-deoxyadenosyl radical is aided by hydrogen bonding of the corrin C19-H to the 3'-O of the ribose moiety of the incipient 5'-deoxyadenosyl radical, which is stabilized by 30 kJ mol(-1) (B. Durbeej et al., Chem. Eur. J. 2009, 15, 8578-8585). The diastereoisomers (R)- and (S)-2,3-dihydroxypropylcobalamin were used as models for coenzyme B(12). A downfield shift of the NMR signal for the C19-H proton was observed for the (R)-isomer (δ=4.45 versus 4.01 ppm for the (S)-isomer) and can be ascribed to an intramolecular hydrogen bond between the C19-H and the oxygen of CHOH. Crystal structures of (R)- and (S)-2,3-dihydroxypropylcobalamin showed C19-H⋅⋅⋅O distances of 3.214(7) Å (R-isomer) and 3.281(11) Å (S-isomer), which suggest weak hydrogen-bond interactions (-ΔG<6 kJ mol(-1)) between the CHOH of the dihydroxypropyl ligand and the C19-H. Exchange of the C19-H, which is dependent on the cobalt redox state, was investigated with cob(I)alamin, cob(II)alamin, and cob(III)alamin by using NMR spectroscopy to monitor the uptake of deuterium from deuterated water in the pH range 3-11. No exchange was found for any of the cobalt oxidation states. 3',5'-Dideoxyadenosylcobalamin, but not the 2',5'-isomer, was found to act as a coenzyme for glutamate mutase, with a 15-fold lower k(cat)/K(M) than 5'-deoxyadenosylcobalamin. This indicates that stabilization of the 5'-deoxyadenosyl radical by a hydrogen bond that involves the C19-H and the 3'-OH group of the cofactor is, at most, 7 kJ mol(-1) (-ΔG). Examination of the crystal structure of glutamate mutase revealed additional stabilizing factors: hydrogen bonds between both the 2'-OH and 3'-OH groups and glutamate 330. The actual strength of a hydrogen bond between the C19-H and the 3'-O of the ribose moiety of the 5'-deoxyadenosyl group is concluded not to exceed 6 kJ mol(-1) (-ΔG).

Adenosylcobalamin enzymes: theory and experiment begin to converge

Adenosylcobalamin (coenzyme B(12)) serves as the cofactor for a group of enzymes that catalyze unusual rearrangement or elimination reactions. The role of the cofactor as the initiator of reactive free radicals needed for these reactions is well established. Less clear is how these enzymes activate the coenzyme towards homolysis and control the radicals once generated. The availability of high resolution X-ray structures combined with detailed kinetic and spectroscopic analyses have allowed several adenosylcobalamin enzymes to be computationally modeled in some detail. Computer simulations have generally obtained good agreement with experimental data and provided valuable insight into the mechanisms of these unusual reactions. Importantly, atomistic modeling of the enzymes has allowed the role of specific interactions between protein, substrate and coenzyme to be explored, leading to mechanistic predictions that can now be tested experimentally. This article is part of a Special Issue entitled: Radical SAM enzymes and Radical Enzymology.

Combined spectroscopic/computational studies of vitamin B12 precursors: geometric and electronic structures of cobinamides

Vitamin B(12) (cyanocobalamin) and its biologically active derivatives, methylcobalamin and adenosylcobalamin, are members of the family of corrinoids, which also includes cobinamides. As biological precursors to cobalamins, cobinamides possess the same structural core, consisting of a low-spin Co(3+) ion that is ligated equatorially by the four nitrogens of a highly substituted tetrapyrrole macrocycle (the corrin ring), but differ with respect to the lower axial ligation. Specifically, cobinamides possess a water molecule instead of the nucleotide loop that coordinates axially to Co(3+)cobalamins via its dimethylbenzimidazole (DMB) base. Compared to the cobalamin species, cobinamides have proven much more difficult to study experimentally, thus far eluding characterization by X-ray crystallography. In this study, we have utilized combined quantum mechanics/molecular mechanics (QM/MM) computations to generate complete structural models of a representative set of cobinamide species with varying upper axial ligands. To validate the use of this approach, analogous QM/MM geometry optimizations were carried out on entire models of the cobalamin counterparts for which high-resolution X-ray structural data are available. The accuracy of the cobinamide structures was assessed further by comparing electronic absorption spectra computed using time-dependent density functional theory to those obtained experimentally. Collectively, the results obtained in this study indicate that the DMB → H(2)O lower axial ligand switch primarily affects the energies of the Co 3d(z(2))-based molecular orbital (MO) and, to a lesser extent, the other Co 3d-based MOs as well as the corrin π-based highest energy MO. Thus, while the energy of the lowest-energy electronic transition of cobalamins changes considerably as a function of the upper axial ligand, it is nearly invariant for the cobinamides.