Search

Modern life is essentially homochiral, containing D-sugars in nucleic acid backbones and L-amino acids in proteins. Since coded proteins are theorized to have developed from a prebiotic RNA World, the homochirality of L-amino acids observed in all known life presumably resulted from chiral transfer from a homochiral D-RNA World. This transfer would have been mediated by aminoacyl-RNAs defining the genetic code. Previous work on aminoacyl transfer using tRNA mimics has suggested that aminoacylation using D-RNA may be inherently biased toward reactivity with L-amino acids, implying a deterministic path from a D-RNA World to L-proteins. Using a model system of self-aminoacylating D-ribozymes and epimerizable activated amino acid analogs, we test the chiral selectivity of 15 ribozymes derived from an exhaustive search of sequence space. All of the ribozymes exhibit detectable selectivity, and a substantial fraction react preferentially to produce the D-enantiomer of the product. Furthermore, chiral preference is conserved within sequence families. These results are consistent with the transfer of chiral information from RNA to proteins but do not support an intrinsic bias of D-RNA for L-amino acids. Different aminoacylation structures result in different directions of chiral selectivity, such that L-proteins need not emerge from a D-RNA World.

RNA structural switches are key regulators of gene expression in bacteria, but their characterization in Metazoa remains limited. Here, we present SwitchSeeker, a comprehensive computational and experimental approach for systematic identification of functional RNA structural switches. We applied SwitchSeeker to the human transcriptome and identified 245 putative RNA switches. To validate our approach, we characterized a previously unknown RNA switch in the 3' untranslated region of the RORC (RAR-related orphan receptor C) transcript. In vivo dimethyl sulfate (DMS) mutational profiling with sequencing (DMS-MaPseq), coupled with cryogenic electron microscopy, confirmed its existence as two alternative structural conformations. Furthermore, we used genome-scale CRISPR screens to identify trans factors that regulate gene expression through this RNA structural switch. We found that nonsense-mediated messenger RNA decay acts on this element in a conformation-specific manner. SwitchSeeker provides an unbiased, experimentally driven method for discovering RNA structural switches that shape the eukaryotic gene expression landscape.

DDX3X regulates the translation of a subset of human transcripts containing complex 5' untranslated regions (5' UTRs). In this study, we developed the helicase activity reporter for translation (HART), which uses DDX3X-sensitive 5' UTRs to measure DDX3X-mediated translational activity in cells. To directly measure RNA structure in DDX3X-dependent mRNAs, we used SHAPE-MaP to determine the secondary structures present in DDX3X-sensitive 5' UTRs and then used HART to investigate how sequence alterations influence DDX3X sensitivity. Additionally, we identified residues 38-44 as potential mediators of DDX3X's interaction with the translational machinery. HART revealed that both DDX3X's association with the translational machinery and its helicase activity are required for its function in promoting the translation of DDX3X-sensitive 5' UTRs. These findings suggest DDX3X plays a crucial role in regulating translation through its interaction with the translational machinery during ribosome scanning and establish the HART reporter as a robust, lentivirally encoded, colorimetric measurement of DDX3X-dependent translation in cells.

Recent discoveries in biology have highlighted the importance of protein and RNA-based condensates as an alternative to classical membrane-bound organelles. Here, we demonstrate the design of pure RNA condensates from nanostructured, star-shaped RNA motifs. We generate condensates using two different RNA nanostar architectures: multi-stranded nanostars whose binding interactions are programmed via linear overhangs, and single-stranded nanostars whose interactions are programmed via kissing loops. Through systematic sequence design, we demonstrate that both architectures can produce orthogonal (distinct and immiscible) condensates, which can be individually tracked via fluorogenic aptamers. We also show that aptamers make it possible to recruit peptides and proteins to the condensates with high specificity. Successful co-transcriptional formation of condensates from single-stranded nanostars suggests that they may be genetically encoded and produced in living cells. We provide a library of orthogonal RNA condensates that can be modularly customized and offer a route toward creating systems of functional artificial organelles for the task of compartmentalizing molecules and biochemical reactions.

The structural diversity of biological macromolecules in different environments contributes complexity to enzymological processes vital for cellular functions. Fluorescence resonance energy transfer and electron microscopy are used to investigate the enzymatic reaction of T4 DNA ligase catalyzing the ligation of nicked DNA. The data show that both the ligase-AMP complex and the ligase-AMP-DNA complex can have four conformations. This finding suggests the parallel occurrence of four ligation reaction pathways, each characterized by specific conformations of the ligase-AMP complex that persist in the ligase-AMP-DNA complex. Notably, these complexes have DNA bending angles of ≈0°, 20°, 60°, or 100°. The mechanism of parallel reactions challenges the conventional notion of simple sequential reaction steps occurring among multiple conformations. The results provide insights into the dynamic conformational changes and the versatile attributes of T4 DNA ligase and suggest that the parallel multiple reaction pathways may correspond to diverse T4 DNA ligase functions. This mechanism may potentially have evolved as an adaptive strategy across evolutionary history to navigate complex environments.

Thermostable clustered regularly interspaced short palindromic repeats (CRISPR) and CRISPR-associated (Cas9) enzymes could improve genome-editing efficiency and delivery due to extended protein lifetimes. However, initial experimentation demonstrated Geobacillus stearothermophilus Cas9 (GeoCas9) to be virtually inactive when used in cultured human cells. Laboratory-evolved variants of GeoCas9 overcome this natural limitation by acquiring mutations in the wedge (WED) domain that produce >100-fold-higher genome-editing levels. Cryoelectron microscopy (cryo-EM) structures of the wild-type and improved GeoCas9 (iGeoCas9) enzymes reveal extended contacts between the WED domain of iGeoCas9 and DNA substrates. Biochemical analysis shows that iGeoCas9 accelerates DNA unwinding to capture substrates under the magnesium-restricted conditions typical of mammalian but not bacterial cells. These findings enabled rational engineering of other Cas9 orthologs to enhance genome-editing levels, pointing to a general strategy for editing enzyme improvement. Together, these results uncover a new role for the Cas9 WED domain in DNA unwinding and demonstrate how accelerated target unwinding dramatically improves Cas9-induced genome-editing activity.

Insertion sequence (IS) elements are the simplest autonomous transposable elements found in prokaryotic genomes1. We recently discovered that IS110 family elements encode a recombinase and a non-coding bridge RNA (bRNA) that confers modular specificity for target DNA and donor DNA through two programmable loops2. Here we report the cryo-electron microscopy structures of the IS110 recombinase in complex with its bRNA, target DNA and donor DNA in three different stages of the recombination reaction cycle. The IS110 synaptic complex comprises two recombinase dimers, one of which houses the target-binding loop of the bRNA and binds to target DNA, whereas the other coordinates the bRNA donor-binding loop and donor DNA. We uncovered the formation of a composite RuvC-Tnp active site that spans the two dimers, positioning the catalytic serine residues adjacent to the recombination sites in both target and donor DNA. A comparison of the three structures revealed that (1) the top strands of target and donor DNA are cleaved at the composite active sites to form covalent 5-phosphoserine intermediates, (2) the cleaved DNA strands are exchanged and religated to create a Holliday junction intermediate, and (3) this intermediate is subsequently resolved by cleavage of the bottom strands. Overall, this study reveals the mechanism by which a bispecific RNA confers target and donor DNA specificity to IS110 recombinases for programmable DNA recombination.

This comprehensive Review delves into the chemical principles governing RNA-mediated crowding events, commonly referred to as granules or biological condensates. We explore the pivotal role played by RNA sequence, structure, and chemical modifications in these processes, uncovering their correlation with crowding phenomena under physiological conditions. Additionally, we investigate instances where crowding deviates from its intended function, leading to pathological consequences. By deepening our understanding of the delicate balance that governs molecular crowding driven by RNA and its implications for cellular homeostasis, we aim to shed light on this intriguing area of research. Our exploration extends to the methodologies employed to decipher the composition and structural intricacies of RNA granules, offering a comprehensive overview of the techniques used to characterize them, including relevant computational approaches. Through two detailed examples highlighting the significance of noncoding RNAs, NEAT1 and XIST, in the formation of phase-separated assemblies and their influence on the cellular landscape, we emphasize their crucial role in cellular organization and function. By elucidating the chemical underpinnings of RNA-mediated molecular crowding, investigating the role of modifications, structures, and composition of RNA granules, and exploring both physiological and aberrant phase separation phenomena, this Review provides a multifaceted understanding of the intriguing world of RNA-mediated biological condensates.

Many experimental and computational efforts have sought to understand DNA origami folding, but the time and length scales of this process pose significant challenges. Here, we present a mesoscopic model that uses a switchable force field to capture the behavior of single- and double-stranded DNA motifs and transitions between them, allowing us to simulate the folding of DNA origami up to several kilobases in size. Brownian dynamics simulations of small structures reveal a hierarchical folding process involving zipping into a partially folded precursor followed by crystallization into the final structure. We elucidate the effects of various design choices on folding order and kinetics. Larger structures are found to exhibit heterogeneous staple incorporation kinetics and frequent trapping in metastable states, as opposed to more accessible structures which exhibit first-order kinetics and virtually defect-free folding. This model opens an avenue to better understand and design DNA nanostructures for improved yield and folding performance.

It is well established that neutrophils adopt malleable polymorphonuclear shapes to migrate through narrow interstitial tissue spaces1-3. However, how polymorphonuclear structures are assembled remains unknown4. Here we show that in neutrophil progenitors, halting loop extrusion-a motor-powered process that generates DNA loops by pulling in chromatin5-leads to the assembly of polymorphonuclear genomes. Specifically, we found that in mononuclear neutrophil progenitors, acute depletion of the loop-extrusion loading factor nipped-B-like protein (NIPBL) induced the assembly of horseshoe, banded, ringed and hypersegmented nuclear structures and led to a reduction in nuclear volume, mirroring what is observed during the differentiation of neutrophils. Depletion of NIPBL also induced cell-cycle arrest, activated a neutrophil-specific gene program and conditioned a loss of interactions across topologically associating domains to generate a chromatin architecture that resembled that of differentiated neutrophils. Removing NIPBL resulted in enrichment for mega-loops and interchromosomal hubs that contain genes associated with neutrophil-specific enhancer repertoires and an inflammatory gene program. On the basis of these observations, we propose that in neutrophil progenitors, loop-extrusion programs produce lineage-specific chromatin architectures that permit the packing of chromosomes into geometrically confined lobular structures. Our data also provide a blueprint for the assembly of polymorphonuclear structures, and point to the possibility of engineering de novo nuclear shapes to facilitate the migration of effector cells in densely populated tumorigenic environments.

The conformational dynamics of nucleosome arrays generate a diverse spectrum of microscopic states, posing challenges to their structural determination. Leveraging cryogenic electron tomography (cryo-ET), we determine the three-dimensional (3D) structures of individual mononucleosomes and arrays comprising di-, tri-, and tetranucleosomes. By slowing the rate of condensation through a reduction in ionic strength, we probe the intra-array structural transitions that precede inter-array interactions and liquid droplet formation. Under these conditions, the arrays exhibite irregular zig-zag conformations with loose packing. Increasing the ionic strength promoted intra-array compaction, yet we do not observe the previously reported regular 30-nanometer fibers. Interestingly, the presence of H1 do not induce array compaction; instead, one-third of the arrays display nucleosomes invaded by foreign DNA, suggesting an alternative role for H1 in chromatin network construction. We also find that the crucial parameter determining the structure adopted by chromatin arrays is the angle between the entry and exit of the DNA and the corresponding tangents to the nucleosomal disc. Our results provide insights into the initial stages of intra-array compaction, a critical precursor to condensation in the regulation of chromatin organization.

Branching is a critical step in RNA splicing that is essential for 5 splice site selection. Recent spliceosome structures have led to competing models for the recognition of the invariant adenosine at the branch point. However, there are no structures of any splicing complex with the adenosine nucleophile docked in the active site and positioned to attack the 5 splice site. Thus we lack a mechanistic understanding of adenosine selection and splice site recognition during RNA splicing. Here we present a cryo-electron microscopy structure of a group II intron that reveals that active site dynamics are coupled to the formation of a base triple within the branch-site helix that positions the 2-OH of the adenosine for nucleophilic attack on the 5 scissile phosphate. This structure, complemented with biochemistry and comparative analyses to splicing complexes, supports a base triple model of adenosine recognition for branching within group II introns and the evolutionarily related spliceosome.

To explore how an early, RNA-based life form could have functioned, in vitro selection experiments have been used to develop catalytic RNAs (ribozymes) with relevant functions. We previously identified ribozymes that use the prebiotically plausible energy source cyclic trimetaphosphate (cTmp) to convert their 5'-hydroxyl group to a 5'-triphosphate. While these ribozymes were developed in the presence of Mg2+, we tested here whether lanthanides could also serve as catalytic cofactors because lanthanides are ideal catalytic cations for this reaction. After an in vitro selection in the presence of Yb3+, several active sequences were isolated, and the most active RNA was analyzed in more detail. This ribozyme required lanthanides for activity, with highest activity at a 10:1 molar ratio of cTmp : Yb3+. Only the four heaviest lanthanides gave detectable signals, indicating a high sensitivity of ribozyme catalysis to the lanthanide ion radius. Potassium and Magnesium did not facilitate catalysis alone but they increased the lanthanide-mediated kOBS by at least 100-fold, with both K+ and Mg2+ modulating the ribozyme's secondary structure. Together, these findings show that RNA is able to use the unique properties of lanthanides as catalytic cofactor. The results are discussed in the context of early life forms.

RNA origami is a method for designing RNA nanostructures that can self-assemble through co-transcriptional folding with applications in nanomedicine and synthetic biology. However, to advance the method further, an improved understanding of RNA structural properties and folding principles is required. Here we use cryogenic electron microscopy to study RNA origami sheets and bundles at sub-nanometre resolution revealing structural parameters of kissing-loop and crossover motifs, which are used to improve designs. In RNA bundle designs, we discover a kinetic folding trap that forms during folding and is only released after 10 h. Exploration of the conformational landscape of several RNA designs reveal the flexibility of helices and structural motifs. Finally, sheets and bundles are combined to construct a multidomain satellite shape, which is characterized by individual-particle cryo-electron tomography to reveal the domain flexibility. Together, the study provides a structural basis for future improvements to the design cycle of genetically encoded RNA nanodevices.

La-related protein 7 (LARP7) are a family of RNA chaperones that protect the 3'-end of RNA and are components of specific ribonucleoprotein complexes (RNP). In Tetrahymena thermophila telomerase, LARP7 protein p65 together with telomerase reverse transcriptase (TERT) and telomerase RNA (TER) form the core RNP. p65 has four known domains-N-terminal domain (NTD), La motif (LaM), RNA recognition motif 1 (RRM1), and C-terminal xRRM2. To date, only the xRRM2 and LaM and their interactions with TER have been structurally characterized. Conformational dynamics leading to low resolution in cryo-EM density maps have limited our understanding of how full-length p65 specifically recognizes and remodels TER for telomerase assembly. Here, we combined focused classification of Tetrahymena telomerase cryo-EM maps with NMR spectroscopy to determine the structure of p65-TER. Three previously unknown helices are identified, one in the otherwise intrinsically disordered NTD that binds the La module, one that extends RRM1, and another preceding xRRM2, that stabilize p65-TER interactions. The extended La module (αN, LaM and RRM1) interacts with the four 3' terminal U nucleotides, while LaM and αN additionally interact with TER pseudoknot, and LaM with stem 1 and 5' end. Our results reveal the extensive p65-TER interactions that promote TER 3'-end protection, TER folding, and core RNP assembly and stabilization. The structure of full-length p65 with TER also sheds light on the biological roles of genuine La and LARP7 proteins as RNA chaperones and core RNP components.

RNA is a key regulator of almost every cellular process, and the structures adopted by RNA molecules are thought to be central to their functions. The recent fast-paced evolution of high-throughput sequencing-based RNA structure mapping methods has enabled the rapid in vivo structural interrogation of entire cellular transcriptomes. Collectively, these studies are shedding new light on the long underestimated complexity of the structural organization of the transcriptome - the RNA structurome. Moreover, recent analyses are challenging the view that the RNA structurome is a static entity by revealing how RNA molecules establish intricate networks of alternative intramolecular and intermolecular interactions and that these ensembles of RNA structures are dynamically regulated to finely tune RNA functions in living cells. This new understanding of how RNA can shape cell phenotypes has important implications for the development of RNA-targeted therapeutic strategies.

R loops (DNA-RNA hybrid) are three-stranded nucleic acid structures that comprise of template DNA strand hybridized with the nascent RNA leaving the displaced non-template strand. Although a programmed R loop formation can serve as powerful regulators of gene expression, these structures can also turn into major sources of genomic instability and contribute to the development of diseases. Therefore, understanding how cells prevent the deleterious consequences of R loops yet allow R loop formation to participate in various physiological processes will help to understand how their homeostasis is maintained. Detection and quantitative measurements of R loops are critical that largely relied on S9.6 antibody. Immunofluorescence methods are frequently used to localize and quantify R loops in the cell but they require specialized tools for analysis and relatively expensive; therefore, they are not always useful for initial assessments of R loop accumulation. Here, we describe an improved slot blot protocol to detect and estimate R loops and show its sensitivity and specificity using the S9.6 antibody. Since specific factors protecting cells from harmful R loop accumulation are expanding, this protocol can be used to determine R loop accumulation in research and clinical settings.

DNA is an attractive material for a range of applications in nanoscience and nanotechnology, and it has recently been demonstrated that the electronic properties of DNA are uniquely sensitive to its sequence and structure, opening new opportunities for the development of electronic DNA biosensors. In this report, we examine the origin of multiple conductance peaks that can occur during single-molecule break-junction (SMBJ)-based conductance measurements on DNA. We demonstrate that these peaks originate from the presence of multiple DNA conformations within the solutions, in particular, double-stranded B-form DNA (dsDNA) and G-quadruplex structures. Using a combination of circular dichroism (CD) spectroscopy, computational approaches, sequence and environmental controls, and single-molecule conductance measurements, we disentangle the conductance information and demonstrate that specific conductance values come from specific conformations of the DNA and that the occurrence of these peaks can be controlled by controlling the local environment. In addition, we demonstrate that conductance measurements are uniquely sensitive to identifying these conformations in solutions and that multiple configurations can be detected in solutions over an extremely large concentration range, opening new possibilities for examining low-probability DNA conformations in solutions.

The Hoogsteen (HG) base pairing conformation, commonly observed in damaged and mutated DNA helices, facilitates DNA repair and DNA recognition. The free energy difference between HG and Watson-Crick (WC) base pairs has been computed in previous studies. However, the mechanism of the conformational transition is not well understood. A detailed understanding of the process of WC to HG base pair transition can provide a deeper understanding of DNA repair and recognition. In an earlier study, we explored the free energy landscape for this process using extensive computer simulation with the CHARMM36 force field. In this work, we study the impact of force field models in describing the WC to HG base pairing transition using meta-eABF enhanced sampling, quasi-harmonic entropy calculation, and nonbonded energy analysis. The secondary structures of both base pairing forms and the topology of the free energy landscapes were consistent over different force field models, although the relative free energy, entropy, and the interaction energies tend to vary. The relative stability of the WC and HG conformations is dictated by a delicate balance between the enthalpic stabilization and the reduced entropy of the structurally rigid HG structure. These findings highlight the impact that subtleties in force field models can have on accurately modeling DNA base pair dynamics and should stimulate further computational investigations into other dynamically important motions in DNA.

ADARs (adenosine deaminases acting on RNA) can be directed to sites in the transcriptome by complementary guide strands allowing for the correction of disease-causing mutations at the RNA level. However, ADARs show bias against editing adenosines with a guanosine 5' nearest neighbor (5'-GA sites), limiting the scope of this approach. Earlier studies suggested this effect arises from a clash in the RNA minor groove involving the 2-amino group of the guanosine adjacent to an editing site. Here we show that nucleosides capable of pairing with guanosine in a syn conformation enhance editing for 5'-GA sites. We describe the crystal structure of a fragment of human ADAR2 bound to RNA bearing a G:G pair adjacent to an editing site. The two guanosines form a Gsyn:Ganti pair solving the steric problem by flipping the 2-amino group of the guanosine adjacent to the editing site into the major groove. Also, duplexes with 2'-deoxyadenosine and 3-deaza-2'-deoxyadenosine displayed increased editing efficiency, suggesting the formation of a Gsyn:AH+anti pair. This was supported by X-ray crystallography of an ADAR complex with RNA bearing a G:3-deaza dA pair. This study shows how non-Watson-Crick pairing in duplex RNA can facilitate ADAR editing enabling the design of next generation guide strands for therapeutic RNA editing.

Xeno-nucleic acids (XNAs) are synthetic genetic polymers with backbone structures composed of non-ribose or non-deoxyribose sugars. Phosphonomethylthreosyl nucleic acid (pTNA), a type of XNA that does not base pair with DNA or RNA, has been suggested as a possible genetic material for storing synthetic biology information in cells. A critical step in this process is the synthesis of XNA episomes using laboratory-evolved polymerases to copy DNA information into XNA. Here, we investigate the polymerase recognition of pTNA nucleotides using X-ray crystallography to capture the post-catalytic complex of engineered polymerases following the sequential addition of two pTNA nucleotides onto the 3'-end of a DNA primer. High-resolution crystal structures reveal that the polymerase mediates Watson-Crick base pairing between the extended pTNA adducts and the DNA template. Comparative analysis studies demonstrate that the sugar conformation and backbone position of pTNA are structurally more similar to threose nucleic acid than DNA even though pTNA and DNA share the same six-atom backbone repeat length. Collectively, these findings provide new insight into the structural determinants that guide the enzymatic synthesis of an orthogonal genetic polymer, and may lead to the discovery of new variants that function with enhanced activity.

Riboswitches are 5'-untranslated regions of mRNA that change their conformation in response to ligand binding, allowing post-transcriptional gene regulation. This ligand-based model of riboswitch function has been expanded with the discovery of a "pH-responsive element" (PRE) riboswitch in Escherichia coli. At neutral pH, the PRE folds into a translationally inactive structure with an occluded ribosome-binding sequence, whereas at alkaline pH, the PRE adopts a translationally active structure. This unique riboswitch does not rely on ligand binding in a traditional sense to modulate its alternative folding outcomes. Rather, pH controls riboswitch folding by two possible modes that are yet to be distinguished; pH either regulates the transcription rate of RNA polymerase (RNAP) or acts on the RNA itself. Previous work suggested that RNAP pausing is prolonged by alkaline pH at two sites, stimulating PRE folding into the active structure. To date, there has been no rigorous exploration into how pH influences RNAP pausing kinetics during PRE synthesis. To provide that understanding and distinguish between pH acting on RNAP versus RNA, we investigated RNAP pausing kinetics at key sites for PRE folding under different pH conditions. We find that pH influences RNAP pausing but not in the manner proposed previously. Rather, alkaline pH either decreases or has no effect on RNAP pause longevity, suggesting that the modulation of RNAP pausing is not the sole mechanism by which pH affects PRE folding. These findings invite the possibility that the RNA itself actively participates in the sensing of pH.

Many large protein-nucleic acid complexes exhibit allosteric regulation. In these systems, the propagation of the allosteric signaling is strongly coupled to conformational dynamics and catalytic function, challenging state-of-the-art analytical methods. Here, we review established and innovative approaches used to elucidate allosteric mechanisms in these complexes. Specifically, we report network models derived from graph theory and centrality analyses in combination with molecular dynamics (MD) simulations, introducing novel schemes that implement the synergistic use of graph theory with enhanced simulations methods and ab-initio MD. Accelerated MD simulations are used to construct "enhanced network models", describing the allosteric response over long timescales and capturing the relation between allostery and conformational changes. "Ab-initio network models" combine graph theory with ab-initio MD and quantum mechanics/molecular mechanics (QM/MM) simulations to describe the allosteric regulation of catalysis by following the step-by-step dynamics of biochemical reactions. This approach characterizes how the allosteric regulation changes from reactants to products and how it affects the transition state, revealing a tense-to-relaxed allosteric regulation along the chemical step. Allosteric models and applications are showcased for three paradigmatic examples of allostery in protein-nucleic acid complexes: (i) the nucleosome core particle, (ii) the CRISPR-Cas9 genome editing system and (iii) the spliceosome. These methods and applications create innovative protocols to determine allosteric mechanisms in protein-nucleic acid complexes that show tremendous promise for medicine and bioengineering.

DNA looping has emerged as a central paradigm of transcriptional regulation, as it is shared across many living systems. One core property of DNA looping-based regulation is its ability to greatly enhance repression or activation of genes with only a few copies of transcriptional regulators. However, this property based on a small number of proteins raises the question of the robustness of such a mechanism with respect to the large intracellular perturbations taking place during growth and division of the cell. Here we address the issue of sensitivity to variations of intracellular parameters of gene regulation by DNA looping. We use the lac system as a prototype to experimentally identify the key features of the robustness of DNA looping in growing Escherichia coli cells. Surprisingly, we observe time intervals of tight repression spanning across division events, which can sometimes exceed 10 generations. Remarkably, the distribution of such long time intervals exhibits memoryless statistics that is mostly insensitive to repressor concentration, cell division events, and the number of distinct loops accessible to the system. By contrast, gene regulation becomes highly sensitive to these perturbations when DNA looping is absent. Using stochastic simulations, we propose that the observed robustness to division emerges from the competition between fast, multiple rebinding events of repressors and slow initiation rate of the RNA polymerase. We argue that fast rebinding events are a direct consequence of DNA looping that ensures robust gene repression across a range of intracellular perturbations.

The ribosomal RNAs, along with their substrates the transfer RNAs, contain the most highly conserved nucleotides in all of biology. We have assembled a database containing structure-based alignments of sequences of the small-subunit rRNAs from organisms that span the entire phylogenetic spectrum, to identify the nucleotides that are universally conserved. In its simplest (bacterial and archaeal) forms, the small-subunit rRNA has ∼1500 nt, of which we identify 140 that are absolutely invariant among the 1961 species in our alignment. We examine the positions and detailed structural and functional interactions of these universal nucleotides in the context of a half century of biochemical and genetic studies and high-resolution structures of ribosome functional complexes. The vast majority of these nucleotides are exposed on the subunit interface surface of the small subunit, where the functional processes of the ribosome take place. However, only 40 of them have been directly implicated in specific ribosomal functions, such as contacting the tRNAs, mRNA, or translation factors. The roles of many other invariant nucleotides may serve to constrain the positions and orientations of those nucleotides that are directly involved in function. Yet others can be rationalized by participation in unusual noncanonical tertiary structures that may uniquely allow correct folding of the rRNA to form a functional ribosome. However, there remain at least 50 nt whose universal conservation is not obvious, serving as a metric for the incompleteness of our understanding of ribosome structure and function.

Metazoan organisms have many tRNA genes responsible for decoding amino acids. The set of all tRNA genes can be grouped in sets of common amino acids and isoacceptor tRNAs that are aminoacylated by corresponding aminoacyl-tRNA synthetases. Analysis of tRNA alignments shows that, despite the high number of tRNA genes, specific tRNA sequence motifs are highly conserved across multicellular eukaryotes. The conservation often extends throughout the isoacceptors and isodecoders with, in some cases, two sets of conserved isodecoders. This study is focused on non-Watson-Crick base pairs in the helical stems, especially GoU pairs. Each of the four helical stems may contain one or more conserved GoU pairs. Some are amino acid specific and could represent identity elements for the cognate aminoacyl tRNA synthetases. Other GoU pairs are found in more than a single amino acid and could be critical for native folding of the tRNAs. Interestingly, some GoU pairs are anticodon-specific, and others are found in phylogenetically-specific clades. Although the distribution of conservation likely reflects a balance between accommodating isotype-specific functions as well as those shared by all tRNAs essential for ribosomal translation, such conservations may indicate the existence of specialized tRNAs for specific translation targets, cellular conditions, or alternative functions.

SignificanceUnderstanding and treating neurological disorders are global priorities. Some of these diseases are engendered by mutations that cause defects in the cellular synthesis of transfer RNAs (tRNAs), which function as adapter molecules that translate messenger RNAs into proteins. During tRNA biogenesis, ribonuclease P catalyzes removal of the transcribed sequence upstream of the mature tRNA. Here, we focus on a cytoplasmic tRNAArgUCU that is expressed specifically in neurons and, when harboring a particular point mutation, contributes to neurodegeneration in mice. Our results suggest that this mutation favors stable alternative structures that are not cleaved by mouse ribonuclease P and motivate a paradigm that may help to understand the molecular basis for disease-associated mutations in other tRNAs.

High-fidelity replication of the large RNA genome of coronaviruses (CoVs) is mediated by a 3'-to-5' exoribonuclease (ExoN) in nonstructural protein 14 (nsp14), which excises nucleotides including antiviral drugs misincorporated by the low-fidelity viral RNA-dependent RNA polymerase (RdRp) and has also been implicated in viral RNA recombination and resistance to innate immunity. Here, we determined a 1.6-Å resolution crystal structure of severe acute respiratory syndrome CoV 2 (SARS-CoV-2) ExoN in complex with its essential cofactor, nsp10. The structure shows a highly basic and concave surface flanking the active site, comprising several Lys residues of nsp14 and the N-terminal amino group of nsp10. Modeling suggests that this basic patch binds to the template strand of double-stranded RNA substrates to position the 3' end of the nascent strand in the ExoN active site, which is corroborated by mutational and computational analyses. We also show that the ExoN activity can rescue a stalled RNA primer poisoned with sofosbuvir and allow RdRp to continue its extension in the presence of the chain-terminating drug, biochemically recapitulating proofreading in SARS-CoV-2 replication. Molecular dynamics simulations further show remarkable flexibility of multidomain nsp14 and suggest that nsp10 stabilizes ExoN for substrate RNA binding to support its exonuclease activity. Our high-resolution structure of the SARS-CoV-2 ExoN-nsp10 complex serves as a platform for future development of anticoronaviral drugs or strategies to attenuate the viral virulence.

For many RNA molecules, the secondary structure is essential for the correct function of the RNA. Predicting RNA secondary structure from nucleotide sequences is a long-standing problem in genomics, but the prediction performance has reached a plateau over time. Traditional RNA secondary structure prediction algorithms are primarily based on thermodynamic models through free energy minimization, which imposes strong prior assumptions and is slow to run. Here, we propose a deep learning-based method, called UFold, for RNA secondary structure prediction, trained directly on annotated data and base-pairing rules. UFold proposes a novel image-like representation of RNA sequences, which can be efficiently processed by Fully Convolutional Networks (FCNs). We benchmark the performance of UFold on both within- and cross-family RNA datasets. It significantly outperforms previous methods on within-family datasets, while achieving a similar performance as the traditional methods when trained and tested on distinct RNA families. UFold is also able to predict pseudoknots accurately. Its prediction is fast with an inference time of about 160 ms per sequence up to 1500 bp in length. An online web server running UFold is available at https://ufold.ics.uci.edu. Code is available at https://github.com/uci-cbcl/UFold.

The binding affinity of therapeutic oligonucleotides (ONs) for their cognate RNA is determined by the rates of association (ka) and dissociation (kd). Single-stranded ONs are highly flexible and can adopt multiple conformations in solution, some of which may not be conducive for hybridization. We investigated if restricting rotation around the sugar-phosphate backbone, by tethering two adjacent backbone phosphonate esters using hydrocarbon bridges, can modulate hybridization kinetics of the modified ONs for complementary RNA. Given the large number of possible analogues with different tether lengths and configurations at the phosphorus atoms, we employed molecular dynamic simulations to optimize the size of the hydrocarbon bridge to guide the synthetic efforts. The backbone-constrained nucleotide trimers with stereodefined configurations at the contiguous backbone phosphorus atoms were assembled using a ring-closing metathesis reaction, then incorporated into oligonucleotides by an in situ synthesis of the phosphoramidites followed by coupling to solid supports. Evaluation of the modified oligonucleotides revealed that 15-membered macrocyclic-constrained analogues displayed similar or slightly improved on-rates but significantly increased off-rates compared to unmodified DNA ONs, resulting in reduced duplex stability. In contrast, LNA ONs with conformationally preorganized furanose rings showed similar on-rates to DNA ONs but very slow off-rates, resulting in net improvement in duplex stability. Furthermore, the experimental data generally supported the molecular dynamics simulation results, suggesting that this strategy can be used as a predictive tool for designing the next generation of constrained backbone ON analogues with improved hybridization properties.

DNA supercoils are generated in genome regulation processes such as transcription and replication and provide mechanical feedback to such processes. Under tension, a DNA supercoil can present a coexistence state of plectonemic and stretched phases. Experiments have revealed the dynamic behaviors of plectonemes, e.g., diffusion, nucleation, and hopping. To represent these dynamics with conformational changes, we demonstrated first the fast dynamics on the DNA to reach torque equilibrium within the plectonemic and stretched phases, and then identified the two-phase boundaries as collective slow variables to describe the essential dynamics. According to the timescale separation demonstrated here, we developed a two-phase model on the dynamics of DNA supercoiling, which can capture physiologically relevant events across timescales of several orders of magnitudes. In this model, we systematically characterized the slow dynamics between the two phases and compared the numerical results with those from the DNA polymer physics-based worm-like chain model. The supercoiling dynamics, including the nucleation, diffusion, and hopping of plectonemes, have been well represented and reproduced, using the two-phase dynamic model, at trivial computational costs. Our current developments, therefore, can be implemented to explore multiscale physical mechanisms of the DNA supercoiling-dependent physiological processes.

Elongating RNA polymerase II (Pol II) can be paused or arrested by a variety of obstacles. These obstacles include DNA lesions, DNA-binding proteins, and small molecules. Hairpin pyrrole-imidazole (Py-Im) polyamides bind to the minor groove of DNA in a sequence-specific manner and induce strong transcriptional arrest. Remarkably, this Py-Im-induced Pol II transcriptional arrest is persistent and cannot be rescued by transcription factor TFIIS. In contrast, TFIIS can effectively rescue the transcriptional arrest induced by a nucleosome barrier. The structural basis of Py-Im-induced transcriptional arrest and why TFIIS cannot rescue this arrest remain elusive. Here we determined the X-ray crystal structures of four distinct Pol II elongation complexes (Pol II ECs) in complex with hairpin Py-Im polyamides as well as of the hairpin Py-Im polyamides-dsDNA complex. We observed that the Py-Im oligomer directly interacts with RNA Pol II residues, introduces compression of the downstream DNA duplex, prevents Pol II forward translocation, and induces Pol II backtracking. These results, together with biochemical studies, provide structural insight into the molecular mechanism by which Py-Im blocks transcription. Our structural study reveals why TFIIS fails to promote Pol II bypass of Py-Im-induced transcriptional arrest.

Riboswitches are structured non-coding RNAs often located upstream of essential genes in bacterial messenger RNAs. Such RNAs regulate expression of downstream genes by recognizing a specific cellular effector. Although nearly 50 riboswitch classes are known, only a handful recognize multiple effectors. Here, we report the 2.60-Å resolution co-crystal structure of a class I type I preQ1-sensing riboswitch that reveals two effectors stacked atop one another in a single binding pocket. These effectors bind with positive cooperativity in vitro and both molecules are necessary for gene regulation in bacterial cells. Stacked effector recognition appears to be a hallmark of the largest subgroup of preQ1 riboswitches, including those from pathogens such as Neisseria gonorrhoeae. We postulate that binding to stacked effectors arose in the RNA World to closely position two substrates for RNA-mediated catalysis. These findings expand known effector recognition capabilities of riboswitches and have implications for antimicrobial development.

During the origin of life, encapsulation of RNA inside vesicles is believed to have been a defining feature of the earliest cells (protocells). The confined biophysical environment provided by membrane encapsulation differs from that of bulk solution and has been shown to increase activity as well as evolutionary rate for functional RNA. However, the structural basis of the effect on RNA has not been clear. Here, we studied how encapsulation of the hairpin ribozyme inside model protocells affects ribozyme kinetics, ribozyme folding into the active conformation, and cleavage and ligation activities. We further examined the effect of encapsulation on the folding of a stem-loop RNA structure and on the formation of a triplex structure in a pH-sensitive DNA switch. The results indicate that encapsulation promotes RNA-RNA association, both intermolecular and intramolecular, and also stabilizes tertiary folding, including the docked conformation characteristic of the active hairpin ribozyme and the triplex structure. The effects of encapsulation were sufficient to rescue the activity of folding-deficient mutants of the hairpin ribozyme. Stabilization of multiple modes of nucleic acid folding and interaction thus enhanced the activity of encapsulated nucleic acids. Increased association between RNA molecules may facilitate the formation of more complex structures and cooperative interactions. These effects could promote the emergence of biological functions in an "RNA world" and may have utility in the construction of minimal synthetic cells.

Recurrent outbreaks of novel zoonotic coronavirus (CoV) diseases in recent years have highlighted the importance of developing therapeutics with broad-spectrum activity against CoVs. Because all CoVs use -1 programmed ribosomal frameshifting (-1 PRF) to control expression of key viral proteins, the frameshift signal in viral mRNA that stimulates -1 PRF provides a promising potential target for such therapeutics. To test the viability of this strategy, we explored whether small-molecule inhibitors of -1 PRF in SARS-CoV-2 also inhibited -1 PRF in a range of bat CoVs-the most likely source of future zoonoses. Six inhibitors identified in new and previous screens against SARS-CoV-2 were evaluated against the frameshift signals from a panel of representative bat CoVs as well as MERS-CoV. Some drugs had strong activity against subsets of these CoV-derived frameshift signals, while having limited to no effect on -1 PRF caused by frameshift signals from other viruses used as negative controls. Notably, the serine protease inhibitor nafamostat suppressed -1 PRF significantly for multiple CoV-derived frameshift signals. These results suggest it is possible to find small-molecule ligands that inhibit -1 PRF specifically in a broad spectrum of CoVs, establishing frameshift signals as a viable target for developing pan-coronaviral therapeutics.

Systems of catalytic RNAs presumably gave rise to important evolutionary innovations, such as the genetic code. Such systems may exhibit particular tolerance to errors (error minimization) as well as coding specificity. While often assumed to result from natural selection, error minimization may instead be an emergent by-product. In an RNA world, a system of self-aminoacylating ribozymes could enforce the mapping of amino acids to anticodons. We measured the activity of thousands of ribozyme mutants on alternative substrates (activated analogs for tryptophan, phenylalanine, leucine, isoleucine, valine, and methionine). Related ribozymes exhibited shared preferences for substrates, indicating that adoption of additional amino acids by existing ribozymes would itself lead to error minimization. Furthermore, ribozyme activity was positively correlated with specificity, indicating that selection for increased activity would also lead to increased specificity. These results demonstrate that by-products of ribozyme evolution could lead to adaptive value in specificity and error tolerance.

Short segments of RNA displace one strand of a DNA duplex during diverse processes including transcription and CRISPR-mediated immunity and genome editing. These strand exchange events involve the intersection of two geometrically distinct helix types-an RNA:DNA hybrid (A-form) and a DNA:DNA homoduplex (B-form). Although previous evidence suggests that these two helices can stack on each other, it is unknown what local geometric adjustments could enable A-on-B stacking. Here we report the X-ray crystal structure of an RNA-5'/DNA-3' strand exchange junction at an anisotropic resolution of 1.6 to 2.2 Å. The structure reveals that the A-to-B helical transition involves a combination of helical axis misalignment, helical axis tilting and compression of the DNA strand within the RNA:DNA helix, where nucleotides exhibit a mixture of A- and B-form geometry. These structural principles explain previous observations of conformational stability in RNA/DNA exchange junctions, enabling a nucleic acid architecture that is repeatedly populated during biological strand exchange events.

The chemical stage of the Monte Carlo track-structure (MCTS) code Geant4-DNA was extended for its use in DNA strand break (SB) simulations and compared against published experimental data. Geant4-DNA simulations were performed using pUC19 plasmids (2686 base pairs) in a buffered solution of DMSO irradiated by60Co or137Csγ-rays. A comprehensive evaluation of SSB yields was performed considering DMSO, DNA concentration, dose and plasmid supercoiling. The latter was measured using the super helix density value used in a Brownian dynamics plasmid generation algorithm. The Geant4-DNA implementation of the independent reaction times method (IRT), developed to simulate the reaction kinetics of radiochemical species, allowed to score the fraction of supercoiled, relaxed and linearized plasmid fractions as a function of the absorbed dose. The percentage of the number of SB after •OH + DNA and H• + DNA reactions, referred as SSB efficiency, obtained using MCTS were 13.77% and 0.74% respectively. This is in reasonable agreement with published values of 12% and 0.8%. The SSB yields as a function of DMSO concentration, DNA concentration and super helix density recreated the expected published experimental behaviors within 5%, one standard deviation. The dose response of SSB and DSB yields agreed with published measurements within 5%, one standard deviation. We demonstrated that the developed extension of IRT in Geant4-DNA, facilitated the reproduction of experimental conditions. Furthermore, its calculations were strongly in agreement with experimental data. These two facts will facilitate the use of this extension in future radiobiological applications, aiding the study of DNA damage mechanisms with a high level of detail.

During spliceosome assembly, interactions that bring the 5' and 3' ends of an intron in proximity are critical for the production of mature mRNA. Here, we report synergistic roles for the stem-loops 3 (SL3) and 4 (SL4) of the human U1 small nuclear RNA (snRNA) in maintaining the optimal U1 snRNP function, and formation of cross-intron contact with the U2 snRNP. We find that SL3 and SL4 bind distinct spliceosomal proteins and combining a U1 snRNA activity assay with siRNA-mediated knockdown, we demonstrate that SL3 and SL4 act through the RNA helicase UAP56 and the U2 protein SF3A1, respectively. In vitro analysis using UV crosslinking and splicing assays indicated that SL3 likely promotes the SL4-SF3A1 interaction leading to enhancement of A complex formation and pre-mRNA splicing. Overall, these results highlight the vital role of the distinct contacts of SL3 and SL4 in bridging the pre-mRNA bound U1 and U2 snRNPs during the early steps of human spliceosome assembly.

RNA molecules are known to fold into specific structures which often play a central role in their functions and regulation. In silico folding of RNA transcripts, especially when assisted with structure profiling (SP) data, is capable of accurately elucidating relevant structural conformations. However, such methods scale poorly to the swaths of SP data generated by transcriptome-wide experiments, which are becoming more commonplace and advancing our understanding of RNA structure and its regulation at global and local levels. This has created a need for tools capable of rapidly deriving structural assessments from SP data in a scalable manner. One such tool we previously introduced that aims to process such data is patteRNA, a statistical learning algorithm capable of rapidly mining big SP datasets for structural elements. Here, we present a reformulation of patteRNA's pattern recognition scheme that sees significantly improved precision without major compromises to computational overhead. Specifically, we developed a data-driven logistic classifier which interprets patteRNA's statistical characterizations of SP data in addition to local sequence properties as measured with a nearest neighbour thermodynamic model. Application of the classifier to human structurome data reveals a marked association between detected stem-loops and RNA binding protein (RBP) footprints. The results of our application demonstrate that upwards of 30% of RBP footprints occur within loops of stable stem-loop elements. Overall, our work arrives at a rapid and accurate method for automatically detecting families of RNA structure motifs and demonstrates the functional relevance of identifying them transcriptome-wide.

Design and preparation of layered composite materials alternating between nucleic acids and proteins has been elusive due to limitations in occurrence and geometry of interaction sites in natural biomolecules. We report the design and kinetically controlled stepwise synthesis of a nano-sandwich composite by programmed noncovalent association of protein, DNA and RNA modules. A homo-tetramer protein core was introduced to control the self-assembly and precise positioning of two RNA-DNA hybrid nanotriangles in a co-parallel sandwich arrangement. Kinetically favored self-assembly of the circularly closed nanostructures at the protein was driven by the intrinsic fast folding ability of RNA corner modules which were added to precursor complex of DNA bound to the protein. The 3D architecture of this first synthetic protein-RNA-DNA complex was confirmed by fluorescence labeling and cryo-electron microscopy studies. The synthesis strategy for the nano-sandwich composite provides a general blueprint for controlled noncovalent assembly of complex supramolecular architectures from protein, DNA and RNA components, which expand the design repertoire for bottom-up preparation of layered biomaterials.

Stochastic fluctuations in gene expression ("noise") are often considered detrimental, but fluctuations can also be exploited for benefit (e.g., dither). We show here that DNA base excision repair amplifies transcriptional noise to facilitate cellular reprogramming. Specifically, the DNA repair protein Apex1, which recognizes both naturally occurring and unnatural base modifications, amplifies expression noise while homeostatically maintaining mean expression levels. This amplified expression noise originates from shorter-duration, higher-intensity transcriptional bursts generated by Apex1-mediated DNA supercoiling. The remodeling of DNA topology first impedes and then accelerates transcription to maintain mean levels. This mechanism, which we refer to as "discordant transcription through repair" ("DiThR," which is pronounced "dither"), potentiates cellular reprogramming and differentiation. Our study reveals a potential functional role for transcriptional fluctuations mediated by DNA base modifications in embryonic development and disease.

Bacteriophages densely pack their long double-stranded DNA genome inside a protein capsid. The conformation of the viral genome inside the capsid is consistent with a hexagonal liquid crystalline structure. Experiments have confirmed that the details of the hexagonal packing depend on the electrochemistry of the capsid and its environment. In this work, we propose a biophysical model that quantifies the relationship between DNA configurations inside bacteriophage capsids and the types and concentrations of ions present in a biological system. We introduce an expression for the free energy that combines the electrostatic energy with contributions from bending of individual segments of DNA and Lennard-Jones-type interactions between these segments. The equilibrium points of this energy solve a partial differential equation that defines the distributions of DNA and the ions inside the capsid. We develop a computational approach that allows us to simulate much larger systems than what is possible using the existing molecular-level methods. In particular, we are able to estimate bending and repulsion between the DNA segments as well as the full electrochemistry of the solution, both inside and outside of the capsid. The numerical results show good agreement with existing experiments and with molecular dynamics simulations for small capsids.

Intron selection during the formation of prespliceosomes is a critical event in pre-mRNA splicing. Chemical modulation of intron selection has emerged as a route for cancer therapy. Splicing modulators alter the splicing patterns in cells by binding to the U2 snRNP (small nuclear ribonucleoprotein)-a complex chaperoning the selection of branch and 3' splice sites. Here we report crystal structures of the SF3B module of the U2 snRNP in complex with spliceostatin and sudemycin FR901464 analogs, and the cryo-electron microscopy structure of a cross-exon prespliceosome-like complex arrested with spliceostatin A. The structures reveal how modulators inactivate the branch site in a sequence-dependent manner and stall an E-to-A prespliceosome intermediate by covalent coupling to a nucleophilic zinc finger belonging to the SF3B subunit PHF5A. These findings support a mechanism of intron recognition by the U2 snRNP as a toehold-mediated strand invasion and advance an unanticipated drug targeting concept.

Telomerase adds telomeric repeats at chromosome ends to compensate for the telomere loss that is caused by incomplete genome end replication1. In humans, telomerase is upregulated during embryogenesis and in cancers, and mutations that compromise the function of telomerase result in disease2. A previous structure of human telomerase at a resolution of 8 Å revealed a vertebrate-specific composition and architecture3, comprising a catalytic core that is flexibly tethered to an H and ACA (hereafter, H/ACA) box ribonucleoprotein (RNP) lobe by telomerase RNA. High-resolution structural information is necessary to develop treatments that can effectively modulate telomerase activity as a therapeutic approach against cancers and disease. Here we used cryo-electron microscopy to determine the structure of human telomerase holoenzyme bound to telomeric DNA at sub-4 Å resolution, which reveals crucial DNA- and RNA-binding interfaces in the active site of telomerase as well as the locations of mutations that alter telomerase activity. We identified a histone H2A-H2B dimer within the holoenzyme that was bound to an essential telomerase RNA motif, which suggests a role for histones in the folding and function of telomerase RNA. Furthermore, this structure of a eukaryotic H/ACA RNP reveals the molecular recognition of conserved RNA and protein motifs, as well as interactions that are crucial for understanding the molecular pathology of many mutations that cause disease. Our findings provide the structural details of the assembly and active site of human telomerase, which paves the way for the development of therapeutic agents that target this enzyme.

Aberrant alternative splicing is a hallmark of cancer, yet the underlying regulatory programs that control this process remain largely unknown. Here, we report a systematic effort to decipher the RNA structural code that shapes pathological splicing during breast cancer metastasis. We discovered a previously unknown structural splicing enhancer that is enriched near cassette exons with increased inclusion in highly metastatic cells. We show that the spliceosomal protein small nuclear ribonucleoprotein polypeptide A' (SNRPA1) interacts with these enhancers to promote cassette exon inclusion. This interaction enhances metastatic lung colonization and cancer cell invasion, in part through SNRPA1-mediated regulation of PLEC alternative splicing, which can be counteracted by splicing modulating morpholinos. Our findings establish a noncanonical regulatory role for SNRPA1 as a prometastatic splicing enhancer in breast cancer.

The mechanism of DNA synthesis has been inferred from static structures, but the absence of temporal information raises longstanding questions about the order of events in one of life's most central processes. Here we follow the reaction pathway of a replicative DNA polymerase using time-resolved X-ray crystallography to elucidate the order and transition between intermediates. In contrast to the canonical model, the structural changes observed in the time-lapsed images reveal a catalytic cycle in which translocation precedes catalysis. The translocation step appears to follow a push-pull mechanism where the O-O1 loop of the finger subdomain acts as a pawl to facilitate unidirectional movement along the template with conserved tyrosine residues 714 and 719 functioning as tandem gatekeepers of DNA synthesis. The structures capture the precise order of critical events that may be a general feature of enzymatic catalysis among replicative DNA polymerases.

Since the arrival of DNA nanotechnology nearly 40 years ago, the field has progressed from its beginnings of envisioning rather simple DNA structures having a branched, multi-strand architecture into creating beautifully complex structures comprising hundreds or even thousands of unique strands, with the possibility to exactly control the positions down to the molecular level. While the earliest construction methodologies, such as simple Holliday junctions or tiles, could reasonably be designed on pen and paper in a short amount of time, the advent of complex techniques, such as DNA origami or DNA bricks, require software to reduce the time required and propensity for human error within the design process. Where available, readily accessible design software catalyzes our ability to bring techniques to researchers in diverse fields and it has helped to speed the penetration of methods, such as DNA origami, into a wide range of applications from biomedicine to photonics. Here, we review the historical and current state of CAD software to enable a variety of methods that are fundamental to using structural DNA technology. Beginning with the first tools for predicting sequence-based secondary structure of nucleotides, we trace the development and significance of different software packages to the current state-of-the-art, with a particular focus on programs that are open source.

RNA molecules can fold into complex two- and three-dimensional shapes that are critical for their function. Chemical probes have long been utilized to interrogate RNA structure and are now considered invaluable resources in the goal of relating structure to function. Recently, the power of deep sequencing and careful chemical probe design have merged, permitting researchers to obtain a holistic understanding of how RNA structure can be utilized to control RNA biology transcriptome-wide. Within this review, we outline the recent advancements in chemical probe design for interrogating RNA structures inside cells and discuss the recent advances in our understanding of RNA biology through the lens of chemical probing.

We present a new class of DNA-based nanoswitches that, upon enzymatic repair, could undergo a conformational change mechanism leading to a change in fluorescent signal. Such folding-upon-repair DNA nanoswitches are synthetic DNA sequences containing O6 -methyl-guanine (O6 -MeG) nucleobases and labelled with a fluorophore/quencher optical pair. The nanoswitches are rationally designed so that only upon enzymatic demethylation of the O6 -MeG nucleobases they can form stable intramolecular Hoogsteen interactions and fold into an optically active triplex DNA structure. We have first characterized the folding mechanism induced by the enzymatic repair activity through fluorescent experiments and Molecular Dynamics simulations. We then demonstrated that the folding-upon-repair DNA nanoswitches are suitable and specific substrates for different methyltransferase enzymes including the human homologue (hMGMT) and they allow the screening of novel potential methyltransferase inhibitors.