Water-mediated recognition of t1-adenosine anchors Argonaute2 to microRNA targets

  1. Nicole T Schirle
  2. Jessica Sheu-Gruttadauria
  3. Stanley D Chandradoss
  4. Chirlmin Joo  Is a corresponding author
  5. Ian J MacRae  Is a corresponding author
  1. The Scripps Research Institute, United States
  2. Delft University of Technology, Netherlands

Abstract

MicroRNAs (miRNAs) direct post-transcriptional regulation of human genes by guiding Argonaute proteins to complementary sites in messenger RNAs (mRNAs) targeted for repression. An enigmatic feature of many conserved mammalian miRNA target sites is that an adenosine (A) nucleotide opposite miRNA nucleotide-1 confers enhanced target repression independently of base pairing potential to the miRNA. In this study, we show that human Argonaute2 (Ago2) possesses a solvated surface pocket that specifically binds adenine nucleobases in the 1 position (t1) of target RNAs. t1A nucleotides are recognized indirectly through a hydrogen-bonding network of water molecules that preferentially interacts with the N6 amine on adenine. t1A nucleotides are not utilized during the initial binding of Ago2 to its target, but instead function by increasing the dwell time on target RNA. We also show that N6 adenosine methylation blocks t1A recognition, revealing a possible mechanism for modulation of miRNA target site potency.

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

eLife digest

Stretches of DNA known as genes provide the instructions to make the proteins and RNA molecules a cell needs to work. To make a protein, the gene is used as a template to make a type of RNA molecule called messenger RNA (mRNA), which is subsequently ‘translated’ into a protein. Most genes do not need to produce proteins all of the time, and so cells have several ways of stopping proteins from being made. For example, the Argonaute family of proteins prevents mRNA molecules from being translated into proteins.

Argonautes are guided to their targets by short RNA molecules called microRNAs. RNA molecules are made up of a sequence of building blocks known as nucleotides, each of which can only bind to one other type of nucleotide. If part of the nucleotide sequence of a microRNA molecule corresponds with part of the nucleotide sequence of the mRNA, the two RNA molecules will bind to each other. This enables the microRNA and the Argonaute protein to prevent the mRNA being translated. If the mRNA has an adenine nucleotide in a particular position (called ‘t1’) near the binding region in the mRNA sequence, Argonaute proteins will prevent translation more effectively. An adenine nucleotide in the t1 position is also known as a t1A nucleotide.

In 2014, researchers revealed the structure of a human Argonaute protein called Argonaute2 when it is bound to a microRNA-mRNA pair. This revealed that t1A nucleotides—but not other nucleotide types in the t1 position—interact with a ‘pocket’ in the Argonaute protein. However, it was not clear how the adenine nucleotide is recognized.

Now, Schirle et al.—including several of the researchers involved in the 2014 work—use a technique called X-ray crystallography to examine how the t1A nucleotide interacts with Argonaute2 in more detail. This revealed that the Argonaute2 pocket contains many water molecules that form an organized network. This network interacts with part of the t1A nucleotide and helps to lock Argonaute2 onto its microRNA target sites. The discovery of the pocket and how t1A is recognized may now be used to design more effective ‘anti-miRs’—synthetic microRNA inhibitors that can treat diseases in which microRNAs work incorrectly, a feature common to many forms of cancer.

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

Introduction

MicroRNAs (miRNAs) are an abundant class of regulatory molecules with diverse biological functions in plants and animals (Lagos-Quintana et al., 2001; Lau et al., 2001; Lee and Ambros, 2001; Reinhart et al., 2002). Over two thousand unique human miRNA sequences have been reported and more than half of all protein-coding human genes are predicted to contain a conserved miRNA recognition site (Friedman et al., 2009; Kozomara and Griffiths-Jones, 2011). miRNAs function as guides for Argonaute proteins, which use the sequence information encoded in each miRNA to identify complementary sties in mRNAs targeted for repression (Liu et al., 2004; Meister et al., 2004). Argonautes then recruit additional silencing factors that mediate translational repression and degradation of the targeted mRNAs (Huntzinger and Izaurralde, 2011).

Pioneering studies in the nematode Caenorhabditis elegans showed that perfect complementarity between miRNAs and their targets is not necessary for silencing (Lee et al., 1993; Wightman et al., 1993). Examination of regulatory elements in fly mRNAs then revealed a striking degree of complementarity to the 5′ ends of subset of conserved miRNAs, suggesting that base paring interactions with nucleotides towards the 5′ end of the miRNA are particularly important for target recognition (Lai, 2002). Indeed, phylogenetic analysis showed that pairing to the miRNA ‘seed region’ (nt 2–7 or 2–8, from the miRNA 5′ end) is the most evolutionarily conserved feature of miRNA targets in animals (Lewis et al., 2003; Brennecke et al., 2005; Krek et al., 2005; Lewis et al., 2005), and complementarity to the miRNA seed is generally sufficient to elicit significant levels of target recognition and repression (Doench and Sharp, 2004; Brennecke et al., 2005; Lim et al., 2005). Careful comparison of miRNA seed-matched sites in the 3′ untranslated regions (UTRs) of human, mouse, rat, dog, and chicken genomes also revealed that an adenosine (A) nucleotide opposite miRNA nucleotide-1 is a conserved feature of many vertebrate miRNA target sites (Lewis et al., 2005). These t1A nucleotides confer enhanced repression of miRNA targets beyond pairing to the seed region alone (Grimson et al., 2007; Nielsen et al., 2007; Baek et al., 2008; Selbach et al., 2008). Curiously, unlike the seed-matched region, t1A nucleotides function independently of the identity of miRNA nucleotide-1, indicating that t1A recognition occurs through a mechanism that is distinct from base pairing with the miRNA guide (Grimson et al., 2007; Nielsen et al., 2007; Baek et al., 2008).

We recently reported crystal structures of human Argonaute2 (Ago2) bound to a guide RNA and short, seed-paired target RNAs (Schirle et al., 2014). These structures revealed that Ago2 cradles the duplex formed by miRNA nucleotides 2–7 and the complementary target RNA, explaining why seed pairing is critical for recognition of miRNA target sites (Figure 1A). We also noted that t1A nucleotides interact with a surface pocket formed at the interface of the L2 and MID domains of Ago2, and proposed that this additional interaction may contribute to the affinity of Ago2 for miRNA target sites. However, the t1-nucleotide binding pocket is large enough to accommodate any of the four natural RNA bases, and the contacts between Ago2 and t1A appear to be almost entirely non-specific (Figure 1B). Therefore, it is not clear how A nucleotides in the t1 position of miRNA targets are recognized and confer enhanced repression. In this study, we show that the t1-nucleotide binding pocket in Ago2 preferentially interacts with A nucleotides through water-mediated contacts to adenosine N6 amine. Adding a methyl group to the t1A N6 amine reduces target affinity, raising the possibility that adenosine methylation could, in principle, lead to partial derepression of miRNA targets containing 7mer-A1 or 8mer sites. We also present data indicating that t1A is not used in the initial search for target sites, but instead provides an anchor that helps retain Ago2 on seed-matched sites on target RNAs.

Structure of the t1 nucleotide binding pocket.

(A) Linear schematic of the Argonaute2 (Ago2) primary structure. (B) Crystal structure of Ago2 bound to a guide RNA (red) and target RNA bearing t1A nucleotide (blue; PDB ID: 4W5O). (C) Close-up view of the t1-binding pocket. Ordered water molecules shown as pink spheres. Protein shown in stick and surface representations. Target RNA shown as sticks. (D) Linear schematic of crystallized guide and target RNAs.

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

Results

An adenosine-specific t1-binding pocket at the MID-L2 interface

We first visualized how Ago2 engages non-A t1 nucleotides by determining crystal structures of Ago2 bound to short target RNAs with cytosine (C), uracil (U), or guanine (G) in the t1 position (Table 1). The overall conformation of Ago2 is nearly identical in the four different t1 structures, indicating that identity of the t1 nucleotide does not affect the structure of the protein. Inspection of target RNA omit maps, in which the target RNAs were excluded from refinement, revealed little or no well-defined electron density for all of the non-A nucleotides, indicating that the non-A t1 nucleotides were mostly disordered (Figure 2). In contrast, the target omit map from the t1A data set had clearly defined density for the t1A inside the t1-nucleotide binding pocket. Moreover, although residual electron density was observed for t1G and t1C nucleotides, binding experiments using the same target RNAs show that these interactions do not contribute substantially to target affinity (Figure 2). We conclude that, although the t1-binding site is large enough to accommodate any of the four natural RNA bases, adenine is the only nucleobase that associates stably enough to contribute to overall target affinity.

Table 1

Crystallographic and refinement statistics for wild-type Ago2-guide-target complexes

https://doi.org/10.7554/eLife.07646.004
Target RNAt1-Ct1-Gt1-Ut1-DAPt1-Inosine
PDB code4Z4C4Z4D4Z4E4Z4F4Z4G
Space groupP1211P1211P1211P1211P1211
Unit cell dimensions
a (Å)55.6955.7455.6455.8655.66
b (Å)116.56117.02116.84116.60117.0
c (Å)69.6169.8769.7470.3870.1
 β (°)92.4392.4392.4392.5292.40
Ago2 per ASU11111
Data collection
 Wavelength (Å)0.979450.979500.979180.979500.97950
 Resolution (Å)38.85–2.30 (2.38–2.30)39.01–1.60 (1.63–1.60)55.60–1.80 (1.90–1.80)38.86–2.80 (2.95–2.80)39.01–2.70 (2.83–2.70)
 Total reflections133,678532,622351,63465,58775,249
 Unique reflections38,61411377582,07121,07824,061
 Completeness (%)98.4 (96.3)96.8 (93.5)99.7 (99.7)95.0 (92.6)97.2 (91.9)
 Redundancy3.5 (3.4)4.7 (4.6)4.3 (3.7)3.1 (3.0)3.1 (2.9)
 I/σI13.1 (2.2)13.7 (2.0)10.3 (1.9)9.9 (2.4)9.6 (2.2)
 Rmerge7.7 (53.0)5.5 (74.9)9.8 (81.9)9.8 (57.7)8.4 (48.2)
 Rpim7.3 (49.2)3.1 (59.0)5.3 (47.4)9.7 (57.1)6.4 (36.3)
Refinement
 Resolution (Å)35.30–2.3039.01–1.6040.27–1.8035.47–2.8035.41–2.70
 R-free/R-factor21.86/16.9318.90/16.3118.54/15.7523.30/18.3922.44/17.86
 R.M.S. deviation
  Bond distances (Å)0.0140.0140.0070.0060.008
  Bond angles (°)1.4541.4601.1310.9030.989
 Number of atoms
  Non-hydrogen, protein64296469642964216404
  Non-hydrogen, RNA568580552571572
  Phenol2828282128
  Isopropanol420800
  Phosphate100500
  Mg33333
  Water27955454516105
 Ramachandran Plot
  Most favored regions95.33%96.79%96.97%94.93%95.17%
  Additionally allowed4.54%3.21%3.03%4.94%4.83%
  Generously llowed0.13%0.00%0.00%0.13%0.00%
The t1 nucleotide binding pocket specifically recognizes adenosine Strong electron density was observed for adenosine in the t1-binding pocket, but not t1C, G, or U.

(A) T1A target RNA omit map contoured at 3σ (grey mesh). (B) T1C target RNA omit map contoured at 3σ (grey mesh). (C) T1G target RNA omit map contoured at 3σ (grey mesh). (D) T1U target RNA omit map contoured at 3σ (grey mesh). (E) Plot of target bound vs Ago2-guide concentration for target RNAs with different t1-nucleotides.

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

An unmodified purine N6 amine is required for t1 recognition

The contacts between Ago2 and t1A appear to be largely non-specific (Figure 1). Therefore, the mechanism by which Ago2 distinguishes the t1A purine ring from t1G is unclear. Two major differences between adenine and guanine are: position-6, where adenine has an exocyclic amine, and guanine has a carbonyl; and position-2, where guanine has an exocyclic amine, and adenine has no substituent group (Figure 3A). We dissected the influence of these two elements on t1-binding by examining interactions with target RNAs with either inosine (I) or 2,6-diaminopurine (DAP) in the t1 position. DAP has an N6 amine (like adenine), and an N2 amine (like guanine). Conversely, inosine has an O6 carbonyl (like guanine), but lacks an N2 amine (like adenine). Ago2 bound the t1I target with an affinity similar to t1G (Table 2) and no electron density was observed for the t1I nucleotide target omit maps, indicating that removal of the guanine N2 amine is insufficient to promote interactions with the t1-binding site (Figure 3B). In contrast, we observed unambiguous electron density for the t1DAP nucleotide, which bound Ago2 in the same position and anti conformation as t1A (Figure 3C). Moreover, the affinity of the t1DAP target is ∼1.7-fold greater than that of the equivalent t1A target, indicating that Ago2 has a modest, though significant, preference for t1DAP over t1A (p-value = 0.0025; two-tailed, unpaired Student's t-test). We conclude that the purine N6 and N2 amines both have positive effects on t1-binding, and that the N6 amine is the major determinant for distinguishing the t1A purine ring from tG. Consistent with this idea, adding a methyl group to the t1A N6 amine reduced target affinity to that of a non-t1A target (Figure 3D).

A purine N6 amine is required for t1 nucleotide recognition.

(A) Chemical structures of adenosine, 2,6-diaminopurine (DAP), guanosine, and inosine (I) nucelobases. Hydrogen bond donors highlighted in blue; hydrogen bond acceptors highlighted in pink. (B) t1I target RNA omit map contoured at 3σ (grey mesh). No electron density was observed for t1I. (C) t1DAP target RNA omit map contoured at 3σ (grey mesh). Clear electron density was observed for t1DAP. (D) Plot of target bound vs Ago2-guide concentration for t1A, t1I, t1DAP, and t1m6A target RNAs.

https://doi.org/10.7554/eLife.07646.006
Table 2

Ago2-target affinities

https://doi.org/10.7554/eLife.07646.007
t1 nucleotideKD (nM)
WT Ago2A481T-Ago2
A0.75 ± 0.041.5 ± 0.13
G1.9 ± 0.091.8 ± 0.11
U1.9 ± 0.10
C1.8 ± 0.12
DAP0.45 ± 0.03
I1.7 ± 0.09
m6A1.8 ± 0.12
  1. Dissociation constants for wild-type (WT) and mutant Ago2 binding short target RNAs with different t1 nucleotides.

An organized water network inside the t1-binding pocket

The only direct contact between the t1A N6 amine and Ago2 is a hydrogen bond to the side chain of Ser-561 (Figure 1). Because the serine alcohol can function as both a hydrogen bond acceptor and donor, it is unlikely that this interaction alone can discriminate a purine N6 amine from an O6 carbonyl. Inspection of the atoms surrounding Ser-561 revealed no set of interactions that could conceivably steer the hydrogen-bonding capacity of the serine alcohol. However, examination of the solvent surrounding the t1A nucleobase identified four ordered water molecules in a hydrogen-bonding network that could selectively interact with purine N6 amines. In this model, both hydrogen atoms in water ‘A’ are engaged in hydrogen bonds with the main chain carbonyls of Lys-440 and Met-437 (and possibly the main chain carbonyl of Asp-480), leaving only lone pair electrons available for interactions with water ‘B’ (Figure 4B). Water B must therefore donate one hydrogen atom to establish a hydrogen bond with water A. The second hydrogen atom in water B is likely involved in a hydrogen bond with the main chain carbonyl of Ile-477. Thus, our model suggests that Ago2 uses main chain carbonyls inside the t1-binding pocket to direct the hydrogen atoms on water B away from the purine 6-position, leaving only lone pair electrons available for establishing interactions with t1 nucleotides. In the proposed orientation, the lone pair electrons on water B are positioned to accept a hydrogen bond from purine N6 amines and to repel O6 carbonyls. Waters ‘C’ and ‘D’ likely provide an additional layer of selectivity through interactions with the unprotonated N1 amine on t1A, and may interact with the purine N2 amine, which does not directly contact Ago2.

Water-mediated recognition of t1A.

(A) Water network within the t1-binding pocket. Protein main chain shown as sticks, with side-chains (except S561) hidden for clarity. Water molecules shown as pink spheres. Water omit map shown contoured at 2.5σ (grey mesh). Potential hydrogen bonds shown as black dashed lines; t1A shown as blue sticks. (B) Flattened cartoon schematic of image shown in (A). Drawing illustrates proposed positions of hydrogen atoms within the t1A recognition network.

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

Disruption of the water network extinguishes t1A recognition

We tested the role of water molecules in t1A recognition by designing a mutation in Ago2 that disrupts the water network inside the t1-binding site. Ala-481, which resides in the back of the t1-binding pocket and is within van der Waals contact distance of water B, was mutated to threonine. Our rational was that the A481T mutation would alter the hydrogen-bonding landscape of the t1-binding pocket and thereby perturb placement of the associated water molecules.

We first assessed the effects of the A481T mutation on Ago2 by determining the crystal structure of the mutant protein bound to a short target RNA bearing a t1A nucleotide (Table 3). The overall structure of A481T Ago2 is nearly identical to the wild type protein, and the structure of the t1-binding pocket is almost completely unaffected by the A481T mutation. The Thr-481 side chain is in the same orientation as that of Ala-481, with the Cγ methyl pointed towards the protein interior and the Oγ alcohol extending into the t1-binding site (Figure 5A). Thus, the A481T mutation adds a hydrogen-bond acceptor/donor to the back of the t1-binding pocket without otherwise perturbing the structure of Ago2.

Table 3

Crystallographic and refinement statistics for mutant (A481T) Ago2-guide-target complexes

https://doi.org/10.7554/eLife.07646.009
Target RNAt1-At1-G
PDB code4Z4H4Z4I
Space groupP1211P1211
Unit cell dimensions
a (Å)55.6955.60
b (Å)116.60116.60
c (Å)70.1069.62
 β (°)92.2992.42
Ago2 per ASU11
Data collection
 Wavelength (Å)0.999990.99999
 Resolution (Å)44.81–2.50 (2.61–2.50)44.69–2.80 (2.95–2.80)
 Total reflections11333285,923
 Unique reflections30,41121,701
 Completeness (%)98.6 (98.3)99.2 (97.6)
 Redundancy3.7 (3.8)4.0 (4.0)
 I/σI9.2 (2.3)8.7 (2.4)
 Rmerge12.7 (69.7)13.1 (60.0)
 Rpim12.0 (65.0)12.0 (55.0)
Refinement
 Resolution (Å)44.81–2.5044.69–2.80
 R-free/R-factor21.13/17.2123.32/19.18
 R.M.S. deviation
  Bond distances (Å)0.0080.005
  Bond angles (°)0.9400.865
 Number of atoms
  Non-hydrogen, protein64326412
  Non-hydrogen, RNA568568
  Phenol2821
  Isopropanol420
  Phosphate100
  Mg33
  Water13523
 Ramachandran plot
  Most favored regions95.19%94.67%
  Additionally allowed4.68%5.33%
  Generously allowed0.13%0.00%
Disruption of the t1 pocket water network abolishes t1A recognition.

(A) Overlay of wild-type and A481T Ago2 structures and close-up view of the region surrounding the A481T mutation (inset). (B, C) Crystal structures of A481T-Ago2 bound to t1A (B) or t1G (C) target RNAs. Target RNA omit maps contoured at 3σ (grey mesh). (D) Plot of bound t1A and t1G target RNAs vs A481T-Ago2-guide concentration. (E) Overlay of A481T-Ago2 and t1A from wild-type (WT) Ago2 (semi-transparent) structures. 4.2 Å distance between the T481 Oγ hydroxyl and the t1A N6 amine indicated as dashed line.

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

Although the structure of the A481T Ago2 is nearly identical to wild type, t1A recognition is severely impaired. Only residual electron density was observed for the t1A nucleotide in the target RNA omit map (Figure 5B). For comparison, we determined the structure of the A481T mutant bound to a target bearing a t1G nucleotide. Again, the structure of the mutant protein was essentially indistinguishable from the wild type and only residual electron density was observed for the t1G (Figure 5C). Moreover, the A481T Ago2 bound t1A and t1G target RNAs with nearly identical affinities (Figure 5D). Importantly, the Oγ hydroxyl of Thr-481 is 4.2 Å from the position of the t1A N6 amine in the wild-type structure, strongly indicating that the mutation affects t1A recognition indirectly (Figure 5E). We conclude that binding specificity in the t1-pocket is achieved by Ago2-directed solvent interactions with the t1 adenosine.

t1A interactions increase the dwell time of Ago2 on miRNA target sites

Binding t1A nucleotides may increase the affinity of Ago2 for target RNAs by increasing the association rate, decreasing the rate of dissociation, or both. To distinguish between these possibilities, we used a single-molecule technique for directly observing Ago2-target binding events that we recently developed (Chandradoss et al., 2015). Briefly, Ago2 is loaded with a Cy3-labeled miRNA and introduced into a microfluidic chamber containing immobilized, Cy5-labled target RNAs. Pairing between the miRNA and the complementary site on the target RNA brings the Cy3 donor and Cy5 acceptor fluorophores into close proximity, leading to high Forster resonance energy transfer (FRET), which is observed at single-molecule resolution by total-internal-reflection microscopy (Figure 6A–C).

t1A nucleotides increase the dwell time on target sites.

(A) Cartoon schematic of a single-molecule FRET assay. (B) Sequences of miRNA and target RNA with base pairs shown. The donor fluorophore (Cy3) is positioned on the ninth nt of miRNA (counting from the 5′ end of miRNA) and the acceptor (Cy5) on target RNA opposite nt 17 of miRNA. (C) Representative time trajectory. Δτ, dwell time of interaction; kon(obs), apparent binding rate. The thin grey box indicates the time of a flow that delivers Ago2 and microRNA into the observation chamber. (D) Accumulated number of first Ago2-miRNA/target binding events vs time for RNA targets bearing t1U, t1A, or t1DAP. The number of events was normalized by the total number of target RNA strands over an imaging area. (E) Binding event dwell times fit to a double exponential decay (blue). The t1A binding events fit two populations (49.8 ± 6.2% that exhibits Δτ1 and 50.2 ± 6.2% that exhibits Δτ2; R2 = 0.998). The t1DAP binding events fit two populations (44.6 ± 4.4% that exhibits Δτ1 and 55.4 ± 4.4% that exhibits Δτ2; R2 = 0.996). Dotted grey lines represent fits to a single exponential decay. The first data columns were excluded to avoid artifacts arising from the time resolution limit. The number of events per bin was normalized by the total amount of binding events per each histogram.

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

We first determined the rates at which Ago2 associates with target RNAs containing six nucleotides of complementarity to the guide seed region (nt 2–7) and either a t1U, t1A, or t1DAP nucleotide. Measuring the time of the first arrival revealed that Ago2 recognizes the three targets at similar rates (kon(obs) = 0.040 ± 0.009, 0.033 ± 0.004, and 0.041 ± 0.008 s−1 nM−1 for t1U, t1A, and t1DAP targets, respectively), demonstrating that t1-nucleotide identity has little, if any, effect on the association of Ago2 with target RNAs (Figure 6D). In contrast, the average dwell time (Δτ) of Ago2 on the t1A and t1DAP targets was more than threefold longer than on the t1U target (Figure 6E). The combined results suggest that t1A nucleotides are not used for initial target site recognition, but instead function primarily by enhancing the stability of the Ago2:target interaction after seed-pairing.

We note that, although the Δτ distribution of t1U is well described by a single-exponential decay (Δτ = 1.64 ± 0.11 s, R2 = 0.994), the Δτ distributions of t1A and t1DAP deviate noticeably from ideal (R2 = 0.972 and 0.961 for t1A and t1DAP, respectively; Figure 6E, gray lines). This observation indicates that t1A and t1DAP binding events consist of at least two discrete populations. Assuming only two populations exist, we fit each data set to a double-exponential decay and calculated that about half of the binding events (49.8% ± 6.2% and 44.6% ± 4.4% for t1A and t1DAP, respectively) had dwell times closely matching that of the t1U target (Δτ1 (t1A) = 1.61 ± 0.29 s; Δτ1 (t1DAP) = 1.97 ± 0.24). In contrast, the other half of the binding events had an average dwell time more than fivefold longer (Δτ2 (t1A) = 8.69 ± 0.45 s; Δτ2 (t1DAP) = 12.87 ± 0.47 s). The simplest explanation for these observations is that, when using limited seed-pairing as in our experimental setup, Ago2 dissociates from the target RNA without stably engaging t1A about half of the time. However, in cases in which Ago2 does engage t1A, the affinity of the interaction is increased substantially. This effect is exaggerated with t1DAP, which stabilizes the association further still. Notably, the t1A nucleobase inserts into the t1-binding pocket at an angle nearly orthogonal to the trajectory of the paired target RNA as it extends into the central cleft of Ago2 (Figure 1). Based on this observation, we suggest that t1A binding increases the dwell time on target RNAs by anchoring Ago2 to seed-paired target sites and inhibiting dissociation.

Discussion

In this study we provide evidence that human Ago2 contains a t1-nucleotide binding pocket that specifically recognizes adenine nucleobases through water-mediated contacts. Our model is reminiscent of interactions in the trp operator/repressor complex, which also uses an ordered network of hydrogen-bonded water molecules to achieve nucleobase recognition (Otwinowski et al., 1988). We also show that the adenosine N6 amine is a key determinant of specificity and that the addition of a methyl group to the N6 amine blocks t1A recognition. Taken with the finding that m6A modifications are enriched in mammalian mRNA 3′ UTRs (Meyer et al., 2012), this result raises the intriguing possibility that, in some cases, miRNA targets may be partially derepressed via adenosine methylation at 7mer-A1 and 8mer sites (Bartel, 2009).

Insights into t1 recognition by Ago2 can also be extended to other eukaryotic members of the Argonaute protein family. The amino acids making up the t1-binding site in Ago2 are conserved in all four human Argonaute proteins, and the structure of the t1-binding pocket in Ago2 is nearly identical to the same region in the structure of human Ago1 (Faehnle et al., 2013; Nakanishi et al., 2013). In fact, two ordered water molecules corresponding to waters A and B were observed in the structure of human Ago1 crystallized in the absence of target RNA (Nakanishi et al., 2013). We therefore suggest that all four human Argonautes likely recognize t1A nucleotides using the same mechanism. Additionally, recent studies indicate that some members of the Piwi clade of Argonaute proteins also display t1-nucleotide preferences. Specifically, the Drosophila protein Aubergine, the silkmoth protein Siwi, and the mouse protein Mili preferentially cleave targets bearing a t1A; while mouse protein Miwi2 prefers targets with a t1 purine (Wang et al., 2014). We note that most Piwi proteins displaying t1 preferences have conserved residues lining the t1-binding site (with the notable exception of Mili, which appears to lack several residues homologous to the t1A-binding residues in Ago2). In contrast, Drosophila and silkmoth Ago3, which do not display t1 preferences (Wang et al., 2014), lack homology in key residues making up the t1-binding pocket (Figure 7). In C. elegans, both t1A and t1U nucleotides are conserved features of many miRNA target sites, indicating that the Argonaute proteins Alg-1 and Alg-2 also have t1 nucleotide preferences (Jan et al., 2011). Curiously, there are only two minor differences between residues making up the human Ago2 t1-binding pocket and the equivalent residues in Alg-1 and Alg-2 (Asn-429 and Arg-479 in Ago2). Conceivably, these substitutions could reorganize the t1-pocket water network to recognize both adenine and uracil nucleobases. Alternatively, we suggest that Alg-1 and Alg-2 may specifically bind t1A, and that conservation of miRNA sites with t1U is connected to recognition by one or more of the other C. elegans Argonaute proteins.

Conservation in the t1A-binding pocket.

(A) Structure of the Ago2 t1A-binding pocket with major structural elements indicated. Ago2 shown as sticks with side chains of non-highlighted residues hidden for clarity. Target RNA shown in blue. (B) Multiple sequence alignment of Ago2 with other members of the extended Argonaute protein family. Conserved structural elements colored as in panel A.

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

Finally, we suggest that understanding the t1-binding site may provide new inroads for development of novel anti-miR oligonucleotides (van Rooij and Kauppinen, 2014). We found that the addition of a t1DAP nucleotide increased target affinity fourfold compared to non-t1A targets (Table 2). More importantly, the t1-binding pocket has a chemically diverse surface and is spacious enough to accommodate nucleotide analogues much more elaborate than DAP. Moreover, displacement of ordered water molecules A and B may lead to entropic gains that translate into higher binding affinity, in a fashion similar to that used in the development of cyclic urea HIV-1 protease inhibitors (Lam et al., 1994). We suggest that t1-nucleotide analogs that make favorable contacts inside the t1-binding pocket could lead to anti-miRs that specifically target the Argonaute-miRNA complex.

Materials and methods

Oligonucleotides

All RNA oligonucleotides were synthesized by Integrated DNA Technologies (IDT, Coralville, IA), with the exception of t1I and t1DAP, which were synthesized by ValueGene, and t1m6A, which was synthesized by GE Healthcare Dharmacon (Lafayette, CO). Prior to binding studies, all target RNAs were 5′ end labeled using γ-32P-ATP and T4 polynucleotide kinase (New England Biolabs, Ipswich, MA) and subsequently purified by denaturing 16% polyacrylamide gel electrophoresis (PAGE) and ethanol precipitation.

RNA nucleotide sequences

Request a detailed protocol

RNAs used in crystal structures and equilibrium binding experiments:

Guide RNA: 5′ p-UUCACAUUGCCCAAGUCUCUU 3′;

Target RNAs (t1 nucleotides in bold):

  • t1A: 5′ CAAUGUGAAAA 3′

  • t1U: 5′ CAAUGUGAUAA 3′

  • t1C: 5′ CAAUGUGACAA 3′

  • t1G: 5′ CAAUGUGAGAA 3′

  • t1I: 5′ CAAUGUGAIAA 3

  • t1DAP: 5′ CAAUGUGA(DAP)AA 3′

  • t1m6A: 5′ CAAUGUGA(m6A)AA 3′

RNAs used in single molecule experiments

Guide RNA (U with C6 amino modifier for Cy3 attachment underlined): 5′ p-UGAGGUAUUUUUUUUUUUUUUU 3′

Target RNAs (U with C6 amino modifier for Cy5 attachment underlined): t1A: 5′ (U)20UUUUUUUUUUUACUACCUCA(U)29-biotin 3′ t1U: 5′ (U)20UUUUUUUUUUUACUACCUCU(U)29-biotin 3′

  • t1DAP: 5′ (U)20UUUUUUUUUUUACUACCUC(DAP)(U)29-biotin 3′

Protein expression and purification

Request a detailed protocol

Full length wild type and A481T Ago2 proteins were expressed in Sf9 cells using a baculovirus system (Schirle and MacRae, 2012), and purified as described previously (Schirle et al., 2014). Briefly, His6-tagged Ago2 was purified by Ni-chelate chromatography and loaded with single-stranded guide RNAs. Loaded Ago2 proteins where then isolated using the Arpón method for purifying Argonaute complexes loaded with a specified guide RNA (Flores-Jasso et al., 2013).

Crystallization and data collection

Request a detailed protocol

Ago2-guide-target complexes were formed by mixing Ago2-guide complexes with target RNAs at a 1:1.2 molar ratio and incubating at room temperature for 10 min. Crystals were grown by hanging drop vapor diffusion at 20°C. Drops contained a 1:1 ratio of Ago2-guide-target to reservoir solution (16% PEG 3350, 0.1 M Tris, pH 8.0, 0.1 M phenol, 12% isopropanol, and 10 mM MgCl2). Crystals typically appeared overnight, were harvested with nylon loops, and flash frozen in liquid N2. Data were collected under cryogenic conditions remotely at beam line 12-2 at the Stanford Synchrotron Radiation Lightsource (SSRL), and beamline 24-ID-E at the Advanced Photon Source (APS) (McPhillips et al., 2002; Soltis et al., 2008). Data were processed using XDS and Scala (Gonzalez and Tsai, 2010; Kabsch, 2010; Winn et al., 2011).

Structure refinement

Request a detailed protocol

All structures were refined using the Ago2-guide-target structure (PDB ID 4W5O), with the t1A nucleotide omitted, as a starting model. Models were built using Coot (Emsley et al., 2010) and were subjected to XYZ coordinate, TLS, and B-factor refinement using PHENIX (Adams et al., 2010). Model building and refinement continued iteratively until all interpretable electron density was modeled. Water molecules were identified automatically in Coot (2Fobs − Fcalc map, above 1.8σ, and greater than 2.4 Å and less than 3.2 Å from hydrogen bond donors or acceptors) and by manual inspection of electron density maps. All structures were refined using an Rfree set identical to that used in refinement of the original 4W5O structure. Structure figures were generated with PyMOL (Schrödinger, LLC, Portland, OR).

Equilibrium binding assays

Request a detailed protocol

Guide-loaded Ago2 samples (0–70 nM) were incubated with 0.1 nM 32P-labeled target RNAs in binding reaction buffer (30 mM Tris pH 8.0, 0.1 M potassium acetate, 2 mM magnesium acetate, 0.5 mM TCEP, 0.005% (vol/vol) NP-40, 0.01 mg/ml baker's yeast tRNA (Sigma, St. Louis, MO)), in a total reaction volume of 100 µl, for 45 min at room temperature. Protein–RNA complexes and free RNA were separated using a dot-blot apparatus (GE Healthcare Life Sciences, Pittsburgh, PA), using Protran nitrocellulose membrane (0.45 µm pore size, Whatman, GE Healthcare Life Sciences) to bind protein complexes, and Hybond Nylon membrane (Amersham, GE Healthcare Life Sciences) to capture free RNA. Samples were applied with vacuum and then washed with 100 µl of ice-cold wash buffer (30 mM Tris pH 8.0, 0.1 M potassium acetate, 2 mM magnesium acetate, 0.5 mM TCEP). Membranes were air-dried and 32P signal visualized by phosphorimaging. Quantification was performed using ImageQuant software (GE Healthcare Life Sciences), and dissociation constants calculated using Prism version 5.0e (GraphPad Software, Inc., La Jolla, CA).

Single molecule FRET

Request a detailed protocol

Ago2 single molecule binding measurements were performed as described elsewhere (Chandradoss et al., 2015). Briefly, target RNAs bearing a Cy5 dye (GE Healthcare) and a 3′ biotin were immobilized on a polymer(PEG)-coated quartz surface in the microfluidic chamber (Chandradoss et al., 2014) of a prism-type total internal reflection fluorescence microscope. Ago2 was loaded with a guide miRNA containing a Cy3 dye (GE Healthcare). The resulting complex was introduced into the microfluidic chamber and Cy3 molecules were excited with a 532 nm diode laser (Compass 215M/50 mW, Coherent). Fluorescence signals of Cy3 and Cy5 were collected through a 60× water immersion objective (UplanSApo, Olympus, Center Valley, PA) with an inverted microscope (IX73, Olympus). Laser scattering was blocked by a 532 nm long pass filter (LPD01-532RU-25, Semrock, Rochester, NY). The Cy3 and Cy5 signals were separated with a dichroic mirror (635 dcxr, Chroma, Bellows Falls, VT) and imaged using a EM-CCD camera (iXon Ultra, DU-897U-CS0-#BV, Andor Technology, United Kingdom) and described previously (Selvin and Ha, 2007).

Data acquisition and analysis

Request a detailed protocol

CCD images of time resolution 0.1 s were recorded, and time traces were extracted from the CCD image series using IDL (ITT Visual Information Solution, Boulder, CO). Colocalization between Cy3 and Cy5 signals was carried out with a custom-made mapping algorithm written in IDL. The extracted time traces were processed using Matlab (MathWorks, Natick, MA) and Origin (Origin Lab, Northampton, MA). The binding rate (kon) was determined by first measuring the time between when Ago2-miRNA was introduced to a microfluidic chamber and when the first Ago2-miRNA was docked to a target RNA; and then fitting the time distribution with a single-exponential growth curve, A(1ekont). The dissociation rate was estimated by measuring the dwell time of a binding event. A dwell time distribution was fitted by either a single-exponential decay curve (Aetτ) or a double-exponential decay curve (A1et/Δτ1+A2et/Δτ2). In case of a double-exponential decay, the percentages of Δτ1 and Δτ2 populations are determined by A1Δτ1/(A1Δτ1 + A2Δτ2) and A2Δτ2/(A1Δτ1 + A2Δτ2), and the average dwell time is determined by (A1Δτ12+A2τ22)/(A1Δτ1+A2Δτ2).

Single molecule sample preparation

Request a detailed protocol

A microfluidic chamber was incubated with 20 µl Streptavidin (0.1 mg/ml, Sigma) for 30 s. Unbound Streptavidin was washed with 100 µl of buffer T50 (10 mM Tris–HCl [pH 8.0], 50 mM NaCl buffer). The 50 μl of 50 pM acceptor-labelled mRNA construct were introduced into the chamber and incubated for 1 min. Unbound labeled constructs were washed with 100 µl of buffer T50. The effector complex was formed by incubating 50 nM purified recombinant hAgo2 with 0.5 nM of donor-labeled hsa-let-7a miRNA in a buffer containing 50 mM Tris–HCl [pH 8.0] (Ambion, Grand Island, NY), 50 mM NaCl (Ambion), and 60 mM KCl (Ambion) at 31°C for 20 min. An imaging buffer for single-molecule FRET was added before the mixture was injected to a microfluidic chamber. The final concentration of the imaging buffer consists of the 0.8% dextrose (Sigma), 0.5 mg/ml glucose oxidase (Sigma), 85 µg/ml Catalase (Merck), and 1 mM Trolox ((±)-6-hydroxy-2,5,7,8-tetramethylchromane-2-carboxylic acid, 238813, Sigma). The experiments were performed at the room temperature (23 ± 2°C).

References

    1. Emsley P
    2. Lohkamp B
    3. Scott WG
    4. Cowtan K
    (2010) Features and development of coot
    Acta Crystallographica. Section D, Biological Crystallography 66:486–501.
    https://doi.org/10.1107/S0907444910007493
    1. Gonzalez A
    2. Tsai Y
    (2010) Autoxds
    In http://smb.slac.stanford.edu/facilities/software/xds/-autoxds_script.
    1. Kabsch W
    (2010) Xds
    Acta Crystallographica. Section D, Biological Crystallography 66:125–132.
    https://doi.org/10.1107/S0907444909047337
  1. Book
    1. Selvin PR
    2. Ha T
    (2007)
    Single-molecule techniques: a laboratory manual (1st edition)
    Cold Spring Harbor, NY: Cold Spring Harbor Laboratory Press.

Decision letter

  1. Phillip D Zamore
    Reviewing Editor; Howard Hughes Medical Institute, University of Massachusetts Medical School, United States

eLife posts the editorial decision letter and author response on a selection of the published articles (subject to the approval of the authors). An edited version of the letter sent to the authors after peer review is shown, indicating the substantive concerns or comments; minor concerns are not usually shown. Reviewers have the opportunity to discuss the decision before the letter is sent (see review process). Similarly, the author response typically shows only responses to the major concerns raised by the reviewers.

Thank you for submitting your work entitled “Water-mediated recognition of t1-adenosine anchors Argonaute2 to microRNA targets” for peer review at eLife. Your submission has been favorably evaluated by John Kuriyan (Senior Editor), and three reviewers, one of whom is a member of our Board of Reviewing Editors.

The following individuals responsible for the peer review of your submission have agreed to reveal their identity: Phillip Zamore (Reviewing Editor), Yukihide Tomari, and David Bartel (peer reviewers).

The reviewers have discussed the reviews with one another and the Reviewing editor has drafted this decision to help you prepare a revised submission.

Despite the presence of a complementary guide g1 nucleotide, the t1 target nucleotide is typically recognized by Argonaute proteins rather than base pairing. In mammals, miRNA target sites often have an adenosine (A) across from the first nucleotide of the guide miRNAs. Such a “t1A” preference is observed also for a subset of PIWI proteins in flies and mice. Recent crystal structures of human Ago2 revealed the existence of a surface pocket at the interface of the L2 linker and the MID domain, which can accommodate the t1 nucleotide. However, the pocket is apparently too large to discriminate among the four nucleotides, so it has been unclear why t1A is specifically preferred. In this manuscript, Schirle et al. report the structures of human Ago2 bound to target RNAs with various nucleotides at the t1 position. MacRae and colleagues find report that a network of water molecules in the t1A pocket provides the specificity for recognizing adenine, especially via its N6 amine. Using single-molecule imaging, the authors suggest that the t1A nucleotide is not used for the initial association with target, but instead slows dissociation of RISC from the target. Finally, the authors identify a base analog, di-aminopurine, that binds more tightly than adenine, leading them to suggest that novel bases could be designed to provide additional specificity to anti-miRNA drugs. The manuscript would be improved by addressing the points below, but otherwise is entirely appropriate for publication in eLife.

1) “[S]liding off of target sites” implies a specific physical model in which Ago2 departs from its seed-matched binding site but not the RNA itself. Does the evidence really support this?

2) In the determination of the koff for t1A, a test statistic or some other analysis (e.g., residuals) should be provided showing the statistical confidence that the data are better explained by a double exponential than a single exponential. If a double-exponential is, in fact, justified, what is the meaning of the t1A single-molecule data fitting a double exponential? That is, what is the physical basis for the 16% of RISC that is slow-to-dissociate? Might it be a non-biological interaction caused by the slide or its coating?

3) Given that the water molecules are critical in the wild-type t1A structure, the omit map for the four water molecules should be shown (for Figure 1C or 4A). The authors should also note whether the four water molecules are observed in the non-t1A structures. In other words, is the hydrogen-bonding network of water molecules formed independently of or in coordination with t1A binding?

4) At the end of the subsection “A purine N6 amine is required for t1 recognition”, the authors state that the N2 amine has a positive effect on t1 binding but do not describe what kind of interactions are made between the N2 amine and Ago2. Do the authors observe the similar interactions for the N2 amine of t1G? Also, it would be timely for the authors to comment on how 6-methyl adenine fits into the t1A pocket compared to adenine.

5) At the end of the subsection “Disruption of the water network extinguishes t1-A recognition“ and in Figure 5A, it is difficult to see the positional relationship between T481 and t1A, which probably could be better shown in a separate figure.

6) It would be good to test by single-molecule analysis that the A481T mutant and/or t1-Inosine do not increase the dwell time.

7) The kon and koff values from the single-molecule data in Figure 6 result in KD values 50-fold higher than those determined from the ensemble assays (Table 2). To prevent confusion, the authors should acknowledge this difference and state whether they think it is due to the different sequences used in the two types of experiments or to the Cy3 on position 9 of the miRNA used for the single-molecule experiments.

8) It is not apparent how the behavior depicted in Figure 6D is used to calculate the kon for RISC-target association. The concentration of RISC in the microfluidic chamber should be provided in the manuscript, as well as a more detailed description of the method by which kon is calculated.

9) The authors should comment on the consistency between the ensemble (Table 2) and the single-molecule data with respect to the double-exponential interpretation of t1A dissociation kinetics. Specifically, the authors should calculate the effects on the ensemble koff for t1A predicted from the single-molecule experiments and verify that the ratio of koff values with t1A and t1U targets is agrees with the ratio of KD values for t1A and t1U targets. (These numbers appear to work out by one Reviewer's calculations, but it would be useful to show that the two methods are in agreement.)

10) Can the authors explain why <40% of target RNAs are bound by Ago2-RISC in Figure 6D? What is the specific activity of Ago2-RISC used in this study?

11) The photobleaching time should be shown for Cy3 alone, Cy5 alone and FRET.

12) A multiple alignment of Ago family proteins would help the discussion of t1A-binding pocket conservation.

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

Author response

Despite the presence of a complementary guide g1 nucleotide, the t1 target nucleotide is typically recognized by Argonaute proteins rather than base pairing. In mammals, miRNA target sites often have an adenosine (A) across from the first nucleotide of the guide miRNAs. Such a “t1A” preference is observed also for a subset of PIWI proteins in flies and mice. Recent crystal structures of human Ago2 revealed the existence of a surface pocket at the interface of the L2 linker and the MID domain, which can accommodate the t1 nucleotide. However, the pocket is apparently too large to discriminate among the four nucleotides, so it has been unclear why t1A is specifically preferred. In this manuscript, Schirle et al. report the structures of human Ago2 bound to target RNAs with various nucleotides at the t1 position. MacRae and colleagues find report that a network of water molecules in the t1A pocket provides the specificity for recognizing adenine, especially via its N6 amine. Using single-molecule imaging, the authors suggest that the t1A nucleotide is not used for the initial association with target, but instead slows dissociation of RISC from the target. Finally, the authors identify a base analog, di-aminopurine, that binds more tightly than adenine, leading them to suggest that novel bases could be designed to provide additional specificity to anti-miRNA drugs. The manuscript would be improved by addressing the points below, but otherwise is entirely appropriate for publication in eLife.

We are grateful to the reviewers for their interest and careful evaluation of our work. We have revised the manuscript to address all of the criticisms and suggestions. The revised paper includes new data showing that adenosine N6 methylation abrogates t1A recognition by Ago2, and that the N2 purine amine extends average dwell time of the Ago2-guide complex on target RNA sites. These findings support and extend conclusions stated in the original submission.

1) “[S]liding off of target sites” implies a specific physical model in which Ago2 departs from its seed-matched binding site but not the RNA itself. Does the evidence really support this?

This point is well taken. We changed “sliding off of target sites” to “dissociating from target sites”.

2) In the determination of the koff for t1A, a test statistic or some other analysis (e.g., residuals) should be provided showing the statistical confidence that the data are better explained by a double exponential than a single exponential.

The calculated coefficients of determination for the koff data (now indicated in the main text) are as follows, (Note: we have included new data for t1DAP binding): t1U: R2 = 0.994 (single exponential) t1A: R2 = 0.972 (single); 0.998 (double) t1DAP: R2 = 0.961 (single); 0.996 (double)

Although all data were collected in the same way, the distribution of t1A and t1DAP dwell times do not fit a single exponential as well as the t1U data do. To better illustrate this point, we displayed dwell time distributions as semi-log plots, in which t1A and t1DAP both clearly deviate from the linear relationship expected of a single exponential function (see updated Figure 6).

If a double-exponential is, in fact, justified, what is the meaning of the t1A single-molecule data fitting a double exponential? That is, what is the physical basis for the 16% of RISC that is slow-to-dissociate? Might it be a non-biological interaction caused by the slide or its coating?

A surface artifact is unlikely because the slow-to-dissociate population was not observed with t1U, which was tested on the same slide (but in a different chamber) as t1A. The simplest explanation is that the two populations represent binding events in which Ago2 either does or does not stably engage t1A before dissociating. This model is supported by: 1) the observation that the dwell time of the short-lived t1A binding events (Δτ1 = 1.61 sec) closely matches the dwell time of t1U (Δτ = 1.64 sec); and, 2) the long-lived dwell time (Δτ2 = 8.69 sec) can be increased further still by replacing t1A with t1DAP (Δτ2 = 12.87 sec); and, 3) t1A and t1DAP do not influence on-rates, demonstrating that their effects on affinity manifest after seed-pairing.

Separately, we note that in the first submission we used incorrect formulas for estimating the percentages of Δτ1 and Δτ2 populations within a double exponential decay. Specifically, in the initial submission, we used wrong equations for the percentages, A1/(A1 + A2) and A2/(A1 + A2), where A1 and A2 are the coefficients of exponentials. This estimate resulted in 16% for the slow-to-dissociate population. With correct equations, A1Δτ1/A1Δτ1 + A2ΔτA and A2τ2/A1Δτ1 + A2Δτ2, we calculate that the slow-to-dissociate population constitutes ∼50% of the binding events. This change does not alter our interpretation of the slow-to-dissociate population. Rather, it indicates that the t1A influences the dissociation pathway to a larger extent than we initially estimated. The correct formulas are described in a new section in the Materials and methods.

3) Given that the water molecules are critical in the wild-type t1A structure, the omit map for the four water molecules should be shown (for Figure 1C or 4A). The authors should also note whether the four water molecules are observed in the non-t1A structures. In other words, is the hydrogen-bonding network of water molecules formed independently of or in coordination with t1A binding?

As suggested, we added a water-omit map calculated from the t1-A dataset to Figure 4A. We also note that ordered water molecules corresponding to waters A, B and D are clearly observed in the t1-U electron density map. Taken with the observation of water molecules equivalent to A and B in the crystal structure of Ago1 in the absence of target RNA (Nakanishi, et al., Cell Reports 2013; PDB ID 4KXT), it appears that Ago organizes t1-binding pocket water molecules, at least to some extent, prior to t1-A binding.

4) At the end of the subsection “A purine N6 amine is required for t1 recognition”, the authors state that the N2 amine has a positive effect on t1 binding but do not describe what kind of interactions are made between the N2 amine and Ago2. Do the authors observe the similar interactions for the N2 amine of t1G?

The N2 amine does not directly contact Ago2 and thus it is difficult to explain the observed increases in dwell time and affinity. Our currently favored explanation is that the N2 amine may make positive interactions with water C, possibly moving it into hydrogen-bonding distance of water B. This idea is based on a positive difference density peak, which may correspond to the repositioned water C, observed in the t1-G data (Author response image 1). We speculate that a similar change in solvent position may occur with t1-DAP binding. Unfortunately, the t1-DAP dataset is relatively low resolution (2.8 Å) and we were not able to model any water molecules in the t1-pocket with confidence. Without direct observation of repositioned water molecules in the t1-DAP structure, we are hesitant to put this model forward and thus simply added the sentence: “Waters “C” and “D” likely provide an additional layer of selectivity through interactions with the unprotonated N1 amine on t1-A, and may interact with the purine N2 amine, which does not directly contact Ago2” to state that N2 does not contact Ago2”.

Author response image 1
Electron density maps surrounding the t1-G nucleotide.

2Fo-Fc map is shown contoured at 1σ(blue mesh). Fo-Fc map (difference map) is shown contoured at 3σ (green mesh) and -3σ (red mesh). A positive difference density peak, possibly representing a repositioned water C, is observed 2.5 Å from water B and 2.7 Å from the t1-G N2 amine. Position of water C in t1-A structure is included for comparison.

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

Also, it would be timely for the authors to comment on how 6-methyl adenine fits into the t1A pocket compared to adenine.

We thank the reviewers for this thoughtful suggestion and decided the issue important enough to warrant determining experimentally how adenosine N6 methylation affects target binding. We found that a t1-m6A target RNA had an affinity similar to that of equivalent non-t1-A targets (see updated Figure 3D and Table 2), indicating that the modification effectively blocks t1-A recognition. This result further supports the idea that the N6 amine is a key determinant in t1-A recognition and raises the intriguing possibility that adenosine methylation could lead to partial derepression of miRNA targets containing 7mer-A1 or 8mer sites.

5) At the end of the subsection “Disruption of the water network extinguishes t1-A recognition“ and in Figure 5A, it is difficult to see the positional relationship between T481 and t1A, which probably could be better shown in a separate figure.

We apologize for this shortcoming and added panel E to Figure 5 to specifically illustrate the spatial relationship between T481 and t1A.

6) It would be good to test by single-molecule analysis that the A481T mutant and/or t1-Inosine do not increase the dwell time.

We thank the reviewers for this thoughtful suggestion. We conducted a slightly different experiment and measured association/dissociation rates of a target RNA with t1-DAP, with dual purposes of further aligning the bulk and single molecule studies, as well as providing insights into why t1-DAP binds better than t1-A. Consistent with the bulk studies, single molecule measurements showed an increase in affinity for t1-DAP over t1-A. The increase in affinity comes entirely from an increase in the average dwell time, with the on-rate unaffected by the additional N2 amine. These results support the notion that t1-nucleotides are not used in the initial target search.

7) The kon and koff values from the single-molecule data in Figure 6 result in KD values 50-fold higher than those determined from the ensemble assays (Table 2). To prevent confusion, the authors should acknowledge this difference and state whether they think it is due to the different sequences used in the two types of experiments or to the Cy3 on position 9 of the miRNA used for the single-molecule experiments.

We thank the reviewers for identifying this potential point of confusion. The target RNAs used in the ensemble and the single molecule experiments are different from each other in several respects. The ensemble assays used very short target RNAs (11 nt) that pair to guide nt 2–9, which provided the best diffracting crystals. Our rational was that using the same RNAs for binding experiments and structure determination would be the best way to correlate structural findings with affinity measurements. The single molecule experiments, on the other hand, required longer target RNAs (69 nt) for reliable photophysics (to establish an optimum distance between donor and acceptor dyes) that paired to guide nt 2–7 (to work in an affinity range where dwell times could be accurately measured). We added a footnote to the subsection “t1A interactions increase the dwell time of Ago2 on miRNA target sites” to direct the reader to these differences.

We also note that, despite substantially different experimental setups and differences in target RNA sequence, length and guide-complementarity, all measured affinities fall within the low nM range. More importantly (for the focus of this study), the relative effects of t1-nucleotide identity are consistent between diverse experiments, and thus appear to be largely independent of these factors.

8) It is not apparent how the behavior depicted in Figure 6D is used to calculate the kon for RISC-target association. The concentration of RISC in the microfluidic chamber should be provided in the manuscript, as well as a more detailed description of the method by which kon is calculated.

We are sorry for this oversight. We added new sections to the Materials and methods that are dedicated to explaining the experimental conditions (Single-molecule sample preparation) and the data analysis (Data acquisition and analysis).

9) The authors should comment on the consistency between the ensemble (Table 2) and the single-molecule data with respect to the double-exponential interpretation of t1A dissociation kinetics. Specifically, the authors should calculate the effects on the ensemble koff for t1A predicted from the single-molecule experiments and verify that the ratio of koff values with t1A and t1U targets is agrees with the ratio of KD values for t1A and t1U targets. (These numbers appear to work out by one Reviewer's calculations, but it would be useful to show that the two methods are in agreement.)

The average dwell times of t1U and t1A are 1.64 and 5.17 sec, respectively (for calculations see equations in Data acquisition and analysis). The ratio between <Δτt1U> and <Δτt1A> is therefore 3.2, which is similar to the ratio between KDt1U and KDt1A, (2.5). Thus, the effects of t1A are relatively consistent between ensemble and single-molecule data.

10) Can the authors explain why <40% of target RNAs are bound by Ago2-RISC in Figure 6D? What is the specific activity of Ago2-RISC used in this study?

We thank the reviewers for this careful observation and have repeated all the single molecule measurements using a new preparation of protein and higher sample concentrations (originally 10 nM Ago2 mixed with 1 nM RNA guide, and now 50 nM Ago2 mixed with 0.5 nM guide RNA) to assemble RISC. Under the new conditions we could observe that the percentage of occupied target RNA reaches nearly 100%. We note that the binding rates we observe also increased, which might be due to the increase of a functional RISC concentration in the new sample preparation.

11) The photobleaching time should be shown for Cy3 alone, Cy5 alone and FRET.

We measured photobleaching of the donor Cy3 dye, which has a half time (394 seconds) > 20x longer than the average dwell times of the observed binding events (Author response image 2). Cy5 and FRET photobleaching time is not necessary to measure since it is only Cy3 signals that are used for measuring kon and koff.

Author response image 2
Photobleaching of the Cy3 donor.
https://doi.org/10.7554/eLife.07646.016

12) A multiple alignment of Ago family proteins would help the discussion of t1A-binding pocket conservation.

We made a sequence alignment of Ago2 with Alg1 and Piwi proteins and related it to major structural elements of the t1-binding pocket in a new Figure 7.

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

Article and author information

Author details

  1. Nicole T Schirle

    Department of Integrative Structural and Computational Biology, The Scripps Research Institute, La Jolla, United States
    Contribution
    NTS, Conception and design, Acquisition of data, Analysis and interpretation of data, Drafting or revising the article
    Competing interests
    The authors declare that no competing interests exist.
  2. Jessica Sheu-Gruttadauria

    Department of Integrative Structural and Computational Biology, The Scripps Research Institute, La Jolla, United States
    Contribution
    JS-G, Conception and design, Acquisition of data, Analysis and interpretation of data
    Competing interests
    The authors declare that no competing interests exist.
  3. Stanley D Chandradoss

    Kavli Institute of NanoScience, Department of BioNanoScience, Delft University of Technology, Delft, Netherlands
    Contribution
    SDC, Conception and design, Acquisition of data, Analysis and interpretation of data
    Competing interests
    The authors declare that no competing interests exist.
  4. Chirlmin Joo

    Kavli Institute of NanoScience, Department of BioNanoScience, Delft University of Technology, Delft, Netherlands
    Contribution
    CJ, Conception and design, Analysis and interpretation of data, Drafting or revising the article
    For correspondence
    C.Joo@tudelft.nl
    Competing interests
    The authors declare that no competing interests exist.
  5. Ian J MacRae

    Department of Integrative Structural and Computational Biology, The Scripps Research Institute, La Jolla, United States
    Contribution
    IJMR, Conception and design, Analysis and interpretation of data, Drafting or revising the article
    For correspondence
    macrae@scripps.edu
    Competing interests
    The authors declare that no competing interests exist.

Funding

National Institute of General Medical Sciences (NIGMS) (R01-GM104475, P41-GM103403)

  • Ian J MacRae

European Research Council (ERC) (309509)

  • Chirlmin Joo

U.S. Department of Energy (DOE) (DE-AC02-76SF00515, DE-AC02-05CH11231, DE-AC02-06CH11357)

  • Ian J MacRae

American Heart Association (AHA) (Graduate Student Fellowship)

  • Nicole T Schirle
  • Jessica Sheu-Gruttadauria

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

Acknowledgements

We are grateful to JR Williamson, DP Bartel, JM Claycomb, E Westhof and JS Oakdale for generous advice and insights. Diffraction data were collected at beam lines: 11-1 and 12-2 at the Stanford Synchrotron Radiation Lightsource, supported by the U.S. Department of Energy, Office of Science, Office of Basic Energy Sciences under Contract No. DE-AC02-76SF00515; 5.0.2 at the Advanced Light Source, supported by the Director, Office of Science, Office of Basic Energy Sciences, of the U.S. Department of Energy under Contract No. DE-AC02-05CH11231; and, 24-ID-E at the Advanced Photon Source, supported by NIGMS grant P41 GM103403 and DOE contract number DE-AC02-06CH11357. NTS and JSG are pre-doctoral fellows of the American Heart Association. NTS is an Achievement Rewards for College Scientists (ARCS) scholar. CJ was funded by the Open Program of the Division for Earth and Life Sciences (822.02.008) of the Netherlands Organization for Scientific Research and European Research Council under the European Union's Seventh Framework Programme [FP7/2007–2013]/ERC grant agreement no. [309509]. The work was supported by NIH grant R01 GM104475 to IJM. Atomic coordinates for the following structures (with PDB IDs) have been deposited in the PDB: wild-type Ago2 bound to target RNAs bearing t1C (4Z4C), t1G (4Z4D), t1U (4Z4E), t1DAP (4Z4F), t1I (4Z4G); and, A481T Ago2 bound to target RNAs with t1A (4Z4H) and t1G (4Z4I).

Reviewing Editor

  1. Phillip D Zamore, Howard Hughes Medical Institute, University of Massachusetts Medical School, United States

Publication history

  1. Received: March 21, 2015
  2. Accepted: September 10, 2015
  3. Accepted Manuscript published: September 11, 2015 (version 1)
  4. Version of Record published: October 15, 2015 (version 2)

Copyright

© 2015, Schirle et al.

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.

Metrics

  • 2,632
    Page views
  • 771
    Downloads
  • 65
    Citations

Article citation count generated by polling the highest count across the following sources: Scopus, Crossref, PubMed Central.

Download links

A two-part list of links to download the article, or parts of the article, in various formats.

Downloads (link to download the article as PDF)

Open citations (links to open the citations from this article in various online reference manager services)

Cite this article (links to download the citations from this article in formats compatible with various reference manager tools)

  1. Nicole T Schirle
  2. Jessica Sheu-Gruttadauria
  3. Stanley D Chandradoss
  4. Chirlmin Joo
  5. Ian J MacRae
(2015)
Water-mediated recognition of t1-adenosine anchors Argonaute2 to microRNA targets
eLife 4:e07646.
https://doi.org/10.7554/eLife.07646

Further reading

    1. Biochemistry and Chemical Biology
    Allyson Li, Rashmi Voleti ... Neel H Shah
    Tools and Resources Updated

    Tyrosine kinases and SH2 (phosphotyrosine recognition) domains have binding specificities that depend on the amino acid sequence surrounding the target (phospho)tyrosine residue. Although the preferred recognition motifs of many kinases and SH2 domains are known, we lack a quantitative description of sequence specificity that could guide predictions about signaling pathways or be used to design sequences for biomedical applications. Here, we present a platform that combines genetically encoded peptide libraries and deep sequencing to profile sequence recognition by tyrosine kinases and SH2 domains. We screened several tyrosine kinases against a million-peptide random library and used the resulting profiles to design high-activity sequences. We also screened several kinases against a library containing thousands of human proteome-derived peptides and their naturally-occurring variants. These screens recapitulated independently measured phosphorylation rates and revealed hundreds of phosphosite-proximal mutations that impact phosphosite recognition by tyrosine kinases. We extended this platform to the analysis of SH2 domains and showed that screens could predict relative binding affinities. Finally, we expanded our method to assess the impact of non-canonical and post-translationally modified amino acids on sequence recognition. This specificity profiling platform will shed new light on phosphotyrosine signaling and could readily be adapted to other protein modification/recognition domains.

    1. Biochemistry and Chemical Biology
    2. Structural Biology and Molecular Biophysics
    Sean M Braet, Theresa SC Buckley ... Ganesh S Anand
    Research Article Updated

    SARS-CoV-2 emergent variants are characterized by increased viral fitness and each shows multiple mutations predominantly localized to the spike (S) protein. Here, amide hydrogen/deuterium exchange mass spectrometry has been applied to track changes in S dynamics from multiple SARS-CoV-2 variants. Our results highlight large differences across variants at two loci with impacts on S dynamics and stability. A significant enhancement in stabilization first occurred with the emergence of D614G S followed by smaller, progressive stabilization in subsequent variants. Stabilization preceded altered dynamics in the N-terminal domain, wherein Omicron BA.1 S showed the largest magnitude increases relative to other preceding variants. Changes in stabilization and dynamics resulting from S mutations detail the evolutionary trajectory of S in emerging variants. These carry major implications for SARS-CoV-2 viral fitness and offer new insights into variant-specific therapeutic development.