Reaction mechanism of soluble epoxide hydrolase: insights from molecular dynamics simulations

Engineering an Epoxide Hydrolase for Chemoenzymatic Asymmetric Synthesis of Chiral Triazole Fungicide (S)- and (R)-Flutriafol

Flutriafol, a globally utilized triazole fungicide in agriculture, is typically applied as a racemic mixture, but its enantiomers differ in bioactivity and environmental impact. The synthesis of flutriafol enantiomers is critically dependent on chiral precursors: 2,2-bisaryl-substituted oxirane [(2-fluorophenyl)-2-(4-fluorophenyl)oxirane, 1a] and 1,2-diol [1-(2-fluorophenyl)-1-(4-fluorophenyl)ethane-1,2-diol, 1b]. Here, we engineered a Rhodotorula paludigensis epoxide hydrolase (RpEH), obtaining mutant Escherichia coli/RpehH336W/L360F with a 6.4-fold enhanced enantiomeric ratio (E) from 5.5 to 35.4. This enabled a gram-scale resolution of rac-1a by E. coli/RpehH336W/L360F, producing (S)-1a (98.2% ees) and (R)-1b (75.0% eep) with 44.3 and 55.7% analytical yields, respectively. As follows, chiral (S)-flutriafol (98.2% ee) and (R)-flutriafol (75.0% ee) were easily synthesized by a one-step chemocatalytic process from (S)-1a and a two-step chemocatalytic process from (R)-1b, respectively. This chemoenzymatic approach offers a superior alternative for the asymmetric synthesis of flutriafol enantiomers. Furthermore, molecular dynamics simulations revealed insight into the enantioselectivity improvement of RpEH toward bulky 2,2-bisaryl-substituted oxirane 1a.

Single-Atom Ce-Doped Metal Hydrides with High Phosphatase-like Activity Amplify Oxidative Stress-Induced Tumor Apoptosis

Phosphates within tumors function as key biomolecules, playing a significant role in sustaining the viability of tumors. To disturb the homeostasis of cancer cells, regulating phosphate within the organism proves to be an effective strategy. Herein, we report single-atom Ce-doped Pt hydrides (Ce/Pt-H) with high phosphatase-like activity for phosphate hydrolysis. The resultant Ce/Pt-H exhibits a 26.90- and 6.25-fold increase in phosphatase-like activity in comparison to Ce/Pt and Pt-H, respectively. Mechanism investigations elucidate that the Ce Lewis acid site facilitates the coordination with phosphate groups, while the surface hydrides enhance the electron density of Pt for promoting catalytic ability in H2O cleavage and subsequent nucleophilic attack of hydroxyl groups. Finally, by leveraging its phosphatase-like activity, Ce/Pt-H can effectively regulate intracellular phosphates to disrupt redox homeostasis and amplify oxidative stress within cancer cells, ultimately leading to tumor apoptosis. This work provides fresh insights into noble-metal-based phosphatase mimics for inducing tumor apoptosis.

Highly sensitive hydrolytic nanozyme-based sensors for colorimetric detection of aluminum ions

Hydrolytic nanozyme-based visual colorimetry has emerged as a promising strategy for the detection of aluminum ions. However, most studies focus on simulating the structure of natural enzymes while neglecting to regulate the rate of hydrolysis-related steps, leading to low enzyme-like activity for hydrolytic nanozymes. Herein, we constructed a ruthenium dioxide (RuO2) in situ embedded cerium oxide (CeO2) nanozyme (RuO2/CeO2) with a Lewis acid-base pair (Ce-O-Ru-OH), which can simulate the catalytic behavior of phosphatase (PPase) and can be quantitatively quenched by Al3+ to achieve accurate and sensitive Al3+ colorimetric sensing detection. The incorporation of Ru into CeO2 nanorods accelerates the dissociation of H2O, followed by subsequent combination of hydroxide species to Lewis acidic Ce-O sites. This synergistic effect facilitates substrate activation and significantly enhances the hydrolysis activity of the nanozyme. The results show that the RuO2/CeO2 nanozyme exhibits a limit of detection as low as 0.5 ng/mL. We also demonstrate their efficacy in detecting Al3+ in various practical food samples. This study offers novel insights into the advancement of highly sensitive hydrolytic nanozyme engineering for sensing applications.

Biomimetic single Al-OH site with high acetylcholinesterase-like activity and self-defense ability for neuroprotection

Neurotoxicity of organophosphate compounds (OPs) can catastrophically cause nervous system injury by inhibiting acetylcholinesterase (AChE) expression. Although artificial systems have been developed for indirect neuroprotection, they are limited to dissociating P-O bonds for eliminating OPs. However, these systems have failed to overcome the deactivation of AChE. Herein, we report our finding that Al3+ is engineered onto the nodes of metal-organic framework to synthesize MOF-808-Al with enhanced Lewis acidity. The resultant MOF-808-Al efficiently mimics the catalytic behavior of AChE and has a self-defense ability to break the activity inhibition by OPs. Mechanism investigations elucidate that Al3+ Lewis acid sites with a strong polarization effect unite the highly electronegative -OH groups to form the enzyme-like catalytic center, resulting in superior substrate activation and nucleophilic attack ability with a 2.7-fold activity improvement. The multifunctional MOF-808-Al, which has satisfactory biosafety, is efficient in reducing neurotoxic effects and preventing neuronal tissue damage.

Evolution of tunnels in α/β-hydrolase fold proteins—What can we learn from studying epoxide hydrolases?

The evolutionary variability of a protein’s residues is highly dependent on protein region and function. Solvent-exposed residues, excluding those at interaction interfaces, are more variable than buried residues whereas active site residues are considered to be conserved. The abovementioned rules apply also to α/β-hydrolase fold proteins—one of the oldest and the biggest superfamily of enzymes with buried active sites equipped with tunnels linking the reaction site with the exterior. We selected soluble epoxide hydrolases as representative of this family to conduct the first systematic study on the evolution of tunnels. We hypothesised that tunnels are lined by mostly conserved residues, and are equipped with a number of specific variable residues that are able to respond to evolutionary pressure. The hypothesis was confirmed, and we suggested a general and detailed way of the tunnels’ evolution analysis based on entropy values calculated for tunnels’ residues. We also found three different cases of entropy distribution among tunnel-lining residues. These observations can be applied for protein reengineering mimicking the natural evolution process. We propose a ‘perforation’ mechanism for new tunnels design via the merging of internal cavities or protein surface perforation. Based on the literature data, such a strategy of new tunnel design could significantly improve the enzyme’s performance and can be applied widely for enzymes with buried active sites.

So far very little is known about proteins tunnels evolution. The goal of this study is to evaluate the evolution of tunnels in the family of soluble epoxide hydrolases—representatives of numerous α/β-hydrolase fold enzymes. As a result two types of tunnels evolution analysis were proposed (a general and a detailed approach), as well as a ‘perforation’ mechanism which can mimic native evolution in proteins and can be used as an additional strategy for enzymes redesign.

Unraveling the mechanism of enantio-controlling switches of an alcohol dehydrogenase toward sterically small ketone

Efficient synthesis of chiral compounds under mild conditions is highly desirable in the chemical and pharmaceutical communities, but it often faces difficulties. Although various enzymes have been harnessed as biocatalysts...

Nearly perfect kinetic resolution of racemic o-nitrostyrene oxide by AuEH2, a microsomal epoxide hydrolase from Aspergillus usamii, with high enantio- and regio-selectivity

Only a few known epoxide hydrolases (EHs) displayed activity towards o-nitrostyrene oxide (4a), presumably owing to the large steric hindrance caused by o-nitro substituent. Therefore, excavating EHs with high activity and enantio- and/or regio-selectivity towards racemic (rac-) 4a is essential but challenging. Here, AuEH2 from Aspergillus usamii was expressed in E. coli BL21(DE3). E. coli/Aueh2, an E. coli transformant expressing AuEH2, possessed EH activities of 16.2-184 U/g wet cell towards rac-styrene oxide (1a) and its derivatives (2a-13a), and the largest enantiomeric ratio of 96 towards rac-4a. The regioselectivity coefficients, β_R and β_S, of AuEH2 were determined to be 99.2% and 98.9%, suggesting that it regiopreferentially attacks the C_β in the oxirane rings of (R)- and (S)-4a. Then, the nearly perfect kinetic resolution of 20 mM rac-4a in pure water was carried out using 20 mg/mL wet cells of E. coli/Aueh2 at 25 °C for 50 min, retaining (S)-4a with over 99% ee_s and 48.9% yield_s, while producing (R)-o-nitrophenyl-1,2-ethanediol (4b) with 95.3% ee_p and 49.8% yield_p. To elucidate the molecular mechanism of AuEH2 with high enantiopreference for (R)-4a, its crystal structure was solved by X-ray diffraction and the molecular docking of AuEH2 with (R)- or (S)-4a was simulated.

Small Molecule Soluble Epoxide Hydrolase Inhibitors in Multitarget and Combination Therapies for Inflammation and Cancer

The enzyme soluble epoxide hydrolase (sEH) plays a central role in metabolism of bioactive lipid signaling molecules. The substrate-specific hydrolase activity of sEH converts epoxyeicosatrienoic acids (EETs) to less bioactive dihydroxyeicosatrienoic acids. EETs exhibit anti-inflammatory, analgesic, antihypertensive, cardio-protective and organ-protective properties. Accordingly, sEH inhibition is a promising therapeutic strategy for addressing a variety of diseases. In this review, we describe small molecule architectures that have been commonly deployed as sEH inhibitors with respect to angiogenesis, inflammation and cancer. We juxtapose commonly used synthetic scaffolds and natural products within the paradigm of a multitarget approach for addressing inflammation and inflammation induced carcinogenesis. Structural insights from the inhibitor complexes and novel strategies for development of sEH-based multitarget inhibitors are also presented. While sEH inhibition is likely to suppress inflammation-induced carcinogenesis, it can also lead to enhanced angiogenesis via increased EET concentrations. In this regard, sEH inhibitors in combination chemotherapy are described. Urea and amide-based architectures feature prominently across multitarget inhibition and combination chemotherapy applications of sEH inhibitors.

Manipulating the regioselectivity of a Solanum lycopersicum epoxide hydrolase for the enantioconvergent synthesis of enantiopure alkane- and alkene-1,2-diols

This work engineered a superior double-site mutant SlEH1^W106T/F189L used for the enantioconvergent biosynthesis of (R)-1b–6b with high ee_p values.

Exploring the origins of selectivity in soluble epoxide hydrolase from Bacillus megaterium

Epoxide hydrolase (EH) enzymes catalyze the hydration of racemic epoxides to yield their corresponding vicinal diols. In this work, the Bacillus megaterium epoxide hydrolase (BmEH)-mediated hydrolysis of racemic styrene oxide (rac-SO) and its para-nitro styrene oxide (rac-p-NSO) derivative are computationally investigated using density functional theory (DFT).

Epoxide hydrolase (EH) enzymes catalyze the hydration of racemic epoxides to yield their corresponding vicinal diols. These enzymes present different enantio- and regioselectivity depending upon either the substrate structure or the substitution pattern of the epoxide ring. In this study, we computationally investigate the Bacillus megaterium epoxide hydrolase (BmEH)-mediated hydrolysis of racemic styrene oxide (rac-SO) and its para-nitro styrene oxide (rac-p-NSO) derivative using density functional theory (DFT) and an active site cluster model consisting of 195 and 197 atoms, respectively. Full reaction mechanisms for epoxide ring opening were evaluated considering the attack at both oxirane carbons and considering two possible orientations of the substrate at the BmEH active site. Our results indicate that for both SO and p-NSO substrates the BmEH enantio- and regioselectivity is opposite to the inherent (R)-BmEH selectivity, the attack at the benzylic position (C1) of the (S)-enantiomer being the most favoured chemical outcome.

Evaluating Ylehd, a recombinant epoxide hydrolase from Yarrowia lipolytica as a potential biocatalyst for the resolution of benzyl glycidyl ether

Glycidyl ethers and their vicinal diols are important building blocks in the organic synthesis of anti-cancer and anti-obesity drugs. Ylehd, an epoxide hydrolase from tropical marine yeast Yarrowia lipolytica, was explored for its enantioselective properties by kinetic, thermodynamic and in silico studies. Kinetic resolution of racemic phenyl glycidyl ether (PGE) yielded (S)-epoxide while for benzyl glycidyl ether (BGE) (R)-epoxide was obtained, with vicinal diols of the opposite configuration. Amongst the enantiomers of PGE and BGE, the (S)-selective conversion of benzyl glycidyl ether to its corresponding diol, (S)-3-benzyloxy-1,2-propanediol while retaining (R)-BGE was most favourable with 95% ee in 20 min. Enantioselective conversion of specific enantiomer of BGE to its corresponding diols was attributed to the favourable kinetic and thermodynamic parameters as well as to the number and proximity of water molecules near the base H325 in the active site pocket. The easily available and highly active Ylehd could be a potential biocatalyst for large scale preparation of pharmaceutically relevant chiral (R)-benzyl glycidyl ether and (S)-3-benzyloxy-1,2-propanediol.

Ylehd, an enantioselective epoxide hydrolase with potential application in resolution of racemic benzyl glycidyl ether.

Active-Site Flexibility and Substrate Specificity in a Bacterial Virulence Factor: Crystallographic Snapshots of an Epoxide Hydrolase

Pseudomonas aeruginosa secretes an epoxide hydrolase with catalytic activity that triggers degradation of the cystic fibrosis transmembrane conductance regulator (CFTR) and perturbs other host defense networks. Targets of this CFTR inhibitory factor (Cif) are largely unknown, but include an epoxy-fatty acid. In this class of signaling molecules, chirality can be an important determinant of physiological output and potency. Here we explore the active-site chemistry of this two-step α/β-hydrolase and its implications for an emerging class of virulence enzymes. In combination with hydrolysis data, crystal structures of 15 trapped hydroxyalkyl-enzyme intermediates reveal the stereochemical basis of Cif's substrate specificity, as well as its regioisomeric and enantiomeric preferences. The structures also reveal distinct sets of conformational changes that enable the active site to expand dramatically in two directions, accommodating a surprising array of potential physiological epoxide targets. These new substrates may contribute to Cif's diverse effects in vivo, and thus to the success of P. aeruginosa and other pathogens during infection.

Consistent Principal Component Modes from Molecular Dynamics Simulations of Proteins

Principal component analysis is a technique widely used for studying the movements of proteins using data collected from molecular dynamics simulations. In spite of its extensive use, the technique has a serious drawback: equivalent simulations do not afford the same PC-modes. In this article, we show that concatenating equivalent trajectories and calculating the PC-modes from the concatenated one significantly enhances the reproducibility of the results. Moreover, the consistency of the modes can be systematically improved by adding more individual trajectories to the concatenated one.

Investigation of the binding mode of 1, 3, 4-oxadiazole derivatives as amide-based inhibitors for soluble epoxide hydrolase (sEH) by molecular docking and MM-GBSA

The soluble epoxide hydrolase (sEH) enzyme plays an important role in the metabolism of endogenous chemical mediators involved in the regulation of blood pressure and inflammation. Inhibition of sEH provides a new approach to the treatment of inflammation, hypertension and atherosclerosis. In this study, the binding modes and inhibition mechanisms of the new oxadiazole-based amide inhibitors of the human soluble epoxide hydrolase were investigated by molecular docking and molecular dynamics (MD) simulation followed by the MM-GBSA method to calculate the binding free energy of each inhibitor to sEH. The results obtained from the binding free energy (ΔG _binding) calculation and normal mode analysis indicate that the major favorable contributors are the van der Waals and electrostatic terms, whereas the polar solvation term opposes binding. In addition, a good agreement between the calculated ΔG _binding and the experimental IC₅₀ was obtained [correlation coefficient, r ² = 0.89 (with) and 0.87 (without) entropy]. Besides, comparison of the enthalpy changes (ΔG _MM-GBSA) with entropy changes (-TΔS) indicates that binding process of all inhibitors to sEH is enthalpy-driven. Based on the ΔG _binding on per residue decomposition, Asp335 and Tyr383 residues from the active site and Trp336, Leu499 and His524 residues from hydrophobic pockets contribute the most to ΔG _binding. Moreover, hydrogen bond analysis reveals that Tyr383, Tyr466 and Asp335 residues have an important role in the binding to inhibitors by forming hydrogen bonds with high occupancies. Our obtained results are useful for the understanding of the sEH-inhibitor interactions and may have great importance in the design of future sEH inhibitors.

Is Promiscuous CALB a Good Scaffold for Designing New Epoxidases?

Candida Antarctica lipase B (CALB) is a well-known enzyme, especially because of its promiscuous activity. Due to its properties, CALB was widely used as a benchmark for designing new catalysts for important organic reactions. The active site of CALB is very similar to that of soluble epoxide hydrolase (sEH) formed by a nucleophile-histidine-acid catalytic triad and an oxyanion hole typical for molecular structures derived from processes of α/β hydrolases. In this work we are exploring these similarities and proposing a Ser105Asp variant of CALB as a new catalyst for epoxide hydrolysis. In particular, the hydrolysis of the trans-diphenylpropene oxide (t-DPPO) is studied by means of quantum cluster models mimicking the active site of both enzymes. Our results, based on semi-empirical and DFT calculations, suggest that mutant Ser105Asp CALB is a good protein scaffold to be used for the bio-synthesis of chiral compounds.

Enantioselective enzymes by computational design and in silico screening

Computational enzyme design holds great promise for providing new biocatalysts for synthetic chemistry. A strategy to design small mutant libraries of complementary enantioselective epoxide hydrolase variants for the production of highly enantioenriched (S,S)-diols and (R,R)-diols is developed. Key features of this strategy (CASCO, catalytic selectivity by computational design) are the design of mutations that favor binding of the substrate in a predefined orientation, the introduction of steric hindrance to prevent unwanted substrate binding modes, and ranking of designs by high-throughput molecular dynamics simulations. Using this strategy we obtained highly stereoselective mutants of limonene epoxide hydrolase after experimental screening of only 37 variants. The results indicate that computational methods can replace a substantial amount of laboratory work when developing enantioselective enzymes.

Engineering the epoxide hydrolase from Agromyces mediolanus for enhanced enantioselectivity and activity in the kinetic resolution of racemic epichlorohydrin

The biocatalytic production of enantiopure epichlorohydrin (ECH) has been steadily attracting more attention.

Predictive model for epoxide hydrolase-generated stereochemistry in the biosynthesis of nine-membered enediyne antitumor antibiotics

Nine-membered enediyne antitumor antibiotics C-1027, neocarzinostatin (NCS), and kedarcidin (KED) possess enediyne cores to which activity-modulating peripheral moieties are attached via (R)- or (S)-vicinal diols. We have previously shown that this stereochemical difference arises from hydrolysis of epoxide precursors by epoxide hydrolases (EHs) with different regioselectivities. The inverting EHs, such as SgcF, hydrolyze an (S)-epoxide substrate to yield an (R)-diol in C-1027 biosynthesis, whereas the retaining EHs, such as NcsF2 and KedF, hydrolyze an (S)-epoxide substrate to yield an (S)-diol in NCS and KED biosynthesis. We now report the characterization of a series of EH mutants and provide a predictive model for EH regioselectivity in the biosynthesis of the nine-membered enediyne antitumor antibiotics. A W236Y mutation in SgcF increased the retaining activity toward (S)-styrene oxide by 3-fold, and a W236Y/Q237M double mutation in SgcF, mimicking NcsF2 and KedF, resulted in a 20-fold increase in the retaining activity. To test the predictive utility of these mutations, two putative enediyne biosynthesis-associated EHs were identified by genome mining and confirmed as inverting enzymes, SpoF from Salinospora tropica CNB-440 and SgrF (SGR_625) from Streptomyces griseus IFO 13350. Finally, phylogenetic analysis of EHs revealed a familial classification according to inverting versus retaining activity. Taken together, these results provide a predictive model for vicinal diol stereochemistry in enediyne biosynthesis and set the stage for further elucidating the origins of EH regioselectivity.

Computational protein engineering: bridging the gap between rational design and laboratory evolution

Enzymes are tremendously proficient catalysts, which can be used as extracellular catalysts for a whole host of processes, from chemical synthesis to the generation of novel biofuels. For them to be more amenable to the needs of biotechnology, however, it is often necessary to be able to manipulate their physico-chemical properties in an efficient and streamlined manner, and, ideally, to be able to train them to catalyze completely new reactions. Recent years have seen an explosion of interest in different approaches to achieve this, both in the laboratory, and in silico. There remains, however, a gap between current approaches to computational enzyme design, which have primarily focused on the early stages of the design process, and laboratory evolution, which is an extremely powerful tool for enzyme redesign, but will always be limited by the vastness of sequence space combined with the low frequency for desirable mutations. This review discusses different approaches towards computational enzyme design and demonstrates how combining newly developed screening approaches that can rapidly predict potential mutation “hotspots” with approaches that can quantitatively and reliably dissect the catalytic step can bridge the gap that currently exists between computational enzyme design and laboratory evolution studies.

Obtaining optical purity for product diols in enzyme-catalyzed epoxide hydrolysis: contributions from changes in both enantio- and regioselectivity

Enzyme variants of the plant epoxide hydrolase StEH1 displaying improved stereoselectivities in the catalyzed hydrolysis of (2,3-epoxypropyl)benzene were generated by directed evolution. The evolution was driven by iterative saturation mutagenesis in combination with enzyme activity screenings where product chirality was the decisive selection criterion. Analysis of the underlying causes of the increased diol product ratios revealed two major contributing factors: increased enantioselectivity for the corresponding epoxide enantiomer(s) and, in some cases, a concomitant change in regioselectivity in the catalyzed epoxide ring-opening half-reaction. Thus, variant enzymes that catalyzed the hydrolysis of racemic (2,3-epoxypropyl)benzene into the R-diol product in an enantioconvergent manner were isolated.

Determinants of reactivity and selectivity in soluble epoxide hydrolase from quantum mechanics/molecular mechanics modeling

Soluble epoxide hydrolase (sEH) is an enzyme involved in drug metabolism that catalyzes the hydrolysis of epoxides to form their corresponding diols. sEH has a broad substrate range and shows high regio- and enantioselectivity for nucleophilic ring opening by Asp333. Epoxide hydrolases therefore have potential synthetic applications. We have used combined quantum mechanics/molecular mechanics (QM/MM) umbrella sampling molecular dynamics (MD) simulations (at the AM1/CHARMM22 level) and high-level ab initio (SCS-MP2) QM/MM calculations to analyze the reactions, and determinants of selectivity, for two substrates: trans-stilbene oxide (t-SO) and trans-diphenylpropene oxide (t-DPPO). The calculated free energy barriers from the QM/MM (AM1/CHARMM22) umbrella sampling MD simulations show a lower barrier for phenyl attack in t-DPPO, compared with that for benzylic attack, in agreement with experiment. Activation barriers in agreement with experimental rate constants are obtained only with the highest level of QM theory (SCS-MP2) used. Our results show that the selectivity of the ring-opening reaction is influenced by several factors, including proximity to the nucleophile, electronic stabilization of the transition state, and hydrogen bonding to two active site tyrosine residues. The protonation state of His523 during nucleophilic attack has also been investigated, and our results show that the protonated form is most consistent with experimental findings. The work presented here illustrates how determinants of selectivity can be identified from QM/MM simulations. These insights may also provide useful information for the design of novel catalysts for use in the synthesis of enantiopure compounds.

Temperature and pH dependence of enzyme-catalyzed hydrolysis of trans-methylstyrene oxide. A unifying kinetic model for observed hysteresis, cooperativity, and regioselectivity

The underlying enzyme kinetics behind the regioselective promiscuity shown by epoxide hydrolases toward certain epoxides has been studied. The effects of temperature and pH on regioselectivity were investigated by analyzing the stereochemistry of hydrolysis products of (1R,2R)-trans-2-methylstyrene oxide between 14-46 degrees C and pH 6.0-9.0, either catalyzed by the potato epoxide hydrolase StEH1 or in the absence of enzyme. In the enzyme-catalyzed reaction, a switch of preferred epoxide carbon that is subjected to nucleophilic attack is observed at pH values above 8. The enzyme also displays cooperativity in substrate saturation plots when assayed at temperatures < or = 30 degrees C and at intermediate pH. The cooperativity is lost at higher assay temperatures. Cooperativity can originate from a kinetic mechanism involving hysteresis and will be dependent on the relationship between k(cat) and the rate of interconversion between two different Michaelis complexes. In the case of the studied reactions, the proposed different Michaelis complexes are enzyme-substrate complexes in which the epoxide substrate is bound in different binding modes, allowing for separate pathways toward product formation. The assumption of separated, but interacting, reaction pathways is supported by that formation of the two product enantiomers also displays distinct pH dependencies of k(cat)/K(M). The thermodynamic parameters describing the differences in activation enthalpy and entropy suggest that (1) regioselectivity is primarily dictated by differences in activation entropy with positive values of both DeltaDeltaH(++) and DeltaDeltaS(++) and (2) the hysteretic behavior is linked to an interconversion between Michaelis complexes with rates increasing with temperature. From the collected data, we propose that hysteresis, regioselectivity, and, when applicable, hysteretic cooperativity are closely linked properties, explained by the kinetic mechanism earlier introduced by our group.

Understanding the mode of action of ThDP in benzoylformate decarboxylase

The mechanism of all elementary steps involved in the catalytic cycle of benzoylformate decarboxylase (BFD, E.C. 4.1.1.7) to generate the acyloin linkage is investigated by extensive molecular dynamics simulations. Models involving different charge states of amino acids and/or mutants of critical residues were constructed to understand the involvement of the catalytically active residues and the reactivity differences between different substrates in this reaction. Our calculations confirm that H70, S26, and H281 are catalytically active amino acids. H281 functions as a base to accept H(donor) in the first nucleophilic attack and as an acid in the second, to donate the proton back to O(acceptor). S26 assists H281 in deprotonation of the donor aldehyde and protonation of the acceptor aldehyde. In both the first and second nucleophilic attacks, H70 interacts with O(aldehyde) and aligns it toward the nucleophilic center. H70 has been found to have an electrostatic effect on the approaching aldehyde whose absence would block the initiation of the reaction. The reactivity difference between benzaldehyde (BA) and acetaldehyde (AA) is mainly explained by the steric interactions of the acceptor aldehyde with the surrounding amino acids in the active center of the enzyme.

Directed evolution of an enantioselective epoxide hydrolase: uncovering the source of enantioselectivity at each evolutionary stage

Directed evolution of enzymes as enantioselective catalysts in organic chemistry is an alternative to traditional asymmetric catalysis using chiral transition-metal complexes or organocatalysts, the different approaches often being complementary. Moreover, directed evolution studies allow us to learn more about how enzymes perform mechanistically. The present study concerns a previously evolved highly enantioselective mutant of the epoxide hydrolase from Aspergillus niger in the hydrolytic kinetic resolution of racemic glycidyl phenyl ether. Kinetic data, molecular dynamics calculations, molecular modeling, inhibition experiments, and X-ray structural work for the wild-type (WT) enzyme and the best mutant reveal the basis of the large increase in enantioselectivity (E = 4.6 versus E = 115). The overall structures of the WT and the mutant are essentially identical, but dramatic differences are observed in the active site as revealed by the X-ray structures. All of the experimental and computational results support a model in which productive positioning of the preferred (S)-glycidyl phenyl ether, but not the (R)-enantiomer, forms the basis of enhanced enantioselectivity. Predictions regarding substrate scope and enantioselectivity of the best mutant are shown to be possible.

A DFT study of solvation effects on the tautomeric equilibrium and catalytic ylide generation of thiamin models

Thiamin diphosphate (ThDP) is the biologically active form of vitamin B1 and an essential cofactor for a number of enzymes. The effect of solvent polarity on the tautomeric equilibria of ThDP using three model systems of the 4'-aminopyrimidine ring is studied by density functional theory calculations (B3LYP/6-311+G(d,p)//B3LYP/6-31G(d)) in the gas phase and selected solvents (cyclohexane, ether, dichloroethane, and water). Solvation effects are investigated using three different schemes: implicit solvation by a continuum model, explicit solvation by inclusion of three water molecules mimicking the first solvation shell of the enzymatic environment, and by a mixed implicit/explicit solvation model. The 4'-aminopyrimidine tautomer is more stable than the 1',4'-iminopyrimidine tautomer in all solvation schemes employed; however, the trend for the stabilities of the 1',4'-iminopyrimidine tautomer in the solvents depends on the specific ThDP-model. Formation of the catalytic important ylide for ThDP-dependent enzymes by deprotonation of ThDP(C2) is also investigated by localization of transition states for two possible pathways. Only the less stable tautomer, 1',4'-iminopyrimidine ThDP, is able to form the catalytic active ylide. Generation of the ylide through a direct intramolecular proton transfer from ThDP(C2) to the ThDP(N4') nitrogen lone pair is favored by 6 kcal/mol in the gas phase, as compared to a water-mediated ylide generation. However, inclusion of a dielectric medium reduces this difference dramatically. Furthermore, inclusion of two water molecules to model the apoenzymatic environment lowers the activation energies of both direct and water-mediated ylide generation.

Proton shuttles and phosphatase activity in soluble epoxide hydrolase

Recently, a novel metal Mg2+-dependent phosphatase activity has been discovered in the N-terminal domain of the soluble epoxide hydrolase (sEH), opening a new branch of fatty acid metabolism and providing an additional site for drug targeting. Importantly, the sEH N-terminal fold belongs to the haloacid dehalogenase (HAD) superfamily, which comprises a vast majority of phosphotransferases. Herein, we present the results of a computational study of the sEH phosphatase activity, which includes classical molecular dynamics (MD) simulations and mixed quantum mechanical/molecular mechanics (QM/MM) calculations. On the basis of experimental results, a two-step mechanism has been proposed and herein investigated: (1) phosphoenzyme intermediate formation and (2) phosphoenzyme intermediate hydrolysis. Building on our earlier work, we now provide a detailed description of the reaction mechanism for the whole catalytic cycle along with its free energy profile. The present computations suggest metaphosphate-like transition states for these phosphoryl transfers. They also reveal that the enzyme promotes water deprotonation and facilitates shuttling of protons via a metal-ligand connecting water bridge (WB). These WB-mediated proton shuttles are crucial for the activation of the solvent nucleophile and for the stabilization of the leaving group. Moreover, due to the conservation of structural features in the N-terminal catalytic site of sEH and other members of the HAD superfamily, we suggest a generalization of our findings to these other metal-dependent phosphatases.

Conformational dynamics of the estrogen receptor alpha: molecular dynamics simulations of the influence of binding site structure on protein dynamics

We present 158 ns of unrestrained all-atom molecular dynamics (MD) simulations of the human estrogen receptor alpha ligand binding domain (ERalpha LBD) sampling the conformational changes upon binding of estradiol. The pivotal role of His524 in maintaining the protein structure in the biologically active agonist conformation is elucidated. With His524 modeled as the epsilon-tautomer, a conserved hydrogen bond to the ligand is found in the active complex. Helices 3 and 11 are held together by a hydrogen-bonding network from His524 to Glu339 via Glu419 and Lys531, arresting the ligand in the binding pocket and creating the "mouse trap" binding site for helix 12 (H12). The simulations reveal how His524 serves as a communication point between the two. When estradiol is bound, His524 is positioned correctly for the hydrogen bond network to be established. H12 is then positioned for interaction with the co-activator protein, leading to the biologically active complex. The conformational dynamics of ERalpha LBD is further investigated from simulations of antagonist and apo conformations of the protein. These simulations suggest a likely sequence of events for the transition from the inactive apo structure to the transcriptionally active conformation of ERalpha LBD. Stable conformations are identified where H12 is placed neither in the "mouse trap" nor in the co-activator binding groove, as is the case for antagonist structures of ERalpha LBD. Finally, the influence of such conformations on the biological function of ERalpha is discussed in relationship to the interaction with selective estrogen receptor modulators and endocrine-disrupting compounds.