RNA Attend a BioResearch Product Faire
HOME HELP FEEDBACK SUBSCRIPTIONS ARCHIVE SEARCH TABLE OF CONTENTS
 QUICK SEARCH:   [advanced]


     


Published online before print August 30, 2005, 10.1261/rna.2990205
RNA (2005), 11:1579-1587. Published by Cold Spring Harbor Laboratory Press. Copyright © 2005 RNA Society.
This Article
Right arrow Abstract Freely available
Right arrow Full Text (PDF)
Right arrow All Versions of this Article:
rna.2990205v1
11/10/1579    most recent
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 DUBEY, A. K.
Right arrow Articles by BABITZKE, P.
Right arrow Search for Related Content
PubMed
Right arrow PubMed Citation
Right arrow Articles by DUBEY, A. K.
Right arrow Articles by BABITZKE, P.
Social Bookmarking
 Add to CiteULike   Add to Connotea   Add to Del.icio.us   Add to Digg   Add to Reddit   Add to Technorati  
What's this?

RNA sequence and secondary structure participate in high-affinity CsrA–RNA interaction

ASHOK K. DUBEY1,3, CAROL S. BAKER1, TONY ROMEO2 and PAUL BABITZKE1

1 Department of Biochemistry and Molecular Biology, The Pennsylvania State University, University Park, Pennsylvania 16802, USA
2 Department of Microbiology and Immunology, Emory University School of Medicine, Atlanta, Georgia 30322, USA

Reprint requests to: Paul Babitzke, Department of Biochemistry and Molecular Biology, The Pennsylvania State University, University Park, PA 16802, USA; e-mail: pxb28{at}psu.edu; fax: (814) 863-7024.


    ABSTRACT
 TOP
 ABSTRACT
 INTRODUCTION
 RESULTS AND DISCUSSION
 MATERIALS AND METHODS
 REFERENCES
 
The global Csr regulatory system controls bacterial gene expression post-transcriptionally. CsrA of Escherichia coli is an RNA binding protein that plays a central role in repressing several stationary phase processes and activating certain exponential phase functions. CsrA regulates translation initiation of several genes by binding to the mRNA leaders and blocking ribosome binding. CsrB and CsrC are noncoding regulatory RNAs that are capable of sequestering CsrA and antagonizing its activity. Each of the known target transcripts contains multiple CsrA binding sites, although considerable sequence variation exists among these RNA targets, with GGA being the most highly conserved element. High-affinity RNA ligands containing single CsrA binding sites were identified from a combinatorial library using systematic evolution of ligands by exponential enrichment (SELEX). The SELEX-derived consensus was determined as RUACARGGAUGU, with the ACA and GGA motifs being 100% conserved and the GU sequence being present in all but one ligand. The majority (51/55) of the RNAs contained GGA in the loop of a hairpin within the most stable predicted structure, an arrangement similar to several natural CsrA binding sites. Strikingly, the identity of several nucleotides that were predicted to form base pairs in each stem were 100% conserved, suggesting that primary sequence information was embedded within the base-paired region. The affinity of CsrA for several selected ligands was measured using quantitative gel mobility shift assays. A mutational analysis of one selected ligand confirmed that the conserved ACA, GGA, and GU residues were critical for CsrA binding and that RNA secondary structure participates in CsrA–RNA recognition.

Keywords: RNA binding protein; translation control; protein–RNA interaction; SELEX


    INTRODUCTION
 TOP
 ABSTRACT
 INTRODUCTION
 RESULTS AND DISCUSSION
 MATERIALS AND METHODS
 REFERENCES
 
The Csr (carbon storage regulation) and homologous Rsm (repressor of stationary phase metabolites) global regulatory systems of several bacterial species control numerous genes and processes post-transcriptionally. Csr (Rsm) systems are characterized by an RNA binding protein that either activates or represses expression of target mRNAs and noncoding RNAs (ncRNAs) that function as antagonists of the RNA binding protein (Romeo 1998Go). For example, in Pseudomonas fluorescens the Rsm system participates in controlling the production of anti-fungal secondary metabolites and exoenzymes (Blumer et al. 1999Go; Reimmann et al. 2005Go), whereas Rsm of Erwinia species regulates a variety of genes involved in soft-rot disease of higher plants (Cui et al. 1999Go). Csr is also known to regulate epithelial cell invasion by Salmonella enterica (Altier et al. 2000Go) and affect swarming and other activities of Proteus mirabilis (Liaw et al. 2003Go) and Legionella pneumophila (Fettes et al. 2001Go).

The Csr system of Escherichia coli is involved in the repression of several stationary phase processes and in the activation of some exponential phase functions (Romeo 1998Go). Three major components of Csr in this organism include the RNA binding protein CsrA and two noncoding RNA (ncRNA) molecules, CsrB and CsrC. CsrA represses various processes such as gluconeogenesis, glycogen metabolism, and biofilm formation (Romeo et al. 1993Go; Sabnis et al. 1995Go; Yang et al. 1996Go; Jackson et al. 2002Go; Wang et al. 2005Go). CsrA also activates glycolysis, acetate metabolism, and flagellum biosynthesis (Sabnis et al. 1995Go; Wei et al. 2000Go, 2001Go). CsrB and CsrC function as antagonists of CsrA by sequestering this protein and preventing its ability to interact with mRNA targets. The multiple imperfect repeat sequences (18 in CsrB and nine in CsrC) in these regulatory RNAs function as CsrA binding sites (Liu et al. 1997Go; Gudapaty et al. 2001Go; Weilbacher et al. 2003Go).

CsrA negatively regulates expression of glgC, a gene involved in glycogen biosynthesis, by binding to four sites in the untranslated leader of the glgCAP operon transcript, one of which overlaps the glgC Shine-Dalgarno (S-D) sequence (Baker et al. 2002Go; A.K. Dubey, T. Romeo, and P. Babitzke, unpubl.). CsrA binding to the glgCAP leader transcript inhibits GlgC synthesis by blocking ribosome binding. Presumably, CsrA-mediated inhibition of glgC translation is responsible for the accelerated rate of glgCAP mRNA decay (Liu et al. 1995Go). CsrA also represses translation of cstA, a carbon starvation-induced gene thought to be involved in peptide transport (Schultz and Matin 1991Go; Dubey et al. 2003Go), as well as the pgaABCD operon, a cluster of genes that are required for the synthesis of the polysaccharide adhesin poly-ß-1,6-N-acetyl-D-glucosamine (PGA) that participates in biofilm formation (Wang et al. 2005Go). CsrA binds to four sites in the cstA transcript and to at least six sites in the pgaA leader transcript. In each case one of the CsrA binding sites overlaps the cognate S-D sequence. Translational repression of these genes proceeds by a mechanism that is similar to that of glgC (Dubey et al. 2003Go; Wang et al. 2005Go).

Considerable sequence variation exists in the known E. coli CsrA binding sites, with GGA being the most conserved element. The GGA motif is predicted to be present in the loop of short RNA hairpins in several of these binding sites, although the importance of this structural arrangement is not known (Liu et al. 1997Go; Baker et al. 2002Go; Weilbacher et al. 2003Go; Wang et al. 2005Go). Systematic evolution of ligands by exponential enrichment (SELEX) is routinely used to isolate high-affinity ligands from a pool of randomized nucleic acid sequences (Ellington and Szostak 1990Go; Tuerk and Gold 1990Go; Ulrich et al. 2002Go). SELEX has been used to identify RNA ligands that either bind to protein (e.g., Schneider et al. 1992Go; Baumann et al. 1997Go) or small metabolites (e.g., Lozupone et al. 2003Go). In this study, SELEX was used to isolate high-affinity CsrA ligands. The results presented herein establish that both the primary sequence and secondary structure of selected RNA ligands are important for high-affinity CsrA interaction.


    RESULTS AND DISCUSSION
 TOP
 ABSTRACT
 INTRODUCTION
 RESULTS AND DISCUSSION
 MATERIALS AND METHODS
 REFERENCES
 
In vitro selection of high-affinity CsrA ligands
The ncRNAs CsrB and CsrC contain 18 and nine apparent CsrA binding sites, respectively. In addition, a total of 14 binding sites have been identified in the glgC, cstA, and pgaA transcripts. Despite the large number of known CsrA binding targets, identification of a reasonable binding site consensus has remained elusive. Since considerable sequence variation exists among the known CsrA binding sites, and because cooperative binding to natural RNAs is common, SELEX was carried out to identify RNA ligands containing single high-affinity CsrA binding sites. The template-primer system used by Ulrich et al. (2002)Go was adapted to produce RNA molecules having a stretch of 15 randomized nucleotides that were flanked by 5' and 3' constant regions. The use of only 15 randomized nucleotides was chosen to minimize the possibility of selecting RNAs containing more than one CsrA binding site. Because the theoretical complexity of 15 randomized positions contains just over 1 x 109 different sequences, our initial randomized RNA pool contained ~1 x 105 molecules of every possible sequence. Since quantitative gel mobility shift assays using native CsrA or C-terminal His-tagged CsrA (CsrA-H6) did not show any significant difference in binding affinities for target transcripts (data not shown), CsrA-H6 was used in SELEX. To remove nonspecific ligands, the RNA pool was subjected to preselection with Ni-NTA beads in the absence of CsrA-H6 protein. RNA molecules that did not bind to the Ni-NTA beads were used for the selection process (see Materials and Methods).

The progress of selection was monitored by gel mobility shift assays using native CsrA and 5' end-labeled RNA pools from rounds 0, 2, 4, 7, and 9 (Fig. 1Go). The binding constants of CsrA–RNA interaction gradually decreased as the number of rounds of selection increased, indicating that subjecting the RNAs to additional rounds of selection enriched for high-affinity ligands. Initially nine rounds of selection were carried out and RT-PCR fragments from rounds 4, 7, and 9 were cloned and sequenced. Three additional rounds of selection were subsequently carried out; however, we did not identify any sequences from round 12 that differed appreciably from those selected in rounds 7 and 9 (see below).



View larger version (63K):
[in this window]
[in a new window]
 
FIGURE 1. Gel mobility shift analysis of CsrA complexed with RNA pools from SELEX rounds 0, 2, 4, 7, and 9. 5' End-labeled RNA was incubated with CsrA at the concentration indicated at the bottom of each lane. Positions of free (F) and bound (B) RNAs are shown. The apparent equilibrium binding constant (Kd) for each RNA pool is shown.

 
Sequence and structural analysis of selected RNA ligands
A total of 55 selected RNA ligands were cloned and sequenced from rounds 4, 7, 9, and 12. Sequence analysis of these RNA molecules identified two ligand classes. Fifty class I RNA molecules contained a single GGA motif, whereas five class II RNAs contained two GGA sequences (Fig. 2Go and data not shown). The RNA ligands from class I were subdivided into two groups. Class I-ARNAs (38 sequences) contained the GGA motif at the extreme 3' end of the randomized sequence, while class I-B RNAs (12 sequences) contained the GGA somewhere in the middle of the randomized sequence. A sequence alignment of all 55 RNA molecules identified a SELEX-derived consensus of RUACARGGAUGU (Fig. 2Go). Strikingly, the ACA and GGA motifs were 100% conserved, while the GU residues were present in all but one RNA ligand. All of the class I-A CsrA binding sites utilized the AGU sequence (following GGA) from the 3' constant region. In all of the class I-B examples, the 3' end of the CsrA binding site was derived from randomized sequence. The presence of UGU (following GGA) in seven out of 12 class I-B RNAs led us to include UGU rather than AGU in the consensus. It is also interesting that four of the class I-B sequences (R9–8, R9–24, R9–54, and R9–56) and one class II sequence (R9–46) derived their AC residues of the conserved ACA element from the 5' constant region (Fig. 2Go). In each of these examples the first A residue originated from a G-to-A substitution that presumably arose in the 5' constant region during PCR. In the case of the class II RNAs, the sequence surrounding only one of the two GGA motifs closely resembles the consensus. When only this GGA motif is considered, four of the class II RNAs are similar to class I-A ligands, with the other being similar to class I-B.



View larger version (59K):
[in this window]
[in a new window]
 
FIGURE 2. Representative selected RNA ligands. The respective class of each RNA ligand is shown (see text for details). The 15 randomized nucleotides are shown in uppercase type, while the constant regions are shown in lowercase type. Residues corresponding to the CsrA binding site are indicated in bold type. The SELEX-derived CsrA binding site consensus is shown at the bottom (R = A or G).

 
The GGA motif of known CsrA binding sites is often found in the loop of short hairpins. Computer modeling of the 55 selected sequences using MFOLD v.3.1 (Mathews et al. 1999Go; Zuker 2003Go) revealed that 51 of the RNAs contained the GGA motif in the loop of a hairpin within the most stable predicted structure ({Delta}G = –17.6 to –6.3 kcal/mol) (Fig. 3Go and not shown). In addition, R9–54 (class I-B, Fig. 2Go) contained its GGA in the loop of a hairpin within the second most stable structure ({Delta}G = –6.6 kcal/mol vs. {Delta}G = –6.2 kcal/mol, not shown). The predicted stabilities of the GGA motif-containing hairpins ranged from –14.2 kcal/mol to –0.7 kcal/mol. The GGA motifs of the remaining three ligands, R9–31 (class II), R9–42 (class II), and R9–49 (class I-B) were also predicted to be present in the loops of suboptimal structures, although in each case their predicted stabilities were substantially lower than the optimal structures (not shown). It is also interesting to note that only three of the hairpins contained more than four contiguous base pairs below the loop, suggesting that particularly stable hairpins might not be favorable for CsrA interaction. Thus, it is possible that CsrA selected all 55 ligands when the GGA motif was present in the loop of a short hairpin.



View larger version (27K):
[in this window]
[in a new window]
 
FIGURE 3. Predicted secondary structures of representative selected RNA ligands. The respective class of each RNA is shown (see text for details). The identity of the purines corresponding to the Rs in the SELEX-derived CsrA binding site consensus (RUACARGGAUGU) is indicated. For R9–44, one of these residues is a C. The apparent CsrA binding site for each transcript is shown in bold type, while the conserved residues predicted to be involved in base-pair formation are boxed. Arrows for R9–31 show a less stable alternative pairing arrangement in which the GGA motif would be present in the loop of a hairpin. The predicted free energy of each structure is shown, while the predicted free energy of the hairpin containing the CsrA binding site is in parentheses.

 
The loops of the hairpins either consisted of four (eight sequences), six (44 sequences), or eight (three sequences) nucleotides. Strikingly, with one exception, the AC residues in the 100% conserved ACA motif were predicted to form base pairs with the conserved GU residues. Moreover, in the cases in which the binding target ended in UGU, the entire ACA motif was predicted to form base pairs with these residues (Fig. 3Go). Finally, it is interesting to note that one class I-B ligand (R9–44) contains a GNRA (GGGA) tetraloop (Fig. 3Go). In this case, the first G in the tetraloop is derived from the second R in the consensus sequence. The observation that the sequence of the seven remaining four-base loops was AGGA, a sequence not known to form tetraloops, suggested that tetraloops are generally unfavorable for CsrA binding. In summary, it is apparent that both primary sequence and RNA secondary structure participate in the formation of high-affinity CsrA binding sites. More striking was the finding that primary sequence information was embedded within the stem.

A selected ligand (R9–43) competes with CsrA–pgaA mRNA interaction
Previously published results demonstrated that CsrA binds specifically to the pgaA mRNA and that this transcript contains six likely CsrA binding sites (Wang et al. 2005Go). Gel mobility shift assays were carried out to determine whether one of the SELEX-derived ligands (R9–43) was capable of competing for CsrA–pgaA mRNA interaction. CsrA binding to the pgaA transcript (+1 to +260 relative to the start of pgaA transcription) was detected as a distinct band in native gels between 2.5 and 40 nM CsrA and an apparent Kd of 25 nM and a cooperativity coefficient of 1.5 (Fig. 4Go, top panel). As the concentration of CsrA was increased further, a shifted complex of slower gel mobility was observed. Competition experiments were carried out with specific (pgaA, +1 to +260 relative to the start of transcription; R9–43) and nonspecific (Bacillus subtilis trp leader, +1 to +49 relative to the start of transcription) RNA competitors. Since the trp leader transcript contained an RNA hairpin (Sudershana et al. 1999Go; Du et al. 2000Go), use of this competitor would allow us to determine whether CsrA was capable of binding to RNA hairpins in general. Both pgaA and R9–43 RNAs were effective competitors, whereas the B. subtilis trp leader RNA did not compete for CsrA–pgaA RNA interaction (Fig. 4Go, bottom panel). These results indicate that CsrA interaction with the R9–43 aptamer and pgaA RNA occurs by a similar mechanism and that CsrA is not capable of interacting with RNA hairpins in general.



View larger version (81K):
[in this window]
[in a new window]
 
FIGURE 4. Ligand R9–43 competes with CsrA–pgaA mRNA interaction. 5' End-labeled pgaA transcript (0.25 nM) was incubated with CsrA at the concentration of CsrA indicated at the bottom of each lane. Gel shift assays were performed in the absence (top) or presence (bottom) of various competitor RNAs. Concentrations of specific (pgaA and R9–43) and nonspecific (trpL of B. subtilis) competitor RNAs are shown at the bottom of each lane. Positions of free (F) and bound (B) RNAs are shown.

 
Binding affinity of CsrA for several selected RNA ligands
Gel mobility shift assays were performed to measure the affinity of CsrA for representative examples of selected RNA ligands (Table 1Go). To determine whether the identity of the purines designated by R in the RUACARGGAUGU consensus sequence influenced CsrA binding, gel shift experiments were carried out with four RNAs from class I-B that contained all possible combinations of A and G at these two positions (Table 1Go; R7–29, R9–9, R9–15, and R9–16). Structural predictions indicated that the GGA motif for each of these RNAs was present in a six-base loop that differed only in the identity of the second R in the consensus sequence (Fig. 3Go). The binding constants for these four RNAs ranged from 8 to 11 nM, indicating that CsrA does not have a significant preference for either purine at these two positions in the context of a six-nucleotide loop. The binding preference of CsrA for A or U following the GGA was also examined. Recall that the A at this position of the class I-A ligands was derived from the 3' constant region, while seven of 12 class B ligands contain a U at this position. In this case, CsrA exhibited a slight preference for U (Table 1Go; R9–24 and R9–54).


View this table:
[in this window]
[in a new window]
 
TABLE 1. Binding affinity of CsrA for selected RNA ligands
 
The binding affinity of CsrA was also determined for each of the class II RNAs. Despite having two GGA sequences, sequence alignments suggested that only one authentic CsrA binding site was present in each transcript (Fig. 2Go). The binding affinity of CsrA for transcripts corresponding to clones R9–22, R9–28, R9–42, and R9–46 were between 7 and 16 nM, values that were similar to the other RNA ligands that were tested (Table 1Go). The finding that only one shifted complex was observed in gel mobility shift assays is consistent with each of these transcripts containing a single CsrA binding site (data not shown). The affinity of CsrA for the R9–31 transcript was considerably lower than for the other RNA ligands that were examined (Table 1Go). One reasonable explanation for the reduced affinity for this RNA molecule is that the GGA motif is sequestered in the stem of a predicted hairpin (Fig. 3Go). The low affinity of CsrA for R9–31 raises the question as to how this RNA was selected in the first place. Perhaps CsrA-H6 bound to this transcript when it was folded into a suboptimal structure in which the GGA motif is present in the loop of a hairpin ({Delta}G = –10.1 kcal/mol vs. {Delta}G = –8.1 kcal/mol) (Fig. 3Go). A possible explanation for the exceptionally high standard error for the binding curve observed with this transcript (Table 1Go) is that CsrA is only capable of binding to this transcript when the suboptimal structure is formed.

Mutational analysis of a selected RNA ligand
The relative contribution of the primary sequence and RNA secondary structure on CsrA binding was examined by altering conserved residues of a high-affinity RNA ligand. R9–43 was chosen as the progenitor wild-type binding site for this analysis because structural predictions using MFOLD indicated that none of the mutations that we intended to introduce would lead to unintended structural rearrangements (Fig. 5Go). The mutations that were introduced altered the primary sequence and/or the predicted secondary structure of R9–43. Gel mobility shift assays were carried out to investigate the effect of each mutation on CsrA binding (Fig. 5Go). The affinity of CsrA for R9–43(WT) was 7 nM (Figs. 5Go, 6Go). The C22U substitution replaced the C–G base pair with a U–G base pair (Fig. 5Go), resulting in a binding site that was identical to the SELEX-derived consensus. As expected, this mutation did not significantly alter the binding affinity. The A23G mutation, which altered a conserved residue and replaced the A–U base pair for a G–U base pair, led to a modest threefold reduction in binding affinity, whereas the A23U mutation, which replaced a conserved residue and disrupted the A–U base pair, led to a four- to fivefold reduction in affinity. The U32A mutation, which altered a conserved residue and disrupted the A–U base pair, resulted in a 10-fold reduction in binding affinity (Figs. 5Go, 6Go). Interestingly, the compensatory A23U:U32A double mutation did not restore binding; the affinity of CsrA for the compensatory mutant was similar to the U32A single mutant. These results indicated that the primary sequence itself, and not just its ability to form a base pair, is critical for high-affinity interaction. A similar pattern was observed for the C24U, C24G, G31C, and C24G:G31C mutations except that the double compensatory change exhibited an additive binding defect. Note that in this case, disruption of the C24–G31 base pair would likely lead to disruption of the A25–U30 closing base pair as well (Fig. 5Go). The A25U substitution, which replaced a conserved residue and disrupted the closing A–U base pair, led to a 15-fold reduction in binding affinity, whereas the U30A substitution did not affect binding (Figs. 5Go, 6Go). Interestingly, all of the SELEX-derived sequences in which the GGA motif was located at the extreme 3' end of the randomized sequence (class I-A) contained an A at this position (see above). The compensatory A25U:U30A change did not restore binding and was similar to the A25U single mutant. Taken together, these data indicate that the primary sequence makes a more substantial contribution to CsrA binding than does the ability to form a hairpin with the GGA motif in the loop. The finding that CsrA can bind to natural RNA targets that are not contained with in RNA hairpins is consistent with this conclusion (Dubey et al. 2003Go).



View larger version (19K):
[in this window]
[in a new window]
 
FIGURE 5. Mutation analysis of a selected high-affinity RNA ligand. Predicted RNA secondary structure of the R9–43 progenitor wild-type (WT) transcript and the position and identity of each mutation are shown. The A23U, U32A, C24G, G31C, A25U, and U30A mutations disrupt base pairing. Equilibrium binding constants (Kd) for CsrA– RNA interaction ± standard deviation, as well as the Kd mutant/Kd R9–43 (WT) ratios, are indicated.

 


View larger version (73K):
[in this window]
[in a new window]
 
FIGURE 6. Gel mobility shift analysis of CsrA complexed with RNA ligand R9–43 and mutant derivatives U32A, U30A, and {Delta}GGA. 5' End-labeled RNA was incubated with CsrA at the concentration indicated at the bottom of each lane. Positions of free (F) and bound (B) RNAs are shown. Binding curves are shown for R9–43 (WT), U32A, and U30A.

 
In addition to changes within the stem, the affect of nucleotide substitutions within the loop was examined. Alteration of either G in the conserved GGA motif (G27C and G28C) led to a 20-fold reduction in binding affinity, whereas substitution of the A (A29C) resulted in a fivefold reduction in Kd (Fig. 5Go). Not surprisingly, deletion of the GGA motif resulted in a severe binding defect (>140-fold) (Figs. 5Go, 6Go). These results confirm the importance of the GGA motif in CsrA binding. Interestingly, the A26G substitution, which generated a GNRA tetraloop (GGGA), resulted in a ninefold reduction in binding affinity (Fig. 5Go). Recall that the SELEX-derived consensus (Fig. 2Go), as well as previous binding data (Table 1Go), suggested that either purine would be equally tolerated at this position. However, as noted above, only one of the selected RNA ligands with a four-base loop contained a G at this position, thereby generating a GNRA tetraloop. The finding that a GNRA tetraloop reduces CsrA binding provides an explanation for the scarcity of selected ligands in which the conserved GGA motif is part of a tetraloop.

Model of CsrA–RNA interaction
CsrA binding sites have been identified in the ncRNAs CsrB (18 apparent sites) and CsrC (nine apparent sites), as well as in glgC (four sites), cstA (four sites), and pgaA (six sites). The majority of the binding sites in CsrB and CsrC are predicted to contain the GGA motif in loops of short hairpins, although the CsrA binding site consensus sequence for these ncRNAs is CAGGAUG, which has a 1-nucleotide (nt) deletion relative to the SELEX-derived consensus (Liu et al. 1997Go; Weilbacher et al. 2003Go). Of the 10 binding sites in the glgC and pgaA transcripts, five are predicted to have a similar structural arrangement. However, none of the GGA motifs within the cstA transcript are predicted to be in the loops of hairpins (Baker et al. 2002Go; Dubey et al. 2003Go; Wang et al. 2005Go). Thus, while it is common for natural CsrA binding sites to contain both conserved sequence and structural arrangements, it is evident that the common structural arrangement is not essential for CsrA binding. It is also important to point out that not every natural CsrA binding site has a GGA motif. In four of the CsrB binding sites, the GGA is replaced with a GGG (Liu et al. 1997Go), while GGA is replaced with AGA in one of the pgaA binding sites (Wang et al. 2005Go).

Of the 55 selected CsrA ligands reported in this paper, 51 have their GGA motif in the loop of a hairpin within the most stable predicted structure. By design, all of these ligands contained a single CsrA binding site. A likely explanation for the high degree of structural conservation among the SELEX-derived CsrA targets is that having the GGA motif in the loop of a hairpin allows for high-affinity interaction. Furthermore, the finding that the AC residues of the 100% conserved ACA motif were predicted to pair with the conserved GU residues in all but one selected ligand strongly suggests that primary sequence information is embedded within the stem of each hairpin. This conclusion is supported by the mutational analysis of the R9–43 ligand, which demonstrated that the identity of the ACA and GU residues was more important for CsrA binding than the simple ability to form base pairs (Fig. 5Go). However, inspection of CsrB and CsrC stems does not reveal conserved primary sequence. Thus, the sites in these ncRNAs do not appear to be optimized for high-affinity binding. Perhaps the absence of the conserved stem sequences in CsrB and CsrC prevent them from competing too effectively for CsrA.

A reasonable model for how CsrA binds to RNA is as follows. CsrA initially interacts with the unpaired loop sequences when the GGA motif is present within the loop of a hairpin. Following this initial interaction, the hairpin is partially melted, leading to additional base-specific contacts. In situations where an RNA secondary structure is not present, CsrA is still able to bind but with lower affinity due to a reduction in the association rate. Consistent with this notion is the fact that of the five naturally occurring RNA targets that have been characterized, CsrA has the lowest affinity for cstA, the only known CsrA target in which none of the GGA motifs are predicted to be present in the loops of hairpins. Finally, one has to take into consideration that CsrA functions as a homodimer (Dubey et al. 2003Go). Thus, the cooperative interaction observed with CsrB, CsrC, cstA, and pgaA likely involves both protein– RNA and dimer–dimer interactions. However, it remains to be determined whether CsrA dimers interact with one or two binding sites. Engineering multiples of the SELEX-derived consensus CsrA binding site into a single transcript will facilitate future studies on cooperative interaction and CsrA–RNA stoichiometry.


    MATERIALS AND METHODS
 TOP
 ABSTRACT
 INTRODUCTION
 RESULTS AND DISCUSSION
 MATERIALS AND METHODS
 REFERENCES
 
Purification of CsrA-H6 protein
Plasmid pCSB12 contains csrA cloned into the NdeI and BamHI sites of pET21a+ (Inovagen). This plasmid was designed to allow production of CsrA containing six additional His residues at the C terminus (CsrA-H6). pCSB12 was transformed into BL21 (DE3) pLysE (Invitrogen) to allow overproduction of CsrA-H6. The resulting strain (PLB328) was grown at 37° C in terrific broth to an OD600 of 0.6, at which time 1 mM IPTG was added to the culture and growth was continued for 20 h. Cells were harvested by centrifugation at 4° C and the cell pellet was frozen at –20° C until further use. Cell pellets (9 g) were thawed and suspended in 45 mL lysis buffer (50 mM NaH2PO4, pH 8.0, 300 mM NaCl, 10 mM imidazole, 10% glycerol). Cell lysates were prepared by sonication, followed by centrifugation at 10,000g for 30 min at 4° C. The supernatant was mixed with 1 mL of prewashed Ni-NTA resin (Qiagen) for 1 h at 4° C and subsequently packed into a 1-mL column. The column was washed once with 20 mL of wash buffer (50 mM NaH2PO4, pH 8.0, 300 mM NaCl, 10% glycerol) containing 20 mM imidazole and twice with wash buffer containing 50 mM imidazole. CsrA-H6 was eluted with 15 mL wash buffer containing 250 mM imidazole. Column fractions were analyzed by 15% SDS-PAGE and Coomassie Blue staining. Fractions containing highly purified CsrA-H6 were pooled and dialyzed in phosphocellulose column buffer (100 mM Tris-HCl, pH 7.0, 100 mM NaCl, 10% glycerol). The dialysate was mixed with phosphocellulose pre-equilibrated with phosphocellulose column buffer and then packed into a 5-mL column. CsrA eluted from the column between 250 and 750 mM NaCl. Column fractions were analyzed by 15% SDS-PAGE and Coomassie Blue staining. Fractions containing pure CsrA-H6 were combined and dialyzed against 10 mM Tris-HCl, pH 7.0, 100 mM KCl, 10 mM MgCl2 and 25% glycerol. The CsrA-H6 concentration was estimated using the Bio-Rad protein assay.

In vitro selection of RNA ligands
SELEX was performed by modifying published procedures (Tuerk and Gold 1990Go; Ulrich et al. 2002Go; Vo et al. 2003Go). The oligonucleotide used in this analysis (5'-ACCGAGTCCAGAAGCTTGTAG TAC(N15)GCCTAGATGGAGTTGAATTCTCCCTATAGTGAGT CGTATTAC-3') contained a 15-nt randomized region (N15) flanked on both sides by constant sequences. To create the initial pool of randomized DNA templates, 10 pmol of the oligonucleotide containing the randomized sequence was used as a template for PCR amplification using 240 pmol each of primer 1 (5'- GTAATACGACTCACTATAGGGAGAATTCAACTCCATCTA- 3') and primer 2 (5'-ACCGAGTCCAGAAGCTTGTAGT-3'). After gel purification of the PCR product, RNA was synthesized using the in vitro Megascript transcription kit (Ambion) in the presence of 0.66 µM [{alpha}32P]ATP. The reaction mixture was then treated with five units of DNase I for 15 min to remove template DNA. Labeled transcripts were gel purified on 8% denaturing polyacrylamide gels and quantified.

RNA suspended in TE was renatured by heating to 85° C followed by slow cooling to room temperature. Binding reactions (500 µL) contained various concentrations of CsrA-H6 and RNA in binding buffer (10 mM Tris-HCl, pH 7.5, 10 mM MgCl2, 100 mM KCl, 0.5 mg of yeast RNA, and 7.5% glycerol). After incubation for 30 min at 37°C, 40 µL of Ni-NTA slurry prewashed with binding buffer was added to the reaction mixture and further incubated for 10 min with repeated mixing. Bound RNAs were recovered by modifying a published procedure (Vo et al. 2003Go). The beads were washed three times with 100 µL of binding buffer and bound RNAs were eluted with 100 µL of binding buffer containing 2 M imidazole. The RNA was then extracted with phenol and phenol-chloroform and ethanol precipitated. RNAs were reverse transcribed using AMV reverse transciptase (Roche) and primer 2 and PCR amplified with primers 1 and 2. In vitro generated transcripts using the RT-PCR templates were subjected to the next round of selection. The CsrA-H6 concentration was 500 nM in the first two rounds of SELEX, 250 nM in rounds 3 and 4, 125 nM in rounds 5 through 7, 62.5 nM in rounds 8 and 9, 31 nM in round 10, 16 nM in round 11, and 8 nM in round 12. In each round the RNA concentration was five- to 10-fold higher than the concentration of CsrA-H6. RT-PCR products from rounds 4, 7, 9, and 12 were cloned and sequenced. A total of 55 clones were sequenced and this information was used to derive a consensus sequence (Fig. 2Go).

Gel mobility shift assays
Quantitative gel mobility shift assays were used to measure the affinity of CsrA-RNA binding by following published procedures (Baker et al. 2002Go; Dubey et al. 2003Go). For individually cloned DNA sequences, RNA was synthesized in vitro with the Ambion MEGAscript kit and linearized plasmid DNA as templates. For gel mobility shift assays with RNA pools, the transcripts from rounds 0, 2, 4, 7, and 9 were used. Gel-purified RNAs and RNA pools from various rounds were dephosphorylated with calf intestinal alkaline phosphatase and subsequently 5' end-labeled using [{gamma}-32P]ATP and polynucleotide kinase. Labeled RNAs were gel purified, suspended in TE (10 mM Tris-HCl, pH 8.0, 1 mM EDTA) and renatured by heating to 85°C and slow cooling to room temperature. Binding reaction mixtures (10 µL) contained 10 mM Tris-HCl, pH 7.5, 100 mM KCl, 10 mM MgCl2, 32.5 ng of yeast RNA, 10% glycerol, 20 mM dithiothreitol, 4 U of RNase inhibitor (Ambion), 5' end-labeled RNAs (30 pM or 250 pM), various concentrations of purified CsrA (Baker et al. 2002Go), and 0.1 mg/mL xylene cyanol. Competition assays also contained unlabeled RNA competitor (see text for details). Reaction mixtures were incubated for 30 min at 37°C to allow CsrA–RNA complex formation. Samples were then fractionated through native polyacrylamide gels. Radioactive bands were visualized with a phosphorimager (Molecular Dynamics). Free and bound RNA species were quantified using ImageQuant software (Molecular Dynamics), and the apparent equilibrium binding constant (Kd) and cooperativity coefficient (n) for CsrA–RNA complex formation was calculated as described (Yakhnin et al. 2000Go).

Site-directed mutagenesis
Effects of mutations on CsrA binding were investigated by mutating various conserved residues of a selected RNA ligand (R9–43) using QuikChange mutagenesis (Stratagene). DNA corresponding to each mutation was cloned into plasmid pTZ18U (Stratagene). The in vitro transcription reactions, gel mobility shift assays, and determination of the binding constants were carried out as described above.


    ACKNOWLEDGMENTS
 
We thank Nam V. Vo for technical advice with SELEX. This work was supported by grant GM059969 from the National Institutes of Health.


    Footnotes
 
3 Present address: Tata Chemicals Ltd., Leela Business Park, Andheri-Kurla Rd., Mumbai 400059, India. Back

Article published online ahead of print. Article and publication date are at http://www.rnajournal.org/cgi/doi/10.1261/rna.2990205.

Received May 2, 2005; accepted July 6, 2005.


    REFERENCES
 TOP
 ABSTRACT
 INTRODUCTION
 RESULTS AND DISCUSSION
 MATERIALS AND METHODS
 REFERENCES
 

Altier, C., Suyemoto, M., and Lawhon, S.D. 2000. Regulation of Salmonella enterica serovar Typhimurium invasion genes by csrA. Infect. Immun. 68: 6790–6797.[Abstract/Free Full Text]

Baker, C.S., Morozov, I., Suzuki, K., Romeo, T., and Babitzke, P. 2002. CsrA regulates glycogen biosynthesis by preventing translation of glgC in Escherichia coli. Mol. Microbiol. 44: 1599–1610.[CrossRef][Medline]

Baumann, C., Xirasagar, S., and Gollinick, P. 1997. The trp RNA-binding attenuation protein (TRAP) from Bacillus subtilis binds to unstacked trp leader RNA. J. Biol. Chem. 272: 19863–19869.[Abstract/Free Full Text]

Blumer, C., Heeb, S., Pessi, G., and Haas, D. 1999. Global GacA-steered control of cyanide and exoprotease production in Pseudomonas fluorescens involves specific ribosome binding sites. Proc. Natl. Acad. Sci. 96: 14073–14078.[Abstract/Free Full Text]

Cui, Y., Mukherjee, A., Dumenyo, C.K., Liu, Y., and Chatterjee, A.K. 1999. rsmC of the soft-rotting bacterium Erwinia carotovora subsp. carotovora negatively controls extracellular enzyme and hairpin (Ecc) production and virulence by modulating levels of regulatory RNA (rsmB) and RNA-binding protein (RsmA). J. Bacteriol. 181: 6042–6052.[Abstract/Free Full Text]

Du, H., Yakhnin, A. V., Dharmaraj, S., and Babitzke, P. 2000. trp RNA-binding attenuation protein (TRAP)-5' stem-loop RNA interaction is required for proper transcription attenuation control of the Bacillus subtilis trpEDCFBA operon. J. Bacteriol. 182: 1819–1827.[Abstract/Free Full Text]

Dubey, A.K., Baker, C.S., Suzuki, K., Jones, A.D., Pandit, P., Romeo, T., and Babitzke, P. 2003. CsrA regulates translation of the Escherichia coli starvation gene, cstA, by blocking ribosome access to the cstA transcript. J. Bacteriol. 185: 4450–4460.[Abstract/Free Full Text]

Ellington, A.D. and Szostak, J.W. 1990. In vitro selection of RNA molecules that bind specific ligands. Nature 346: 818–822.[CrossRef][Medline]

Fettes, P.S., Forsbach-Birk, V., Lynch, D., and Marre, R. 2001. Over-expression of a Legionella pneumophila homologue of the E. coli regulator CsrA affects cell size, flagellation, and pigmentation. Int. J. Med. Microbiol. 291: 353–360.[CrossRef][Medline]

Gudapaty, S., Suzuki, K., Wang, X., Babitzke, P., and Romeo, T. 2001. Regulatory interactions of Csr components: The RNA binding protein CsrA activates csrB transcription in Escherichia coli. J. Bacteriol. 183: 6017–6027.[Abstract/Free Full Text]

Jackson, D.W., Suzuki, K., Oakford, L., Simecka, J.W., Hart, M.E., and Romeo, T. 2002. Biofilm formation and dispersal is under the influence of the global regulator CsrA of Escherichia coli. J. Bacteriol. 184: 290–301.[Abstract/Free Full Text]

Liaw, S.J., Lai, H.C., Ho, S.W., Luh, K.T., and Wang, W.B. 2003. Role of RsmA in the regulation of swarming motility and virulence factor expression in Proteus mirabilis. J. Med. Microbiol. 52: 19–28[Abstract/Free Full Text]

Liu, M.Y., Yang, H., and Romeo, T. 1995. The product of the pleiotropic Escherichia coli gene csrA modulates glycogen biosynthesis via effects on mRNA stability. J. Bacteriol. 177: 2663–2672.[Abstract/Free Full Text]

Liu, M.Y., Gui, G., Wei, B., Preston, J.F. 3rd, Oakford, L., Yuksel, U., Giedroc, D.P., and Romeo, T. 1997. The RNA molecule CsrB binds to the global regulatory protein CsrA and antagonizes its activity in Escherichia coli. J. Biol. Chem. 272: 17502–17510.[Abstract/Free Full Text]

Lozupone, C., Changayil, S., Majerfeld, I., and Yarus, M. 2003. Selection of the simplest RNA that binds isoluecine. RNA 9: 1315–1322.[Abstract/Free Full Text]

Mathews, D.H., Sabina, J., Zuker, M., and Turner, D.H. 1999. Expanded sequence dependence of thermodynamic parameters improves prediction of RNA secondary structure. J. Mol. Biol. 288: 911–940.[CrossRef][Medline]

Reimmann, C., Valverde, C., Kay, E., and Haas, D. 2005. Post-transcriptional repression of GacS/GacA-controlled genes by the RNA-binding protein RsmE acting together with RsmA in the biocontrol strain Pseudomonas fluorescens CHA0. J. Bacteriol. 187: 276–285.[Abstract/Free Full Text]

Romeo, T. 1998. Global regulation by the small RNA-binding protein CsrA and the non-coding RNA molecule CsrB. Mol. Microbiol. 29: 1321–1330.[CrossRef][Medline]

Romeo, T., Gong, M., Liu, M.Y., and Brun-Zinkernagel, A.M. 1993. Identification and molecular characterization of csrA, a pleiotropic gene from Escherichia coli that affects glycogen biosynthesis, gluconeogenesis, cell size, and surface properties. J. Bacteriol. 175: 4744–4755.[Abstract/Free Full Text]

Sabnis, N.A., Yang, H., and Romeo, T. 1995. Pleiotropic regulation of central carbohydrate metabolism in Escherichia coli via the gene csrA. J. Biol. Chem. 270: 29096–29104.[Abstract/Free Full Text]

Schneider, D., Tuerk, C., and Gold, L. 1992. Selection of high affinity RNA ligands to the bacteriophage R17 coat protein. J. Mol. Biol. 228: 862–869.[CrossRef][Medline]

Schultz, J.E. and Matin, A. 1991. Molecular and functional characterization of a carbon starvation gene of Escherichia coli. J. Mol. Biol. 218: 129–140.[CrossRef][Medline]

Sudershana, S., Du, H., Mahalanabis, M., and Babitzke, P. 1999. A 5' stem-loop participates in the transcription attenuation mechanism that controls expression of the Bacillus subtilis trpEDCFBA operon. J. Bacteriol. 181: 5742–5749.[Abstract/Free Full Text]

Tuerk, C. and Gold, L. 1990. Systematic evolution of ligands by exponential enrichment: RNA ligands to bacteriophage T4 DNA polymerase. Science 249: 505–510.[Abstract/Free Full Text]

Ulrich, H., Magdesian M.H., Alves M.J., and Colli, W. 2002. In vitro selection of RNA aptamers that bind to cell adhesion receptors of Trypanosoma cruzi and inhibit cell invasion. J. Biol. Chem. 277: 20756–20762.[Abstract/Free Full Text]

Vo, N.V., Oh, J.W., and Lai, M.M.C. 2003. Identification of RNA ligands that bind hepatitis C virus polymerase selectively and inhibit its RNA synthesis from the natural viral RNA templates. Virology 307: 301–316.[CrossRef][Medline]

Wang, X., Dubey, A.K., Suzuki, K., Baker, C.S., Babitzke, P., and Romeo, T. 2005. CsrA post-transcriptionally represses pgaABCD, responsible for synthesis of a biofilm polysaccharide adhesin of Escherichia coli. Mol. Microbiol. 56: 1648–1663.[Medline]

Wei, B., Shin, S., LaPorte, D., Wolfe, A.J., and Romeo, T. 2000. Global regulatory mutations in csrA and rpoS cause severe central carbon stress in Escherichia coli in the presence of acetate. J. Bacteriol. 182: 1632–1640.[Abstract/Free Full Text]

Wei, B.L., Brun-Zinkernagel, A.M., Simecka, J.W., Pruss, B.M., Babitzke, P., and Romeo, T. 2001. Positive regulation of motility and flhDC expression by the RNA-binding protein CsrA of Escherichia coli. Mol. Microbiol. 40: 245–256.[CrossRef][Medline]

Weilbacher, T., Suzuki, K., Dubey, A.K., Wang, X., Gudapaty, S., Morozov, I., Baker, C.S., Georgellis, D., Babitzke, P., and Romeo, T. 2003. A novel sRNA component of the carbon storage regulatory system of Escherichia coli. Mol. Microbiol. 48: 657–670.[CrossRef][Medline]

Yakhnin, A.V., Trimble, J.J., Chiaro, C.R., and Babitzke, P. 2000. Effects of mutations in the L-tryptophan binding pocket of the trp RNA-binding attenuation protein of Bacillus subtilis. J. Biol. Chem. 275: 4519–4524.[Abstract/Free Full Text]

Yang, H., Liu, M.Y., and Romeo, T. 1996. Coordinate genetic regulation of glycogen catabolism and biosynthesis in Escherichia coli via the CsrA gene product. J. Bacteriol. 178: 1012–1017.[Abstract/Free Full Text]

Zuker, M. 2003. Mfold web server for nucleic acid folding and hybridization prediction. Nucleic Acids Res. 31: 3406–3415.[Abstract/Free Full Text]
Add to CiteULike CiteULike   Add to Connotea Connotea   Add to Del.icio.us Del.icio.us   Add to Digg Digg   Add to Reddit Reddit   Add to Technorati Technorati    What's this?


This article has been cited by other articles:


Home page
Infect. Immun.Home page
J. H. Carlson, W. M. Whitmire, D. D. Crane, L. Wicke, K. Virtaneva, D. E. Sturdevant, J. J. Kupko III, S. F. Porcella, N. Martinez-Orengo, R. A. Heinzen, et al.
The Chlamydia trachomatis Plasmid Is a Transcriptional Regulator of Chromosomal Genes and a Virulence Factor
Infect. Immun., June 1, 2008; 76(6): 2273 - 2283.
[Abstract] [Full Text] [PDF]


Home page
J. Bacteriol.Home page
T. L. Campbell and E. D. Brown
Genetic Interaction Screens with Ordered Overexpression and Deletion Clone Sets Implicate the Escherichia coli GTPase YjeQ in Late Ribosome Biogenesis
J. Bacteriol., April 1, 2008; 190(7): 2537 - 2545.
[Abstract] [Full Text] [PDF]


Home page
MicrobiologyHome page
C. Lucchetti-Miganeh, E. Burrowes, C. Baysse, and G. Ermel
The post-transcriptional regulator CsrA plays a central role in the adaptation of bacterial pathogens to different stages of infection in animal hosts
Microbiology, January 1, 2008; 154(1): 16 - 29.
[Abstract] [Full Text] [PDF]


Home page
J. Bacteriol.Home page
C. S. Baker, L. A. Eory, H. Yakhnin, J. Mercante, T. Romeo, and P. Babitzke
CsrA Inhibits Translation Initiation of Escherichia coli hfq by Binding to a Single Site Overlapping the Shine-Dalgarno Sequence
J. Bacteriol., August 1, 2007; 189(15): 5472 - 5481.
[Abstract] [Full Text] [PDF]


Home page
J. Bacteriol.Home page
Y. Sakuragi and R. Kolter
Quorum-Sensing Regulation of the Biofilm Matrix Genes (pel) of Pseudomonas aeruginosa
J. Bacteriol., July 15, 2007; 189(14): 5383 - 5386.
[Abstract] [Full Text] [PDF]


Home page
J. Bacteriol.Home page
H. Chouayekh, H. Nothaft, S. Delaunay, M. Linder, B. Payrastre, N. Seghezzi, F. Titgemeyer, and M. J. Virolle
Phosphoinositides Are Involved in Control of the Glucose-Dependent Growth Resumption That Follows the Transition Phase in Streptomyces lividans
J. Bacteriol., February 1, 2007; 189(3): 741 - 749.
[Abstract] [Full Text] [PDF]


Home page
J. Biol. Chem.Home page
J. Mercante, K. Suzuki, X. Cheng, P. Babitzke, and T. Romeo
Comprehensive Alanine-scanning Mutagenesis of Escherichia coli CsrA Defines Two Subdomains of Critical Functional Importance
J. Biol. Chem., October 20, 2006; 281(42): 31832 - 31842.
[Abstract] [Full Text] [PDF]


Home page
Nucleic Acids ResHome page
M. Hiller, R. Pudimat, A. Busch, and R. Backofen
Using RNA secondary structures to guide sequence motif finding towards single-stranded regions
Nucleic Acids Res., October 18, 2006; 34(17): e117 - e117.
[Abstract] [Full Text] [PDF]


Home page
Genes Dev.Home page
K. Suzuki, P. Babitzke, S. R. Kushner, and T. Romeo
Identification of a novel regulatory protein (CsrD) that targets the global regulatory RNAs CsrB and CsrC for degradation by RNase E
Genes & Dev., September 15, 2006; 20(18): 2605 - 2617.
[Abstract] [Full Text] [PDF]


Home page
J. Bacteriol.Home page
E. Kay, B. Humair, V. Denervaud, K. Riedel, S. Spahr, L. Eberl, C. Valverde, and D. Haas
Two GacA-Dependent Small RNAs Modulate the Quorum-Sensing Response in Pseudomonas aeruginosa.
J. Bacteriol., August 1, 2006; 188(16): 6026 - 6033.
[Abstract] [Full Text] [PDF]


Home page
Nucleic Acids ResHome page
P. R. Kulkarni, X. Cui, J. W. Williams, A. M. Stevens, and R. V. Kulkarni
Prediction of CsrA-regulating small RNAs in bacteria and their experimental verification in Vibrio fischeri
Nucleic Acids Res., July 5, 2006; 34(11): 3361 - 3369.
[Abstract] [Full Text] [PDF]


This Article
Right arrow Abstract Freely available
Right arrow Full Text (PDF)