Active-site engineering of biphenyl dioxygenase: effect of substituted amino acids on substrate specificity and regiospecificity

Engineering Burkholderia xenovorans LB400 BphA through Site-Directed Mutagenesis at Position 283

Biphenyl dioxygenase (BPDO), which is a Rieske-type oxygenase (RO), catalyzes the initial dioxygenation of biphenyl and some polychlorinated biphenyls (PCBs). In order to enhance the degradation ability of BPDO in terms of a broader substrate range, the BphAE_S283M, BphAE_p4-S283M, and BphAE_RR41-S283M variants were created from the parent enzymes BphAE_LB400, BphAE_p4, and BphAE_RR41, respectively, by a substitution at one residue, Ser283Met. The results of steady-state kinetic parameters show that for biphenyl, the k _cat/K_m values of BphAE_S283M, BphAE_p4-S283M, and BphAE_RR41-S283M were significantly increased compared to those of their parent enzymes. Meanwhile, we determined the steady-state kinetics of BphAEs toward highly chlorinated biphenyls. The results suggested that the Ser283Met substitution enhanced the catalytic activity of BphAEs toward 2,3',4,4'-tetrachlorobiphenyl (2,3',4,4'-CB), 2,2',6,6'-tetrachlorobiphenyl (2,2',6,6'-CB), and 2,3',4,4',5-pentachlorobiphenyl (2,3',4,4',5-CB). We compared the catalytic reactions of BphAE_LB400 and its variants toward 2,2'-dichlorobiphenyl (2,2'-CB), 2,5-dichlorobiphenyl (2,5-CB), and 2,6-dichlorobiphenyl (2,6-CB). The biochemical data indicate that the Ser283Met substitution alters the orientation of the substrate inside the catalytic site and, thereby, its site of hydroxylation, and this was confirmed by docking experiments. We also assessed the substrate ranges of BphAE_LB400 and its variants with degradation activity. BphAE_S283M and BphAE_p4-S283M were clearly improved in oxidizing some of the 3-6-chlorinated biphenyls, which are generally very poorly oxidized by most dioxygenases. Collectively, the present work showed a significant effect of mutation Ser283Met on substrate specificity/regiospecificity in BPDO. These will certainly be meaningful elements for understanding the effect of the residue corresponding to position 283 in other Rieske oxygenase enzymes.IMPORTANCE The segment from positions 280 to 283 in BphAEs is located at the entrance of the catalytic pocket, and it shows variation in conformation. In previous works, results have suggested but never proved that residue Ser283 of BphAE_LB400 might play a role in substrate specificity. In the present paper, we found that the Ser283Met substitution significantly increased the specificity of the reaction of BphAE toward biphenyl, 2,3',4,4'-CB, 2,2',6,6'-CB, and 2,3',4,4',5-CB. Meanwhile, the Ser283Met substitution altered the regiospecificity of BphAE toward 2,2'-dichlorobiphenyl and 2,6-dichlorobiphenyl. Additionally, this substitution extended the range of PCBs metabolized by the mutated BphAE. BphAE_S283M and BphAE_p4-S283M were clearly improved in oxidizing some of the more highly chlorinated biphenyls (3 to 6 chlorines), which are generally very poorly oxidized by most dioxygenases. We used modeled and docked enzymes to identify some of the structural features that explain the new properties of the mutant enzymes. Altogether, the results of this study provide better insights into the mechanisms by which BPDO evolves to change and/or expand its substrate range and its regiospecificity.

Estructura y función de las oxigenasas tipo Rieske/mononuclear

Las oxigenasas Rieske/mononuclear son un grupo de metaloenzimas que catalizan la oxidación de una variedad de compuestos, destaca su participación en la degradación de compuestos xenobióticos contaminantes; estas enzimas también participan en la biosíntesis de algunos compuestos de interés comercial. Poseen una amplia especificidad por el sustrato, convirtiéndolas en un grupo de enzimas con un alto potencial de aplicación en procesos biotecnológicos que hasta el momento no ha sido explotado. La presente revisión aborda aspectos generales acerca de la función y estructura de este importante grupo de enzimas. 

Relationship between the binding free energy and PCBs' migration, persistence, toxicity and bioaccumulation using a combination of the molecular docking method and 3D-QSAR

The molecular docking method was used to calculate the binding free energies between biphenyl dioxygenase and 209 polychlorinated biphenyl (PCB) congeners. The relationships between the calculated binding free energies and migration (octanol-air partition coefficients, KOA), persistence (half-life, t1/2), toxicity (half maximal inhibitory concentration, IC50), and bioaccumulation (bioconcentration factor, BCF) values for the PCBs were used to gain insight into the degradation of PCBs in the presence of biphenyl dioxygenase. The relationships between the calculated binding free energies and the molecular weights, KOA, BCF, and t1/2 values for the PCBs were statistically significant (P < 0.01), whereas the relationship between the calculated binding free energies and the IC50 for the PCBs was not statistically significant (P > 0.05). The electrostatic field, derived from three-dimensional quantitative structure-activity relationship studies, was a primary factor governing the binding free energy, which agreed with literature findings for KOA, t1/2, and BCF. Comparative molecular field analysis and comparative molecular similarity indices analysis contour maps showed that the binding free energies, KOA, t1/2, and BCF values for the PCBs decreased simultaneously when substituents with electropositive groups at the 3-position or electronegative groups at the 3'-position were introduced. This indicated the binding free energy was correlated with the persistent organic pollutant characteristics of PCBs. Furthermore, low binding free energies improved the degradation of the PCBs and simultaneously decreased the KOA, t1/2, and BCF values, thereby reducing the persistent organic pollutant characteristics of PCBs in the environment. These results are expected to be beneficial in providing a theoretical foundation for further elucidation of the degradation and molecular modification of PCBs.

Structural basis of the divergent oxygenation reactions catalyzed by the rieske nonheme iron oxygenase carbazole 1,9a-dioxygenase

Carbazole 1,9a-dioxygenase (CARDO), a Rieske nonheme iron oxygenase (RO), is a three-component system composed of a terminal oxygenase (Oxy), ferredoxin, and a ferredoxin reductase. Oxy has angular dioxygenation activity against carbazole. Previously, site-directed mutagenesis of the Oxy-encoding gene from Janthinobacterium sp. strain J3 generated the I262V, F275W, Q282N, and Q282Y Oxy derivatives, which showed oxygenation capabilities different from those of the wild-type enzyme. To understand the structural features resulting in the different oxidation reactions, we determined the crystal structures of the derivatives, both free and complexed with substrates. The I262V, F275W, and Q282Y derivatives catalyze the lateral dioxygenation of carbazole with higher yields than the wild type. A previous study determined the crystal structure of Oxy complexed with carbazole and revealed that the carbonyl oxygen of Gly178 hydrogen bonds with the imino nitrogen of carbazole. In these derivatives, the carbazole was rotated approximately 15, 25, and 25°, respectively, compared to the wild type, creating space for a water molecule, which hydrogen bonds with the carbonyl oxygen of Gly178 and the imino nitrogen of carbazole. In the crystal structure of the F275W derivative complexed with fluorene, C-9 of fluorene, which corresponds to the imino nitrogen of carbazole, was oriented close to the mutated residue Trp275, which is on the opposite side of the binding pocket from the carbonyl oxygen of Gly178. Our structural analyses demonstrate that the fine-tuning of hydrophobic residues on the surface of the substrate-binding pocket in ROs causes a slight shift in the substrate-binding position that, in turn, favors specific oxygenation reactions toward various substrates.

Retuning Rieske-type oxygenases to expand substrate range

Rieske-type oxygenases are promising biocatalysts for the destruction of persistent pollutants or for the synthesis of fine chemicals. In this work, we explored pathways through which Rieske-type oxygenases evolve to expand their substrate range. BphAE(p4), a variant biphenyl dioxygenase generated from Burkholderia xenovorans LB400 BphAE(LB400) by the double substitution T335A/F336M, and BphAE(RR41), obtained by changing Asn(338), Ile(341), and Leu(409) of BphAE(p4) to Gln(338), Val(341), and Phe(409), metabolize dibenzofuran two and three times faster than BphAE(LB400), respectively. Steady-state kinetic measurements of single- and multiple-substitution mutants of BphAE(LB400) showed that the single T335A and the double N338Q/L409F substitutions contribute significantly to enhanced catalytic activity toward dibenzofuran. Analysis of crystal structures showed that the T335A substitution relieves constraints on a segment lining the catalytic cavity, allowing a significant displacement in response to dibenzofuran binding. The combined N338Q/L409F substitutions alter substrate-induced conformational changes of protein groups involved in subunit assembly and in the chemical steps of the reaction. This suggests a responsive induced fit mechanism that retunes the alignment of protein atoms involved in the chemical steps of the reaction. These enzymes can thus expand their substrate range through mutations that alter the constraints or plasticity of the catalytic cavity to accommodate new substrates or that alter the induced fit mechanism required to achieve proper alignment of reaction-critical atoms or groups.

Expression in Escherichia coli of biphenyl 2,3-dioxygenase genes from a Gram-positive polychlorinated biphenyl degrader, Rhodococcus jostii RHA1

Rhodococcus jostii RHA1 is a polychlorinated biphenyl degrader. Multi-component biphenyl 2,3-dioxygenase (BphA) genes of RHA1 encode large and small subunits of oxygenase component and ferredoxin and reductase components. They did not express enzyme activity in Escherichia coli. To obtain BphA activity in E. coli, hybrid BphA gene derivatives were constructed by replacing ferredoxin and/or reductase component genes of RHA1 with those of Pseudomonas pseudoalcaligenes KF707. The results obtained indicate a lack of catalytic activity of the RHA1 ferredoxin component gene, bphAc in E. coli. To determine the cause of inability of RHA1 bphAc to express in E. coli, the bphAc gene was introduced into Rosetta (DE3) pLacI, which has extra tRNA genes for rare codons in E. coli. The resulting strain abundantly produced the bphAc product, and showed activity. These results suggest that codon usage bias is involved in inability of RHA1 bphAc to express its catalytic activity in E. coli.

Toward protein engineering for phytoremediation: possibilities and challenges

The combination of rational protein engineering and directed evolution techniques allow for the redesign of enzymes with tailored properties for use in environmental remediation. This review summarizes current molecular methods for either altering or improving protein function and highlights examples of how these methods can address bioremediation problems. Although much of the protein engineering applied to environmental clean-up employs microbial systems, there is great potential for and significant challenges to translating these approaches to plant systems for phytoremediation purposes. Protein engineering technologies combined with genomic information and metabolic engineering strategies hold promise for the design of plants and microbes to remediate organic and inorganic pollutants.

Structural insight into the expanded PCB-degrading abilities of a biphenyl dioxygenase obtained by directed evolution

The biphenyl dioxygenase of Burkholderia xenovorans LB400 is a multicomponent Rieske-type oxygenase that catalyzes the dihydroxylation of biphenyl and many polychlorinated biphenyls (PCBs). The structural bases for the substrate specificity of the enzyme's oxygenase component (BphAE(LB400)) are largely unknown. BphAE(p4), a variant previously obtained through directed evolution, transforms several chlorobiphenyls, including 2,6-dichlorobiphenyl, more efficiently than BphAE(LB400), yet differs from the parent oxygenase at only two positions: T335A/F336M. Here, we compare the structures of BphAE(LB400) and BphAE(p4) and examine the biochemical properties of two BphAE(LB400) variants with single substitutions, T335A or F336M. Our data show that residue 336 contacts the biphenyl and influences the regiospecificity of the reaction, but does not enhance the enzyme's reactivity toward 2,6-dichlorobiphenyl. By contrast, residue 335 does not contact biphenyl but contributes significantly to expansion of the enzyme's substrate range. Crystal structures indicate that Thr335 imposes constraints through hydrogen bonds and nonbonded contacts to the segment from Val320 to Gln322. These contacts are lost when Thr is replaced by Ala, relieving intramolecular constraints and allowing for significant movement of this segment during binding of 2,6-dichlorobiphenyl, which increases the space available to accommodate the doubly ortho-chlorinated congener 2,6-dichlorobiphenyl. This study provides important insight about how Rieske-type oxygenases can expand substrate range through mutations that increase the plasticity and/or mobility of protein segments lining the catalytic cavity.

Hybrid pseudomonads engineered by two-step homologous recombination acquire novel degradation abilities toward aromatics and polychlorinated biphenyls

Pseudomonas pseudoalcaligenes KF707 possesses a chromosomally encoded bph gene cluster responsible for the catabolism of biphenyl and polychlorinated biphenyls. Previously, we constructed chimeric versions of the bphA1 gene, which encodes a large subunit of biphenyl dioxygenase, by using DNA shuffling between bphA1 genes from P. pseudoalcaligenes KF707 and Burkholderia xenovorans LB400. In this study, we demonstrate replacement of the bphA1 gene with chimeric bphA1 sequence within the chromosomal bph gene cluster by two-step homologous recombination. Notably, some of the hybrid strains acquired enhanced and/or expanded degradation capabilities for specific aromatic compounds, including single aromatic hydrocarbons and polychlorinated biphenyls.

Transgenic plants to improve rhizoremediation of polychlorinated biphenyls (PCBs)

Recent investigations have shown that the three components of the biphenyl dioxygenase and the 2,3-dihydroxybiphenyl dioxygenase can be produced actively in transgenic plants. Both enzymes catalyze critical steps of the bacterial polychlorinated biphenyl (PCB) degrading pathway. On the basis of these observations, optimized plant-microbe bioremediation processes in which transgenic plants would initiate PCB metabolism and release the metabolites for further degradation by rhizobacteria has been proposed. Since this is still a relatively new approach for PCB remediation, its successful application will require efforts first, to engineer improved PCB-degrading enzymes; second, to co-ordinately express these enzymes' components in plants; and third, to better understand the mechanisms by which plants and rhizobacteria interact to degrade organic pollutants.

Family shuffling of soil DNA to change the regiospecificity of Burkholderia xenovorans LB400 biphenyl dioxygenase

Previous work has shown that the C-terminal portion of BphA, especially two amino acid segments designated region III and region IV, influence the regiospecificity of the biphenyl dioxygenase (BPDO) toward 2,2'-dichlorobiphenyl (2,2'-CB). In this work, we evolved BPDO by shuffling bphA genes amplified from polychlorinated biphenyl-contaminated soil DNA. Sets of approximately 1-kb DNA fragments were amplified with degenerate primers designed to amplify the C-terminal portion of bphA. These fragments were shuffled, and the resulting library was used to replace the corresponding fragment of Burkholderia xenovorans LB400 bphA. Variants were screened for their ability to oxygenate 2,2'-CB onto carbons 5 and 6, which are positions that LB400 BPDO is unable to attack. Variants S100, S149, and S151 were obtained and exhibited this feature. Variant S100 BPDO produced exclusively cis-5,6-dihydro-5,6-dihydroxy-2,2'-dichlorobiphenyl from 2,2'-CB. Moreover, unlike LB400 BPDO, S100 BphA catalyzed the oxygenation of 2,2',3,3'-tetrachlorobiphenyl onto carbons 5 and 6 exclusively and it was unable to oxygenate 2,2',5,5'-tetrachlorobiphenyl. Based on oxygen consumption measurements, variant S100 oxygenated 2,2'-CB at a rate of 16 +/- 1 nmol min(-1) per nmol enzyme, which was similar to the value observed for LB400 BPDO. cis-5,6-Dihydro-5,6-dihydroxy-2,2'-dichlorobiphenyl was further oxidized by 2,3-dihydro-2,3-dihydroxybiphenyl dehydrogenase (BphB) and 2,3-dihydroxybiphenyl dioxygenase (BphC). Variant S100 was, in addition, able to oxygenate benzene, toluene, and ethyl benzene. Sequence analysis identified amino acid residues M237 S238 and S283 outside regions III and IV that influence the activity toward doubly ortho-substituted chlorobiphenyls.