QM/MM studies of the enzyme-catalyzed dechlorination of 4-chlorobenzoyl-CoA provide insight into reaction energetics

Microbial degradation of halogenated aromatics: molecular mechanisms and enzymatic reactions

Halogenated aromatics are used widely in various industrial, agricultural and household applications. However, due to their stability, most of these compounds persist for a long time, leading to accumulation in the environment. Biological degradation of halogenated aromatics provides sustainable, low-cost and environmentally friendly technologies for removing these toxicants from the environment. This minireview discusses the molecular mechanisms of the enzymatic reactions for degrading halogenated aromatics which naturally occur in various microorganisms. In general, the biodegradation process (especially for aerobic degradation) can be divided into three main steps: upper, middle and lower metabolic pathways which successively convert the toxic halogenated aromatics to common metabolites in cells. The most difficult step in the degradation of halogenated aromatics is the dehalogenation step in the middle pathway. Although a variety of enzymes are involved in the degradation of halogenated aromatics, these various pathways all share the common feature of eventually generating metabolites for utilizing in the energy-producing metabolic pathways in cells. An in-depth understanding of how microbes employ various enzymes in biodegradation can lead to the development of new biotechnologies via enzyme/cell/metabolic engineering or synthetic biology for sustainable biodegradation processes.

Catalytic mechanism of phenylacetone monooxygenases for non-native linear substrates

Phenylacetone monooxygenase (PAMO) is the most stable and thermo-tolerant member of the Baeyer-Villiger monooxygenase family, and therefore it is an ideal candidate for the synthesis of industrially relevant compounds. However, its limited substrate scope has largely limited its industrial applications. In the present work, we provide, for the first time, the catalytic mechanism of PAMO for the native substrate phenylacetone as well as for a linear non-native substrate 2-octanone, using molecular dynamics simulations, quantum mechanics and quantum mechanics/molecular mechanics calculations. We provide a theoretical basis for the preference of the enzyme for the native aromatic substrate over non-native linear substrates. Our study provides fundamental atomic-level insights that can be employed in the rational engineering of PAMO for wide applications in industrial biocatalysis, in particular, in the biotransformation of long-chain aliphatic oils into potential biodiesels.

Catalytic promiscuity of the non-native FPP substrate in the TEAS enzyme: non-negligible flexibility of the carbocation intermediate

By QM(DFT)/MM MD simulations, it has been revealed that the non-native substrate catalytic promiscuity of TEAS (one of the sesquiterpene cyclases) is mostly attributable to its notable conformational flexibility of the branching intermediate bisabolyl cation.

Molecular perspectives and recent advances in microbial remediation of persistent organic pollutants

Nutrition and pollution stress stimulate genetic adaptation in microorganisms and assist in evolution of diverse metabolic pathways for their survival on several complex organic compounds. Persistent organic pollutants (POPs) are highly lipophilic in nature and cause adverse effects to the environment and human health by biomagnification through the food chain. Diverse microorganisms, harboring numerous plasmids and catabolic genes, acclimatize to these environmentally unfavorable conditions by gene duplication, mutational drift, hypermutation, and recombination. Genetic aspects of some major POP catabolic genes such as biphenyl dioxygenase (bph), DDT 2,3-dioxygenase, and angular dioxygenase assist in degradation of biphenyl, organochlorine pesticides, and dioxins/furans, respectively. Microbial metagenome constitutes the largest genetic reservoir with miscellaneous enzymatic activities implicated in degradation. To tap the metabolic potential of microorganisms, recent techniques like sequence and function-based screening and substrate-induced gene expression are proficient in tracing out novel catabolic genes from the entire metagenome for utilization in enhanced biodegradation. The major endeavor of today's scientific world is to characterize the exact genetic mechanisms of microbes for bioremediation of these toxic compounds by excavating into the uncultured plethora. This review entails the effect of POPs on the environment and involvement of microbial catabolic genes for their removal with the advanced techniques of bioremediation.

Mechanistic insights into cobalt(ii/iii)-catalyzed C-H oxidation: a combined theoretical and experimental study

Cobalt-mediated C-H functionalization has been the subject of extensive interest in synthetic chemistry, but the mechanisms of many of these reactions (such as the cobalt-catalyzed C-H oxidation) are poorly understood. In this paper, possible mechanisms including single electron transfer (SET) and the concerted metalation-deprotonation (CMD) pathways of the CoII/CoIII-catalyzed alkoxylation of C(sp2)-H bonds have been investigated for the first time using the DFT method. CoII(OAc)2 has been employed as an efficient catalyst in our previous experimental study, but the calculated results unexpectedly indicated that the intermolecular SET pathway with CoIII as the actual catalyst might be the most favorable pathway. To support this theoretical prediction, we have explored a series of Cp*CoIII(CO)I2 catalyzed C(sp2)-H bond alkoxylations, extending the application of cobalt-catalyzed functionalization of C-H bonds. Furthermore, kinetic isotope effect (KIE) data, electron paramagnetic resonance (EPR) data, and TEMPO inhibition experiments also support the SET mechanism in both the Co-catalyzed alkoxylation reactions. Thus, this work should support an understanding of the possible mechanisms of the CoII/CoIII-catalyzed C(sp2)-H functionalization, and also provide an example of the rational design of novel catalytic reactions guided by theoretical calculations.

Investigating the role of a backbone to substrate hydrogen bond in OMP decarboxylase using a site-specific amide to ester substitution

Hydrogen bonds between backbone amide groups of enzymes and their substrates are often observed, but their importance in substrate binding and/or catalysis is not easy to investigate experimentally. We describe the generation and kinetic characterization of a backbone amide to ester substitution in the orotidine 5'-monophosphate (OMP) decarboxylase from Methanobacter thermoautotrophicum (MtOMPDC) to determine the importance of a backbone amide-substrate hydrogen bond. The MtOMPDC-catalyzed reaction is characterized by a rate enhancement (∼10(17)) that is among the largest for enzyme-catalyzed reactions. The reaction proceeds through a vinyl anion intermediate that may be stabilized by hydrogen bonding interaction between the backbone amide of a conserved active site serine residue (Ser-127) and oxygen (O4) of the pyrimidine moiety and/or electrostatic interactions with the conserved general acidic lysine (Lys-72). In vitro translation in conjunction with amber suppression using an orthogonal amber tRNA charged with L-glycerate ((HO)S) was used to generate the ester backbone substitution (S127(HO)S). With 5-fluoro OMP (FOMP) as substrate, the amide to ester substitution increased the value of Km by ∼1.5-fold and decreased the value of kcat by ∼50-fold. We conclude that (i) the hydrogen bond between the backbone amide of Ser-127 and O4 of the pyrimidine moiety contributes a modest factor (∼10(2)) to the 10(17) rate enhancement and (ii) the stabilization of the anionic intermediate is accomplished by electrostatic interactions, including its proximity of Lys-72. These conclusions are in good agreement with predictions obtained from hybrid quantum mechanical/molecular mechanical calculations.

Dehalogenation of aromatics by nucleophilic aromatic substitution

Nucleophilic aromatic substitution has been implicated as a mechanism for both the biotic and abiotic hydrodehalogenation of aromatics. Two mechanisms for the aqueous dehalogenation of aromatics involving nucleophilic aromatic substitution with hydride as a nucleophile are investigated using a validated density functional and continuum solvation protocol. For chlorinated and brominated aromatics, nucleophilic addition ortho to carbon-halogen bonds via an anionic intermediate is predicted to be the preferred mechanism in the majority of cases, while concerted substitution is predicted to be preferred for most fluorinated aromatics. Nucleophilic aromatic substitution reactions with the hydroxide and hydrosulfide anions as nucleophiles are also investigated and compared.

A pre-steady state kinetic analysis of the αY60W mutant of trans-3-chloroacrylic acid dehalogenase: implications for the mechanism of the wild-type enzyme

The bacterial degradation of the nematicide 1,3-dichloropropene, an isomeric mixture, requires the action of trans- and cis-3-chloroacrylic acid dehalogenase (CaaD and cis-CaaD, respectively). Both enzymes are tautomerase superfamily members and share a core catalytic mechanism for the hydrolytic dehalogenation of the respective isomer of 3-haloacrylate. The observation that cis-CaaD requires two additional residues raises the question of how CaaD conducts a comparable reaction with fewer catalytic residues. As part of an effort to determine the basis for the apparently simpler CaaD-catalyzed reaction, the kinetic mechanism was determined by stopped-flow and chemical-quench techniques using a fluorescent mutant form of the enzyme, αY60W-CaaD, and trans-3-bromoacrylate as the substrate. The data from these experiments as well as bromide inhibition studies are best accommodated by a six-step model that provides individual rate constants for substrate binding, chemistry, and a proposed conformational change occurring after chemistry followed by release of malonate semialdehyde and bromide. The conformational change and product release rates are comparable, and together they limit the rate of turnover. The kinetic analysis and modeling studies validate the αY60W-CaaD mutant as an accurate reporter of active site events during the course of the enzyme-catalyzed reaction. The kinetic mechanism for the αY60W-CaaD-catalyzed reaction is comparable to that obtained for the cis-CaaD-catalyzed reaction. The kinetic model and the validated αY60W-CaaD mutant set the stage for an analysis of active site mutants to explore the contributions of individual catalytic residues and the basis for the simplicity of the reaction.

Hydrogen-bond stabilization in oxyanion holes: grand jeté to three dimensions

We recently reported crystallographic evidence that the hydrogen bonds which can stabilize oxygen-centered negative charge within enzyme oxyanion holes are rarely found in the place they should be expected on the basis of the analysis of small-molecule crystal structures. We investigated this phenomenon using calculations on simplified active site models. A recent paper suggested that several aspects of the analysis required further exploration. In this paper we: (i) review the results of our crystallographic study; (ii) report molecular dynamics studies which investigate the effect of protein movement; (iii) report ONIOM calculations which trace the reaction coordinate for an oxyanion hole reaction in the presence of a complete enzyme active site. These results show that the limitations of gas phase calculations on simplified models do not invalidate our comparison of competing active site geometries. These new results reaffirm the conclusion that oxyanion holes are not usually stabilized by planar arrangements of H-bonds, and that this sub-optimal transition state stabilization leads to better overall catalysis.

Rate-determining factors in nucleophilic aromatic substitution reactions

Quantum chemical calculations (OPBE/6-311++G(d,p)) have been performed to uncover the electronic factors that govern reactivity in the prototypical S(N)Ar reaction. It was found that intrinsic nucleophilicity--expressed as the critical energy (the energy required for forming the Meisenheimer structure Ph(X)(2)(-)) in the identity substitution reaction X(-) + PhX --> X(-) + PhX (Ph = phenyl)--shows the following approximate trend: NH(2)(-) approximately OH(-) approximately F(-) >> PH(2)(-) approximately SH(-) approximately Cl(-) > AsH(2)(-) approximately SeH(-) approximately Br(-). The periodic trends are discussed in terms of molecular properties (proton affinity of X(-) expressing Lewis basicity of the nucleophile and C(1s) orbital energy expressing Lewis acidity of the substrate) based on a dative bonding model. Furthermore, the stepwise progress of the reactions and the critical structures are analyzed applying energy decomposition analysis. Increased stability, and thereby increased intrinsic nucleophilicity, correlates with decreasing aromatic character of the Meisenheimer structure. This apparent contradiction is explained in consistency with the other observations using the same model.

Pre-steady-state kinetic analysis of cis-3-chloroacrylic acid dehalogenase: analysis and implications

Isomer-specific 3-chloroacrylic acid dehalogenases catalyze the hydrolytic dehalogenation of the cis- and trans-isomers of 3-chloroacrylate to yield malonate semialdehyde. These reactions represent key steps in the degradation of the nematocide, 1,3-dichloropropene. The kinetic mechanism of cis-3-chloroacrylic acid dehalogenase (cis-CaaD) has now been examined using stopped-flow and chemical-quench techniques. Stopped-flow analysis of the reaction, following the fluorescence of an active site tryptophan, is consistent with a minimal three-step model involving substrate binding, chemistry, and product release. Chemical-quench experiments show burst kinetics, indicating that product release is at least partially rate limiting. Global fitting of all of the kinetic results by simulation is best accommodated by a four-step mechanism. In the final kinetic model, the enzyme binds substrate with an immediate isomerization to an alternate fluorescent form and chemistry occurs, followed by the ordered release of two products, with the release of the first product as the rate-limiting step. Bromide ion is a competitive inhibitor of the reaction indicating that it binds to the free enzyme rather than to the enzyme with one product still bound. This observation suggests that malonate semialdehyde is the first product released by the enzyme (rate limiting), followed by halide. A comparison of the unliganded cis-CaaD crystal structure with that of an inactivated cis-CaaD where the prolyl nitrogen of Pro-1 is covalently attached to (R)-2-hydroxypropanoate provides a possible explanation for the isomerization step. The structure of the covalently modified enzyme shows that a seven-residue loop comprised of residues 32-38 is closed down on the active site cavity where the backbone amides of two residues (Phe-37 and Leu-38) interact with the carboxylate group of the adduct. In the unliganded form, the same loop points away from the active site cavity. Similarly, substrate binding may cause this loop to close down on the active site and sequester the reaction from the external environment.

Quantum mechanical/molecular mechanical and density functional theory studies of a prototypical zinc peptidase (carboxypeptidase A) suggest a general acid-general base mechanism

Carboxypeptidase A is a zinc-containing enzyme that cleaves the C-terminal residue in a polypeptide substrate. Despite much experimental work, there is still a significant controversy concerning its catalytic mechanism. In this study, the carboxypeptidase A-catalyzed hydrolysis of the hippuryl-L-Phe molecule (k(cat) = 17.7 +/- 0.7 s(-1)) is investigated using both density functional theory and a hybrid quantum mechanical/molecular mechanical approach. The enzymatic reaction was found to proceed via a promoted-water pathway with Glu270 serving as the general base and general acid. Free-energy calculations indicate that the first nucleophilic addition step is rate-limiting, with a barrier of 17.9 kcal/mol. Besides activating the zinc-bound water nucleophile, the zinc cofactor also serves as an electrophilic catalyst that stabilizes the substrate carbonyl oxygen during the formation of the tetrahedral intermediate. In the Michaelis complex, Arg127, rather than Zn(II), is responsible for the polarization of the substrate carbonyl and it also serves as the oxyanion hole. As a result, its mutation leads to a higher free-energy barrier, in agreement with experimental observations.

Aromatic C-F activation by complexes containing the {Pt2S2} core via nucleophilic substitution: a combined experimental and theoretical study

The C-F bond activation of perfluorobenzene and perfluoropyridine have been achieved by means of the complex [Pt(2)(mu-S)(2)(dppp)(2)], where dppp denotes 1,3-bis(diphenylphosphino)propane. The reaction with the first substrate requires a long time (five days) and high temperature (reflux in toluene) to yield [Pt(o-S(2)C(6)F(4))(dppp)] and [Pt(3)(mu(3)-S)(2)(dppp)(3)]F(2), and involves replacement of two fluorides in the ortho position. In contrast, the reaction with perfluoropyridine is much faster (15 min at 0 degrees C) yielding [Pt(2)(mu-S){mu-(p-SC(5)F(4)N)}(dppp)(2)]F, which implies the C-F activation in the para position with respect the pyridine nitrogen. The mechanism of both reactions has been studied computationally and the geometries of the transition states are consistent with an S(N)Ar mechanism where a sulfido bridging ligand replaces the fluoride anion. The energy barriers corresponding to the first and the second fluoride substitution are 131.7 and 137.1 kJ mol(-1) for perfluorobenzene and 85.9 and 142.7 kJ mol(-1) for perfluoropyridine, respectively. The different energy barrier of the first substitution explains the different experimental conditions required and the various products obtained for these reactions.

Kinetic isotope effects on dehalogenations at an aromatic carbon

In order to interpret the observed isotopic fractionation it is necessaryto understand its relationship with the isotope effect(s) on steps that occur during the conversion of the initial reactant to the final product. We examine this relationship from the biochemical point of view and elaborate on the consequences of the assumptions that it is based on. We illustrate the discrepancies between theoretical and experimental interpretation of kinetic isotope effects on examples of dehalogenation reactions that occur at an aromatic carbon atom. The examples include 4-chlorobenzoyl-CoA dehalogenase-catalyzed conversion of 4-chlorobenzoyl-CoA to 4-hydroxybenzoyl-CoA, dehaloperoxidase-catalyzed conversion of 2,4,6-trichlorophenol to 2,6-dichloroquinone, and spontaneous hydrolysis of atrazine at pH 12. For this latter reaction we have measured the chlorine kinetic isotope effect and estimated its value theoretically at the DFT level of theory. Results of chlorine kinetic isotope effects suggest that the studied dechlorination reactions proceed in a single step with significant weakening of the carbon-chlorine bond in the transition state.

Theoretical study of general base-catalyzed hydrolysis of aryl esters and implications for enzymatic reactions

In this work, the mechanism of general base-catalyzed hydrolysis of aryl esters is investigated in vacuo with density functional theory and in solutions with a polarized continuum model. The hydrolysis is found to proceed via a concerted mechanism featuring simultaneous addition and elimination steps accompanied by proton transfers, consistent with experimental evidence. Reasonable agreement with measured kinetic isotope effects provides additional validation. It is found that solvation substantially lowers the transition state energy, but has a small effect on the reaction exothermicity. An enzyme oxyanion hole, modeled by an ammonia molecule hydrogen bonded to the acyl carbonyl oxygen, is found to stabilize the near-tetrahedral transition state. Implications of these findings for the hydrolysis step of the dehalogenation reaction catalyzed by 4-chlorobenzoyl-CoA dehalogenase are discussed.

Theoretical Study of the Catalytic Mechanism and Metal-Ion Dependence of Peptide Deformylase

The reaction pathway of deformylation catalyzed by E. coli peptide deformylase (PDF) has been investigated by the density functional theory method of PBE1PBE on a small model and by a two-layer ONIOM method on a realistic protein model. The deformylation proceeds in sequential steps involving nucleophilic addition of metal-coordinated water/hydroxide to the carbonyl carbon of the formyl group, proton transfer, and cleavage of the C-N bond. The first step is rate-determining for the deformylation, which occurs through a pentacoordinated metal center. The estimated activation energies with the ONIOM method are about 23.0, 15.0, and 14.9 kcal/mol for Zn-, Ni-, and Fe-PDFs, respectively. These calculated barriers are in close agreement with experimental observations. Our results demonstrate that the preference for metal coordination geometry exerts a significant influence on the catalytic activity of PDFs by affecting the activation of the carbonyl group of the substrate, the deprotonation of the metal-coordinated water, and the stabilization of the transition state. This preference for coordination geometry is mainly determined by the ligand environment and the intrinsic electronic structures of the metal center in the active site of the PDFs.

A computational study on the amine-oxidation mechanism of monoamine oxidase: Insight into the polar nucleophilic mechanism

The proposed polar nucleophilic mechanism of MAO was investigated using quantum chemical calculations employing the semi-empirical PM3 method. In order to mimic the reaction at the enzyme's active site, the reactions between the flavin and the p-substituted benzylamine substrate analogs were modeled. Activation energies and rate constants of all the reactions were calculated and compared with the published experimental data. The results showed that electron-withdrawing groups at the para position of benzylamine increase the reaction rate. A good correlation between the log of the calculated rate constants and the electronic parameter (sigma) of the substituent was obtained. These results agree with the previous kinetic experiments on the effect of p-substituents on the reduction of MAO-A by benzylamine analogs. In addition, the calculated rate constants showed a correlation with the rate of reduction of the flavin in MAO-A. In order to verify the results obtained from the PM3 method single-point B3LYP/6-31G*//PM3 calculations were performed. These results demonstrated a strong reduction in the activation energy for the reaction of benzylamine derivatives having electron-withdrawing substituents, which is in agreement with the PM3 calculations and the previous experimental QSAR study. PM3 and B3LYP/6-31G* energy surfaces were obtained for the overall reaction of benzylamine with flavin. Results suggest that PM3 is a reasonable method for studying this kind of reaction. These theoretical findings support the proposed polar nucleophilic mechanism for MAO-A.

Modelling enzyme reaction mechanisms, specificity and catalysis

Modern modelling methods can now give uniquely detailed understanding of enzyme-catalyzed reactions, including the analysis of mechanisms and the identification of determinants of specificity and catalytic efficiency. A new field of computational enzymology has emerged that has the potential to contribute significantly to structure-based design and to develop predictive models of drug metabolism and, for example, of the effects of genetic polymorphisms. This review outlines important techniques in this area, including quantum-chemical model studies and combined quantum-mechanics and molecular-mechanics (QM/MM) methods. Some recent applications to enzymes of pharmacological interest are also covered, showing the types of problems that can be tackled and the insight they can give.

Electrostatic influence of active-site waters on the nucleophilic aromatic substitution catalyzed by 4-chlorobenzoyl-CoA dehalogenase

The 4-chlorobenzoyl-CoA dehalogenase catalyzes the hydrolytic dechlorination of 4-chlorobenzoyl-CoA via a two-step mechanism, namely nucleophilic aromatic substitution and ester hydrolysis. The mutation of an active-site Histidine residue has been shown to reduce the catalytic activity in both the substitution and subsequent hydrolysis steps. In this communication, we report a quantum mechanical/molecular mechanical simulation of the potential of mean force for the substitution step, which confirms the increased barrier height in the H90Q mutant and provides evidence on the electrostatic influence of two active-site waters on the rate-limiting barrier.

2-Haloacrylate reductase, a novel enzyme of the medium chain dehydrogenase/reductase superfamily that catalyzes the reduction of a carbon-carbon double bond of unsaturated organohalogen compounds

A soil bacterium, Burkholderia sp. WS, grows on 2-chloroacrylate as the sole carbon source. To identify the enzymes metabolizing 2-chloroacrylate, we carried out comparative two-dimensional gel electrophoresis of the proteins from 2-chloroacrylate- and lactate-grown bacterial cells. As a result, we found that a protein named CAA43 was inducibly synthesized when the cells were grown on 2-chloroacrylate. The CAA43 gene was cloned and shown to encode a protein of 333 amino acid residues (M(r) 35,788) that shared a significant sequence similarity with NADPH-dependent quinone oxidoreductase from Escherichia coli (38.2% identity). CAA43 was overproduced in E. coli and purified to homogeneity. The purified protein catalyzed the NADPH-dependent reduction of the carbon-carbon double bond of 2-chloroacrylate to produce (S)-2-chloropropionate, which is probably further metabolized to (R)-lactate by (S)-2-haloacid dehalogenase in Burkholderia sp. WS. NADH did not serve as a reductant. Despite the sequence similarity to quinone oxidoreductases, CAA43 did not act on 1,4-benzoquinone and 1,4-naphthoquinone. 2-Chloroacrylate analogs, such as acrylate and methacrylate, were also inert as the substrates. In contrast, 2-bromoacrylate served as the substrate. Thus, we named this novel enzyme 2-haloacrylate reductase. This study revealed a new pathway for the degradation of unsaturated organohalogen compounds. It is also notable that the enzyme is useful for the production of (S)-2-chloropropionate, which is used for the industrial production of aryloxyphenoxypropionic acid herbicides.