Structural basis for terminal loop recognition and stimulation of pri-miRNA-18a processing by hnRNP A1 AbstractPost-transcriptional mechanisms play a predominant role in the control of microRNA (miRNA) production. Recognition of the terminal loop of precursor miRNAs by RNA-binding proteins (RBPs) influences their processing; however, the mechanistic basis for how levels of individual or subsets of miRNAs are regulated is mostly unexplored. We previously showed that hnRNP A1, an RBP implicated in many aspects of RNA processing, acts as an auxiliary factor that promotes the Microprocessor-mediated processing of pri-mir-18a. Here, by using an integrative structural biology approach, we show that hnRNP A1 forms a 1:1 complex with pri-mir-18a where both RNA recognition motifs (RRMs) bind to cognate RNA sequence motifs in the terminal loop of pri-mir-18a. Terminal loop binding induces an allosteric destabilization of base-pairing in the pri-mir-18a stem that promotes its downstream processing. Our results highlight terminal loop RNA recognition by RBPs as a potential general principle of miRNA biogenesis and regulation. IntroductionMicroRNAs (miRNAs, miRs) are a class of highly conserved small non-coding RNAs that play a crucial role in the regulation of gene expression. They are involved in a variety of biological processes including cell growth, proliferation, and differentiation1. Mature miRNAs are generated by two RNA cleavage steps involving nuclear and cytoplasmic RNase III enzymes (Drosha and Dicer, respectively). Primary miRNA (pri-miRNA, pri-mir) transcripts are cropped by the Microprocessor complex (comprising Drosha and DGCR8) in the nucleus giving rise to ~70 nucleotide (nt) stem-loop precursor miRNAs (pre-miRNAs, pre-mir), which following export to the cytoplasm, are further processed by Dicer into mature miRNAs2. In higher organisms, many miRNA genes are transcribed together as a cluster3. A prototypical example is the miR-17鈥?2 cluster that is encoded as an intronic polycistron comprising six individual miRNAs that are highly conserved in vertebrates (miR-17, miR-18a, miR-19a, miR-20a, miR-19b-1, miR-92a-1)4,5. This cluster is frequently amplified and overexpressed in human cancers; hence, it is also referred to as oncomiR-1.The biogenesis of miRNAs is tightly regulated and results in tissue-specific and developmental-specific expression patterns of miRNAs6. A number of specific RNA-binding proteins (RBPs) have recently emerged as important post-transcriptional regulators of miRNA processing. However, very little is known about their mechanism of action. Previously, we identified heterogenous nuclear ribonucleoprotein A1 (hnRNP A1) as an auxiliary factor聽that positively regulates the processing of miRNA-18a primary transcript (pri-mir-18a) by making specific contacts to the terminal loop of the RNA7,8 (Fig.聽1a, b). HnRNP A1 is a highly abundant RBP that has been implicated in diverse cellular functions related to RNA processing, including alternative splicing regulation9,10, mRNA export11, IRES (internal ribosome entry site)-mediated translation12, mRNA stability13, and telomere maintenance14,15.Fig. 1The tandem RRMs of hnRNP A1 promote pri-mir-18a biogenesis in living cells. a Domain structure of human hnRNP A1. The sequences of conserved RNP-1 and RNP-2 motifs in the RRM domains are indicated. M9 is a transport signal linked with both nuclear import and export of this protein. b Secondary structure of pri-mir-18a RNA based on footprinting analysis7. Regions corresponding to the terminal loop and stem are boxed. The cleavage sites for Microprocessor (Drosha/DGCR8) are indicated by arrowheads. c Effect of transiently transfected epitope-tagged T7-hnRNP A1, UP1, UP1-M9, RRM1-M9, and RRM2-M9 in the processing of pri-mir-18a in HeLa cells in culture. Processing of pri-mir-16 was included as a control (white bars). The upper panel shows the level of expression of T7 epitope-tagged hnRNP A1 WT or constructs expressing individual domains (RRM1 or RRM2) or the UP1 fragment (tandem RRM1-RRM2). An M9 sequence was included to direct the nuclear localization of the UP1, RRM1, and RRM2 constructs. d Electrophoretic mobility shift assay (EMSA) of UP1 in complex with pri-mir-18a loop and stem RNAsFull size imageHere, we have combined an integrative structural biology approach with biochemical and functional assays to provide mechanistic insights into the role of hnRNP A1 in stimulating pri-mir-18a processing. We show that hnRNP A1 forms a 1:1 complex with pri-mir-18a in solution, with the recognition of two UAG motifs by the tandem RRM domains of hnRNP A1 revealed by a high-resolution crystal structure. NMR and biophysical data show that high-affinity binding involves recognition of two UAG motifs in the pri-mir-18a terminal loop and the proximal stem region. Notably, binding to the terminal loop induces an allosteric destabilization of base-pairing in the pri-mir-18a stem that promote its processing. These findings may serve as a paradigm for the regulation of miRNA processing by the recognition of the terminal loop by RBPs.ResultsUP1 stimulates pri-mir-18a processingWe have previously shown that hnRNP A1 acts as an auxiliary factor that promotes miR-18a biogenesis; however, the underlying molecular mechanisms are unknown7. HnRNP A1 has two RNA recognition motif (RRM) domains, each harboring conserved RNP-1 and RNP-2 submotifs, and a C-terminal flexible glycine-rich tail, which includes the M9 sequence responsible for nuclear import and export11,16. The RNA-binding region of hnRNP A1, comprising the tandem RRM1-RRM2 domains, is referred to as UP1 (Unwinding Protein 1)17 (Fig.聽1a). To identify which regions in hnRNP A1 are required for stimulating pri-mir-18a processing in living cells we used an in vivo processing assay. For this, several N-terminal T7-tagged hnRNP A1 constructs were transiently overexpressed in HeLa cells and the level of mature miR-18a was analyzed by qRT-PCR. We found that overexpression of full-length hnRNP A1 results in a ~two-fold increase in the levels of mature miR-18a, whereas UP1, comprising both RRMs but lacking the M9 sequence, has no effect鈥攄ue to its cytoplasmic localization (Fig.聽1c; Supplementary Fig.聽1a, b; Supplementary Table聽1). This was confirmed by transient expression of UP1-M9, where UP1 is fused to the M9 sequence that directs nuclear localization of hnRNP A116. UP1-M9 localized exclusively to the nucleus (Supplementary Fig.聽1b; Supplementary Table聽1) and, importantly, displayed similar activity as full-length hnRNP A1 in stimulating miR-18a production in vivo (Fig.聽1c). By contrast, RRM1-M9 and RRM2-M9 that partially localize to the nucleus (Supplementary Fig.聽1b; Supplementary Table聽1) did not increase miR-18a levels, showing that nuclear localized UP1 is necessary and sufficient to enhance pri-mir-18a processing (Fig.聽1c). To assess the effect of RNA-binding in the functional activity of UP1, we mutated conserved Phe residues within RNP motifs of the two RRM domains that directly contact the RNA and are required for the RNA-binding activity in hnRNP A117. Substitution of Phe with Asp within RRM2 (UP1-M9_FD2) or with Asp or Ala in both RRM1 and RRM2 domains (UP1-M9_FD12, UP1-M9_FD12a, or UP1-M9_FA12b) is sufficient to abolish the activity of UP1-M9 in our in vivo pri-mir-18a processing assay without affecting the nuclear localization of the protein constructs (Supplementary Fig.聽1a, b; Supplementary Table聽1). This indicates that the RNA-binding activity of hnRNP A1 is essential for enhancing miR-18a biogenesis. The in vivo functional data further confirm that the stimulatory function of hnRNP A1 in processing of pri-mir-18a requires both RRM domains, as mutations that affect RNA binding in one domain (or deletion of one domain) abolish the activity of hnRNP A1.UP1 specifically recognizes the loop region of pri-mir-18aTo determine the regions of pri-mir-18a that are recognized by hnRNP A1 we performed electro-mobility shift assays (EMSA) with RNA variants corresponding to the terminal loop and stem of pri-mir-18a (Fig.聽1b) and UP1, which has been shown to recapitulate most of the functions of full-length hnRNP A1 in vitro17 (Fig.聽1d). We observed that UP1 specifically binds to the terminal loop RNA, whereas no binding to the stem RNA was detected even at higher protein-to-RNA ratios. Single RRM1 and RRM2 domains do not show any detectable RNA-binding activity in this assay (Supplementary Fig.聽1c). Altogether, these data show that (i) UP1 binds specifically to the terminal loop region of pri-mir-18a, and (ii) both RRM domains of UP1 are required for high-affinity RNA binding, indicating that they bind cooperatively.Minimal UP1-binding element in the pri-mir-18a terminal loopTo provide a quantitative analysis of binding affinities of UP1 with pri-mir-18a, we performed isothermal titration calorimetry (ITC) experiments. Pri-mir-18a 71-mer RNA binds to UP1 with a dissociation constant (KD) of 147鈥塶M (Fig.聽2a, right panel). To identify a minimal RNA region required for efficient UP1 binding, we tested RNA fragments of various sizes for binding to UP1 and individual RRMs by ITC (Table聽1; Fig.聽2a; Supplementary Fig.聽2a). A 7-mer oligonucleotide (5鈥?AGUAGAU 3鈥? corresponding to the terminal loop of pri-mir-18a harbors an UAG motif, which is known to be recognized by hnRNP A118. This 7-mer single-stranded RNA binds to the individual RRM1 and RRM2 domains with low micromolar affinity (KD鈥?鈥?0.4鈥壜礛 and 6.9鈥壜礛, respectively) forming 1:1 complexes (Table聽1; Supplementary Fig.聽2a). Binding of the 7-mer to UP1 has a KD of 3.4 碌M (Table聽1; Fig.聽2a, left panel), suggesting that each RRM domain in UP1 can recognize the 7-mer RNA. Notably, a single-stranded 12-mer oligonucleotide (5鈥?AGUAGAUUAGCA 3鈥? derived from the pri-mir-18a terminal loop and flanking sequences (Fig.聽1b) shows more than 200-fold higher affinity (KD鈥?鈥?5.5鈥塶M with a 1:1 stoichiometry) to UP1 (Fig.聽2a, middle panel). The 12-mer RNA harbors two UAG motifs suggesting that each of these can be recognized by one of the RRMs in a cooperative manner. This is evident from the very large increase in binding affinity compared to the binding of the 7-mer to UP1. It is remarkable that binding of UP1 to this single-stranded 12-mer RNA is at least 10-fold stronger than binding to full-length pri-mir-18a (KD鈥?鈥?5.5鈥塶M vs. 147鈥塶M, respectively) (Fig.聽2a, middle and right panel, respectively). This may be explained by considering that in the pri-mir-18a stem-loop the second UAG motif is predicted to be base-paired in the loop-proximal stem and thus not freely accessible (Fig.聽1b). Binding of UP1 to pri-mir-18a may thus require opening (melting) of these base pairs, whereas in the single-stranded 12-mer RNA both UAG binding sites are readily available for interaction with UP1 (see below), hence the higher affinity. We conclude that the 12-mer RNA is a high-affinity UP1-binding sequence where both RRM domains of UP1 recognize UAG motifs. The RNA recognition features in the 12-mer are expected to represent the interaction of UP1 with the pri-mir-18a. A summary of thermodynamic parameters of the hnRNP A1鈥揜NA interactions is given in Table聽1.Fig. 2Biophysical characterization of UP1鈥揜NA interactions. a ITC of the binding of UP1 to 7-mer, 12-mer, and pri-mir-18a RNAs. KD values are indicated. The sequences of 7-mer and 12-mer oligonucleotides derived from the terminal loop of pri-mir-18a are shown on the left. b Combined 1H and 15N chemical shift perturbations (CSPs) of RRM1/7-mer, RRM2/7-mer, UP1/7-mer, and UP1/12-mer are plotted against the residue number. Secondary structure elements are shown above the plot. The gaps in the graph are proline residues; negative gray bars represent residues that could not be assigned due to line-broadening. c, d 1H, 15N correlation spectra of UP1, free (black) and in the presence of either the (c) 7-mer or (d) the 12-mer RNA (red). Selected residues experiencing large chemical shift perturbations and line-broadening upon RNA binding are labeled in green and black, respectively. e The CSPs are mapped onto the structure of UP1 (gradient of white to blue indicates weak to strong CSPs). Residues corresponding to amide signals that are exchange-broadened in the RNA-bound spectra are colored redFull size imageTable 1 Isothermal titration calorimetry of protein鈥揜NA interactionsFull size tableNMR analysis of UP1鈥揜NA interactionsWe next characterized the binding interface of UP1 with various RNA ligands derived from pri-mir-18a using NMR titration experiments. Addition of the 7-mer RNA harboring one UAG motif causes extensive chemical shift perturbations (CSPs) in RRM1 and RRM2 constructs, thus demonstrating that each of the RRMs can interact with the 7-mer RNA (Fig.聽2b). The CSP pattern obtained upon titration of the tandem RRM domains in UP1 with the 7-mer RNA is very similar to that obtained for individual RRM domains (Fig.聽2b, c). This shows that recognition of the 7-mer RNA in the isolated RRM domains and in the context of the UP1 construct is very similar. Similarly, the 12-mer RNA harboring two UAG motifs induces large CSPs in both RRM domains of UP1 (Fig.聽2d) that are comparable to those obtained at saturating levels of the 7-mer. The CSPs map to the canonical RNA-binding surface on the 尾-sheets of the two RRM domains (Fig.聽2e). In addition, strong CSPs are observed for residues in the C-terminal region of RRM2. This region is flexible in free UP1 but upon RNA-binding forms an additional helix (伪3), which is not present in the free protein19 (Supplementary Fig.聽3a, b), suggesting that helix 伪3 is induced and stabilized upon RNA binding. Interestingly, a number of NMR signals that correspond to residues in RRM2 and the RRM1-RRM2 linker are severely broadened in the RNA-bound spectrum, suggesting dynamics on the 碌s-ms time-scale. The affected residues map to the interface between the RRM1 and RRM2 domains (Fig.聽2e), suggesting that some conformational dynamics and adaptation of this interface is associated with RNA binding.Recognition of the pri-mir-18a terminal loop by UP1To gain insight into the molecular details of pri-mir-18a recognition by UP1, we determined the crystal structure of the UP1/12-mer RNA complex at 2.5鈥壝?resolution (Table聽2; Fig.聽3a; Supplementary Fig.聽4a). Surprisingly, the crystal structure of the UP1/12-mer RNA complex exhibits two molecules of UP1 and two RNA chains in the asymmetric unit in a 2:2 stoichiometry (Supplementary Fig.聽4a), similar to a previously reported structure of UP1 with single-stranded telomeric DNA20. As this peculiar 2:2 stoichiometry most likely does not represent the UP1:RNA complex in solution, we analyzed UP1, RNA and the protein鈥揜NA complexes by static light scattering (SLS) (Fig.聽3b). Both UP1 and pri-mir-18a alone behave as single species with a molecular weight corresponding to respective monomeric conformations. Importantly, the molecular weight obtained for the UP1/12-mer complex (26.9鈥塳Da) indicates a 1:1 complex (theoretical molecular weight: 26.3鈥塳Da) and demonstrates that the 2:2 stoichiometry observed in the crystal structure is an artifact induced by the crystal environment. Notably, also the molecular weight obtained for the UP1/pri-mir-18a complex (45鈥塳Da) is fully consistent with the formation of a 1:1 complex (Fig.聽3a, b).Table 2 Crystallographic data collection and refinement statisticsFull size tableFig. 3Structure of the UP1/12-mer complex. a Schematic representation of the 2:2 UP1/12-mer RNA complex observed in the crystal structure (Supplementary Fig.聽4a) and the proposed 1:1 complex in solution. b Static light scattering (SLS) profiles of UP1, pri-mir-18a, UP1/12-mer, and UP1/pri-mir-18a. The molecular weight (MW) obtained from SLS is indicated. UP1 forms a 1:1 complex with both pri-mir-18a and the 12-mer oligonucleotide. c Structural model of the 1:1 UP1/12-mer complex, where each RRM domain recognizes a UAG motif in the 12-mer RNA (magenta) forming a 1:1 complex. d Structural details of recognition of the UAG motif by RRM1 and RRM2 domains. Side-chain of selected residues involved in the interaction are indicated. e, f Correlation of experimental 1H-15N residual dipolar couplings (RDCs) (e) and small angle X-ray scattering (SAXS) data (f) vs. those calculated from the structural model shown in cFull size imageThe crystal structure reveals that each RRM domain in UP1 specifically recognizes one UAG motif in the RNA. Although the stoichiometry does not reflect the solution conformation, the RNA-binding interface is expected to be conserved. This is consistent with the NMR chemical shift perturbations which map to the canonical binding interface, although some minor differences may exist regarding specific contacts considering the overall different topology. The crystal structure shows that the central AG dinucleotide in each UAG motif is recognized by contacts mainly through conserved RNP motif residues on the 尾-sheets (Fig.聽3c, d), which resembles the recognition of TAG in the UP1-telomeric DNA complex20 and recently reported structures with RNA21,22. Two conserved aromatic residues in RRM1, Phe17 (RNP-2 motif residue located on 尾1) and Phe59 (RNP-1 motif residue located on 尾3), are involved in stacking interactions with the bases of A4 and G5, respectively (Fig.聽3d). Similarly, in RRM2 Phe108 (RNP-2 motif residue located on 尾1) and Phe150 (RNP-1 motif residue located on 尾3) stack with the bases of A9 and G10, respectively. In addition to the stacking interactions with RNP residues (vide infra), the adenosine is specifically recognized by hydrogen bonds of its exocyclic NH2 group with the main chain carbonyl oxygen of residues Arg88 and Lys179 in RRM1 and RRM2, respectively. A positively charged residue in each domain, Arg55 in RRM1 and Arg146 in RRM2, makes electrostatic interactions with the phosphate backbone of the AG dinucleotide. Two charged residues in each domain, Glu85 and Lys87 in RRM1 and Glu176 and Arg178 in RRM2, contact the uridines in the UAG motifs (U3 and U8), while another charged residue, Lys15 in RRM1 and Lys106 in RRM2, interacts with G5 and G10, respectively (Fig.聽3d), thereby recognizing the U and G residues in the UAG motif. Overall, the mode of RNA recognition by RRM1 and RRM2 is very similar; in each domain an AG dinucleotide is sandwiched between the 尾-sheet surface and a C-terminal helix (Fig.聽3d).The two RRM domains of hnRNP A1 are connected by an approximately 17-residue linker, which is evolutionary conserved both in terms of sequence and length23 (Supplementary Fig.聽3e). Residues in this linker are highly flexible as determined by NMR relaxation data (Supplementary Fig.聽3a). However, the two RRM domains interact with each other and this intramolecular domain interface is stabilized by two conserved salt bridges (Arg75鈥揂sp155 and Arg88鈥揂sp157) as well as a small cluster of hydrophobic residues in the interface. Virtually the same domain arrangement and salt bridges are observed in the monomeric solution structure of UP1 free19 and a crystal structure bound to a short RNA22 (Supplementary Fig.聽4f), indicating that this arrangement of RRM1 and RRM2 stabilized by a number of interactions is preformed and may be retained in free and ligand-bound structures of UP1.To confirm the domain arrangement of RNA-bound UP1 in solution, we measured NMR and small angle X-ray scattering (SAXS) (Supplementary Fig.聽4b鈥揹). NMR paramagnetic relaxation enhancements (PRE) were recorded for UP1 spin-labeled at position 66 (UP1 Glu66Cys mutant). The PRE data provide long-range (up to ~20鈥壝? distance information and can thus report on domain/domain arrangements24,25,26 (Supplementary Fig.聽4c). The PRE profiles of free and RNA-bound form of UP1 are similar, suggesting that the domain arrangement does not change significantly in the presence of the 12-mer RNA.To determine a structural model of the 1:1 UP1/12-mer RNA complex in solution a protein monomer from our crystal structure was used as template. The RNA was initially placed by constraining the recognition of the two UAG motifs in the 12-mer RNA (taken from two different RNA molecules in the 2:2 crystal structure) (Fig.聽3a). SAXS and NMR residual dipolar coupling (RDC) data were then used as restraints in a molecular dynamics simulation together with the distances of the Arg75鈥揂sp155 and Arg88鈥揂sp157 salt bridges as observed in the crystal structure. The refinement shows that a 1:1 complex is fully compatible with the data (Fig.聽3e, f; Supplementary Table聽2). The arrangement of RRM1 and RRM2 is similar to the corresponding interface in the crystal structure (coordinate rmsd of 2.8鈥壝?calculated over all the backbone atoms excluding linker residues) with only minor conformational changes in the domain/domain interface (Fig.聽3c; Supplementary Fig.聽4a鈥揺).Validation of the UP1/12-mer RNA structural modelThe structural model of the 1:1 UP1/12-mer RNA complex was confirmed by mutational analysis of protein and RNA. ITC data with 12-mer RNAs where the first or second UAG motif has been replaced by UUU, 12-mer-mut1: AGUUUAUUAGCA and 12-mer-mut2: AGUAGAUUUUCA show 10-fold and 20-fold (KD鈥?鈥?54鈥塶M and 330鈥塶M) reduced binding affinity, respectively, compared to the wild-type sequence (Table聽1; Supplementary Fig.聽2b). This confirms that both motifs are recognized by the protein. Additional RNA variants with an AG鈫扷U mutation or lacking the initial AG dinucleotide, that precedes the two 鈥渃anonical鈥?UAG motifs, have the same binding affinity as the wild-type 12-mer RNA (Table聽1; Supplementary Fig.聽2b). As mutation or deletion of the AG dinucleotide at the very 5鈥?end of the 12-mer RNA has essentially no effect, this establishes that UP1 has a preference for the recognition of two closely spaced UAG motifs.Further, the domain interface in UP1 was probed by mutations that are expected to disrupt the two salt bridges Arg75鈥揂sp155 and Arg88鈥揂sp157 in the RRM1-RRM2 interface, which have been observed in all reported structures of free and nucleic acid bound forms of UP119,20,21,22. Introducing charge clashes (UP1-Arg75Glu/Arg88Glu) in this interface, which is remote from the RNA-binding surface, decreases the binding affinity to the 12-mer RNA by 鈮?-fold (KD ~15.5鈥塶M vs. ~40鈥塶M) (Supplementary Fig.聽2c). This suggests that the salt bridges play an indirect role for RNA binding by stabilizing the arrangement of the two RRM domains.Functionally, we established that mutations of conserved Phe residues in the RNP motifs of RRM2 or both RRMs that affect RNA binding, consequently abrogate the stimulatory effect of UP1-M9 in the processing of pri-mir-18a in HeLa cells in culture (Supplementary Fig.聽1a, b; Supplementary Table聽1). Collectively, the structural model of the 1:1 UP1/12-mer RNA complex is fully consistent with our biochemical and functional data regarding the requirement for two RRM domains and two UAG motifs for high-affinity interaction.UP1 binding destabilizes pri-mir-18a RNA near the terminal loopThe structure of UP1 bound to the 12-mer RNA in a 1:1 stoichiometry implies that the loop-proximal region of the pri-mir-18a stem should be destabilized to enable recognition of two UAG motifs, one accessible in the terminal loop and the second one as part of the pri-mir-18a duplex. To study this further, we first analyzed the structure of pri-mir-18a alone. For this, a model of the RNA was prepared using the MC-Fold/MC-Sym server27 (Fig.聽4a) and assessed with experimental SAXS and NMR data. The predicted secondary structure and elongated shape of pri-mir-18a 71-mer is supported by SAXS data of the free RNA (Fig.聽4b). We then used NMR to analyze the base-pairing in the pri-mir-18a 71-mer using imino NOESY spectra (Fig.聽4c), as imino proton NMR signals probe the presence and stability of base pairs28. We could unambiguously assign imino-imino cross-peaks corresponding to the stem region of the 71-mer RNA. However, no imino correlation was observed for the upper part of the stem loop (Fig.聽4c, see secondary structure on the right), suggesting that this region of the RNA helix is dynamic with less stable base-pairing.Fig. 4Structural analysis of pri-mir-18a 71-mer RNA. a A structural model of pri-mir-18a 71-mer RNA obtained using the MC-Fold/MC-Sym server. b Experimental and predicted SAXS data for the pri-mir-18a model are shown in black and red, respectively. c NMR analysis of the 71-mer pri-mir-18a RNA. 1H, 15N-HSQC of 13C, 15N-labeled pri-mir-18a in the absence (black) and presence (red) of UP1. Imino signals observed indicating stable base-pairing are shown on the secondary structure of pri-mir-18a (magenta boxes). The signal marked with an asterisk has not been assigned. Nucleotides undergoing large perturbations and/or line broadening upon UP1 binding are highlighted with green boxes (left panel) or green letters and arrows (right panel). d 2D-imino NOESY and e H5 correlated HNN experiment of the 17-mer RNA derived from the terminal loop and flanking stem of pri-mir-18a. Base pairs confirmed by these experiments are boxed in blue and red on the secondary structure of the 17-mer RNA. The first G:C base pair, which is not part of the native sequence is shown in lower case. f 2D-imino NOESY of the 17-mer(A35C) RNA. The mutated nucleotide is shown in red. Base pairs confirmed by the detection of the imino correlation are boxedFull size imageThis was further confirmed by analyzing a 17-mer stem-loop construct (Fig.聽4d, e), which represents the loop-proximal stem region in the full-length pri-mir-18a (Fig.聽1b). The upper stem region in this RNA harbors the second UAG motif that is also present in the single-stranded 12-mer RNA and recognized by UP1. An imino NOESY spectrum of the 17-mer RNA shows imino correlations only for the G:C and G:U base pairs at the stem of this RNA. The fact that the imino protons of the predicted A:U base pairs adjacent to the terminal loop are not detectable indicates that they are dynamic (Fig.聽4d). Nevertheless, the presence of these A:U base pairs is confirmed by the (Py) H(CC)NN-COSY experiment, which allows detection of transient and weak-pairs, even in the absence of observable imino proteins29 (Fig.聽4e). The absence of detectable imino protons for many predicted base pairs in the 17-mer and in the upper part of pri-mir-18a indicates that the corresponding stem region exhibits weak base pairs and is dynamic.The observation of dynamic base-pairing in the loop-proximal region is consistent with melting of the stem region to enable recognition of both UAG motifs by UP1. In fact, NMR titrations of UP1 with the 17-mer stem-loop RNA show virtually identical CSPs to those seen with the single-stranded 12-mer RNA (Fig.聽2d, Supplementary Fig.聽3c). As NMR chemical shifts are sensitive indicators of the three-dimensional structure, this finding confirms that the protein鈥揜NA interactions observed with the 12-mer RNA also reflect the RNA recognition within the 17-mer RNA. In both cases, two UAG motifs are involved in the interaction, thus requiring melting of the 17-mer RNA helical stem. This is also consistent with the ITC binding curve (Supplementary Fig.聽2d), which shows unusual concentration-dependent enthalpy changes that cannot be fitted to a simple binding event, but may reflect coupled RNA unfolding and protein binding. To further support this, we studied an A35C variant of the 17-mer RNA, which introduces complete Watson鈥揅rick complementarity in the RNA stem. In contrast to the 17-mer, the 17-mer (A35C) exhibits all expected base pairs as evidenced by detectable imino protons throughout the stem region (Fig.聽4f). NMR titrations show much smaller chemical shift perturbations in UP1 for the mutant 17-mer RNA (Supplementary Fig.聽3d), which binds with a KD of ~3鈥壜礛, i.e., corresponding to a 20鈥?00-fold reduced binding affinity compared to the full-length pri-mir-18a and the single-stranded 12-mer RNA (Table聽1; Supplementary Fig.聽2d). These findings can be rationalized by considering that in the A35C mutant only one single-stranded UAG motif is available for binding to UP1, as the second motif is now part of a highly stable RNA helix and thus not accessible. Consistently, the 17-mer (A35C) exhibits micromolar affinity to UP1, i.e., comparable to the affinity observed for the 7-mer RNA with a single UAG motif. On the other hand, the fact that the wild-type pri-mir-18a has one UAG motif in the terminal loop and a second one in the dynamic stem region suggests that binding of UP1 requires melting of the stem region flanking the terminal loop. Taken together, these data suggest that melting of the loop-proximal stem region in pri-mir-18a enables recognition of the second UAG motif and thus provides high-affinity (low nanomolar KD) binding of UP1 to both UAG motifs.Next, we compared the 1H, 15N imino correlations of the pri-mir-18a 71-mer RNA in the free form and bound to UP1 (Fig.聽4c). Notably, signals for the U10:G52 base pair in the middle of the stem region are not detectable in the complex but readily observed in the free form, indicative of destabilization and partial melting of this part of the duplex. Also, other residues especially next to mismatches and other less stable regions of the RNA exhibit reduced intensity or chemical shift perturbation such as G23, U49, and U59 (Fig.聽4c, nucleotides and arrows in green). This is consistent with allosteric effects that lead to destabilization of the complete pri-mir-18a stem-loop induced by binding of UP1 to the terminal loop.Structural model of the UP1/pri-mir-18a complexTo derive a structural model of UP1 bound to pri-mir-18a, we performed molecular dynamics simulations restrained by the structural information obtained for the UP1/12-mer RNA complex and experimental SAXS data of the UP1/pri-mir-18a complex (Fig.聽5a; see Methods for details). The resulting structure was then scored against a combination of small angle neutron scattering (SANS) data with contrast matching (Fig.聽5a; Supplementary Fig.聽4g, h). SAXS data show that the dimensions of the complex correspond to a radius of gyration (Rg) and maximum pairwise distance (Dmax) values of 37.7鈥壝?and 130鈥壝? respectively, consistent with a 1:1 complex (Supplementary Table聽2). The structural model shows very good agreement with all experimental data and is also consistent with ab initio SAXS derived models of the protein鈥揜NA complex (Supplementary Fig.聽4g, h). The UP1/pri-mir-18a complex shows the recognition of two UAG motifs in the terminal loop and a partially melted upper stem of pri-mir-18a (Fig.聽5a).Fig. 5hnRNP A1 recognizes the terminal loop of pri-mir-18a. a Structural model of the 71-mer pri-mir-18a/UP1 complex. The region corresponding to the 12-mer RNA in the terminal loop of the pri-mir-18a is shown in dark magenta. b Footprint analysis of the pri-mir-18a/UP1 complex. Cleavage patterns were obtained for 5鈥?32P-labeled pri-mir-18a transcript (100鈥壝椻€?03 c.p.m.) incubated in the presence of recombinant UP1 protein (+200鈥塶g, 500鈥塶M), treated with Ribonuclease T1 at 1.5鈥塙/渭L. F and T identify nucleotide residues subjected to partial digest with formamide (every nucleotide) or ribonuclease T1 (G-specific cleavage), respectively. The cleavages intensities generated by Ribonuclease T1 are indicated on the pri-mir-18a secondary structure. The region of the major UP1 footprints is indicated by a blue oval shape. Residues that show enhanced reactivity upon UP1 binding are indicated by green arrows. c A schematic model of the effects of hnRNP A1 on pri-mir-18a. Binding of UP1 to the terminal loop, where each RRM domain recognizes a UAG motif, leads to melting of the terminal loop that then spreads down leading to destabilization and enhanced dynamics in the RNA stem, thereby facilitating Drosha cleavage by unknown mechanisms. d Averaged SHAPE and NMR data were used to calculate the secondary structure of pri-mir-18a free and when bound to UP1. SHAPE data were implemented in RNAstructure31 as a soft constraint, while base pairs identified by NMR were used as hard constraints. The structures shown correspond to the lowest free energy models. Color coding corresponds to SHAPE reactivity as indicated with red/orange/yellow corresponding to high/semi/low reactivityFull size imageTo validate the structural model described and to assess changes in accessibility of the pri-mir-18a induced by UP1 binding, we performed foot-printing analysis of the complete pri-mir-18a RNA in the absence and presence of UP1 (Fig.聽5b). This revealed that in the free RNA the terminal loop and flanking stem region comprising the two UAG motifs are accessible and dynamic (Fig.聽5b), which is consistent with the NMR data (Fig.聽4c, d). Binding of UP1 protects this region in the RNA, in full agreement with the structural model of the UP1/pri-mir-18a complex. The significant reduction in accessibility observed for residues in the terminal loop is consistent with the simultaneous binding of both RRM domains to the RNA (Figs.聽3c and 5a). Interestingly, residues at the bottom part of the stem of pri-mir-18a become more accessible for nuclease cleavage upon binding of UP1 to the terminal loop (Fig.聽5b, green nucleotides). This is fully consistent with the NMR analysis which indicates that UP1 binding leads to destabilization of the pri-mir-18a stem (Fig.聽4c).SHAPE analysis of pri-mir-18aTo assess the RNA structure of pri-mir-18a and potential effects of flanking regions in the context of the pri-mir-17-18a-19a cluster, we performed selective 2鈥?hydroxyl acylation analyzed by primer extension (SHAPE, Supplementary Fig.聽5a, b)30. In vitro-transcribed RNAs comprising either pri-mir-18a or the pri-mir-17-18a-19a cluster were incubated with increasing amounts of purified UP1 or full-length hnRNP A1 proteins, prior to treatment with N-methylisatoic anhydride (NMIA). The SHAPE reactivity reflects the average intensity of the NMIA-treated RNA primer extension products, normalized to the corresponding untreated RNA, in the presence or absence of UP1/hnRNP A1 proteins.SHAPE analysis of pri-mir-18a reveals a region of medium to high nucleotide flexibility, corresponding to the terminal loop (nts 170鈥?86). As indicated by NMR (compare Fig.聽4d, e), this region of the pri-mir-18a stem loop is dynamic, consistent with varying degrees of SHAPE reactivity for the loop-proximal base pairs, thus supporting the dynamic nature of the upper stem loop. The secondary structure of the free pri-mir-18a was predicted using both SHAPE reactivity and NMR-identified base pairs with RNAstructure31 (Fig.聽5d). The structure is in excellent agreement with the SAXS, enzymatic probing, and NMR data.SHAPE analysis of pri-mir-17-18a-19a revealed few differences between SHAPE reactivity of pri-mir-18a in the cluster and isolated pri-mir-18a, indicating that the pri-mir-18a hairpin folds largely independently of the rest of the cluster (Supplementary Fig.聽5b). Despite the overall similar reactivity pattern, differences observed from nucleotide 210 through 220 are presumably induced by the presence of pri-mir-17 and pri-mir-19a in the whole transcript, which may stabilize the basal AU base pairs of pri-mir-18a. These results are corroborated by a recent SHAPE and chemical probing assays performed on the complete oncomiR-1 cluster32.Importantly, upon addition of UP1 and hnRNP A1, distinct changes in SHAPE reactivities are observed compared to the free pri-mir-18a indicating that residues around the terminal loop and flanking stem region (nts 170鈥?86) become protected by protein binding (Fig.聽5d; Supplementary Fig.聽5a). This region includes the two UAG motifs that bind to RRM1 and RRM2 of hnRNP A1/UP1. Analysis of individual replicate SHAPE data sets reveals a modest increase in SHAPE reactivity down the RNA stem in the presence of protein (nts 50, 51, 53鈥?5). However, these changes are not visible after averaging the three SHAPE replicate data sets, and are thus not statistically significant. While SHAPE may not be sensitive enough to detect the destabilization of the central stem upon hnRNP A1 binding, the SHAPE reactivity observed supports the proposed secondary structure and a plethora of evidence indicating that hnRNP A1 binds to and opens the terminal loop.Changes in SHAPE reactivity observed for the pri-mir-17-18a-19a cluster in the presence of hnRNP A1/UP1 also revealed the protection of nucleotides spanning the terminal loop of pri-mir-18a (nts 170鈥?86). Thus, most of the highly reactive residues in pri-mir-18a display a similar behavior as the pri-mir-17-18a-19a transcript upon addition of UP1/hnRNP A1 proteins (Fig.聽5d; Supplementary Fig.聽5). Importantly, both UP1 and hnRNP A1 induce comparable changes in SHAPE reactivity, confirming that UP1 is sufficient to bind to the pri-miRNA, in agreement with the EMSA and functional assays (Fig.聽1c, d).In summary, relative to the free RNA the presence of hnRNP A1/UP1 leads to a decreased SHAPE reactivity around the terminal loop of pri-mir-18a, indicating that this region is the major binding site for the RRMs of hnRNP A1.Mechanism of hnRNP A1 stimulation of pri-mir-18a processingFinally, we attempted to address the mechanism by which hnRNP A1 activates the processing of pri-mir-18a. One possible scenario is that binding of hnRNP A1 (or UP1) leads to partial opening/melting of the terminal loop, which can lead to destabilization of the stem region and thus render it more accessible for processing by Drosha7. Indeed, footprinting and site-directed mutagenesis of pri-mir-18a suggested that hnRNP A1 alters the local conformation of the stem in the vicinity of Drosha cleavage sites7, although the molecular mechanism was unclear. Indeed, UP1 (unwinding protein 1) can unwind secondary and higher order structures of DNA and RNA15,33. We constructed a series of mutants in the terminal loop region of pri-mir-18a. These include single and double nucleotide mutants, where an A residue within the UAG motif in the terminal loop or within the second UAG motif, were mutated to C (UCG in pri-mir-18a[A30C] and pri-mir-18a[A30C/A35C], respectively), and a triple mutant (pri-mir-18a[U21A/U29A/U34A]), where all UAG motifs were mutated to AAG. We also designed a mutant, in which the terminal loop was stabilized by five G:C base pairs (pri-mir-18a[5GC]). The wild-type sequence was efficiently processed in the context of the pri-mir-17-18a-19 cluster (Fig.聽6a) and isolated pri-mir-18a (Supplementary Fig.聽6). The single (pri-mir-18a [A30C]) and double nucleotide (pri-mir-18a[A30C/A35C]) terminal loop mutant RNAs retained hnRNP A1 binding, although with lower affinity. These mutant RNAs were still efficiently processed in the context of the pri-mir-17-18a-19 cluster (Fig.聽6a) and the isolated pri-mir-18a (Supplementary Fig.聽6). The triple mutant (pri-mir-18a[U21A/U29A/U34A]) showed binding to hnRNP A1, although with reduced affinity, and retained efficient processing in the pri-mir-17-18a-19 cluster (Fig.聽6b). However, it displayed decreased processing in the isolated pri-mir-18a (Supplementary Fig.聽6). The pri-mir-18a[5GC] mutant does not bind to hnRNP A1 in the RNA pull-down assay and consequently is not processed (Fig.聽6a; Supplementary Fig.聽6). This lack of processing could result from disruption of UP1 binding and/or conformational changes that inhibit Microprocessor activity, for example, by stabilizing the loop-proximal stem region. Based on these data we conclude that hnRNP A1 binding is essential for pri-mir-18a processing, although some reduction in binding affinity can still support miRNA processing, as shown by the fact that even low hnRNP A1 levels are sufficient to stimulate processing activity.Fig. 6Mechanism of hnRNP A1 stimulation of pri-miRNA processing. a Pri-mir-18a with single or double A鈫扖 point mutations still bind hnRNP A1 in RNA pull-down assays and are processed by Drosha (in vitro processing assay with the pri-mir-17-18a-19a cluster, wildtype and mutants), whereas pri-mir-18a mutant with a 5GC clamp does not bind hnRNP A1 (RNA pull-down assay is shown on the right) and is not processed by Drosha. b Pri-mir-18a with a 5GC_internal clamp and wild-type terminal loop binds hnRNP A1 but is not processed by Drosha. Pri-mir-18a with triple mutations [U21A/U29A/U34A] binds hnRNP A1 with lower affinity than the wild-type pri-mir-18a but is still efficiently processed by Drosha. c In vitro processing assay of pri-mir-16 with 5GC_internal clamp shows efficient processing by Drosha, similar to pri-mir-16 wildtypeFull size imageNext, we hypothesized that destabilization and unfolding of the upper stem region close to the terminal loop in pri-mir-18a induced by hnRNP A1 binding can spread along the RNA stem and thus lead to destabilization of the RNA near the Drosha cleavage site. To examine this hypothesis, we tested a different mutant pri-mir-18a, in which the terminal loop was clamped by 5鈥塆:C base pairs (pri-mir-18a[5GC_internal]). Note, that this RNA variant is different from the pri-mir-18a[5GC] mutant described above, in which the terminal loop region was stabilized by 5鈥塆:C base pairs. Remarkably, the processing of this RNA was impaired in the context of pri-mir-17-18a-19 cluster (Fig.聽6b) and isolated pri-mir-18a (Supplementary Fig.聽6), despite retaining full binding to hnRNP A1 in the RNA pull-down assay (Fig.聽6b). This is highly important, because it demonstrates that abrogating the allosteric change induced by hnRNP A1 binding leads to lack of stimulatory processing activity, even in the presence of high-affinity binding.To rule out the possibility that the internal 5GC clamp affects Drosha processing irrespective of the requirement for hnRNP A1 as an auxiliary factor, we examined the processing of pri-mir-16 with and without the 5GC clamp. Both wild-type pri-mir-16 and pri-mir-16[5GC_internal] were processed by Drosha, indicating that the internal 5GC clamp is not sufficient to impair Drosha processing (Fig.聽6c). Furthermore, to control for undesired effects of our mutations on Drosha cleavage we generated a set of pri-mir-16 chimerae that harbor wild-type or mutated pri-mir-18a terminal loops (Fig.聽7). Drosha efficiently processes all but one of the pri-mir-16/18a chimerae analyzed here, where only the processing of pri-mir-16/18a[5GC] was reduced. This shows that conformational changes in the apical part of the pri-mir stem can directly affect Microprocessor activity. By contrast, the processing of pri-mir-16/18a[5GC_internal]) was efficient in the context of this chimera (Fig.聽7), unlike what we observed with the processing of pri-mir-18a[5GC_internal], which was abolished (Fig.聽6b). Therefore, these data strongly support our hypothesis that binding of hnRNP A1 to the terminal loop of pri-mir-18a results in a partial unwinding of this pri-miRNA that spreads from the terminal loop towards the stem and is essential for the stimulatory role of hnRNP A1 in miR-18a biogenesis.Fig. 7In vitro processing of pri-mir-16 and 18a chimerae. Wild-type pri-mir-16 or chimerae with 18a terminal loop sequences (in green) with point mutations introduced in Supplementary Fig.聽6 (in red) were subjected to in vitro processing assay. All transcripts, except for pri-mir-16/18a[5GC], are efficiently processed by Drosha. Pri-mir-16/18a[5GC] processing efficiency is reduced, indicating that conformational changes to the pri-mir stem can directly affect Microprocessor activityFull size imageDiscussionHere, we have used a multi-disciplinary approach to reveal the molecular mechanism by which hnRNP A1 binds to pri-mir-18a and facilitates its processing. We show that hnRNP A1 specifically binds to pri-mir-18a through interactions involving both RRM domains in UP1 and a region comprising two UAG RNA sequence motifs in the terminal loop of pri-mir-18a and a flanking stem region. Cooperative binding of both domains to cognate RNA motifs results in substantially increased binding affinity and allows the unwinding of the target stem-loop RNA. This mode of binding is distinct from the recognition of a viral RNA, where apparently only RRM1 is involved in the recognition of a single AGU motif22, but reminiscent of recognition of two RNA motifs in single-stranded cis regulatory elements in alternative splicing regulation by hnRNP A121. A common feature for all structures is that the overall domain arrangement and interface in free and nucleic acid bound hnRNP A1 tandem RRMs is conserved (Supplementary Fig.聽4f). Nevertheless, some adaptation and fine-tuning of the domain arrangement occurs upon RNA binding. This is consistent with the line-broadening of amide signals in the RRM1/RRM2 interface upon binding to the 12-mer RNA (Fig.聽2e). A further common feature is the formation of 1:1 complexes in solution, very distinct from the 2:2 stoichiometry observed in crystal structures. The various distinct modes of nucleic acid recognition by hnRNP A1 are intriguing and may reflect how it can play important roles in the regulation of many distinct biological activities by its RNA interactions.hnRNP A1 can also act as a negative regulator of let-7a processing, by competing with the activator protein KSRP for binding to the pri-let-7a terminal loop34,35. In addition to hnRNP A1, several other RBPs recognize the terminal loop of miRNA precursors and influence, either positively or negatively, their biogenesis at the post-transcriptional level36. For example, binding of Lin28 to the terminal loop of pre-let-7 was shown to inhibit processing by both Drosha and Dicer37. Another example is the RBP Rbfox, which binds to the terminal loop of a subset of pri-miRNAs and suppresses their nuclear processing38.How does binding of a regulator to the terminal loop of pri-miRNA affect Drosha cleavage at the opposite end of the RNA? We propose that binding of hnRNP A1 to the terminal loop of pri-mir-18a leads to destabilization of base pairs and (partial) melting of the loop-proximal stem region, with subsequent spreading of destabilization of the stem region towards the Drosha cleavage site (Fig.聽5c). In support of this model, processing of a mutant pri-mir-18a, in which the terminal loop was clamped by 5鈥塆:C base pairs (5GC_internal) was abolished, despite strong hnRNP A1 binding (Fig.聽6b). This strongly argues that the effect of hnRNP A1 binding at the terminal loop is somehow propagated and leads to a stimulatory effect of Drosha processing. Indeed, pri-mir-18b, which is part of the homologous miR-106a-363 cluster, does not require hnRNP A1 for efficient processing8. Mechanistically, this can be explained by the fact that the conformation of the stem in pri-mir-18b resembles the more open stem structure comprising a bulge in the stem, which is only observed in pri-mir-18a upon binding of hnRNP A17,8. This affects Drosha processing as simply introducing this bulge in the pri-miR-18a stem (UC鈫扜U) made its processing more efficient and completely independent of hnRNP A17,8.We propose three non-mutually exclusive mechanisms by which hnRNP A1 may enhance the processing of pri-mir-18a. hnRNP A1 binding may 1) promote Drosha processing by stem destabilization and thereby facilitate the recruitment and/or assembly of an active Microprocessor substrate complex, 2) enhance the assembly of pri-mir-18a with the correct orientation onto the Microprocessor by generating a thermodynamic or conformational asymmetry in the stem, 3) increase the accessibility of the Microprocessor to the pri-mir-18a hairpin in the context of the complete oncomiR-1 cluster. Options 1 and 2 are consistent with recent biochemical, structural and computational studies that have shown that the Microprocessor recognizes two regions at either end of pri-miRNA, i.e., near the ssRNA/stem basal junctions at the bottom and near the terminal loop junction39,40,41,42,43,44. Recognition of specific sequences (a UG motif 5鈥?or CNNC motif in single-stranded region 3鈥?of the basal junction, and/or a UGUG motif near the terminal loop) and/or structural features at these two opposite regions of a pri-miRNA are proposed to guide the binding orientation and conformation of the Drosha/DGCR8 complex and define the cleavage site. Thus, various RBPs that recognize specific sequence motifs in terminal loops may further regulate and enhance processing of specific pri-miRNAs.Here, we found that strengthening the upper part of pri-mir-18a stem by GC base pairs blocks miRNA processing, whereas disruption of the base pairs enhances Microprocessor cleavage efficiency. It is noteworthy, that the destabilization of the apical RNA helix by binding of hnRNP A1 induces a structural and thermodynamic asymmetry in the stem region (stably base-paired vs. dynamic RNA helix) that may help define, in which orientation pri-mir-18a is recognized by the Microprocessor. It was recently proposed that heme binding by DGCR8 is required to ensure the correct Microprocessor assembly41. A similar role may thus extend to RBPs such as hnRNP A1 that bind to terminal loops of specific pri-miRNAs.It is also likely that the processing efficiency of pri-mir-18a is context-dependent, suggesting that the sequence and/or structure of pri-mir-18a as part of the miR-17鈥?2 cluster is not optimal for Drosha processing8. Interestingly, several studies have shown that the miR-17鈥?2 cluster adopts a compact tertiary structure, in which individual miRNAs have different expression levels depending on their accessibility32,45,46. Notably, a recent SHAPE analysis of the miR-17鈥?2 cluster revealed that the terminal loop of pri-mir-18a in the cluster is partially聽solvent inaccessible32. As the terminal loop corresponds to the sequence that we have identified as the main hnRNP A1-binding site in pri-mir-18a, it is tempting to speculate that binding of hnRNP A1 to the miR-17鈥?2 cluster is associated with a conformational change in the RNA, which can facilitate Microprocessor assembly and Drosha cleavage. We propose that the tertiary structure of pri-mir-18a in the context of the miR-17鈥?2 cluster, as well as sequences in the stem and loop region, are important determinants of miRNA processing by Drosha. The recognition of the pri-mir-18a stem-loop by hnRNP A1 thus adds an additional layer for the regulation of pri-miRNA processing by an RBP, beyond features that have been recently identified39,43,47.In conclusion, our data demonstrate that recognition of a conserved terminal loop RNA sequence in pri-miRNAs by an RBP can strongly modulate miRNA biogenesis by conformational changes and dynamic destabilization induced by RNA binding. This suggests that recognition of pri-miRNAs by RBPs is a general paradigm for context-dependent regulation of miRNA biogenesis and function.MethodsProtein expression and purificationThe sequences encoding RRM1 (amino acids 1鈥?7), RRM2 (amino acids 94鈥?96), and UP1 (amino acids 1鈥?96) were subcloned into the pETM11 vector (EMBL) encoding a His tag followed by TEV cleavage site (primers are listed in Supplementary Table聽3). Recombinant proteins were expressed in E. coli BL21(DE3) cells (Novagen) in standard media or minimal M9 media supplemented with 1鈥塯/L 15N-ammonium chloride and 2鈥塯/L 13C-glucose. Protein expression and purification was done as described previously10. After growth of bacterial cells up to an OD600 of 0.7鈥?.8, protein expression was induced by 0.5鈥塵M IPTG and continued overnight at 20鈥壜癈. Cells were resuspended in buffer A (30鈥塵M Tris/HCl pH 7.5, 500鈥塵M NaCl, 10鈥塵M imidazole, 1鈥塵M TCEP, 5% (v/v) glycerol) supplemented with protease inhibitors, and lysed by sonication. The cleared lysate was loaded on Ni鈥怤TA resin and after several washing steps with buffer A and buffer A containing 25鈥塵M imidazole the protein was eluted with elution buffer containing 300鈥塵M imidazole. After cleavage of the tag by His鈥恡agged TEV protease at 4鈥壜癈 overnight, samples were reloaded on Ni鈥怤TA resin to remove the tag, TEV protease and uncleaved protein. All protein samples were further purified by size-exclusion chromatography on a HiLoad 16/60 Superdex 75 column (GE Healthcare) equilibrated with NMR buffer (20鈥塵M sodium phosphate pH 6.5, 100鈥塵M NaCl, 2鈥塵M DTT).Expression vectorsThe plasmid pCGT7 hnRNP A1 and pCG T7 UP1 have been previously described48. In brief, these expression vectors are under the control of the CMV enhancer/promoter and include an N-terminal T7 epitope tag that corresponds to the first eleven residues of the bacteriophage T7 gene 10 capsid protein. The construct pCG T7 UP1-M9 harbors the M9 import/export sequence fused downstream of both RRM domains of UP1. The UP1-M9 mutants were synthesized by Invitrogen including the corresponding mutations and subcloned into pCG T7 (Supplementary Table聽1). The mutated residues correspond to F148D/F150D (UP1-M9_FD2); F57D/F59D/F148D/F150D (UP1-M9_FD12); F17D/F57D/F59D/F108D/F148D/F150D (UP1-M9_FD12a); F17A/F57A/F59A/F108A/F148A/F150A (UP1-M9_FA12b).Indirect immunofluorescenceCells were fixed 24鈥塰 after transfection with 4% p-formaldehyde in PBS for 15鈥?0鈥塵in at room temperature and permeabilized for 10鈥塵in in 0.2 % Triton X-100. Cells were later incubated with 1:1000 anti-T7 monoclonal antibody (Merck, 69522) for 1鈥塰 at room temperature, washed with PBS, and further incubated with 1:200 fluorescein-conjugated goat anti-mouse IgG (Cappel Laboratories) (1鈥塰 at room temperature). Samples were observed on a Zeiss Axioskop microscope and the images were acquired with a Photometrics CH250 cooled CCD camera using Digital Scientific Smartcapture extensions (IP Lab Spectrum software). The immunofluorescence figures show representative data, and each experiment was reproduced in multiple independent transfections.Analysis of miRNA levels in living cellsHeLa cells were cultured in standard DMEM medium (Life Technologies) supplemented with 10% FBS (Life Technologies). Plasmids were transfected into HeLa cells using Lipofectamine 2000 reagent, as described before49. Mouse monoclonal anti-T7 tag HRP conjugate (1:10000, 69048, RRID鈥擜B10807495, Novagen) was used to detect T7-tagged proteins. Mouse monoclonal anti-tubulin (1:10000, T6199, RRID鈥擜B_477583, Sigma-Aldrich) was used as loading control. miScript qRT-PCR kit (Qiagen) was used to detect miRNAs in total RNA isolated with TRIzol reagent (Life Technologies). Each sample was run in duplicate. To measure the levels of the corresponding miRNAs, values were normalized to 5鈥塖 RNA. For each measurement, three independent experiments were performed.RNA samplesUnlabeled 7-mer, 10-mer, 12-mer, and 17-mer RNA oligonucleotides were purchased from IBA in double-desalted form (G枚ttingen, Germany). All other RNA samples were made by in vitro transcription and purified by denaturing PAGE. RNA samples were heated to 95鈥壜癈 for 3鈥塵in and snap-cooled on ice before use to promote proper folding.EMSAElectrophoretic mobility shift assays (EMSA) were performed with internally labeled RNAs and proteins produced in E. coli. Gel-purified RNAs (50鈥壝椻€?03 c.p.m. (counts per minute), ~20鈥塸mol) were incubated in 15鈥壩糒 reaction mixtures containing the specified amounts of proteins in Roeder D buffer (100鈥塵M KCl, 20% (v/v) glycerol, 0.2鈥塵M EDTA, 100鈥塵M Tris at pH鈥?鈥?.0, 0.5鈥塵M DTT, 0.2鈥塵M PMSF) complemented with 0.5鈥塵M ATP, 20鈥塵M creatine phosphate, and 3.2鈥塵M MgCl2. Reactions mixtures were incubated at 4鈥壜癈 for 1鈥塰 followed by electrophoresis on a 6% (w/v) non-denaturing gel. The radioactive signal was captured with radiographic film or was exposed to a phosphoimaging screen and scanned on a FLA-5100 scanner (Fujifilm).In vitro processing assaysPri-miRNA substrates were obtained by in vitro transcription with [伪-32P]-UTP. Gel-purified substrates (20鈥壝椻€?03 c.p.m. (counts per minute), ~20鈥塸mol) were incubated in 30鈥壩糒 reaction mixtures containing 50% HeLa cell extract in Roeder D buffer, 0.5鈥塵M ATP, 20鈥塵M creatine phosphate, and 3.2鈥塵M MgCl2. Then the reactions were incubated at 37鈥壜癈 for 30鈥塵in. The reactions were stopped with 2鈥壝椻€塙rea Dye and followed by 8% (w/v) denaturing gel electrophoresis. The signal was registered with a radiographic film or by exposure to a phosphoimaging screen and scanning on a FLA-5100 scanner (Fujifilm).RNA pull-downRNA pull-down was performed as previously described49. In brief, protein extracts from HeLa cells were incubated with in vitro-transcribed RNAs, coupled to adipic acid dihydrazide-agarose beads (Sigma). Next, the beads were washed three times with buffer G (20鈥塵M Tris pH 7.5, 135鈥塵M NaCl, 1.5鈥塵M MgCl2, 10% (v/v) glycerol, 1鈥塵M EDTA, 1鈥塵M DTT, and 0.2鈥塵M PMSF). After the final wash, the proteins associated with the beads were resolved and identified by SDS鈥揚AGE and western blotting, respectively. The following antibodies were used: rabbit polyclonal anti-DHX9 (1:1000, 17721-1-AP, RRID鈥擜B_2092506, Protein-Tech); rabbit polyclonal anti-hnRNP A1 antibody (1:1000, PA5-19431, Invitrogen).Footprinting assaysThe assays were performed as described before50. In summary, the RNA substrates were synthesized by in vitro transcription and were 5鈥?labeled with PKA. A formamide ladder and ribonuclease T1 ladder were generated to assign the nucleotide position. Ribonuclease T1 at 0.5鈥塙/渭L was added to 1鈥壩糒 of RNA (50鈥壝椻€?03 c.p.m.) and 7鈥壩糒 of 1鈥壝椻€塻tructure buffer (12鈥塵M Tris-HCl at pH鈥?鈥?.5, 48鈥塵M NaCl, 1.2鈥塵M MgCl2). RNAs substrates were unfolded at 90鈥壜癈 for 1鈥塵in and left at RT for 5鈥塵in to refold. The footprinting reactions were incubated at 37鈥壜癈 for 10鈥塵in. Reactions were run in the presence and absence of the recombinant UP1 protein and ribonuclease T1. Reactions were run on 10% polyacrylamide gel. The signal was registered with a radiographic film or via exposure to a phosphoimaging screen and then scanned on a FLA-5100 scanner (Fujifilm).NMR spectroscopyNMR experiments were recorded at 298鈥塊 on 900, 800, and 600鈥塎Hz Bruker Avance NMR spectrometers, equipped with cryogenic triple resonance gradient probes. NMR spectra were processed by NMRPipe51 and analyzed using Sparky (T. D. Goddard and D. G. Kneller, SPARKY 3, University of California, San Francisco). Samples were measured in NMR buffer with 10% D2O added as lock signal. Backbone resonance assignments of UP1 alone and in complex with RNA were obtained from a uniformly 15N,13C-labeled UP1 (with random fractional deuteration) in the absence and presence of saturating concentrations of RNA. Standard triple resonance experiments HNCA, HNCACB and CBCA(CO)NH52 were recorded at 600鈥塎Hz.15N relaxation experiments were recorded on a 600-MHz spectrometer at 25鈥壜癈. 15N T1 and T1蟻 relaxation times were obtained from pseudo-3D HSQC-based experiments recorded in an interleaved fashion with 12 different relaxation delays (21.6, 86.4, 162, 248.4, 345.6, 432, 518.4, 669.6, 885.6, 1144.8, 1404, and 1782鈥塵s) for T1 and 8 different relaxation delays (5, 7, 10, 15, 20, 25, 30, and 40鈥塵s) for T1蟻. Two delays in each experiment were recorded in duplicates for error estimation. Relaxation rates were extracted by fitting the data to an exponential function using the relaxation module in NMRViewJ53.Paramagnetic relaxation enhancements (PREs) were recorded at 600鈥塎Hz using a sample with concentration of ~200鈥壜礛. UP1(Glu66Cys) was spin-labeled by adding 10-fold excess of 3-(2-Iodoacetamido)-PROXYL to the protein solution (in 100鈥塵M Tris pH 8.0, 100鈥塵M NaCl) and incubating at 4鈥壜癈 overnight. Unreacted spin-label was removed by size-exclusion chromatography. Completion of the spin-labeling reaction and attachment of a single spin-label on the protein was confirmed by LC/MS. Control spin-labeling reactions using wild-type UP1 protein showed that the native cysteines at positions 43 and 175 are inaccessible. PRE data were recorded and analyzed as described previously26,54.Residual dipolar couplings (RDCs) for 1HN-15N bonds were measured on a ~200鈥壜礛 sample of the UP1/12-mer RNA complex aligned in Pf1 phage (~10鈥塵g/mL) using the IPAP experiment55. RDCs were best-fitted to the structure by singular value decomposition (SVD) using PALES56.X-ray crystallographyScreening for crystallization conditions was done using commercial screens (QIAGEN and Hampton Research) at 20鈥壜癈. UP1 (in 20鈥塵M HEPES pH 7.0, 100鈥塵M NaCl, and 1鈥塵M TCEP) was mixed with 12-mer RNA in a 1:1 molar ratio (0.95鈥塵M protein鈥揜NA complex concentration) and incubated on ice for 1鈥塰. Crystallization drops were set-up by mixing 100鈥塶L of complex and 100鈥塶L of reservoir using the sitting drop vapor diffusion method. Crystals with hexagonal plate morphology were obtained in 0.2鈥塎 sodium citrate and 20% PEG 3350 after one week. The crystals were cryoprotected in reservoir solution supplemented with 20% (v/v) ethylene glycol and flash frozen in liquid nitrogen. X-ray diffraction data were collected at the European Synchrotron Radiation Facility (ESRF). All diffraction images were processed by XDS57. The CCP4 package58 was used for all subsequent data analysis. The structure of the UP1/12-mer RNA complex was solved by molecular replacement with the program Phaser59 using the coordinates of UP1 bound to modified telomeric DNA (PDB code: 1U1R)60 as search model. Model building was performed manually in Coot61 and refinement was done using Phenix62.Restrained molecular dynamics simulationsMD simulations were performed using GROMACS63 and PLUMED64. The systems were prepared using the AMBER99SB-ILDN force field65 in tandem with the TIP3P water model66. After solvation in a cubic box (100鈥壝?side including 21,300 water molecules and 250鈥壝?side including 531,000 water molecules for the UP1/12-mer RNA and the UP1/pri-mir-18a complex, respectively), a short initial energy minimization of 100 steps was performed to resolve steric clashes between solvent atoms and the complex. Structures were refined using simulated annealing simulations. Bonds were constrained using LINCS67 and water was kept rigid using SETTLE.The temperature was adjusted with a period of 20鈥塸s between 300鈥塊 and 100鈥塊. The SAXS data set of the UP1/12-mer RNA complex encompassed data points with scattering wavenumbers between 0.025鈥壝?sup>鈭? and 0.685鈥壝?sup>鈭?. Intensities up to 0.5鈥壝?sup>鈭? were fitted with a polynomial of 16th degree. From this fit, 43 intensities were calculated for wavenumbers between 0.03 and 0.45鈥壝?sup>鈭? in steps of 0.01鈥壝?sup>鈭?. These representative intensities were utilized as restraints. For the UP1/pri-mir-18a complex representative intensities between 0.08鈥壝?sup>鈭? and 0.22鈥壝?sup>鈭? were calculated and used as restraints as described above. An approximate scaling factor relating calculated and measured SAXS intensities was estimated by the average ratio between the experimental SAXS intensities as taken from the fit and the first round of calculated SAXS intensities. SAXS intensities were calculated using the Debye formula and standard atomistic structure factors corrected for the effect of the solvent68.Metainference69 with a Gaussian likelihood per data point on the representative SAXS intensities was applied every 10th step. Note, that metainference in the approximation of the absence of dynamics as used here (without replicas and a standard error of the mean set to zero) is equivalent to the Inferential Structure Determination approach70. The attributes of the uncertainty parameter were initially set to large values to allow a slow increase of the restrain force. Additionally, an additional scaling factor between the experimental data was sampled using a flat prior between 0.9 and 1.1.The protein鈥揜NA interface was restrained by harmonic upper-wall potentials centered聽at 3.5鈥壝?applied on the distances between the centers of the respective rings (Phe17鈥揂4, Phe59鈥揋5, Phe108鈥揂9, and Phe150鈥揋10) with energy constants of 1000鈥塳J/mol. Furthermore, two crystallographic salt bridges (Arg75鈥揂sp155 and Arg88鈥揂sp157, distance between side-chain nitrogens and oxygens restrained to be below 5鈥壝? between the two RRM domains were restrained with similar potentials centered at a distance of 4鈥壝?applied on the distances between their charged groups. Secondary structures identified by STRIDE71 from the crystal structure were restrained using an upper-wall potential on the rmsd of the backbone atoms of residues involved with an energy constant of 10000鈥塳J/mol, centered at 0鈥壝?For the UP1/12-mer complex RDCs restraints were applied using the 胃-method72. To take into account for the multiple possible alignments of the molecule with the phage, RDCs were calculated as averages over two replicas and a linear restraint with a slope of 鈭?0,000鈥塳J/mol was applied on the correlation between the average and experimental RDCs. Each replica was independently restrained with all the formerly introduced restraints.After 300 preliminary annealing cycles (for a nominal simulation time of 6鈥塶s per refinement), the refined structure was chosen as the one with the lowest metainference energy among those sampled at 300鈥塊 in the latest 30 cycles. The quality of the structure was then further assessed using PROCHECK73 and the SAXS/SANS profiles and RDCs were independently confirmed using CRYSOL74 and PALES56, respectively.Isothermal titration calorimetryITC (isothermal titration calorimetry) measurements were carried out at 25鈥壜癈 using an iTC200 calorimeter (GE Healthcare). Both protein and ligand were exchanged into NMR buffer without DTT. Protein concentrations in the range of 300鈥?000鈥壜礛, depending on the affinity of the interaction, were injected into the sample cell containing RNA with a concentration of 20鈥?00鈥壜礛. Titrations consisted of 20 injections of 2鈥壜礚 or 26 injections of 1.5鈥壜礚 with a 3-min spacing between each injection. After correction for heat of dilution, data were fitted to a one-site binding model using the Microcal Origin 7.0 software. Each measurement was repeated at least three times.Static light scatteringStatic light scattering (SLS) experiments were performed on a S75 10/300 size-exclusion column (GE Healthcare) connected to a Viscotek Tetra Detector Array (TDA) instrument equipped with refractive index (RI), light scattering, viscosity and photo diode array (PDA) detectors (Malvern Instruments). Sample volume of 100鈥壜礚 was injected onto the column pre-equilibrated with NMR buffer and the flow rate was set to 0.5鈥塵L/min. Calibration was done with BSA (bovine serum albumin) at a concentration of 4鈥?鈥塵g/mL. Data were analyzed by the OmniSEC software using refractive index increment (dn/dc) values of 0.185鈥塵L/g and 0.17鈥塵L/g for protein and RNA samples, respectively75.Small angle X-ray/neutron scatteringSAXS data of UP1, pri-mir-18a, and UP1/pri-mir-18a complex were recorded at the X33 beamline of the European Molecular Biology Laboratory (EMBL) at Deutsches Elektronen Synchrotron (DESY, Hamburg) at 15鈥壜癈. The scattering curves were measured with 120-s exposure times (8 frames, 15鈥塻 each) for concentrations in the range of 1鈥?0鈥塵g/mL. The scattering intensity was measured covering the momentum transfer range 0.007鈥?lt;鈥?i>q鈥?lt;鈥?.63鈥壝?sup>鈭?. Individual frames collected during the exposure time were compared to check for radiation damage before averaging. Scattering of buffer measured before and after each sample was averaged and subtracted from the scattering of the sample. SAXS data of the UP1/12-mer complex were recorded on a BioSAXS-1000 instrument (Rigaku) equipped with a Pilatus detector using 20鈥壜礚 of sample in a capillary tube. All SAXS data were processed using the ATSAS software package76. The radius of gyration (Rg) and the maximum dimension (Dmax) values were obtained from the GNOM program, which evaluates the pair-distance distribution function, P(r)77. For ab initio modeling of the protein鈥揜NA complex, three scattering profiles corresponding to individual components and the complex were used as input for the multiphase modeling program MONSA78. Theoretical scattering curves were calculated using CRYSOL74.SANS data were recorded at the large dynamic range diffractometer D22 at the Institut Laue-Langevin (ILL) Grenoble, France, using a neutron wavelength of 6鈥壝?and a detector-collimator set-up of 2鈥塵/2鈥塵. The scattering intensity was measured covering the momentum transfer range 0.02鈥?lt;鈥?i>q鈥?lt;鈥?.35鈥壝?sup>鈭?. Sample and buffer volumes were 200鈥壩糒 and exposed for 60鈥塵in. The 2D detector signals of the samples were corrected for detector efficiency, empty cell scattering, directly calibrated (against water), and azimuthally averaged using ILL in-house software (Gosh, R. E., Egelhaaf, S. U. Rennie, A. R. A Computing Guide for Small Angle Scattering Experiments. Technical Report ILL06GH05T, 2006, Institut Laue-Langevin, Grenoble, France). The final 1D scattering curves were further analyzed using the ATSAS software.SHAPE analysisPri-mir-17-18a-19a (or pri-mir-18a) in complex with purified hnRNP A1 or UP1 proteins (125鈥塶M) were assembled in folding buffer (100鈥塵M HEPES pH 8, 6鈥塵M MgCl2) using 160鈥塶M RNA prior to treatment with N-methylisatoic anhydride NMIA (Invitrogen), as the modifying agent, as recently described79. RNA was phenol-extracted and ethanol-precipitated and then subjected to primer extension analysis. For primer extension, equal amounts of NMIA-treated and untreated RNAs (10鈥壜礚) were incubated with 0.5鈥壜礚 of antisense 5鈥?end 32P-labeled primers (SHAPE1-3; Supplementary Table聽3), which were used for the analysis of regions comprising miR-17, miR-18a, and miR-19a, respectively. Primer extension and SHAPE data processing was conducted, as previously described80.Data availabilityAtomic coordinates and structure files for the UP1/12-mer RNA crystal structure have been deposited in the Protein Data Bank (http://www.pdb.org/) with accession code 6DCL. Coordinates for the structural models of the 1:1 complexes of UP1/12-mer and UP1/pri-mir-18a are available upon request. Uncropped gels are presented in Supplementary Fig.聽7. Other data are available from the corresponding authors upon reasonable request. References1.Ebert, M. S. Sharp, P. A. Roles for microRNAs in conferring robustness to biological processes. Cell 149, 515鈥?24 (2012).Article聽 PubMed聽 PubMed Central聽 CAS聽Google Scholar聽 2.Ha, M. Kim, V. N. Regulation of microRNA biogenesis. Nat. Rev. Mol. Cell Biol. 15, 509鈥?24 (2014).Article聽 PubMed聽 CAS聽Google Scholar聽 3.Rodriguez, A., Griffiths-Jones, S., Ashurst, J. L. Bradley, A. Identification of mammalian microRNA host genes and transcription units. Genome Res. 14, 1902鈥?910 (2004).Article聽 PubMed聽 PubMed Central聽 CAS聽Google Scholar聽 4.Concepcion, C. P., Bonetti, C. Ventura, A. The microRNA-17-92 family of microRNA clusters in development and disease. Cancer J. 18, 262鈥?67 (2012).Article聽 PubMed聽 PubMed Central聽 CAS聽Google Scholar聽 5.He, L. et al. A microRNA polycistron as a potential human oncogene. Nature 435, 828鈥?33 (2005).ADS聽 Article聽 PubMed聽 PubMed Central聽 CAS聽Google Scholar聽 6.Krol, J., Loedige, I. Filipowicz, W. The widespread regulation of microRNA biogenesis, function and decay. Nat. Rev. Genet. 11, 597鈥?10 (2010).Article聽 PubMed聽 CAS聽Google Scholar聽 7.Michlewski, G., Guil, S., Semple, C. A. Caceres, J. F. Posttranscriptional regulation of miRNAs harboring conserved terminal loops. Mol. Cell 32, 383鈥?93 (2008).Article聽 PubMed聽 PubMed Central聽 CAS聽Google Scholar聽 8.Guil, S. Caceres, J. F. The multifunctional RNA-binding protein hnRNP A1 is required for processing of miR-18a. Nat. Struct. Mol. Biol. 14, 591鈥?96 (2007).Article聽 PubMed聽 CAS聽Google Scholar聽 9.Mayeda, A. Krainer, A. R. Regulation of alternative pre-mRNA splicing by hnRNP A1 and splicing factor SF2. Cell 68, 365鈥?75 (1992).Article聽 PubMed聽 CAS聽Google Scholar聽 10.Tavanez, J. P., Madl, T., Kooshapur, H., Sattler, M. Valcarcel, J. hnRNP A1 proofreads 3鈥?splice site recognition by U2AF. Mol. Cell 45, 314鈥?29 (2012).Article聽 PubMed聽 CAS聽Google Scholar聽 11.Izaurralde, E. et al. A role for the M9 transport signal of hnRNP A1 in mRNA nuclear export. J. Cell Biol. 137, 27鈥?5 (1997).Article聽 PubMed聽 PubMed Central聽 CAS聽Google Scholar聽 12.Bonnal, S. et al. Heterogeneous nuclear ribonucleoprotein A1 is a novel internal ribosome entry site trans-acting factor that modulates alternative initiation of translation of the fibroblast growth factor 2 mRNA. J. Biol. Chem. 280, 4144鈥?153 (2005).Article聽 PubMed聽 CAS聽Google Scholar聽 13.Hamilton, B. J., Burns, C. M., Nichols, R. C. Rigby, W. F. Modulation of AUUUA response element binding by heterogeneous nuclear ribonucleoprotein A1 in human T lymphocytes. The roles of cytoplasmic location, transcription, and phosphorylation. J. Biol. Chem. 272, 28732鈥?8741 (1997).Article聽 PubMed聽 CAS聽Google Scholar聽 14.LaBranche, H. et al. Telomere elongation by hnRNP A1 and a derivative that interacts with telomeric repeats and telomerase. Nat. Genet. 19, 199鈥?02 (1998).Article聽 PubMed聽 CAS聽Google Scholar聽 15.Zhang, Q. S., Manche, L., Xu, R. M. Krainer, A. R. hnRNP A1 associates with telomere ends and stimulates telomerase activity. RNA 12, 1116鈥?128 (2006).Article聽 PubMed聽 PubMed Central聽 CAS聽Google Scholar聽 16.Siomi, H. Dreyfuss, G. A nuclear localization domain in the hnRNP A1 protein. J. Cell Biol. 129, 551鈥?60 (1995).Article聽 PubMed聽 CAS聽Google Scholar聽 17.Mayeda, A., Munroe, S. H., Caceres, J. F. Krainer, A. R. Function of conserved domains of hnRNP A1 and other hnRNP A/B proteins. EMBO J. 13, 5483鈥?495 (1994).PubMed聽 PubMed Central聽 CAS聽 Article聽Google Scholar聽 18.Burd, C. G. Dreyfuss, G. RNA binding specificity of hnRNP A1: significance of hnRNP A1 high-affinity binding sites in pre-mRNA splicing. EMBO J. 13, 1197鈥?204 (1994).PubMed聽 PubMed Central聽 CAS聽 Article聽Google Scholar聽 19.Barraud, P. Allain, F. H. Solution structure of the two RNA recognition motifs of hnRNP A1 using segmental isotope labeling: how the relative orientation between RRMs influences the nucleic acid binding topology. J. Biomol. NMR 55, 119鈥?38 (2013).Article聽 PubMed聽 CAS聽Google Scholar聽 20.Ding, J. et al. Crystal structure of the two-RRM domain of hnRNP A1 (UP1) complexed with single-stranded telomeric DNA. Genes. Dev. 13, 1102鈥?115 (1999).Article聽 PubMed聽 PubMed Central聽 CAS聽Google Scholar聽 21.Beusch, I., Barraud, P., Moursy, A., Clery, A., Allain, F. H. Tandem hnRNP A1 RNA recognition motifs act in concert to repress the splicing of survival motor neuron exon 7. eLife 6, e25736 (2017).22.Morgan, C. E. et al. The first crystal structure of the UP1 domain of hnRNP A1 bound to RNA reveals a new look for an old RNA binding protein. J. Mol. Biol. 427, 3241鈥?257 (2015).Article聽 PubMed聽 PubMed Central聽 CAS聽Google Scholar聽 23.Mayeda, A., Munroe, S. H., Xu, R. M. Krainer, A. R. Distinct functions of the closely related tandem RNA-recognition motifs of hnRNP A1. RNA 4, 1111鈥?123 (1998).Article聽 PubMed聽 PubMed Central聽 CAS聽Google Scholar聽 24.Gobl, C., Madl, T., Simon, B. Sattler, M. NMR approaches for structural analysis of multidomain proteins and complexes in solution. Prog. NMR Spectrosc. 80, 26鈥?3 (2014).Article聽 CAS聽Google Scholar聽 25.Hennig, J. et al. Structural analysis of protein-RNA complexes in solution using NMR paramagnetic relaxation enhancements. Methods Enzymol. 558, 333鈥?62 (2015).Article聽 PubMed聽 CAS聽Google Scholar聽 26.Simon, B., Madl, T., Mackereth, C. D., Nilges, M. Sattler, M. An efficient protocol for NMR-spectroscopy-based structure determination of protein complexes in solution. Angew. Chem. Int. Ed. Engl. 49, 1967鈥?970 (2010).Article聽 PubMed聽 CAS聽Google Scholar聽 27.Parisien, M. Major, F. The MC-Fold and MC-Sym pipeline infers RNA structure from sequence data. Nature 452, 51鈥?5 (2008).ADS聽 Article聽 PubMed聽 CAS聽Google Scholar聽 28.Dallmann, A. Sattler, M. Detection of hydrogen bonds in dynamic regions of RNA by NMR spectroscopy. Curent Protoc. Nucleic Acid Chem. 59, 7 22 21鈥?9 (2014). Google Scholar聽 29.Dallmann, A. et al. Efficient detection of hydrogen bonds in dynamic regions of RNA by sensitivity-optimized NMR pulse sequences. Angew. Chem. Int. Ed. Engl. 52, 10487鈥?0490 (2013).Article聽 PubMed聽 CAS聽Google Scholar聽 30.Wilkinson, K. A., Merino, E. J. Weeks, K. M. Selective 2鈥?hydroxyl acylation analyzed by primer extension (SHAPE): quantitative RNA structure analysis at single nucleotide resolution. Nat. Protoc. 1, 1610鈥?616 (2006).Article聽 PubMed聽 CAS聽Google Scholar聽 31.Reuter, J. S. Mathews, D. H. RNAstructure: software for RNA secondary structure prediction and analysis. BMC Bioinformatics. 11, 129 (2010).Article聽 PubMed聽 PubMed Central聽 CAS聽Google Scholar聽 32.Chakraborty, S. Krishnan, Y. A structural map of oncomiR-1 at single-nucleotide resolution. Nucleic Acids Res. 45, 9694鈥?705 (2017).Article聽 PubMed聽 PubMed Central聽 CAS聽Google Scholar聽 33.Herrick, G., Delius, H. Alberts, B. Single-stranded DNA structure and DNA polymerase activity in the presence of nucleic acid helix-unwinding proteins from calf thymus. J. Biol. Chem. 251, 2142鈥?146 (1976).PubMed聽 CAS聽Google Scholar聽 34.Michlewski, G. Caceres, J. F. Antagonistic role of hnRNP A1 and KSRP in the regulation of let-7a biogenesis. Nat. Struct. Mol. Biol. 17, 1011鈥?018 (2010).Article聽 PubMed聽 PubMed Central聽 CAS聽Google Scholar聽 35.Trabucchi, M. et al. The RNA-binding protein KSRP promotes the biogenesis of a subset of microRNAs. Nature 459, 1010鈥?014 (2009).ADS聽 Article聽 PubMed聽 PubMed Central聽 CAS聽Google Scholar聽 36.Castilla-Llorente, V., Nicastro, G. Ramos, A. Terminal loop-mediated regulation of miRNA biogenesis: selectivity and mechanisms. Biochem. Soc. Trans. 41, 861鈥?65 (2013).Article聽 PubMed聽 PubMed Central聽 CAS聽Google Scholar聽 37.Nam, Y., Chen, C., Gregory, R. I., Chou, J. J. Sliz, P. Molecular basis for interaction of let-7 microRNAs with Lin28. Cell 147, 1080鈥?091 (2011).Article聽 PubMed聽 PubMed Central聽 CAS聽Google Scholar聽 38.Chen, Y. et al. Rbfox proteins regulate microRNA biogenesis by sequence-specific binding to their precursors and target downstream Dicer. Nucleic Acids Res. 44, 4381鈥?395 (2016).Article聽 PubMed聽 PubMed Central聽 CAS聽Google Scholar聽 39.Auyeung, V. C., Ulitsky, I., McGeary, S. E. Bartel, D. P. Beyond secondary structure: primary-sequence determinants license pri-miRNA hairpins for processing. Cell 152, 844鈥?58 (2013).Article聽 PubMed聽 PubMed Central聽 CAS聽Google Scholar聽 40.Nguyen, T. A. et al. Functional anatomy of the human microprocessor. Cell 161, 1374鈥?387 (2015).Article聽 PubMed聽 CAS聽Google Scholar聽 41.Partin, A. C. et al. Heme enables proper positioning of Drosha and DGCR8 on primary microRNAs. Nat. Commun. 8, 1737 (2017).ADS聽 Article聽 PubMed聽 PubMed Central聽 CAS聽Google Scholar聽 42.Roden, C. et al. Novel determinants of mammalian primary microRNA processing revealed by systematic evaluation of hairpin-containing transcripts and human genetic variation. Genome Res. 27, 374鈥?84 (2017).Article聽 PubMed聽 PubMed Central聽 CAS聽Google Scholar聽 43.Fang, W. Bartel, D. P. The menu of features that define primary microRNAs and enable de novo design of MicroRNA genes. Mol. Cell 60, 131鈥?45 (2015).Article聽 PubMed聽 PubMed Central聽 CAS聽Google Scholar聽 44.Kim, K., Nguyen, T. D., Li, S. Nguyen, T. A. SRSF3 recruits DROSHA to the basal junction of primary microRNAs. RNA 24,聽https://doi.org/10.1261/rna.065862.118聽(2018).45.Chaulk, S. G., Xu, Z., Glover, M. J. Fahlman, R. P. MicroRNA miR-92a-1 biogenesis and mRNA targeting is modulated by a tertiary contact within the miR-17~92 microRNA cluster. Nucleic Acids Res. 42, 5234鈥?244 (2014).Article聽 PubMed聽 PubMed Central聽 CAS聽Google Scholar聽 46.Chakraborty, S., Mehtab, S., Patwardhan, A. Krishnan, Y. Pri-miR-17-92a transcript folds into a tertiary structure and autoregulates its processing. RNA 18, 1014鈥?028 (2012).Article聽 PubMed聽 PubMed Central聽 CAS聽Google Scholar聽 47.Du, P., Wang, L., Sliz, P. Gregory, R. I. A biogenesis step upstream of Microprocessor controls miR-17 approximately 92 expression. Cell 162, 885鈥?99 (2015).Article聽 PubMed聽 PubMed Central聽 CAS聽Google Scholar聽 48.Guil, S., Long, J. C. Caceres, J. F. hnRNP A1 relocalization to the stress granules reflects a role in the stress response. Mol. Cell. Biol. 26, 5744鈥?758 (2006).Article聽 PubMed聽 PubMed Central聽 CAS聽Google Scholar聽 49.Choudhury, N. R. et al. Trim25 is an RNA-specific activator of Lin28a/TuT4-mediated uridylation. Cell Rep. 9, 1265鈥?272 (2014).Article聽 PubMed聽 PubMed Central聽 CAS聽Google Scholar聽 50.Nowak, J. S. et al. Lin28a uses distinct mechanisms of binding to RNA and affects miRNA levels positively and negatively. RNA 23, 317鈥?32 (2017).Article聽 PubMed聽 PubMed Central聽 CAS聽Google Scholar聽 51.Delaglio, F. et al. NMRPipe: a multidimensional spectral processing system based on UNIX pipes. J. Biomol. NMR 6, 277鈥?93 (1995).Article聽 PubMed聽 CAS聽Google Scholar聽 52.Sattler, M., J, S. C, G. Heteronuclear multidimensional NMR experiments for the structure determination of proteins in solution employing pulsed field gradients. Prog. NMR Spectrosc. 34, 93鈥?58 (1999).Article聽 CAS聽Google Scholar聽 53.Johnson, B. A. Blevins, R. A. NMRView: a computer program for the visualization and analysis of NMR data. J. Biomol. NMR 4, 603鈥?14 (1994).Article聽 PubMed聽 CAS聽Google Scholar聽 54.Mackereth, C. D. et al. Multi-domain conformational selection underlies pre-mRNA splicing regulation by U2AF. Nature 475, 408鈥?11 (2011).Article聽 PubMed聽 CAS聽Google Scholar聽 55.Cordier, F., Dingley, A. J. Grzesiek, S. A doublet-separated sensitivity-enhanced HSQC for the determination of scalar and dipolar one-bond J-couplings. J. Biomol. NMR 13, 175鈥?80 (1999).Article聽 PubMed聽 CAS聽Google Scholar聽 56.Zweckstetter, M. NMR: prediction of molecular alignment from structure using the PALES software. Nat. Protoc. 3, 679鈥?90 (2008).Article聽 PubMed聽 CAS聽Google Scholar聽 57.Kabsch, W. XDS. Acta Cryst. D Biol. Cryst. 66, 125鈥?32 (2010).Article聽 CAS聽Google Scholar聽 58.Winn, M. D. et al. Overview of the CCP4 suite and current developments. Acta Cryst. D Biol. Cryst. 67, 235鈥?42 (2011).Article聽 CAS聽Google Scholar聽 59.McCoy, A. J. et al. Phaser crystallographic software. J. Appl. Cryst. 40, 658鈥?74 (2007).Article聽 CAS聽Google Scholar聽 60.Myers, J. C. Shamoo, Y. Human UP1 as a model for understanding purine recognition in the family of proteins containing the RNA recognition motif (RRM). J. Mol. Biol. 342, 743鈥?56 (2004).Article聽 PubMed聽 CAS聽Google Scholar聽 61.Emsley, P., Lohkamp, B., Scott, W. G. Cowtan, K. Features and development of Coot. Acta Cryst. D Biol. Cryst. 66, 486鈥?01 (2010).Article聽 CAS聽Google Scholar聽 62.Adams, P. D. et al. PHENIX: a comprehensive Python-based system for macromolecular structure solution. Acta Cryst. D Biol. Cryst. 66, 213鈥?21 (2010).Article聽 CAS聽Google Scholar聽 63.Abraham, M. J. et al. GROMACS: High performance molecular simulations through multi-level parallelism from laptops to supercomputers. SoftwareX 1-2, 19鈥?5 (2015).ADS聽 Article聽Google Scholar聽 64.Tribello, G. A., Bonomi, M., Branduardi, D., Camilloni, C. Bussi, G. PLUMED 2: New feathers for an old bird. Comput. Phys. Commun. 185, 604鈥?13 (2014).ADS聽 Article聽 CAS聽Google Scholar聽 65.Lindorff-Larsen, K. et al. Improved side-chain torsion potentials for the Amber ff99SB protein force field. Proteins 78, 1950鈥?958 (2010).PubMed聽 PubMed Central聽 CAS聽Google Scholar聽 66.Jorgensen, W. L. Quantum and statistical mechanical studies of liquids. 10. Transferable intermolecular potential functions for water, alcohols, and ethers. Application to liquid water. J. Am. Chem. Soc. 103, 335鈥?40 (1981).Article聽 CAS聽Google Scholar聽 67.Hess, B., Bekker, H., Berendsen, H. J. Fraaije, J. G. LINCS: a linear constraint solver for molecular simulations. J. Comput. Chem. 18, 1463鈥?472 (1997).Article聽 CAS聽Google Scholar聽 68.Fraser, R. D. B., MacRae, T. P. Suzuki, E. An improved method for calculating the contribution of solvent to the X-ray diffraction pattern of biological molecules. J. Appl. Cryst. 11, 693鈥?94 (1978).Article聽 CAS聽Google Scholar聽 69.Bonomi, M., Camilloni, C., Cavalli, A. Vendruscolo, M. Metainference: a Bayesian inference method for heterogeneous systems. Sci. Adv. 2, e1501177 (2016).ADS聽 Article聽 PubMed聽 PubMed Central聽Google Scholar聽 70.Rieping, W., Habeck, M. Nilges, M. Inferential structure determination. Science 309, 303鈥?06 (2005).ADS聽 Article聽 PubMed聽 CAS聽Google Scholar聽 71.Frishman, D. Argos, P. Knowledge-based protein secondary structure assignment. Proteins 23, 566鈥?79 (1995).Article聽 PubMed聽 CAS聽Google Scholar聽 72.Camilloni, C. Vendruscolo, M. A tensor-free method for the structural and dynamical refinement of proteins using residual dipolar couplings. J. Phys. Chem. B 119, 653鈥?61 (2015).Article聽 PubMed聽 CAS聽Google Scholar聽 73.Laskowski, R. A., Macarthur, M. W., Moss, D. S. Thornton, J. M. {PROCHECK}: a program to check the stereochemical quality of protein structures. J. Appl. Cryst. 26, 283鈥?91 (1993).Article聽 CAS聽Google Scholar聽 74.Svergun, D., Barberato, C. Koch, M. H. J. CRYSOL鈥?a program to evaluate X-ray solution scattering of biological macromolecules from atomic coordinates. J. Appl. Cryst. 28, 768鈥?73 (1995).Article聽 CAS聽Google Scholar聽 75.Rambo, R. P. Doudna, J. A. Assembly of an active group II intron-maturase complex by protein dimerization. Biochemistry 43, 6486鈥?497 (2004).Article聽 PubMed聽 CAS聽Google Scholar聽 76.Petoukhov, M. V. et al. New developments in the ATSAS program package for small-angle scattering data analysis. J. Appl. Cryst. 45, 342鈥?50 (2012).Article聽 CAS聽Google Scholar聽 77.Svergun, D. I. Determination of the regularization parameter in indirect-transform methods using perceptual criteria. J. Appl. Cryst. 25, 495鈥?03 (1992).Article聽 CAS聽Google Scholar聽 78.Svergun, D. I. Restoring low resolution structure of biological macromolecules from solution scattering using simulated annealing. Biophys. J. 76, 2879鈥?886 (1999).ADS聽 Article聽 PubMed聽 PubMed Central聽 CAS聽Google Scholar聽 79.Fernandez, N. et al. Genetic variation and RNA structure regulate microRNA biogenesis. Nat. Commun. 8, 15114 (2017).ADS聽 Article聽 PubMed聽 PubMed Central聽Google Scholar聽 80.Fernandez, N. et al. Structural basis for the biological relevance of the invariant apical stem in IRES-mediated translation. Nucleic Acids Res. 39, 8572鈥?585 (2011).Article聽 PubMed聽 PubMed Central聽 CAS聽Google Scholar聽 Download referencesAcknowledgementsWe thank Marianne Keith (MRC HGU, University of Edinburgh) for Immunofluorescence data, Robert Janowski for collecting the X-ray diffraction data, and Ralf Stehle for help with in-house SAXS measurements. We are grateful to the beamline staff at DESY (Hamburg) and Anne Martel at ILL (Grenoble) for support with SAXS and SANS measurements, respectively. We acknowledge the use of the X-ray crystallography platform at the Helmholtz Zentrum M眉nchen, NMR access at the Bavarian NMR Center and SAXS measurements at the facility of the SFB1035 at Department Chemie, Technical University of Munich. This work was supported by the Deutsche Forschungsgemeinschaft through grants SFB1035 and GRK1721 (to M.S.), MRC core funding and the Wellcome Trust (Grant 095518/Z/11/Z to J.F.C.), Medical Research Council Career Development Award (G10000564 to G.M.). Wellcome Trust Center Core Grants (077707 and 092076 to G.M.). M.M., A.J., and C.C. acknowledge the support of the Technical University of Munich鈥擨nstitute for Advanced Study, funded by the German Excellence Initiative and the European Union Seventh Framework Programme under grant agreement n掳 291763.Author informationAffiliationsInstitute of Structural Biology, Helmholtz Zentrum M眉nchen, Neuherberg, 85764, GermanyHamed Kooshapur,聽Alisha N. Jones,聽Andre Dallmann聽 聽Michael SattlerCenter for Integrated Protein Science Munich at Department Chemie, Technical University of Munich, Garching, 85747, GermanyHamed Kooshapur,聽Max M眉hlbauer,聽Alexander Jussupow,聽Alisha N. Jones,聽Andre Dallmann,聽Carlo Camilloni聽 聽Michael SattlerWellcome Trust Centre for Cell Biology, University of Edinburgh, Michael Swann Building, Edinburgh, EH9 3BF, UKNila Roy Choudhury聽 聽Gracjan MichlewskiStructural and Computational Biology Unit, European Molecular Biology Laboratory, Heidelberg, 69117, GermanyBernd SimonInstitute for Advanced Study, Technical University of Munich, Garching, 85747, GermanyMax M眉hlbauer,聽Alexander Jussupow聽 聽Carlo CamilloniMRC Human Genetics Unit, Institute of Genetics and Molecular Medicine, University of Edinburgh, Western General Hospital, Edinburgh, EH4 2XU, UKNoemi Fernandez,聽Gracjan Michlewski聽 聽Javier F. CaceresInstitut Laue-Langevin, Grenoble, 38044, FranceFrank GabelUniv. Grenoble Alpes, CEA, CNRS, IBS, Grenoble, 38042, FranceFrank GabelDivision of Infection and Pathway Medicine, University of Edinburgh, The Chancellor鈥檚 Building, 49 Little France Crescent, Edinburgh, EH16 4SB, UKGracjan MichlewskiZJU-UoE Institute, Zhejiang University, 718 East Haizhou Road, Haining, Zhejiang, 314400, ChinaGracjan MichlewskiAuthorsHamed KooshapurView author publicationsYou can also search for this author in PubMed聽Google ScholarNila Roy ChoudhuryView author publicationsYou can also search for this author in PubMed聽Google ScholarBernd SimonView author publicationsYou can also search for this author in PubMed聽Google ScholarMax M眉hlbauerView author publicationsYou can also search for this author in PubMed聽Google ScholarAlexander JussupowView author publicationsYou can also search for this author in PubMed聽Google ScholarNoemi FernandezView author publicationsYou can also search for this author in PubMed聽Google ScholarAlisha N. JonesView author publicationsYou can also search for this author in PubMed聽Google ScholarAndre DallmannView author publicationsYou can also search for this author in PubMed聽Google ScholarFrank GabelView author publicationsYou can also search for this author in PubMed聽Google ScholarCarlo CamilloniView author publicationsYou can also search for this author in PubMed聽Google ScholarGracjan MichlewskiView author publicationsYou can also search for this author in PubMed聽Google ScholarJavier F. CaceresView author publicationsYou can also search for this author in PubMed聽Google ScholarMichael SattlerView author publicationsYou can also search for this author in PubMed聽Google ScholarContributionsH.K. performed molecular biology, protein/RNA preparation, ITC, NMR, crystallographic structure determination, RNA modeling, and SAXS analysis; N.R.C. and G.M. performed footprinting and miRNA processing assays both in vitro and in cells in culture, N.F. performed SHAPE experiments, N.F. and A.N.J. analyzed SHAPE data, B.S., M.M., A.J., and C.C. performed structural modeling and restrained molecular dynamics, F.G. analyzed SAXS and SANS data, A.D. performed NMR of RNA samples. H.K., G.M., J.F.C., and M.S. conceived and designed the project. H.K., G.M., J.F.C., and M.S. wrote the paper. All authors discussed the results and commented on the manuscript.Corresponding authorsCorrespondence to Gracjan Michlewski or Javier F. Caceres or Michael Sattler.Ethics declarations Competing interests The authors declare no competing interests. Additional informationPublisher\'s note: Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.Electronic supplementary material CommentsBy submitting a comment you agree to abide by our Terms and Community Guidelines. If you find something abusive or that does not comply with our terms or guidelines please flag it as inappropriate. Sign up for the Nature Briefing newsletter 鈥?what matters in science, free to your inbox daily.