Theoretical Study of the Full Reaction Mechanism of Human Soluble Epoxide Hydrolase

Investigation of iron release from the N- and C-lobes of human serum transferrin by quantum chemical calculations

Cluster models of iron binding sites of the N- and C-lobes highlights the inequivalence of each lobe in iron release.

Computer-aided understanding and engineering of enzymatic selectivity

Enzymes offering chemo-, regio-, and stereoselectivity enable the asymmetric synthesis of high-value chiral molecules. Unfortunately, the drawback that naturally occurring enzymes are often inefficient or have undesired selectivity toward non-native substrates hinders the broadening of biocatalytic applications. To match the demands of specific selectivity in asymmetric synthesis, biochemists have implemented various computer-aided strategies in understanding and engineering enzymatic selectivity, diversifying the available repository of artificial enzymes. Here, given that the entire asymmetric catalytic cycle, involving precise interactions within the active pocket and substrate transport in the enzyme channel, could affect the enzymatic efficiency and selectivity, we presented a comprehensive overview of the computer-aided workflow for enzymatic selectivity. This review includes a mechanistic understanding of enzymatic selectivity based on quantum mechanical calculations, rational design of enzymatic selectivity guided by enzyme-substrate interactions, and enzymatic selectivity regulation via enzyme channel engineering. Finally, we discussed the computational paradigm for designing enzyme selectivity in silico to facilitate the advancement of asymmetric biosynthesis.

Computational insights into the known inhibitors of human soluble epoxide hydrolase

Human soluble epoxide hydrolase (hsEH) is involved in the hydrolysis of epoxyeicosatrienoic acids (EETs), which have potent anti-inflammatory properties. Given that EET conversion generates nonbioactive molecules, inhibition of this enzyme would be beneficial. Past decades of work on hsEH inhibitors resulted in numerous potential compounds, of which a hundred hsEH-ligand complexes were crystallized and deposited in the Protein Data Bank (PDB). We analyzed all deposited hsEH-ligand complexes to gain insight into the binding of inhibitors and to provide feedback on the future drug design processes. We also reviewed computationally driven strategies that were used to propose novel hsEH inhibitors.

Identification and catalytic properties of new epoxide hydrolases from the genomic data of soil bacteria

Epoxide hydrolases (EHs) catalyse the conversion of epoxides into vicinal diols. These enzymes have extensive value in biocatalysis as they can generate enantiopure epoxides and diols which are important and versatile synthetic intermediates for the fine chemical and pharmaceutical industries. Despite these benefits, they have seen limited use in the bioindustry and novel EHs continue to be reported in the literature. We identified twenty-nine putative EHs within the genomes of soil bacteria. Eight of these EHs were explored in terms of their activity. Two limonene epoxide hydrolases (LEHs) and one ⍺/β EH were active on a model compound styrene oxide and its ring-substituted derivatives, with low to good percentage conversions of 18-86%. Further exploration of the substrate scope with enantiopure (R)-styrene oxide and (S)-styrene oxide, showed different epoxide ring opening regioselectivities. Two enzymes, expressed from plasmids pQR1984 and pQR1990 de-symmetrised the meso-epoxide cyclohexene oxide, forming the (R,R)-diol with high enantioselectivity. Two LEHs, from plasmids pQR1980 and pQR1982 catalysed the hydrolysis of (+) and (-) limonene oxide, with diastereomeric preference for the (1S,2S,4R)- and (1R,2R,4S)-diol products, respectively. The enzyme from plasmid pQR1982 had a good substrate scope for a LEH, being active towards styrene oxide, its analogues, cyclohexene oxide and 1,2-epoxyhexane in addition to (±)-limonene oxide. The enzymes from plasmids pQR1982 and pQR1984 had good substrate scopes and their enzymatic properties were characterised with respect to styrene oxide. They had comparable temperature optima and pQR1984 had 70% activity in the presence of 40% of the green solvent MeOH, a useful property for bio-industrial applications. Overall, this study has provided novel EHs with potential value in industrial biocatalysis.

Applications of water molecules for analysis of macromolecule properties

Water molecules maintain proteins' structures, functions, stabilities and dynamics. They can occupy certain positions or pass quickly via a protein's interior. Regardless of their behaviour, water molecules can be used for the analysis of proteins' structural features and biochemical properties. Here, we present a list of several software programs that use the information provided by water molecules to: i) analyse protein structures and provide the optimal positions of water molecules for protein hydration, ii) identify high-occupancy water sites in order to analyse ligand binding modes, and iii) detect and describe tunnels and cavities. The analysis of water molecules' distribution and trajectories sheds a light on proteins' interactions with small molecules, on the dynamics of tunnels and cavities, on protein composition and also on the functionality, transportation network and location of functionally relevant residues. Finally, the correct placement of water molecules in protein crystal structures can significantly improve the reliability of molecular dynamics simulations.

15-deoxy-Δ12,14-Prostaglandin J2 inhibits human soluble epoxide hydrolase by a dual orthosteric and allosteric mechanism

Human soluble epoxide hydrolase (hsEH) is an enzyme responsible for the inactivation of bioactive epoxy fatty acids, and its inhibition is emerging as a promising therapeutical strategy to target hypertension, cardiovascular disease, pain and insulin sensitivity. Here, we uncover the molecular bases of hsEH inhibition mediated by the endogenous 15-deoxy-Δ12,14-Prostaglandin J2 (15d-PGJ2). Our data reveal a dual inhibitory mechanism, whereby hsEH can be inhibited by reversible docking of 15d-PGJ2 in the catalytic pocket, as well as by covalent locking of the same compound onto cysteine residues C423 and C522, remote to the active site. Biophysical characterisations allied with in silico investigations indicate that the covalent modification of the reactive cysteines may be part of a hitherto undiscovered allosteric regulatory mechanism of the enzyme. This study provides insights into the molecular modes of inhibition of hsEH epoxy-hydrolytic activity and paves the way for the development of new allosteric inhibitors.

Expression, purification, and characterisation of human soluble Epoxide Hydrolase (hsEH) and of its functional C-terminal domain

The human soluble Epoxide Hydrolase (hsEH) is an enzyme involved in the hydrolysis of endogenous anti-inflammatory and cardio-protective signalling mediators known as epoxyeicosatrienoic acids (EETs). EETs’ conversion into the corresponding diols by hsEH generates non-bioactive molecules, thereby the enzyme inhibition would be expected to enhance the EETs bioavailability, and their beneficial properties. Numerous inhibitors have been developed to target the enzyme, some of which are showing promising antihypertensive and anti-inflammatory properties in vivo. Thus far, the preparation of the recombinant enzyme for enzymatic and structural in vitro studies has been performed mainly using a baculovirus expression system. More recently, it was reported that the enzyme could be exogenously expressed and isolated from E. coli, although limited amounts of active protein were obtained. We herein describe two novel methods to yield pure recombinant enzyme. The first describes the expression and purification of the full-length enzyme from eukaryotic cells HEK293-F, whilst the second concerns the C-terminal domain of hsEH obtained from the cost-effective and rapid E. coli prokaryotic system. The two methods successfully generated satisfactory amounts of functional enzyme, with virtually identical enzymatic activity. Overall, the protocols described in this paper can be employed for the recombinant expression and purification of active hsEH, to be used in future biomedical investigations and for high-throughput screening of inhibitors for potential use in the treatment of cardiovascular disease.

•

hsEH is a key regulator of cardiovascular homeostasis.

•

A HEK293-F mammalian expression system for hsEH full-length (FL) was developed.

•

An E. coli expression system for the hsEH C-terminal Domain (CTD) was established.

•

Both proteins exhibited the same enzymatic specific activity in vitro.

•

The CTD preparation provides benefits of easy operation, and high yield and purity.

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.

Probing the orientation of inhibitor and epoxy-eicosatrienoic acid binding in the active site of soluble epoxide hydrolase

Soluble epoxide hydrolase (sEH) is an important therapeutic target of many diseases, such as chronic obstructive pulmonary disease (COPD) and diabetic neuropathic pain. It acts by hydrolyzing and thus regulating specific bioactive long chain polyunsaturated fatty acid epoxides (lcPUFA), like epoxyeicosatrienoic acids (EETs). To better predict which epoxides could be hydrolyzed by sEH, one needs to dissect the important factors and structural requirements that govern the binding of the substrates to sEH. This knowledge allows further exploration of the physiological role played by sEH. Unfortunately, a crystal structure of sEH with a substrate bound has not yet been reported. In this report, new photoaffinity mimics of a sEH inhibitor and EET regioisomers were prepared and used in combination with peptide sequencing and computational modeling, to identify the binding orientation of different regioisomers and enantiomers of EETs into the catalytic cavity of sEH. Results indicate that the stereochemistry of the epoxide plays a crucial role in dictating the binding orientation of the substrate.

Conformational diversity and enantioconvergence in potato epoxide hydrolase 1

Potato epoxide hydrolase 1 (StEH1) is a biocatalytically important enzyme that exhibits rich enantio- and regioselectivity in the hydrolysis of chiral epoxide substrates. In particular, StEH1 has been demonstrated to enantioconvergently hydrolyze racemic mixes of styrene oxide (SO) to yield (R)-1-phenylethanediol. This work combines computational, crystallographic and biochemical analyses to understand both the origins of the enantioconvergent behavior of the wild-type enzyme, as well as shifts in activities and substrate binding preferences in an engineered StEH1 variant, R-C1B1, which contains four active site substitutions (W106L, L109Y, V141K and I155V). Our calculations are able to reproduce both the enantio- and regioselectivities of StEH1, and demonstrate a clear link between different substrate binding modes and the corresponding selectivity, with the preferred binding modes being shifted between the wild-type enzyme and the R-C1B1 variant. Additionally, we demonstrate that the observed changes in selectivity and the corresponding enantioconvergent behavior are due to a combination of steric and electrostatic effects that modulate both the accessibility of the different carbon atoms to the nucleophilic side chain of D105, as well as the interactions between the substrate and protein amino acid side chains and active site water molecules. Being able to computationally predict such subtle effects for different substrate enantiomers, as well as to understand their origin and how they are affected by mutations, is an important advance towards the computational design of improved biocatalysts for enantioselective synthesis.

Quantum Chemical Modeling of the Dehalogenation Reaction of Haloalcohol Dehalogenase

The dehalogenation reaction of haloalcohol dehalogenase HheC from Agrobacterium radiobacter AD1 was investigated theoretically using hybrid density functional theory methods. HheC catalyzes the enantioselective conversion of halohydrins into their corresponding epoxides. The reaction is proposed to be mediated by a catalytic Ser132-Tyr145-Arg149 triad, and a distinct halide binding site is suggested to facilitate halide displacement by stabilizing the free ion. We investigated the HheC-mediated dehalogenation of (R)-2-chloro-1-phenylethanol using three quantum chemical models of various sizes. The calculated barriers and reaction energies give support to the suggested reaction mechanism. The dehalogenation occurs in a single concerted step, in which Tyr145 abstracts a proton from the halohydrin substrate and the substrate oxyanion displaces the chloride ion, forming the epoxide. Characterization of the involved stationary points is provided. Furthermore, by using three different models of the halide binding site, we are able to assess the adopted modeling methodology.

Expanding the Catalytic Triad in Epoxide Hydrolases and Related Enzymes

Potato epoxide hydrolase 1 exhibits rich enantio- and regioselectivity in the hydrolysis of a broad range of substrates. The enzyme can be engineered to increase the yield of optically pure products as a result of changes in both enantio- and regioselectivity. It is thus highly attractive in biocatalysis, particularly for the generation of enantiopure fine chemicals and pharmaceuticals. The present work aims to establish the principles underlying the activity and selectivity of the enzyme through a combined computational, structural, and kinetic study using the substrate trans-stilbene oxide as a model system. Extensive empirical valence bond simulations have been performed on the wild-type enzyme together with several experimentally characterized mutants. We are able to computationally reproduce the differences between the activities of different stereoisomers of the substrate and the effects of mutations of several active-site residues. In addition, our results indicate the involvement of a previously neglected residue, H104, which is electrostatically linked to the general base H300. We find that this residue, which is highly conserved in epoxide hydrolases and related hydrolytic enzymes, needs to be in its protonated form in order to provide charge balance in an otherwise negatively charged active site. Our data show that unless the active-site charge balance is correctly treated in simulations, it is not possible to generate a physically meaningful model for the enzyme that can accurately reproduce activity and selectivity trends. We also expand our understanding of other catalytic residues, demonstrating in particular the role of a noncanonical residue, E35, as a "backup base" in the absence of H300. Our results provide a detailed view of the main factors driving catalysis and regioselectivity in this enzyme and identify targets for subsequent enzyme design efforts.

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.

Epoxide hydrolase activities and epoxy fatty acids in the mosquito Culex quinquefasciatus

Culex mosquitoes have emerged as important model organisms for mosquito biology, and are disease vectors for multiple mosquito-borne pathogens, including West Nile virus. We characterized epoxide hydrolase activities in the mosquito Culex quinquefasciatus, which suggested multiple forms of epoxide hydrolases were present. We found EH activities on epoxy eicosatrienoic acids (EETs). EETs and other eicosanoids are well-established lipid signaling molecules in vertebrates. We showed EETs can be synthesized in vitro from arachidonic acids by mosquito lysate, and EETs were also detected in vivo both in larvae and adult mosquitoes by LC-MS/MS. The EH activities on EETs can be induced by blood feeding, and the highest activity was observed in the midgut of female mosquitoes. The enzyme activities on EETs can be inhibited by urea-based inhibitors designed for mammalian soluble epoxide hydrolases (sEH). The sEH inhibitors have been shown to play diverse biological roles in mammalian systems, and they can be useful tools to study the function of EETs in mosquitoes. Besides juvenile hormone metabolism and detoxification, insect epoxide hydrolases may also play a role in regulating lipid signaling molecules, such as EETs and other epoxy fatty acids, synthesized in vivo or obtained from blood feeding by female mosquitoes.

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.

Electronic and chelation effects on the unusual C2-methylation of N-(para-substituted)phenylaziridines with lithium organocuprates

Density functional theory calculations with the B3LYP functional were performed for the title ring-opening reaction to understand the intrinsic activating and directing effects of the N-substituents, as well as the electron donating effect of the para-substituted (Y = Cl, H, Me) phenyl group at the more hindered benzylic C2 atom. The N-tosyl group (i.e., N-Tos) or the N-(2-pyridyl)sulfonyl group (i.e., N-Py) was introduced to activate the ring nitrogen atom (N1) and the para-substituted (Y = Cl, H, Me) phenyl group for the activation of the C2 atom. Conformational searches and geometry optimizations were performed for the N-(para-substituted)phenylaziridines (1∼6). Calculations indicate that the aziridine 6 (i.e., Py/Me) has the most elongated C2-N1 bond intrinsically due to the electronic activating effects, implying the aziridine 6 to be the most potent candidate for the more-hindered C2 opening. Transition states (TSs) were investigated for the prospective ring-opening paths (I∼IV), considering the types of intermolecular push-pull interactions between the N-activated phenylaziridines and the cuprate. The N-Py group provides an unique C2-favored TS along the path IV, which the N-Tos group cannot afford, due to the less charge transfer from the nucleophilic CH 3δ- of the cuprate into the electrophilic C2 atom. Furthermore, the e-donating effect of the para-substituents (Y = Cl, H, Me) enhances the C2 opening for the path IV. This study enables us to understand the unusual ring-opening phenomena in terms of electronic and directing effects and hence may serve as a tool to design substrates for highly regioselective ring openings.

Quantum Chemical Modeling of Enzymatic Reactions: The Case of Decarboxylation

We present a systematic study of the decarboxylation step of the enzyme aspartate decarboxylase with the purpose of assessing the quantum chemical cluster approach for modeling this important class of decarboxylase enzymes. Active site models ranging in size from 27 to 220 atoms are designed, and the barrier and reaction energy of this step are evaluated. To model the enzyme surrounding, homogeneous polarizable medium techniques are used with several dielectric constants. The main conclusion is that when the active site model reaches a certain size, the solvation effects from the surroundings saturate. Similar results have previously been obtained from systematic studies of other classes of enzymes, suggesting that they are of a quite general nature.

Substituted pyrazoles as novel sEH antagonist: investigation of key binding interactions within the catalytic domain

A novel series of pyrazole sEH inhibitors is reported. Lead optimization efforts to replace the aniline core are also described. In particular, 2-pyridine, 3-pyridine and pyridazine analogs are potent sEH inhibitors with favorable CYP3A4 inhibitory and microsomal stability profiles.

Quantum chemical modeling of enzymatic reactions: the case of histone lysine methyltransferase

Quantum chemical cluster models of enzyme active sites are today an important and powerful tool in the study of various aspects of enzymatic reactivity. This methodology has been applied to a wide spectrum of reactions and many important mechanistic problems have been solved. Herein, we report a systematic study of the reaction mechanism of the histone lysine methyltransferase (HKMT) SET7/9 enzyme, which catalyzes the methylation of the N-terminal histone tail of the chromatin structure. In this study, HKMT SET7/9 serves as a representative case to examine the modeling approach for the important class of methyl transfer enzymes. Active site models of different sizes are used to evaluate the methodology. In particular, the dependence of the calculated energies on the model size, the influence of the dielectric medium, and the particular choice of the dielectric constant are discussed. In addition, we examine the validity of some technical aspects, such as geometry optimization in solvent or with a large basis set, and the use of different density functional methods.

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.

Biosynthesis of ethers: unusual or common natural events?

Ether bonds are found in a wide variety of natural products--mainly secondary metabolites--including lipids, oxiranes, terpenoids, flavonoids, polyketides, and carbohydrate derivatives, to name some representative examples. To furnish such a biodiversity of structures, a large number of different enzymes are involved in several different biosynthetic pathways. Depending on the compound and on the (micro) environment in which the reaction is performed, ethers are produced by very different (enzymatic) reactions, thus providing an impressive display of how Nature has combined evolution and thermodynamics to be able to produce a vast number of compounds. In addition, many of these compounds possess different biological activities of pharmacological interest. Moreover, some of these ethers (i.e., epoxides) have high chemical reactivity, and can be useful starting materials for further synthetic processes. This review aims to provide an overview of the different strategies that are found in Nature for the formation of these "bioethers". Both fundamental and practical insights of the biosynthetic processes will be discussed.

Reaction mechanism of cis-3-chloroacrylic acid dehalogenase: a theoretical study

The reaction mechanism of cis-3-chloroacrylic acid dehalogenase (cis-CaaD) is studied using the B3LYP density functional theory method. This enzyme catalyzes the hydrolytic dehalogenation of cis-3-chloroacrylic acid to yield malonate semialdehyde and HCl. The uncatalyzed reaction is first considered, and excellent agreement is found between the calculated barrier and the measured rate constant. The enzymatic reaction is then studied with an active site model consisting of 159 atoms. The results suggest an alternative mechanism for cis-CaaD catalysis and different roles for some active site residues in this mechanism.

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.

Stereochemical preference of yeast epoxide hydrolase for the O-axial C3 epimers of 1-oxaspiro[2.5]octanes

The 1-oxaspiro[2.5]octane moiety is a common motif in many biologically active spiroepoxide compounds. Stereochemistry plays an important role in the action of these spiroepoxides, since the O-axial C3 epimers are predominantly responsible for biological activity. In view of this, the reactivity of the yeast epoxide hydrolase (YEH) from Rhodotorula glutinis towards both O-axial and O-equatorial C3 epimers of various 1-oxaspiro[2.5]octanes was investigated. O-axial C3 Epimers were hydrolyzed faster than the O-equatorial C3 epimers. The stereochemical preference was greatly dependent on the type of substitution on the cyclohexane ring. The preference of YEH for O-axial C3 epimers, found throughout this study, illustrates the effectiveness of YEH in enzymatic detoxification of spiroepoxides.