Crystal structures of substrate-free and biphenyl-bound forms of BphAEII9.
The crystal structures of the substrate-free and biphenyl-bound forms of BphAE
II9 were determined to resolutions of 2.5 Å and 1.9 Å, respectively. Representative crystallographic data and refinement statistics are presented in
Table 1.
The crystal structure of the substrate-free form was refined at a resolution of 2.5 Å to final
Rcryst and
Rfree values of 18.1% and 23.4%, respectively, with four αβ-heterodimers in the asymmetric unit of space group H3. Three dimers form one biological hexamer (α
3β
3) associated with the asymmetric unit, and the fourth forms a hexamer with symmetry-related heterodimers. The final model includes residues 18 to 143 and 153 to 459 for all α-subunits and residues 6 to 188 for all β-subunits. Interpretable electron density was not observed for residues 144 to 152 of the α-subunit. These residues were also absent in the wild-type BphAE
LB400 structure (
11). Similar to other Rieske-type dioxygenases, each α-subunit binds a [2Fe-2S] Rieske-type cluster and a mononuclear Fe
2+ ion at the active site, which is coordinated to two His residues (His233 and His239), an Asp residue (Asp388), and one or more water molecules.
The biphenyl-bound structure was refined to a final Rcryst of 21.6% and an Rfree of 25.4% at a resolution of 1.9 Å in space group P1. The structure includes 12 crystallographically independent αβ-heterodimers. α-subunits are labeled chains A, C, E, and G, etc., and β-subunits are labeled chains B, D, F, and H, etc. The final model includes residues 18 to 143 and 153 to 459 for all α-subunits except chain U (missing residues 143 and 459) and chain W (missing residue 143). For the β-subunits, all chains extend to residue 188, but the N-terminal content varies among the chains: chain F begins at residue 5, chain J begins at residue 14, and the remaining chains begin within the range of residues 6 to 13.
In eight of the α-subunits (chains C, E, I, K, M, O, Q, and W), initial unbiased (
Fo −
Fc) Fourier maps showed bulky electron density at the active site, consistent with bound biphenyl (
Fig. 1B), whereas such density was not observed for the other four α-subunits (chains A, G, S, and U). Density consistent with a partially occupied biphenyl was observed for chain S at a later stage.
The absence of biphenyl in some active sites appears to be a consequence of crystal packing contacts associated with structural variation in the N-terminal segments of the β-subunits. For most β-subunits, the N-terminal segment interacts with residues of the same subunit or an adjacent β-subunit of the same α3β3 hexamer. However, for three β-subunits (chains D, F, and J), the N-terminal segment extends away from the rest of the β-subunit and interacts with a loop formed by α-subunit residues 247 to 263 at the mouth of one of the actives sites of a neighboring α3β3 hexamer. This contact apparently interferes with the binding of biphenyl. In fact, for chains D and F, the side chain of β-Phe9 extends just into the mouth of the active site of a neighboring α-subunit (chains A and G, respectively) and would prevent the entry of the substrate into the active site, as shown in Fig. S1 in the supplemental material. In the case of chain J, β-Phe9 was not modeled because interpretable electron density begins at residue 14. Nevertheless, the N-terminal segment in chain J extends toward chain U, and density for biphenyl was not observed in the active site of U.
Comparison among α-subunits of the substrate-free and biphenyl-bound structures of BphAEII9.
The C-α atoms of the α-subunits from the substrate-free and biphenyl-bound BphAEII9 structures, 16 chains in total, can be superposed with a C-α RMSD of 0.5 Å or less, demonstrating close agreement. In some regions, however, variations in the local structure are observed.
Examination of superposed αβ-dimers reveals that the BphAE
II9:biphenyl complex shows greater variability in the position of key active-site elements, as shown in
Fig. 2B, than either the substrate-free BphAE
II9 structure (
Fig. 2C) or the biphenyl-bound form of its progenitor, BphAE
LB400 (
Fig. 2A). For example, in the latter two crystal structures, the superposed positions of the Fe
2+ ions lie in single, tight clusters. In contrast, in the BphAE
II9:biphenyl complex, five Fe
2+ ions (chains A, C, E, G, and U) are in one cluster, and six (chains I, K, M, O, Q, and W) are found in a second cluster ∼1.3 Å distant (average of all intercluster Fe
2+-Fe
2+ distances). Although the first cluster includes all chains for which biphenyl is not modeled (chains A, G, and U), two active sites with biphenyl present (chains C and E) are also members. The Fe
2+ ion from chain S lies approximately between the two clusters and nearly equidistant (0.8 Å) from both.
Comparatively greater positional variability is also observed in the BphAE
II9:biphenyl complex for two of the Fe
2+-coordinating residues, Asp388 and His239, as well as Ser283, Val287, Gly321, and Gln322. The case of Asp388 is highlighted in
Fig. 2A to
C. Although the carboxylate groups of Asp338 for all chains lie close to a common plane, the C-δ atoms are as much as 1.5 Å apart. The variations for Asp388 extend to the backbone atoms: the average (0.8 Å) and maximum (1.5 Å) distances among C-α atoms are indistinguishable from the average and maximum distances for C-γ and Oδ1, which binds to the Fe
2+ ion. Ligands would be expected to track with the Fe
2+ ion, and thus, clusters are again observed for the Asp388 carboxylates, although they are not as distinct or as tightly clustered.
To understand the variations in Fe2+ and Asp388 locations, it is useful to consider them with reference to a vector from the C-α atom of Gly321 to the C-α atom of Asp388, which passes close to the center of biphenyl and close to the Fe2+ ion. The length of the vector is significantly shorter for chains without density for biphenyl (ŕ = 15.0 Å; range = 15.0 Å to 15.1 Å) than for chains with biphenyl (ŕ = 16.5 Å; range = 15.9 Å to 17.3 Å). Qualitatively, expansion of the active site along this vector moves Fe2+ and Asp388 away from Gly321, opening space for biphenyl.
The orientation of the peptide plane between Gly321 and Gln322 and the interactions of the carbonyl of Gly321 also differ between the substrate-free and biphenyl complexes and vary among the chains of the latter (
Fig. 3A). In all four chains of the BphAE
II9 substrate-free structure, ψ is in the range of 12° to 37° such that the carbonyl of Gly321 points into the active site and interacts with C-ε of Met231 through a distance of 3.0 Å to 3.2 Å. For the biphenyl-soaked structure, the carbonyl of Gly321 is similarly placed in six of the chains, including all three without density for biphenyl (chains A, G, and U) and three with biphenyl modeled (chains M, Q, and W). In the other six chains, all with biphenyl modeled (chains C, E, I, K, O, and S), ψ is in the range of 95° to 115°, such that the direction of the carbonyl of Gly321 differs by ∼90°. This allows the carbonyl to form a hydrogen bond with the amide NH of Tyr277 through a distance of 3.0 Å to 3.1 Å. As in the case of the Fe
2+ location, the local structures of chains without biphenyl modeled are consistent, but the chains with biphenyl bound disperse between two groups, one that is consistent with the substrate-free structure and one that is distinct.
A difference in the extent of variability between the biphenyl complex and the substrate-free enzyme is also seen for the short α-helix (residues Pro281 to Met288) that forms one side of the active-site portal and one wall of the biphenyl binding site, where the side chains of Ser283 and Val287 are in contact with the distal ring of biphenyl. In the BphAEII9 biphenyl complex, the maximum difference in the positions of Ser283 C-α atoms for any pair of chains is 1.9 Å, but among the chains in the BphAEII9 substrate-free structure, the maximum value of the same measure is 0.6 Å.
Finally, a comparable extent of variability is seen in both structures for the backbone and side-chain atoms from residues Pro249 to Thr260. Here, maximum C-α displacements of up to 0.8 Å and 2 Å in the BphAEII9 biphenyl-bound structure and of up to 3.8 Å and 1.1 Å in the BphAEII9 biphenyl-free structure are seen at residues Ser254 and Ile258, respectively. A comparable variation is observed for side-chain atom positions and conformations. This segment contains a large fraction of solvent-exposed residues and does not contribute directly to biphenyl binding.
Transformation of 2,3,4′-trichlorobiphenyl, 2,2′,5,5′-tetrachlorobiphenyl, and DDT and docking studies.
To provide more insight into the structural features of BphAEII9, BphAELB400, and BphAEB356 responsible for their differential abilities to metabolize biphenyl analogs, we examined the biochemistry and the structural interaction of these enzymes toward 2,2′,5,5′-tetrachlorobiphenyl, 2,3,4′-trichlorobiphenyl, and DDT.
Biochemical studies were performed to assess the ability of BphAEII9 to transform DDT and two PCB congeners, 2,3,4′-trichlorobiphenyl, 2,2′,5,5′-tetrachlorobiphenyl, and the regiospecificity of the reactions.
BphAE
LB400 generated one major metabolite from 2,3,4′-trichlorobiphenyl as determined by examination of the GC-MS profile of the assay solution (
Fig. 5A). The mass spectral features of this metabolite's
nBuB derivative were consistent with those of a dihydro-dihydroxy-trichlorobiphenyl, namely, a fragmentation pattern exhibiting a molecular ion at
m/z 356 and diagnostically important ions at
m/z 340 (M
+ − O),
m/z 321 (M
+ − Cl),
m/z 286 (M
+ − 2Cl), and
m/z 256 (M
+ −
nBuBO
2). A second metabolite exhibiting mass spectral features of a dihydro-dihydroxy-trichlorobiphenyl was also detected, but these features are not observable in
Fig. 5A, as this metabolite was present in only trace amounts. BphAE
II9 transformed 2,3,4′-trichlorobiphenyl into the same major metabolite, as shown in
Fig. 5A; however, the amount produced was significantly larger. BphAE
B356 produced the same metabolite from 2,3,4′-trichlorobiphenyl, and the amount produced was similar to that produced by BphAE
LB400 (not shown). Therefore, our assay demonstrated that the major product produced by three enzymes is the same metabolite and that BphAE
II9 has a greater ability to degrade 2,3,4′-trichlorobiphenyl, as demonstrated by the much larger peak for the metabolite observed in the chromatogram. The substrate was docked into the active site of the three enzymes.
For BphAELB400, 2,3,4′-trichlorobiphenyl docked with the C-3 and C-3–chlorine atoms in approximately the same location as the C-2 and C-4 atoms of biphenyl in its complex with BphAELB400, and the proximal ring is in approximately the same plane as found in the biphenyl complex. This places the dichlorinated ring in a location and orientation consistent with dioxygenation across the C-4—C-5 bond. In this pose, as shown in Fig. S2B in the supplemental material, C-4 and C-5 are 4.3 Å and 4.4 Å from the Fe2+ ion, respectively, and 2.8 Å and 3.0 Å from a water molecule that occupies the presumed binding site for dioxygen, respectively.
For BphAE
II9, docking of 2,3,4′-trichlorobiphenyl produced a pose consistent with the same regio- and stereospecificity as those for BphAE
LB400, but the binding mode is remarkably different (
Fig. 6A and
B; see also Fig. S2B in the supplemental material). The dichlorinated ring does not penetrate as deeply into the active site, such that its location is comparable to that of the distal ring in the BphAE
II9:biphenyl complex. In addition, relative to the biphenyl complex, the orientation of the dichlorinated ring is shifted by ∼70°, and the torsion angle between the rings is shifted by +110°. Nevertheless, the distances from C-4 and C-5 to the water molecule are 3.3 Å and 3.0 Å, respectively, and the distances to the Fe
2+ ion are 5.1 Å and 4.5 Å, respectively. For BphAE
B356, 2,3,4′-trichlorobiphenyl docked in a manner similar to that of BphAE
II9, as shown in Fig. S2A in the supplemental material, but with reduced docking scores.
Our docking study showed that the preferred docking pose for all three enzymes was one that positioned the substrate to be converted to a 3,4-dihydro-3,4-dihydroxy-2,3,4′-trichlorobiphenyl. This proposed product would not be further degraded into a chlorobenzoate, which may explain why Seeger et al. did not detect any chlorobenzoate produced from 2,3,4′-trichlorobiphenyl in their experiment with BphAE
LB400, since BphAE
LB400 metabolizes this substrate through 3,4-dioxygenation (
42). Furthermore, the docking study showed that the preferred pose for BphAE
LB400 oriented the substrate much farther into the active site than the preferred pose for BphAE
II9 (
Fig. 6B; see also Fig. S2B in the supplemental material). This deeper-binding mode may explain why BphAE
LB400 exhibits lower activity toward 2,3,4′-trichlorobiphenyl than does BphAE
II9.
BphAE
II9 and BphAE
LB400 produced dihydro-dihydroxy-tetrachlorinated as the only metabolite of 2,2′,5,5′-tetrachlorobiphenyl, and on the basis of the peak area, the amounts produced were approximately the same for both enzymes (
Fig. 5B). Consistent with data from a previous report (
43), the dihydro-dihydroxy-tetrachlorobiphenyl produced from 2,2′,5,5′-tetrachlorobiphenyl by BphAE
II9 and BphAE
LB400 must have been 3,4-dihydro-3,4-dihydroxy-2,2′,5,5′-tetrachlorobiphenyl, and GC-MS analysis showed that the metabolite produced by BphAE
II9 was the same as the one produced by BphAE
LB400.
The bulkier 2,2′,5,5′-tetrachlorobiphenyl docked at the active sites of BphAE
LB400, BphAE
II9, and BphAE
B356 in distinct locations and orientations (
Fig. 7B). Once again, the deepest penetration was found for BphAE
LB400, where docking of 2,2′,5,5′-tetrachlorobiphenyl gave a single pose competent for dioxygenation at positions C-3 and C-4, consistent with available biochemical data (
Table 2). C-3 and C-4 are almost equidistant from the key water molecule at 3.3 Å and from the Fe
2+ ion at 4.2 Å (see Fig. S3B in the supplemental material). The C-5–chlorine faces toward His323, the C-2–chlorine points toward Phe384, and the proximal ring of docked 2,2′,5,5′-tetrachlorobiphenyl is nearly in the same orientation as that observed for the proximal ring of biphenyl in its complex with BphAE
LB400.
In the active site of BphAE
II9, 2,2′,5,5′-tetrachlorobiphenyl docked in multiple poses that would allow
ortho-meta dioxygenation across the C-5—C-6 or C-2—C-3 bonds as well as
meta-para dioxygenation across the C-3—C-4 bond. Since dechlorinated products were not observed in the biochemical analysis, poses consistent with
ortho-meta dioxygenation were ruled out for the purposes of this study. The orientation most competent for
meta-para dioxygenation, consistent with the biochemical data, was chosen for further analysis, although its docking metrics were lower. In this pose, as shown in
Fig. 7A, carbons 3 and 4 are 3.0 Å and 4.2 Å away from the water molecule and 4.1 Å and 5.4 Å away from the Fe
2+ ion, respectively. Unlike BphAE
LB400, the chlorine at the 2-C position of 2,2′,5,5′-tetrachlorobiphenyl is positioned toward His239 in BphAE
II9. The chlorine at the 5-C position faces toward His323 in BphAE
II9 but in a different orientation than that in BphAE
LB400.
Interestingly, the docking experiments with the two enzymes showed very different preferred poses for the binding of this substrate. The binding pose of 2,2′,5,5′-tetrachlorobiphenyl in BphAEII9 would likely generate a steric conflict between the substrate and Phe336 if adopted in BphAELB400. The pose adopted in BphAELB400 likewise appears to be inaccessible to BphAEII9, as it would create a steric conflict with Phe378. However, as discussed above, Phe378 is less constrained in BphAEII9 than in BphAELB400. Therefore, Phe378II9 may be more flexible in reality than was allowed in our docking experiments, which may explain the differences between the results of the biochemical and docking experiments.
For BphAE
B356, the best-scoring pose places 2,2′,5,5′-tetrachlorobiphenyl much farther from the water molecule (at distances of 4.1 Å and 4.3 Å for C-3 and C-4, respectively) and Fe
2+ (at a distance of 5.6 Å for both carbons) than in BphAE
LB400 and BphAE
II9 (
Table 2; see also Fig. S3A in the supplemental material). C-3 and C-4 are closest to the Fe
2+ ion, but the increased distances reflect a binding mode that is likely incompetent for the reaction, consistent with biochemical data that show that BphAE
B356 has poor activity toward 2,2′,5,5′-tetrachlorobiphenyl.
BphAE
II9 retains activity toward 2,2′,5,5′-tetrachlorobiphenyl, while BphAE
B356 shows very little activity toward this particular congener. This apparent contradiction is resolved when one considers the role of an adjacent residue, Thr375
B356, in the position of Phe376
B356. Thr375
B356 Oγ1 makes a hydrogen bond to the backbone carbonyl of Gln371
B356. This hydrogen bond forces Phe376
B356 into a position closer to the Fe
2+ ion (
Fig. 8). The residue corresponding to Thr375
B356 is Asn377
LB400/II9. Asn377
LB400/II9 also forms a hydrogen bond with the backbone carbonyl of His373
LB400/II9 (corresponding to Gln371
B356) and with the backbone carbonyl of Val287
LB400/II9 (
Fig. 8). The second interaction (between residues Thr375
B356 and Gln371
B356) is missing in BphAE
B356, appears to anchor Phe378
II9, restricts how far toward the Fe
2+ ion the residue can move, and thus seems responsible for preserving BphAE
II9's activity toward 2,2′,5,5′-tetrachlorobiphenyl.
Based on this analysis, one would expect that if Thr375
B356 was mutated to an Asn residue, BphAE
B356 could potentially gain the ability to oxygenate 2,2′,5,5′-tetrachlorobiphenyl. Indeed, BphAE
KF707, from a strain that exhibits no activity toward 2,2′,5,5′-tetrachlorobiphenyl, gains the ability to oxygenate this particular congener by a Thr376Asn substitution (corresponding to Thr375
B356 and Asn377
LB400/II9) (
44,
45). Likewise, the reverse is true: when region III of BphAE
LB400 is replaced by region III of BphAE
KF707 along with mutation of Asn377
LB400 to Thr (as in BphAE
KF707/B356), the enzyme loses its activity against
ortho-substituted congeners (
46). However, it was also reported previously by Mondello et al. that the Asn377
LB400Thr substitution alone (i.e., without swapping of region III) did not affect the substrate specificity of the mutated enzyme (
23). Thus, neither the single Asn377Thr substitution nor the swapping of region III alone appears to restructure BphAE
LB400 enough to disrupt its ability to act against
ortho-substituted congeners.
As reported previously (
37), on the basis of product formation, the
Km and
kcat values of BphAE
B356 for DDT were 174 ± 5 μM and 0.15 ± 0.08 s
−1, respectively. During the current work, we found that the activity of BphAE
LB400 toward DDT was too poor to calculate accurate steady-state kinetic values. On the other hand, BphAE
II9 was able to metabolize DDT, exhibiting
Km and
kcat values of 82 ± 2 μM and 0.18 ± 0.04 s
−1, respectively. The GC-MS chromatogram produced from BphAE
B356 against DDT shows two metabolites identified as stereoisomers of 1,1,1-trichloro-2,(4-chlorophenyl-2,3-dihydroxy-4,6-cyclohexadiene)-2-(4′-chlorophenyl)ethane, as shown in
Fig. 5C. BphAE
II9 produced the same two stereoisomers but in an inverse ratio to that of BphAE
B356 (
Fig. 5C).
The active sites of BphAE
II9, BphAE
LB400, and BphAE
B356 all accommodated DDT with corresponding 4-Cl-phenyl rings in similar locations. For the proximal ring, the distances between C-4 atoms are 1.1 Å to 1.4 Å, and the orientations of the rings are also similar (
Fig. 9B). Nevertheless, the placement of the trichloromethyl group distinguishes the BphAE
LB400 complex from the others, and this complex had the poorest docking metrics (
Table 2).
For BphAE
II9, the top-ranked pose is consistent with dioxygenation across the C-2—C-3 bond, and the trichloromethyl group is directed away from the Fe
2+ ion toward Gly321 (
Fig. 9A). The distances from C-2 and C-3 of the proximal ring to the water molecule that marks the dioxygen binding site are 2.6 Å and 2.3 Å, respectively, and the distances to Fe
2+ are 4.5 Å and 4.1 Å, respectively. The best pose for BphAE
B356 is similar (see Fig. S4A in the supplemental material), but the C-2 and C-3 atoms are more distant from the water atom, 3.3 Å and 3.1 Å, respectively, and from the Fe
2+ ion, 4.9 Å and 4.3 Å, respectively. The trichloromethyl group points in a similar direction toward Gly319 (aligns with Gly321 of BphAE
II9).
In the best pose for the BphAELB400 complex, the C-2 and C-3 atoms lie 3.4 Å from the water molecule and 4.0 Å and 4.1 Å from the Fe2+ ion, respectively; however, the trichloromethyl group lies near Phe378 and the Fe2+ ion, as shown in Fig. S4B in the supplemental material, such that the shortest Cl-Fe2+ distance is 4.1 Å. Moreover, although the C-water and C-Fe2+ distances are similar in the three docked complexes, the geometric relationships between the ring, the water molecule, and the Fe2+ ion in BphAELB400 differ markedly. In the BphAEII9 and BphAEB356 complexes, the plane defined by C-2, C-3, and the Fe2+ ion is nearly orthogonal to the plane of the ring, and the water molecule is only ∼0.8 Å out of the C-2–C-3–Fe2+ plane but >2.0 Å out of the plane of the ring. In contrast, in the BphAELB400 complex, the angle between the planes is ∼45°, and the water molecule is ∼1.5 Å away from the C-2–C-3–Fe2+ plane and within 1.0 Å of the plane of the ring.
Therefore, our study demonstrates that BphAEII9 and BphAEB356 have similar levels of activity toward DDT. This finding is further supported by docking studies that showed that BphAEB356 and BphAEII9 have preferred binding poses for DDT that are quite similar. Unlike the other two enzymes, BphAELB400 metabolized DDT very poorly, and the docking experiments generated a binding pose that was completely different from that of the other two enzymes. On the basis of biochemical data, the preferred binding pose modeled in BphAELB400 is not a productive orientation. However, if DDT was oriented in BphAELB400 in a way similar to what was observed for the docking with BphAEII9, the substrate would interfere with Gly321-Gln322 and Phe336. This interference is not present in BphAEII9 due to the substitution of Phe336 with Ile and the elimination of a hydrogen bond that constrains the position of Gly321-Gln322.
In conclusion, our study of the structure of BphAEII9 suggests an enzyme with increased flexibility in and around its active site compared to its two parent enzymes. This increased flexibility would allow it to better accommodate a wide variety of potential substrates, and this is reflected in its enhanced substrate profile compared to its parent enzymes.