Influence of dimethylsulfoxide on RNA structure and ligand binding

Parageobacillus thermoglucosidasius Strain Engineering Using a Theophylline Responsive RiboCas for Controlled Gene Expression

The relentless increase in atmospheric greenhouse gas concentrations as a consequence of the exploitation of fossil resources compels the adoption of sustainable routes to chemical and fuel manufacture based on biological fermentation processes. The use of thermophilic chassis in such processes is an attractive proposition; however, their effective exploitation will require improved genome editing tools. In the case of the industrially relevant chassis Parageobacillus thermoglucosidasius, CRISPR/Cas9-based gene editing has been demonstrated. The constitutive promoter used, however, accentuates the deleterious nature of Cas9, causing decreased transformation and low editing efficiencies, together with an increased likelihood of off-target effects or alternative mutations. Here, we rectify this issue by controlling the expression of Cas9 through the use of a synthetic riboswitch that is dependent on the nonmetabolized, nontoxic, and cheap inducer, theophylline. We demonstrate that the riboswitches are dose-dependent, allowing for controlled expression of the target gene. Through their use, we were then able to address the deleterious nature of Cas9 and produce an inducible system, RiboCas93. The benefits of RiboCas93 were demonstrated by increased transformation efficiency of the editing vectors, improved efficiency in mutant generation (100%), and a reduction of Cas9 toxicity, as indicated by a reduction in the number of single nucleotide polymorphisms (SNPs) observed. This new system provides a quick and efficient way to produce mutants in P. thermoglucosidasius.

Fluorescent Ligand Equilibrium Displacement: A High-Throughput Method for Identification of FMN Riboswitch-Binding Small Molecules

Antibiotic resistance remains a pressing global concern, with most antibiotics targeting the bacterial ribosome or a limited range of proteins. One class of underexplored antibiotic targets is bacterial riboswitches, structured RNA elements that regulate key biosynthetic pathways by binding a specific ligand. We developed a methodology termed Fluorescent Ligand Equilibrium Displacement (FLED) to rapidly discover small molecules that bind the flavin mononucleotide (FMN) riboswitch. FLED leverages intrinsically fluorescent FMN and the quenching effect on RNA binding to create a label-free, in vitro method to identify compounds that can bind the apo population of riboswitch in a system at equilibrium. The response difference between known riboswitch ligands and controls demonstrates the robustness of the method for high-throughput screening. An existing drug discovery library that was screened using FLED resulted in a final hit rate of 0.67%. The concentration response of each hit was determined and revealed a variety of approximate effective concentration values. Our preliminary screening data support the use of FLED to identify small molecules for medicinal chemistry development as FMN riboswitch-targeted antibiotic compounds. This robust, label-free, and cell-free method offers a strong alternative to other riboswitch screening methods and can be adapted to a variety of laboratory setups.

High-throughput, fluorescent-aptamer-based measurements of steady-state transcription rates for Mycobacterium tuberculosis RNA polymerase

The first step in gene expression is the transcription of DNA sequences into RNA. Regulation at the level of transcription leads to changes in steady-state concentrations of RNA transcripts, affecting the flux of downstream functions and ultimately cellular phenotypes. Changes in transcript levels are routinely followed in cellular contexts via genome-wide sequencing techniques. However, in vitro mechanistic studies of transcription have lagged with respect to throughput. Here, we describe the use of a real-time, fluorescent-aptamer-based method to quantitate steady-state transcription rates of the Mycobacterium tuberculosis RNA polymerase. We present clear controls to show that the assay specifically reports on promoter-dependent, full-length RNA transcription rates that are in good agreement with the kinetics determined by gel-resolved, α- 32 P NTP incorporation experiments. We illustrate how the time-dependent changes in fluorescence can be used to measure regulatory effects of nucleotide concentrations and identity, RNAP and DNA concentrations, transcription factors, and antibiotics. Our data showcase the ability to easily perform hundreds of parallel steady-state measurements across varying conditions with high precision and reproducibility to facilitate the study of the molecular mechanisms of bacterial transcription.

Significance Statement

RNA polymerase transcription mechanisms have largely been determined from in vitro kinetic and structural biology methods. In contrast to the limited throughput of these approaches, in vivo RNA sequencing provides genome-wide measurements but lacks the ability to dissect direct biochemical from indirect genetic mechanisms. Here, we present a method that bridges this gap, permitting high-throughput fluorescence-based measurements of in vitro steady-state transcription kinetics. We illustrate how an RNA-aptamer-based detection system can be used to generate quantitative information on direct mechanisms of transcriptional regulation and discuss the far-reaching implications for future applications.

RNA:RNA interaction in ternary complexes resolved by chemical probing

RNA regulation can be performed by a second targeting RNA molecule, such as in the microRNA regulation mechanism. Selective 2'-hydroxyl acylation analyzed by primer extension (SHAPE) probes the structure of RNA molecules and can resolve RNA:protein interactions, but RNA:RNA interactions have not yet been addressed with this technique. Here, we apply SHAPE to investigate RNA-mediated binding processes in RNA:RNA and RNA:RNA-RBP complexes. We use RNA:RNA binding by SHAPE (RABS) to investigate microRNA-34a (miR-34a) binding its mRNA target, the silent information regulator 1 (mSIRT1), both with and without the Argonaute protein, constituting the RNA-induced silencing complex (RISC). We show that the seed of the mRNA target must be bound to the microRNA loaded into RISC to enable further binding of the compensatory region by RISC, while the naked miR-34a is able to bind the compensatory region without seed interaction. The method presented here provides complementary structural evidence for the commonly performed luciferase-assay-based evaluation of microRNA binding-site efficiency and specificity on the mRNA target site and could therefore be used in conjunction with it. The method can be applied to any nucleic acid-mediated RNA- or RBP-binding process, such as splicing, antisense RNA binding, or regulation by RISC, providing important insight into the targeted RNA structure.

A toolkit for the identification of NEAT1_2/paraspeckle modulators

Paraspeckles are ribonucleoprotein granules assembled by NEAT1_2 lncRNA, an isoform of Nuclear Paraspeckle Assembly Transcript 1 (NEAT1). Dysregulation of NEAT1_2/paraspeckles has been linked to multiple human diseases making them an attractive drug target. However currently NEAT1_2/paraspeckle-focused translational research and drug discovery are hindered by a limited toolkit. To fill this gap, we developed and validated a set of tools for the identification of NEAT1_2 binders and modulators comprised of biochemical and cell-based assays. The NEAT1_2 triple helix stability element was utilized as the target in the biochemical assays, and the cellular assay ('ParaQuant') was based on high-content imaging of NEAT1_2 in fixed cells. As a proof of principle, these assays were used to screen a 1,200-compound FDA-approved drug library and a 170-compound kinase inhibitor library and to confirm the screening hits. The assays are simple to establish, use only commercially-available reagents and are scalable for higher throughput. In particular, ParaQuant is a cost-efficient assay suitable for any cells growing in adherent culture and amenable to multiplexing. Using ParaQuant, we identified dual PI3K/mTOR inhibitors as potent negative modulators of paraspeckles. The tools we describe herein should boost paraspeckle studies and help guide the search, validation and optimization of NEAT1_2/paraspeckle-targeted small molecules.

Binding Properties of RNA Quadruplex of SARS-CoV-2 to Berberine Compared to Telomeric DNA Quadruplex

Previous studies suggest that berberine, an isoquinoline alkaloid, has antiviral potential and is a possible therapeutic candidate against SARS-CoV-2. The molecular underpinnings of its action are still unknown. Potential targets include quadruplexes (G4Q) in the viral genome as they play a key role in modulating the biological activity of viruses. While several DNA-G4Q structures and their binding properties have been elucidated, RNA-G4Qs such as RG-1 of the N-gene of SARS-CoV-2 are less explored. Using biophysical techniques, the berberine binding thermodynamics and the associated conformational and hydration changes of RG-1 could be characterized and compared with human telomeric DNA-G4Q 22AG. Berberine can interact with both quadruplexes. Substantial changes were observed in the interaction of berberine with 22AG and RG-1, which adopt different topologies that can also change upon ligand binding. The strength of interaction and the thermodynamic signatures were found to dependent not only on the initial conformation of the quadruplex, but also on the type of salt present in solution. Since berberine has shown promise as a G-quadruplex stabilizer that can modulate viral gene expression, this study may also contribute to the development of optimized ligands that can discriminate between binding to DNA and RNA G-quadruplexes.

Using DMS-MaPseq to Uncover the Roles of DEAD-box Proteins in Ribosome Assembly

DMS (dimethylsulfate) is a time-tested chemical probe for nucleic acid secondary structure that has recently re-emerged as a powerful tool to study RNA structure and structural changes, by coupling it to high throughput sequencing techniques. This variant, termed DMS-MaPseq, allows for mapping of all RNAs in a cell at the same time. However, if an RNA adopts different structures, for example during the assembly of an RNA-protein complex, or as part of its functional cycle, then DMS-MaPseq cannot differentiate between these structures, and an ensemble average will be produced. This is especially challenging for long-lived RNAs, such as ribosomes, whose steady-state abundance far exceeds that of any assembly intermediates, rendering those inaccessible to DMS-MaPseq on total RNAs. These challenges can be overcome by purification of assembly intermediates stalled at specific assembly steps (or steps in the functional cycle), via a combination of affinity tags and mutants stalled at defined steps, and subsequent DMS probing of these intermediates. Interpretation of the differences in DMS accessibility is facilitated by additional structural information, e.g. from cryo-EM experiments, available for many functional RNAs. While this approach is generally useful for studying RNA folding or conformational changes within RNA-protein complexes, it can be particularly valuable for studying the role(s) of DEAD-box proteins, as these tend to lead to larger conformational rearrangements, often resulting from the release of an RNA-binding protein from a bound RNA. Here we provide an adaptation of the DMS-MaPseq protocol to study RNA conformational transitions during ribosome assembly, which addresses the challenges arising from the presence of many assembly intermediates, all at concentrations far below that of mature ribosomes. While this protocol was developed for the yeast S. cerevisiae, we anticipate that it should be readily transferable to other model organisms for which affinity purification has been established.

A high-throughput assay for directly monitoring nucleolar rRNA biogenesis

Studies of the regulation of nucleolar function are critical for ascertaining clearer insights into the basic biological underpinnings of ribosome biogenesis (RB), and for future development of therapeutics to treat cancer and ribosomopathies. A number of high-throughput primary assays based on morphological alterations of the nucleolus can indirectly identify hits affecting RB. However, there is a need for a more direct high-throughput assay for a nucleolar function to further evaluate hits. Previous reports have monitored nucleolar rRNA biogenesis using 5-ethynyl uridine (5-EU) in low-throughput. We report a miniaturized, high-throughput 5-EU assay that enables specific calculation of nucleolar rRNA biogenesis inhibition, based on co-staining of the nucleolar protein fibrillarin (FBL). The assay uses two siRNA controls: a negative non-targeting siRNA control and a positive siRNA control targeting RNA Polymerase 1 (RNAP1; POLR1A), and specifically quantifies median 5-EU signal within nucleoli. Maximum nuclear 5-EU signal can also be used to monitor the effects of putative small-molecule inhibitors of RNAP1, like BMH-21, or other treatment conditions that cause FBL dispersion. We validate the 5-EU assay on 68 predominately nucleolar hits from a high-throughput primary screen, showing that 58/68 hits significantly inhibit nucleolar rRNA biogenesis. Our new method establishes direct quantification of nucleolar function in high-throughput, facilitating closer study of RB in health and disease.

Assembly of long L-RNA by native RNA ligation

Due to their intrinsic nuclease resistance, L-oligonucleotides are being increasingly utilized in the development of molecular tools and sensors. Yet, it remains challenging to synthesize long L-oligonucleotides, potential limiting future applications. Herein, we report straightforward and versitile approach to assemble long L-RNAs from two or more shorter fragments using T4 RNA ligase 1. We show that this approach is compatible with the assembly of several classes of functional L-RNA, which we highlight by generating a 124 nt L-RNA biosensor that functions in serum.

A New Specific and Sensitive RT-qPCR Method Based on Splinted 5' Ligation for the Quantitative Detection of RNA Species Shorter than microRNAs

Recently, we discovered a new family of unusually short RNAs mapping to 5.8S ribosomal RNA (rRNA) and which we named dodecaRNAs (doRNAs), according to the number of core nucleotides (12 nt) their members contain. To confirm these small RNA-sequencing (RNA-Seq) data, validate the existence of the two overly abundant doRNAs-the minimal core 12-nt doRNA sequence and its + 1-nt variant bearing a 5' Cytosine, C-doRNA-and streamline their analysis, we developed a new specific and sensitive splinted 5' ligation reverse transcription (RT)-quantitative polymerase chain reaction (qPCR) method. This method is based on a splint-assisted ligation of an adapter to the 5' end of doRNAs, followed by RT-qPCR amplification and quantitation. Our optimized protocol, which may discriminate between doRNA, C-doRNA, mutated and precursor sequences, can accurately detect as low as 240 copies and is quantitatively linear over a range of 7 logs. This method provides a unique tool to expand and facilitate studies exploring the molecular and cellular biology of RNA species shorter than microRNAs.

Regulation of reversible conformational change, size switching, and immunomodulation of RNA nanocubes

In biological systems, conformational changes and allosteric modulation play pivotal roles in regulating biological functions, such as the dynamic change of protein molecules, in response to binding or interacting with other factors such as pH, voltage, salt, light, or ligand. RNA can be manipulated and tuned with a level of simplicity that is characteristic of DNA or polymers, while displaying versatility in structure, diversity in function, and adaptability in a configuration similar to proteins. In the past, the work on the investigation of conformational change mainly focused on protein. The induced-fit and conformational capture in RNA have also been explored, such as in the study of riboswitches. Herein, we report the engineering of three-dimensional RNA nanocubes and demonstrated the operation and regulation for its configuration. We demonstrate the operation of reconfigurable RNA nanocubes whose shapes change precisely and reversibly in response to a specific trigger strand. The shape, size, and conformation can be regulated precisely and reversibly in response to the specific triggering signals. The shape and conformational conversion were observed by cryo-EM and gel electrophoresis, respectively. Harnessing the size, shape, conformation, and self-assembly capabilities of the RNA nanocube can provide a new potential use of this technology as nanocarriers for the treatment of various diseases.

Molecular diagnostic of toxigenic Corynebacterium diphtheriae strain by DNA sensor potentially suitable for electrochemical point-of-care diagnostic

The presented study is focused on the development of electrochemical genosensor for detection of tox gene fragment of toxigenic Corynebacterium diphtheriae strain. Together with our previous studies it fulfils the whole procedure for fast and accurate diagnostic of diphtheria at its early stage of infection with the use of electrochemical methods. The developed DNA sensor potentially can be used in more sophisticated portable device. After the electrochemical stem-loop probe structure optimization the conditions for real asymmetric PCR (aPCR) product detection were selected. As was shown it was crucial to optimize the magnesium and organic solvent concentrations in detection buffer. Under optimal conditions it was possible to selectively detect as low as 20.8 nM of complementary stand in 5 min or 0.5 nM in 30 min with sensitivity of 12.81 and 0.24 1⋅μM^-1 respectively. The unspecific biosensor response was elucidated with the use of new electrode blocking agent, diethyldithiocarbamate. Its application in electrochemical genosensors lead to significant higher current values and the biosensor response even in conditions with magnesium ion depletion. The developed biosensor selectivity was examined using samples containing genetic material originated from a number of non-target bacterial species which potentially can be present in the human upper respiratory tract.

Optimization and evaluation of fluorescence in situ hybridization chain reaction in cleared fresh-frozen brain tissues

Transcript labeling in intact tissues using in situ hybridization chain reaction has potential to provide vital spatiotemporal information for molecular characterization of heterogeneous neuronal populations. However, large tissue labeling in non-perfused or fresh-frozen rodent and postmortem human samples, which provide more flexible utilization than perfused tissues, is largely unexplored. In the present study, we optimized the combination of in situ hybridization chain reaction in fresh-frozen rodent brains and then evaluated the uniformity of neuronal labeling between two clearing methods, CLARITY and iDISCO⁺. We found that CLARITY yielded higher signal-to-noise ratios but more limited imaging depth and required longer clearing times, whereas, iDISCO⁺ resulted in better tissue clearing, greater imaging depth and a more uniform labeling of larger samples. Based on these results, we used iDISCO⁺-cleared fresh-frozen rodent brains to further validate this combination and map the expression of a few genes of interest pertaining to mood disorders. We then examined the potential of in situ hybridization chain reaction to label transcripts in cleared postmortem human brain tissues. The combination failed to produce adequate mRNA labeling in postmortem human cortical slices but produced visually adequate labeling in the cerebellum tissues. We next, investigated the multiplexing ability of in situ hybridization chain reaction in cleared tissues which revealed inconsistent fluorescence output depending upon the fluorophore conjugated to the hairpins. Finally, we applied our optimized protocol to assess the effect of glucocorticoid receptor overexpression on basal somatostatin expression in the mouse cortex. The constitutive glucocorticoid receptor overexpression resulted in lower number density of somatostatin-expressing neurons compared to wild type. Overall, the combination of in situ hybridization chain reaction with clearing methods, especially iDISCO⁺, may find broad application in the transcript analysis in rodent studies, but its limited use in postmortem human tissues can be improved by further optimizations.

Template-guided selection of RNA ligands using imine-based dynamic combinatorial chemistry

This study establishes the applicability of imine-based dynamic combinatorial chemistry to discover non-covalent ligands for RNA targets. We elucidate properties underlying the reactivity of arylamines and demonstrate target-guided amplification of tight binders in an amiloride-based dynamic library.

2-Amino-1,3-benzothiazole-6-carboxamide Preferentially Binds the Tandem Mismatch Motif r(UY:GA)

RNA helices are often punctuated with non-Watson-Crick features that may be targeted by chemical compounds, but progress toward identifying such compounds has been slow. We embedded a tandem UU:GA mismatch motif (5'-UG-3':5'-AU-3') within an RNA hairpin stem to identify compounds that bind the motif specifically. The three-dimensional structure of the RNA hairpin and its interaction with a small molecule identified through virtual screening are presented. The G-A mismatch forms a sheared pair upon which the U-U base pair stacks. The hydrogen bond configuration of the U-U pair involves O2 of the U adjacent to the G and O4 of the U adjacent to the A. The G-A and U-U pairs are flanked by A-U and G-C base pairs, respectively, and the stability of the mismatch is greater than when the motif is within the context of other flanking base pairs or when the 5'-3' orientation of the G-A and U-U pairs is swapped. Residual dipolar coupling constants were used to generate an ensemble of structures against which a virtual screen of 64480 small molecules was performed. The tandem mismatch was found to be specific for one compound, 2-amino-1,3-benzothiazole-6-carboxamide, which binds with moderate affinity but extends the motif to include the flanking A-U and G-C base pairs. The finding that the affinity for the UU:GA mismatch is dependent on flanking sequence emphasizes the importance of the motif context and potentially increases the number of small noncanonical features within RNA that can be specifically targeted by small molecules.

Challenges and perspectives for structural biology of lncRNAs-the example of the Xist lncRNA A-repeats

Following the discovery of numerous long non-coding RNA (lncRNA) transcripts in the human genome, their important roles in biology and human disease are emerging. Recent progress in experimental methods has enabled the identification of structural features of lncRNAs. However, determining high-resolution structures is challenging as lncRNAs are expected to be dynamic and adopt multiple conformations, which may be modulated by interaction with protein binding partners. The X-inactive specific transcript (Xist) is necessary for X inactivation during dosage compensation in female placental mammals and one of the best-studied lncRNAs. Recent progress has provided new insights into the domain organization, molecular features, and RNA binding proteins that interact with distinct regions of Xist. The A-repeats located at the 5′ end of the transcript are of particular interest as they are essential for mediating silencing of the inactive X chromosome. Here, we discuss recent progress with elucidating structural features of the Xist lncRNA, focusing on the A-repeats. We discuss the experimental and computational approaches employed that have led to distinct structural models, likely reflecting the intrinsic dynamics of this RNA. The presence of multiple dynamic conformations may also play an important role in the formation of the associated RNPs, thus influencing the molecular mechanism underlying the biological function of the Xist A-repeats. We propose that integrative approaches that combine biochemical experiments and high-resolution structural biology in vitro with chemical probing and functional studies in vivo are required to unravel the molecular mechanisms of lncRNAs.

Hidden Structural Modules in a Cooperative RNA Folding Transition

Large-scale, cooperative rearrangements underlie the functions of RNA in RNA-protein machines and gene regulation. To understand how such rearrangements are orchestrated, we used high-throughput chemical footprinting to dissect a seemingly concerted rearrangement in P5abc RNA, a paradigm of RNA folding studies. With mutations that systematically disrupt or restore putative structural elements, we found that this transition reflects local folding of structural modules, with modest and incremental cooperativity that results in concerted behavior. First, two distant secondary structure changes are coupled through a bridging three-way junction and Mg²⁺-dependent tertiary structure. Second, long-range contacts are formed between modules, resulting in additional cooperativity. Given the sparseness of RNA tertiary contacts after secondary structure formation, we expect that modular folding and incremental cooperativity are generally important for specifying functional structures while also providing productive kinetic paths to these structures. Additionally, we expect our approach to be useful for uncovering modularity in other complex RNAs.

Ligand-observed NMR techniques to probe RNA-small molecule interactions

A growing understanding of the structure and function of RNA has revealed it as a key regulator of gene expression and disease. A multitude of noncoding functions apart from the central roles of RNA in coding for and facilitating protein biogenesis has stimulated research into RNA as a pharmacological target. Despite many exciting advances, RNA remains an understudied target for small molecules, and techniques to investigate RNA-binding molecules are still emerging. A key stumbling block in this area has been validation of RNA-small molecule interactions. Our laboratory has recently used multiple ligand-observed NMR techniques in this regard, including CPMG and WaterLOGSY. This work describes methods to use these techniques in the context of studying RNA-ligand interactions.

Understanding the Contributions of Conformational Changes, Thermodynamics, and Kinetics of RNA-Small Molecule Interactions

The implication of RNA in multiple cellular processes beyond protein coding has revitalized interest in the development of small molecules for therapeutically targeting RNA and for further probing its cellular biology. However, the process of rationally designing such small molecule probes is hampered by the paucity of information about fundamental molecular recognition principles of RNA. In this Review, we summarize two important and often underappreciated aspects of RNA-small molecule recognition: RNA conformational dynamics and the biophysical properties of interactions of small molecules with RNA, specifically thermodynamics and kinetics. While conformational flexibility is often said to impede RNA ligand development, the ability of small molecules to influence the RNA conformational landscape can have a significant effect on the cellular functions of RNA. An analysis of the conformational landscape of RNA and the interactions of individual conformations with ligands can thus guide the development of new small molecule probes, which needs to be investigated further. Additionally, while it is common practice to quantify the binding affinities ( K_a or K_d) of small molecules for biomacromolecules as a measure of their activity, further biophysical characterization of their interaction can provide a deeper understanding. Studies that focus on the thermodynamic and kinetic parameters for interaction between RNA and ligands are next discussed. Finally, this Review provides the reader with a perspective on how such in-depth analysis of biophysical characteristics of the interaction of RNA and small molecules can impact our understanding of these interactions and how they will benefit the future design of small molecule probes.

Low dose dimethyl sulfoxide driven gross molecular changes have the potential to interfere with various cellular processes

Dimethyl sulfoxide (DMSO) is a small molecule with polar, aprotic and amphiphilic properties. It serves as a solvent for many polar and nonpolar molecules and continues to be one of the most used solvents (vehicle) in medical applications and scientific research. To better understand the cellular effects of DMSO within the concentration range commonly used as a vehicle (0.1-1.5%, v/v) for cellular treatments, we applied Attenuated Total Reflectance (ATR) Fourier Transform Infrared (FT-IR) spectroscopy to DMSO treated and untreated epithelial colon cancer cells. Both unsupervised (Principal Component Analysis-PCA) and supervised (Linear Discriminant Analysis-LDA) pattern recognition/modelling algorithms applied to the IR data revealed total segregation and prominent differences between DMSO treated and untreated cells at whole, lipid and nucleic acid regions. Several of these data were supported by other independent techniques. Further IR data analyses of macromolecular profile indicated comprehensive alterations especially in proteins and nucleic acids. Protein secondary structure analysis showed predominance of β-sheet over α-helix in DMSO treated cells. We also observed for the first time, a reduction in nucleic acid level upon DMSO treatment accompanied by the formation of Z-DNA. Molecular docking and binding free energy studies indicated a stabilization of Z-DNA in the presence of DMSO. This alternate DNA form may be related with the specific actions of DMSO on gene expression, differentiation, and epigenetic alterations. Using analytical tools combined with molecular and cellular biology techniques, our data indicate that even at very low concentrations, DMSO induces a number of changes in all macromolecules, which may affect experimental outcomes where DMSO is used as a solvent.

RNA as a Complex Polymer with Coupled Dynamics of Ions and Water in the Outer Solvation Sphere

We unravel the internal and collective modes of a widely studied 58-nucleotide rRNA fragment in solvent using atomically detailed molecular dynamics simulations. The variation of lifetimes for water hydrogen bonds with nucleotide groups indicates heterogeneity of water dynamics on the RNA surface. The time scales of interactions of the discrete water molecules with RNA nucleotides extend from several hundred picoseconds to a few nanoseconds. We determine all of the association sites and the spatial distribution of residence times for Mg²⁺, K⁺, and water molecules in those sites. We provide insights into the population of Mg²⁺ and K⁺ ions and water molecules in the outer sphere and how their fluctuations are intricately linked with the kinetics of the 58-mer. We find that many of the long-lived Mg²⁺ sites identified from the simulations agree with the locations of ions in the X-ray structure. We determine the excess ion atmosphere around the rRNA and compare it with experimental data. We investigate the collective behavior of RNA, ions, and water, by performing a joint principle component analysis for the Cartesian coordinates of the RNA phosphorus atoms and for the occupation counts of the association sites. Our results indicate that the 58-mer system is a complex polymer, composed of RNA that is encased by a fluctuating network of associated counterions, co-ions, and water.

Multi-environment fitness landscapes of a tRNA gene

A fitness landscape (FL) describes the genotype-fitness relationship in a given environment. To explain and predict evolution, it is imperative to measure the FL in multiple environments because the natural environment changes frequently. Using a high-throughput method that combines precise gene replacement with next-generation sequencing, we determine the in vivo FL of a yeast tRNA gene comprising over 23,000 genotypes in four environments. Although genotype-by-environment interaction (G×E) is abundantly detected, its pattern is so simple that we can transform an existing FL to that in a new environment with fitness measures of only a few genotypes in the new environment. Under each environment, we observe prevalent, negatively biased epistasis between mutations (G×G). Epistasis-by-environment interaction (G×G×E) is also prevalent, but trends in epistasis difference between environments are predictable. Our study thus reveals simple rules underlying seemingly complex FLs, opening the door to understanding and predicting FLs in general.

An anti-PDGFRβ aptamer for selective delivery of small therapeutic peptide to cardiac cells

Small therapeutic peptides represent a promising field for the treatment of pathologies such as cardiac diseases. However, the lack of proper target-selective carriers hampers their translation towards a potential clinical application. Aptamers are cell-specific carriers that bind with high affinity to their specific target. However, some limitations on their conjugation to small peptides and the functionality of the resulting aptamer-peptide chimera exist. Here, we generated a novel aptamer-peptide chimera through conjugation of the PDGFRβ-targeting Gint4.T aptamer to MP, a small mimetic peptide that via targeting of the Cavβ2 subunit of the L-type calcium channel (LTCC) can recover myocardial function in pathological heart conditions associated with defective LTCC function. The conjugation reaction was performed by click chemistry in the presence of N,N,N',N',N"-pentamethyldiethylenetriamine as a Cu (I) stabilizing agent in a DMSO-free aqueous buffer. When administered to cardiac cells, the Gint4.T-MP aptamer-peptide chimera was successfully internalized in cells, allowing the functional targeting of MP to LTCC. This approach represents the first example of the use of an internalizing aptamer for selective delivery of a small therapeutic peptide to cardiac cells.

Effect of Dimethyl Sulfoxide on the Binding of 1-Adamantane Carboxylic Acid to β- and γ-Cyclodextrins

Most therapeutic targets are proteins whose binding sites are hydrophobic cavities. For this reason, the majority of drugs under development are hydrophobic molecules exhibiting low solubility in water. To tackle this issue, a few percent of cosolvent, such as dimethyl sulfoxide (DMSO), is usually employed to increase drug solubility during the drug screening process. However, the few published studies dealing with the effect of adding DMSO showed that the affinity of hydrophobic ligands is systematically underestimated. To better understand the effect of DMSO, there is a need of studying its effect on a large range of systems. In this work, we used β- and γ-cyclodextrins (made of 6 and 7 α-d-glucopyranoside units, respectively) as models of hydrophobic cavities to investigate the effect of the addition 5% DMSO on the affinity of 1-adamantane carboxylic acid (ADA) to these cyclodextrins. The two systems differ by the size of the cyclodextrin cavity. The evaluation of binding constants was performed using ultrasound velocimetry, nuclear magnetic resonance spectroscopy, and molecular simulations. All techniques show that the presence of 5% DMSO does not significantly modify the affinity of ADA for γ-cyclodextrin, while the affinity is dramatically reduced for β-cyclodextrin. The bias induced by the presence of DMSO is thus more important when the ligand volume better fits the cyclodextrin cavity. Our work also suggests that free energy calculations provide a sound alternative to experimental techniques when dealing with poorly water-soluble drugs.

Evaluation of global HIV/SIV envelope gp120 RNA structure and evolution within and among infected hosts

Lentiviral RNA genomes contain structural elements that play critical roles in viral replication. Although structural features of 5'-untranslated regions have been well characterized, attempts to identify important structures in other genomic regions by Selective 2'-Hydroxyl Acylation analyzed by Primer Extension (SHAPE) have led to conflicting structural and mechanistic conclusions. Previous approaches accounted neither for sequence heterogeneity that is ubiquitous in viral populations, nor for selective constraints operating at the protein level. We developed an approach that augments SHAPE with phylogenetic analyses and applied it to investigate structure in coding regions (cRNA) within the HIV and SIV envelope genes. Analysis of single-genome SHAPE data with phylogenetic information from diverse lentiviral sequences argues against the conservation of a putative global gp120 RNA structure but points to the existence of core RNA sub-structures. Our findings establish a framework for considering sequence heterogeneity and protein function in de novo RNA structure inference approaches.

An in vivo high-throughput screening for riboswitch ligands using a reverse reporter gene system

Riboswitches are bacterial RNA elements that regulate gene expression in response to metabolite or ion abundance and are considered as potential drug targets. In recent years a number of methods to find non-natural riboswitch ligands have been described. Here we report a high-throughput in vivo screening system that allows identifying OFF-riboswitch modulators in a 384 well bioluminescence assay format. We use a reverse reporter gene setup in Bacillus subtilis, consisting of a primary screening assay, a secondary assay as well as counter assays to detect compounds in a library of 1,280 molecules that act on the guanine-responsive xpt riboswitch from B. anthracis. With this in vivo high-throughput approach we identified several hit compounds and could validate the impact of one of them on riboswitch-mediated gene regulation, albeit this might not be due to direct binding to the riboswitch. However, our data demonstrate the capability of our screening assay for bigger high-throughput screening campaigns. Furthermore, the screening system described here can not only be generally employed to detect non-natural ligands or compounds influencing riboswitches acting as genetic OFF switches, but it can also be used to investigate natural ligands of orphan OFF-riboswitches.

Protein Synthesis in Coupled and Uncoupled Cell-Free Prokaryotic Gene Expression Systems

Secondary structure formation of mRNA, caused by desynchronization of transcription and translation, is known to impact gene expression in vivo. Yet, inactivation of mRNA by secondary structures in cell-free protein expression is frequently overlooked. Transcription and translation rates are often not highly synchronized in cell-free expression systems, leading to a temporal mismatch between the processes and a drop in efficiency of protein production. By devising a cell-free gene expression platform in which transcriptional and translational elongation are successfully performed independently, we determine that sequence-dependent mRNA secondary structures are the main cause of mRNA inactivation in in vitro gene expression.

Shortening the HIV-1 TAR RNA Bulge by a Single Nucleotide Preserves Motional Modes over a Broad Range of Time Scales

Helix-junction-helix (HJH) motifs are flexible building blocks of RNA architecture that help define the orientation and dynamics of helical domains. They are also frequently involved in adaptive recognition of proteins and small molecules and in the formation of tertiary contacts. Here, we use a battery of nuclear magnetic resonance techniques to examine how deleting a single bulge residue (C24) from the human immunodeficiency virus type 1 (HIV-1) transactivation response element (TAR) trinucleotide bulge (U23-C24-U25) affects dynamics over a broad range of time scales. Shortening the bulge has an effect on picosecond-to-nanosecond interhelical and local bulge dynamics similar to that casued by increasing the Mg(2+) and Na(+) concentration, whereby a preexisting two-state equilibrium in TAR is shifted away from a bent flexible conformation toward a coaxial conformation, in which all three bulge residues are flipped out and flexible. Surprisingly, the point deletion minimally affects microsecond-to-millisecond conformational exchange directed toward two low-populated and short-lived excited conformational states that form through reshuffling of bases pairs throughout TAR. The mutant does, however, adopt a slightly different excited conformational state on the millisecond time scale, in which U23 is intrahelical, mimicking the expected conformation of residue C24 in the excited conformational state of wild-type TAR. Thus, minor changes in HJH topology preserve motional modes in RNA occurring over the picosecond-to-millisecond time scales but alter the relative populations of the sampled states or cause subtle changes in their conformational features.

The structural stability and catalytic activity of DNA and RNA oligonucleotides in the presence of organic solvents

Organic solvents and apolar media are used in the studies of nucleic acids to modify the conformation and function of nucleic acids, to improve solubility of hydrophobic ligands, to construct molecular scaffolds for organic synthesis, and to study molecular crowding effects. Understanding how organic solvents affect nucleic acid interactions and identifying the factors that dominate solvent effects are important for the creation of oligonucleotide-based technologies. This review describes the structural and catalytic properties of DNA and RNA oligonucleotides in organic solutions and in aqueous solutions with organic cosolvents. There are several possible mechanisms underlying the effects of organic solvents on nucleic acid interactions. The reported results emphasize the significance of the osmotic pressure effect and the dielectric constant effect in addition to specific interactions with nucleic acid strands. This review will serve as a guide for the selection of solvent systems based on the purpose of the nucleic acid-based experiments.

Rapid Exchange Between Free and Bound States in RNA-Dendrimer Polyplexes: Implications on the Mechanism of Delivery and Release

A combination of solution NMR, dynamic light scattering (DLS), and fluorescence quenching assays were employed to obtain insights into the dynamics and structural features of a polyplex system consisting of HIV-1 transactivation response element (TAR) and PEGylated generation 5 poly(amidoamine) dendrimer (G5-PEG). NMR chemical shift mapping and (13)C spin relaxation based dynamics measurements depict the polyplex system as a highly dynamic assembly where the RNA, with its local structure and dynamics preserved, rapidly exchanges ( Collapse

Key Words

• therapeutic rna delivery
• pamam
• pegylation
• tar
• dls
• nmr spin relaxation

Collapse

MESH Headings

• Dendrimers/chemistry
• HIV-1/chemistry
• Magnetic Resonance Spectroscopy/methods
• Polyamines/chemistry
• Polyethylene Glycols/chemistry
• RNA/chemistry
• Response Elements/genetics
• Transfection/methods

Collapse

Grants

• R01 AI066975 NIAID NIH HHS
• P50 GM103297 NIGMS NIH HHS
• R01 AI087463-01 NIAID NIH HHS
• R01 EB005028 NIBIB NIH HHS
• R01EB005028 NIBIB NIH HHS
• R01 GM089846 NIGMS NIH HHS

Collapse

Affiliation(s)

Anisha Shakya Department of Chemistry, University of Michigan, Ann Arbor, MI 48109-1055, USA
Casey A. Dougherty Department of Chemistry, University of Michigan, Ann Arbor, MI 48109-1055, USA
Yi Xue Department of Biochemistry and Chemistry, Duke University Medical Center, Durham, NC 27710, USA
Hashim M. Al-Hashimi Department of Biochemistry and Chemistry, Duke University Medical Center, Durham, NC 27710, USA
Mark M. Banaszak Holl Department of Chemistry, University of Michigan, Ann Arbor, MI 48109-1055, USA

Collapse

Fluorescent DNA-Protected Silver Nanoclusters for Ligand-HIV RNA Interaction Assay

Studying ligand-biomacromolecule interactions provides opportunities for creating new compounds that can efficiently regulate specific biological processes. Ribonucleic acid (RNA) molecules have become attractive drug targets since the discovery of their roles in modulating gene expression, while only a limited number of studies have investigated interactions between ligands and functional RNA molecules, especially those based on nanotechnology. DNA-protected silver nanoclusters (AgNCs) were used to investigate ligand-RNA interactions for the first time in this study. The anthracycline anticancer drug mitoxantrone (MTX) was found to quench the fluorescence of AgNCs. After adding human immunodeficiency virus trans-activation responsive region (TAR) RNA or Rev-response element (RRE) RNA to AgNCs-MTX mixtures, the fluorescence of the AgNCs recovered due to interactions between MTX with RNAs. The binding constants and number of binding sites of MTX to TAR and RRE RNA were determined through theoretical calculations. MTX-RNA interactions were further confirmed in fluorescence polarization and mass spectrometry experiments. The mechanism of MTX-based fluorescence quenching of the AgNCs was also explored. This study provides a new strategy for ligand-RNA binding interaction assay.

Factors that influence T box riboswitch efficacy and tRNA affinity

The T box riboswitch is an intriguing potential target for antibacterial drug discovery. Found primarily in Gram-positive bacteria, the riboswitch regulates gene expression by selectively responding to uncharged tRNA to control transcription readthrough. Polyamines and molecular crowding are known to specifically affect RNA function, but their effect on T box riboswitch efficacy and tRNA affinity have not been fully characterized. A fluorescence-monitored in vitro transcription assay was developed to readily quantify these molecular interactions and to provide a moderate-throughput functional assay for a comprehensive drug discovery screening cascade. The polyamine spermidine specifically enhanced T box riboswitch readthrough efficacy with an EC50 = 0.58 mM independent of tRNA binding. Molecular crowding, simulated by the addition of polyethylene glycol, had no effect on tRNA affinity for the riboswitch, but did reduce the efficacy of tRNA-induced readthrough. These results indicate that the T box riboswitch tRNA affinity and readthrough efficacy are intricately modulated by environmental factors.

Fluorescence assays for monitoring RNA-ligand interactions and riboswitch-targeted drug discovery screening

Riboswitches and other noncoding regulatory RNA are intriguing targets for the development of therapeutic agents. A significant challenge in the drug discovery process, however, is the identification of potent compounds that bind the target RNA specifically and disrupt its function. Essential to this process is an effectively designed cascade of screening assays. A screening cascade for identifying compounds that target the T box riboswitch antiterminator element is described. In the primary assays, moderate to higher throughput screening of compound libraries is achieved by combining the sensitivity of fluorescence techniques with functionally relevant assays. Active compounds are then validated and the binding to target RNA further characterized in secondary assays. The cascade of assays monitor ligand-induced changes in the steady-state fluorescence of an attached dye or internally incorporated 2-aminopurine; the fluorescence anisotropy of an RNA complex; and, the thermal denaturation fluorescence profile of a fluorophore-quencher labeled RNA. While the assays described have been developed for T box riboswitch-targeted drug discovery, the fluorescence methods and screening cascade design principles can be applied to drug discovery efforts targeted toward other medicinally relevant noncoding RNA.