Investigation of substrate activation by 4-chlorobenzoyl-coenzyme A dehalogenase

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.

Structure and denaturation of 4-chlorobenzoyl coenzyme A dehalogenase from Arthrobacter sp. strain TM-1

The secondary structure of the trimeric protein 4-chlorobenzoyl coenzyme A dehalogenase from Arthrobacter sp. strain TM-1, the second of three enzymes involved in the dechlorination of 4-chlorobenzoate to form 4-hydroxybenzoate, has been examined. E(mM) for the enzyme was 12.59. Analysis by circular dichroism spectrometry in the far uv indicated that 4-chlorobenzoyl coenzyme A dehalogenase was composed mostly of alpha-helix (56%) with lesser amounts of random coil (21%), beta-turn (13%) and beta-sheet (9%). These data are in close agreement with a computational prediction of secondary structure from the primary amino acid sequence, which indicated 55.8% alpha-helix, 33.7% random coil and 10.5% beta-sheet; the enzyme is, therefore, similar to the 4-chlorobenzoyl coenzyme A dehalogenase from Pseudomonas sp. CBS-3. The three-dimensional structure, including that of the presumed active site, predicted by computational analysis, is also closely similar to that of the Pseudomonas dehalogenase. Study of the stability and physicochemical properties revealed that at room temperature, the enzyme was stable for 24 h but was completely inactivated by heating to 60 degrees C for 5 min; thereafter by cooling at 1 degrees C min(-1) to 45 degrees C, 20.6% of the activity could be recovered. Mildly acidic (pH 5.2) or alkaline (pH 10.1) conditions caused complete inactivation, but activity was fully recovered on returning the enzyme to pH 7.4. Circular dichroism studies also indicated that secondary structure was little altered by heating to 60 degrees C, or by changing the pH from 7.4 to 6.0 or 9.2. Complete, irreversible destruction of, and maximal decrease in the fluorescence yield of the protein at 330-350 nm were brought about by 4.5 M urea or 1.1 M guanidinium chloride. Evidence was obtained to support the hypothetical three-dimensional model, that residues W140 and W167 are buried in a non-polar environment, whereas W182 appears at or close to the surface of the protein. At least one of the enzymes of the dehalogenase system (the combined 4-chlorobenzoate:CoA ligase, the dehalogenase and 4-hydroxybenzoyl coenzyme A thioesterase) appears to be capable of association with the cell membrane.

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.

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

The conversion of 4-chlorobenzoyl-CoA to 4-hydroxybenzoyl-CoA catalyzed by 4-chlorobenzoyl-CoA dehalogenase is investigated using combined QM/MM approaches. The calculated potential of mean force at the PM3/CHARMM level supports the proposed nucleophilic aromatic substitution mechanism. In particular, a Meisenheimer intermediate was found, stabilized by hydrogen bonds between the benzoyl carbonyl of the ligand and two backbone amide NHs at positions 64 and 114. Mutation of Gly113 to Ala significantly increases the barrier by disrupting the hydrogen bond with the Gly114 backbone. The formation of the Meisenheimer complex is accompanied by significant charge redistribution and structural changes in the substrate benzoyl moiety, consistent with experimental observations. Theoretical results suggest that the reaction rate is limited by the formation of the Meisenheimer complex, rather than by its decomposition. A kinetic model based on the calculated free energy profile is found to be consistent with the experimental time course data.

Stereoselectivity of enoyl-CoA hydratase results from preferential activation of one of two bound substrate conformers

Enoyl-CoA hydratase catalyzes the hydration of trans-2-crotonyl-CoA to 3(S)- and 3(R)-hydroxybutyryl-CoA with a stereoselectivity (3(S)/3(R)) of 400,000 to 1. Importantly, Raman spectroscopy reveals that both the s-cis and s-trans conformers of the substrate analog hexadienoyl-CoA are bound to the enzyme, but that only the s-cis conformer is polarized. This selective polarization is an example of ground state strain, indicating the existence of catalytically relevant ground state destabilization arising from the selective complementarity of the enzyme toward the transition state rather than the ground state. Consequently, the stereoselectivity of the enzyme-catalyzed reaction results from the selective activation of one of two bound substrate conformers rather than from selective binding of a single conformer. These findings have important implications for inhibitor design and the role of ground state interactions in enzyme catalysis.

Divergent evolution of enzymatic function: mechanistically diverse superfamilies and functionally distinct suprafamilies

The protein sequence and structure databases are now sufficiently representative that strategies nature uses to evolve new catalytic functions can be identified. Groups of divergently related enzymes whose members catalyze different reactions but share a common partial reaction, intermediate, or transition state (mechanistically diverse superfamilies) have been discovered, including the enolase, amidohydrolase, thiyl radical, crotonase, vicinal-oxygen-chelate, and Fe-dependent oxidase superfamilies. Other groups of divergently related enzymes whose members catalyze different overall reactions that do not share a common mechanistic strategy (functionally distinct suprafamilies) have also been identified: (a) functionally distinct suprafamilies whose members catalyze successive transformations in the tryptophan and histidine biosynthetic pathways and (b) functionally distinct suprafamilies whose members catalyze different reactions in different metabolic pathways. An understanding of the structural bases for the catalytic diversity observed in super- and suprafamilies may provide the basis for discovering the functions of proteins and enzymes in new genomes as well as provide guidance for in vitro evolution/engineering of new enzymes.

The active site dynamics of 4-chlorobenzoyl-CoA dehalogenase

A molecular dynamics study was performed to compare the differences in the active-site dynamics of the wild-type and W137F mutant enzymes of 4-chlorobenzoyl-CoA dehalogenase. Only in the wild-type simulation are conformations formed between the catalytic Asp-145 and 4-chlorobenzoyl-CoA, which resemble the ab initio calculated gas-phase transition-state geometry. In the W137F simulation, the hydrogen bond formed between His-90 and Asp-145 persisted throughout the simulation, causing the carboxylate of Asp-145 to be distant from the benzoyl ring of 4-chlorobenzoyl-CoA. In both simulations, water molecules were able to diffuse into the active site of the enzymes. The trajectories provide insight into the routes that water may use to get into position for the hydrolysis portion of the dehalogenation reaction. In both simulations, the water molecule entering the active site forms a hydrogen bond with Asp-145.

Chlorine kinetic isotope effects on the haloalkane dehalogenase reaction

We have found chlorine kinetic isotope effects on the dehalogenation catalyzed by haloalkane dehalogenase from Xanthobacter autotrophicus GJ10 to be 1.0045 +/- 0.0004 for 1,2-dichloroethane and 1.0066 +/- 0.0004 for 1-chlorobutane. The latter isotope effect approaches the intrinsic chlorine kinetic isotope effect for the dehalogenation step. The intrinsic isotope effect has been modeled using semiempirical and DFT theory levels using the ONIOM QM/QM scheme. Our results indicate that the dehalogenation step is reversible; the overall irreversibility of the enzyme-catalyzed reaction is brought about by a step following the dehalogenation.

Product catalyzes the deamidation of D145N dehalogenase to produce the wild-type enzyme

Aspartate 145 plays an essential role in the active site of 4-chlorobenzoyl-CoA dehalogenase, forming a transient covalent link at the 4-position of the benzoate during the conversion of the substrate to 4-hydroxybenzoyl-CoA. Replacement of Asp 145 by residues such as alanine or serine results in total inactivation, and stable complexes can be formed with either substrate or product. The Raman spectroscopic characterization of some of the latter is described in the preceding publication (Dong et al.). The present work investigates complexes formed by D145N dehalogenase and substrate or product. Time-resolved absorption and Raman difference spectroscopic data show that these systems evolve rapidly with time. For the substrate complex, initially the absorption and Raman spectra show the signatures of the substrate bound in the active site of the asparagine 145 form of the enzyme but these signatures are accompanied by those for the ionized product. After several minutes these signatures disappear to be replaced with those closely resembling the un-ionized product in the active site of wild-type dehalogenase. Similarly, for the product complex, the absorption and Raman spectra initially show evidence for ionized product in the active site of D145N, but these are rapidly replaced by signatures closely resembling the un-ionized product bound to wild-type enzyme. It is proposed that product bound to the active site of asparagine 145 dehalogenase catalyzes the deamidation of the asparagine side chain to produce the wild-type aspartate 145. For the complexes involving substrate, the asparagine 145 enzyme population contains a small amount of the WT enzyme, formed by spontaneous deamidation, that produces product. In turn, these product molecules catalyze the deamidation of Asn 145 in the major enzyme population. Thus, conversions of substrate to product and of D145N to D145D dehalogenase go on simultaneously. The spontaneous deamidation of asparagine 145 has been characterized by allowing the enzyme to stand at RT in Hepes buffer at pH 7.5. Under these conditions deamidation occurs with a rate constant of 0.0024 h-1. The rate of product-catalyzed deamidation in Hepes buffer at 22 degrees C was measured by stopped-flow kinetics to be 0.024 s-1, 36000 times faster than the spontaneous process. A feature near 1570 cm-1 could be observed in the early Raman spectra of both substrate and product-enzyme complexes. This band is not associated with either substrate or product and is tentatively assigned to an ester-like species formed by the attack of the product's 4-O- group on the carbonyl of asparagine's side chain and the subsequent release of ammonia. A reaction scheme is proposed, incorporating these observations.

Modulating electron density in the bound product, 4-hydroxybenzoyl-CoA, by mutations in 4-chlorobenzoyl-CoA dehalogenase near the 4-hydroxy group

The enzyme 4-chlorobenzoyl-CoA dehalogenase hydrolyzes 4-chlorobenzoyl-CoA (4-CBA-CoA) to 4-hydroxybenzoyl-CoA (4-HBA-CoA). Biochemical and crystallographic studies have identified a critical role for the dehalogenase residue Asp 145 in close proximity to the ligand's 4-hydroxy group in the structure of the product-enzyme complex. In the present study the effects of site selective mutations at Asp 145 on the product complex are explored by Raman spectroscopy. The spectral signatures of the WT-product complex, the large red shift in lambdamax, and the complete reorganization of the benzoyl ring modes in Raman data are absent for the D145E complex. The major spectral perturbations in the WT complex are brought about by strong electron "pull" at the benzoyl carbonyl and electron "push" by the side chain of Asp 145 near the 4-OH group. Acting in concert, these factors polarize the benzoyl's pi-electrons. Since the Raman data show that very strong electron pull occurs at the benzoyl's carbonyl in the D145E complex, it is apparent that the needed electron push near the benzoyl's 4-OH group is missing. Thus, very precise positioning of Asp 145's side chain near the benzoyl's 4-position is needed to bring about the dramatic electron reorganization seen in the WT complex, and this criterion cannot be met by the glutamate side chain with its additional CH2 group. For two other Asp145 mutants D145A and D145S that lack catalytic activity, Raman difference spectroscopic data for product complexes demonstrate the presence of a population of ionized product (i.e., 4-O-) in the active sites. The presence of the ionized phenolate form explains the observation that these complexes have highly red-shifted absorbance maxima with lambdamaxs near 400 nm. For the WT complex only the 4-OH form is seen, ionization being energetically expensive with the presence of the proximal negative charge on the Asp 145 side chain. Semiquantitative estimates of the pKa for the bound product in D145S and D145A indicate that this ionization lies in the pH 6.5-7.0 range. This is approximately 2 pH units below the pKa for the free product. The Raman spectrum of 4-dimethylaminobenzoyl-CoA undergoes major changes upon binding to dehalogenase. The bound form has two features near 1562 and 1529 cm-1 and therefore closely resembles the spectrum of product bound to wild-type enzyme, which underlines the quinonoid nature in these complexes. The use of a newly developed Raman system allowed us to obtain normal (nonresonance) Raman data for the dehalogenase complexes in the 100-300 microM range and heralds an important advance in the application of Raman spectroscopy to dilute solutions of macromolecules.

The crystal structure of dienoyl-CoA isomerase at 1.5 A resolution reveals the importance of aspartate and glutamate sidechains for catalysis

BACKGROUND

The degradation of unsaturated fatty acids is vital to all living organisms. Certain unsaturated fatty acids must be catabolized via a pathway auxiliary to the main beta-oxidation pathway. Dienoyl-coenzyme A (dienoyl-CoA) isomerase catalyzes one step of this auxiliary pathway, the isomerization of 3-trans,5-cis-dienoyl-CoA to 2-trans,4-trans-dienoyl-CoA, and is imported into both mitochondria and peroxisomes. Dienoyl-CoA isomerase belongs to a family of CoA-binding proteins that share the enoyl-CoA hydratase/isomerase sequence motif.

RESULTS

The crystal structure of rat dienoyl-CoA isomerase has been determined at 1.5 A resolution. The fold closely resembles that of enoyl-CoA hydratase and 4-chlorobenzoyl-CoA dehalogenase. Dienoyl-CoA isomerase forms hexamers made up of two trimers. The structure contains a well ordered peroxisomal targeting signal type-1 which is mostly buried in the inter-trimer space. The active-site pocket is deeply buried and entirely hydrophobic, with the exception of the acidic residues Asp176, Glu196 and Asp204. Site-directed mutagenesis of Asp204 revealed that this residue is essential for catalysis. In a molecular modeling simulation, a molecule of 3-trans,5-cis-octadienoyl-CoA was docked into the active site.

CONCLUSIONS

The structural data, supported by the mutagenesis data, suggest a reaction mechanism where Glu196 acts as a proton acceptor and Asp204 acts as a proton donor. Asp176 is paired with Glu196 and is important for optimizing the catalytic proton transfer properties of Glu196. In the predicted mode of substrate binding, an oxyanion hole stabilizes the transition state by binding the thioester oxygen. The presence of a buried peroxisomal targeting signal suggests that dienoyl-CoA isomerase is prevented from reaching its hexameric structure in the cytosol.

Three-dimensional structure of L-2-haloacid dehalogenase from Xanthobacter autotrophicus GJ10 complexed with the substrate-analogue formate

The L-2-haloacid dehalogenase from the 1,2-dichloroethane degrading bacterium Xanthobacter autotrophicus GJ10 catalyzes the hydrolytic dehalogenation of small L-2-haloalkanoic acids to yield the corresponding D-2-hydroxyalkanoic acids. Its crystal structure was solved by the method of multiple isomorphous replacement with incorporation of anomalous scattering information and solvent flattening, and was refined at 1.95-A resolution to an R factor of 21.3%. The three-dimensional structure is similar to that of the homologous L-2-haloacid dehalogenase from Pseudomonas sp. YL (1), but the X. autotrophicus enzyme has an extra dimerization domain, an active site cavity that is completely shielded from the solvent, and a different orientation of several catalytically important amino acid residues. Moreover, under the conditions used, a formate ion is bound in the active site. The position of this substrate-analogue provides valuable information on the reaction mechanism and explains the limited substrate specificity of the Xanthobacter L-2-haloacid dehalogenase.

Raman study of the polarizing forces promoting catalysis in 4-chlorobenzoate-CoA dehalogenase

The enzyme 4-chlorobenzoate-CoA dehalogenase catalyzes the hydrolysis of 4-chlorobenzoate-CoA (4-CBA-CoA) to 4-hydroxybenzoyl-CoA (4-HBA-CoA). In order to facilitate electrophilic catalysis, the dehalogenase utilizes a strong polarizing interaction between the active site residues and the benzoyl portion of the substrate [Taylor, K. L., et al. (1995) Biochemistry 34, 13881]. As a result of this interaction, the normal modes of the benzoyl moiety of the bound 4-HBA-CoA undergo a drastic rearrangement as shown by Raman spectroscopy. Here, we present Raman difference spectroscopic data on the product-enzyme complex where the product's benzoyl carbonyl is labeled with 18O (C=18O) or 13C (13C=O) or where the 4-OH group is labeled with 18O. The data demonstrate that the carbonyl group participates in the most intense normal modes occurring in the Raman spectrum in the 1520-1560 cm-1 region. The substrate analog 4-methylbenzoate-CoA (4-MeBA-CoA) has also been characterized by Raman difference spectroscopy in its free form and bound to the dehalogenase. Upon binding, the 4-MeBA-CoA shows evidence of polarization within the delocalized pi-electrons, but to a lesser extent compared to that seen for the product. The use of 4-MeBA-CoA labeled with 18O at the carbonyl enables us to estimate the degree of electron polarization within the C=O group of the bound 4-MeBA-CoA. The C=O stretching frequency occurs near 1663 cm-1 in non-hydrogen bonding solvents such as CCl4, near 1650 cm-1 in aqueous solution, and near 1610 cm-1 in the active site of dehalogenase. From model studies, we can estimate that in the active site the carbonyl group behaves as though it is being polarized by hydrogen bonds approximately 57 kJ mol-1 in strength. Major contributions to this polarization come from hydrogen bonds from the peptide NHs of Gly114 and Phe64. However, an additional contribution, which may account for up to half of the observed shift in nuC=O, originates in the electrostatic field due to the alpha-helix dipole from residues 121-114. The helix which terminates at Gly114, near the C=O group of the bound benzoyl, provides a dipolar electrostatic component which contributes to the polarization of the C=O bond and to the polarization of the entire benzoyl moiety. The effect of both the helix dipole and the hydrogen bonds on the C=O is a "pull" of electrons onto the carbonyl oxygen, which, in turn, polarizes the electron distribution within the benzoyl pi-electron system. The ability of these two factors to polarize the electrons within the benzoyl moiety is increased by the environment about the benzoyl ring; it is surrounded by hydrophobic residues which provide a low-dielectric constant microenvironment. Electron polarization promotes catalysis by reducing electron density at the C4 position of the benzoyl ring, thereby assisting attack by the side chain of Asp145. An FTIR study on the model compound 4-methylbenzoyl S-ethyl thioester, binding to a number of hydrogen bonding donors in CCl4, is described and is used to relate the observed shift of the C=O stretching mode of 4-MeBA-CoA in the active site to the hydrogen bonding strength value. Since the shift of the C=O frequency upon binding is due to hydrogen bonding and helix dipole effects, we refer to this bonding strength as the effective hydrogen bonding strength.