Acidic C-terminal domains autoregulate the RNA chaperone Hfq

Abstract
Introduction
Results
Discussion
Materials and methods
References
Article and author information
Metrics

Abstract

The RNA chaperone Hfq is an Sm protein that facilitates base pairing between bacterial small RNAs (sRNAs) and mRNAs involved in stress response and pathogenesis. Hfq possesses an intrinsically disordered C-terminal domain (CTD) that may tune the function of the Sm domain in different organisms. In Escherichia coli, the Hfq CTD increases kinetic competition between sRNAs and recycles Hfq from the sRNA-mRNA duplex. Here, de novo Rosetta modeling and competitive binding experiments show that the acidic tip of the E. coli Hfq CTD transiently binds the basic Sm core residues necessary for RNA annealing. The CTD tip competes against non-specific RNA binding, facilitates dsRNA release, and prevents indiscriminate DNA aggregation, suggesting that this acidic peptide mimics nucleic acid to auto-regulate RNA binding to the Sm ring. The mechanism of CTD auto-inhibition predicts the chaperone function of Hfq in bacterial genera and illuminates how Sm proteins may evolve new functions.

https://doi.org/10.7554/eLife.27049.001

Introduction

Host factor for RNA phage Qβ replication (Hfq) is found in most sequenced bacterial genomes (Sun et al., 2002) and plays a well characterized role in post-transcriptional regulation by small non-coding RNA (sRNA) (Gottesman et al., 2006; Storz et al., 2011). Regulation by Hfq and sRNAs is important for controlling the expression of metabolic, stress-response and virulence genes in many genera (Feliciano et al., 2016). Hfq binds sRNA and facilitates interactions between sRNAs and their mRNA targets (Zhang et al., 2002; Moller et al., 2002). To chaperone sRNA target recognition, Hfq must select its substrates from a large pool of nucleic acid in the cell and efficiently dissociate from its products at the end of each RNA annealing cycle (Rajkowitsch et al., 2007).

E. coli Hfq contains an Sm-like domain (residues 7–65) that oligomerizes into a homohexameric ring with two sequence-specific RNA-binding faces. The proximal face of the ring is highly conserved and binds to uridines (Zhang et al., 2002; Schumacher et al., 2002) at the 3’-ends of bacterial small non-coding RNA (sRNA) that resemble a classic Sm binding site (Zhou et al., 2014). In E. coli and many Gram negative bacteria, the distal face of Hfq binds to AAN triplet repeats (Mikulecky et al., 2004; Link et al., 2009) found in mRNA leaders (Link et al., 2009; Soper et al., 2011) and certain sRNAs (Schu et al., 2015; Małecka et al., 2015). In addition to these sequence-specific RNA binding sites, arginine-rich basic patches at the rim of the E. coli Hfq hexamer interact with the sRNA body (Zhang et al., 2002; Otaka et al., 2011; Sauer et al., 2012; Ishikawa et al., 2012; Zhang et al., 2013) and facilitate annealing with target mRNAs (Panja et al., 2013; Zheng et al., 2016).

Like many RNA binding proteins, Hfq also possesses intrinsically disordered domains that have the potential to modulate the function of the core Sm ring. The E. coli Hfq Sm domain is flanked by a short, disordered, N-terminal domain (NTD; residues 1–6), which protrudes from the proximal face of the hexamer, and a longer disordered C-terminal domain (CTD; residues 66–102), which extends from the rim (Beich-Frandsen et al., 2011a; Vincent et al., 2012). NMR chemical shift perturbations from a comparison of full-length Hfq (Hfq102) and a truncated variant lacking the CTD (Hfq65) suggested that some part of the CTD contacts residues on the rim of the hexamer, although the specificity of these proposed contacts was uncertain since they occur near where the CTD protrudes from the ring (Beich-Frandsen et al., 2011a; Vincent et al., 2012).

The functional importance of the CTD for sRNA regulation has also been unclear, owing to the conflicting results of different studies (Sonnleitner et al., 2004; Olsen et al., 2010; Večerek et al., 2008; Salim et al., 2012). Using a combination of biophysical and genetic approaches, we recently showed that the CTD displaces RNA from the rim and proximal face of Hfq (Santiago-Frangos et al., 2016), with two important consequences. First, release of annealed dsRNA from the arginine-rich rim is accelerated, increasing Hfq turnover. Second, kinetic competition between sRNAs is increased, which allows dominant sRNAs to bind to Hfq and accumulate in the cell, while weaker competitors are degraded (Santiago-Frangos et al., 2016). The latter creates a hierarchy of sRNA regulation that depends on the CTD.

The mechanism by which the CTD displaces RNA from the core (Sm domain and NTD) of Hfq is unknown. No common sequence motifs have been identified in the CTD (Sun et al., 2002; Vincent et al., 2012; Weichenrieder, 2014; Sobrero and Valverde, 2012; Fortas et al., 2015; Updegrove et al., 2016), which varies in length and composition across bacteria (Attia et al., 2008; Schilling and Gerischer, 2009; Baba et al., 2010). This diversity is characteristic of disordered peptides, which rapidly evolve via non-conservative substitutions and indels (Liu et al., 2008; Brown et al., 2010; Light et al., 2013). Two models could explain the displacement of RNA by CTDs in E. coli Hfq. The ‘polymer brush’ model suggests the CTDs passively obstruct RNA binding sites. This model is attractive since it depends only on the length and flexibility of the CTD. In contrast, the ‘nucleic acid mimic’ model suggests that the CTDs specifically bind to basic core residues and actively compete against nucleic acids. Given the divergence of CTD and core sequences, this model predicts that CTD auto-regulation is likely in some Hfq clades but not others.

In this study, we use de novo modeling and biophysical experiments to determine the mechanism by which the CTD regulates Hfq activity. We propose that the acidic CTD tip transiently binds the rim of E. coli Hfq and makes distributed interactions with basic residues, thereby modulating RNA and DNA binding and RNA annealing. Applying our modeling procedure to Hfqs from other bacteria demonstrates that stable interactions between the acidic CTD and the basic rim correlate with the importance of Hfq for sRNA regulation in that host. Thus, our proposed mechanism of CTD auto-inhibition provides a basis for predicting the function of Hfq in different bacteria. Our approach may be useful for predicting the sequence–function relationship of disordered domains in other partially disordered proteins.

Results

C-terminus of Hfq is enriched for acidic residues

To search for conserved features or amino acid motifs amongst the highly heterogeneous Hfq CTDs, we first built a phylogenetic tree (Figure 1—figure supplement 1) from the multiple sequence alignment of nearly 1000 non-redundant sequences (see Methods). The cluster containing E. coli Hfq contained many other Hfq variants previously identified as functional in RNA annealing (Zheng et al., 2016) or sRNA regulation (Gottesman et al., 2006). Therefore, we examined the sequence logo of this cluster of 222 Hfqs in more detail (Figure 1A).

Figure 1 with 3 supplements see all

Download asset Open asset

Acidic residues in the Hfq CTD are predicted to bind basic core.

(A) Sequence logo (Crooks et al., 2004) of the CTD generated from gapped alignment of Hfq sequences that clustered with *Escherichia coli* Hfq (Group 1 in Figure 1—figure supplement 1A,B), numbered according to the *E. coli* sequence. Regions of interest are denoted above. The gapped *E. coli* CTD sequence is shown below. Eukaryotic Sm proteins cluster separately (Figure 1—figure supplement 1D). (B) (Top) Average number of times a given core residue favorably interacts (E < −1.0 Rosetta Energy Units) with at least one acidic CTD residue, per low energy model. Acidic CTD residues most frequently interact with basic Hfq core residues. (Bottom) Mutation of acidic CTD residues 97, 99, 100 and 102 to basic or polar residues decreases the number of predicted core interactions. Error bars represent ± 1 s.d. as computed by bootstrap resampling of the computational models (see Methods and Figure 1—figure supplement 2). Of 36 core residues not predicted to interact with the CTD, 14 had accessible surface area < 2.0 A², computed in PyMOL. (C) (Left) Example low-energy model of wildtype *E. coli* Hfq; top-down proximal view. Light grey, NTD; cyan, Hfq core; pink-purple, CTD; red, CTD tip. (Center) Side view of rim of the same Hfq model. (Inset) Example hydrogen bonding network at the CTD–core binding interface showing interactions between the acidic CTD residues (red) and core residues as indicated. Additional models in Figure 1—figure supplement 3.

https://doi.org/10.7554/eLife.27049.002

The start of the CTD region is delineated by a proline at position 64 of E. coli Hfq that is strongly conserved across all clades. Additionally, an arginine at the beginning of the CTD (position 66 in E. coli) that packs against the lateral edge of the Hfq hexamer (Beich-Frandsen et al., 2011a; Sauter et al., 2003; Dimastrogiovanni et al., 2014) is strongly conserved. Although the middle linker region of the CTD lacks conserved motifs (Figure 1—figure supplement 1), the C-terminus is rich in acidic residues, corresponding to the sequence DSEETE in E. coli. Noting that most Hfq clusters containing a basic patch on the rim also end in acidic residues, we hypothesized that the CTD tip binds the rim. Because the basic patch is essential for sRNA binding and annealing, direct interaction between the CTD tip and the Hfq core could explain the previously observed auto-inhibition of the CTD (Santiago-Frangos et al., 2016).

De novo modeling of CTD interactions in the Hfq hexamer

To determine whether the acidic tip of the E. coli Hfq CTD could interact with basic residues in the core, we used Rosetta FloppyTail (Kleiger et al., 2009), a de novo modeling approach for disordered regions of proteins. We updated the original FloppyTail algorithm to model multiple disordered regions simultaneously and to ensure adequate sampling of backbone degrees of freedom (see Materials and methods and Figure 1—figure supplement 2). Then, we generated and analyzed ~30,000 models of the full-length E. coli Hfq hexamer. In the lowest energy (1%) subset of models, the acidic CTD residues (D97, E99, E100, and E102) frequently interact with basic residues on the rim (R16, R17, R19 and K47) and in the NTD (K3) (Figure 1B, top). By contrast, K31 on the distal face is not predicted to be contacted by the CTD, although K31 is highly accessible. This bias accords with prior observations that the CTD does not displace RNA from Hfq’s distal face (Santiago-Frangos et al., 2016). As anticipated for a disordered domain, no single conformation dominated the ensemble of models (Figure 1—figure supplement 3). Rather, the acidic CTD tip was found to bind to various combinations of residues in the basic patch (Figure 1C, inset).

To confirm we were not simply observing the non-specific collapse of the disordered CTD onto the core, or enriching interactions between highly solvent-accessible polar residues, we repeated our simulations using a mutant Hfq in which the acidic CTD residues were replaced with polar or basic side chains (D97R-E99N-E100K-E102N). These mutations drastically decreased the frequency of predicted interactions between the basic core residues and CTD residues at positions 97, 99, 100 and 102 in our simulations (Figure 1B, bottom), without increasing predicted interactions between this mutant CTD and solvent-accessible acidic residues on the Hfq core (D9, E18, E37 and D40).

Acidic CTD specifically binds Hfq rim

To determine whether the CTD interacts with the rim as predicted by our models, we used fluorescence anisotropy to measure the affinity of core Hfq (Hfq65) for a fluorescently-labeled CTD peptide, CTD-FITC (Figure 2A and Figure 2—figure supplement 1). CTD-FITC lacks residues 65–72 to avoid contributions to binding from this region, which packs against the Sm domain as one strand of the β-sheet (Arluison et al., 2004). Hfq65 bound to CTD-FITC with a K_d of 2.9 µM Hfq monomer in low salt buffer (cyan in Figure 2B) and 22 µM in a higher salt buffer (Figure 2—figure supplement 2). These interaction strengths are meaningful even at higher ionic strength, because the effective concentration of each individual acidic CTD tip is roughly 350 μM in the full-length protein (see Materials and methods).

Figure 2 with 3 supplements see all

Download asset Open asset

CTD and RNA occupy the same rim binding site.

(A) Scheme for in vitro binding of fluorescent CTD-FITC peptide by Hfq core. The CTD linker is shown in purple, the acidic tip in red and the N-terminal FITC as a yellow star. Hfq core is shown in cyan, with basic rim patches in dark blue. (B) Binding of CTD-FITC to variants of Hfq65 core at 30 ˚C. 45 nM CTD-FITC was titrated with 0–100 µM Hfq monomer in duplicate, and the average (±s.d.) was fit to Equation 5 (Materials and methods). (C) Reaction scheme for annealing an RNA molecular beacon to a target RNA (open bar) (Hopkins et al., 2011; Panja et al., 2015). (D) Progress curves for annealing 50 nM molecular beacon and 100 nM Target by 50 nM Hfq65 hexamer at 30°C, measured by stopped-flow fluorescence. See Figure 2—figure supplement 3 for further data. (E) Contribution of core residues to CTD binding. Interaction energy (Expected Energetic Contribution; EEC) in silico for a core residue in the Rosetta models (solid symbols and solid line; adjusted R² = 0.77) or the average annealing rates for Target and Target-A18 relative to Hfq65 (open symbols and dashed line; adjusted R² = 0.94) *versus* experimental CTD binding energy (∆∆G°) for each Hfq65 variant. The binding energy, ∆∆G° = –RTln(*K_d*^MUT/*K_d*), reflects the perturbation to CTD binding by a mutation in Hfq65. The interaction energy in silico or EEC is defined as the average energy of a tail–core interaction multiplied by the average number of tail–core interactions per model (Figure 1B, top and Equation 3). The relative annealing rate for Hfq65 variants, k_rel = k_obs^MUT/k_obs^WT, is <1 if the mutated residue is important for RNA annealing.

https://doi.org/10.7554/eLife.27049.006

Binding of the CTD-FITC peptide to Hfq65 was weakened by mutations in the basic rim residues R16A, R19D and K47A (Figure 2B), which frequently interact with the CTD in our computational models (Figure 1B). Although we were not able to test the predicted interactions between NTD K3 and the CTD (see Materials and methods), K3 is also known to form electrostatic interactions with the RNA backbone (Dimastrogiovanni et al., 2014). In contrast, mutation of a surface-accessible polar residue (Q35A) close to the binding interface (Figure 1C, inset), slightly enhanced CTD binding (Figure 2B). Intriguingly, A35 is common in Hfq from γ-proteobacteria. Finally, a CTD peptide containing the mutated C-terminal tip (RSNKTN) was not able to bind Hfq65, confirming that the acidic residues on the CTD peptide are necessary for this interaction (grey in Figure 2B).

CTD-bound core residues play a role in RNA annealing

To determine how much core residues that bind the CTD contribute to Hfq’s RNA annealing activity, we compared the effect of rim mutations on the rate of base pairing between an RNA molecular beacon and the 16 nt Target RNA by stopped-flow fluorescence spectroscopy (Figure 2C) (Hopkins et al., 2011). In the absence of competition from the CTD, the rate of annealing in this assay depends only on interactions between the two RNAs and the Hfq core. As previously observed (Santiago-Frangos et al., 2016), Hfq65 is highly active in single-turnover annealing assays (Figure 2—figure supplement 3A). The observed annealing rate was most diminished by the loss of basic residues, especially the conserved R16A, and relatively unaffected by the mutation Q35A (Figure 2D). Similar results were obtained with Target-A18, which anchors to the distal face (Figure 2—figure supplement 3B). The average relative annealing rates of Hfq65 variants correlated well with the importance of each residue for CTD binding in vitro (Figure 2E), suggesting that the CTD peptide and the RNA interact with the same residues on the rim of Hfq.

The predictive value of our computational approach was validated by a direct correlation between the experimentally measured contribution (∆∆G°) of each core residue for CTD binding with the predicted Expected Energetic Contribution (EEC) of that core residue to interactions with the acidic CTD in silico (solid symbols and solid line, Figure 2E). EEC is defined as the average energy of a tail–core interaction multiplied by the average number of tail–core interactions per model. The absolute binding and simulated interaction energies cannot be directly compared because the peptide binding assay is performed in trans rather than in cis, and the Rosetta Energy does not account for entropic contributions to binding. Nevertheless, amino acids that most strongly impacted the free energy of CTD binding when mutated, also had larger contributions to CTD binding in silico (solid symbols and solid line, Figure 2E; linear regression p-value=0.078), and had stronger effects on Hfq65 RNA annealing activity in vitro (open symbols and dashed line, Figure 2E; linear regression p-value=0.020).

Nucleic acids compete with the CTD for binding the Hfq core

If the CTD peptide and RNA interact with the same basic Hfq residues, direct competition between the two could explain how the CTD triggers the release of annealed RNAs from Hfq (Santiago-Frangos et al., 2016) and why it increases the stringency of RNA or DNA binding. To examine whether nucleic acids are in competition against the CTD for binding to the Hfq core, we compared the ability of different nucleic acids to displace CTD-FITC from a preformed Hfq65·CTD-FITC complex (Figure 3A). The strength of competition was expressed as the concentration of nucleic acid needed to displace 50% of CTD-FITC from Hfq65, IC₅₀.

Figure 3

Download asset Open asset

Nucleic acids displace the CTD from the Hfq core.

(A) Scheme for *in vitro* competition of the CTD-FITC·Hfq65 interaction by nucleic acids. (B) Example titrations of the CTD-FITC·Hfq65 complex with sRNA. Titrations were done in duplicate, and the averages (±s.d.) fit to Equation 6 (Materials and methods) to determine the IC_50.(C) IC₅₀ values for nucleic acids of different lengths (blue, RNA oligomers; orange, DNA oligomers; green, sRNAs). From shortest to longest: Target, A18, minRCRB, Target-U6, Target-A18, DNA1, DNA2, DNA2c, dsDNA2, ChiX, RyhB, DsrA, RprA). Dashed circles indicate nucleic acids that deviate from the linear relationship between log IC₅₀ and length (adjusted R² = 0.85).

https://doi.org/10.7554/eLife.27049.010

Natural sRNAs strongly competed against the CTD-FITC peptide, in keeping with their strong (~10 nM) affinity for Hfq (Figure 3B). Short RNA and DNA oligomers that bind Hfq weakly were poorer competitors than natural sRNAs, as expected. In general, competition against the CTD peptide correlated with nucleic acid length, suggesting little sequence specificity, or that longer nucleic acids may occupy more than one basic patch (Figure 3C). However, minRCRB RNA, DNA1 and dsDNA2 deviated from this linear trend (dashed circles; Figure 3C). minRCRB, a stronger than expected competitor, consists of a stem-loop with a 5’-overhang, and has been shown to specifically bind to the rim (Santiago-Frangos et al., 2016; Dimastrogiovanni et al., 2014). Similarly, DNA1 possesses a stable minRCRB-like motif at its 5-end and was also a stronger competitor than expected based on its length. By contrast, the completely double-stranded dsDNA2 was a weaker competitor than expected, consistent with the CTD’s ability to displace annealed dsRNAs from the rim of Hfq (Santiago-Frangos et al., 2016).

Higher local concentration of acidic residues increases CTD autoinhibition

The results above indicate that an interaction between the acidic CTD tip and the basic rim inhibits RNA binding to the basic patch and stimulates release of dsRNA. If this model is correct, shortening the CTD should increase the local concentration of the acidic CTD tip and exacerbate autoinhibition. Alternatively, if the CTD acts as a polymer brush, a shorter CTD should exhibit less autoinhibition because it will exclude less volume around the Hfq core. To test these predictions, we generated the mutant Hfq-sCTD, which lacks residues 86–96 (inclusive). Truncation of this non-conserved ‘linker’ region is predicted to increase the local concentration of the acidic tip around the Hfq core roughly three-fold, from roughly 350 μM per CTD in Hfq102, to ~1220 μM per CTD in Hfq-sCTD (Equation 4 in Materials and methods).

As previously shown (Santiago-Frangos et al., 2016), Hfq65, which lacks the CTD entirely, anneals the 16 nt Target and Target-U6 RNA about five times faster than full-length Hfq102 (100–60-fold vs. 20–10-fold; cyan and black in Figure 4A). In our model, this is because the CTD sweeps RNAs from the rim and proximal face of Hfq. By contrast, both proteins accelerate Target-A18 annealing roughly 100-fold compared to no Hfq, because this RNA remains anchored to the distal face of Hfq and resists CTD displacement (Santiago-Frangos et al., 2016). Importantly, the shortened CTD linker (Hfq-sCTD) decreased annealing rates relative to Hfq102 for all RNA targets (Figure 4A and Figure 4—figure supplement 1), suggesting that access to the basic patch was more restricted.

Figure 4 with 2 supplements see all

Download asset Open asset

Shortening the CTD increases auto-inhibition by acidic residues.

(A) Observed annealing rate constants for 100 nM target RNA and 50 nM molecular beacon with 0–200 nM Hfq hexamer at 30°C, measured by stopped-flow fluorescence. Black, Hfq102; red, Hfq-sCTD; cyan, Hfq65. (B) Fluorescence anisotropy assay for RNA binding and release. In stopped-flow FRET experiments, the dsRNA is released after pairing with its complementary strand (Panja et al., 2013; Santiago-Frangos et al., 2016). D16-FAM RNA (50 nM) was allowed to bind 50 nM Hfq hexamer. The Hfq·D16-FAM complex was challenged with 50 nM complementary R16 RNA. Most of the D16-FAM·R16 product is released from Hfq102 and Hfq102-sCTD, but not Hfq65. Remaining RNA was displaced from Hfq by excess ssDNA (400 nM DNA2). The averages and standard deviations for three trials are plotted for each Hfq variant. (C) Molar fractions of D16-FAM·R16 product released (light grey), remaining D16-FAM·Hfq·R16 ternary complex (hatched), and Hfq·D16-FAM binary complex (dark grey) calculated from Equation 7 (Materials and methods), based on the anisotropies at the end of the annealing phase and the maximum anisotropies of ternary complexes from equilibrium binding experiments (Table 3). (D) Annealing of 5 nM ³²P-ChiX sRNA with 30 nM *chiP* mRNA without Hfq, or with Hfq65, Hfq102 or Hfq-sCTD, as indicated above each lane. Samples were loaded on a native polyacrylamide gel 60 min or 20 s after the components were mixed. Full gel images are in Figure 4—figure supplement 2. (E) Fractions of free ChiX (black), ChiX·*chiP* dsRNA (light grey), ChiX·Hfq·*chiP* ternary complex (hatched) and ChiX·Hfq binary complex (dark grey) after 4 min of annealing, as analyzed by EMSA (Figure 4—figure supplement 2).

https://doi.org/10.7554/eLife.27049.011

Next, we used steady-state fluorescence anisotropy to examine how a shorter CTD affects the release of annealed dsRNA from the Hfq core (Figure 4B). Binding of Hfq102, Hfq65 or Hfq-sCTD to FAM-labeled D16 RNA increased the anisotropy of FAM fluorescence, as expected (Figure 4D). The smaller anisotropy of the Hfq65·D16-FAM complex is due to its smaller hydrodynamic drag, since all Hfqs have similar affinities for D16-FAM (Table 1). When complementary R16 RNA was added to the Hfq102·D16-FAM complex, the anisotropy decreased since most of the annealed dsRNA dissociated from Hfq102 (Santiago-Frangos et al., 2016) (black, Figure 4B). Whereas, when complementary RNA was added to the Hfq65·D16-FAM complex, the anisotropy increased (Santiago-Frangos et al., 2016), because a large proportion of the dsRNA remained bound to Hfq65 (cyan, Figure 4C). When the same experiment was done for Hfq-sCTD, the anisotropy decreased even further than for Hfq102, suggesting that more dsRNA was released when the CTD is shortened. This smaller anisotropy cannot be explained by the slightly smaller molecular weight of Hfq-sCTD, since the maximum anisotropies of Hfq102 and Hfq-sCTD ternary complexes during equilibrium binding experiments were very similar.

Table 1

Equilibrium dissociation constants for Hfq.

Values are the mean ± SD of three independent experiments. K_d values were determined by fluorescence anisotropy (see Methods). *Values were previously determined (Santiago-Frangos et al., 2016).

https://doi.org/10.7554/eLife.27049.014

	K_d (nM hexamer)			Hill coefficient
RNA	Hfq₁₀₂	Hfq₆₅	Hfq-sCTD	Hfq₁₀₂	Hfq₆₅	Hfq-sCTD
D16-FAM	15.5 ± 0.9*	20.0 ± 1.3*	12.9 ± 2.1	0.7 ± 0.1*	0.6 ± 0.1*	0.8 ± 0.1
D16-FAM·R16	117 ± 12	45.5 ± 2.4	181 ± 12	1.2 ± 0.1	1.0 ± 0.1	1.0 ± 0.1
minRCRB	13.9 ± 1.5*	6.46 ± 0.7*	20.1 ± 1.9	1.4 ± 0.2*	2.5 ± 0.6*	1.1 ± 0.1

As a control, we also measured the relative affinities of each Hfq for the dsRNA product (D16-FAM·R16) versus ssRNA substrate (D16-FAM), K_rel = K_d(P)/K_d(S) (Table 1). The K_rel of each protein corresponded to the efficiency of product release in the anisotropy experiment (Figure 4C and Table 3): Hfq102 had a high K_rel = 7.5, Hfq65 had a low K_rel = 2.3, and Hfq-sCTD had the highest K_rel = 14.0. Thus, the number of non-conserved residues between the acidic tip and the core of Hfq dictates the stringency of CTD autoinhibition and the efficiency of dsRNA displacement, presumably by controlling the effective concentration of the acidic tip around the Hfq core.

CTD limits sRNA-mRNA association

To determine whether the above results on RNA oligomers apply to natural RNA substrates for Hfq, we examined the annealing of the Class II sRNA ChiX to the mRNA chiP via electrophoretic mobility shift assays (EMSAs). Because the mRNA targets of class II sRNAs interact with the rim of Hfq (Schu et al., 2015; Małecka et al., 2015), we reasoned that this sRNA-mRNA pair would be sensitive to displacement by the CTD. Low nanomolar amounts of ChiX and chiP anneal very slowly at 10°C, in the absence of Hfq (Figure 4—figure supplement 2). Whereas, Hfq102, Hfq65 and Hfq-sCTD all form a ternary complex with ChiX and chiP within 20 s (Figure 4D), reaching equilibrium in a few minutes (Figure 4—figure supplement 2). These results demonstrate that the CTD is not necessarily required for annealing longer natural RNAs. In addition, Hfq65 formed the most stable ChiX·Hfq·chiP ternary complex, with almost no ChiX-Hfq65 binary complex remaining after 20 s (Figure 4E). Less ternary complex was formed by Hfq102, and the least by Hfq-sCTD. These results are consistent with the idea that the CTD limits access of chiP mRNA to the rim of Hfq.

DNA binding is regulated by the CTD

Hfq binds dsDNA and has been reported to associate with the bacterial chromosome (Kajitani et al., 1994; Takada et al., 1997; Azam and Ishihama, 1999; Updegrove et al., 2010; Jiang et al., 2015). Since the CTD peptide competes with DNA for binding to the Hfq core (Figure 3B) and inhibits binding of dsRNA and dsDNA more strongly than ssRNA (Table 1 and Figure 3), we asked whether the CTD modulates binding of Hfq to DNA. A change in DNA binding could alter the distribution of Hfq within the cell.

We quantified Hfq102 binding to linearized pUC19 DNA by measuring the change in DNA electrophoretic mobility in 1.5% agarose (Figure 5A,B). The mobility of the DNA-Hfq complexes decreased with added Hfq102, indicating increasing numbers of Hfq hexamers bound per DNA. The complexes exhibited uniform mobility at each protein concentration, however, except for a smear below the main band that may arise from dissociation of Hfq during electrophoresis (Figure 5C, top). This pattern is consistent with an equal (non-cooperative) distribution of Hfq102 between DNA molecules, although Hfq was previously suggested to bind DNA cooperatively (Cech et al., 2016). Cooperative binding to neighboring sites on the DNA would result in distinct bound and unbound populations of pUC19 (Tapias et al., 2000; Kozlov et al., 2010), which we do not observe here. The rim mutation R16A lowered the maximum gel retardation and increased the intensity of the ‘smear’, consistent with a reduced affinity of this mutant for dsDNA (Updegrove et al., 2010).

Figure 5 with 2 supplements see all

Download asset Open asset

CTD and RNAs inhibit DNA binding to Hfq by obstructing the basic rim patch.

(**A, B**) Agarose gel electrophoretic mobility shift assays of 0–3.3 μM hexamer Hfq102, Hfq102-R16A, Hfq65 or Hfq65-R16A binding to 6 nM linear pUC19 DNA (2635 nts) at 25°C, stained with SYBR Gold. Hfq65 forms large aggregates that fail to enter the gel (Figure 5—figure supplement 1). (**C, D**) Line densitometry of DNA migration in (**A, B**). Free pUC19 (no Hfq) is shown in light grey. Samples with increasing Hfq concentration are shown in darker shades of grey. (**E, F**) RNA competition. Complexes of 0.5 μM Hfq102 hexamer and 6 nM linear pUC19 (black) were challenged by 0–2 μM RNA or DNA competitor (darker shades of red). pUC19 control in the absence of Hfq is shown in light grey. Icons indicate the binding surface for each competitor. See Figure 5—figure supplement 2 for gel images.

https://doi.org/10.7554/eLife.27049.015

In contrast to the results with the full-length protein, we observed two behaviors when Hfq65 core interacted with pUC19 DNA. At low protein concentrations, the DNA was sparsely bound by Hfq65, resulting in a small mobility shift (Figure 5B). At higher protein concentrations, the DNA formed aggregates that were too large to enter the gel, resulting in a loss of signal (Figure 5B and Figure 5D, top). Aggregation of the DNA was confirmed by solubilization with Proteinase K or by pelleting assays (Figure 5—figure supplement 1). The rim mutation Hfq65-R16A rescued the formation of insoluble aggregates, confirming that dsDNA interacts with the basic rim and not the CTD itself, as previously proposed (Jiang et al., 2015). Moreover, the Hfq65-R16A complexes migrated more slowly than those formed by similar concentrations of Hfq102 (Figure 5C and Figure 5D, bottom), suggesting that Hfq65 and Hfq65-R16A bind pUC19 at higher densities than Hfq102. Thus, the CTD appears to limit Hfq102 binding to DNA, perhaps by maintaining a regular spacing between Hfq hexamers or by enforcing a dynamic equilibrium between bound and free protein. By contrast, when the Hfq core is exposed by deletion of the CTD, Hfq65 binds and aggregates dsDNA indiscriminately.

To examine which Hfq surfaces bind dsDNA, we challenged complexes of 0.5 μM Hfq102 hexamer and pUC19 DNA with 0–2 μM RNA or DNA oligomers that interact with different sites on Hfq (Figure 5E and Figure 5—figure supplement 2). A18 RNA that binds the distal face of Hfq did not compete for DNA binding (Figure 5E, top). However, when the 16 nt Target RNA that weakly interacts with the rim is appended to A18 (Target-A18), the oligomer strongly competed against DNA for binding to Hfq (Figure 5E). Target-U6 was also a good competitor, whereas the 16 nt Target sequence alone, which binds Hfq102 weakly at the rim (1 µM), was a poor competitor. DNA1 and minRCRB, which bind to the rim of Hfq (Santiago-Frangos et al., 2016; Dimastrogiovanni et al., 2014) (Figure 3C), but have low affinities for Hfq102, were weak competitors for DNA (Figure 5F). Finally, the sRNAs ChiX and RyhB were the strongest competitors, with competition saturating at a 1:1 ratio of sRNA to Hfq102 hexamer (Figure 5F).

Overall these data indicate that the CTD does not directly bind DNA as previously suggested (Updegrove et al., 2010; Jiang et al., 2015), but rather modulates the ability of the Hfq core to bind DNA, so that the extent of binding is limited and Hfq-DNA complexes remain soluble. Our data do not conflict with a low-resolution SANS model that suggests Hfq binds perpendicularly to the DNA duplex with a slight tilt (Jiang et al., 2015), but suggest that this occurs when the basic rim of Hfq interacts with the phosphate backbone of the DNA duplex. We note that the potent competition from sRNAs, which are more numerous than Hfq in the cell (Wagner, 2013), calls into question the hypothesis that E. coli Hfq regulates cellular processes via DNA binding (Sobrero and Valverde, 2012; Cech et al., 2016).

CTD–core interactions in other bacterial Hfq’s

Our results on E. coli Hfq show that the strength and frequency of CTD–core interactions depend on the number of basic residues in the core, the acidic residues in the CTD, and the linker length. Thus, our proposed mechanism for CTD-core interactions can be used to predict how the degree of CTD autoinhibition may vary among bacterial Hfq’s. We applied our de novo modeling procedure to estimate the CTD–core interactions in four other bacterial Hfqs (Figure 6A) for which the genetic function and in vitro annealing activity have been previously characterized (Zheng et al., 2016; Bohn et al., 2007; Liu et al., 2010; Rochat et al., 2015; Christiansen et al., 2004; Oglesby-Sherrouse and Vasil, 2010). We examined low energy models of each Hfq hexamer, and compared how frequently acidic CTD residues interact with basic rim and NTD residues (‘on-target’) versus other core residues (‘off-target’) (see Materials and methods) (Figure 6B,C). This comparison was quantitatively expressed as the difference in the EEC of on-target and off-target interactions (∆EEC).

Figure 6 with 3 supplements see all

Download asset Open asset

Acidic CTD and basic patch correlate with the chaperone activity of bacterial Hfq’s.

(A) Alignment of modeled Hfq sequences in order of decreasing *in vitro* RNA annealing activity. Residues are numbered according to the *E. coli* sequence. Yellow stars, residues mutated in this study; red hexagon, last residue in Hfq65; grey box, linker removed in Hfq-sCTD. (B) Average number of favorable interactions per model for each core residue with at least one acidic CTD residue (Equation 1) in the lowest energy models (≤1%). As in Figure 1B. Number of residues with < 2.0 A² accessible surface area: *P. aeruginosa,* 10; *L. monocytogenes,* 25; *B. subtilis,* 19; *S. aureus,* 12. (C) Top-down (proximal face) and side (rim) views of example low-energy models for each Hfq, as in Figure 1C. (D) Relative RNA beacon annealing rate for Target-U6 (boxes) and Target-A18 (circles) in Hfq vs. no Hfq (relative k_obs) versus the specificity of predicted CTD–core interactions (∆EEC) for *B. subtilis*, *S. aureus*, *L. monocytogenes* and *E. coli* Hfq (blue), and *P. aeruginosa* Hfq, which is more active *in vitro* than predicted by its ΔEEC (red). Annealing data are from Zheng et al. (2016).

https://doi.org/10.7554/eLife.27049.018

For E. coli Hfq, an active chaperone with a basic rim patch and long CTD, the CTD tip tended to interact with basic residues on the rim and NTD more often and more strongly than with other residues, resulting in ΔEEC = −1.11 ± 0.20 REU. This was also true for Listeria monocytogenes Hfq (ΔEEC = −1.08 ± 0.22 REU). In contrast, Bacillus subtilis (ΔEEC = 0.51 ± 0.14 REU) and Staphylococcus aureus (ΔEEC = 0.19 ± 0.05 REU) Hfq, which are inactive in our in vitro annealing assay (Zheng et al., 2016), did not exhibit specific CTD-core interactions. Finally, in models of full-length Pseudomonas aeruginosa Hfq, the CTD adopts an extended β conformation that wraps over the rim of the hexamer and places the C-terminal acidic residues near the weakly basic NTDs (ΔEEC = −0.30 ± 0.04 REU) (Figure 6C). In the absence of the NTD, however, the CTDs dock with R16 and K17 on the rim (ΔEEC = −0.50 ± 0.18 REU). Thus, Hfqs that do not anneal RNA in vitro tend to possess shorter, less acidic CTDs that form weaker and less frequent interactions with the basic rim and NTD in silico (Figure 6D). There is a similar trend between ∆EEC and the importance of Hfq for sRNA regulation in each bacterium (Zhang et al., 2013; Bohn et al., 2007; Liu et al., 2010; Rochat et al., 2015; Oglesby-Sherrouse and Vasil, 2010; Tsui et al., 1994; Nielsen et al., 2010; Rochat et al., 2012).

In the above examples, both the CTD and the core co-vary between different species. We next asked whether the CTD conferred specificity or strength to CTD-core interactions. We modeled an Hfq chimera consisting of the highly basic E. coli Sm core, fused to the shorter and slightly less acidic B. subtilis CTD. In our models, the B. subtilis CTD contacted K3 in the NTD and R17 on the rim more frequently than E. coli CTD (Figure 6—figure supplement 1A), but contacted R16 and R19, which are functionally very important (Figure 1 and Figure 2), less frequently than E. coli CTD (Figure 6—figure supplement 1A). This was corroborated with fluorescence anisotropy results showing that E. coli Hfq65 binds a BsCTD-FITC peptide about three times more weakly than its own CTD (8.7 µM vs. 2.9 µM; Figure 6—figure supplement 1B). Although we have shown that a foreign CTD can bind the core of E. coli Hfq, the ‘specificity’ of this interaction may have been lost.

Discussion

We previously found that the flexible CTD of E. coli Hfq sweeps RNAs from the proximal and rim surfaces of the Hfq ring by an unknown mechanism (Santiago-Frangos et al., 2016). Because the mechanism was not known, it was not possible to predict whether other bacterial Hfq CTDs, which are highly variable in sequence composition and length, would perform similar functions. Here, we have used computational models and experiments to show that the acidic tip of the CTD directly displaces RNA from basic patches on the rim of Hfq. The CTD’s mimicry of nucleic acids is supported by direct competition between nucleic acids and the CTD for binding the Hfq core, and stronger autoinhibition when the linker connecting the acidic tip to the core is shortened (Hfq-sCTD). The good agreement between the modeled CTD–core contacts and the contributions of individual residues to the measured CTD binding energies and to RNA annealing validates our modeling approach, and further suggests that nucleic acids and the acidic tip of the CTD interact with the same residues in the Hfq core. As expected for a nucleic acid mimic, CTD·core interactions are dominated by electrostatics and exhibit a salt dependence (Figure 2 and Figure 2—figure supplement 2) similar to that seen for the autoregulatory CTD of HTLV-1 nucleocapsid (Qualley et al., 2010).

Our results show that competition between the CTD and RNA improves the efficiency of E. coli Hfq’s chaperone activity while increasing the stringency of substrate selection. In our model, sRNA and mRNA substrates are recruited through specific interactions with the proximal or distal face of the Sm ring. When complementary RNA segments engage one or more basic patches on the rim of Hfq, these interactions favor nucleation and zippering of the double helix (Panja et al., 2013; Panja et al., 2015). Transient interactions between the CTD and the rim leads to the displacement of the dsRNA product, preventing strand dissociation and recycling Hfq. In support of the model, we observe that Hfq102 binds dsRNA less strongly than single-stranded RNA, whereas Hfq65 binds them more similarly (Table 1).

We propose that the CTD makes Hfq a more selective RNA binding protein by inhibiting access to its rim. Single-stranded nucleic acids compete with the CTD approximately in proportion to length, suggesting that nucleic acid binding to the Hfq rim has low sequence-specificity. Consequently, short RNAs that only bind the rim, such as the 16 nt Target, weakly compete with the CTD and are poor substrates for annealing. By contrast, sRNAs or mRNAs that specifically bind the proximal or distal face of Hfq strongly compete against the CTD, and gain access to the basic patch on the rim. We previously showed that the CTD increases competition among E. coli sRNAs, resulting in different levels of sRNA accumulation in the cell (Santiago-Frangos et al., 2016). It remains to be shown whether the CTD also increases the stringency of target site selection. The six CTDs, which are disordered and mobile (Beich-Frandsen et al., 2011b), exclude a large volume around the core of E. coli Hfq (Figure 1C). This excluded volume is expected to limit the number of sRNAs that may bind E. coli Hfq at any one time, perhaps further increasing the stringency of RNA selection.

The hyper-variability of the Hfq CTD among different bacteria points to a continuous optimization of autoinhibition and binding selectivity, possibly in response to the acquisition or evolution of novel sRNA-mRNA regulatory pairs (Peer and Margalit, 2014). A balance of interactions at the rim of Hfq is needed, since a CTD that inhibits RNA binding too strongly may adversely affect interactions with genuine RNA substrates (Figure 4). Conversely, our DNA binding results show that an exposed basic patch can bind DNA (and RNA) indiscriminately, hinting that a basic patch necessitates co-evolution with a ‘protective’ CTD. Intriguingly, most Hfqs contain acidic sequences at the end of the CTD (Figure 1—figure supplement 1), despite a general bias toward basic residues at protein C-termini (Berezovsky et al., 1999) and within intrinsically disordered domains (Williams et al., 2001; Lise and Jones, 2005). Additionally, the CTDs of E. coli Hfq are long enough to contact the rims of neighboring monomers (Figure 6—figure supplements 2–3), which may explain the contribution of the CTD to hexamer stability (Vincent et al., 2012), and inter-hexamer contacts (Figure 2—figure supplement 1).

Computational modeling provided atomic-scale insight to the accessible conformations of the disordered N- and C- termini of E. coli Hfq (Figure 1C). Our experimentally validated EEC metric (Figure 2E) defines residue–residue interactions more accurately with respect to experimental data than commonly used distance cutoffs (Cα–Cα and Cβ–Cβ) (Kleiger et al., 2009; Fischer et al., 2006). This de novo modeling strategy was able to identify frequent and specific CTD–rim interactions in Lm and Pa Hfq, which act in sRNA regulation and annealing (Panja et al., 2013; Zheng et al., 2016; Bohn et al., 2007; Liu et al., 2010; Rochat et al., 2015; Rochat et al., 2012), but not for Bs and Sa Hfq (Figure 6D), in agreement with in vitro experiments. This suggests that the FloppyTail algorithm could be generally useful for predicting the interactions of disordered regions with ordered domains.

Many RNA and DNA binding proteins contain disordered or flexible domains that have been implicated in cooperativity, autoinhibition and liquid phase separation (Trudeau et al., 2013; Varadi et al., 2015; Järvelin et al., 2016). Hfq is an example of an emerging paradigm of autoregulation of nucleic acid binding by nucleic acid mimic peptides. Other examples in which a disordered CTD autoinhibits RNA or DNA binding include HTLV-1 NC (Qualley et al., 2010), E. coli gyrase (Tretter and Berger, 2012), E. coli ssDNA binding protein (Kozlov et al., 2010) and mammalian high-mobility group B1 (Watson et al., 2007). Unlike HTLV-1 NC, which also remodels RNA, the Hfq CTD gives rise to dynamic cycling of bound RNAs needed to chaperone sRNA-mRNA interactions. Our modeling procedure could be utilized to screen disordered domains found in kinases, such as myosin light chain kinases and protein kinase C (Kobe and Kemp, 1999), and nucleic acid binding proteins from all kingdoms of life (Trudeau et al., 2013; Adams, 2003; Ward et al., 2004; Wang et al., 2016).

Materials and methods

Hfq alignments and sequence logos

Request a detailed protocol

All Hfq gene sequences were taken from Uniprot (UniProt Consortium, 2015). 5359 sequences were aligned using the G-INS-1 algorithm on MAFFT webservers (Yamada et al., 2016). This alignment was reduced using CD-HIT (Li and Godzik, 2006) and Max-Align (Gouveia-Oliveira et al., 2007). An unrooted, neighbor-joining tree of the remaining 985 non-redundant, representative, sequences was made on MAFFT webservers (Yamada et al., 2016). Sequence logos of re-aligned sequences from chosen clusters were generated using WebLogo (Crooks et al., 2004).

Computational modeling of the intrinsically disordered regions

Structure preparation

Request a detailed protocol

The crystal structures of E. coli (1HK9) (Sauter et al., 2003), P. aeruginosa (1U1S) (Nikulin et al., 2005), L. monocytogenes (4NL2) (Kovach et al., 2014), B. subtilis (3HSB) (Someya et al., 2012), and S. aureus (1KQ1) (Schumacher et al., 2002) Hfqs were used as starting points for the computational modeling. All crystal structures contained the hexameric form of Hfq, except for 4NL2, for which we generated the biologically relevant hexamer using the reported symmetry operations. Missing residues were appended or prepended to the crystal structures in the following manner. First, on a single subunit, absent N-terminal residues were prepended and all N-terminal residues predicted to be disordered (Buchan et al., 2013; Jones and Cozzetto, 2015) were initialized in an extended conformation, with backbone dihedral angles set to: $ϕ = - 135 °$ and $ψ = 135 °$ . Since the Hfq hexamer is C6 symmetric, the modified subunit could be symmetrized to all other subunits. The same process was repeated to append C-terminal residues, except the base of the tail (residues 64–69 in 1HK9, 1U1S, and 1KQ1, 66–71 in 4NL2, and 62–67 in 3HSB) was ‘kinked’ to point proximally as in the 1HK9 structure. For the RSNKTN tail mutant, side chains were mutated using the PyMOL ‘mutate’ function. The structures with extended termini were ‘relaxed’ with constraints, using Rosetta (Conway et al., 2014), to eliminate energetically unfavorable atomic clashes, before modeling.

Modeling

Request a detailed protocol

A modified version of the FloppyTail algorithm (Kleiger et al., 2009) was used to model the disordered termini (see Appendix 1 for step-by-step protocol). The source code is freely available to academic users through the RosettaCommons: www.rosettacommons.org. The FloppyTail algorithm generates hypothetical, low-energy conformations of disordered regions through two stages of modeling: (i) low-resolution modeling, where side chains are represented as single pseudo-atom centroids, with aggressive sampling of backbone conformational space and gradient-based minimization, and (ii) all-atom modeling, where all side-chain atoms are restored, with fine sampling of backbone conformational space, side-chain optimization, and minimization. We adapted the original algorithm to permit simultaneous modeling of multiple disordered termini and to more extensively sample conformation space. In our simulations, any Hfq region predicted to be disordered was allowed to move and approximately 500 backbone moves (changes in ϕ/ψ angles) were attempted per disordered residue (Kleiger et al., 2009) attempted ~25 backbone moves per disordered residue). Non-disordered residues had no backbone motion, but were permitted to sample side-chain conformations. In total, simulations were used to generate ~30,000 hypothetical structures for each species’ Hfq.

Analysis

PyRosetta (Chaudhury et al., 2010) was used to evaluate the energies of pairwise residue–residue interactions. Pairwise energies were computed with the talaris2014 energy function (O’Meara et al., 2015), comprised of terms capturing van der Waals, solvation, hydrogen bonding and electrostatic interactions. If a pairwise energy was unfavorable (0 or greater), we did not consider it for further analysis.

In our analysis, we considered a set of core residues $𝒞$ and a set of tail residues $𝒯$ . We calculated the average number of tail interactions for a single core residue, $x \in 𝒞$ , by counting the number of pairwise interactions, with a lower energy than the threshold, between x and every residue in the $𝒯$ and dividing by the total number of CTDs:

⟨ N_{x} ⟩ = \sum_{m o d e l s} \sum_{s u b u n i t s} \sum_{y \in 𝒯} δ (x, y) / (N_{s u b u n i t s} N_{m o d e l s}),

where

δ (x, y) = {\begin{cases} 1, & i f E (x, y) < - 1 R E U \\ 0, & i f E (x, y) \geq - 1 R E U \end{cases} a n d E (x, y) i s t h e p a i r w i s e e n e r g y .

To compute a standard deviation for the average number of interactions with core residue $x$ , we used bootstrap resampling as described in Chaudhury et al. (2011). We resampled, with replacement, our set of models for $B = 1,000$ times and re-computed $N_{x}^{'}$ (same as Equation 1, but using the resampled set of models), acquiring a standard deviation: $σ_{N}^{2} = \frac{1}{B} \sum_{B} {(N_{x}^{'} - ⟨N_{x}^{'}⟩)}^{2}$ .

In addition, we calculated the average energy for each interaction above the threshold (between one core residue, $x \in 𝒞$ , and a set of tail residues, $𝒯$ ):

⟨ E_{x : 𝒯} ⟩ = \frac{1}{N_{s u b u n i t s} N_{m o d e l s}} \sum_{m o d e l s} \sum_{s u b u n i t s} \sum_{y \in 𝒯} E (x, y) δ (x, y),

The standard deviation for the interaction energy was computed without bootstrap resampling; the energy has a distribution within a set of models, whereas the presence of an interaction is binary and only varies when the models are resampled. We compute the standard deviation as:

σ_{E}^{2} = \frac{1}{N_{s u b u n i t s} N_{m o d e l s}} \sum_{m o d e l s} \sum_{s u b u n i t s} {(\sum_{y \in 𝒯} E (x, y) δ (x, y) - ⟨ E_{x : 𝒯} ⟩)}^{2} .

Finally, EEC was computed over a set of basic core residues, by multiplying the average tail–core interaction energy by the average number of interactions per model and summing:

EEC = \sum_{ℬ} ⟨ N_{x} ⟩ ⟨ E_{x : 𝒯} ⟩ .

Standard deviation for EEC was computed by assuming that the standard deviations of the above values are independent: $σ_{E E C}^{2} = σ_{E}^{2} σ_{N}^{2} + σ_{E}^{2} {⟨ N_{x} ⟩}^{2} + σ_{N}^{2} {⟨ E_{x : 𝒯} ⟩}^{2}$ .

Tail/Core Selections (E. coli numbering)

Request a detailed protocol

Species	Core ( $𝒞$ )	Basic core ( $ℬ$ ,for EEC)	Acidic tail ( $𝒯$ )
E. coli	1–65	3, 16, 17, 19, 47	97, 99, 100, 102
P. aeruginosa	1–65	3, 5, 16, 17, 19, 47	94, 97
L. monocytogenes	1–65	2, 16, 17, 19, 35	100, 102
B. subtilis	1–65	2, 16, 17, 37	100, 102
S. aureus	1–65	10, 16, 41	65, 67, 99, 101, 102

Local concentration of the acidic CTD tip in E. coli Hfq

Request a detailed protocol

The intrinsically disordered CTD linker was represented by a worm-like chain model (Kratky and Porod, 1949) with a statistical chain segment of 35 Å or 10 residues, which is twice the persistence length of 15–20 Å for a random coil polypeptide chain (Krigbaum and Hsu, 1975; Damaschun et al., 1991; Damaschun et al., 1993; Kellermayer et al., 1997). The disordered linker region (Schumacher et al., 2002; Beich-Frandsen et al., 2011a; Dimastrogiovanni et al., 2014), was assumed to begin at residue 71 because the first few residues of the CTD tend to pack along the core (Arluison et al., 2004). The last five residues of the CTD constitute the acidic tip. The local concentration of the acidic tip on a single CTD was calculated for full-length Hfq (Hfq102) and for Hfq-sCTD:

C = \frac{(\frac{1}{V_{t a i l} - V_{c o r e}})}{N_{A}} .

In which C is the concentration of the acidic tip, V_tail is the total volume the acidic tip can access around the center of mass of a single Hfq hexamer, V_core is the inaccessible volume of the Hfq core, and N_A is Avogadro’s number. V_tail is estimated as a sphere with radius of 105 Å for Hfq102, and a radius of 70 Å for Hfq-sCTD. V_core is estimated as a cylinder with radius 31.5 Å and height 25 Å (Sauter et al., 2003).

Hfq purification

Request a detailed protocol

Untagged E. coli Hfq102, Hfq-sCTD, Hfq65, Hfq65-Q35A, Hfq65-K47A, Hfq65-R19D and Hfq65-R16A were over-expressed in E. coli BL21(DE3)Δhfq::cat-sacB cells grown in 1 L LB-Miller media (10 g/L Tryptone, 10 g/L NaCl, 5 g/L yeast extract) supplemented with 100 µg/mL ampicillin. Plasmids for over-expression of mutant Hfq proteins were created by site-directed mutagenesis of pET21b-Hfq (Zhang et al., 2002). The purification method has been previously described (Santiago-Frangos et al., 2016). In brief, resuspended cell lysates of Hfq102 and Hfq-sCTD variants were clarified by heat denaturation and untagged Hfq was purified via Ni²⁺-affinity. Lysates of Hfq65 variants were further clarified by ammonium sulfate precipitation after heat treatment, and the protein purified by hydrophobic interaction chromatography. Finally, all Hfq variants were purified by cation-exchange chromatography to remove nucleic acids (Figure 2—figure supplement 1).

Nucleic acid preparation

Request a detailed protocol

The sequences of RNA and DNA substrates are listed in Table 2. Synthetic Target RNAs, molecular beacon (Panja and Woodson, 2012a), A18, D16-FAM and R16 have been previously described (Hopkins et al., 2009). minRCRB RNA (IDT) was reduced with TCEP (tris(2-carboxyethyl)phosphine) and purified by denaturing PAGE before labeling with Cy3-maleimide (GE Healthcare), as previously described (Santiago-Frangos et al., 2016). The extent of labeling was estimated from the absorbance at 260 and 552 nm. The sRNAs ChiX, RprA, DsrA and RyhB and mRNA chiP were transcribed in vitro as previously described (Lease and Woodson, 2004). pUC19 plasmid (NEB) for Hfq binding assays was isolated from transformed DH5α cells (NEB) using Plasmid Maxi kit (QIAGEN) and digested with EcoRI (NEB) and HindIII (NEB) and purified by phenol-chloroform extraction followed by ethanol precipitation.

CTD binding and displacement

Request a detailed protocol

To measure binding of CTD-FITC, CTDpos-FITC, or BsCTD-FITC peptides (Table 2) to Hfq65 or Hfq65 mutants, the fluorescence polarization of FITC-labeled peptide was measured 3 min after the addition of 0–33.3 µM Hfq65. Anisotropy measurements were normalized to the average anisotropy in the absence of Hfq. Samples were prepared in a 100 µL cuvette containing 100 µL 50 mM Tris·HCl pH 7.5, 45 nM CTD-FITC or CTDpos-FITC, at 30°C. Fluorescence polarization with grating correction factor was measured using a Horiba Fluorolog-3 (L-format) with single excitation and emission monochromators at 495 nm and 515 nm respectively (5 nm slit widths). Titrations were performed in duplicate and the curves were fit to a single-site binding isotherm:

y = \frac{K_{a} * x}{1 + (K_{a} * x)},

in which K_a is the association constant.

Although residue K3 in the NTD was observed to bind the acidic CTD in silico, the contribution of K3 to in vitro binding could not be determined because neither Hfq65-K3S nor Hfq65-K3Q formed stable proteins.

To measure the displacement of CTD-FITC from Hfq65·CTD-FITC complexes by nucleic acids, samples were prepared in a 100 µL cuvette containing 1 µM CTD-FITC and 1.67 µM Hfq65 hexamer so that roughly 50% of CTD-FITC peptides were bound at the start of the experiment. The polarization of CTD-FITC was measured 3 min after the addition of increasing amounts of RNA or DNA as above. Competition curves were fit to:

y = \frac{m i n Y + (m a x Y - m i n Y)}{1 + {(\frac{x}{{I C}_{50}})}^{n}},

where minY and maxY are the minimum and maximum anisotropy values measured, and IC₅₀ is the concentration of nucleic acid which displaced 50% of the bound CTD-FITC from Hfq65.

RNA binding and annealing

Request a detailed protocol

Binding constants for D16-FAM or minRCRB-Cy3 (5 nM) were measured in TNK buffer (10 mM Tris·HCl, pH 7.5, 50 mM NaCl, 50 mM KCl) at 30°C by FAM fluorescence anisotropy as described before (Hopkins et al., 2009). To measure the affinity of Hfq for D16-FAM·R16 dsRNA complex, 50 nM of both RNAs were mixed and allowed to equilibrate at 30°C for 10 min before titration with Hfq. Annealing kinetics of molecular beacon (50 nM) to either Target or Target-A18 RNA (100 nM) by 0–200 nM Hfq hexamer, in 1X TNK (10 mM Tris·HCl pH 7.5, 50 mM NaCl, 50 mM KCl) buffer at 30°C was measured by stopped-flow fluorescence spectroscopy as described previously (Soper et al., 2011; Panja and Woodson, 2012b). Annealing progress curves were fit to single or double-exponential rate equations.

Anisotropy time-course

Request a detailed protocol

To measure RNA binding and release from unlabeled Hfq102, Hfq-sCTD and Hfq65 by anisotropy, the polarization of D16-FAM was recorded every 20 s for ≥ 3 min after each addition, as previously described (Santiago-Frangos et al., 2016). Samples were prepared in a 500 µL cuvette containing 50 nM D16-FAM in TNK buffer at 30°C, with additions of 50 nM Hfq102, 50 nM Hfq-sCTD, or 50 nM Hfq65 (binding phase), followed by 50 nM R16 RNA (annealing and release phase), and finally 400 nM of the ssDNA competitor, DNA2 (stimulated release phase). RNA binding and release experiments were done in triplicate for each Hfq variant. The molar fractions of released dsRNA product D16-FAM·R16 (χ_dr), remaining ternary complex D16-FAM·Hfq·R16 (χ_hdr) and binary complex Hfq·D16-FAM (χ_hd) at the end of the ‘annealing and release phase’ of the experiment were calculated from

r_{A P} = (χ_{d r} * r_{M d r}) + (χ_{h d r} * r_{M h d r}) + (χ_{h d} * r_{M h d}),

where r_AP is the anisotropy measured at the end of the ‘annealing and release’ phase in the above annealing experiments, r_Mdr is the average anisotropy of D16-FAM·R16 complex during the ‘stimulated release’ phase that is indistinguishable from its anisotropy without Hfq, r_Mhdr is the maximum anisotropy of Hfq·D16-FAM·R16 complex from equilibrium binding experiments, and r_Mhd is the maximum anisotropy of Hfq·D16-FAM complex from equilibrium binding experiments. Using the conservation of mass, χ_dr + χ_hdr + χ_hd=1, and the relative K_d values for Hfq binding to D16Fl and D16Fl·R16, χ_hd=K_rel * χ_hdr, yields an expression for the molar fraction of ternary complex:

χ_{h d r} = \frac{(r_{A P} - r_{M d r})}{(r_{M d r} * (1 + K_{d r e l})) + r_{M h d r} + (r_{h d} * K_{d r e l})}

Hfq–plasmid DNA binding assays

Request a detailed protocol

Samples (10 µL) containing linear pUC19 (0.145 nmol bp), 0–3.333 µM Hfq hexamer, in 40 mM Tris-HCl pH 7.5, 0.14 mM EDTA, 35 mM NH₄Cl, 3.7% (v/v) glycerol, 0.05% (w/v) bromophenol blue were incubated at 25°C for 30 min. 2 µL of each reaction was loaded into a sample well of a 15 × 8 cm agarose gel (1.5% (w/v) Seakem LE agarose (Lonza) in 1X TAE (40 mM Tris, 20 mM acetate, 1 mM EDTA, pH 8.0). Electrophoresis was carried out in the cold room (4°C) at 4 V/cm for 6.5 hr. Hfq was dissociated from the bound complexes by soaking agarose gels in 150 mL TBE (89 mM Tris, 89 mM borate, 2 mM EDTA, pH 8.3) and 1 M NaCl, for 30 min at 25°C, at 85 rpm. The gels were washed twice with 150 mL TBE for 10 min, stained with 1X SYBR Gold (Invitrogen) in 150 mL TBE for 45 min, and washed twice with 150 mL TBE for 10 min before imaging on a Typhoon 9410 (GE Healthcare) via excitation at 488 nm and using a 555 nm bandpass 30 emission filter. The fluorescence intensity was measured on a line from the bottom edge of the well through the middle of the lane to visualize the migration of pUC19. For RNA competition experiments, 10 µL samples were prepared as above with 0–2 µM RNA or DNA competitor and 0 or 0.5 µM Hfq102 hexamer.

Table 2

Sequences of oligomers and sRNAs.

https://doi.org/10.7554/eLife.27049.022

RNA or DNA oligomers	Sequences (5’ to 3’)
Target	GUGGUCAGUCGAGUGG
Target-U6	GUGGUCAGUCGAGUGGUUUUUU
Target-A18	GUGGUCAGUCGAGUGGAAAAAAAAAAAAAAAAAA
A18	AAAAAAAAAAAAAAAAAA
R16	GCACUUAAAAAAUUCG
Molecular beacon	FAM-GGUCCCCCACUCGACUCACCACCGGACC-DABCYL
D16-FAM	FAM-CGAAUUUUUUAAGUGC
minRCRB	Thiol-C6-CUUCCGUCCAUUUCGGACG
DNA1	TATCCGTATGACGTTCCGGACTATGCGGCTAAGGGGCAATCTTTAC
DNA2	TTTTTCAAACTGCGGATGAGACCACATATGTATATCTCCTTCTTAAAGTTAAAC
DNA2c	CAAATTGAAATTCTTCCTCTATATGTATACACCAGAGTAGGCGTCAAACTTTTT

Transcribed RNAs
RprA	gggACGGUUAUAAAUCAACAUAUUGAUUUAUAAGCAUGGAAAUC CCCUGAGUGAAACAACGAAUUGCUGUGUGUAGUCUUUGCCCAU CUCCCACGAUGGGCUUUUUUUU
DsrA	gggAACACAUCAGAUUUCCUGGUGUAACGAAUUUU UUAAGUGCUUCUUGCUUAAGCAAGUUUCAUCCCGA CCCCCUCAGGGUCGGGAUUUUUUU
RyhB	ggGCGAUCAGGAAGACCCUCGAGGAGAACCUGAAAGCA CGACAUUGCUCACAUUGCUUCCAGUAUUACUU AGCCAGCCGGGUGCUGGCUUUUU
ChiX	gggACACCGUCGCUUAAAGUGACGGCAUAAUA AUAAAAAAAUGAAAUUCCUCUUUGACGGGC CAAUAGCGAUAUUGGCCAUUUUUUU
chiP	GUAGUCAGCGAGACUUUUCUCAACGCUACUU UUUUAAUUUUUAUUUUUUCGCUGUUCACCUUUG GUGCAGCAAUUUAUACGUCAAAGAGG AUUAACCCAUGCGUACGUUUAGUGGC

Peptides
CTD-FITC	FITC-Ahx-NNAGGGTSSNYHHGSSAQNTSAQQDSEETE-COOH
CTDpos-FITC	FITC-Ahx-NNAGGGTSSNYHHGSSAQNTSAQQRSNKTN-COOH
BsCTD-FITC	FITC-Ahx-QLELE-COOH

Table 3

Parameters used to calculate molar fractions of D16-FAM complexes.

Values are the mean of at least two independent experiments determined by fluorescence anisotropy. A further description of the parameters and their usage is provided in the Methods.

https://doi.org/10.7554/eLife.27049.023

Parameter	Anisotropy	Complex/complexes
r_Mdr	0.0400	[D16-FAM·R16]
Hfq₁₀₂ r_Mhdr	0.2227	[Hfq102·D16-FAM·R16]
Hfq-sCTD r_Mhdr	0.2203	[HfqsCTD·D16-FAM·R16]
Hfq₆₅ r_Mhdr	0.2161	[Hfq65·D16-FAM·R16]
Hfq₁₀₂ r_AP	0.0732	[Hfq102·D16-FAM·R16] + [D16-FAM·R16] + [Hfq102·D16-FAM]
Hfq-sCTD r_AP	0.0663	[Hfq-sCTD·D16-FAM·R16] + [D16-FAM·R16] + [Hfq-sCTD·D16-FAM]
Hfq₆₅ r_AP	0.1099	[Hfq65·D16-FAM·R16] + [D16-FAM·R16] + [Hfq65·D16-FAM]
Hfq₁₀₂ r_Mhd	0.2195	[Hfq102·D16-FAM]
Hfq-sCTD r_Mhd	0.1985	[Hfq-sCTD·D16-FAM]
Hfq₆₅ r_Mhd	0.1935	[Hfq65·D16-FAM]
Hfq₁₀₂ r_hd	0.1332	[Hfq102·D16-FAM]
Hfq-sCTD r_hd	0.1353	[Hfq-sCTD·D16-FAM]
Hfq₆₅ r_hd	0.0969	[Hfq65·D16-FAM]

References

1. Adams JA
(2003) Activation loop phosphorylation and catalysis in protein kinases: is there functional evidence for the autoinhibitor model?
Biochemistry 42:601–607.

https://doi.org/10.1021/bi020617o
- PubMed
- Google Scholar
1. Arluison V
2. Folichon M
3. Marco S
4. Derreumaux P
5. Pellegrini O
6. Seguin J
7. Hajnsdorf E
8. Regnier P
(2004) The C-terminal domain of Escherichia coli Hfq increases the stability of the hexamer
European Journal of Biochemistry 271:1258–1265.

https://doi.org/10.1111/j.1432-1033.2004.04026.x
- PubMed
- Google Scholar
1. Attia AS
2. Sedillo JL
3. Wang W
4. Liu W
5. Brautigam CA
6. Winkler W
7. Hansen EJ
(2008) Moraxella catarrhalis expresses an unusual Hfq protein
Infection and Immunity 76:2520–2530.

https://doi.org/10.1128/IAI.01652-07
- PubMed
- Google Scholar
1. Azam TA
2. Ishihama A
(1999)
Twelve species of the nucleoid-associated protein from Escherichia coli. Sequence recognition specificity and DNA binding affinity

The Journal of Biological Chemistry 274:33105–33113.
- PubMed
- Google Scholar
1. Baba S
2. Someya T
3. Kawai G
4. Nakamura K
5. Kumasaka T
6. Expression KT
(2010) Expression, crystallization and preliminary crystallographic analysis of RNA-binding protein Hfq (YmaH) from Bacillus subtilis in complex with an RNA aptamer
Acta Crystallographica Section F Structural Biology and Crystallization Communications 66:563–566.

https://doi.org/10.1107/S1744309110009942
- PubMed
- Google Scholar
(2011a) Structural analysis of full-length Hfq from Escherichia coli
Acta Crystallographica Section F Structural Biology and Crystallization Communications 67:536–540.

https://doi.org/10.1107/S174430911100786X
- PubMed
- Google Scholar
(2011b) Structural insights into the dynamics and function of the C-terminus of the E. coli RNA chaperone Hfq
Nucleic Acids Research 39:4900–4915.

https://doi.org/10.1093/nar/gkq1346
- PubMed
- Google Scholar
(1999) Amino acid composition of protein termini are biased in different manners
Protein Engineering, Design and Selection 12:23–30.

https://doi.org/10.1093/protein/12.1.23
- PubMed
- Google Scholar
(2007) No detectable effect of RNA-binding protein Hfq absence in Staphylococcus aureus
BMC Microbiology 7:10.

https://doi.org/10.1186/1471-2180-7-10
- PubMed
- Google Scholar
(2010) Comparing models of evolution for ordered and disordered proteins
Molecular Biology and Evolution 27:609–621.

https://doi.org/10.1093/molbev/msp277
- PubMed
- Google Scholar
1. Buchan DW
2. Minneci F
3. Nugent TC
4. Bryson K
5. Jones DT
(2013) Scalable web services for the PSIPRED Protein Analysis Workbench
Nucleic Acids Research 41:W349–W357.

https://doi.org/10.1093/nar/gkt381
- PubMed
- Google Scholar
(2016) The Escherichia Coli Hfq Protein: An Unattended DNA-Transactions Regulator
Frontiers in Molecular Biosciences 3:36.

https://doi.org/10.3389/fmolb.2016.00036
- PubMed
- Google Scholar
(2010) PyRosetta: a script-based interface for implementing molecular modeling algorithms using Rosetta
Bioinformatics 26:689–691.

https://doi.org/10.1093/bioinformatics/btq007
- PubMed
- Google Scholar
1. Chaudhury S
2. Berrondo M
3. Weitzner BD
4. Muthu P
5. Bergman H
6. Gray JJ
(2011) Benchmarking and analysis of protein docking performance in Rosetta v3.2
PLoS One 6:e22477.

https://doi.org/10.1371/journal.pone.0022477
- PubMed
- Google Scholar
(2004) The RNA-binding protein Hfq of Listeria monocytogenes: role in stress tolerance and virulence
Journal of Bacteriology 186:3355–3362.

https://doi.org/10.1128/JB.186.11.3355-3362.2004
- PubMed
- Google Scholar
1. Conway P
2. Tyka MD
3. DiMaio F
4. Konerding DE
5. Baker D
(2014) Relaxation of backbone bond geometry improves protein energy landscape modeling
Protein Science 23:47–55.

https://doi.org/10.1002/pro.2389
- PubMed
- Google Scholar
(2004) WebLogo: a sequence logo generator
Genome Research 14:1188–1190.

https://doi.org/10.1101/gr.849004
- PubMed
- Google Scholar
1. Damaschun G
2. Damaschun H
3. Gast K
4. Gernat C
5. Zirwer D
(1991) Acid denatured apo-cytochrome c is a random coil: evidence from small-angle X-ray scattering and dynamic light scattering
Biochimica et Biophysica Acta (BBA) - Protein Structure and Molecular Enzymology 1078:289–295.

https://doi.org/10.1016/0167-4838(91)90571-G
- PubMed
- Google Scholar
1. Damaschun G
2. Damaschun H
3. Gast K
4. Misselwitz R
5. Müller JJ
6. Pfeil W
7. Zirwer D
(1993) Cold denaturation-induced conformational changes in phosphoglycerate kinase from yeast
Biochemistry 32:7739–7746.

https://doi.org/10.1021/bi00081a019
- PubMed
- Google Scholar
(2014) Recognition of the small regulatory RNA RydC by the bacterial Hfq protein
eLife 3:e05375.

https://doi.org/10.7554/eLife.05375
- Google Scholar
(2016) Hfq: a multifaceted RNA chaperone involved in virulence
Future Microbiology 11:137–151.

https://doi.org/10.2217/fmb.15.128
- PubMed
- Google Scholar
1. Fischer TB
2. Holmes JB
3. Miller IR
4. Parsons JR
5. Tung L
6. Hu JC
7. Tsai J
8. Jc H
(2006) Assessing methods for identifying pair-wise atomic contacts across binding interfaces
Journal of Structural Biology 153:103–112.

https://doi.org/10.1016/j.jsb.2005.11.005
- PubMed
- Google Scholar
(2015) New insight into the structure and function of Hfq C-terminus
Bioscience Reports 35:1–9.

https://doi.org/10.1042/BSR20140128
- PubMed
- Google Scholar
(2006) Small RNA regulators and the bacterial response to stress
Cold Spring Harbor Symposia on Quantitative Biology 71:1–11.

https://doi.org/10.1101/sqb.2006.71.016
- PubMed
- Google Scholar
(2007) MaxAlign: maximizing usable data in an alignment
BMC Bioinformatics 8:312.

https://doi.org/10.1186/1471-2105-8-312
- PubMed
- Google Scholar
(2009) Effect of salt and RNA structure on annealing and strand displacement by Hfq
Nucleic Acids Research 37:6205–6213.

https://doi.org/10.1093/nar/gkp646
- PubMed
- Google Scholar
(2011) Rapid binding and release of Hfq from ternary complexes during RNA annealing
Nucleic Acids Research 39:5193–5202.

https://doi.org/10.1093/nar/gkr062
- PubMed
- Google Scholar
1. Ishikawa H
2. Otaka H
3. Maki K
4. Morita T
5. Aiba H
(2012) The functional Hfq-binding module of bacterial sRNAs consists of a double or single hairpin preceded by a U-rich sequence and followed by a 3' poly(U) tail
RNA 18:1062–1074.

https://doi.org/10.1261/rna.031575.111
- PubMed
- Google Scholar
1. Jiang K
2. Zhang C
3. Guttula D
4. Liu F
5. van Kan JA
6. Lavelle C
7. Kubiak K
8. Malabirade A
9. Lapp A
10. Arluison V
11. van der Maarel JR
(2015) Effects of Hfq on the conformation and compaction of DNA
Nucleic Acids Research 43:4332–4341.

https://doi.org/10.1093/nar/gkv268
- PubMed
- Google Scholar
1. Jones DT
2. Cozzetto D
(2015) DISOPRED3: precise disordered region predictions with annotated protein-binding activity
Bioinformatics 31:857–863.

https://doi.org/10.1093/bioinformatics/btu744
- PubMed
- Google Scholar
(2016) The new (dis)order in RNA regulation
Cell Communication and Signaling 14:3.

https://doi.org/10.1186/s12964-016-0132-3
- PubMed
- Google Scholar
1. Kajitani M
2. Kato A
3. Wada A
4. Inokuchi Y
5. Ishihama A
(1994) Regulation of the Escherichia coli hfq gene encoding the host factor for phage Q beta
Journal of Bacteriology 176:531–534.

https://doi.org/10.1128/jb.176.2.531-534.1994
- PubMed
- Google Scholar
(1997) Folding-unfolding transitions in single titin molecules characterized with laser tweezers
Science 276:1112–1116.

https://doi.org/10.1126/science.276.5315.1112
- PubMed
- Google Scholar
1. Kleiger G
2. Saha A
3. Lewis S
4. Kuhlman B
5. Deshaies RJ
(2009) Rapid E2-E3 assembly and disassembly enable processive ubiquitylation of cullin-RING ubiquitin ligase substrates
Cell 139:957–968.

https://doi.org/10.1016/j.cell.2009.10.030
- PubMed
- Google Scholar
1. Kobe B
2. Kemp BE
(1999) Active site-directed protein regulation
Nature 402:373–376.

https://doi.org/10.1038/46478
- PubMed
- Google Scholar
1. Kovach AR
2. Hoff KE
3. Canty JT
4. Orans J
5. Brennan RG
(2014) Recognition of U-rich RNA by Hfq from the Gram-positive pathogen Listeria monocytogenes
RNA 20:1548–1559.

https://doi.org/10.1261/rna.044032.113
- PubMed
- Google Scholar
(2010) Regulation of single-stranded DNA binding by the C termini of Escherichia coli single-stranded DNA-binding (SSB) protein
Journal of Biological Chemistry 285:17246–17252.

https://doi.org/10.1074/jbc.M110.118273
- PubMed
- Google Scholar
1. Kratky O
2. Porod G
(1949) Diffuse small-angle scattering of X-rays in colloid systems
Journal of Colloid Science 4:35–70.

https://doi.org/10.1016/0095-8522(49)90032-X
- PubMed
- Google Scholar
1. Krigbaum WR
2. Hsu TS
(1975) Molecular conformation of bovine A1 basic protein, a coiling macromolecule in aqueous solution
Biochemistry 14:2542–2546.

https://doi.org/10.1021/bi00682a038
- PubMed
- Google Scholar
1. Lease RA
2. Woodson SA
(2004) Cycling of the Sm-like protein Hfq on the DsrA small regulatory RNA
Journal of Molecular Biology 344:1211–1223.

https://doi.org/10.1016/j.jmb.2004.10.006
- PubMed
- Google Scholar
1. Letunic I
2. Bork P
(2016) Interactive tree of life (iTOL) v3: an online tool for the display and annotation of phylogenetic and other trees
Nucleic Acids Research 44:W242–W245.

https://doi.org/10.1093/nar/gkw290
- PubMed
- Google Scholar
1. Li W
2. Godzik A
(2006) Cd-hit: a fast program for clustering and comparing large sets of protein or nucleotide sequences
Bioinformatics 22:1658–1659.

https://doi.org/10.1093/bioinformatics/btl158
- PubMed
- Google Scholar
1. Light S
2. Sagit R
3. Sachenkova O
4. Ekman D
5. Elofsson A
(2013) Protein expansion is primarily due to indels in intrinsically disordered regions
Molecular Biology and Evolution 30:2645–2653.

https://doi.org/10.1093/molbev/mst157
- PubMed
- Google Scholar
(2009) Structure of Escherichia coli Hfq bound to polyriboadenylate RNA
PNAS 106:19292–19297.

https://doi.org/10.1073/pnas.0908744106
- PubMed
- Google Scholar
1. Lise S
2. Jones DT
(2005) Sequence patterns associated with disordered regions in proteins
Proteins: Structure, Function, and Bioinformatics 58:144–150.

https://doi.org/10.1002/prot.20279
- PubMed
- Google Scholar
1. Liu J
2. Zhang Y
3. Lei X
4. Zhang Z
(2008) Natural selection of protein structural and functional properties: a single nucleotide polymorphism perspective
Genome Biology 9:r69.

https://doi.org/10.1186/gb-2008-9-4-r69
- PubMed
- Google Scholar
1. Liu Y
2. Wu N
3. Dong J
4. Gao Y
5. Zhang X
6. Mu C
7. Shao N
8. Yang G
(2010) Hfq is a global regulator that controls the pathogenicity of Staphylococcus aureus
PLoS One 5:e13069.

https://doi.org/10.1371/journal.pone.0013069
- PubMed
- Google Scholar
(2015) Structure of bacterial regulatory RNAs determines their performance in competition for the chaperone protein Hfq
Biochemistry 54:1157–1170.

https://doi.org/10.1021/bi500741d
- PubMed
- Google Scholar
1. Mikulecky PJ
2. Kaw MK
3. Brescia CC
4. Takach JC
5. Sledjeski DD
6. Feig AL
(2004) Escherichia coli Hfq has distinct interaction surfaces for DsrA, rpoS and poly(A) RNAs
Nature Structural & Molecular Biology 11:1206–1214.

https://doi.org/10.1038/nsmb858
- PubMed
- Google Scholar
(2013) Archaeal and eukaryotic homologs of Hfq: A structural and evolutionary perspective on Sm function
RNA biology 10:636–651.

https://doi.org/10.4161/rna.24538
- PubMed
- Google Scholar
(2002)
Hfq: a bacterial Sm-like protein that mediates RNA-RNA interaction

Molecular Cell 9:23–30.
- PubMed
- Google Scholar
(2010) Defining a role for Hfq in Gram-positive bacteria: evidence for Hfq-dependent antisense regulation in Listeria monocytogenes
Nucleic Acids Research 38:907–919.

https://doi.org/10.1093/nar/gkp1081
- PubMed
- Google Scholar
1. Nikulin A
2. Stolboushkina E
3. Perederina A
4. Vassilieva I
5. Blaesi U
6. Moll I
7. Kachalova G
8. Yokoyama S
9. Vassylyev D
10. Garber M
11. Nikonov S
(2005) Structure of Pseudomonas aeruginosa Hfq protein
Acta Crystallographica Section D Biological Crystallography 61:141–146.

https://doi.org/10.1107/S0907444904030008
- PubMed
- Google Scholar
1. O'Meara MJ
2. Leaver-Fay A
3. Tyka MD
4. Stein A
5. Houlihan K
6. DiMaio F
7. Bradley P
8. Kortemme T
9. Baker D
10. Snoeyink J
11. Kuhlman B
(2015) Combined covalent-electrostatic model of hydrogen bonding improves structure prediction with Rosetta
Journal of Chemical Theory and Computation 11:609–622.

https://doi.org/10.1021/ct500864r
- PubMed
- Google Scholar
1. Oglesby-Sherrouse AG
2. Vasil ML
(2010) Characterization of a heme-regulated non-coding RNA encoded by the prrF locus of Pseudomonas aeruginosa
PLoS One 5:e9930.

https://doi.org/10.1371/journal.pone.0009930
- PubMed
- Google Scholar
(2010) C-terminally truncated derivatives of Escherichia coli Hfq are proficient in riboregulation
Journal of Molecular Biology 404:173–182.

https://doi.org/10.1016/j.jmb.2010.09.038
- PubMed
- Google Scholar
1. Otaka H
2. Ishikawa H
3. Morita T
4. Aiba H
(2011) PolyU tail of rho-independent terminator of bacterial small RNAs is essential for Hfq action
PNAS 108:13059–13064.

https://doi.org/10.1073/pnas.1107050108
- PubMed
- Google Scholar
1. Panja S
2. Woodson SA
(2012a) Hfq proximity and orientation controls RNA annealing
Nucleic Acids Research 40:8690–8697.

https://doi.org/10.1093/nar/gks618
- PubMed
- Google Scholar
1. Panja S
2. Woodson SA
(2012b) Hexamer to monomer equilibrium of E. coli Hfq in solution and its impact on RNA annealing
Journal of Molecular Biology 417:406–412.

https://doi.org/10.1016/j.jmb.2012.02.009
- PubMed
- Google Scholar
(2013) Conserved arginines on the rim of Hfq catalyze base pair formation and exchange
Nucleic Acids Research 41:7536–7546.

https://doi.org/10.1093/nar/gkt521
- PubMed
- Google Scholar
(2015) Light-Triggered RNA Annealing by an RNA Chaperone
Angewandte Chemie International Edition 54:7281–7284.

https://doi.org/10.1002/anie.201501658
- PubMed
- Google Scholar
1. Peer A
2. Margalit H
(2014) Evolutionary patterns of Escherichia coli small RNAs and their regulatory interactions
RNA 20:994–1003.

https://doi.org/10.1261/rna.043133.113
- PubMed
- Google Scholar
(2010) C-terminal domain modulates the nucleic acid chaperone activity of human T-cell leukemia virus type 1 nucleocapsid protein via an electrostatic mechanism
Journal of Biological Chemistry 285:295–307.

https://doi.org/10.1074/jbc.M109.051334
- PubMed
- Google Scholar
1. Rajkowitsch L
2. Chen D
3. Stampfl S
4. Semrad K
5. Waldsich C
6. Mayer O
7. Jantsch MF
8. Konrat R
9. Bläsi U
10. Schroeder R
(2007) RNA chaperones, RNA annealers and RNA helicases
RNA Biology 4:118–130.

https://doi.org/10.4161/rna.4.3.5445
- PubMed
- Google Scholar
(2009) A conserved small RNA promotes silencing of the outer membrane protein YbfM
Molecular Microbiology 72:566–577.

https://doi.org/10.1111/j.1365-2958.2009.06688.x
- PubMed
- Google Scholar
1. Rath A
2. Glibowicka M
3. Nadeau VG
4. Chen G
5. Deber CM
(2009) Detergent binding explains anomalous SDS-PAGE migration of membrane proteins
PNAS 106:1760–1765.

https://doi.org/10.1073/pnas.0813167106
- PubMed
- Google Scholar
1. Rochat T
2. Bouloc P
3. Yang Q
4. Bossi L
5. Figueroa-Bossi N
(2012) Lack of interchangeability of Hfq-like proteins
Biochimie 94:1554–1559.

https://doi.org/10.1016/j.biochi.2012.01.016
- PubMed
- Google Scholar
1. Rochat T
2. Delumeau O
3. Figueroa-Bossi N
4. Noirot P
5. Bossi L
6. Dervyn E
7. Bouloc P
(2015) Tracking the Elusive Function of Bacillus subtilis Hfq
PLoS One 10:e0124977.

https://doi.org/10.1371/journal.pone.0124977
- PubMed
- Google Scholar
1. Sagawa S
2. Shin JE
3. Hussein R
4. Lim HN
(2015) Paradoxical suppression of small RNA activity at high Hfq concentrations due to random-order binding
Nucleic Acids Research 43:8502–8515.

https://doi.org/10.1093/nar/gkv777
- PubMed
- Google Scholar
1. Salim NN
2. Faner MA
3. Philip JA
4. Feig AL
(2012) Requirement of upstream Hfq-binding (ARN)x elements in glmS and the Hfq C-terminal region for GlmS upregulation by sRNAs GlmZ and GlmY
Nucleic Acids Research 40:8021–8032.

https://doi.org/10.1093/nar/gks392
- PubMed
- Google Scholar
(2016) C-terminal domain of the RNA chaperone Hfq drives sRNA competition and release of target RNA
PNAS 113:E6089–E6096.

https://doi.org/10.1073/pnas.1613053113
- PubMed
- Google Scholar
(2012) Small RNA binding to the lateral surface of Hfq hexamers and structural rearrangements upon mRNA target recognition
PNAS 109:9396–9401.

https://doi.org/10.1073/pnas.1202521109
- PubMed
- Google Scholar
(2003) Sm-like proteins in Eubacteria: the crystal structure of the Hfq protein from Escherichia coli
Nucleic Acids Research 31:4091–4098.

https://doi.org/10.1093/nar/gkg480
- PubMed
- Google Scholar
1. Schilling D
2. Gerischer U
(2009) The Acinetobacter baylyi Hfq gene encodes a large protein with an unusual C terminus
Journal of Bacteriology 191:5553–5562.

https://doi.org/10.1128/JB.00490-09
- PubMed
- Google Scholar
1. Schu DJ
2. Zhang A
3. Gottesman S
4. Storz G
(2015) Alternative Hfq-sRNA interaction modes dictate alternative mRNA recognition
The EMBO Journal 34:2557–2573.

https://doi.org/10.15252/embj.201591569
- PubMed
- Google Scholar
(2002) Structures of the pleiotropic translational regulator Hfq and an Hfq-RNA complex: a bacterial Sm-like protein
The EMBO Journal 21:3546–3556.

https://doi.org/10.1093/emboj/cdf322
- PubMed
- Google Scholar
1. Shi Y
2. Mowery RA
3. Ashley J
4. Hentz M
5. Ramirez AJ
6. Bilgicer B
7. Slunt-Brown H
8. Borchelt DR
9. Shaw BF
(2012) Abnormal SDS-PAGE migration of cytosolic proteins can identify domains and mechanisms that control surfactant binding
Protein Science 21:1197–1209.

https://doi.org/10.1002/pro.2107
- PubMed
- Google Scholar
1. Sobrero P
2. Valverde C
(2012) The bacterial protein Hfq: much more than a mere RNA-binding factor
Critical Reviews in Microbiology 38:276–299.

https://doi.org/10.3109/1040841X.2012.664540
- PubMed
- Google Scholar
1. Someya T
2. Baba S
3. Fujimoto M
4. Kawai G
5. Kumasaka T
6. Nakamura K
(2012) Crystal structure of Hfq from Bacillus subtilis in complex with SELEX-derived RNA aptamer: insight into RNA-binding properties of bacterial Hfq
Nucleic Acids Research 40:1856–1867.

https://doi.org/10.1093/nar/gkr892
- PubMed
- Google Scholar
(2004) Functional effects of variants of the RNA chaperone Hfq
Biochemical and Biophysical Research Communications 323:1017–1023.

https://doi.org/10.1016/j.bbrc.2004.08.190
- PubMed
- Google Scholar
(2011) Major role for mRNA binding and restructuring in sRNA recruitment by Hfq
RNA 17:1544–1550.

https://doi.org/10.1261/rna.2767211
- PubMed
- Google Scholar
(2011) Regulation by small RNAs in bacteria: expanding frontiers
Molecular Cell 43:880–891.

https://doi.org/10.1016/j.molcel.2011.08.022
- PubMed
- Google Scholar
(2002) Predicted structure and phyletic distribution of the RNA-binding protein Hfq
Nucleic Acids Research 30:3662–3671.

https://doi.org/10.1093/nar/gkf508
- PubMed
- Google Scholar
1. Takada A
2. Wachi M
3. Kaidow A
4. Takamura M
5. Nagai K
(1997) DNA binding properties of the hfq gene product of Escherichia coli
Biochemical and Biophysical Research Communications 236:576–579.

https://doi.org/10.1006/bbrc.1997.7013
- PubMed
- Google Scholar
(2000) Bacillus subtilis LrpC is a sequence-independent DNA-binding and DNA-bending protein which bridges DNA
Nucleic Acids Research 28:552–559.

https://doi.org/10.1093/nar/28.2.552
- PubMed
- Google Scholar
1. Tretter EM
2. Berger JM
(2012) Mechanisms for defining supercoiling set point of DNA gyrase orthologs: I. A nonconserved acidic C-terminal tail modulates Escherichia coli gyrase activity
The Journal of Biological Chemistry 287:18636–18644.

https://doi.org/10.1074/jbc.M112.345678
- PubMed
- Google Scholar
1. Trudeau T
2. Nassar R
3. Cumberworth A
4. Wong ET
5. Woollard G
6. Gsponer J
(2013) Structure and intrinsic disorder in protein autoinhibition
Structure 21:332–341.

https://doi.org/10.1016/j.str.2012.12.013
- PubMed
- Google Scholar
(1994) Characterization of broadly pleiotropic phenotypes caused by an hfq insertion mutation in Escherichia coli K-12
Molecular Microbiology 13:35–49.

https://doi.org/10.1111/j.1365-2958.1994.tb00400.x
- PubMed
- Google Scholar
1. UniProt Consortium
(2015)
UniProt: a hub for protein information

Nucleic Acids Research 43:204.
- Google Scholar
(2010) E. coli DNA associated with isolated Hfq interacts with Hfq's distal surface and C-terminal domain
Biochimica et Biophysica Acta (BBA) - Gene Regulatory Mechanisms 1799:588–596.

https://doi.org/10.1016/j.bbagrm.2010.06.007
- PubMed
- Google Scholar
(2016) Hfq: the flexible RNA matchmaker
Current Opinion in Microbiology 30:133–138.

https://doi.org/10.1016/j.mib.2016.02.003
- PubMed
- Google Scholar
1. Varadi M
2. Zsolyomi F
3. Guharoy M
4. Tompa P
(2015) Functional Advantages of Conserved Intrinsic Disorder in RNA-Binding Proteins
PLoS One 10:e0139731.

https://doi.org/10.1371/journal.pone.0139731
- PubMed
- Google Scholar
(2008) The C-terminal domain of Escherichia coli Hfq is required for regulation
Nucleic Acids Research 36:133–143.

https://doi.org/10.1093/nar/gkm985
- PubMed
- Google Scholar
1. Vincent HA
2. Henderson CA
3. Ragan TJ
4. Garza-Garcia A
5. Cary PD
6. Gowers DM
7. Malfois M
8. Driscoll PC
9. Sobott F
10. Callaghan AJ
(2012) Characterization of Vibrio cholerae Hfq provides novel insights into the role of the Hfq C-terminal region
Journal of Molecular Biology 420:56–69.

https://doi.org/10.1016/j.jmb.2012.03.028
- PubMed
- Google Scholar
1. Wagner EG
(2013) Cycling of RNAs on Hfq
RNA Biology 10:619–626.

https://doi.org/10.4161/rna.24044
- PubMed
- Google Scholar
1. Wang W
2. Wang L
3. Zou Y
4. Zhang J
5. Gong Q
6. Wu J
7. Shi Y
(2011) Cooperation of Escherichia coli Hfq hexamers in DsrA binding
Genes & Development 25:2106–2117.

https://doi.org/10.1101/gad.16746011
- PubMed
- Google Scholar
(2016) Disordered nucleiome: Abundance of intrinsic disorder in the DNA- and RNA-binding proteins in 1121 species from Eukaryota, Bacteria and Archaea
Proteomics 16:1486–1498.

https://doi.org/10.1002/pmic.201500177
- PubMed
- Google Scholar
1. Ward JJ
2. Sodhi JS
3. McGuffin LJ
4. Buxton BF
5. Jones DT
(2004) Prediction and functional analysis of native disorder in proteins from the three kingdoms of life
Journal of Molecular Biology 337:635–645.

https://doi.org/10.1016/j.jmb.2004.02.002
- PubMed
- Google Scholar
(2007) Mapping intramolecular interactions between domains in HMGB1 using a tail-truncation approach
Journal of Molecular Biology 374:1286–1297.

https://doi.org/10.1016/j.jmb.2007.09.075
- PubMed
- Google Scholar
1. Weichenrieder O
(2014) RNA binding by Hfq and ring-forming (L)Sm proteins: a trade-off between optimal sequence readout and RNA backbone conformation
RNA Biology 11:537–549.

https://doi.org/10.4161/rna.29144
- PubMed
- Google Scholar
1. Williams RM
2. Obradovi Z
3. Mathura V
4. Braun W
5. Garner EC
6. Young J
7. Takayama S
8. Brown CJ
9. Dunker AK
(2001)
The protein non-folding problem: amino acid determinants of intrinsic order and disorder

Pacific Symposium on Biocomputing. Pacific Symposium on Biocomputing pp. 89–100.
- PubMed
- Google Scholar
(2016) Application of the MAFFT sequence alignment program to large data-reexamination of the usefulness of chained guide trees
Bioinformatics 32:3246–3251.

https://doi.org/10.1093/bioinformatics/btw412
- PubMed
- Google Scholar
1. Zhang A
2. Wassarman KM
3. Ortega J
4. Steven AC
5. Storz G
(2002) The Sm-like Hfq protein increases OxyS RNA interaction with target mRNAs
Molecular Cell 9:11–22.

https://doi.org/10.1016/S1097-2765(01)00437-3
- PubMed
- Google Scholar
1. Zhang A
2. Schu DJ
3. Tjaden BC
4. Storz G
5. Gottesman S
(2013) Mutations in interaction surfaces differentially impact E. coli Hfq association with small RNAs and their mRNA targets
Journal of Molecular Biology 425:3678–3697.

https://doi.org/10.1016/j.jmb.2013.01.006
- PubMed
- Google Scholar
(2016) Arginine patch predicts the RNA annealing activity of Hfq from gram-negative and gram-positive bacteria
Journal of Molecular Biology 428:2259–2264.

https://doi.org/10.1016/j.jmb.2016.03.027
- PubMed
- Google Scholar
1. Zhou L
2. Hang J
3. Zhou Y
4. Wan R
5. Lu G
6. Yin P
7. Yan C
8. Shi Y
(2014) Crystal structures of the Lsm complex bound to the 3' end sequence of U6 small nuclear RNA
Nature 506:116–120.

https://doi.org/10.1038/nature12803
- PubMed
- Google Scholar

Article and author information

Author details

Andrew Santiago-Frangos

Cell, Molecular and Developmental Biology and Biophysics Program, Johns Hopkins University, Baltimore, United States

Contribution
Conceptualization, Validation, Investigation, Visualization, Writing—original draft, Writing—review and editing

Competing interests
No competing interests declared

"This ORCID iD identifies the author of this article:" 0000-0001-9615-065X
Jeliazko R Jeliazkov

Program in Molecular Biophysics, Johns Hopkins University, Baltimore, United States

Contribution
Conceptualization, Software, Validation, Investigation, Visualization, Methodology, Writing—review and editing

Competing interests
No competing interests declared

"This ORCID iD identifies the author of this article:" 0000-0003-4249-1955
Jeffrey J Gray

Department of Chemical and Biomolecular Engineering, Johns Hopkins University, Baltimore, United States

Contribution
Conceptualization, Software, Supervision, Funding acquisition, Methodology, Writing—review and editing

Competing interests
No competing interests declared

"This ORCID iD identifies the author of this article:" 0000-0001-6380-2324
Sarah A Woodson

T.C. Jenkins Department of Biophysics, Johns Hopkins University, Baltimore, United States

Contribution
Conceptualization, Supervision, Funding acquisition, Visualization, Project administration, Writing—review and editing

For correspondence
swoodson@jhu.edu

Competing interests
No competing interests declared

"This ORCID iD identifies the author of this article:" 0000-0003-0170-1987

Funding

National Institute of General Medical Sciences (R01 GM120425-01)

Sarah A Woodson

National Institute of General Medical Sciences (R01 GM078221)

Jeffrey J Gray

National Institute of General Medical Sciences (T32 GM008403-25)

Jeliazko R Jeliazkov

National Institute of General Medical Sciences (T32 GM007231-40)

Andrew Santiago-Frangos

The funders had no role in study design, data collection and interpretation, or the decision to submit the work for publication.

Acknowledgements

The authors thank Steven Lewis (Cyrus Bio), Susan Gottesman (NCI) and Agata Groszewska (Adam Mickiewicz University, Poznan) for helpful discussion, and Nadim Majdalani (NCI) for the gift of E. coli strain NM694 hfq^–.

Copyright

This article is distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use and redistribution provided that the original author and source are credited.