|
Advertisement | ||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
J. Biol. Chem., Vol. 280, Issue 11, 10683-10694, March 18, 2005
BzdR, a Repressor That Controls the Anaerobic Catabolism of Benzoate in Azoarcus sp. CIB, Is the First Member of a New Subfamily of Transcriptional Regulators*![]() ![]() ![]() ![]() ![]() ![]() ||**
From the
Received for publication, October 29, 2004 , and in revised form, December 13, 2004.
In this work, we have studied the transcriptional regulation of the bzd operon involved in the anaerobic catabolism of benzoate in the denitrifying Azoarcus sp. strain CIB. The transcription start site of the PN promoter running the expression of the bzd catabolic genes was identified. Gel retardation assays and PN::lacZ translational fusion experiments performed both in Azoarcus sp. CIB and Escherichia coli cells have shown that bzdR encodes a specific repressor that controls the inducible expression of the adjacent bzd catabolic operon, being the first intermediate of the catabolic pathway (i.e. benzoyl-CoA, the actual inducer molecule). This is the first report of a transcriptional repressor and a CoA-derived aromatic inducer controlling gene expression in the anaerobic catabolism of aromatic compounds. DNase I footprinting experiments revealed that BzdR protected three regions (operators) at the PN promoter. The three operators contain direct repetitions of a TGCA sequence that forms part of longer palindromic structures. In agreement with the repressor role of BzdR, operator region I spans the transcription initiation site as well as the -10 sequence for recognition of the RNA polymerase. Primary sequence analyses of BzdR showed an unusual modular organization with an N-terminal region homologous to members of the HTH-XRE family of transcriptional regulators and a C-terminal region similar to shikimate kinases. A three-dimensional model of the N-terminal and C-terminal regions of BzdR, generated by comparison with the crystal structures of the SinR regulator from Bacillus subtilis and the shikimate kinase I protein from E. coli, strongly suggests that they contain the helix-turn-helix DNA-binding motif and the benzoyl-CoA binding groove, respectively. The BzdR protein constitutes, therefore, the prototype of a new subfamily of transcriptional regulators.
In the last 3 decades, some insights have been added to the knowledge of the anaerobic catabolism of aromatic acids. Most of the studies have focused on the biochemical steps that lead to the mineralization of the aromatic ring in a few microorganisms such as certain strains of Rhodopseudomonas (photosynthetic -Proteobacteria), Thauera, and Azoarcus (denitrifying -Proteobacteria), in some Fe(III) and sulfate reducers, and in fermentative bacteria (14).
Benzoate is a good model compound by which to study the anaerobic catabolism of aromatic compounds (1). A general strategy used by microbes to degrade anaerobically benzoate is its activation to benzoyl-CoA that is then degraded to central biosynthetic intermediates through a series of reactions that involve aromatic ring reduction, The bzd cluster of Azoarcus sp. CIB codes for the set of proteins involved in the anaerobic catabolism of benzoate (2) (Fig. 1A). The bzd catabolic genes are organized in a single operon (bzdNOPQMSTUVWXYZA) encoding the enzymes that activate, dearomatize, and cleave the aromatic ring of benzoate yielding CoA derivative molecules that are channeled into the central metabolism of the cell (2). The bzd catabolic operon is driven by the PN promoter (Fig. 1A). Upstream of the bzd operon is located the bzdR gene that is expressed from the PR promoter (Fig. 1A) (2). It has been suggested previously that the bzdR gene might encode the specific transcriptional regulator (BzdR) of the bzd cluster (2).
By using different genetic and biochemical approaches, we show in this work that bzdR codes for a specific regulator of the bzd catabolic genes, being the first example of a benzoyl-CoA-dependent transcriptional repressor that is described in anaerobic catabolism of aromatic compounds. Moreover, primary sequence analyses and three-dimensional modeling of BzdR reveal an unusual modular structure that allows us to consider BzdR as the first member of a new subfamily of transcriptional regulators.
Strains, Plasmids, and Growth ConditionsThe Escherichia coli and Azoarcus strains as well as the plasmids used in this work are listed in Table I. To construct plasmid pCK01BzdR, a 1.6-kb EcoRI/BamHI DNA fragment containing the bzdR gene was PCR-amplified from plasmid pECOR7 (Table I) by using oligonucleotides 5R (5'-GGAATTCAGGTTCATTCGGCAGCGAGACA-3'; an engineered EcoRI site is underlined) and 3R (5'-GCGGATCCTTCGACATACTCGGCCCGAGCGGTATCTAC-3'; an engineered BamHI is underlined); the resulting fragments were then cloned into EcoRI/BamHI double-digested pCK01 plasmid (Table I). Plasmid pQE32-His6BzdR was constructed by cloning into double-digested pQE32 plasmid a 1-kb BamHI/HindIII fragment harboring the bzdR gene that was PCR-amplified from the pECOR7 plasmid by using oligonucleotides 5HISReg (5'-CGGGATCCTTTCCAACGATGAGAACTCATCAC-3'; an engineered BamHI is underlined) and 3HISReg (5'-GGGAAGCTTTCAGCGTGCCAGGACTTCGAGG-3'; an engineered HindIII is underlined). Unless indicated otherwise, E. coli cells were grown at 37 °C in Luria-Bertani (LB) medium (15). When required, E. coli cells were grown anaerobically in M63 minimal medium (15) at 30 °C using the corresponding necessary nutritional supplements and 20 mM glycerol, as carbon source, and 10 mM nitrate, as terminal electron acceptor. Azoarcus strains were grown anaerobically in MC medium as described previously (2). Where appropriate, antibiotics were added at the following concentrations: ampicillin (100 µg/ml), chloramphenicol (30 µg/ml), gentamycin (7.5 µg/ml), kanamycin (50 µg/ml), and rifampicin (50 µg/ml).
Molecular Biology TechniquesRecombinant DNA techniques were carried out by published methods (8). Plasmid DNA was prepared with a High Pure plasmid isolation kit (Roche Applied Science). DNA fragments were purified with Gene-Clean Turbo (Q-BIOgene). Oligonucleotides were synthesized on an Oligo-1000 M nucleotide synthesizer (Beckman Instruments, Inc.). All cloned inserts and DNA fragments were confirmed by DNA sequencing through an ABI Prism 377 automated DNA sequencer (Applied Biosystems Inc.). Transformation of E. coli cells was carried out by using the RbCl method or by electroporation (Gene Pulser; Bio-Rad) (8). Plasmids were transferred from E. coli S171 pir (donor strain) into Azoarcus sp. recipient strains by biparental filter mating as described previously (2). Proteins were analyzed by SDS-PAGE as described previously (8). The protein concentration in cell extracts was determined by the method of Bradford (16) by using bovine serum albumin as the standard.
N-Terminal Amino Acid Sequence DeterminationThe N-terminal sequence of BzdR was determined by Edman degradation with a 477A automated protein sequencer (Applied Biosystems Inc.). A crude extract of E. coli DH5
Sequence Data AnalysesThe amino acid sequences of open reading frames were compared with those present in finished and unfinished microbial genome data bases using the TBLAST algorithm (18) at the National Center for Biotechnology Information server (available on the World Wide Web at www.ncbi.nlm.nih.gov/blast/blast.cgi). Multiple protein sequence alignments were made with the ClustalW (19) program at the INFOBIOGEN server (available on the World Wide Web at www.infobiogen.fr/services/menuserv.html). Phylogenetic analysis of the BzdR-like proteins was carried out according to the neighbor-joining method of the PHYLIP program (20, 21) at the TreeTop-GeneBee server (available on the World Wide Web at www.genebee.msu.su/service/phtree_reduced.html).
Construction of E. coli Strains Harboring Chromosomal PN::lacZ Translational FusionsTwo different 'lacZ translational fusions were accomplished. Thus, a 598-bp DNA fragment spanning from the 3'-end of the bzdR gene to the translation start codon of the bzdN gene was PCR-amplified using plasmid pECOR7 (Table I) as template and oligonucleotides 3REG (5'-GGGGTACCCGTGCATCAATGATCCGGCAAG-3'; an engineered KpnI site is underlined) and 5BZN (5'-GCTCTAGACCCATCGAACTATCTCCTCTGATG-3'; the bzdN start codon is indicated in boldface type, and an engineered XbaI site is underlined) (2). On the other hand, a 2070-bp DNA fragment spanning from the 5'-end of the bzdR gene to the translation start codon of bzdN was PCR-amplified using plasmid pECOR7 as template and oligonucleotides 5REG (5'-GGGGTACCGGTTTCGTCGCAGGTGCTGTCTGGC-3'; an engineered KpnI site is underlined) and 5BZN. The PCR-amplified fragments were digested with KpnI and XbaI restriction endonucleases and ligated to the KpnI/XbaI double-digested pSJ3 promoter probe vector, giving rise to plasmids pSJ3PN (PN::lacZ) and pSJ3RPN (Pr-bzdR/PN::lacZ) (Table I, Fig. 4). The correct lacZ translational fusions were confirmed by nucleotide sequence analysis. Plasmids pUTminiTn5PN and pUTminiTn5RPN were constructed by subcloning the NotI cassettes from pSJ3PN and pSJ3RPN into the mini-Tn5 delivery plasmid pUTminiTn5Km2 (Table I, Fig. 4), and they were transferred from E. coli S171
Construction of Azoarcus sp. Strain CIBdbzdRFor disruption of the bzdR gene through single homologous recombination, a 499-bp internal fragment of bzdR was PCR-amplified by using primers BzdRInt5' (5'-GCGAGCGCAGCGGGGAATG-3') and N-INV-III (5'-GCTTGCGACCACCGACTCC-3'), and it was cloned into the SmaI-digested pK18mob (a mobilizable plasmid that does not replicate in Azoarcus) (Table I). The resulting construct, pK18mobbzdR (Table I), was transferred from E. coli S171 pir (donor strain) into Azoarcus sp. strain CIB (recipient strain) by biparental filter mating (2). Exconjugants harboring the disrupted bzdR gene by insertion of the suicide plasmid, namely Azoarcus sp. strain CIBdbzdR, were isolated aerobically on kanamycin-containing MC medium lacking nitrate and containing 0.4% citrate as the sole carbon source for counterselection of donor cells. The mutant strain was analyzed by PCR to confirm the disruption of the target gene.
Overproduction and Purification of His6-BzdRThe recombinant plasmid pQE32-His6BzdR (Table I) carries the bzdR gene without the ATG start codon and with a His6 tag coding sequence at its 5'-end, under control of the T5 promoter and two lac operator boxes. The His tail added 13 amino acids (MCGSHHHHHHGIL) to the N-terminal end of His6-BzdR-protein (34.842 Da). The His-tagged protein was overproduced in E. coli M15 strain harboring plasmid pREP4 (Table I) that produces the LacI repressor to strictly control gene expression from pQE32 derivatives in the presence of isopropyl-1-thio-
Gel Retardation AssaysA 268-bp DNA fragment spanning from position -174 to +79 with respect to the transcription start site of PN promoter (Fig. 1B) was PCR-amplified from plasmid pSJ3PN (Table I) by using oligonucleotides 5BZN (see above) and RetIII (5'-GGGAATTCCGAGCCTCGCGTTTTACTGC-3'; an engineered EcoRI site is underlined). This fragment was digested with XbaI and EcoRI restriction enzymes, and it was cloned in XbaI/EcoRI double-digested pUC18 vector, giving rise to plasmid pUC18FP4 (Table I). The DNA fragment was then excised from pUC18FP4 by digestion with EcoRI and HincII restriction enzymes, and it was labeled by filling in the overhanging EcoRI-digested end with [
DNase I Footprinting AssaysThe DNA fragment used for DNase I footprinting assays was PCR-amplified from plasmid pECOR7 by using oligonucleotides 5IVTPN (5'-CGGAATTCCGTGCATCAATGATCCGGCAAG-3'), which hybridizes with the 3'-end of bzdR, and 3IVTPN (5'-CGGAATTCCATCGAACTATCTCCTCTGATG-3'; an engineered EcoRI site is underlined, and the bzdN start codon is indicated in boldface type). The amplified DNA fragment was then digested with PvuII and EcoRI restriction enzymes, and the resulting 376-bp substitution was singly 3'-end-labeled by filling in the overhanging EcoRI-digested end with [
Primer Extension AnalysisAzoarcus sp. CIB cells harboring plasmid pBBR2PN (PN::lacZ) were grown anaerobically on benzoate-containing MC medium until the culture reached an A600 of 0.3. Total RNA was isolated by using the RNeasy Mini kit (Qiagen) according to the instructions of the supplier. Primer extension reactions were carried out with the avian myeloblastosis virus reverse transcriptase (Promega) and 10 µg of total RNA as described previously (24), using oligonucleotide 5BZN (which hybridizes with the coding strand of the bzdN gene; see above) labeled at its 5'-end with phage T4 polynucleotide kinase and [ Modeling of BzdRThe threading methods used to obtain the three-dimensional models for the N- and C-terminal regions of BzdR were BIOINBGU (25), FFAS (26), FUGUE (27), HMM (28), and 3D-PSSM (29). These programs provide a list of hits with known three-dimensional structures, which presumably are similar to that of the query protein, and thus can be used as modeling templates (Table II). Sequence alignment results between the BzdR sequences and those of the templates were analyzed with ClustalX (19). Once a modeling template was chosen, a full-atom three-dimensional model was obtained by using the SWISS-MODEL Protein Modeling Server (30, 31). The models were energy-minimized with CNS (32), and their structures were validated using PROCHECK (33) and the ERRAT Protein Verification Server (34). The final models after energy minimization with CNS showed a correct stereochemistry, with 100% of the predicted residues in allowed regions in the Ramachandran plot (not shown). Additionally, the results obtained with ERRAT (34) indicate that the percentage below the 95% confidence limit is 77.3 and 86.3 for N-terminal and C-terminal BzdR regions, respectively. Manual docking of the interaction between BzdR and benzoyl-CoA was accomplished by using the program O (35).
BzdR Is the Prototype of a New Subfamily of Transcriptional RegulatorsDetermination of the N-terminal amino acid sequence (MSNDENSSRL) of the bzdR gene product (see "Experimental Procedures"), suggests that BzdR is a protein of 298 amino acids (33,488 Da). Primary sequence analyses and amino acid sequence comparisons revealed that the BzdR protein exhibits two distinct regions. The N-terminal region (residues 1130) shows a significant similarity with members of the helix-turn-helix (HTH)1-XRE family of transcriptional regulators (SMART release; available on the World Wide Web at smart.embl-heidelberg.de/smart) that includes more than 1300 proteins from eukaryota, archaea, bacteriophages, and bacteria with an HTH DNA-binding motif similar to that of the well characterized Cro protein of phage (36, 37). The highest amino acid sequence identity (34%) was found with SinR (Fig. 2), a 14-kDa multifunctional protein that, besides repressing sporulation, regulates genes for motility, alkaline protease expression (aprE), and competence development in Bacillus subtilis (38). The C-terminal region of BzdR (residues 131298) is homologous to shikimate kinases (PFAM accession number PF01202), enzymes that catalyze the conversion of shikimate to shikimate 3-phosphate using ATP as co-substrate, and it shows a sequence identity of 23% with the E. coli shikimate kinase I (aroK gene product) (Fig. 2) (39).
BzdR-like proteins are found encoded in the genomes of other denitrifying bacteria such as T. aromatica (51% identity; AAN32624 [GenBank] , A. evansii (50% identity; AAN39374 [GenBank] , Burkholderia fungorum LB400 (two proteins with 51 and 48% identity; ZP_00279512 and ZP_00283918, respectively), Ralstonia metallidurans (53% identity; ZP_00274120), Ralstonia eutropha (53% identity; ZP_001700685), M. magnetotacticum MS-1 (two proteins with 47 and 44% identity; ZP_00056533 and ZP_00052337, respectively), Bradyrhizobium japonicum (47% identity; NP_767716 [GenBank] ), and the photosynthetic bacterium R. palustris (45% identity; CAE30227 [GenBank] (Fig. 2). With the exception of the BzdR-like protein of M. magnetotacticum (ZP_00056533), there is a good correlation between the taxonomical position of the organism and the level of identity among BzdR homologues (Fig. 2). Thus, BzdR-like proteins from and -Proteobacteria cluster as different branches of the phylogenetic tree (data not shown). Whereas in most of the bacteria, the bzdR-like genes are in the vicinity of or integrated into the aerobic benzoate degradation cluster (40), Azoarcus sp. CIB is the only case described so far with the bzdR gene located in the anaerobic benzoate degradation cluster. All of these data strongly suggest that BzdR is the specific regulator that controls expression of the bzd catabolic operon in Azoarcus sp. CIB. Since the modular organization of BzdR has not been reported so far in other regulators of the HTH-XRE family, BzdR might constitute the prototype of a new subfamily of transcriptional regulators.
BzdR Is the Specific Repressor of the bzd Catabolic Operon To determine whether bzdR is indeed involved in the regulation of the anaerobic benzoate degradation, a bzdR mutant of Azoarcus sp. strain CIB was constructed by gene disruption as described under "Experimental Procedures." The mutant, Azoarcus sp. CIBdbzdR (Table I), grew as the wild-type Azoarcus sp. CIB strain in benzoate-containing minimal medium (data not shown), which suggests that the BzdR protein, if playing any regulatory role, is not an activator but a repressor of the bzd catabolic operon. When we measured expression of the last gene of the bzd operon (i.e. the bzdA gene encoding benzoate-CoA ligase) (Fig. 1A), we observed that whereas Azoarcus sp. CIB cells showed a benzoate-inducible expression of bzdA, a high level of benzoate-CoA ligase activity was detected in Azoarcus sp. CIBdbzdR cells growing either in the presence of benzoate (inducing conditions) or pyruvate (noninducing conditions) (Fig. 3A). This result indicates a constitutive expression of bzdA in Azoarcus CIB sp. cells that lack the bzdR gene, which is in agreement with BzdR being a putative repressor of the bzd catabolic operon. To confirm that BzdR controls the expression of the bzd catabolic operon at the level of the PN promoter, plasmid pBBR5PN, which harbors a PN::lacZ translational fusion (Table I), was introduced both in wild-type Azoarcus sp. CIB cells and in Azoarcus sp. CIBdbzdR cells. Whereas a significant
To further investigate the regulatory role of the BzdR protein on the expression of the PN promoter, E. coli AFMC cells were transformed with plasmids pSJ3PN (PN::lacZ), pSJ3RPN (PR-bzdR/PN::lacZ), which harbors the bzdR gene under control of its own promoter located adjacent to the PN::lacZ reporter fusion (Fig. 4), and pCK01BzdR (Plac::bzdR), which encodes the bzdR gene under control of the heterologous Plac promoter (Table I). Whereas -galactosidase assays of permeabilized E. coli AFMC (pSJ3PN) cells grown anaerobically in glycerol-containing mineral medium revealed expression of PN (Fig. 5A), the presence of the BzdR protein both in cis (E. coli AFMC (pSJ3RPN) cells) or in trans (E. coli AFMC (pSJ3PN, pCK01BzdR) cells) repressed the expression of the PN::lacZ translational fusion (Fig. 5A). As expected, -galactosidase levels in E. coli AFMC cells harboring plasmid pSJ3PN and the control plasmid pCK01 were similar to those observed in E. coli AFMC (pSJ3PN) cells (Fig. 5A). To analyze faithfully the bzd regulatory system at monocopy dosage, the PN::lacZ and PR-bzdR/PN::lacZ reporter fusions were subcloned as NotI-DNA cassettes within mini-Tn5 vectors, rendering plasmids pUTminiTn5PN and pUTminiTn5RPN, respectively (Fig. 4), which were then used to deliver by transposition the corresponding translational fusions into the chromosome of E. coli AFMC. The presence of a strong T7 phage transcriptional terminator downstream of the lacZ fusion and the orientation of such fusions within the mini-Tn5 elements (Fig. 4) prevented read-through transcription from nearby promoters after insertion into the chromosome of the resulting E. coli AFMCPN and AFMCRPN strains (Table I). Despite the fact that -galactosidase levels of E. coli AFMCPN cells harboring the control plasmid pCK01 were about 1 order of magnitude lower than those observed with E. coli AFMC (pSJ3PN) cells (Fig. 5B) due to the decrease in the genetic dosage, the repressor role of BzdR on the PN promoter was confirmed by monitoring a significant reduction of -galactosidase activity both in E. coli AFMCRPN and E. coli AFMCPN (pCK01BzdR) cells (Fig. 5B).
It should be noted that all transcriptional regulators so far described in the anaerobic catabolism of aromatic compounds (i.e. BadR, AadR, and HbaR controlling degradation of benzoate and p-hydroxybenzoate in R. palustris (6, 7) and TutBC and TdiRS controlling degradation of toluene in Thauera sp. T1 and T. aromatica K172, respectively (41, 42)) were shown (BadR, AadR, and HbaR) or were suggested to be (TutBC and TdiRS) activators of their corresponding promoters. The BzdR protein from Azoarcus sp. CIB constitutes, therefore, the first transcriptional repressor that is shown to control the expression of genes involved in the anaerobic catabolism of aromatic compounds.
Benzoyl-CoA Is the Actual Inducer of the Benzoate Degradation PathwayTo elucidate whether the substrate of the bzd pathway from Azoarcus sp. CIB (i.e. benzoate) was the inducer of the bzd catabolic operon, E. coli AFMCRPN cells were grown anaerobically on minimal medium supplemented with glycerol in the presence of benzoate and the expression of the PN::lacZ fusion was monitored through Although previous reports have suggested that the inducer of the genes involved in anaerobic catabolism of benzoate in R. palustris and T. aromatica might be benzoyl-CoA rather than benzoate (1, 43), no direct experimental evidence supporting this assumption has been provided. Here we present for the first time experimental evidence that benzoyl-CoA is the actual inducer of the pathway in Azoacus sp. CIB. It should be noted that another aromatic CoA derivative, phenylacetyl-CoA, was shown to be the inducer of an aerobic hybrid pathway, the phenylacetic acid degradation pathway, whose intermediates are also CoA-derived compounds (11). Whether transcriptional regulators of both anaerobic and aerobic hybrid pathways for catabolism of aromatic acids have evolved to recognize CoA-derived aromatic compounds rather than the free acids is a hypothesis that requires further confirmation when additional regulatory systems of such type of degradative routes become characterized.
In Vitro Binding of BzdR to the PN PromoterTo study in vitro the interaction of the BzdR protein with the PN promoter, we first mapped the transcription start site of this promoter and overproduced the regulatory protein in recombinant E. coli cells. Primer extension analysis was performed with total RNA isolated from Azoarcus sp. strain CIB cells containing plasmid pBBR2PN (Table I). The transcription initiation site in the PN promoter was mapped 75 nucleotides upstream of the ATG translation initiation codon of the bzdN gene (Fig. 6), showing putative -10 (TAACAT) and -35 (TCAACA) boxes typical of
To overproduce BzdR, we have engineered plasmid pQE32-His6BzdR (Table I) that expresses a His-tagged BzdR protein (see "Experimental Procedures"). As expected, the His6-BzdR protein was a repressor of the PN promoter in E. coli AFMCPN (pQE32-His6BzdR) cells growing anaerobically in glycerol-containing minimal medium (Fig. 7A). Induction by benzoyl-CoA of the His6-BzdR mediated repression of PN was confirmed in E. coli AFMCPN (pQE32-His6BzdR, pCK01BzdA) cells growing anaerobically in glycerol-containing minimal medium in the presence (induction) or absence of benzoate (Fig. 7A). These results indicate that the functionality of His-tagged BzdR is similar to that of the wild-type protein, and, thus, His6-BzdR can be used to study the interaction of this regulator with the cognate PN promoter. The purification of the His6-BzdR protein was carried out through affinity chromatography as described under "Experimental Procedures," and it was checked by SDS-PAGE (Fig. 7B). To demonstrate the interaction of the BzdR regulatory protein with the PN promoter, the purified His6-BzdR protein was subjected to gel retardation assays, using as probe a 268-bp DNA fragment that carries the PN promoter region from position -174 to +79 (PN probe) (Fig. 1B). The His6-BzdR protein was able to retard the migration of the PN probe in a protein concentration-dependent manner (Fig. 8A). Moreover, binding of BzdR to the PN promoter was highly specific, because whereas it was inhibited by adding an unlabeled PN probe to the retardation assays, it was not affected by adding an unlabeled heterologous probe (data not shown). Depending on the amount of His6-BzdR protein used in the retardation assay, two distinct complexes, complex 1 (higher mobility) and complex 2 (lower mobility), were observed (Fig. 8A). Thus, whereas complex I appeared preferentially at low amounts (2.510 nM) of His6-BzdR protein, 5-fold higher amounts of His6-BzdR gave rise to complex II (Fig. 8A). These results might suggest that BzdR binds to different regions at the PN promoter (see below). Gel retardation assays with His6-BzdR were also carried out in the presence of different concentrations of the benzoyl-CoA (Sigma) inducer molecule. As shown in Fig. 8B, increasing concentrations of benzoyl-CoA ranging from 250 µM to 2 mM significantly reduced the amount of bound PN probe, suggesting that benzoyl-CoA inhibits binding of His6-BzdR to the PN promoter. Although a percentage of PN probe (about 20%) remained free when 2 mM benzoate was added to the retardation assay, this percentage increased to 95% when 2 mM benzoyl-CoA was used (Fig. 8B). All of these data taken together are in agreement with the lacZ-reporter fusion experiments reported above, and they confirm BzdR as the repressor of the PN promoter and benzoyl-CoA as the actual inducer molecule. Nevertheless, the molecular mechanism of the benzoyl-CoA-mediated induction of the PN promoter is still unknown, and it will be subject of future research.
Identification of the BzdR Binding Site(s)The purified His6-BzdR protein was also used to determine the BzdR-binding site(s) (operator) in the PN promoter region by DNase I footprinting. As shown in Fig. 9, the His6-BzdR protein protected three different regions: region I (63 bp) spanning from position -32 to +31, region II (21 bp) spanning from position -83 to -63, and region III (21 bp) spanning from position -146 to -126 (Fig. 1B). The three protected regions contain direct repetitions of a sequence, TGCA, that forms a part of longer palindromic structures (Fig. 1B). Whereas the TGCA sequences are separated by 6 nucleotides in regions II and III, the longer region I presents a pair of TGCA sequences separated by 1 nucleotide and another pair of TGCA sequences separated by 15 nucleotides (Fig. 1B). Other members of the HTH-XRE family of transcriptional regulators (see above) such as the SinR regulator from B. subtilis and the Cro and 434 repressors from phages and 434, respectively, also bind to short repeated sequences that, in most cases, are located within palindromic regions that span the promoters (38, 4446). Binding of such regulators to the multiple operators of their cognate promoters generates different protein-DNA complexes in gel retardation assays (38, 47), a pattern of binding that was also observed with BzdR at the PN promoter (Fig. 8). Moreover, binding of BzdR induces changes in the DNA structure of PN as revealed by several phosphodiester bonds that become hypersensitive to DNase I cleavage (Fig. 9). The fact that BzdR-binding region I spans the transcription initiation site as well as the -10 sequence for recognition of the 70-RNA polymerase (Fig. 1B) is in agreement with the observed repressor role of BzdR at the PN promoter. Nevertheless, to precisely understand the molecular basis of the BzdR-mediated repression at the PN promoter, more work needs to be done.
Modeling of BzdR and Docking Simulations of Benzoyl-CoA BindingWith the aim of building three-dimensional models of BzdR, we carried out a -fold recognition approach for the N-terminal region, from residue 25 to 87 (N-BzdR), and the C-terminal region, residues 131298 (C-BzdR), of the protein. The five threading methods used in this work found potential templates with high scores, always above the 95% certainty level (see "Experimental Procedures"). The threading results for N-BzdR revealed that all of the templates found belong to the same SCOP superfamily, the repressor-like DNA-binding domains. Pairwise comparisons of the C atoms of the crystal structures of all of the templates found, as calculated with the program O (35), yield a maximum root mean square deviation of 1.3 and correspond to a five-helix closed bundle folding (Fig. 10A). The template selected for building a three-dimensional model of N-BzdR was SinR of B. subtilis (1b0n) (Table II), since it shares the highest identity. The presence within this fold of an HTH (residues 3858 of BzdR) (Figs. 2 and 10A), a motif found in numerous DNA-binding proteins (51), provides a reasonable structural basis for suggesting N-BzdR as the region directly interacting with the target DNA.
The threading results obtained for C-BzdR also provided consistent results, since all of the hits found share the same SCOP fold, the P-loop-containing nucleoside triphosphate hydrolase fold, almost all of them belonging to the shikimate kinase family (see "Experimental Procedures"). The template used to build a full-atom three-dimensional model of C-BzdR was the shikimate kinase I of E. coli (52, 53) (Fig. 2). Interestingly, the amino acid sequence identity shared between C-BzdR and the selected template (23%) is similar to that observed between the two different isoenzymes of E. coli, shikimate kinases I and II (30% identity) (39). This observation indicates that the protein fold of shikimate kinases is highly versatile, since it can accommodate a wide array of sequences without significant structural departures. As shown in Fig. 10B, C-BzdR presents a canonical mononucleotide-binding fold (54) found in a number of structurally diverse proteins (55), which is constituted by a five-stranded parallel -sheet, strand order 23145, flanked by eight -helices. Within this highly conserved fold, a phosphate-binding loop (P-loop) or Walker-A motif (56, 57) can be observed between 1 and 6 (Fig. 10B). The Walker-A motif (consensus sequence: GLRGAGK(T/S)) is highly conserved in all BzdR orthologues (Fig. 2), and, in fact, the suggested three-dimensional structure of C-BzdR reveals that the environment around the P-loop is also essentially conserved in all of these proteins. Moreover, the strictly conserved Gly present in the Walker B-motif of purine nucleotide-binding proteins (56), usually located around the 3 strand (55), is also present in BzdR (Gly-209) and its orthologues (Fig. 2).
The three-dimensional model of C-BzdR permits us to envisage a putative three-dimensional model for the complex between C-BzdR and the effector molecule benzoyl-CoA. The feasibility of this model is accentuated by the striking structural similarity between both ends of the benzoyl-CoA molecule (namely the ADP and the benzoyl moieties) and the two substrates of shikimate kinases (ATP and shikimate, respectively). Manual docking of the ADP moiety of benzoyl-CoA into the C-BzdR structure is immediate, since the C-BzdR fold has a conserved nucleotide-binding site (Fig. 10C). On the other hand, the pantothenate and
* This work was supported by Comunidad Autónoma de Madrid Grant 07M/0076/2002 and Comisión Interministerial de Ciencia y Tecnología Grants BIO2003-01482 and VEM2003-20075-CO2-02. The costs of publication of this article were defrayed in part by the payment of page charges. This article must therefore be hereby marked "advertisement" in accordance with 18 U.S.C. Section 1734 solely to indicate this fact.
|| Holder of the Ramón y Cajal Program of the Spanish Ministerio de Educación y Ciencia. ** To whom correspondence should be addressed: Dept. de Microbiología Molecular, Centro de Investigaciones Biológicas-CSIC, Ramiro de Maeztu 9, 28040 Madrid, Spain. Tel.: 34-91-8373112; E-mail: mcarmona{at}cib.csic.es.
1 The abbreviations used are: HTH, helix-turn-helix; LB, Luria-Bertani.
The technical work of I. Alonso is greatly appreciated. We thank E. Aporta and J. Varela for help with oligonucleotide synthesis and N-terminal protein sequencing, respectively. The help of A. Díaz, S. Carbajo, M. Cayuela, and G. Porras with sequencing is gratefully acknowledged.
This article has been cited by other articles:
|
| |||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
|
Advertisement | ||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||