Search

The discoveries of myriad non-coding RNA molecules, each transiting through multiple flexible states in cells or virions, present major challenges for structure determination. Advances in high-throughput chemical mapping give new routes for characterizing entire transcriptomes in vivo, but the resulting one-dimensional data generally remain too information-poor to allow accurate de novo structure determination. Multidimensional chemical mapping (MCM) methods seek to address this challenge. Mutate-and-map (M2 ), RNA interaction groups by mutational profiling (RING-MaP and MaP-2D analysis) and multiplexed •OH cleavage analysis (MOHCA) measure how the chemical reactivities of every nucleotide in an RNA molecule change in response to modifications at every other nucleotide. A growing body of in vitro blind tests and compensatory mutation/rescue experiments indicate that MCM methods give consistently accurate secondary structures and global tertiary structures for ribozymes, ribosomal domains and ligand-bound riboswitch aptamers up to 200 nucleotides in length. Importantly, MCM analyses provide detailed information on structurally heterogeneous RNA states, such as ligand-free riboswitches that are functionally important but difficult to resolve with other approaches. The sequencing requirements of currently available MCM protocols scale at least quadratically with RNA length, precluding general application to transcriptomes or viral genomes at present. We propose a modify-cross-link-map (MXM) expansion to overcome this and other current limitations to resolving the in vivo ‘RNA structurome’. [ABSTRACT FROM AUTHOR]

Landscapes exhibiting multiple secondary structures arise in natural RNA molecules that modulate gene expression, protein synthesis, and viral. We report herein that high-throughput chemical experiments can isolate an RNA’s multiple alternative secondary structures as they are stabilized by systematic mutagenesis (mutate-and-map, M2) and that a computational algorithm, REEFFIT, enables unbiased reconstruction of these states’ structures and populations. In an in silico benchmark on non-coding RNAs with complex landscapes, M2-REEFFIT recovers 95% of RNA helices present with at least 25% population while maintaining a low false discovery rate (10%) and conservative error estimates. In experimental benchmarks, M2-REEFFIT recovers the structure landscapes of a 35-nt MedLoop hairpin, a 110-nt 16S rRNA four-way junction with an excited state, a 25-nt bistable hairpin, and a 112-nt three-state adenine riboswitch with its expression platform, molecules whose characterization previously required expert mutational analysis and specialized NMR or chemical mapping experiments. With this validation, M2-REEFFIT enabled tests of whether artificial RNA sequences might exhibit complex landscapes in the absence of explicit design. An artificial flavin mononucleotide riboswitch and a randomly generated RNA sequence are found to interconvert between three or more states, including structures for which there was no design, but that could be stabilized through mutations. These results highlight the likely pervasiveness of rich landscapes with multiple secondary structures in both natural and artificial RNAs and demonstrate an automated chemical/computational route for their empirical characterization. [ABSTRACT FROM AUTHOR]

Consistently predicting biopolymer structure at atomic resolution from sequence alone remains a difficult problem, even for small sub-segments of large proteins. Such loop prediction challenges, which arise frequently in comparative modeling and protein design, can become intractable as loop lengths exceed 10 residues and if surrounding side-chain conformations are erased. Current approaches, such as the protein local optimization protocol or kinematic inversion closure (KIC) Monte Carlo, involve stages that coarse-grain proteins, simplifying modeling but precluding a systematic search of all-atom configurations. This article introduces an alternative modeling strategy based on a ‘stepwise ansatz’, recently developed for RNA modeling, which posits that any realistic all-atom molecular conformation can be built up by residue-by-residue stepwise enumeration. When harnessed to a dynamic-programming-like recursion in the Rosetta framework, the resulting stepwise assembly (SWA) protocol enables enumerative sampling of a 12 residue loop at a significant but achievable cost of thousands of CPU-hours. In a previously established benchmark, SWA recovers crystallographic conformations with sub-Angstrom accuracy for 19 of 20 loops, compared to 14 of 20 by KIC modeling with a comparable expenditure of computational power. Furthermore, SWA gives high accuracy results on an additional set of 15 loops highlighted in the biological literature for their irregularity or unusual length. Successes include cis-Pro touch turns, loops that pass through tunnels of other side-chains, and loops of lengths up to 24 residues. Remaining problem cases are traced to inaccuracies in the Rosetta all-atom energy function. In five additional blind tests, SWA achieves sub-Angstrom accuracy models, including the first such success in a protein/RNA binding interface, the YbxF/kink-turn interaction in the fourth ‘RNA-puzzle’ competition. These results establish all-atom enumeration as an unusually systematic approach to ab initio protein structure modeling that can leverage high performance computing and physically realistic energy functions to more consistently achieve atomic accuracy. [ABSTRACT FROM AUTHOR]

A complete macromolecule modeling package must be able to solve the simplest structure prediction problems. Despite recent successes in high resolution structure modeling and design, the Rosetta software suite fares poorly on small protein and RNA puzzles, some as small as four residues. To illustrate these problems, this manuscript presents Rosetta results for four well-defined test cases: the 20-residue mini-protein Trp cage, an even smaller disulfide-stabilized conotoxin, the reactive loop of a serine protease inhibitor, and a UUCG RNA tetraloop. In contrast to previous Rosetta studies, several lines of evidence indicate that conformational sampling is not the major bottleneck in modeling these small systems. Instead, approximations and omissions in the Rosetta all-atom energy function currently preclude discriminating experimentally observed conformations from de novo models at atomic resolution. These molecular ''puzzles'' should serve as useful model systems for developers wishing to make foundational improvements to this powerful modeling suite. [ABSTRACT FROM AUTHOR]

We present fragment assembly of RNA with full-atom refinement (FARFAR), a Rosetta framework for predicting and designing noncanonical motifs that define RNA tertiary structure. In a test set of thirty-two 6–20-nucleotide motifs, FARFAR recapitulated 50% of the experimental structures at near-atomic accuracy. Sequence redesign calculations recovered native bases at 65% of residues engaged in noncanonical interactions, and we experimentally validated mutations predicted to stabilize a signal recognition particle domain. [ABSTRACT FROM AUTHOR]

The prospect of phasing diffraction data sets 'de novo' for proteins with previously unseen folds is appealing but largely untested. In a first systematic exploration of phasing with Rosetta de novo models, it is shown that all-atom refinement of coarse-grained models significantly improves both the model quality and performance in molecular replacement with the Phaser software. 15 new cases of diffraction data sets that are unambiguously phased with de novo models are presented. These diffraction data sets represent nine space groups and span a large range of solvent contents (33-79%) and asymmetric unit copy numbers (1-4). No correlation is observed between the ease of phasing and the solvent content or asymmetric unit copy number. Instead, a weak correlation is found with the length of the modeled protein: larger proteins required somewhat less accurate models to give successful molecular replacement. Overall, the results of this survey suggest that de novo models can phase diffraction data for approximately one sixth of proteins with sizes of 100 residues or less. However, for many of these cases, 'de novo phasing with de novo models' requires significant investment of computational power, much greater than 103 CPU days per target. Improvements in conformational search methods will be necessary if molecular replacement with de novo models is to become a practical tool for targets without homology to previously solved protein structures. [ABSTRACT FROM AUTHOR]

DNA is thought to behave as a stiff elastic rod with respect to the ubiquitous mechanical deformations inherent to its biology. To test this model at short DNA lengths, we measured the mean and variance of end-to-end length for a series of DNA double helices in solution, using small-angle x-ray scattering interference between gold nanocrystal labels. In the absence of applied tension, DNA is at least one order of magnitude softer than measured by single-molecule stretching experiments. Further, the data rule out the conventional elastic rod model. The variance in end-to-end length follows a quadratic dependence on the number of base pairs rather than the expected linear dependence, indicating that DNA stretching is cooperative over more than two turns of the DNA double helix. Our observations support the idea of long-range allosteric communication through DNA structure. [ABSTRACT FROM AUTHOR]

We report a novel molecular ruler for measurement of distances and distance distributions with accurate external calibration. Using solution X-ray scattering we determine the scattering interference between two gold nanocrystal probes attached site-specifically to a macromolecule of interest. Fourier transformation of the interference pattern provides a model-independent probability distribution for the distances between the probe centers-of-mass. To test the approach, we measure end-to-end distances for a variety of DNA structures. We demonstrate that measurements with independently prepared samples and using different X-ray sources are highly reproducible, we demonstrate the quantitative accuracy of the first and second moments of the distance distributions, and we demonstrate that the technique recovers complex distribution shapes. Distances measured with the solution scattering-interference ruler match the corresponding crystallographic values, but differ from distances measured previously with alternate ruler techniques. The X-ray scattering interference ruler should be a powerful tool for relating crystal structures to solution structures and for studying molecular fluctuations. [ABSTRACT FROM AUTHOR]

Advances over the past few years have begun to enable prediction and design of macromolecular structures at near-atomic accuracy. Progress has stemmed from the development of reasonably accurate and efficiently computed all-atom potential functions as well as effective conformational sampling strategies appropriate for searching a highly rugged energy landscape, both driven by feedback from structure prediction and design tests. A unified energetic and kinematic framework in the Rosetta program allows a wide range of molecular modeling problems, from fibril structure prediction to RNA folding to the design of new protein interfaces, to be readily investigated and highlights areas for improvement. The methodology enables the creation of novel molecules with useful functions and holds promise for accelerating experimental structural inference. Emerging connections to crystallographic phasing, NMR modeling, and lower-resolution approaches are described and critically assessed. [ABSTRACT FROM AUTHOR]

The biological behaviors of ribozymes, riboswitches, and numerous other functional RNA molecules are critically dependent on their tertiary folding and their ability to sample multiple functional states. The conformational heterogeneity and partially folded nature of most of these states has rendered their characterization by high-resolution structural approaches difficult or even intractable. Here we introduce a method to rapidly infer the tertiary helical arrangements of large RNA molecules in their native and non-native solution states. Multiplexed hydroxyl radical (·OH) cleavage analysis (MOHCA) enables the high-throughput detection of numerous pairs of contacting residues via random incorporation of radical cleavage agents followed by two-dimensional gel electrophoresis. We validated this technology by recapitulating the unfolded and native states of a well studied model RNA, the P4-P6 domain of the Tetrahymena ribozyme, at subhelical resolution. We then applied MOHCA to a recently discovered third state of the P4-P6 RNA that is stabilized by high concentrations of monovalent salt and whose partial order precludes conventional techniques for structure determination. The three-dimensional portrait of a corn- pact, non-native RNA state reveals a well ordered subset of native tertiary contacts, in contrast to the dynamic but otherwise similar molten globule states of proteins. With its applicability to nearly any solution state, we expect MOHCA to be a powerful tool for illuminating the many functional structures of large RNA molecules and RNA/protein complexes. [ABSTRACT FROM AUTHOR]

RNA tertiary structure prediction has been based almost entirely on base-pairing constraints derived from phylogenetic covariation analysis. We describe here a complementary approach, inspired by the Rosetta low-resolution protein structure prediction method, that seeks the lowest energy tertiary structure for a given RNA sequence without using evolutionary information. In a benchmark test of 20 RNA sequences with known structure and lengths of ≈30 nt, the new method reproduces better than 90% of Watson-Crick base pairs, comparable with the accuracy of secondary structure prediction methods. In more than half the cases, at least one of the top five models agrees with the native structure to better than 4 Å rmsd over the backbone. Most importantly, the method recapitulates more than one-third of non-Watson-Crick base pairs seen in the native structures. Tandem stacks of "sheared" base pairs, base triplets, and pseudoknots are among the noncanonical features reproduced in the models. In the cases in which none of the top five models were native-like, higher energy conformations similar to the native structures are still sampled frequently but not assigned low energies. These results suggest that modest improvements in the energy function, together with the incorporation of information from phylogenetic covariance, may allow confident and accurate structure prediction for larger and more complex RNA chains. [ABSTRACT FROM AUTHOR]

Abstract: Riboswitches are complex folded RNA domains found in noncoding regions of mRNA that regulate gene expression upon small molecule binding. Recently, Breaker and coworkers reported a tandem aptamer riboswitch (VCI-II) that binds glycine cooperatively. Here, we use hydroxyl radical footprinting and small-angle X-ray scattering (SAXS) to study the conformations of this tandem aptamer as a function of Mg2+ and glycine concentration. We fit a simple three-state thermodynamic model that describes the energetic coupling between magnesium-induced folding and glycine binding. Furthermore, we characterize the structural conformations of each of the three states: In low salt with no magnesium present, the VCI-II construct has an extended overall conformation, presumably representing unfolded structures. Addition of millimolar concentrations of Mg2+ in the absence of glycine leads to a significant compaction and partial folding as judged by hydroxyl radical protections. In the presence of millimolar Mg2+ concentrations, the tandem aptamer binds glycine cooperatively. The glycine binding transition involves a further compaction, additional tertiary packing interactions and further uptake of magnesium ions relative to the state in high Mg2+ but no glycine. Employing density reconstruction algorithms, we obtain low resolution 3-D structures for all three states from the SAXS measurements. These data provide a first glimpse into the structural conformations of the VCI-II aptamer, establish rigorous constraints for further modeling, and provide a framework for future mechanistic studies. [Copyright &y& Elsevier]

Abstract: Like many structured RNAs, the Tetrahymena group I ribozyme is prone to misfolding. Here we probe a long-lived misfolded species, referred to as M, and uncover paradoxical aspects of its structure and folding. Previous work indicated that a non-native local secondary structure, termed alt P3, led to formation of M during folding in vitro. Surprisingly, hydroxyl radical footprinting, fluorescence measurements with site-specifically incorporated 2-aminopurine, and functional assays indicate that the native P3, not alt P3, is present in the M state. The paradoxical behavior of alt P3 presumably arises because alt P3 biases folding toward M, but, after commitment to this folding pathway and before formation of M, alt P3 is replaced by P3. Further, structural and functional probes demonstrate that the misfolded ribozyme contains extensive native structure, with only local differences between the two states, and the misfolded structure even possesses partial catalytic activity. Despite the similarity of these structures, re-folding of M to the native state is very slow and is strongly accelerated by urea, Na+, and increased temperature and strongly impeded by Mg2+ and the presence of native peripheral contacts. The paradoxical observations of extensive native structure within the misfolded species but slow conversion of this species to the native state are readily reconciled by a model in which the misfolded state is a topological isomer of the native state, and computational results support the feasibility of this model. We speculate that the complex topology of RNA secondary structures and the inherent rigidity of RNA helices render kinetic traps due to topological isomers considerably more common for RNA than for proteins. [Copyright &y& Elsevier]

Understanding biological and physical processes involving nucleic acids, such as the binding of proteins to DNA and RNA, DNA condensation, and RNA folding, requires an understanding of the ion atmosphere that surrounds nucleic acids. We have used a simple model DNA system to determine how the ion atmosphere modulates interactions between duplexes in the absence of specific metal ion-binding sites and other complicated interactions. In particular, we have tested whether the Coulomb repulsion between nucleic acids can be reversed by counterions to give a net attraction, as has been proposed recently for the rapid collapse observed early in KNA folding. The conformation of two DNA duplexes tethered by a flexible neutral linker was determined in the presence of a series of cations by small angle x-ray scattering. The small angle x-ray scattering profiles of two control molecules with distinct shapes (a continuous duplex and a mimic of the compact DNA) were in good agreement with predictions, establishing the applicability of this approach. Under low-salt conditions (20 mM Nat), the tethered duplexes are extended because of a Coulombic repulsion estimated to be 2-5 kT/bp. Addition of high concentrations of Na+ (1.2 M), Mg2+(0.6 M), and spermidine3+ (75 mM) resulted in electrostatic relaxation to a random state. These results indicate that a counterion-induced attractive force between nucleic acid duplexes is not significant under physiological conditions. An upper limit on the magnitude of the attractive potential under all tested ionic conditions is estimated. [ABSTRACT FROM AUTHOR]

Counterions are required for RNA folding, and divalent metal ions such as Mg2+ are often critical. To dissect the role of counterions, we have compared global and local folding of wild-type and mutant variants of P4-P6 RNA derived from the Tetrahymena group I ribozyme in monovalent and in divalent metal ions. A remarkably simple picture of the folding thermodynamics emerges. The equilibrium folding pathway in monovalent ions displays two phases. In the first phase, RNA molecules that are initially in an extended conformation enforced by charge–charge repulsion are relaxed by electrostatic screening to a state with increased flexibility but without formation of long-range tertiary contacts. At higher concentrations of monovalent ions, a state that is nearly identical to the native folded state in the presence of Mg2+ is formed, with tertiary contacts that involve base and backbone interactions but without the subset of interactions that involve specific divalent metal ion-binding sites. The folding model derived from these and previous results provides a robust framework for understanding the equilibrium and kinetic folding of RNA. [Copyright &y& Elsevier]

Large RNAs can collapse into compact conformations well before the stable formation of the tertiary contacts that define their final folds. This study identifies likely physical mechanisms driving these early compaction events in RNA folding. We have employed time-resolved small-angle X-ray scattering to monitor the fastest global shape changes of the Tetrahymena ribozyme under different ionic conditions and with RNA mutations that remove long-range tertiary contacts. A partial collapse in each of the folding time-courses occurs within tens of milliseconds with either monovalent or divalent cations. Combined with comparison to predictions from structural models, this observation suggests a relaxation of the RNA to a more compact but denatured conformational ensemble in response to enhanced electrostatic screening at higher ionic concentrations. Further, the results provide evidence against counterion-correlation-mediated attraction between RNA double helices, a recently proposed model for early collapse. A previous study revealed a second 100 ms phase of collapse to a globular state. Surprisingly, we find that progression to this second early folding intermediate requires RNA sequence motifs that eventually mediate native long-range tertiary interactions, even though these regions of the RNA were observed to be solvent-accessible in previous footprinting studies under similar conditions. These results help delineate an analogy between the early conformational changes in RNA folding and the “burst phase” changes and molten globule formation in protein folding. [Copyright &y& Elsevier]

Describes a maximum-likelihood approach to removing radio sources from observations of the Sunyaev-Zel'dovich effect. Removal of contaminating radio sources; Cosmology; Galaxies.

Hypothesis-driven experimentation – the scientific method – can be subverted by fraud, irreproducibility, and lack of rigorous predictive tests. A robust solution to these problems may be the ‘massive open laboratory’ model, recently embodied in the internet-scale videogame EteRNA. Deploying similar platforms throughout biology could enforce the scientific method more broadly. [ABSTRACT FROM AUTHOR]

We propose a rapid chemical strategy for identifying base pairs in structured nucleic acid systems. The approach goes beyond traditional chemical mapping approaches by monitoring perturbations of each residue's chemical accessibility in response to systematic mutagenesis of residues that are distant in sequence but nearby in three dimensions. As a proof of concept, we present high-throughput dimethyl sulfate accessibility data for a chimeric DNA/RNA system in which every possible sequence variation and deletion in a 20 bp region has been synthesized and tested. The data demonstrate that 88% of the system's base pairs can be robustly inferred, with A/A and T/C DNA/RNA mismatches giving the strongest signals. These results point to the feasibility of rapid base pair inference in larger and more complex nucleic acid systems with unknown structure. [ABSTRACT FROM AUTHOR]

Coronavirus genomes sequester their start codons within stem-loop 5 (SL5), a structured, 5' genomic RNA element. In most alpha- and betacoronaviruses, the secondary structure of SL5 is predicted to contain a four-way junction of helical stems, some of which are capped with UUYYGU hexaloops. Here, using cryogenic electron micros-copy (cryo-EM) and computational modeling with biochemically determined secondary structures, we present three-dimensional structures of SL5 from six coronaviruses. The SL5 domain of betacoronavirus severe-acute-respiratory-syndrome-related coronavirus 2 (SARS-CoV-2), resolved at 4.7 Å resolution, exhibits a T-shaped structure, with its UUYYGU hexaloops at opposing ends of a coaxial stack, the T's "arms." Further analysis of SL5 domains from SARS-CoV-1 and MERS (7.1 and 6.4 to 6.9 Å resolution, respectively) indicate that the junction geometry and inter-hexaloop distances are conserved features across these human-infecting betacoronaviruses. The MERS SL5 domain displays an additional tertiary interaction, which is also observed in the non-human-infecting betacoronavirus BtCoV-HKU5 (5.9 to 8.0 Å resolution). SL5s from human-infecting alphacoronaviruses, HCoV-229E and HCoV-NL63 (6.5 and 8.4 to 9.0 A resolution, respectively), exhibit the same coaxial stacks, including the UUYYGU-capped arms, but with a phylogenetically distinct crossing angle, an X-shape. As such, all SL5 domains studied herein fold into stable tertiary structures with cross-genus similarities and notable differences, with implications for potential protein-binding modes and therapeutic targets. [ABSTRACT FROM AUTHOR]

Rapid RNA synthesis of comprehensive single mutant libraries and targeted multiple mutant libraries is enabling new multidimensional chemical approaches to solve RNA structures. PCR assembly of DNA templates and in vitro transcription allow synthesis and purification of hundreds of RNA mutants in a cost-effective manner, with sharing of primers across constructs allowing significant reductions in expense. However, these protocols require organization of primer locations across numerous 96 well plates and guidance for pipetting, non-trivial tasks for which informatics and visualization tools can prevent costly errors. We report here an online tool to accelerate synthesis of large libraries of desired mutants through design and efficient organization of primers. The underlying program and graphical interface have been experimentally tested in our laboratory for RNA domains with lengths up to 300 nucleotides and libraries encompassing up to 960 variants. In addition to the freely available Primerize-2D server, the primer design code is available as a stand-alone Python package for broader applications. [ABSTRACT FROM AUTHOR]

Summary RNA-protein complexes underlie numerous cellular processes including translation, splicing, and posttranscriptional regulation of gene expression. The structures of these complexes are crucial to their functions but often elude high-resolution structure determination. Computational methods are needed that can integrate low-resolution data for RNA-protein complexes while modeling de novo the large conformational changes of RNA components upon complex formation. To address this challenge, we describe RNP-denovo , a Rosetta method to simultaneously fold-and-dock RNA to a protein surface. On a benchmark set of diverse RNA-protein complexes not solvable with prior strategies, RNP-denovo consistently sampled native-like structures with better than nucleotide resolution. We revisited three past blind modeling challenges involving the spliceosome, telomerase, and a methyltransferase-ribosomal RNA complex in which previous methods gave poor results. When coupled with the same sparse FRET, crosslinking, and functional data used previously, RNP-denovo gave models with significantly improved accuracy. These results open a route to modeling global folds of RNA-protein complexes from low-resolution data. Graphical Abstract Highlights • A Rosetta method, RNP-denovo , folds and docks RNA to protein surfaces • RNP-denovo consistently samples native-like structures for a diverse benchmark set • Applying RNP-denovo to past modeling problems gives models with improved accuracy RNA and protein molecules interact intimately during ribosome translation, pre-mRNA splicing, and other fundamental cellular processes. Structures of the relevant RNA-protein complexes often elude high-resolution determination. Kappel and Das report a computational Rosetta method RNP-denovo that can incorporate sparse experimental data and sample native-like structures of RNA-protein complexes. [ABSTRACT FROM AUTHOR]

The article reports that in living organisms and in test tubes, structured RNA molecules make intimate interactions with divalent metal cations, especially magnesium. Numerous thermodynamic, structural, and biochemical studies have probed important aspects of metal ion association with nucleic acids, and there is a rich literature on the theory of ion-polyelectrolyte interactions. One widely used approach to estimate metal ion stoichiometry for RNA folding events has relied on fitting equilibrium data to the Hill equation.

The structural interconversions that mediate the gene regulatory functions of RNA molecules may be different from classic models of allostery, but the relevant structural correlations have remained elusive in even intensively studied systems. Here, we present a four-dimensional expansion of chemical mapping called lock-mutate-map-rescue (LM²R), which integrates multiple layers of mutation with nucleotide-resolution chemical mapping. This technique resolves the core mechanism of the adenine-responsive V. vulnificus add riboswitch, a paradigmatic system for which both Monod-Wyman-Changeux (MWC) conformational selection models and non-MWC alternatives have been proposed. To discriminate amongst these models, we locked each functionally important helix through designed mutations and assessed formation or depletion of other helices via compensatory rescue evaluated by chemical mapping. These LM²R measurements give strong support to the pre-existing correlations predicted by MWC models, disfavor alternative models, and suggest additional structural heterogeneities that may be general across ligand-free riboswitches. [ABSTRACT FROM AUTHOR]

The article presents a response by researcher Rebecca S. Mathew-Fenn and colleagues regarding comments made by Nils B. Becker and Ralf Everaers on their report "Remeasuring the Double Helix." The authors cite their remeasurement of the mean and variance of the length of short DNA fragments in solution and reiterate the presentation of three experimental tests explaining their explanation for DNA bending.

Internet-based scientific communities promise a means to apply distributed, diverse human intelligence toward previously intractable scientific problems. However, current implementations have not allowed communities to propose experiments to test all emerging hypotheses at scale or to modify hypotheses in response to experiments. We report high-throughput methods for molecular characterization of nucleic acids that enable the large-scale video game–based crowdsourcing of RNA sensor design, followed by high-throughput functional characterization. Iterative design testing of thousands of crowdsourced RNA sensor designs produced near–thermodynamically optimal and reversible RNA switches that act as self-contained molecular sensors and couple five distinct small molecule inputs to three distinct protein binding and fluorogenic outputs. This work suggests a paradigm for widely distributed experimental bioscience. [ABSTRACT FROM AUTHOR]

Single-molecule tweezers measurements of double-stranded nucleic acids (dsDNA and dsRNA) provide unprecedented opportunities to dissect how these fundamental molecules respond to forces and torques analogous to those applied by topoisomerases, viral capsids, and other biological partners. However, tweezers data are still most commonly interpreted post facto in the framework of simple analytical models. Testing falsifiable predictions of state-of-the-art nucleic acid models would be more illuminating but has not been performed. Here we describe a blind challenge in which numerical predictions of nucleic acid mechanical properties were compared to experimental data obtained recently for dsRNA under applied force and torque. The predictions were enabled by the HelixMC package, first presented in this paper. HelixMC advances crystallography-derived base-pair level models (BPLMs) to simulate kilobase-length dsDNAs and dsRNAs under external forces and torques, including their global linking numbers. These calculations recovered the experimental bending persistence length of dsRNA within the error of the simulations and accurately predicted that dsRNA's “spring-like” conformation would give a two-fold decrease of stretch modulus relative to dsDNA. Further blind predictions of helix torsional properties, however, exposed inaccuracies in current BPLM theory, including three-fold discrepancies in torsional persistence length at the high force limit and the incorrect sign of dsRNA link-extension (twist-stretch) coupling. Beyond these experiments, HelixMC predicted that ‘nucleosome-excluding’ poly(A)/poly(T) is at least two-fold stiffer than random-sequence dsDNA in bending, stretching, and torsional behaviors; Z-DNA to be at least three-fold stiffer than random-sequence dsDNA, with a near-zero link-extension coupling; and non-negligible effects from base pair step correlations. We propose that experimentally testing these predictions should be powerful next steps for understanding the flexibility of dsDNA and dsRNA in sequence contexts and under mechanical stresses relevant to their biology. [ABSTRACT FROM AUTHOR]

Ions surround nucleic acids in what is referred to as an ion atmosphere. As a result, the folding and dynamics of RNA and DNA and their complexes with proteins and with each other cannot be understood without a reasonably sophisticated appreciation of these ions' electrostatic interactions. However, the underlying behavior of the ion atmosphere follows physical rules that are distinct from the rules of site binding that biochemists are most familiar and comfortable with. The main goal of this review is to familiarize nucleic acid experimentalists with the physical concepts that underlie nucleic acid-ion interactions. Throughout, we provide practical strategies for interpreting and analyzing nucleic acid experiments that avoid pitfalls from oversimplified or incorrect models. We briefly review the status of theories that predict or simulate nucleic acid-ion interactions and experiments that test these theories. Finally, we describe opportunities for going beyond phenomenological fits to a next-generation, truly predictive understanding of nucleic acid-ion interactions. [ABSTRACT FROM AUTHOR]

The tertiary structures of functional RNA molecules remain difficult to decipher. A new generation of automated RNA structure prediction methods may help address these challenges but have not yet been experimentally validated. Here we apply four prediction tools to a class of double glycine riboswitches that can bind two ligands cooperatively. A novel method (BPPalign), RMdetect, JAR3D, and Rosetta 3D modeling give consistent predictions for a new stem P0 and a kink-turn motif. These elements structure the linker between the RNAs' double aptamers. Chemical mapping on the Fusobacterium nucleatum riboswitch with N-methylisatoic anhydride, dimethyl sulfate and 1-cyclohexyl-3-(2-morpholinoethyl)carbodiimide metho-p-toluenesulfonate probing, mutate-and-map studies, and mutation/rescue experiments all provide strong evidence for the structured linker. Under solution conditions that permit rigorous thermodynamic analysis, disrupting this helix-junction-helix structure gives 120- and 6-30-fold poorer dissociation constants for the RNA's two glycine-binding transitions, corresponding to an overall energetic impact of 4.3 ± 0.5 kcal/mol. Prior biochemical and crystallography studies did not include this critical element due to over-truncation of the RNA. We speculate that several further undiscovered elements are likely to exist in the flanking regions of this and other functional RNAs, and automated prediction tools can play a useful role in their detection and dissection. [ABSTRACT FROM AUTHOR]

Atomic-accuracy structure prediction of macromolecules should be achievable by optimizing a physically realistic energy function but is presently precluded by incomplete sampling of a biopolymer's many degrees of freedom. We present herein a working hypothesis, called the "stepwise ansatz," for recursively constructing well-packed atomic-detail models in small steps, enumerating several million conformations for each monomer, and covering all build-up paths. By making use of high-performance computing and the Rosetta framework, we provide first tests of this hypothesis on a benchmark of 15 RNA loop-modeling problems drawn from riboswitches, ribozymes, and the ribosome, including 10 cases that are not solvable by current knowledge-based modeling approaches. For each loop problem, this deterministic stepwise assembly method either reaches atomic accuracy or exposes flaws in Rosetta's all-atom energy function, indicating the resolution of the conformational sampling bottleneck. As a further rigorous test, we have carried out a blind all-atom prediction for a noncanonical RNA motif, the C7.2 tetraloop/receptor, and validated this model through nucleotide-resolution chemical mapping experiments. Stepwise assembly is an enumerative, ab initio build-up method that systematically outperforms existing Monte Carlo and knowledge-based methods for 3D structure prediction. [ABSTRACT FROM AUTHOR]

As the fastest folding protein, the villin headpiece (HP35) serves as an important bridge between simulation and experimental studies of protein folding. Despite the simplicity of this system, experiments continue to reveal a number of surprises, including structure in the unfolded state and complex equilibrium dynamics near the native state. Using 2.5 ms of molecular dynamics and Markov state models, we connect to current experimental results in three ways. First, we present and validate a novel method for the quantitative prediction of triplet-triplet energy transfer experiments. Second, we construct a many-state model for HP35 that is consistent with previous experiments. Finally, we predict contact-formation time traces for all 1,225 possible triplet-triplet energy transfer experiments on HP35. [ABSTRACT FROM AUTHOR]

RNA molecules adopt three-dimensional structures that are critical to their function and of interest in drug discovery. Few RNA structures are known, however, and predicting them computationally has proven challenging. We introduce a machine learning approach that enables identification of accurate structural models without assumptions about their defining characteristics, despite being trained with only 18 known RNA structures. The resulting scoring function, the Atomic Rotationally Equivariant Scorer (ARES), substantially outperforms previous methods and consistently produces the best results in community-wide blind RNA structure prediction challenges. By learning effectively even from a small amount of data, our approach overcomes a major limitation of standard deep neural networks. Because it uses only atomic coordinates as inputs and incorporates no RNA-specific information, this approach is applicable to diverse problems in structural biology, chemistry, materials science, and beyond. [ABSTRACT FROM AUTHOR]

Docking of the P1 duplex into the pre-folded core of the Tetrahymena group I ribozyme exemplifies the formation of tertiary interactions in the context of a complex, structured RNA. We have applied Φ-analysis to P1 docking, which compares the effects of modifications on the rate constant for docking (kdock) with the effects on the docking equilibrium (Kdock). To accomplish this we used a single molecule fluorescence resonance energy transfer assay that allows direct determination of the rate constants for formation of thermodynamically favorable, as well as unfavorable, states. Modification of the eight groups of the P1 duplex that make tertiary interactions with the core and changes in solution conditions decrease Kdock up to 500-fold, whereas kdock changes by ≤2-fold. The absence of effects on kdock, both from atomic modifications and global perturbations, strongly suggests that the transition state for docking is early and does not closely resemble the docked state. These results, the slow rate of docking of 3 s−1, and the observation that a modification that is expected to increase the degrees of freedom between the P1 duplex and the ribozyme core accelerates docking, suggest a model in which a kinetic trap(s) slows docking substantially. Nonetheless, urea does not increase kdock, suggesting that there is little change in the exposed surface area between the trapped, undocked state and the transition state. The findings highlight that urea and temperature dependencies can be inadequate to diagnose the presence of kinetic traps in a folding process. The results described here, combined with previous work, provide an in-depth view of an RNA tertiary structure formation event and suggest that large, highly structured RNAs may have local regions that are misordered. [Copyright &y& Elsevier]

Predicting RNA three-dimensional structures from sequence could accelerate understanding of the growing number of RNA molecules being discovered across biology. Rosetta's Fragment Assembly of RNA with Full-Atom Refinement (FARFAR) has shown promise in community-wide blind RNA-Puzzle trials, but lack of a systematic and automated benchmark has left unclear what limits FARFAR performance. Here, we benchmark FARFAR2, an algorithm integrating RNA-Puzzle-inspired innovations with updated fragment libraries and helix modeling. In 16 of 21 RNA-Puzzles revisited without experimental data or expert intervention, FARFAR2 recovers native-like structures more accurate than models submitted during the RNA-Puzzles trials. Remaining bottlenecks include conformational sampling for >80-nucleotide problems and scoring function limitations more generally. Supporting these conclusions, preregistered blind models for adenovirus VA-I RNA and five riboswitch complexes predicted native-like folds with 3- to 14 Å root-mean-square deviation accuracies. We present a FARFAR2 webserver and three large model archives (FARFAR2-Classics, FARFAR2-Motifs, and FARFAR2-Puzzles) to guide future applications and advances. • A new RNA fragment assembly method, FARFAR2, uses RNA-Puzzle-inspired innovations • FARFAR2 recovers native-like structures on six new blind challenges • A FARFAR2 webserver is available at https://rosie.rosettacommons.org/farfar2 • A dataset of millions of FARFAR2-predicted models is available for further study Watkins et al. benchmark FARFAR2, an algorithm for RNA structure prediction that consolidates a decade of ad hoc strategies. FARFAR2 recovers native-like structures more accurately than original challenge submissions for 16 of 21 RNA-Puzzles and achieves native-like preregistered blind predictions for adenovirus VA-I RNA and five riboswitch complexes. [ABSTRACT FROM AUTHOR]

Summary: We have established an RNA mapping database (RMDB) to enable structural, thermodynamic and kinetic comparisons across single-nucleotide-resolution RNA structure mapping experiments. The volume of structure mapping data has greatly increased since the development of high-throughput sequencing techniques, accelerated software pipelines and large-scale mutagenesis. For scientists wishing to infer relationships between RNA sequence/structure and these mapping data, there is a need for a database that is curated, tagged with error estimates and interfaced with tools for sharing, visualization, search and meta-analysis. Through its on-line front-end, the RMDB allows users to explore single-nucleotide-resolution mapping data in heat-map, bar-graph and colored secondary structure graphics; to leverage these data to generate secondary structure hypotheses; and to download the data in standardized and computer-friendly files, including the RDAT and community-consensus SNRNASM formats. At the time of writing, the database houses 53 entries, describing more than 2848 experiments of 1098 RNA constructs in several solution conditions and is growing rapidly.Availability: Freely available on the web at http://rmdb.stanford.eduContact: rhiju@stanford.eduSupplementary information: Supplementary data are available at Bioinformatics Online. [ABSTRACT FROM AUTHOR]

Background: The understanding of the importance of RNA has dramatically changed over recent years. As in the case of proteins, the function of an RNA molecule is encoded in its tertiary structure, which in turn is determined by the molecule's sequence. The prediction of tertiary structures of complex RNAs is still a challenging task. Results: Using the observation that RNA sequences from the same RNA family fold into conserved structure, we test herein whether parallel modeling of RNA homologs can improve ab initio RNA structure prediction. EvoClustRNA is a multi-step modeling process, in which homologous sequences for the target sequence are selected using the Rfam database. Subsequently, independent folding simulations using Rosetta FARFAR and SimRNA are carried out. The model of the target sequence is selected based on the most common structural arrangement of the common helical fragments. As a test, on two blind RNA-Puzzles challenges, EvoClustRNA predictions ranked as the first of all submissions for the L-glutamine riboswitch and as the second for the ZMP riboswitch. Moreover, through a benchmark of known structures, we discovered several cases in which particular homologs were unusually amenable to structure recovery in folding simulations compared to the single original target sequence. Conclusion: This work, for the first time to our knowledge, demonstrates the importance of the selection of the target sequence from an alignment of an RNA family for the success of RNA 3D structure prediction. These observations prompt investigations into a new direction of research for checking 3D structure "foldability" or "predictability" of related RNA sequences to obtain accurate predictions. To support new research in this area, we provide all relevant scripts in a documented and ready-to-use form. By exploring new ideas and identifying limitations of the current RNA 3D structure prediction methods, this work is bringing us closer to the near-native computational RNA 3D models. [ABSTRACT FROM AUTHOR]

Structured RNAs and RNA complexes underlie biological processes ranging from control of gene expression to protein translation. Approximately 50% of nucleotides within known structured RNAs are folded intoWatson–Crick (WC) base pairs, and sequence changes that preserve these pairs are typically assumed to preserve higherorder RNA structure and binding of macromolecule partners. Here, we report that indirect effects of the helix sequence on RNA tertiary stability are, in fact, significant but are nevertheless predictable from a simple computational model called RNAMake-ΔΔG. When tested through the RNA on a massively parallel array (RNA-MaP) experimental platform, blind predictions for >1500 variants of the tectoRNA heterodimer model system achieve high accuracy (rmsd 0.34 and 0.77 kcal/mol for sequence and length changes, respectively). Detailed comparison of predictions to experiments support a microscopic picture of how helix sequence changes subtly modulate conformational fluctuations at each base-pair step, which accumulate to impact RNA tertiary structure stability. Our study reveals a previously overlooked phenomenon in RNA structure formation and provides a framework of computation and experiment for understanding helix conformational preferences and their impact across biological RNA and RNAprotein assemblies. [ABSTRACT FROM AUTHOR]

Emerging RNA-based approaches to disease detection and gene therapy require RNA sequences that fold into specific base-pairing patterns, but computational algorithms generally remain inadequate for these secondary structure design tasks. The Eterna project has crowdsourced RNA design to human video game players in the form of puzzles that reach extraordinary difficulty. Here, we demonstrate that Eterna participants’ moves and strategies can be leveraged to improve automated computational RNA design. We present an eternamoves-large repository consisting of 1.8 million of player moves on 12 of the most-played Eterna puzzles as well as an eternamoves-select repository of 30,477 moves from the top 72 players on a select set of more advanced puzzles. On eternamoves-select, we present a multilayer convolutional neural network (CNN) EternaBrain that achieves test accuracies of 51% and 34% in base prediction and location prediction, respectively, suggesting that top players’ moves are partially stereotyped. Pipelining this CNN’s move predictions with single-action-playout (SAP) of six strategies compiled by human players solves 61 out of 100 independent puzzles in the Eterna100 benchmark. EternaBrain-SAP outperforms previously published RNA design algorithms and achieves similar or better performance than a newer generation of deep learning methods, while being largely orthogonal to these other methods. Our study provides useful lessons for future efforts to achieve human-competitive performance with automated RNA design algorithms. [ABSTRACT FROM AUTHOR]

Interactions between RNA and proteins are pervasive in biology, driving fundamental processes such as protein translation and participating in the regulation of gene expression. Modeling the energies of RNA-protein interactions is therefore critical for understanding and repurposing living systems but has been hindered by complexities unique to RNA-protein binding. Here, we bring together several advances to complete a calculation framework for RNA-protein binding affinities, including a unified free energy function for bound complexes, automated Rosetta modeling of mutations, and use of secondary structure-based energetic calculations to model unbound RNA states. The resulting Rosetta-Vienna RNP-ΔΔG method achieves root-mean-squared errors (RMSEs) of 1.3 kcal/mol on highthroughput MS2 coat protein-RNA measurements and 1.5 kcal/mol on an independent test set involving the signal recognition particle, human U1A, PUM1, and FOX-1. As a stringent test, the method achieves RMSE accuracy of 1.4 kcal/mol in blind predictions of hundreds of human PUM2-RNA relative binding affinities. Overall, these RMSE accuracies are significantly better than those attained by prior structure-based approaches applied to the same systems. Importantly, Rosetta-Vienna RNP-ΔΔG establishes a framework for further improvements in modeling RNA-protein binding that can be tested by prospective high-throughput measurements on new systems. [ABSTRACT FROM AUTHOR]

Phase separation of multivalent protein and RNA molecules underlies the biogenesis of biomolecular condensates such as membraneless organelles. In vivo, these condensates encompass hundreds of distinct types of molecules that typically organize into multilayered structures supporting the differential partitioning of molecules into distinct regions with distinct material properties. The interplay between driven (active) versus spontaneous (passive) processes that are required for enabling the formation of condensates with coexisting layers of distinct material properties remains unclear. Here, we deploy systematic experiments and simulations based on coarse-grained models to show that the collective interactions among the simplest, biologically relevant proteins and archetypal RNA molecules are sufficient for driving the spontaneous emergence of multilayered condensates with distinct material properties. These studies yield a set of rules regarding homotypic and heterotypic interactions that are likely to be relevant for understanding the interplay between active and passive processes that control the formation of functional biomolecular condensates. [ABSTRACT FROM AUTHOR]

The article focuses on the prediction of RNA structure from nucleotide sequence. It focuses on modeling of noncanonical base pairs of new RNA structural motifs. It offers information on stepwise Monte Carlo (SWM)method with add-and-delete move set that allows predictions of noncanonical base pairs of complex RNA structures.

Despite the critical roles RNA structures play in regulating gene expression, sequencing-based methods for experimentally determining RNA base pairs have remained inaccurate. Here, we describe a multidimensional chemical-mapping method called "mutate-and-map read out through next-generation sequencing" (M2-seq) that takes advantage of sparsely mutated nucleotides to induce structural perturbations at partner nucleotides and then detects these events through dimethyl sulfate (DMS) probing and mutational profiling. In special cases, fortuitous errors introduced during DNA template preparation and RNA transcription are sufficient to give M2-seq helix signatures; these signals were previously overlooked or mistaken for correlated double-DMS events. When mutations are enhanced through error-prone PCR, in vitro M2-seq experimentally resolves 33 of 68 helices in diverse structured RNAs including ribozyme domains, riboswitch aptamers, and viral RNA domains with a single false positive. These inferences do not require energy minimization algorithms and can be made by either direct visual inspection or by a neural-network-inspired algorithm called M2-net. Measurements on the P4-P6 domain of the Tetrahymena group I ribozyme embedded in Xenopus egg extract demonstrate the ability of M2-seq to detect RNA helices in a complex biological environment. [ABSTRACT FROM AUTHOR]

The predictive modeling and design of biologically active RNA molecules requires understanding the energetic balance among their basic components. Rapid developments in computer simulation promise increasingly accurate recovery of RNA's nearest-neighbor (NN) freeenergy parameters, but these methods have not been tested in predictive trials or on nonstandard nucleotides. Here, we present, to our knowledge, the first such tests through a RECCES-Rosetta (reweighting of energy-function collection with conformational ensemble sampling in Rosetta) framework that rigorously models conformational entropy, predicts previously unmeasured NN parameters, and estimates these values' systematic uncertainties. RECCES-Rosetta recovers the 10 NN parameters for Watson-Crick stacked base pairs and 32 single-nucleotide dangling-end parameters with unprecedented accuracies: rmsd of 0.28 kcal/mol and 0.41 kcal/mol, respectively. For setaside test sets, RECCES-Rosetta gives rmsd values of 0.32 kcal/mol on eight stacked pairs involving G-U wobble pairs and 0.99 kcal/mol on seven stacked pairs involving nonstandard isocytidine-isoguanosine pairs. To more rigorously assess RECCES-Rosetta,we carried out four blind predictions for stacked pairs involving 2,6-diaminopurine-U pairs, which achieved 0.64 kcal/mol rmsd accuracy when tested by subsequent experiments. Overall, these results establish that computational methods can now blindly predict energetics of basic RNA motifs, including chemically modified variants, with consistently better than 1 kcal/mol accuracy. Systematic tests indicate that resolving the remaining discrepancies will require energy function improvements beyond simply reweighting component terms, and we propose further blind trials to test such efforts. [ABSTRACT FROM AUTHOR]

Designing RNAs that form specific secondary structures is enabling better understanding and control of living systems through RNA-guided silencing, genome editing and protein organization. Little is known, however, about which RNA secondary structures might be tractable for downstream sequence design, increasing the time and expense of design efforts due to inefficient secondary structure choices. Here, we present insights into specific structural features that increase the difficulty of finding sequences that fold into a target RNA secondary structure, summarizing the design efforts of tens of thousands of human participants and three automated algorithms (RNAInverse, INFO-RNA and RNA-SSD) in the Eterna massive open laboratory. Subsequent tests through three independent RNA design algorithms (NUPACK, DSS-Opt and MODENA) confirmed the hypothesized importance of several features in determining design difficulty, including sequence length, mean stem length, symmetry and specific difficult-to-design motifs such as zigzags. Based on these results, we have compiled an Eterna100 benchmark of 100 secondary structure design challenges that span a large range in design difficulty to help test future efforts. Our in silico results suggest new routes for improving computational RNA design methods and for extending these insights to assess “designability” of single RNA structures, as well as of switches for in vitro and in vivo applications. [ABSTRACT FROM AUTHOR]

I130 --> 115.3 --> Linear eukaryotic chromosome ends are capped with telomeres, which play key roles in maintaining genome integrity. Telomeres are, however, progressively shortened because of incomplete replication by the replisome. Critically short telomere length leads to proliferative senescence and cell death. To balance this loss, telomerase ribonucleoprotein (RNP) synthesizes telomeric repeats through its telomerase reverse transcriptase (TERT) and an integral RNA subunit (hTR) carrying the template for repeat synthesis. In addition to TERT and hTR, human telomerase holoenzyme consists of a number of other protein factors required for RNP assembly and localization. We used cryo‐EM to determine the structure of DNA substrate‐bound human telomerase holoenzyme first at medium resolution and more recently at sub‐4Å resolution. The structure revealed not only detailed insight into telomerase active site and assembly but also novel histone subunits. Here I will discuss our new findings and its implications. [ABSTRACT FROM AUTHOR]

The article presents a study that determines the role of multidimensional chemical mapping (MCM) pipeline in providing sequencing-based route to uncover non-coding RNA (ncRNA) 3D structure which is essential in the heterogenous state of chemical mapping

Emerging evidence suggests that the ribosome has a regulatory function in directing how the genome is translated in time and space. However, how this regulation is encoded in the messenger RNA sequence remains largely unknown. Here we uncover unique RNA regulons embedded in homeobox (Hox) 5′ untranslated regions (UTRs) that confer ribosome-mediated control of gene expression. These structured RNA elements, resembling viral internal ribosome entry sites (IRESs), are found in subsets of Hox mRNAs. They facilitate ribosome recruitment and require the ribosomal protein RPL38 for their activity. Despite numerous layers of Hox gene regulation, these IRES elements are essential for converting Hox transcripts into proteins to pattern the mammalian body plan. This specialized mode of IRES-dependent translation is enabled by an additional regulatory element that we term the translation inhibitory element (TIE), which blocks cap-dependent translation of transcripts. Together, these data uncover a new paradigm for ribosome-mediated control of gene expression and organismal development. [ABSTRACT FROM AUTHOR]

RNA plays myriad roles in the transmission and regulation of genetic information that are fundamentally constrained by its mechanical properties, including the elasticity and conformational transitions of the double-stranded (dsRNA) form. Although double-stranded DNA (dsDNA) mechanics have been dissected with exquisite precision, much less is known about dsRNA. Here we present a comprehensive characterization of dsRNA under external forces and torques using magnetic tweezers. We find that dsRNA has a force–torque phase diagram similar to that of dsDNA, including plectoneme formation, melting of the double helix induced by torque, a highly overwound state termed "P-RNA," and a highly underwound, left-handed state denoted "L-RNA." Beyond these similarities, our experiments reveal two unexpected behaviors of dsRNA: Unlike dsDNA, dsRNA shortens upon overwinding, and its characteristic transition rate at the plectonemic buckling transition is two orders of magnitude slower than for dsDNA. Our results challenge current models of nucleic acid mechanics, provide a baseline for modeling RNAs in biological contexts, and pave the way for new classes of magnetic tweezers experiments to dissect the role of twist and torque for RNA–protein interactions at the single-molecule level. [ABSTRACT FROM AUTHOR]

Structured noncoding RNAs underlie fundamental cellular processes, but determining their three-dimensional structures remains challenging. We demonstrate that integrating 1H NMR chemical shift data with Rosetta de novo modeling can be used to consistently determine high-resolution RNA structures. On a benchmark set of 23 noncanonical RNA motifs, including 11 'blind' targets, chemical-shift Rosetta for RNA (CS-Rosetta-RNA) recovered experimental structures with high accuracy (0.6-2.0 Å all-heavy-atom r.m.s. deviation) in 18 cases. [ABSTRACT FROM AUTHOR]

Self-assembling RNA molecules present compelling substrates for the rational interrogation and control of living systems. However, imperfect in silico models—even at the secondary structure level— hinder the design of new RNAs that function properly when synthesized. Here, we present a unique and potentially general approach to such empirical problems: the Massive Open Laboratory. The EteRNA project connects 37,000 enthusiasts to RNA design puzzles through an online interface. Uniquely, EteRNA participants not only manipulate simulated molecules but also control a remote experimental pipeline for high-throughput RNA synthesis and structure mapping. We show herein that the EteRNA community leveraged dozens of cycles of continuous wet laboratory feedback to learn strategies for solving in vitro RNA design problems on which automated methods fail. The top strategies—including several previously unrecognized negative design rules—were distilled by machine learning into an algorithm, EteRNABot. Over a rigorous 1-y testing phase, both the EteRNA community and EteRNABot significantly outperformed prior algorithms in a dozen RNA secondary structure design tests, including the creation of dendrimer-like structures and scaffolds for small molecule sensors. These results show that an online community can carry out largescale experiments, hypothesis generation, and algorithm design to create practical advances in empirical science. [ABSTRACT FROM AUTHOR]