The role of halogen substituents and substrate pKa in defining the substrate specificity of 2,6-dichlorohydroquinone 1,2-dioxygenase (PcpA).

2,6-Dichlorohydroquinone 1,2-dioxygenase (PcpA) is a non-heme Fe(II) enzyme that is specific for ortho-dihalohydroquinones. Here we deconvolute the role of halogen polarizability vs. substrate pKa in defining this specificity, and show how substrate binding compares to the structurally homologous catechol extradiol dioxygenases. The substrates 2,6-dichloro- and 2,6-dibromohydroquinone (polarizable halogens, pKa1 = 7.3), 2,6-difluorohydroquinone (nonpolarizable halogens, pKa1 = 7.5), and 2-chloro-6-methylhydroquinone (polarizable halogen, pKa1 = 9.0) were examined through spectrophotometric titrations and steady-state kinetics. The results show that binding of the substrates to the enzyme decreased [Formula: see text] by about 0.5, except for 2,6-difluorohydroquinone, which showed no change. Additionally, the Kd values of 2,6-dichloro- and 2,6-dibromohydroquinone are about equal to their respective [Formula: see text]. For comparison, with catechol 2,3-dioxygenase (XylE), the substrates 4-methyl- and 3-bromocatechol are bound to the enzyme exclusively in the monoanion form over a wide pH range, indicating a [Formula: see text] of at least - 2.9 and - 1.2, respectively. The steady-state kinetic studies showed that 2,6-difluorohydroquinone is a poor substrate, with [Formula: see text] approximately 40-fold lower and [Formula: see text] 20-fold higher than 2,6-dichlorohydroquinone, despite its similar pKa1. Likewise, the pH dependence of [Formula: see text] for 2-chloro-6-methylhydroquinone is nearly identical to that of 2,6-dichlorohydroquinone, despite its very different pKa1. These results show that (1) it is clearly the halogen polarizability and not the lower substrate pKa that determines the substrate specificity of PcpA, and (2) that PcpA, unlike the catechol extradiol dioxygenases, lacks an active site base that assists with substrate deprotonation, highlighting a key functional difference in what are otherwise similar active sites that defines their different reactivity.

Structural and spectroscopic characterization of iron(II), cobalt(II), and nickel(II) ortho-dihalophenolate complexes: insights into metal-halogen secondary bonding.

Metal complexes incorporating the tris(3,5-diphenylpyrazolyl)borate ligand (Tp(Ph2)) and ortho-dihalophenolates were synthesized and characterized in order to explore metal-halogen secondary bonding in biorelevant model complexes. The complexes Tp(Ph2)ML were synthesized and structurally characterized, where M was Fe(II), Co(II), or Ni(II) and L was either 2,6-dichloro- or 2,6-dibromophenolate. All six complexes exhibited metal-halogen secondary bonds in the solid state, with distances ranging from 2.56 Å for the Tp(Ph2)Ni(2,6-dichlorophenolate) complex to 2.88 Å for the Tp(Ph2)Fe(2,6-dibromophenolate) complex. Variable temperature NMR spectra of the Tp(Ph2)Co(2,6-dichlorophenolate) and Tp(Ph2)Ni(2,6-dichlorophenolate) complexes showed that rotation of the phenolate, which requires loss of the secondary bond, has an activation barrier of ~30 and ~37 kJ/mol, respectively. Density functional theory calculations support the presence of a barrier for disruption of the metal-halogen interaction during rotation of the phenolate. On the other hand, calculations using the spectroscopically calibrated angular overlap method suggest essentially no contribution of the halogen to the ligand-field splitting. Overall, these results provide the first quantitative measure of the strength of a metal-halogen secondary bond and demonstrate that it is a weak noncovalent interaction comparable in strength to a hydrogen bond. These results provide insight into the origin of the specificity of the enzyme 2,6-dichlorohydroquinone 1,2-dioxygenase (PcpA), which is specific for ortho-dihalohydroquinone substrates and phenol inhibitors.

Substrate specificity of Sphingobium chlorophenolicum 2,6-dichlorohydroquinone 1,2-dioxygenase.

PcpA is an aromatic ring-cleaving dioxygenase that is homologous to the well-characterized Fe(II)-dependent catechol extradiol dioxygenases. This enzyme catalyzes the oxidative cleavage of 2,6-dichlorohydroquinone in the catabolism of pentachlorophenol by Sphingobium chlorophenolicum ATCC 39723. (1)H NMR and steady-state kinetics were used to determine the regiospecificity of ring cleavage and the substrate specificity of the enzyme. PcpA exhibits a high degree of substrate specificity for 2,6-disubstituted hydroquinones, with halogens greatly preferred at those positions. Notably, the k(cat)(app)/K(mA)(app) of 2,6-dichlorohydroquinone is ~40-fold higher than that of 2,6-dimethylhydroquinone. The asymmetric substrate 2-chloro-6-methylhydroquinone yields a mixture of 1,2- and 1,6-cleavage products. These two modes of cleavage have different K(mO(2))(app) values (21 and 260 µM, respectively), consistent with a mechanism in which the substrate binds in two catalytically productive orientations. In contrast, monosubstituted hydroquinones show a limited amount of ring cleavage but rapidly inactivate the enzyme in an O(2)-dependent fashion, suggesting that oxidation of the Fe(II) may be the cause. Potent inhibitors of PcpA include ortho-disubstituted phenols and 3-bromocatechol. 2,6-Dibromophenol is the strongest competitive inhibitor, consistent with PcpA's substrate specificity. Several factors that could yield this specificity for halogen substituents are discussed. Interestingly, 3-bromocatechol also inactivates the enzyme, while 2,6-dihalophenols do not, indicating a requirement for two hydroxyl groups for ring cleavage and for enzyme inactivation. These results provide mechanistic insights into the hydroquinone dioxygenases.

Solution and structural characterization of iron(II) complexes with ortho-halogenated phenolates: insights into potential substrate binding modes in hydroquinone dioxygenases.

The new ligand cis,cis-1,3,5-tris-(E)-(tolylideneimino)cyclohexane (TACH-o-tolyl) forms a 1:1 complex with iron(II). Addition of substituted phenolates forms 1:1:1 ligand:iron:phenolate complexes, which have been characterized both in the solid state and in solution. There is complete binding of the phenolate to the complex only when there are ortho-halogens on the phenolate. The tertiary complexes with ortho-halo-substituted phenolates exhibit short Fe-halogen distances, and the complex containing a non-coordinating but similarly sized ortho-methyl phenolate has a significantly different conformation and coordination geometry. Therefore, it is likely that the metal-halogen interaction stabilizes the complexes. The iron(II)-halogen interaction in these complexes may explain the substrate specificity of PcpA and LinE, enzymes that preferentially bind phenols and hydroquinones containing halogen substituents in ortho positions.

Determination of the active site of Sphingobium chlorophenolicum 2,6-dichlorohydroquinone dioxygenase (PcpA).

2,6-Dichlorohydroquinone 1,2-dioxygenase (PcpA) from Sphingobium chlorophenolicum ATCC 39723 is a member of a class of Fe(II)-containing hydroquinone dioxygenases that is involved in the mineralization of the pollutant pentachlorophenol. This enzyme has not been extensively characterized, despite its interesting ring-cleaving activity and use of Fe(II), which are reminiscent of the well-known extradiol catechol dioxygenases. On the basis of limited sequence homology to the extradiol catechol dioxygenases, the residues ligating the Fe(II) center were originally proposed to be H159, H227, and E276 (Xu et al. in Biochemistry 38:7659-7669, 1999). However, PcpA has higher sequence homology to a newly reported, crystallographically characterized zinc metalloenzyme that has a similar predicted fold. We generated a homology model of the structure of PcpA based upon the structure of this zinc metalloenzyme. The homology model predicts that the tertiary structure of PcpA differs significantly from that of the extradiol dioxygenases, and that the residues ligating the Fe(II) are H11, H227, and E276. This structural model was tested by mutating each of H11, H159, H227, and E276 to alanine. An additional residue that is predicted to lie near the active site and is conserved among PcpA, its closest homologues, and the extradiol dioxygenases, Y266, was mutated to phenylalanine. Of these mutants, only H159A retained significant activity, thus confirming the active-site location predicted by the homology-based structural model. The model provides an important basis for understanding the origin of the unique function of PcpA.

Hyperfine-shifted (13)C and (15)N NMR signals from Clostridium pasteurianum rubredoxin: extensive assignments and quantum chemical verification.

Stable isotope-labeling methods, coupled with novel techniques for detecting fast-relaxing NMR signals, now permit detailed investigations of paramagnetic centers of metalloproteins. We have utilized these advances to carry out comprehensive assignments of the hyperfine-shifted (13)C and (15)N signals of the rubredoxin from Clostridium pasteurianum (CpRd) in both its oxidized and reduced states. We used residue-specific labeling (by chemical synthesis) and residue-type-selective labeling (by biosynthesis) to assign signals detected by one-dimensional (15)N NMR spectroscopy, to nitrogen atoms near the iron center. We refined and extended these (15)N assignments to the adjacent carbonyl carbons by means of one-dimensional (13)C[(15)N] decoupling difference experiments. We collected paramagnetic-optimized SuperWEFT (13)C[(13)C] constant time COSY (SW-CT-COSY) data to complete the assignment of (13)C signals of reduced CpRd. By following these (13)C signals as the protein was gradually oxidized, we transferred these assignments to carbons in the oxidized state. We have compared these assignments with hyperfine chemical shifts calculated from available X-ray structures of CpRd in its oxidized and reduced forms. The results allow the evaluation of the X-ray structural models as representative of the solution structure of the protein, and they provide a framework for future investigation of the active site of this protein. The methods developed here should be applicable to other proteins that contain a paramagnetic center with high spin and slow electron exchange.

Changes in hydrogen-bond strengths explain reduction potentials in 10 rubredoxin variants.

The rubredoxin from Clostridium pasteurianum (CpRd) provides an excellent system for investigating how the protein sequence modulates the reduction potential of the active site in an iron-sulfur protein. (15)N NMR spectroscopy has allowed us to determine with unprecedented accuracy the strengths of all six key hydrogen bonds between protein backbone amides and the sulfur atoms of the four cysteine residues that ligate the iron in the oxidized (Fe(III)) and reduced (Fe(II)) forms of wild-type CpRd and nine mutants (V44G, V44A, V44I, V44L, V8G, V8A, V8I, V8L, and V8G/V44G). The length (or strength) of each hydrogen bond was inferred from the magnitude of electron spin delocalized across the hydrogen bond from the iron atom onto the nitrogen. The aggregate lengths of these six hydrogen bonds are shorter in both oxidation states in variants with higher reduction potential than in those with lower reduction potential. Differences in aggregate hydrogen bonding upon reduction correlate linearly with the published reduction potentials for the 10 CpRd variants, which span 126 mV. Sequence effects on the reduction potential can be explained fully by their influence on hydrogen-bond strengths.

Paramagnetic NMR spectroscopy and density functional calculations in the analysis of the geometric and electronic structures of iron-sulfur proteins.

Paramagnetic NMR spectroscopy has been underutilized in the study of metalloproteins. One difficulty of the technique is that paramagnetic relaxation broadens signals from nuclei near paramagnetic centers. In systems with low electronic relaxation rates, this makes such signals difficult to observe or impossible to assign by traditional methods. We show how the challenges of detecting and assigning signals from nuclei near the metal center can be overcome through the combination of uniform and selective 2H, 13C, and 15N isotopic labeling with NMR experiments that utilize direct one-dimensional (2H, 13C, and 15N) and two-dimensional (13C-X) detection. We have developed methods for calculating NMR chemical shifts and relaxation rates by density functional theory (DFT) approaches. We use the correspondence between experimental NMR parameters and those calculated from structural models of iron-sulfur clusters derived from X-ray crystallography to validate the computational approach and to investigate how structural differences are manifested in these values. We have applied this strategy to three iron-sulfur proteins: Clostridium pasteurianum rubredoxin, Anabaena [2Fe-2S] ferredoxin, and human [2Fe-2S] ferredoxin. Provided that an accurate structural model of the iron-sulfur cluster and surrounding residues is available from diffraction data, our results show that DFT calculations can return NMR observables with excellent accuracy. This suggests that it might be possible to use calculations to refine structures or to generate structural models of active sites when crystal structures are unavailable. The approach has yielded insights into the electronic structures of these iron-sulfur proteins. In rubredoxin, the results show that substantial unpaired electron spin is delocalized across NH...S hydrogen bonds and that the reduction potential can be changed by 77 mV simply by altering the strength of one of these hydrogen bonds. In reduced [2Fe-2S] ferredoxins, hyperfine shift data have provided quantitative information on the degree of valence trapping. The approach described here for iron-sulfur proteins offers new avenues for detailed studies of these and other metalloprotein systems.

Crystallographic studies of V44 mutants of Clostridium pasteurianum rubredoxin: effects of side-chain size on reduction potential.

Understanding the structural origins of differences in reduction potentials is crucial to understanding how various electron transfer proteins modulate their reduction potentials and how they evolve for diverse functional roles. Here, the high-resolution structures of several Clostridium pasteurianum rubredoxin (Cp Rd) variants with changes in the vicinity of the redox site are reported in order to increase this understanding. Our crystal structures of [V44L] (at 1.8 A resolution), [V44A] (1.6 A), [V44G] (2.0 A) and [V44A, G45P] (1.5 A) Rd (all in their oxidized states) show that there is a gradual decrease in the distance between Fe and the amide nitrogen of residue 44 upon reduction in the size of the side chain of residue 44; the decrease occurs from leucine to valine, alanine or glycine and is accompanied by a gradual increase in their reduction potentials. Mutation of Cp Rd at position 44 also changes the hydrogen-bond distance between the amide nitrogen of residue 44 and the sulfur of cysteine 42 in a size-dependent manner. Our results suggest that residue 44 is an important determinant of Rd reduction potential in a manner dictated by side-chain size. Along with the electric dipole moment of the 43-44 peptide bond and the 44-42 NH--S type hydrogen bond, a modulation mechanism for solvent accessibility through residue 41 might regulate the redox reaction of the Rds.

Strategy for the study of paramagnetic proteins with slow electronic relaxation rates by nmr spectroscopy: application to oxidized human [2Fe-2S] ferredoxin.

NMR studies of paramagnetic proteins are hampered by the rapid relaxation of nuclei near the paramagnetic center, which prevents the application of conventional methods to investigations of the most interesting regions of such molecules. This problem is particularly acute in systems with slow electronic relaxation rates. We present a strategy that can be used with a protein with slow electronic relaxation to identify and assign resonances from nuclei near the paramagnetic center. Oxidized human [2Fe-2S] ferredoxin (adrenodoxin) was used to test the approach. The strategy involves six steps: (1) NMR signals from (1)H, (13)C, and (15)N nuclei unaffected or minimally affected by paramagnetic effects are assigned by standard multinuclear two- and three-dimensional (2D and 3D) spectroscopic methods with protein samples labeled uniformly with (13)C and (15)N. (2) The very broad, hyperfine-shifted signals from carbons in the residues that ligate the metal center are classified by amino acid and atom type by selective (13)C labeling and one-dimensional (1D) (13)C NMR spectroscopy. (3) Spin systems involving carbons near the paramagnetic center that are broadened but not hyperfine-shifted are elucidated by (13)C[(13)C] constant time correlation spectroscopy (CT-COSY). (4) Signals from amide nitrogens affected by the paramagnetic center are assigned to amino acid type by selective (15)N labeling and 1D (15)N NMR spectroscopy. (5) Sequence-specific assignments of these carbon and nitrogen signals are determined by 1D (13)C[(15)N] difference decoupling experiments. (6) Signals from (1)H nuclei in these spin systems are assigned by paramagnetic-optimized 2D and 3D (1)H[(13)C] experiments. For oxidized human ferredoxin, this strategy led to assignments (to amino acid and atom type) for 88% of the carbons in the [2Fe-2S] cluster-binding loops (residues 43-58 and 89-94). These included complete carbon spin-system assignments for eight of the 22 residues and partial assignments for each of the others. Sequence-specific assignments were determined for the backbone (15)N signals from nine of the 22 residues and ambiguous assignments for five of the others.

Correlation between hydrogen bond lengths and reduction potentials in Clostridium pasteurianum rubredoxin.

15N NMR hyperfine-shift data were collected for wild-type and site-specific mutant (V44I, V44A, and V44G) Clostridium pasteurianum rubredoxins in the oxidized state. Whereas most of the (15)N NMR signals did not exhibit large systematic changes upon mutation of residue 44, the signal from the backbone nitrogen of residue 44 itself (arrows) shifted by approximately 400 ppm. These shifts were used to determine the lengths of the hydrogen bond between the backbone amide of residue 44 and the side-chain sulfur of cysteine-44, which is covalently ligated to the iron of the metal center. The results, which demonstrated that this hydrogen bond is shorter in mutants with higher reduction potential, point to the importance of hydrogen bonds in modulating the reduction potential of iron-sulfur proteins.

(13)C[(13)C] 2D NMR: a novel strategy for the study of paramagnetic proteins with slow electronic relaxation rates.

Oxidized human [2Fe-2S] ferredoxin has a notably slow electronic relaxation rate, which precludes the observation of signals from nuclei near the iron-sulfur cluster by conventional 2D or 3D methods that utilize proton detection. We have demonstrated the utility of (13)C[(13)C]CT-COSY in identifying connectivity information from fast relaxing carbon nuclei near the paramagnetic center, including those from residues that ligate the metal center.