JBC

HOME HELP FEEDBACK SUBSCRIPTIONS ARCHIVE SEARCH TABLE OF CONTENTS
 QUICK SEARCH:   [advanced]


     


Originally published In Press as doi:10.1074/jbc.M610838200 on April 10, 2007

J. Biol. Chem., Vol. 282, Issue 22, 16476-16491, June 1, 2007
This Article
Right arrow Abstract Freely available
Right arrow Full Text (PDF)
Right arrow All Versions of this Article:
282/22/16476    most recent
M610838200v1
Right arrow Alert me when this article is cited
Right arrow Alert me if a correction is posted
Right arrow Citation Map
Services
Right arrow Email this article to a friend
Right arrow Similar articles in this journal
Right arrow Similar articles in PubMed
Right arrow Alert me to new issues of the journal
Right arrow Download to citation manager
Right arrow reprints & permissions
Citing Articles
Right arrow Citing Articles via HighWire
Right arrow Citing Articles via Google Scholar
Google Scholar
Right arrow Articles by Lorca, G. L.
Right arrow Articles by Savchenko, A.
Right arrow Search for Related Content
PubMed
Right arrow PubMed Citation
Right arrow Articles by Lorca, G. L.
Right arrow Articles by Savchenko, A.
Social Bookmarking
 Add to CiteULike   Add to Complore   Add to Connotea   Add to Del.icio.us   Add to Digg   Add to Reddit   Add to Technorati  
What's this?

Glyoxylate and Pyruvate Are Antagonistic Effectors of the Escherichia coli IclR Transcriptional Regulator*

Graciela L. Lorca{ddagger}1, Alexandra Ezersky{ddagger}§, Vladimir V. Lunin, John R. Walker||, Svetlana Altamentova, Elena Evdokimova, Masoud Vedadi||, Alexey Bochkarev{ddagger}||, and Alexei Savchenko{ddagger}2

From the {ddagger}Banting and Best Department of Medical Research, Toronto, Ontario M5G 1L6, Canada, §Department of Medical Biophysics, University of Toronto, Toronto, Ontario M5G 1L6, Canada, Ontario Center for Structural Proteomics, Toronto, Ontario M5G 1L6, Canada, and ||Canada Structural Genomics Consortium, Toronto, Ontario M5G 1L6, Canada

Received for publication, November 22, 2006 , and in revised form, March 28, 2007.


    ABSTRACT
 TOP
 ABSTRACT
 INTRODUCTION
 EXPERIMENTAL PROCEDURES
 RESULTS
 DISCUSSION
 REFERENCES
 
The Escherichia coli isocitrate lyase regulator (IclR) regulates the expression of the glyoxylate bypass operon (aceBAK). Founding member of a large family of common fold transcriptional regulators, IclR comprises a DNA binding domain that interacts with the operator sequence and a C-terminal domain (C-IclR) that binds a hitherto unknown small molecule. We screened a chemical library of more than 150 metabolic scaffolds using a high-throughput protein stability assay to identify molecules that bind IclR and then tested the active compounds in in vitro assays of operator binding. Glyoxylate and pyruvate, identified by this method, bound the C-IclR domain with KD values of 0.9 ± 0.2 and 156.2 ± 7.9 µM, as defined by isothermal titration calorimetry. Both compounds altered IclR interactions with operator DNA in electrophoretic mobility shift assays but showed an antagonistic effect. Glyoxylate disrupted the formation of the IclR/operator complex in vitro by favoring the inactive dimeric state of the protein, whereas pyruvate increased the binding of IclR to the aceBAK promoter by stabilizing the active tetrameric form of the protein. Structures of the C-IclR domain alone and in complex with each effector were determined at 2.3 Å, confirming the binding of both molecules in the effector recognition site previously characterized for the other representative of the family, the E. coli AllR regulator. Site-directed mutagenesis demonstrated the importance of hydrophobic patch formed by Met-146, Leu-154, Leu-220, and Leu-143 in interactions with effector molecules. In general, our strategy of combining chemical screens with functional assays and structural studies has uncovered two small molecules with antagonistic effects that regulate the IclR-dependent transcription of the aceBAK operon.


    INTRODUCTION
 TOP
 ABSTRACT
 INTRODUCTION
 EXPERIMENTAL PROCEDURES
 RESULTS
 DISCUSSION
 REFERENCES
 
A classical bacterial repressor, termed IclR, participates in the regulation of the aceBAK operon in Escherichia coli (1, 2). The aceBAK operon encodes for the enzymes of the glyoxylate bypass (3) that are required during growth on acetate since it bypasses the two CO2-evolving steps of the Krebs cycle (4). The expression of these enzymes is induced during growth on minimal medium supplemented with acetate or fatty acids (3) as well as in rich medium as a result of the acetate accumulation during exponential phase. According to early genetic studies (1, 2) IclR represses the expression of aceBAK operon as well as its own gene by binding to specific operator sequences in the promoter region. Two separate operator sequences, termed IclR-boxes, have been identified in the promoter region upstream of the aceB gene (5). The primary operator sequence (IclR box II) has been mapped between -52 and -19 bases of the aceB promoter, and binding of IclR to this site prevents the RNA polymerase-promoter interaction (5). Binding to the second site, IclR box I, which is located between -125 and -99 bases of the aceB promoter, disassembles the open complex through a protein-protein interaction between IclR and the RNA polymerase {alpha}-subunit (5). A third IclR-box was located between +14 and -21 bases of the iclR promoter (6). All three IclR-boxes identified in the promoters of aceB and iclR genes contain near-perfect palindromes. Although not identical, their sequences share a 15-mer consensus motif, 5'-TGGAAATNATTTCCA, which was identified by in vitro random oligonucleotide selection experiments (7).

In growth media containing acetate, IclR is believed to release the DNA as a result of binding a small molecule. Although this molecule is likely related to a metabolite within the glyoxylate bypass pathway, the precise nature of this IclR effector remains nevertheless in doubt. Phosphoenolpyruvate was reported to alter IclR binding to DNA in vitro (8) but recently was shown to have no effect on IclR activity by Yamamoto and Ishihama (5).

IclR is the founding member of a large family (Pfam 01614, COG1414) of sequence-related microbial transcriptional regulators with more than 500 members identified in bacterial and archaeal genomes. Numerous bacterial genomes such as E. coli K12 (9) and Pseudomonas aeruginosa PA01 (10) contain multiple members of this family (8 and 9, respectively) with the largest number (40 members) identified so far in Bordetella bronchiseptica RB50 (11). The characterized members of this family are involved in the regulation of diverse catabolic pathways ranging from the degradation of plant cell polysaccharides in the plant pathogen Erwinia chrysanthemi (12) to the metabolism of aromatic acids in E. coli, Acinetobacter (13, 14), and Pseudomonas (15).

The characteristic structural features for members of this family, gained primarily from the crystal structure of TM0065 from Thermotoga maritime (PDB code 1MKM), are fully consistent with the role of transcription regulators responding to chemical stimuli (16). The TM0065 structure confirmed the presence of N-terminal classic-winged helix-turn-helix DNA binding domain as well as demonstrated that the effector binding C-terminal domain has significant structural homology to the PAS/GAF domain, a known small molecule binding motif (17, 18). The N- and C-terminal domains in the structure are connected by a linker helix, which together with N-terminal domain participates in protein dimerization. In the crystal lattice TM0065 dimers further oligomerize into tetramers, which appear to be the functional unit common for a number of IclR regulators (16). The C-terminal domain retains its structure and small molecule binding function when expressed without the DNA binding domain and was structurally characterized for three E. coli representatives of IclR family, KdgR (PDB code 1YSP), YiaJ (PDB code 1YSQ), and AllR3 (Ref. 20, PDB code 1TF1). Although sharing relatively low sequence identity between them and with previously characterized TM0065, all four structures were very similar, confirming a common structural framework for most of the IclR family members.

AllR is involved in transcriptional regulation of the allantoin catabolism (19) and shares 42% sequence identity with IclR. The AllR C-terminal domain structure was recently determined in the presence and absence of its chemical effector glyoxylate (20), providing the first detailed view of the molecular interactions of an IclR family representative with a small molecule. According to the AllR-glyoxylate complex structure, the location of the effector binding pocket in IclR regulators corresponded to the co-factor-accommodating cavity in PAS/GAF domains; however, the binding mechanisms were not conserved between these structurally similar domains (17, 18, 20).

Of the hundreds of IclR family members, only a few of the sensor chemicals have been identified, which apart from glyoxylate, include 3-(3-hydroxyphenyl) propionic acid for E. coli MhpR (21), 2-keto-2-deoxyglyconate for E. chrysanthemi KdgR (12), p-hydroxybenzoate and protocatechuate for Acinetobacter PobR and PcaU, respectively (13, 22), N-3-oxo-octanoyl homoserine lactone (3OC8HSL) for Agrobacterium tumefaciens AttJ (23), and 1-naphthol and 2,3-dihydroxynaphthalene for Pseudomonas putida TtgV (24). Thus, the effector specificity of hundreds IclR regulators, including the founding member of the family, remains uncertain.

To address these issues we used IclR as a test case and screened it against a large number of intermediates associated with glyoxylate bypass including phosphoenolpyruvate as part of a library of more than 150 metabolites by a high-throughput protein stability assay (25, 26). This methodology takes advantage of the well documented phenomenon of enhanced thermal stability of the protein in the presence of a specifically bound ligand (27). The molecules that had a stabilizing effect were further assessed for the effector role for IclR using biophysical and functional assays.


    EXPERIMENTAL PROCEDURES
 TOP
 ABSTRACT
 INTRODUCTION
 EXPERIMENTAL PROCEDURES
 RESULTS
 DISCUSSION
 REFERENCES
 
DNA Manipulations and Gene Cloning—Standard methods were used for site-directed mutagenesis, chromosomal DNA isolation, restriction enzyme digestion, agarose gel electrophoresis, ligation, and transformation (28). Plasmids were isolated using spin miniprep kits (Qiagen), and PCR products were purified using Qiaquick purification kits (Qiagen). Mutagenesis was performed using the QuikChange site-directed mutagenesis kit (Stratagene).

For protein expression and purification the iclR gene was amplified from wild type E. coli BW25113 chromosomal DNA by PCR. To identify the boundaries of the ligand binding domain of IclR, a multiple sequence alignment of IclR and its homologues was performed. The link between the DNA and ligand binding domains in E. coli IclR was predicted to lie between amino acids Ser 98 and Leu 101. Two fragments corresponding to the 176 and 173 C-terminal amino acids were amplified and tested for expression of the soluble domain as previously described (20). The expression of a fragment of iclR corresponding to amino acids 98–274 yielded a soluble domain of IclR (C-IclR) that was used for further experiments. All amplified fragments were cloned into the NdeI and BamHI sites of a modified form of pET15b (EMD Biosciences) in which a TEV protease cleavage site replaced the thrombin cleavage site and a double stop codon was introduced downstream from the BamHI site. This construct provides for an N-terminal hexahistidine tag separated from the protein by a tobacco etch virus protease recognition site (ENLYFQ{downarrow}GS).

Construction of Deletion Mutants—An IclR deletion mutant was generated by the methods described by Datsenko and Wanner (29). To prepare competent cells for transformation, BW25113 (lacIq rrnBT14 {Delta}lacZWJ16 hsdR514 {Delta}araBADAH33 {Delta}rhaBADLD78; 29) containing pKD46 was cultured at 30 °C in SOB broth (28) containing 100 µg of ampicillin/ml. When the A600 reached 0.5, the culture was centrifuged at 4000 rpm for 5 min, and the cells were washed 3 times with cold water before being resuspended in 1% of the original culture volume of water. PCR methods were used to amplify the kanamycin resistance gene (km) from pKD4 by using the primers iclRFw (5'-GTTCAACATTAACTCATCGGATCAGTTCAGTAACTATTGCATTAGCTAACGTGTAGGCTGGAGCTGCTTC-3') and iclRrv (5'-GCGATTAACAGACACCCTTATTCTATTGCCACTCAGGTATGATCATATGAATATCCTCCTTA-3'). The PCR products were purified with a Qiagen kit, treated with DpnI, and repurified by electrophoresis. The km gene was transformed into BW25113-competent cells by electroporation (Gene Pulser; pulse controller at 200 ohms, capacitance at 250 microfarads, and voltage at 25 kV). After electroporation, the cells were grown with shaking in 1 ml of SOC (28) medium at 37 °C for 1 h, and the cultures were plated onto Luria-Bertani (LB) agar containing 25 µg of kanamycin/ml. The Kmr transformants were purified on new kanamycin-LB plates. The mutants in which the target genes were replaced by the km gene were verified by PCR using the primers iclRcfw (5'-CTCATCGGATCAGTTCAGTAAC-3') and iclRcrv (5'-CTTATTCTATTGCCACTCAGG-3'). To delete the km gene from the chromosome, pKD46 was removed from the cells by growing the bacteria at 37 °C, and then pCP20, expressing the Flip recombinase, was introduced by transformation. The transformants containing pCP20 were grown overnight with shaking at 42 °C, and the cultures were plated on LB agar without antibiotics. Colonies were tested for sensitivity to kanamycin and ampicillin. crp mutants were constructed by P1 transduction (30) using the strain JW3320 as donor (National BioResource Project, National Institute of Genetics, Mishima, Japan (31)).

Real-time Quantitative PCR Studies—Bacterial cells were cultured in MOPS minimal media (32) with 10 mM pyruvate as carbon source, and when required, 10 mM glyoxylate was added. When the A600 reached 0.5, cells were collected by centrifugation at 4 °C. Total RNA was subsequently isolated with RiboPureTM-Bacteria (Ambion) in accordance with the manufacturer's protocol. cDNAs were synthesized with the superscript first-strand synthesis kit (Invitrogen) in accordance with the manufacturer's instructions and stored at -80 °C before use. Real-time quantitative PCR was carried out on the Applied Biosystems 7300 apparatus (Applied Biosystems) using Platinum® SYBR® Green qPCR SuperMix UDG (Invitrogen) in accordance with the manufacturer's recommended protocol. Primers used for the real-time PCR were as follows: for rrsC, 5'-CAGCCACACTGGAACTGAGA-3' and 5'-GTTAGCCGGTGCTTCTTCTG-3'; for aceB, 5'-CTGCGTGACCATATTGTTGG-3' and 5'-CAGGCGTGAGTAAGCATTCA-3'; for aceK, 5'-TTCGTGCCTGCTATCAACTG-3' and 5'-AGATATTGAGCGGCACCATC-3'.

Protein Purification—The His-tagged fusion proteins were overexpressed in E. coli BL21-Star(DE3) cells (Stratagene) harboring an extra plasmid encoding three rare tRNAs (AGG and AGA for Arg, ATA for Ile). The cells were grown in LB at 37 °C to an A600 ~ 0.6, and expression was induced with 0.4 mM isopropyl 1-thio-beta-D-galactopyranoside. After the addition of isopropyl 1-thio-beta-D-galactopyranoside, the cells were incubated with shaking at 15 °C overnight. The cells were harvested, resuspended in binding buffer (500 mM NaCl, 5% glycerol, 50 mM HEPES, pH 7.5, 5 mM imidazole), flash-frozen in liquid N2, and stored at -70 °C. The thawed cells were lysed by sonication after the addition of 0.5% Nonidet P-40 and 1 mM each of PMSF and benzamidine. The lysate was clarified by centrifugation (30 min at 17,000 rpm) and passed through a DE52 column preequilibrated in binding buffer. The flow-through fraction was then applied to a metal chelate affinity column charged with Ni2+. After the column was washed, the protein was eluted from the column in elution buffer (binding buffer with 500 mM imidazole). The hexahistidine tag was then cleaved from the protein by treatment with recombinant His-tagged tobacco etch virus protease. The cleaved protein was then resolved from the cleaved His tag and the His-tagged protease by flowing the mixture through a second Ni2+ column.

The purified proteins were dialyzed against 10 mM HEPES, pH 7.5, 500 mM NaCl and concentrated using a BioMax concentrator (Millipore). Before crystallization, any particulate matter was removed from the sample by passage through a 0.2-µm Ultrafree-MC centrifugal filtration device (Millipore). Selenomethionine-labeled proteins were expressed using the same vector and host strain but in supplemented M9 media (33). The sample was prepared under the same conditions as the native protein except for the addition of 5 mM beta-mercapto-ethanol to the purification buffers.

Size-exclusion Chromatography—100 µl of protein sample contained 10 mM HEPES, pH 7.5, 500 mM NaCl, 25 µM IclR wild type or its mutant derivatives, and when indicated, 1 mM glyoxylate or 10 mM pyruvate. After 20 min of incubation on ice, samples were injected onto a Superose 12 10/300 GL gel filtration column (Amersham Biosciences) installed on anÁkta system (Amersham Biosciences) equilibrated with 10 mM HEPES, pH 7.5, 500 mM NaCl. Filtration was performed at 4 °C at a flow rate of 0.5 ml/min, and the protein concentration was monitored by measuring the absorbance at 280 nm. Blue dextran 2000 was used to determine the void volume. A mixture of protein molecular mass standards containingbeta-amylase (200 kDa), bovine serum albumin (66 kDa), carbonic anhydrase (29 kDa), and cytochrome c (12.4 kDa) was applied to the column under similar conditions. The elution volumes and molecular masses of the protein standards were used to generate a standard curve from which the apparent molecular mass was determined.

Crystallization and Data Collection—The primary crystallization condition was determined by a sparse crystallization matrix (Hampton research kits Crystal Screen 1 and PEG/Ion Screen) at room temperature using the sitting drop vapor diffusion technique in 96-well plates. This condition was optimized by varying the pH and the concentration of the solutes. For apoIclR the best crystals were obtained using hanging drops in 0.2 M potassium acetate, pH 7.0, 20% polyethylene glycol 3350 for 2–5 days at room temperature. For diffraction studies, the crystals were flash-frozen with the crystallization buffer plus 30% ethylene glycol as cryoprotectant. For the complex structures the crystals obtained in the crystallization condition described above were soaked in 10 mM solution of either glyoxylate or pyruvate for 1 h. For diffraction studies, the crystals were flash-frozen with the crystallization buffer plus 25% ethylene glycol as the cryoprotectant.

Diffraction data of crystals of C-IclR were collected at 100 K at the 19ID beamline of the Structural Biology Center at the Advanced Photon Source, Argonne National Laboratory. The absorption peak and rising inflection point were determined by calculating and plotting f' and f'' values against energy (34) from the fluorescence spectrum. The three-wavelength anomalous diffraction data were collected from a Se-Met-substituted protein crystal using an inverse beam strategy. All crystallographic data were measured with the custom-built 3 x 3-tiled CCD (charge-coupled device) detector (35). The experiment, data collection, visualization, and processing were controlled with the HKL2000 program package (36). Statistics of data collection and processing are provided in Table 3. Diffraction data for glyoxylate and pyruvate complexes were collected at home source (Rigaku FR-E generator with Raxis4++ detector), and the data were integrated and scaled with *TREK (37).


View this table:
[in this window]
[in a new window]

 
TABLE 3
Data collection and refinement statistics

 
Structure Solution and Refinement—Data from each of the three wavelengths was processed with the HKL2000 (36) and then input into the program SOLVE (38). SOLVE located 24 of a possible 28 selenium sites and gave a mean figure of merit of 0.66 (0.51 in the highest resolution bin). Density modification and automatic model building using RESOLVE (39, 40) resulted in excellent maps and partial models for each of the four subunits. The models were completed with the aid of 4-fold noncrystallographic symmetry and the graphics program O (41). Refinement of atomic positions and individual b-factors using the program CNS-1.1 resulted in an Rfactor of 23.0% (Rfree 29.9%) for data from 19.91 to 2.30 Å. The final model comprises 5499 protein atoms (4 C-IclR domains) and 335 water molecules. The model has excellent stereochemistry as judged by PROCHECK (42), with no Ramachandran violations. All residues of the ligand binding domain except for three C-terminal amino acids (two in the case of subunit D) were located in the experiment. Both glyoxylate and pyruvate complexes structures were solved by molecular replacement using apo-structure as starting model. Structures were refined using REFMAC5 (43, 44).

Electrophoretic Mobility Shift Assays (EMSA)—Electrophoretic mobility shift assays for IclR and mutant proteins were performed using proteins purified and concentrated according the procedures described above. A fragment of the aceBAK promoter region, from position -117 to +17 (containing IclR box II), was generated by PCR using biotin-prelabeled (5'-end) primers and purified using QIAquick spin columns (Qiagen). Incubation mixtures for EMSA (20 µl) contained 2.5 nM 5'-labeled DNA fragment, 50 mM Tris-HCl, pH 7.5, 150 mM KCl, 10 mM MgCl2, 0.01% Triton X-100, 50 ng/µl poly(dI-dC)-nonspecific competitor DNA, purified IclR protein (0–100 nM), and ligand (0–1 mM) where indicated.

After incubation for 20 min at 37 °C, samples were separated on 5% acrylamide-bisacrylamide nondenaturing gels in 0.5x Tris borate-EDTA buffer, pH 8.3 (TBE). Electrophoresis was performed at 100 V using ice-cold 0.5x TBE as a running buffer, and DNA was transferred from the polyacrylamide gel to the Biodyne B Positive Nylon Membrane (Pierce) by electroblotting at 380 mA for 40 min in 0.5x TBE. Transferred DNA was cross-linked for 15 min using UV cross-linker equipped with 312-nm bulbs. Biotin-labeled DNA was detected by horseradish peroxidase/Super Signal detection system (Pierce). Membranes were exposed to Eastman Kodak Co. x-ray film. For quantitative EMSA, the distribution of bound and unbound DNA was quantified using ImageJ (National Institutes of Health).

In Vitro Transcription Runoff Assays—For in vitro transcription runoff the complete aceB promoter was amplified from position -297 to +434 and cloned in a derivative of the pGF-Puv vector (Clontech) in which the gfp gene was replaced with the aceB promotor in the SapI-EcoRI position. The expected specific transcript was 434 bases long. The recombinant plasmids were amplified in E. coli DH5{alpha}, cut with EcoRI, and gel-purified. The reaction mixtures (20 µl) contained template DNA (5 nM), purified IclR or its mutant forms (100 nM), and glyoxylate (1–1000 µM) in a buffer consisting of 40 mM Tris-HCl, pH 8.0, 50 mM KCl, 10 mM MgCl2, 0.01% Triton X-100, and 1 unit of RNase inhibitor. The reaction mix was preincubated for 10 min at 37 °C, and 20 nM RNA polymerase (Epicenter) was added. The reaction mix was further incubated for 5 min at 37 °C, and transcription was initiated by adding 2 µl of a mixture containing 2 mM each ATP, GTP, and CTP, 1.5 mM UTP, and 0.5 mM biotin-11-UTP (PerkinElmer Life Sciences). After 20 min at 37 °C, the reactions were terminated with 10 mM EDTA. The ethanol-precipitated transcripts were analyzed on 6% acrylamide, 7 M urea gels. Electrophoresis was performed at 100 V using ice-cold 0.5x TBE as a running buffer, and RNA was transferred from the polyacrylamide gel to the Biodyne B positive nylon membrane (Pierce) by electroblotting at 380 mA for 40 min in 0.5x TBE. Transferred RNA was cross-linked for 15 min using UV cross-linker equipped with 312-nm bulbs. Detection was performed as described above. Preheated biotinylated 2-log DNA ladder (BioLabs) was used as molecular weight marker.

Screening Using Static Light Scattering—Purified C-IclR protein was screened against a library of 160 compounds using static light scattering (StarGazer) (45). Protein samples were diluted to a final concentration of 0.4 mg/ml in 100 mM Hepes, pH 7, 150 mM NaCl (26). 50-µl aliquots of protein solution containing the chemical compounds at 1 mM were placed in duplicate into clear-bottom 384-well plates (NuncTM, Nalgene Nunc International) and heated from 27 to 80 °C at the rate of 1 °C/min. The formation of protein aggregates was monitored by static light scattering. Images of scattered light were captured every 30 s, and the light intensities were translated to arbitrary numbers using StarGazerTM proprietary software. Intensities were plotted against temperature for each sample well, and transition curves were fitted using the Boltzmann equation. The midpoint of each transition was calculated and compared with the one calculated for the reference sample. If the difference between them was greater than 1.5 °C, the corresponding compound was considered to be a "hit," and the experiment was repeated to confirm the effect.

Circular Dichroism (CD) Spectroscopy—CD spectra were collected at 25 °C on a AVIV 62D CD spectrophotometer from 200 to 260 nm using a 0.1-mm path length cell with a scan rate of 100 nm/min, a time constant of 1.0 s, a bandwidth of 1 nm, and a sensitivity of 100 millidegrees. Each spectrum was the average of 10 scans. After background subtraction, the spectrum was expressed in molar ellipticity. The proteins were analyzed at a final concentration of 1 mg/ml, and the assays were performed in 10 mM HEPES, pH 7.5, 0.5 M NaCl. When necessary, glyoxylate or pyruvate was added at 1 mM.

Isothermal Titration Calorimetry—Measurements were performed on a VP-Microcalorimeter (MicroCal, Northampton, MA) at 25 °C. The protein was thoroughly dialyzed against 10 mM HEPES, pH 7.5, and 500 mM NaCl. Solutions of glyoxylate (0.5 mM) and pyruvate (1.5 mM) were directly prepared in dialysis buffer. Each titration involved a series of 4-µl injections of effector molecules into the protein solution. The mean enthalpies measured from injection of the ligand in the buffer were subtracted from raw titration data before data analysis with ORIGIN software (MicroCal). Titration curves were fitted by a nonlinear least squares method to a function for the binding of a ligand to a macromolecule (46). From the curve thus fitted, the parameters {Delta}H (reaction enthalpy), KA (binding constant, KA = 1/KD), and n (reaction stoichiometry), were determined. From the values of KA and {Delta}H, the changes in free energy ({Delta}G) and in entropy ({Delta}S) were calculated with the equation {Delta}G = -RT ln KA ={Delta}H - T{Delta}S, where R is the universal molar gas constant, and T is the absolute temperature.


    RESULTS
 TOP
 ABSTRACT
 INTRODUCTION
 EXPERIMENTAL PROCEDURES
 RESULTS
 DISCUSSION
 REFERENCES
 
Screening of a Small Molecule Library for Binding the IclR Significant sequence similarity between IclR and previously characterized AllR regulators, particularly among amino acids involved in glyoxylate binding (20), indicated that the IclR effector may have chemical features similar to glyoxylate rather than phosphoenolpyruvate, which was proposed as the IclR effector by Cortay et al. (8). To clarify the nature of the IclR effector molecule, we decided to perform an unbiased screen of IclR against a variety of glyoxylate-associated compounds along with phosphoenolpyruvate as part of large set of metabolic scaffolds.


Figure 1
View larger version (10K):
[in this window]
[in a new window]

 
FIGURE 1.
Circular dichroism of C-IclR with (squares) and without (circles) 1 mM glyoxylate. The proteins were used at a final concentration 1 mg/ml, and the assays were performed in 10 mM HEPES, pH 7.5, 0.5 M NaCl.

 
The expression of the IclR Ser-97 to Arg-287 fragment produced a soluble polypeptide (see "Experimental Procedures" for details) corresponding to the effector binding domain. This protein, termed C-IclR, was tested for binding against 158 metabolic scaffolds (47) by high-throughput protein stability assay using static light scattering technology (25, 26). The library contained intermediates of the glyoxylate bypass and the three carboxylic acids cycle, including glyoxylate, malate, oxaloacetate, citrate, isocitrate, {alpha}-ketoglutarate, succinate, and fumarate as well as phosphoenolpyruvate. The glyoxylate binding domain of AllR regulator, C-AllR (20), was also tested against the same library as a control. The thermal melt conditions for C-IclR and C-AllR were established to generate interpretable unfolding data (described under "Experimental Procedures"). The compounds that induced a shift in the midpoint transition temperature ({Delta}Tm) of proteins by more than 1.5 °C are listed in Table 1.


View this table:
[in this window]
[in a new window]

 
TABLE 1
Stabilization effect of ligand binding The thermal stabilization of each protein by 1 mM ligand concentration was evaluated using differential static light scattering.

 
Glyoxylate was identified as the strongest thermostabilizing compound for C-IclR. The shift in C-IclR {Delta}Tm in the presence of glyoxylate (11.1 °C) was close to the effect of this compound (9.9 °C) on C-AllR. Additional compounds inducing significant stabilization of C-IclR included allantoate, pyruvate, and phosphoenolpyruvate. Apart from glyoxylate, allantoate, DL-glyceraldehyde 3-phosphate, and phosphoenolpyruvate, C-AllR was stabilized by allantoin but not by pyruvate (Table 1). The thermal stabilizing effect of glyoxylate on C-IclR was confirmed by circular dichroism spectroscopy (Fig. 1). In agreement with the differential light-scattering results, the circular dichroism analysis showed an increase in 12 °C in the melting temperature of C-IclR in the presence of glyoxylate.

All five compounds identified for C-IclR contained a common chemical motif in the form of adjacent carbonyl or hydroxyl groups (Table 1) with the strongest thermal stabilizing agent, glyoxylate, being the simplest representation of this chemical scaffold. The intermediate role played by glyoxylate in the pathway regulated by IclR also strengthens its case as a candidate for the native IclR effector molecule. At the same time identification of several metabolites causing IclR thermal stabilization and having a common chemical signature indicates that this sensor protein might be reactive to more than one effector molecule.

Effect of Small Molecule Modulation of IclR Binding in Vitro The effect of glyoxylate, allantoate, DL-glyceraldehyde-3-phosphate, pyruvate, and phosphoenolpyruvate on IclR activity was tested in vitro by EMSA. The binding of IclR to the DNA fragment corresponding to the promoter region of the aceBAK operon and containing IclR box II (5) was analyzed in the presence and the absence of 100 µM or 1 mM concentrations of each compound. According to EMSA results, glyoxylate caused dissociation of the IclR regulator from DNA, thus increasing the fraction of unbound DNA (Fig. 2). On the contrary, pyruvate and, to a lesser extent, DL-glyceraldehyde 3-phosphate and phosphoenolpyruvate stabilized the IclR-DNA complex and resulted in a decrease of unbound DNA (Fig. 2). Allantoate did not affect the IclR/aceB complex (data not shown).


Figure 2
View larger version (36K):
[in this window]
[in a new window]

 
FIGURE 2.
Effect of different ligands on the binding of IclR to the aceB promoter. EMSA results using 2.5 nM biotin-labeled aceB promoter and 25 nM (A) or 50 nM (B) IclR with different inducer molecules. The full binding conditions are described under "Experimental procedures." Lanes 1 and 11 show the migration of the target DNA fragment; lanes 2–10, IclR 25 nM; lanes 12–14, IclR 50 nM. Glyoxylate (lanes 3, 4, and 13), pyruvate (lanes 5, 6, and 14), phosphoenolpyruvate (lanes 7 and 8), and glyceraldehyde 3-phosphate (lanes 9 and 10) were added at a concentration of 0.1 mM (3, 5, 7, 9, and 13) or 1 mM (4, 6, 8, 10, and 14).

 


Figure 3
View larger version (52K):
[in this window]
[in a new window]

 
FIGURE 3.
In vitro transcription runoff experiments using the aceB promoter, IclR, and different concentrations of glyoxylate. A, in vitro transcription assays were performed as described under "Experimental Procedures." The assays were performed at 37 °C for 30 min in the absence (lane 1) or in the presence of 100 nM IclR added before 20 nM RNA polymerase (lanes 2–6). Glyoxylate was added at 1 µM (lane 3), 5 µM (lane 4), 10 µM (lane 5), and 50 µM (lane 6). The arrows denote the aceB (434 nucleotides) and the bla (860 nucleotides) transcripts. B, quantitative analysis of the in vitro transcription runoff experiment presented on A. For each lane the amount of aceB transcript was normalized (average of three experiments) to the expression level of the bla transcript. Expression of aceB in the absence of IclR represented 10% of the expression of bla transcript (internal control). Each bar represents the percentage of the aceB transcript relative to the transcription level in the absence of IclR protein (lane 1).

 


Figure 4
View larger version (43K):
[in this window]
[in a new window]

 
FIGURE 4.
Effect of glyoxylate and pyruvate on the binding of IclR to its operator site. A, EMSA results using 2.5 nM biotin-label aceB promoter, 50 nM of IclR, and increasing concentrations of glyoxylate (0–500 µM) as shown on the top of the picture. Pyruvate 500 µM was added to a second set of samples that also contained increasing concentrations of glyoxylate as shown on the top of the panel. No protein was added to the first lane. B, EMSA results using 25 nM IclR and increasing concentrations of pyruvate (0–500 µM) as shown on the top of the picture. Glyoxylate 500 µM was added to a second set of samples that also contained increasing concentrations of pyruvate as shown on the top of the panel. C and D, EMSA results from A and B, respectively, quantified by image analysis. Each graph represents the percentage (average of three experiments) of the bound DNA in the presence of the ligand relative to the amount of the DNA bound to the same protein in the absence of the ligand.

 
In vitro transcription run-off (Fig. 3) and EMSA (Fig. 4) experiments were performed at different concentrations of glyoxylate. In the in vitro transcription run-off experiment the level of the aceB transcript was normalized to the concentration of the bla transcript used as internal control. The presence of IclR reduced the level of the aceB transcript to 30% of the level measured in the absence of this protein (Fig. 3, lanes 1 and 2). The addition of glyoxylate resulted in gradual increase of the aceB transcription, reaching up to 74% of the original value at 50 µM glyoxylate (Fig. 3, lane 6). In EMSA experiments more than half of the bound operator DNA was released by IclR in the presence of 5 µM glyoxylate (Fig. 4A), whereas 50 µM pyruvate was required to achieve 100% DNA shift (Fig. 4B). Although a substantial difference in the concentrations of glyoxylate needed to decrease IclR binding to DNA was observed between the two techniques tested, they both confirmed the dissociation of the IclR-operator complex in a concentration-dependent manner. This apparent discrepancy could be partially explained by the IclR interactions with the C-terminal domain of the RNA polymerase {alpha} subunit, as described earlier (5).

Phosphoenolpyruvate and DL-glyceraldehyde 3-phosphate were able to stabilize the IclR-DNA complex (Fig. 2); however, the concentrations required to obtain this effect (over 10 mM, data not shown) were out of the physiological range for these compounds (4750). Thus, phosphoenolpyruvate and DL-glyceraldehyde 3-phosphate were not further considered for an IclR effector role.

Because glyoxylate and pyruvate had distinctive and antagonist effects on the formation of the IclR-DNA complex, in the next series of EMSA experiments we analyzed the combined effect of glyoxylate and pyruvate on the stability of the IclR-DNA complex. In this setting one of compounds was added in increasing amounts, whereas the concentration of the other was kept constant. The results obtained (Fig. 4A) demonstrated that the glyoxylate-induced dissociation of the IclR-operator complex was overturned by a 50-fold excess of pyruvate (500 µM pyruvate with 10 µM glyoxylate). These results were corroborated by doing the reciprocal experiment (Fig. 4B). The common chemical features of these compounds suggested that this effect might be due to the competition for the same binding site in IclR, with glyoxylate having a significantly higher binding affinity than pyruvate.

To further characterize IclR interactions with potential effector molecules, the thermodynamic properties of IclR interactions with glyoxylate and pyruvate were determined using isothermal titration calorimetry. The titration of C-IclR with each compound followed an exothermal heat change profile, giving rise to a sigmoidal binding curve with glyoxylate or hyperbolic with pyruvate (Fig. 5). The calculated thermodynamic parameters are summarized in Table 2. The stoichiometry of the reaction, 0.5, is consistent with the binding of one ligand molecule per IclR dimer. In accordance with EMSA results, the C-IclR dissociation constant (KD) for glyoxylate (0.9 ± 0.2 µM) was significantly lower (~150 times) than that for pyruvate (156.2 ± 7.9 µM). Although the KD values obtained for binding of glyoxylate and pyruvate to the full-length IclR were higher (14.2 ± 1.2 and 321.7 ± 20 µM, respectivelytively), their relative ratio was within the same range as obtained for C-IclR.


View this table:
[in this window]
[in a new window]

 
TABLE 2
Thermodynamic parameters derived from the calorimetric titration of C-Iclr and C-AllR with effector molecules Protein solutions at 50 µM in 10 mM KH2PO4/K2HPO4, pH 7.4, and 500 mM NaCl were titrated with 0.5-1.5 mM solutions of effectors. Experiments were carried out at 25 °C. Further assay conditions are given under "Experimental Procedures."

 
In conclusion, our EMSA results show glyoxylate as the strongest candidate for the IclR effector molecule, which destabilized IclR interactions with DNA in vitro. By contrast, pyruvate affected IclR interactions by stabilizing the IclR-DNA complex.

Modulation of IclR Binding by Glyoxylate in Vivo—To assess the effect of glyoxylate on IclR activity in vivo, the expression level of aceBAK operon was measured during growth on minimal media in the presence and absence of this compound. Although early genetic studies identified IclR as the main regulator of aceBAK operon (51), recent E. coli transcriptome analysis indicated that this operon is also under control of the global catabolite regulation protein, Crp. The Crp potential binding sequence was found overlapping the IclR box II (52). In light of this new data, the IclR activity was monitored in the presence and the absence of a functional Crp regulator.

Single {Delta}iclR and {Delta}crp as well as double ({Delta}iclR{Delta}crp) deletion strains were prepared in the E. coli BW25113 strain background. Deletion strains were grown to exponential phase on MOPS minimal medium with pyruvate (10 mM) as the carbon source in the presence and in the absence of 5 mM glyoxylate. The expression levels of aceB and aceK genes were measured using real-time quantitative PCR. The expression profiles of both genes followed the same trend; thus, only the results for aceB gene will be discussed further on.


Figure 5
View larger version (32K):
[in this window]
[in a new window]

 
FIGURE 5.
Isothermal titration calorimetry data for the binding of glyoxylate (A) or pyruvate (B) to C-IclR. Heat changes (upper panel) and integrated peak areas (lower panel) for the injection of a series of 4-µl aliquots of 500 µM glyoxylate or 1.5 mM pyruvate in a solution of 50 µM protein. Data were fitted with ORIGIN.

 
The expression of aceB in the {Delta}iclR and {Delta}crp strains was increased by 3.3- and 24-fold, respectively, compared with the wild type strain. The double {Delta}iclR{Delta}crp deletion caused a 31-fold increase in aceB expression corresponding to the combined effect of each gene deletion. The presence of glyoxylate in the media had no significant effect on the level of aceB expression in {Delta}iclR and {Delta}iclR{Delta}crp strains lacking the functional IclR regulator, whereas it triggered a 1.4-fold increase when the IclR regulator was present in the cells ({Delta}crp strain) (Fig. 6). The glyoxylate-mediated induction was dose-dependent since in the presence of 10 mM glyoxylate a 2.7-fold induction in aceB expression was observed (data not shown). In the case of the wild type strain containing both IclR and Crp regulators, the presence of glyoxylate in the media did not result in a detectable increase of aceBAK, expression probably due to effector-independent tight repression by the functional Crp regulator.

The substantial derepression effect of the Crp deletion confirmed the major role played by this global regulator on the expression of the aceBAK operon. The effect of the IclR deletion was significantly smaller yet consistent with its role as a specific regulator of this operon. The derepression effect in the presence of glyoxylate observed in the {Delta}crp strain appears to be specific to the presence of functional IclR, indicating that glyoxylate can affect the activity of the IclR regulator in a manner expected from the native effector molecule. Unfortunately we were unable to study the corepressor effect of pyruvate in vivo since the crp deletion strain is unable to grow on the variety of carbon sources (52) required for these kinds of studies.


Figure 6
View larger version (13K):
[in this window]
[in a new window]

 
FIGURE 6.
Effect of glyoxylate on the expression of aceB on different strains measured by quantitative real-time PCR. The strains BW25113 (wild type), {Delta}iclR, {Delta}crp, or {Delta}iclR{Delta}crp were grown under aerobic conditions in MOPS with 10 mM pyruvate as carbon source with (light gray bars) or without 5 mM glyoxylate (dark gray bars). RNA extractions and real-time quantitative PCR was performed as described under "Experimental Procedures." The amplification values obtained were corrected for those obtained using rrsC as internal control. The values shown are relative to those observed for the same gene in the wild type strain grown with pyruvate.

 
Structure of the IclR Effector Binding Domain—To provide a structural framework for IclR interactions with the newly identified small molecules, the crystal structures of the IclR effector binding domain alone (apo structure) and in complex with glyoxylate and pyruvate were determined at 2.3, 1.8, and 2.3 Å, respectively. Final R/Rfree for the apo structure was 0.23/0.299 (PDB code 1TD5). Final R/Rfree for complex structures was 0.215/0.262 (PDB code 2O99) for the glyoxylate and 0.177/0.233 (PDB code 2O9A) for the pyruvate complexes. All three models possessed excellent stereochemistry as judged by PROCHECK (42), with no Ramachandran violations. Statistics for all three structures are shown in Table 3. To ease the description and comparison to other structures, the assigned numbers for beta sheets and {alpha}-helices are given in relation to the full-length homolog TM0065 (16).

The C-IclR domain displays a similar architecture to the two other IclR family members whose structures have been determined, TM0065 (16) and AllR (20). The C-IclR structure comprises a centrally located six-stranded anti-parallel beta-sheet surrounded on one side by two long {alpha}-helices ({alpha}5 and {alpha}9) and on the other by three shorter {alpha}-helices ({alpha}6–{alpha}8). The beta-sheet is strongly curved, with the shortest helix ({alpha}6) fitting on the inside of the half-barrel (Fig. 7).

In complex structures both the glyoxylate and pyruvate molecules occupied the cavity corresponding to the effector binding pocket in the AllR structure (20). The map quality was assessed by constructing an omit Fo - Fc electron density map. The map was calculated after 10 cycles of refinement without the ligand in the model (Fig. 8). Similar to the C-AllR-glyoxylate complex structure, no significant conformational changes were observed in C-IclR complexes with pyruvate or glyoxylate when compared with C-IclR apo structure. The ligands were bound in the same orientation with almost identical positions for atoms C1, C2, O1, O2, and O3 (Fig. 8). Oxygen atoms of glyoxylate and pyruvate were coordinated by an extensive network of hydrogen bonds with the side chains and backbone of the protein. For each ligand the O1 atom is hydrogen-bound to the Asp-212O{delta}2 and Gly-160N, the O2 atom to the Ser-239O{gamma} and Ser-241O{gamma}, and the O3 to the Ala-161N. The C1-Sg distance varies from 1.86 to 2.09 Å for glyoxylate complex and from 1.90 to 1.97 Å for the pyruvate complex. Both complex models were refined with the same restraints with no covalent bond between ligand C1 and Cys-222S{gamma} in the dictionary. In agreement with previous sequence analysis, all but two amino acids, Ile-134 and Leu-143, in the effector binding area are conserved between IclR and AllR (20). The corresponding positions in AllR are represented by leucine (Leu-129) and methionine (Met-138), respectively (Fig. 9).


Figure 7
View larger version (34K):
[in this window]
[in a new window]

 
FIGURE 7.
Overall structure of C-IclR. A, schematic diagram of a monomer of the IclR ligand binding domain in complex with glyoxylate (PDB code 2O99). B, C-IclR interdomain interface between two monomers in the (PDB code 2O99) structure. Loops are colored green, {alpha} helices are red, and beta-strands are yellow. Glyoxylate is represented as a stick figure in cyan. The N and C termini are labeled.

 
The presence of an isoleucine residue at position 134 in the IclR structure results in the shortening of the distances between Leu-154 and glyoxylate atoms by2Åinthe IclR-glyoxylate complex structure compared with the same region in C-AllR-glyoxylate structure (contact distance C{delta}1 to C2, 3.68 versus 4.44 Å; C{delta}2 to C1, 4.17 versus 6.16 Å). In addition, the average distances between protein atoms are shorter by an average of 1 Å in C-IclR-glyoxylate versus C-IclR structures (10.8 and 11.6 Å, respectively), whereas there is no such difference between C-AllR and C-AllR-glyoxylate structures (11.96 and 12.1 Å, respectively) for the same parameters. These structural dissimilarities between binding cavities in C-IclR and C-AllR may contribute toward close to a 10-fold difference in KD for glyoxylate as determined by isothermal titration calorimetry (Table 2).

Taking into account the almost identical position of glyoxylate and pyruvate in the C-IclR structures, the interactions involving the methyl group of pyruvate, which is absent in glyoxylate, become critical in understanding the different effects of these compounds on IclR activity. The methyl group is well defined in the pyruvate electron density map (Fig. 8B) and is oriented toward the side chain of the Leu-143 residue belonging to the other molecule of the C-IclR dimer. The distances between C3 of the pyruvate and C{delta}2 atom of Leu-143 are on average 4.24 ± 0.4 Å for the four molecules in the C-IclR asymmetric unit. Leu-143 makes part of the hydrophobic patch formed by side chains of Met-146, Leu-154, His-216, and Leu-220, which are in the range of van de Waals interactions with the pyruvate methyl group. The presence of methionine instead of leucine in the corresponding position in AllR would explain why pyruvate had no stabilization effect (Table 1).

In C-IclR the interdomain interactions between C-IclR monomers are similar to those observed in the C-AllR and TM0065 structures. The C-IclR and C-AllR dimers are superimposable to a root mean square deviation (r.m.s.d.) of 1.75 Å, compared with an r.m.s.d. of 1.01 Å for individual subunits. This common arrangement of the C-terminal domains of three different IclR family members in three different unit cells supports the physiological relevance of these interactions for tetramerization of IclR regulators (Fig. 7B). The additional intermolecular interactions between the pyruvate methyl group and hydrophobic patch observed in the C-IclR-pyruvate complex structure would, thus, stabilize the tetramerization state of IclR protein.

Effect of Glyoxylate and Pyruvate on the Oligomerization of IclR—To establish the effect of glyoxylate or pyruvate on the oligomeric state of IclR protein, this protein was preincubated with 1 mM glyoxylate or 10 mM pyruvate and analyzed by gel filtration. The untreated IclR protein was used as the control.

In the absence of effector molecules IclR was eluted in two peaks corresponding to the dimeric (60 ± 5 kDa) or tetrameric (120 ± 12 kDa) states. The abundance of both fractions (54/46% relative percentage, respectively) indicated that under our experimental conditions the IclR protein is in a dynamic equilibrium between these two oligomeric states (Fig. 10, Table 4). After preincubation with glyoxylate, IclR was eluted predominantly in the form of the lower (60 kDa) molecular mass fraction, indicating a dramatic shift of equilibrium toward the dimeric state. On the other hand, preincubation with 10 mM pyruvate resulted in elution of IclR in the form of the higher (120 kDa) molecular mass fraction corresponding to the tetramer. Similar results were obtained if glyoxylate was included in the elution buffer; however, the tetrameric state of IclR was more favorable if pyruvate was included in the elution buffer. These results can be explained by the difference in affinity between the ligands (Table 2). Overall, size exclusion experiments demonstrated that glyoxylate and pyruvate have a dramatic effect on the oligomeric state of the IclR regulator. This suggests that the glyoxylate effect on IclR activity is due to the destabilization of tetramers, which is required for efficient binding to the operator DNA, whereas the co-repressor effect observed for pyruvate might be due to the stabilization of the IclR tetramer on the DNA.


View this table:
[in this window]
[in a new window]

 
TABLE 4
Summary of the effect of glyoxylate and pyruvate on DNA binding and on the oligomeric state of IcIR or its mutant variants

 


Figure 8
View larger version (55K):
[in this window]
[in a new window]

 
FIGURE 8.
Molecular details of C-IclR interactions with glyoxylate (A and C) or pyruvate (B and D). In A and B the side chains of amino acids within a distance of 8 Å are shown and labeled in the figure. Protein residues (yellow), glyoxylate (pink), and pyruvate (green) are represented as stick figures. The backbone of the protein is shown as a thin yellow line. Leu-143 from another monomer is represented in cyan.In C and D a3{sigma} Fo - Fc omit map was built using the REFMAC5 program (43, 44) for the glyoxylate and pyruvate complex structures, respectively. The map was calculated after 10 cycles of refinement without the ligand in the model. The omit Fo - Fc electron density map covering the ligands is shown in stick representation (blue for C or cyan for D). Neighboring residues are also shown in stick representation (cyan in C or yellow in D).

 


Figure 9
View larger version (45K):
[in this window]
[in a new window]

 
FIGURE 9.
Ligand binding pocket superimposition of AllR and IclR. Superimposition of the active site of C-IclR with that of C-AllR. Protein residues and glyoxylate are represented as stick figures. Components of C-AllR are represented in green and light pink, whereas components of IclR are colored in yellow and cyan.

 
Functional Analysis of Key Amino Acids Involved in Interdomain Interactions and Effector Binding—In light of the important effect of glyoxylate and pyruvate on the IclR oligomeric state, we investigated the specific role of amino acids Leu-143, Met-146, Leu-154, and Leu-220 in the interdomain interactions and effector binding. In the course of this analysis each of these amino acids was replaced by alanine or aspartic acid. The mutant proteins were purified and submitted to gel filtration analysis and EMSA with the IclR box II fragment in the presence and absence of glyoxylate or pyruvate (Table 4, Figs. 11 and 12).

The alanine substitution of Leu-143 did not significantly affect the gel filtration profile of IclR (Fig. 11). However, according to the EMSA experiments a higher concentration (100 nM) of the L143A mutant protein was required to achieve the DNA binding comparable with the wild type IclR (Table 4). The addition of pyruvate stabilized complex of this mutant with operator DNA similarly to the wild type IclR, whereas the glyoxylate had a significantly lower effect on the L143A mutant (Fig. 12). In the presence of the latter ligand an additional band was observed on the EMSA gels. This additional band migrated faster than the one corresponding to the wild type IclR-DNA complex but slower that the free DNA (Fig. 12). This band would correspond to an IclR-DNA complex with a lower molecular weight than in case of the wild type IclR. Taking into account that IclR might bind DNA as a tetramer (5), this lower molecular weight complex would correspond to the association of operator DNA with the dimers of L143A as an intermediate in the releasing of IclR from the protein-DNA complex.

The gel filtration profiles of the alanine substitutions of Met-146, Leu-154, and Leu-220 showed them predominantly as dimers in solution (70/30 relative percentage, respectively, Table 4); however, they were still able to tetramerize over DNA as evidenced by a high molecular weight migrating complex (Fig. 12). The mutant proteins showed a significant variation of responses to glyoxylate and pyruvate (Table 4). Although the M146A mutant response to these ligands was similar to the wild type IclR protein, the L220A mutant showed no response to glyoxylate but was stabilized by pyruvate. The L154A mutant, to the contrary, responded to glyoxylate but not to pyruvate. These results put in evidence the different involvement of these amino acids in ligand binding.


Figure 10
View larger version (24K):
[in this window]
[in a new window]

 
FIGURE 10.
Gel filtration analysis of IclR oligomerization in presence of glyoxylate or pyruvate. IclR (thick line), IclR preincubated with pyruvate (continuous line), or with glyoxylate (dotted line) was passed through a Superose 12 column as described under "Experimental Procedures." The inset shows the calibration curve. The standards used were beta-amylase (200 kDa), bovine serum albumin (66 kDa), carbonic anhydrase (29 kDa), and cytochrome C (12.4 kDa). mAU, milliabsorbance units.

 
In the next series of mutagenesis the hydrophobic region was altered by single replacement of Leu-143, Met-146, Leu-154, or Leu-220 residues to the charged residue, aspartate. The mutations had a dramatic effect on IclR oligomerization, with all mutant proteins being dimers as judged by gel filtration analysis (Table 4). According to EMSA results, L143D and L220D bound the operator DNA at a higher protein concentration (100–150 nM) if compared with the wild type regulator (50 nM) and formed only low molecular weight complexes with DNA. Both mutants were not able to respond to either glyoxylate or to pyruvate (Fig. 12).

EMSA experiments on M146D showed low and high molecular weight protein-DNA complexes (Fig. 12). In agreement with the mutations to alanine, M146D does not interfere in the binding of glyoxylate or pyruvate. Moreover, glyoxylate disrupted the formation of the high molecular weight protein-DNA complex, increasing the amount of the lower molecular weight complex (Fig. 12). The presence of pyruvate, on the contrary, increases the formation of the high molecular weight complex. These results show that although the Met-146 is important for the stabilization of tetramers in solution, it is not directly involved in the binding of the ligands, although it might contribute to the overall hydrophobicity of the area.

In the case of L154D, the EMSA profile was similar to that of the wild type IclR; however, the response to glyoxylate was weaker (Fig. 12). Pyruvate had no stabilization effect on this mutant protein (Fig. 12), confirming the Leu-154 role in the specific interactions with pyruvate.


    DISCUSSION
 TOP
 ABSTRACT
 INTRODUCTION
 EXPERIMENTAL PROCEDURES
 RESULTS
 DISCUSSION
 REFERENCES
 
The IclR family of transcriptional regulators shares common structural features with their C-terminal domains representing small molecule binding modules adapted to recognize a wide variety of cellular metabolites as effector molecules. The identity of the effector molecules remains unknown for a vast majority of these proteins. This reflects a general situation within the large transcriptional factor families such as LysR (53), LuxR (54), TetR (55), and GntR (56), accounting for thousands of transcriptional regulators that respond to unidentified chemical stimuli which need to be addressed.


Figure 11
View larger version (21K):
[in this window]
[in a new window]

 
FIGURE 11.
Analysis of IclR L143A and M146A oligomerization in presence of glyoxylate and pyruvate. L143A (A) and M146A (B)(thick line) were preincubated with pyruvate (continuous line), or with glyoxylate (dotted line) and passed through a Superose 12 column as described under "Experimental Procedures." mAU, milliabsorbance units.

 
Using the founding member of the IclR family as a test case, we explored the possibility of screening a small molecule library using a high throughput protein stability assay as a general strategy to identify putative ligands for such families' representatives. Although the phenomenon of protein thermostabilization by a bound small molecule was known for decades (5759), only recently it became applicable in screening chemical libraries (25, 26, 60). The breakthrough came from the development of technologies enabling the thermo-denaturation of a large number of protein samples in parallel monitoring by static light scattering or fluorimetry adapted to microplate format (25, 26). This kind of technology enables the high-throughput screening of protein-ligand interactions independent of the knowledge of the protein function. A similar technology named ThermoFluor was previously used to identify and elaborate upon the hypothetical function of a pyridoxal phosphate-dependent aminotransferase enzyme (61). In this study we screened a small molecule library containing representatives of the main chemical scaffolds identified in the E. coli metabolome (47) as well as a large number of chemicals associated with the three carboxylic acids cycle and the glyoxylate bypass. Glyoxylate was identified as the strongest thermo-stabilizing compound for IclR. In EMSA and in vitro transcription runoff experiments the presence of glyoxylate disrupted IclR complex formation with operator DNA, an effect that would be expected from a physiologically relevant effector. Allantoate, DL-glyceraldehyde 3-phosphate, pyruvate, and phosphoenolpyruvate also had a stabilizing effect on IclR. Surprisingly, three of these compounds stabilize the binding of IclR to DNA. All five compounds demonstrated a common chemical moiety in the form of vicinal carboxyl groups. Structural analysis of the C-IclR complexes with glyoxylate and pyruvate clearly showed that this chemical substructure is the key element in interactions with IclR. This observation indicates that the chemical screen-based approach can provide the necessary information for identification of chemical moieties specifically recognized by IclR family members as well as for other effector-binding proteins in general.


Figure 12
View larger version (66K):
[in this window]
[in a new window]

 
FIGURE 12.
Functionalanalysisofkeyaminoacidsinvolvedininterdomaininteractionsandeffectorbinding. The glyoxylate and pyruvate binding on mutant IclR proteins was tested by EMSA. The name of the corresponding mutant protein and its concentration used for binding is indicated above each image. Glyoxylate and pyruvate were tested at 1 mM. The full binding conditions are described under "Experimental Procedures."

 
Only two compounds, glyoxylate and pyruvate, altered the binding of IclR in vitro at a physiologically relevant concentration. The role of glyoxylate as IclR effector is supported by its role as an intermediate in the glyoxylate bypass, which is regulated by IclR. The importance of this pathway for cellular metabolism is reflected in the complex regulation of this bypass at transcriptional and post-transcriptional levels. At the transcriptional level it is tightly regulated by the binding of at least six different transcription factors, and it is extremely dependent on the metabolic state of the cell. The aceBAK operon is under negative regulation by IclR, FadR (indirectly by regulation of the iclR expression (62), Crp (52), and ArcA (63). Additionally, it is under positive regulation by integration host factor (IHF) (64) and FruR (65, 66). Once induced at transcriptional level, the flow through the pathway is regulated in part by the inhibitory phosphorylation of isocitrate dehydrogenase (67, 68). Gluconeogenesis also plays an essential role during growth on acetate by synthesizing glucose from nonhexose precursors. The genes pckA and ppsA belong to two parallel pathways for the gluconeogenic conversions of the three carboxylic acid cycle intermediates to phosphoenolpyruvate with the malic enzymes (sfcA and maeB) and ppsA forming one path for the conversion of malate to phosphoenolpyruvate and pckA forming the other for the conversion of oxaloacetate to pyruvate (6971).

Taking into account the physiological data, we could hypothesize that the cell could maximize the efficiency of expression/repression of the aceBAK operon by using two antagonistic effectors on the same transcription factor. By EMSA and biophysical tests we were able to establish that IclR binds glyoxylate and pyruvate. Pyruvate is the first non-phosphorylated intermediate in glycolysis and a key intermediate in catabolic and biosynthetic reaction pathways. There are several pathways that deliver this compound (48, 72, 73). Siddiquee et al. (48, 49) have shown that the intracellular concentration of pyruvate can be from 0.15 to 0.4 mM depending if samples were taken from batch or continuous cultures. These values are within the affinity constants defined for IclR binding. Glyoxylate, on the other hand, is an important intermediate of the central microbial metabolism in the glyoxylate bypass. Glyoxylate is also generated from glycolate or purine degradation (74, 75). The intracellular concentration for this metabolite in E. coli has not been reported; however, there are studies in the closely related bacterium Salmonella typhimurium (76) which suggest that the concentration of glyoxylate should be very low (~10 µM). Thus, the big difference in affinity between glyoxylate and pyruvate binding to IclR mig