Lawrence P. Wackett,1,2 Lawrence Que Jr.,2,3 John D. Lipscomb,1,2 and Douglas H. Ohlendorf1,2*
Department of Biochemistry, Molecular Biology and Biophysics,1 Department of Chemistry,3 Center for Metals in Biocatalysis, University of Minnesota, Minneapolis, Minnesota 554552
Received 25 July 2003/ Accepted 4 December 2003
| ABSTRACT |
|---|
|
|
|---|
ßßß modules forming two structurally homologous N-terminal and C-terminal barrel-shaped domains. The active-site metal is located in the C-terminal barrel and is ligated by two equatorial ligands, H214NE1 and E267OE1; one axial ligand, H155NE1; and two to three water molecules. The first and second coordination spheres of these enzymes are virtually identical (root mean square difference over all atoms, 0.19 Å), suggesting that the metal selectivity must be due to changes at a significant distance from the metal and/or changes that occur during folding. The substrate (2,3-dihydroxyphenylacetate [HPCA]) chelates the metal asymmetrically at sites trans to the two imidazole ligands and interacts with a unique, mobile C-terminal loop. The loop closes over the bound substrate, presumably to seal the active site as the oxygen activation process commences. An "open" coordination site trans to E267 is the likely binding site for O2. The geometry of the enzyme-substrate complexes suggests that if a transiently formed metal-superoxide complex attacks the substrate without dissociation from the metal, it must do so at the C-3 position. Second-sphere active-site residues that are positioned to interact with the HPCA and/or bound O2 during catalysis are identified and discussed in the context of current mechanistic hypotheses. | INTRODUCTION |
|---|
|
|
|---|
Dioxygenases that act on ortho-dihydroxylated aromatic compounds are divided into two classes, termed intradiol and extradiol, which differ in their mode of ring cleavage and the oxidation state of the active-site metal (44, 63). Intradiol dioxygenases utilize Fe3+ and cleave the bond between the carbons bearing the two hydroxyls; extradiol dioxygenases utilize Fe2+ or, rarely, Mn2+ and cleave one of the carbon-carbon bonds adjacent to the ortho-hydroxyl substituents. Although intradiol and extradiol dioxygenases oxidize an overlapping set of substrates, they are not structurally related and are proposed to have fundamentally different mechanisms (44, 45, 61, 62). Briefly, the mechanism of the intradiol dioxygenases is thought to involve activation of the aromatic substrate for electrophilic attack by molecular oxygen. In contrast, in the extradiol dioxygenase mechanism, the substrate activates the metal center for dioxygen binding and subsequent nucleophilic attack on the aromatic substrate.
On the basis of sequence alignments, the extradiol group of dioxygenases has been divided into three families (20, 22, 73, 80). As recently reviewed by Vaillancourt et al. (80), type I includes the 2,3-dihydroxybiphenyl 1,2-dioxygenases, such as the single-domain enzyme from Rhodococcus globerulus (5) and the two-domain enzymes from Pseudomonas putida KF715 (29), Pseudomonas sp. strain KKS102 (BPHC_PS102) (39, 68), Pseudomonas cepacia LB400 (BPHC_LB400) (27, 32), and Bacillus sp. strain JF8 (BPHC_JF8) (28). Type I also includes most of the catechol 2,3-dioxygenases (2,3-CTDs), such as the xylE gene product from the TOL plasmid of P. putida mt-2 (51). A consensus sequence for the metal binding region of type I enzymes has been defined (7, 22). The type II extradiol dioxygenases include such enzymes as protocatechuate 4,5-dioxygenase (LigAB) from Pseudomonas paucimobilis (53), which has two different types of subunits, and 2,3-CTD from Alcaligenes eutrophus JMP222 (38). The type III extradiol dioxygenases include gentistate dioxygenase from Haloferax sp. strain D1227 (24), human homogentisate dioxygenase (77), and 1-hydroxy-2-naphthoate dioxygenase from Nocardioides sp. strain KP7 (35). The type III enzymes have been identified as part of the cupin superfamily (20, 21).
Extradiol dioxygenases of each type have no intense visible chromophore or electron paramagnetic resonance (EPR) spectrum, which has slowed mechanistic and structural investigations. However, early studies using type I and II enzymes showed that NO, used as an O2 mimic, binds to the active-site Fe2+ to convert it to an EPR-active form (2, 44). Hyperfine interactions with 17O-labeled solvents and substrates, detected as a broadening of this EPR spectrum, showed that at least three sites in the iron coordination can be occupied by exogenous molecules (1, 3). Solvent bound initially to the iron is displaced by the substrate as both hydroxyl groups bind. Later, magnetic circular dichroism-circular dichroism (MCD/CD) (47) and extended X-ray absorption fine structure (EXAFS) (71) studies showed that the Fe2+ ions in both the resting enzyme and the substrate complex are pentacoordinate in the absence of NO, suggesting that ligands (probably the solvents) are displaced as the substrate binds and that there is an open site that could be occupied by O2 (or NO) as part of the activation mechanism. This would bring the substrate and oxygen into close spatial and electronic proximity in a prescribed orientation that could be controlled by the enzyme. Similar mechanistic strategies (12) have recently been invoked for a wide range of nonheme iron-containing enzymes of the "2-His, 1-carboxylate facial triad" group (30) which all share the property that one face of the active-site Fe2+ coordination is bound by solvents which can be displaced by substrates as part of the catalytic cycle. The crystal structures (27, 68) of 2,3-dihydroxybiphenyl 1,2-dioxygenase bore out all of the conclusions deduced from the earlier spectroscopic studies of related extradiol dioxygenases. One important finding from the EXAFS and EPR studies of 2,3-CTD was that the substrate chelates the iron asymmetrically, probably due to the deprotonation of only one hydroxyl function (71). This could serve to direct the position of oxygen attack. Also, it may differentiate extradiol from intradiol ring cleavage mechanisms, because the substrates bind as dianionic chelate complexes in the latter enzyme class (56).
Insight into the structural details of the extradiol dioxygenases has recently been gained through additional X-ray crystallographic studies. The majority of extradiol enzymes that have been structurally characterized belong to type I. Type I enzymes can be further subdivided based on their preference for substrates that have one or two aromatic rings. The crystal structure of at least one member of each subclass has been solved: P. putida mt-2 2,3-CTD (40), which catalyzes the cleavage of monocyclic aromatics, and BPHC_LB400 (27) and BPHC_PS102 (68), which catalyze the cleavage of bicyclic aromatics. These structures demonstrate that these enzymes have a common fold, though different oligomeric structures. The different folds found in the type II extradiol dioxygenase LigAB (76) and the type III extradiol dioxygenase homogentisate dioxygenase (77) suggest that these types arose independently of the type I enzymes. However, all three types appear to share similar active-site features and probably similar catalytic mechanisms.
Although the spectroscopic and crystallographic structural studies conducted thus far are generally in agreement, several questions remain that can potentially be addressed through structural studies of additional extradiol enzymes. Among these are the origin of the mechanistically relevant asymmetry in substrate binding in the active site, the identities and roles of the second-sphere active-site residues that may participate in catalysis, a rationale for the strict specificity of the ring cleavage site for extradiol and intradiol dioxygenases, and the origin of the remarkable ability of the homologous extradiol dioxygenases to select between iron and manganese as an active-site metal. In this study, we examine these aspects of extradiol dioxygenase catalysis by using two enzymes that both catalyze the oxygen-dependent extradiol cleavage of 3,4-dihydroxyphenylacetate (also known as homoprotocatechuate [HPCA]) to produce 5-carboxymethyl-2-hydroxymuconic semialdehyde. Interestingly, homoprotocatechuate 2,3-dioxygenase (2,3-HPCD) from Arthrobacter globiformis (MndD) utilizes Mn2+ (7, 86), whereas the enzyme isolated from Brevibacterium fuscum (Bf 2,3-HPCD) utilizes Fe2+ (48, 84). Primary sequence analysis indicates that MndD has 14 of the 18 amino acid residues previously identified as conserved in the type I family (7, 32). Despite this, MndD was found to have on average only 23% sequence identity with the type I family (7). In contrast, the sequence identity between MndD and Bf 2,3-HPCD is 82%. In comparison, the Mn2+-dependent enzyme BPHC_JP8 has less than 25% sequence identity with other 2,3-dihydroxybiphenyl 2,3-dioxygenases (28) and only 30% sequence identity with MndD. A multiple sequence alignment of MndD and Bf 2,3-HPCD with several type I extradiol dioxygenases indicates low but significant sequence homology, without the need for large gaps or insertions. However, MndD and Bf 2,3-HPCD contain an additional 20 to 30 residues at the C terminus beyond those present in the majority of type I extradiol dioxygenases. Here we report the refined structures of MndD, Bf 2,3-HPCD, and their substrate complexes. Active-site residues and their possible roles are identified through comparison of these structures with those of other extradiol dioxygenases. Although some insight is gained into the origin of metal ion specificity, the true basis of this important process is shown not to lie in the final folded active-site region, which is nearly identical in the two 2,3-HPCD enzymes.
| MATERIALS AND METHODS |
|---|
|
|
|---|
Prior to some crystallization trials, the proteins were modified proteolytically. Bovine pancreas trypsin (Sigma-Aldrich, St. Louis, Mo.) was added (final concentration, 2 mg/ml) to a concentrated solution (75 to 100 mg/ml) of Bf 2,3-HPCD or MndD. This mixture was then incubated at 18°C under an anaerobic hood for several days. Proteolysis was monitored by native polyacrylamide gel electrophoresis (PAGE) until the desired degree of cleavage was achieved (see Results for details).
Crystallization. Bf 2,3-HPCD crystallized in two forms. The P3221 form grew as hexagonal rods (0.15 by 0.15 by 0.5 mm) in hanging drops at 18°C when equal volumes of protein (20 mg/ml) were combined with a solution containing 1.5 M (NH4)2SO4 and 100 mM morpholinepropanesulfonic acid (MOPS) (pH 7.5). The I41 form grew as bricks (0.4 by 0.4 by 0.8 mm) at 18°C. Here 10 µl of proteolytically cleaved Bf 2,3-HPCD (100 mg/ml) was combined with 90 µl of a stock solution consisting of 10 to 15% polyethylene glycol 5000 monomethyl ether (MME), 0.2 M magnesium acetate, and 100 mM MOPS (pH 7.5) in a 96-well cell culture plate.
Thin plates of MndD (dimensions, 0.1 by 0.5 by 0.5 mm) grew from batch crystallization experiments with incubation at 18°C. Typically 10 µl of proteolytically cleaved MndD (100 mg/ml) was combined with 90 µl of a stock solution consisting of 6 to 10% polyethylene glycol 8000, 0.2 M magnesium acetate, and 100 mM sodium cacodylate at pH 6.5 in a 96-well cell culture plate. These crystals were found to belong to the monoclinic space group C2 with one dimer per asymmetric unit.
Data collection. A number of the original data sets were collected at room temperature on a Siemens HI-STAR area detector and processed by using XENGEN (33). X rays were generated by a Rigaku RU200 X-ray generator equipped with a copper anode operating at 45 kV and 200 mA. These included the MndD room temperature data set (Native 1), the Bf 2,3-HPCD trigonal data set (Native 3), and all of the Bf 2,3-HPCD tetragonal data sets used in heavy atom phasing.
High-resolution, low-temperature data sets of Bf 2,3-HPCD (Native 2) and MndD (Native 2) were collected at the Argonne National Laboratory at beamline 19-ID. Data were produced by using monochromatic radiation and were collected with a 3-by-3 charge-coupled device detector. Diffraction intensities were indexed with DENZO and scaled with SCALEPACK (58). Data were collected at -80°C by using a mother liquor containing 20% glycerol.
Crystalline enzyme-substrate (ES) complexes were prepared and mounted under anaerobic conditions under an argon atmosphere in a glove box (Coy Systems). Crystals were incubated for 1 h in anaerobic solutions containing the substrate (1 mM HPCA) and their respective mother liquors and were then mounted in quartz capillaries. A layer of mother liquor and then a layer of mineral oil were added to the top of the capillaries prior to their removal from the anaerobic chamber. The capillaries were then sealed with low-melt glue to maintain anaerobic conditions during the experiment. Diffraction patterns for the anaerobic ES complexes were measured at room temperature on an R-AXIS IV++ imaging plate detector and were processed by using CRYSTALCLEAR (Molecular Structure Corporation) (59). X rays were generated by utilizing a Rigaku RU200 X-ray generator operating at 45 kV and 100 mA, equipped with a copper anode and confocal MaxFlux Blue focusing optics (Molecular Structure Corporation).
Crystals of the ES complex remained colorless during incubation with the substrate and subsequent data collection, indicating that anaerobic conditions were maintained. After data collection, the capillaries were opened to expose the crystals to O2. The crystals immediately turned an intense yellow, the color of the product of the extradiol dioxygenase reaction. This suggested that both MndD and Bf 2,3-HPCD are active within the crystalline lattice and that the metal remained in its active state.
Phasing. Molecular replacement attempts using the structural homologues 2,3-CTD (Protein Data Bank [PDB] identifier, 1MPY) and BPHC_LB400 (PDB identifier, 1HAN) were unsuccessful. Therefore, phases for the Bf 2,3-HPCD I41 crystal form were obtained by multiple isomorphous replacement (MIR) utilizing K2HgI4, UO2(CH3COO)2, and an iodinated substrate analog (3-iodo, 4-hydroxyphenylacetate) to produce derivatives. Heavy atoms were located by difference Patterson and difference Fourier analysis using the CCP4 program suite (15); heavy atom positions were refined by using the MLPHARE program (57). Data and phasing statistics are shown in Table 1. The K2HgI4 derivative had excellent phasing statistics, especially at low resolution, producing Fourier maps with clear boundaries for solvent and protein. Although UO2(CH3COO)2 and the iodinated substrate analog were very weak derivatives, they facilitated the identification of the noncrystallographic symmetry operators. The MIR phases were refined and extended to 2.5 Å by solvent flattening, histogram matching, and fourfold noncrystallographic symmetry averaging using the DM program (16). A Fourier map based on the DM-refined MIR phases was sufficient to position 96% of the Bf 2,3-HPCD sequence.
|
and 1
were found in Fo-Fc and 2Fo-Fc maps, respectively, and when the solvent molecules formed good-geometry hydrogen bonds with the protein or other solvents. The van der Waals interaction energy was turned off for all ligands to the active-site metal so that the refinement of ligand-to-metal distances was driven purely by the electron density. The Bf 2,3-HPCD trigonal crystal form and the MndD structure were solved by molecular replacement using the CNS and AMORE (52) programs, respectively, with a molecular replacement search model from a partially refined monomer of Bf 2,3-HPCD (I41, Native 2) that had been stripped of all heteroatoms. Unambiguous solutions were obtained in both cases and were subsequently submitted to simulated annealing torsional dynamics. This reduced the initial Rfactor from 34.2 to 26.8% for the Bf 2,3-HPCD trigonal crystal form and from 37 to 28.9% for MndD. A high-resolution, low-temperature MndD data set (MndD C2, Native 2) became available at a later date and was used to complete the refinement.
Sequence correction. Two regions of the previously published MndD sequence (7) were discovered to be in error. The new sequences for residues 230 to 255 and 335 to 360 are DKMGALRISDRIERGPGRHGVSNAFY and GADGFSYTRKDETGEAAEGFKLGAQV, respectively. The resequenced portion from residue 230 to residue 255 is consistent with the electron density found for the MndD structure for this region. No electron density was available for the second resequenced region, because these residues were removed upon trypsin digestion. However, the new sequence is more consistent (i.e., homologous) with the sequence of the C terminus of Bf 2,3-HPCD.
Accession numbers. The atomic coordinates for the structures of Bf 2,3-HPCD and MndD have been deposited in the Protein Data Bank. Accession numbers are given in Table 2.
|
| RESULTS |
|---|
|
|
|---|
It was noted that a proteolytic fragment of MndD occasionally appeared during purification and that when this occurred, better crystals were obtained. Controlled digestion by trypsin of Bf 2,3-HPCD and MndD produced more-homogeneous enzymes as judged by native PAGE. Trypsin-treated Bf 2,3-HPCD grew in a new crystal form (I41). Both cleaved proteins were active and had no substantial change in Vmax or Km.
Overall model quality. All models presented in this study were refined against their respective data sets to Rfactors less than 18.0% and Rfrees less than 22.0%. Structural analysis using the PROCHECK program (42) demonstrated that all models conformed to typical protein stereochemistry and exhibited >90% of their residues in the most favored or allowed regions of a Ramachandran plot. Refinement statistics for the various models are shown in Table 2.
In general, the first two to three residues of Bf 2,3-HPCD and MndD are disordered and therefore were not modeled. Also, the data sets typically showed no density beyond residue 322 or 323. The exception was the P3221 crystal form of Bf 2,3-HPCD, for which the data did yield sufficient electron density to model 37 of these C-terminal residues.
Oligomeric structure. Bf 2,3-HPCD and MndD form very similar tetramers, each with 222 point symmetry and overall dimensions of approximately 75 by 75 by 110 Å. Within the tetramer, each monomer has distinct interactions with the other three monomers. The most extensive interaction is that between monomers A and B (or C and D), where a 3,369-Å2 surface area is buried upon formation of the dimer. The interaction of dimer AC with dimer BD is composed almost entirely of solvent-mediated polar interactions. In this interface there are three solvent pockets: one at the center of the A-B interface, one at the center of the C-D interface, and one around the 222 symmetry point. Conversely, the interactions between monomers A and C (or B and D) and between monomers A and D (or B and C) bury only 1,733 and 1,006 Å2, respectively. The interaction of the AB dimer with the CD dimer includes both hydrophobic and solvent-mediated polar interactions. There are no biochemical data to suggest that either enzyme exists as anything but tetramers in solution (48, 86). The active-site cavities are located approximately along the long edges of the rectangular tetramers. The active-site metal of monomer A is 31, 33, and 44 Å away from the corresponding metals of monomers B, C, and D, respectively. Biochemical studies of Bf 2,3-HPCD and MndD do not indicate any allostery between monomers (48, 86).
The C-terminal 43 residues of Bf 2,3-HPCD run transverse from the corners of the tetramer, down the long axis, and occupy space between the A-C and B-D monomer contacts. As stated earlier, these residues can be removed upon incubation of Bf 2,3-HPCD or MndD with trypsin. Inclusion of the C terminus, as seen in the P3221 crystal form of Bf 2,3-HPCD, increases the width of the tetramer to 80 Å. This small increase in volume is not consistent with the large decrease in volume apparent upon treatment with trypsin, as viewed by native PAGE, suggesting that the C termini are mobile and can take extended conformations.
Monomer structure.
The monomer structures of Bf 2,3-HPCD and MndD are composed of two structurally homologous barrel-shaped domains: an N-terminal domain (residues 1 to 145), a C-terminal domain (residues 145 to 290), and a lid domain (residues 291 to 353) (Fig. 1). The substrate-binding determinants and the active-site metal are located within the C-terminal domain barrel. Each barrel can be further subdivided into two modules, each of which has the same ß
ßßß motif: modules A and B (N-terminal domain) and modules C and D (C-terminal domain). The four ß strands of each module combine to form a mixed-topology ß sheet with an amphipathic helix lying against the ß sheet. The basic fold of the individual modules has been described for a number of other structures (4, 6, 25) including several that are not type I extradiol dioxygenases, i.e., human glyoxylase I (14), bleomycin resistance protein (19), and Pseudomonas fluorescens 4-hydroxyphenylpyruvate dioxygenase (69). This superfamily has been labeled the vicinal oxygen chelate superfamily (4).
|
The monomer structures of Bf 2,3-HPCD and MndD are essentially identical (RMS difference over 319 structurally equivalent C
s, 0.250 Å), as expected from their high sequence identity. MndD and Bf 2,3-HPCD also have high structural homology with the individual subunits of two other type I extradiol dioxygenases, 2,3-CTD (RMS difference over 285 C
s, 1.65 Å) and BPHC_LB400 (RMS difference over 240 C
s, 1.79 Å), despite sequence identities of <25%. BPHC_LB400 is a homooctamer, while the other proteins assemble as homotetramers with similar structural features. Additionally, 2,3-CTD and BPHC_LB400 do not have H2L, S3L, and S4L of the lid domain, because this region is absent in the these enzymes.
Active-site cavity. We believe that the trypsin-cleaved structures are representative of the conformation the structure assumes before binding the substrate. Evidence for this will be presented in the following sections. The active-site metal is located approximately 15 to 20 Å within the C-terminal barrel. The outer rim of the active site has a very broad, flat opening constructed by residues 290 to 305 of the lid domain, the loop between strands S2C and S3C, and the loop between strands S4C and S1D. An approaching substrate would first encounter residues Y189, V287, R292, Y305, V300, G206, and G207. Farther in, the substrate would encounter a circular opening, approximately 8 to 10 Å in diameter, defined by the side chains of R293, I181, A190, W304, R243, and H213. The side chain of R293 is disordered in the substrate-free, trypsin-cleaved forms of MndD and Bf 2,3-HPCD. There is ample room for the substrate to enter if R293 and H2L of the lid domain move out of the way. Finally, upon interacting with the active-site metal, the substrate would settle in a pocket formed by the side chains of W192, H248, V250, S251, and V257.
Native metal coordination. The electron densities for the ligands to the active-site metals of Bf 2,3-HPCD and MndD are shown in Fig. 2A and C, respectively. All the ligands to the metal are provided by the C-terminal domain and arise from ß strands. The active-site metal is coordinated by two histidines, one glutamic acid, and two to three solvent molecules, depending on the crystal form, in octahedral geometry. This confirms the previous mutagenesis work, which showed that these three residues are critical to both activity and the ability to bind manganese in MndD (8). In our terminology we identify H214NE2 and E267OE1 as equatorial ligands and H155NE2 as an axial ligand. The three potentially open ligand positions (A, B, and C) are trans to H214NE2, E267OE1, and H155NE2, respectively. In two of the native structures (Bf 2,3-HPCD, Native 2, and MndD, Native 1), the three open positions are, in fact, occupied by solvent molecules. The solvent in site A appears to be the weakest ligand to the iron, because it typically has the highest B-factor and a significantly longer bond length, and it is absent in one of the native structures (Bf 2,3-HPCD, Native 3). At this resolution, changes in bond lengths of 0.25 to 0.3 Å are deemed significant. In the other sites, the average ligand distance is 2.1 to 2.2 Å depending on the native structure examined (Table 3). Having a weakly bound solvent or nothing in site A is consistent with the MCD/CD (47) and X-ray absorption spectroscopy (71) studies of 2,3-CTD, which show a pentacoordinate metal center with five N or O ligands at an average ligand-to-metal distance of 2.09 Å. It is assumed that, on average, one of the solvent molecules exists as a hydroxide, resulting in charge neutrality at the active site. The solvent molecule at site C has the shortest bond lengths, suggesting that it has the most hydroxide character. A superposition of the MndD and Bf 2,3-HPCD active sites does not indicate any discernible differences in the metal coordination geometry or the positions of critical residues (RMS deviation over 64 atoms, 0.12 Å).
|
|
The Bf 2,3-HPCD-substrate complex demonstrates that Fe2+ ligation also occurs through both hydroxyls of the substrate. The most significant difference between the two complexes is the absence of the solvent molecule in site B, leading to a square pyramidal Fe2+ geometry leaving one open coordination site. In both complexes the binding is asymmetric, with HPCAO3 exhibiting the shorter bond to the iron (Table 3).
Several residues in the active-site cavity appear to be critical to the positioning of the substrate (Fig. 3). The substrate is wedged into a planar crevice. One wall is formed by Y257 and H248, both of which are conserved among all members of the type I extradiol dioxygenase family. Y257OH forms a hydrogen bond with HPCAO3, while the center of the imidazole ring of H248 stacks with HPCAC5 and HPCAC6. On the opposite side of the substrate, the side chain of W192 is approximately perpendicular to the plane of the substrate where W192NH interacts with the
-electron cloud of HPCA. The corresponding residue in 2,3-CTD and BPHC_LB400 is a phenylalanine, which is similarly oriented. At the bottom of the crevice, the edge of the aromatic ring (HPCAC5, HPCAC6) is within van der Waals contact with the side chains of V250 and S251 and with the backbone carbonyl of V250. The close proximity of these residues would severely restrict the abilities of MndD and Bf 2,3-HPCD to catalyze the cleavage of any substrates with bulky substituents at these positions. On the opposite side of the ring, HPCAC2 has more distant contacts; it is approximately 4.5 Å from the side chain of W304. W304 is therefore an important determinant of what substituents could be placed at C2 and still retain functional binding. In BPHC_LB400, the benzyl substituent at this position is accommodated by having a much abbreviated lid domain.
|
Major movements upon substrate binding. The binding of the substrate is facilitated by and/or facilitates the movement of several residues in the active-site cavity. The movements appear to propagate outward from several residues directly adjacent to the substrate.
In the native structure, a hydrogen bonding network links Y257OH with H248NE2 and H248ND1 with R293O. When a substrate binds, these interactions are disrupted, as also noted for the analogous residues in BPHC_LB400 (79). Y257OH forms a new hydrogen bond with HPCAO3. The plane of the imidazole ring of H248 slides approximately 1.4 Å along the plane of the substrate ring. The hydrogen bond of Y257OH with H248NE2 is broken, and H248NE2 forms a new hydrogen bond with E267OE2, one of the ligands to the metal. E267OE2 shifts slightly to optimize this interaction with H248NE2. The movement of H248 also breaks the hydrogen bond of H248ND1 with R293O. This interaction is replaced by one with the acetate group of HPCA, which forms a hydrogen bond with both H248ND1 and R293O. These movements permit the acetate group of the substrate to interact with the guanidinium group of R243. A smaller (0.6-Å) but similar shift of the homologous H241 of BPHC_LB400 occurs when a substrate binds that enzyme (79).
In the native structure, D294 forms a hydrogen bond with R243NE. Upon substrate binding, the conformation of D294 switches from gauche- (
1 = -90°) to trans (
1 = -168°), away from R243 and toward the 293-to-302 loop. This appears to be caused by the movement of R293 (due to the ordering of the side chain) and the loss of the hydrogen bond of R293O with H248ND1. In addition, the interaction of the acetate group of HPCA with R293 means that the electrostatic component of the interaction between R293 and D294 is no longer required. The ordering of the side chain of R293 causes the side chain of Y305 to move outward toward the solvent by approximately 1 to 2 Å. These displacements permit V300CB to move approximately 2.9 Å inward toward R243. There are also large changes in the positions of P299, V301, and P302, as this whole loop packs closer to the protein.
A second distinct set of displacements propagates from near the active-site metal to the interface between the A and B (or C and D) subunits. In the native structure, a solvent-filled pocket directly adjacent to the active-site metal is separated from a second solvent-filled pocket in the AB (or CD) dimer interface by the side chains of Y269 and D272. In the ES complex, Y269 and D272 have moved away from the dimer interface toward the substrate while maintaining the hydrogen bond of Y269OH with D272OE1. The side chain of Y269 fills a large portion of the solvent pocket adjacent to the iron, where Y269CE2 is only
3.6 Å from HPCAO4. These displacements cause the side chain of Q271, a residue in the AB (or CD) dimer interface, to change conformations (from a
1 of -150° to a
1 of -76°). The trigger for the second group of conformational changes is unclear, but they may be due to a change in the conformational stability in this area when solvents are expelled upon substrate binding.
A superposition of the Bf 2,3-HPCD and MndD substrate-bound structures indicates that all of the conformational changes upon substrate binding are conserved between the two structures. The most telling conformational changes are the ordering of the side chain of R293, the movement of the 293-to-302 loop (most specifically the side chain of V300), and the conformational change of Y269. The structures, therefore, suggest that the active sites of MndD and Bf 2,3-HPCD can be in an open or a closed conformation depending on the binding of the substrate.
P3221 crystal form. Analysis of native and cleaved Bf 2,3-HPCD and MndD by native PAGE and by their propensities to crystallize indicated that there was most likely a flexible amino or carboxy terminus on these proteins. Sequence alignment of several enzymes from the extradiol family with a diverse substrate range indicated that Bf 2,3-HPCD and MndD have unique carboxy termini. The function of these residues is unknown. We were able to obtain one data set from a crystal form of Bf 2,3-HPCD that did not require trypsin cleavage. Interestingly, 36 of the missing 42 C-terminal residues became visible in the electron density. The average B-factor for these residues is 25.5 Å2. As discussed earlier, the portion of the lid domain that remained after trypsin cleavage reduced the accessibility of the active site to small molecules. However, since the side chain of R293 was disordered, the mouth of the active site still permitted access to a molecule the size of HPCA. The structure of the P3221 crystal form demonstrates that a helix from the trypsin-sensitive portion of the lid domain is capable of covering the remaining opening of the active-site cavity. Figure 4 shows the electron density around R293 in the P3221 crystal form of Bf 2,3-HPCD. The size and shape of the electron density in the coordination sphere of the Fe2+ available to exogenous substrates suggested a partially occupied carboxylate, e.g., acetate. However, since acetate was not present in the crystallization solutions, we chose to model this density as two solvent molecules.
|
| DISCUSSION |
|---|
|
|
|---|
Effect of the carboxy terminus on crystallization. It was observed in this study that the space group and quality of the crystals obtained correlated with the presence or absence of the C-terminal fragment. The I41 space group of cleaved Bf 2,3-HPCD has one tetramer per asymmetric unit, with each monomer exhibiting unique contacts with crystallographically related molecules. In each case, if the flexible portion of the C terminus remained bound in the closed conformation, it would result in direct collisions with neighboring molecules and would prevent the formation of this crystal form. A similar comparison with the C2 crystal form of MndD suggested that the C-terminal fragment could fit comfortably within the solvent space of the crystal and not disrupt crystal contacts. However, the presence of the flexible C-terminal fragment would most likely interfere with the conformation of the immediately preceding residues P318 to V323, which form the most intimate crystal contacts. This postulate is consistent with the fact that disordered MndD crystals were obtained prior to treatment of the protein with trypsin. In the P3221 crystal form of Bf 2,3-HPCD, there are no specific contacts between the flexible C-terminal fragment and neighboring molecules in the lattice, indicating that order is not imposed by interactions with neighboring molecules. Thus, not only could a nearly complete lid domain be seen, but these observations also suggested that the lid domain has the probable function of closing over the active site after the substrate binds.
Metal ion specificities of MndD and Bf 2,3-HPCD. MndD and Bf 2,3-HPCD display a high degree of sequence identity, which fails to reveal any obvious basis for the observed metal ion specificity. The structures presented here permit an analysis of residues in and around the active site that might influence the binding and reactivity of a particular metal. Surprisingly, there are no significant differences in the positions of atoms within the first (H155, H214, and E267) or the second (H200, N157, H213, Y257, D154, and W192) sphere of the active-site metal. However, there are a number of amino acid differences at the interface between the two ß barrels of the C-terminal domain (10 to 15 Å from the metal). The amino acid differences are complementary exchanges of hydrophobic residues so that the character of the interface is maintained. The reordering of this interface could have a long-range effect on the metal selectivity, but no clear rationale for this emerges. There are also a number of differences in the sequence or position of residues in the solvent pocket between the A and B (or C and D) subunits. One of these residues is R152, which in MndD forms a hydrogen bond with the side chain of D154OD1, while in Bf 2,3-HPCD it extends out into the cavity, forming hydrogen bonds with D272'O (symmetry-related monomer in tetramer), E57OE1, and several solvent molecules (Fig. 5). Since the side chain of D154 also forms a hydrogen bond to the imidazole of H155, a metal ligand, the position of R152 could have an indirect effect on the electronic character of the metal site. It is not clear if these differences are significant enough to have any impact on metal selectivity.
|
Enzyme mechanism. The high-resolution structures of the substrate-free and -bound forms of MndD and Bf 2,3-HPCD allow us to evaluate and expand the proposals that have been made for the mechanism of the extradiol dioxygenases (3, 13, 48, 71). Key questions relate to the mechanism of oxygen activation, the attack on the substrate, and the specificity of the ring cleavage site. A mechanistic proposal based on our structural insight is shown in Fig. 6.
|
2.7 Å). A similar role has been proposed for the analogous tyrosine residues in BPHC_PS102 and BPHC_LB400 (66, 81). The notion of a monoanionic ES complex for extradiol dioxygenases has been suggested previously by X-ray absorption spectroscopy studies of 2,3-CTD (71), UV resonance Raman and crystallographic studies of BPHC_LB400 (79), and studies of 4-nitrocatechol binding to these enzymes (64, 78). The asymmetry of the complex offers a possible means by which the site of attack of activated oxygen on the aromatic ring and, from this, the regiospecificity of ring cleavage are determined. Oxygen binding. Studies utilizing NO as an oxygen analog showed that NO coordinates to the Fe2+ center of extradiol dioxygenases with an affinity that is dramatically enhanced by substrate binding (3). By extension, O2 is proposed to bind to the metal center in the solvent binding site (site C) that remains after the substrate is bound. This assignment is strongly supported by the crystal structure of BphC with NO bound in the analogous site (66). It may be necessary to stabilize the M2+-O2 bond to prevent adventitious oxidation of the metal center. One residue that that might carry out this function is H200, a conserved residue in the active sites of extradiol dioxygenases (27, 65, 68). The present structures show that H200NE2 is well positioned to interact with a metal-bound dioxygen and stabilize it, much like the function of the distal histidine in myoglobin and hemoglobin (49, 50, 54, 60, 70, 72, 85). H200 could also carry out a second function, since it is approximately 3.5 Å from HPCAO4. If this hydroxyl must be deprotonated as part of the oxygen insertion reaction of the reaction mechanism, H200 could serve as a general base for this reaction, as suggested for the equivalent reaction in BphC (68). A protonated H200 would also be more effective in stabilizing an M2+-O2 bond, so oxygen binding may promote removal of the second proton from HPCA. However, it must also be stated that the H200NE2 geometry is not ideal for the formation of a reasonable hydrogen bond with HPCAO4; reorganization would be required. Groce and Lipscomb have recently directly demonstrated the importance of H200 by generating the H200F mutant and showing that the rate of turnover is greatly decreased despite little effect on the substrate-binding reaction itself (26). This suggests that the primary effects of this residue are exerted on later steps in the reaction cycle, including oxygen binding, activation, and insertion. Interestingly, the H200F mutant also switches the ring cleavage specificity for some alternative substrates from extradiol to intradiol, suggesting that it has a profound effect on directing either the reactivity or the relative geometry of the substrate and oxygen in the active site.
Oxygen activation and attack on substrate. It has been proposed that, subsequent to O2 binding, the metal center acts as a conduit for the flow of electron density from the electron-rich substrate to the bound dioxygen (3, 48), generating a semiquinone-M2+-superoxide intermediate. The oxygen activated in this way would then attack the adjacent substrate at C-2 or C-3 to form a peroxy intermediate. The proximity of Y257 to the axial binding site forces HPCAC3 closer to the putative oxygen binding site, as illustrated by the 15 to 20° deviation of HCPAO3 from octahedral geometry. The present structures show that the substrate would have to rotate about the HPCAC3/C5 axis in order to bring HPCAC2 close enough to form the M2+-O-O-HPCAC2 peroxy adduct favored by earlier mechanistic postulates (3, 48, 71). This rotation is unlikely because of the many interactions the acetate group makes with the protein and the potential for collisions between the substrate and W192. On the other hand, the structure shows that an M2+-O-O-HPCAC3 peroxy adduct can be formed easily without much distortion or dissociation of the incipient superoxide. A similar conclusion was reached for the BphC-substrate complex (66), suggesting that this may be a consistent aspect of the extradiol dioxygenase mechanism.
Ring-opening reaction and its regiospecificity. The ring-opening reaction is envisioned as a Criegee rearrangement in which the O-O bond is broken and one atom is inserted into the ring in a concerted process. This reaction would be promoted by removal of the proton from the second hydroxyl group of HPCA and proton donation to the peroxy intermediate undergoing O-O bond cleavage by one or more active-site acid or base catalysts as discussed above. Following the ring insertion chemistry, hydrolysis of the intermediate by the remaining oxygen atom from O2 bound to the metal center would complete the reaction. This sequence would ensure the overall dioxygenase stoichiometry observed for these enzymes. The proposed Criegee-type oxygen insertion reaction could occur in one of two ways, yielding either the C2-C3 lactone species that would lead to extradiol cleavage or the C3-C4 anhydride species that would yield the intradiol product after hydrolysis, as shown in Fig. 7.
|
It has also been suggested that the precise stereochemistry of the metal-peroxo-substrate unit determines the sites of cleavage (13). In this proposal, the axial orientation of the peroxo-substrate bond (see Fig. 6) that results from the presence and relative positions of the three available binding sites on the metal, together with the ferrous oxidation state of the iron and the presence of a proton donor to promote the Criegee rearrangement, favors alkenyl group migration to form the C2-C3 lactone intermediate. The present structural studies show that each of these features is present in the active sites of the 2,3-HPCD enzymes. Mutational studies that examine the effects of different metal ligands and active-site residues are needed to test these mechanistic proposals. The structural studies described here provide a starting point for such mutational studies, which are currently under way.
| ACKNOWLEDGMENTS |
|---|
Special thanks go to Stephan L. Ginell, Norma E. C. Duke, Rongguang Zhang, and Andrzej Joachimiak, who assisted in the collection of data sets obtained on beamline ID-19 at the Argonne National Laboratory.
| FOOTNOTES |
|---|
Present address: Department of Biochemistry, Albert Einstein College of Medicine, Bronx, NY 10461-1602. ![]()
| REFERENCES |
|---|
|
|
|---|