How the Conformational Movement of the Substrate Drives the Regioselective C-N Bond Formation in P450 TleB: Insights from Molecular Dynamics Simulations and Quantum Mechanical/Molecular Mechanical Calculations.

P450 TleB catalyzes the oxidative cyclization of the dipeptide N-methylvalyl-tryptophanol into indolactam V through selective intramolecular C-H bond amination at the indole C4 position. Understanding its catalytic mechanism is instrumental for the engineering or design of P450-catalyzed C-H amination reactions. Using multiscale computational methods, we show that the reaction proceeds through a diradical pathway, involving a hydrogen atom transfer (HAT) from N1-H to Cpd I, a conformational transformation of the substrate radical species, and a second HAT from N13-H to Cpd II. Intriguingly, the conformational transformation is found to be the key to enabling efficient and selective C-N coupling between N13 and C4 in the subsequent diradical coupling reaction. The underlined conformational transformation is triggered by the first HAT, which proceeds with an energy-demanding indole ring flip and is followed by the facile approach of the N13-H group to Cpd II. Detailed analysis shows that the internal electric field (IEF) from the protein environment plays key roles in the transformation process, which not only provides the driving force but also stabilizes the flipped conformation of the indole radical. Our simulations provide a clear picture of how the P450 enzyme can smartly modulate the selective C-N coupling reaction. The present findings are in line with all available experimental data, highlighting the crucial role of substrate dynamics in controlling this highly valuable reaction.

Biodegradation of 2,5-Dihydroxypyridine by 2,5-Dihydroxypyridine Dioxygenase and Its Mutants: Insights into O-O Bond Activation and Flexible Reaction Mechanisms from QM/MM Simulations.

2,5-Dihydroxypyridine dioxygenase (NicX) from Pseudomonas putida KT2440 is a mononuclear non-heme iron oxygenase responsible for the biodegradation of 2,5-dihydroxypyridine (DHP) to N-formylmaleamic acid (NFM). Here, extensive quantum mechanical-molecular mechanical (QM/MM) calculations and molecular dynamics (MD) simulations are used to elucidate the degradation mechanism of DHP by wild-type NicX and its H105F variant (NicXH105F) and the roles of key residues. In particular, NicX and NicXH105F can catalyze the ring opening degradation of DHP to NFM, but flexible mechanisms are adopted therein. Both reactions of NicX and NicXH105F are initiated by the attack of FeIII superoxide species onto the substrate, during which a proton-coupled electron transfer (PCET) process is involved. For wild-type NicX, the PCET reaction is mediated by the adjacent His105, while the further proton transfer from His105 to the peroxo species can remarkably enhance the following O-O cleavage. However, for the NicXH105F mutant, a water molecule replaces the role of residue His105, which not only stabilizes the substrate binding via a H bonding network but also functions as a base to mediate the PCET process. For the NicXH105A mutant, MD simulations show that the disruption of the H bonding network can displace the substrate binding, leading to the loss of enzyme activity. These findings can expand our understanding of the PCET-mediated O-O bond activation and the flexible catalytic routes in various mutants, which have general implications on enzyme catalysis.

Degradation of pesticides diazinon and diazoxon by phosphotriesterase: insight into divergent mechanisms from QM/MM and MD simulations.

Enzymatic hydrolysis by phosphotriesterase (PTE) is one of the most effective ways of degrading organophosphorus pesticides, but the catalytic efficiency depends on the structural features of substrates. Here the enzymatic degradation of diazinon (DIN) and diazoxon (DON), characterized by PS and PO, respectively, have been investigated by QM/MM calculations and MM MD simulations. Our calculations demonstrate that the hydrolysis of DON (with PO) is inevitably initiated by the nucleophilic attack of the bridging-OH- on the phosphorus center, while for DIN (with PS), we proposed a new degradation mechanism, initiated by the nucleophilic attack of the Znα-bound water molecule, for its low-energy pathway. For both DIN and DON, the hydrolytic reaction is predicted to be the rate-limiting step, with energy barriers of 18.5 and 17.7 kcal mol-1, respectively. The transportation of substrates to the active site, the release of the leaving group and the degraded product are generally verified to be favorable by MD simulations via umbrella sampling, both thermodynamically and dynamically. The side-chain residues Phe132, Leu271 and Tyr309 play the gate-switching role to manipulate substrate delivery and product release. In comparison with the DON-enzyme system, the degraded product of DIN is more easily released from the active site. These new findings will contribute to the comprehensive understanding of the enzymatic degradation of toxic organophosphorus compounds by PTE.

QM/MM and MM MD simulations on decontamination of the V-type nerve agent VX by phosphotriesterase: toward a comprehensive understanding of steroselectivity and activity.

Due to deadly toxicity and high environmental stability of the nerve agent VX, an efficient decontamination approach is desperately needed in tackling its severe threat to human security. The enzymatic destruction of nerve agents has been generally considered as one of the most effective ways, and here the hydrolysis of VX by phosphotriesterase (PTE) was investigated by extensive QM/MM and MM MD simulations. The hydrolytic cleavage of P-S by PTE is a two-step process with the free energy spans of 15.8 and 26.0 kcal mol-1 for the RP- and SP-enantiomer VX, respectively, and such remarkable stereospecificity of VX enantiomers in the enzymatic degradation is attributed to their conformational compatibility with the active pocket. The structurally less adaptive SP-enantiomer allows one additional water molecule to enter the binuclear zinc center and remarkably facilitates the release of the degraded product. Overall, the rate-limiting steps in the enzymatic degradation of VX by PTE involve the degraded product release of the RP-enantiomer and the enzymatic P-S cleavage of the SP-enantiomer. Further computational analysis on the mutation of selected residues also revealed that H257Y, H257D, H254Q-H257F, and L7ep-3a variants allow more water molecules to enter the active site, which improves the catalytic efficiency of PTE, as observed experimentally. The present work provides mechanistic insights into the stereoselective hydrolysis of VX by PTE and the activity manipulation through the active-site accessibility of water molecules, which can be used for the enzyme engineering to defeat chemical warfare agents.

Hydrogen- and Halogen-Bonds between Ions of like Charges: Are They Anti-Electrostatic in Nature?

Recent theoretical studies suggested that hydrogen bonds between ions of like charges are of a covalent nature due to the dominating nD →σ*H-A charge-transfer (CT) interaction. In this work, energy profiles of typical hydrogen (H) and halogen (X) bonding systems formed from ions of like charges are explored using the block-localized wavefunction (BLW) method, which can derive optimal geometries and wave functions with the CT interaction "turned off." The results demonstrate that the kinetic stability, albeit reduced, is maintained for most investigated systems even after the intermolecular CT interaction is quenched. Further energy decomposition analyses based on the BLW method reveal that, despite a net repulsive Coulomb repulsion, a stabilizing component exists due to the polarization effect that plays significant role in the kinetic stability of all systems. Moreover, the fingerprints of the augmented electrostatic interaction due to polarization are apparent in the variation patterns of the electron density. All in all, much like in standard H- and X-bonds, the stability of such bonds between ions of like charges is governed by the competition between the stabilizing electrostatic and charge transfer interactions and the destabilizing deformation energy and Pauli exchange repulsion. While in most cases of "anti-electrostatic" bonds the CT interaction is of a secondary importance, we also find cases where CT is decisive. As such, this work validates the existence of anti-electrostatic H- and X-bonds. © 2017 Wiley Periodicals, Inc.

Elucidating the Enzymatic Mechanism of Dihydrocoumarin Degradation: Insight into the Functional Evolution of Methyl-Parathion Hydrolase from QM/MM and MM MD Simulations.

Methyl-parathion hydrolase (MPH), which evolved from dihydrocoumarin hydrolase, offers one of the most efficient enzymes for the hydrolysis of methyl-parathion. Interestingly, the substrate preference of MPH shifts from the methyl-parathion to the lactone dihydrocoumarin (DHC) after its mutation of five specific residues (R72L, L273F, L258H, T271I, and S193Δ, m5-MPH). Here, extensive QM/MM calculations and MM MD simulations have been used to delve into the structure-function relationship of MPH enzymes and plausible mechanisms for the chemical and nonchemical steps, including the transportation and binding of the substrate DHC to the active site, the hydrolysis reaction, and the product release. The results reveal that the five mutations remodel the active pocket and reposition DHC within the active site, leading to stronger enzyme-substrate interactions. The MM/GBSA-estimated binding free energies are about -20.7 kcal/mol for m5-MPH and -17.1 kcal/mol for wild-type MPH. Furthermore, this conformational adjustment of the protein may facilitate the chemical step of DHC hydrolysis and the product release, although there is a certain influence on the substrate transport. The hydrolytic reaction begins with the nucleophilic attack of the bridging OH- with the energy barriers of 22.0 and 18.0 kcal/mol for the wild-type and m5-MPH enzymes, respectively, which is rate-determining for the entire process. Unraveling these mechanistic intricacies may help in the understanding of the natural evolution of enzymes for diverse substrates and establish the enzyme structure-function relationship.

QM/MM and MM MD Simulations on Enzymatic Degradation of the Nerve Agent VR by Phosphotriesterase.

V-type nerve agents are hardly degraded by phosphotriesterase (PTE). Interestingly, the PTE variant of BHR-73MNW can effectively improve the hydrolytic efficiency of VR, especially for its Sp-enantiomer. Here, the whole enzymatic degradation of both Sp and Rp enantiomers of VR by the wild-type PTE and its variant BHR-73MNW was investigated by quantum mechanics/molecular mechanics (QM/MM) calculations and MM molecular dynamics simulations. Present results indicate that the degradation of VR can be initiated by the nucleophilic attack of the bridging OH- and the zinc-bound water molecule. The QM/MM-predicted energy barriers for the hydrolytic process of Sp-VR are 19.8 kcal mol-1 by the variant with water as a nucleophile and 22.0 kcal mol-1 by the wild-type PTE with OH- as a nucleophile, and corresponding degraded products are bound to the dinuclear metal site in monodentate and bidentate coordination modes, respectively. The variant effectively increases the volume of the large pocket, allowing more water molecules to enter the active pocket and resulting in the improvement of the degradation efficiency of Sp-VR. The hydrolysis of Rp-VR is triggered only by the hydroxide with an energy span of 20.6 kcal mol-1 for the wild-type PTE and 20.7 kcal mol-1 for the variant BHR-73-MNW PTE. Such mechanistic insights into the stereoselective degradation of VR by PTE and the role of water may inspire further studies to improve the catalytic efficiency of PTE toward the detoxification of nerve agents.

Water-Regulated Mechanisms for Degradation of Pesticides Paraoxon and Parathion by Phosphotriesterase: Insight from QM/MM and MD Simulations.

The enzymatic degradation of pesticides paraoxon (PON) and parathion (PIN) by phosphotriesterase (PTE) has been investigated by QM/MM calculations and MD simulations. In the PTE-PON complex, Znα and Znß in the active site are five- and six-coordinated, respectively, while both zinc ions are six coordinated with the Znα -bound water molecule (WT1) for the PTE-PIN system. The hydrolytic reactions for PON and PIN are respectively driven by the nucleophilic attack of the bridging-OH- and the Znα -bound water molecule on the phosphorus center of substrate, and the two-step hydrolytic process is predicted to be the rate-limiting step with the energy spans of 13.8 and 14.4 kcal/mol for PON and PIN, respectively. The computational studies reveal that the presence of the Znα -bound water molecule depends on the structural feature of substrates characterized by P=O and P=S, which determines the hydrolytic mechanism and efficiency for the degradation of organophosphorus pesticides by PTE.

Conformational Change of H64 and Substrate Transportation: Insight Into a Full Picture of Enzymatic Hydration of CO2 by Carbonic Anhydrase.

The enzymatic hydration of CO2 into HCO3 - by carbonic anhydrase (CA) is highly efficient and environment-friendly measure for CO2 sequestration. Here extensive MM MD and QM/MM MD simulations were used to explore the whole enzymatic process, and a full picture of the enzymatic hydration of CO2 by CA was achieved. Prior to CO2 hydration, the proton transfer from the water molecule (WT1) to H64 is the rate-limiting step with the free energy barrier of 10.4 kcal/mol, which leads to the ready state with the Zn-bound OH-. The nucleophilic attack of OH- on CO2 produces HCO3 - with the free energy barrier of 4.4 kcal/mol and the free energy release of about 8.0 kcal/mol. Q92 as the key residue manipulates both CO2 transportation to the active site and release of HCO3 -. The unprotonated H64 in CA prefers in an inward orientation, while the outward conformation is favorable energetically for its protonated counterpart. The conformational transition of H64 between inward and outward correlates with its protonation state, which is mediated by the proton transfer and the product release. The whole enzymatic cycle has the free energy span of 10.4 kcal/mol for the initial proton transfer step and the free energy change of -6.5 kcal/mol. The mechanistic details provide a comprehensive understanding of the entire reversible conversion of CO2 into bicarbonate and roles of key residues in chemical and nonchemical steps for the enzymatic hydration of CO2.