Skip to main content

Thank you for visiting nature.com. You are using a browser version with limited support for CSS. To obtain the best experience, we recommend you use a more up to date browser (or turn off compatibility mode in Internet Explorer). In the meantime, to ensure continued support, we are displaying the site without styles and JavaScript.

  • Article
  • Published:

Quantitative analysis of RNA-protein interactions on a massively parallel array reveals biophysical and evolutionary landscapes

Nature Biotechnology volume 32, pages 562–568 (2014)Cite this article

Subjects

Abstract

RNA-protein interactions drive fundamental biological processes and are targets for molecular engineering, yet quantitative and comprehensive understanding of the sequence determinants of affinity remains limited. Here we repurpose a high-throughput sequencing instrument to quantitatively measure binding and dissociation of a fluorescently labeled protein to >107 RNA targets generated on a flow cell surface by in situ transcription and intermolecular tethering of RNA to DNA. Studying the MS2 coat protein, we decompose the binding energy contributions from primary and secondary RNA structure, and observe that differences in affinity are often driven by sequence-specific changes in both association and dissociation rates. By analyzing the biophysical constraints and modeling mutational paths describing the molecular evolution of MS2 from low- to high-affinity hairpins, we quantify widespread molecular epistasis and a long-hypothesized, structure-dependent preference for G:U base pairs over C:A intermediates in evolutionary trajectories. Our results suggest that quantitative analysis of RNA on a massively parallel array (RNA-MaP) provides generalizable insight into the biophysical basis and evolutionary consequences of sequence-function relationships.

This is a preview of subscription content, access via your institution

Access options

Buy this article

  • Purchase on SpringerLink
  • Instant access to full article PDF
Buy now

Prices may be subject to local taxes which are calculated during checkout

Figure 1: A massively parallel RNA array for quantitative, high-throughput biochemistry.
Figure 2: A quantitative map of MS2 binding across RNA sequence variants.
Figure 3: Decomposition of primary and secondary RNA structure determinants of binding affinity.
Figure 4: Sequence-specific contributions of association and dissociation rates to binding affinity.
Figure 5: Evolutionary landscapes are highly constrained by biophysical requirements.

Similar content being viewed by others

Accession codes

Primary accessions

Sequence Read Archive

Referenced accessions

Protein Data Bank

References

  1. Keene, J.D. RNA regulons: coordination of post-transcriptional events. Nat. Rev. Genet. 8, 533–543 (2007).

    Article  CAS  Google Scholar 

  2. Carey, J., Cameron, V., De Haseth, P.L. & Uhlenbeck, O.C. Sequence-specific interaction of R17 coat protein with its ribonucleic acid binding site. Biochemistry 22, 2601–2610 (1983).

    Article  CAS  Google Scholar 

  3. Tsvetanova, N.G., Klass, D.M., Salzman, J. & Brown, P.O. Proteome-wide search reveals unexpected RNA-binding proteins in Saccharomyces cerevisiae . PLoS ONE 5, e12671 (2010).

    Article  Google Scholar 

  4. Scherrer, T., Mittal, N., Janga, S.C. & Gerber, A.P. A screen for RNA-binding proteins in yeast indicates dual functions for many enzymes. PLoS ONE 5, e15499 (2010).

    Article  CAS  Google Scholar 

  5. Butter, F., Scheibe, M., Morl, M. & Mann, M. Unbiased RNA-protein interaction screen by quantitative proteomics. Proc. Natl. Acad. Sci. USA 106, 10626–10631 (2009).

    Article  CAS  Google Scholar 

  6. Castello, A. et al. Insights into RNA biology from an atlas of mammalian mRNA-binding proteins. Cell 149, 1393–1406 (2012).

    Article  CAS  Google Scholar 

  7. Wang, K.C. et al. A long noncoding RNA maintains active chromatin to coordinate homeotic gene expression. Nature 472, 120–124 (2011).

    Article  CAS  Google Scholar 

  8. Tsai, M.C. et al. Long noncoding RNA as modular scaffold of histone modification complexes. Science 329, 689–693 (2010).

    Article  CAS  Google Scholar 

  9. Guttman, M. & Rinn, J.L. Modular regulatory principles of large non-coding RNAs. Nature 482, 339–346 (2012).

    Article  CAS  Google Scholar 

  10. Culler, S.J., Hoff, K.G. & Smolke, C.D. Reprogramming cellular behavior with RNA controllers responsive to endogenous proteins. Science 330, 1251–1255 (2010).

    Article  CAS  Google Scholar 

  11. Ausländer, S., Ausländer, D., Müller, M., Wieland, M. & Fussenegger, M. Programmable single-cell mammalian biocomputers. Nature 487, 123–127 (2012).

    Article  Google Scholar 

  12. SantaLucia, J. & Turner, D.H. Measuring the thermodynamics of RNA secondary structure formation. Biopolymers 44, 309–319 (1997).

    Article  CAS  Google Scholar 

  13. Kertesz, M. et al. Genome-wide measurement of RNA secondary structure in yeast. Nature 467, 103–107 (2010).

    Article  CAS  Google Scholar 

  14. Ding, Y. et al. In vivo genome-wide profiling of RNA secondary structure reveals novel regulatory features. Nature 505, 696–700 (2014).

    Article  CAS  Google Scholar 

  15. Rouskin, S., Zubradt, M., Washietl, S., Kellis, M. & Weissman, J. Genome-wide probing of RNA structure reveals active unfolding of mRNA structures in vivo . Nature 505, 701–705 (2014).

    Article  CAS  Google Scholar 

  16. Ban, N., Nissen, P., Hansen, J., Moore, P.B. & Steitz, T.A. The complete atomic structure of the large ribosomal subunit at 2.4 Å resolution. Science 289, 905–920 (2000).

    Article  CAS  Google Scholar 

  17. Wan, Y., Kertesz, M., Spitale, R.C., Segal, E. & Chang, H.Y. Understanding the transcriptome through RNA structure. Nat. Rev. Genet. 12, 641–655 (2011).

    Article  CAS  Google Scholar 

  18. Martin, L. et al. Systematic reconstruction of RNA functional motifs with high-throughput microfluidics. Nat. Methods 9, 1192–1194 (2012).

    Article  CAS  Google Scholar 

  19. Ray, D. et al. Rapid and systematic analysis of the RNA recognition specificities of RNA-binding proteins. Nat. Biotechnol. 27, 667–670 (2009).

    Article  CAS  Google Scholar 

  20. Ray, D. et al. A compendium of RNA-binding motifs for decoding gene regulation. Nature 499, 172–177 (2013).

    Article  CAS  Google Scholar 

  21. Araya, C.L. et al. A fundamental protein property, thermodynamic stability, revealed solely from large-scale measurements of protein function. Proc. Natl. Acad. Sci. USA 109, 16858–16863 (2012).

    Article  CAS  Google Scholar 

  22. Pitt, J.N. & Ferre-D'Amare, A.R. Rapid construction of empirical RNA fitness landscapes. Science 330, 376–379 (2010).

    Article  CAS  Google Scholar 

  23. Guenther, U.-P. et al. Hidden specificity in an apparently nonspecific RNA-binding protein. Nature 502, 385–388 (2013).

    Article  CAS  Google Scholar 

  24. Matzas, M. et al. High-fidelity gene synthesis by retrieval of sequence-verified DNA identified using high-throughput pyrosequencing. Nat. Biotechnol. 28, 1291–1294 (2010).

    Article  CAS  Google Scholar 

  25. Myllykangas, S., Buenrostro, J.D., Natsoulis, G., Bell, J.M. & Ji, H.P. Efficient targeted resequencing of human germline and cancer genomes by oligonucleotide-selective sequencing. Nat. Biotechnol. 29, 1024–1027 (2011).

    Article  CAS  Google Scholar 

  26. Uemura, S. et al. Real-time tRNA transit on single translating ribosomes at codon resolution. Nature 464, 1012–1017 (2010).

    Article  CAS  Google Scholar 

  27. Nutiu, R. et al. Direct measurement of DNA affinity landscapes on a high-throughput sequencing instrument. Nat. Biotechnol. 29, 659–664 (2011).

    Article  CAS  Google Scholar 

  28. Bentley, D.R. et al. Accurate whole human genome sequencing using reversible terminator chemistry. Nature 456, 53–59 (2008).

    Article  CAS  Google Scholar 

  29. Carey, J., Lowary, P.T. & Uhlenbeck, O.C. Interaction of R17 coat protein with synthetic variants of its ribonucleic acid binding site. Biochemistry 22, 4723–4730 (1983).

    Article  CAS  Google Scholar 

  30. Lim, F. & David, S.P. Mutations that increase the affinity of a translational repressor for RNA. Nucleic Acids Res. 22, 3748–3752 (1994).

    Article  CAS  Google Scholar 

  31. Valegård, K., Murray, J.B., Stockley, P.G., Stonehouse, N.J. & Liljas, L. Crystal structure of an RNA bacteriophage coat protein-operator complex. Nature 371, 623–626 (1994).

    Article  Google Scholar 

  32. Romaniuk, P.J., Lowary, P., Wu, H.N., Stormo, G. & Uhlenbeck, O.C. RNA binding site of R17 coat protein. Biochemistry 26, 1563–1568 (1987).

    Article  CAS  Google Scholar 

  33. Grahn, E. et al. Structural basis of pyrimidine specificity in the MS2 RNA hairpin-coat-protein complex. RNA 7, 1616–1627 (2001).

    CAS  PubMed  PubMed Central  Google Scholar 

  34. Bardwell, V.J. & Wickens, M. Purification of RNA and RNA-protein complexes by an R17 coat protein affinity method. Nucleic Acids Res. 18, 6587–6594 (1990).

    Article  CAS  Google Scholar 

  35. Bertrand, E. et al. Localization of ASH1 mRNA particles in living yeast. Mol. Cell 2, 437–445 (1998).

    Article  CAS  Google Scholar 

  36. Kivioja, T. et al. Counting absolute numbers of molecules using unique molecular identifiers. Nat. Methods 9, 72–74 (2012).

    Article  CAS  Google Scholar 

  37. Greenleaf, W.J., Frieda, K.L., Foster, D.A., Woodside, M.T. & Block, S.M. Direct observation of hierarchical folding in single riboswitch aptamers. Science 319, 630–633 (2008).

    Article  CAS  Google Scholar 

  38. Hobson, D. & Uhlenbeck, O.C. Alanine scanning of MS2 coat protein reveals protein–phosphate contacts involved in thermodynamic hot spots. J. Mol. Biol. 356, 613–624 (2006).

    Article  CAS  Google Scholar 

  39. Varani, G. & McClain, W.H. The G·U wobble base pair. EMBO Rep. 1, 18–23 (2000).

    Article  CAS  Google Scholar 

  40. Valegård, K. et al. The three-dimensional structures of two complexes between recombinant MS2 capsids and RNA operator fragments reveal sequence-specific protein-RNA interactions. J. Mol. Biol. 270, 724–738 (1997).

    Article  Google Scholar 

  41. Breen, M.S., Kemena, C., Vlasov, P.K., Notredame, C. & Kondrashov, F.A. Epistasis as the primary factor in molecular evolution. Nature 490, 535–538 (2012).

    Article  CAS  Google Scholar 

  42. McCandlish, D.M., Rajon, E., Shah, P., Ding, Y. & Plotkin, J.B. The role of epistasis in protein evolution. Nature 497, E1–E2 (2013).

    Article  CAS  Google Scholar 

  43. Weinreich, D.M. Darwinian evolution can follow only very few mutational paths to fitter proteins. Science 312, 111–114 (2006).

    Article  CAS  Google Scholar 

  44. Bridgham, J.T., Ortlund, E.A. & Thornton, J.W. An epistatic ratchet constrains the direction of glucocorticoid receptor evolution. Nature 461, 515–519 (2009).

    Article  CAS  Google Scholar 

  45. Natarajan, C. et al. Epistasis among adaptive mutations in deer mouse hemoglobin. Science 340, 1324–1327 (2013).

    Article  CAS  Google Scholar 

  46. Rousset, F., Pélandakis, M. & Solignac, M. Evolution of compensatory substitutions through G.U intermediate state in Drosophila rRNA. Proc. Natl. Acad. Sci. USA 88, 10032–10036 (1991).

    Article  CAS  Google Scholar 

  47. Gell, C. et al. Single-molecule fluorescence resonance energy transfer assays reveal heterogeneous folding ensembles in a simple RNA stem–loop. J. Mol. Biol. 384, 264–278 (2008).

    Article  CAS  Google Scholar 

  48. Licatalosi, D.D. et al. HITS-CLIP yields genome-wide insights into brain alternative RNA processing. Nature 456, 464–469 (2008).

    Article  CAS  Google Scholar 

  49. Zhao, J. et al. Genome-wide identification of polycomb-associated RNAs by RIP-seq. Mol. Cell 40, 939–953 (2010).

    Article  CAS  Google Scholar 

  50. Buenrostro, J.D., Giresi, P.G., Zaba, L.C., Chang, H.Y. & Greenleaf, W.J. Transposition of native chromatin for fast and sensitive epigenomic profiling of open chromatin, DNA-binding proteins and nucleosome position. Nat. Methods 10, 1213–1218 (2013).

    Article  CAS  Google Scholar 

Download references

Acknowledgements

This work was supported by National Institutes of Health (NIH) NIH R01-HG004361 (to H.Y.C.); H.Y.C. is an Early Career Scientist of the Howard Hughes Medical Institute. J.D.B. and L.M.C. acknowledge support from the National Science Foundation Graduate Research Fellowships. J.D.B. also acknowledges support from NIH training grant T32HG000044. L.M.C. acknowledges support from NIH training grant T32GM067586. J.D.B. and L.M.C. each acknowledge support of a National Science Foundation graduate research fellowship. M.P.S. and C.L.A. acknowledge the NIH and the National Human Genome Research Institute (NHGRI) for funding through 5U54HG00455805. We thank D. Herschlag for feedback and advice throughout the methods development, and G. Sherlock and D. Herschlag for critical readings of earlier versions of this manuscript. We thank R. Landick for discussions regarding the design of our synthetic DNA library and K. Bajaj for discussions regarding quantification of cluster fluorescence. We also thank O.D. Phanstiel, M. Sikora and O. Cornejo for discussions on the modeling and evolutionary analyses.

Author information

Author notes

  1. Jason D Buenrostro and Carlos L Araya: These authors contributed equally to this work.

Authors and Affiliations

  1. Department of Genetics, Stanford University School of Medicine, Stanford, California, USA

    Jason D Buenrostro, Carlos L Araya, Lauren M Chircus, Curtis J Layton, Michael P Snyder & William J Greenleaf

  2. Program in Epithelial Biology and the Howard Hughes Medical Institute, Stanford University School of Medicine, Stanford, California, USA

    Jason D Buenrostro & Howard Y Chang

  3. Department of Chemical and Systems Biology, Stanford University School of Medicine, Stanford, California, USA

    Lauren M Chircus

Authors
  1. Jason D Buenrostro
    View author publications

    You can also search for this author in PubMed Google Scholar

  2. Carlos L Araya
    View author publications

    You can also search for this author in PubMed Google Scholar

  3. Lauren M Chircus
    View author publications

    You can also search for this author in PubMed Google Scholar

  4. Curtis J Layton
    View author publications

    You can also search for this author in PubMed Google Scholar

  5. Howard Y Chang
    View author publications

    You can also search for this author in PubMed Google Scholar

  6. Michael P Snyder
    View author publications

    You can also search for this author in PubMed Google Scholar

  7. William J Greenleaf
    View author publications

    You can also search for this author in PubMed Google Scholar

Contributions

W.J.G., J.D.B. and C.J.L. conceived of the method. J.D.B. developed the RNA display protocol. J.D.B. and L.M.C. designed and performed on-chip assays. L.M.C. designed and performed the protein purification and in vitro binding assays. J.D.B. wrote the image analysis algorithm with input from W.J.G. and C.J.L. C.L.A. developed and implemented the structural (epistatic), functional (modeling, kinetic) and evolutionary analyses. All authors interpreted the data and wrote the manuscript. W.J.G. supervised all aspects of this work.

Corresponding author

Correspondence to William J Greenleaf.

Ethics declarations

Competing interests

The authors declare no competing financial interests.

Supplementary information

Supplementary Text and Figures

Supplementary Figures 1–13 and Supplementary Discussion (PDF 4153 kb)

Supplementary Table 1

Oligonucleotide sequences used in this study. (XLSX 24 kb)

Supplementary Table 2

Measured binding energies and quality metrics for 129,248 MS2 RNA hairpin sequences. Note: Position indexes in -15,+3 indexing whereby. "NA:NA" indicates the consensus sequence (with zero mutations). (XLSX 21708 kb)

Supplementary Table 3

Measured dissociation and inferred association rates for 3,029 MS2 RNA hairpin sequences. Note: Position indexes in -15,+3 indexing whereby. "NA:NA" indicates the consensus sequence (with zero mutations). (XLSX 793 kb)

Supplementary Table 4

Summary of evolutionary path probabilities and constraint in 1,997 tesseracts. Note: Position indexes in -15,+3 indexing whereby. "NA:NA" indicates the consensus sequence (with zero mutations). (XLSX 217 kb)

Supplementary Data

Image analysis software. (ZIP 97 kb)

Equilibrium binding and dissociation rate measurements on a sequenced flow cell

A small region of the flow cell is shown as fluorescently tagged MS2 coat protein is bound at increasing concentrations to the RNA clusters and is then removed from solution to determine dissociation constants. In the first frame, the fluorescence signal from a complementary oligonucleotide annealed to all RNAs is shown. A single cluster is circled (blue), and quantified fluorescence is shown in the inset. Subsequent frames show the quantified fluorescence signal at increasing concentrations of labeled MS2 coat protein. At the end of the binding experiment, the fit for that individual cluster is shown in the inset. The following frames show dissociation of labeled MS2, measured by replacing labeled with unlabeled MS2 in solution, and observing the decay of fluorescence over time. The inset dissociation curve depicts data from all clusters that share the sequence with the circled cluster (-5C variant). (MOV 996 kb)

Rights and permissions

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Buenrostro, J., Araya, C., Chircus, L. et al. Quantitative analysis of RNA-protein interactions on a massively parallel array reveals biophysical and evolutionary landscapes. Nat Biotechnol 32, 562–568 (2014). https://doi.org/10.1038/nbt.2880

Download citation

  • Received:

  • Accepted:

  • Published:

  • Issue Date:

  • DOI: https://doi.org/10.1038/nbt.2880

Search

Advanced search

Quick links

Nature Briefing

Sign up for the Nature Briefing newsletter — what matters in science, free to your inbox daily.

Get the most important science stories of the day, free in your inbox. Sign up for Nature Briefing