Anion binding to nucleic acids

Principles of ion binding to RNA inferred from the analysis of a 1.55 Å resolution bacterial ribosome structure - Part I: Mg2

The importance of Mg2+ ions for RNA structure and function cannot be overstated. Several attempts were made to establish a comprehensive Mg2+ binding site classification. However, such descriptions were hampered by poorly modelled ion binding sites as observed in a recent cryo-EM 1.55 Å Escherichia coli ribosome structure where incomplete ion assignments blurred our understanding of their binding patterns. We revisited this model to establish general binding principles applicable to any RNA of sufficient resolution. These principles rely on the 2.9 Å distance separating two water molecules bound in cis to Mg2+. By applying these rules, we could assign all Mg2+ ions bound with 2-4 non-water oxygens. We also uncovered unanticipated motifs where up to five adjacent nucleotides wrap around a single ion. The formation of such motifs involves a hierarchical Mg2+ ion dehydration process that plays a significant role in ribosome biogenesis and in the folding of large RNAs. Besides, we established a classification of the Mg2+…Mg2+ and Mg2+…K+ ion pairs observed in this ribosome. Overall, the uncovered binding principles enhance our understanding of the roles of ions in RNA structure and will help refining the solvation shell of other RNA systems.

Strategies for targeted gene delivery using lipid nanoparticles and cell-derived nanovesicles

Gene therapy is a promising approach for the treatment of many diseases. However, the effective delivery of the cargo without degradation in vivo is one of the major hurdles. With the advent of lipid nanoparticles (LNPs) and cell-derived nanovesicles (CDNs), gene delivery holds a very promising future. The targeting of these nanosystems is a prerequisite for effective transfection with minimal side-effects. In this review, we highlight the emerging strategies utilized for the effective targeting of LNPs and CDNs, and we summarize the preparation methodologies for LNPs and CDNs. We have also highlighted the non-ligand targeting of LNPs toward certain organs based on their composition. It is highly expected that continuing the developments in the targeting approaches of LNPs and CDNs for the delivery system will further promote them in clinical translation.

Anionic G•U pairs in bacterial ribosomal rRNAs

Wobble GU pairs (or G•U) occur frequently within double-stranded RNA helices interspersed between standard G=C and A-U Watson-Crick pairs. Another type of G•U pair interacting via their Watson-Crick edges has been observed in the A site of ribosome structures between a modified U34 in the tRNA anticodon triplet and G + 3 in the mRNA. In such pairs, the electronic structure of the U is changed with a negative charge on N3(U), resulting in two H-bonds between N1(G)…O4(U) and N2(G)…N3(U). Here, we report that such pairs occur in other highly conserved positions in ribosomal RNAs of bacteria in the absence of U modification. An anionic cis Watson-Crick G•G pair is also observed and well conserved in the small subunit. These pairs are observed in tightly folded regions.

Oligonucleotide nanoassemblies with allyl bromide scaffold-based small molecules

The development of oligonucleotide nanoassemblies with small molecules has shown great potential in bio-medical applications. However, the interaction of negatively charged oligonucleotides with halogenated small molecules represents a scientific challenge. Here, we introduced a distinct allyl bromide halogenated scaffold, which exhibits specific interaction with adenine nucleic bases of the oligonucleotides, thus leading to the formation of self-assembled nanostructures.

Structural, Dynamical, and Entropic Differences between SARS-CoV and SARS-CoV-2 s2m Elements Using Molecular Dynamics Simulations

The functional role of the highly conserved stem-loop II motif (s2m) in SARS-CoV and SARS-CoV-2 in the viral lifecycle remains enigmatic and an intense area of research. Structure and dynamics of the s2m are key to establishing a structure-function connection, yet a full set of atomistic resolution coordinates is not available for SARS-CoV-2 s2m. Our work constructs three-dimensional coordinates consistent with NMR solution phase data for SARS-CoV-2 s2m and provides a comparative analysis with its counterpart SARS-CoV s2m. We employed initial coordinates based on PDB ID 1XJR for SARS-CoV s2m and two models for SARS-CoV-2 s2m: one based on 1XJR in which we introduced the mutations present in SARS-CoV-2 s2m and the second based on the available SARS-CoV-2 NMR NOE data supplemented with knowledge-based methods. For each of the three systems, 3.5 μs molecular dynamics simulations were used to sample the structure and dynamics, and principal component analysis (PCA) reduced the ensembles to hierarchal conformational substates for detailed analysis. Dilute solution simulations of SARS-CoV s2m demonstrate that the GNRA-like terminal pentaloop is rigidly defined by base stacking uniquely positioned for possible kissing dimer formation. However, the SARS-CoV-2 s2m simulation did not retain the reported crystallographic SARS-CoV motifs and the terminal loop expands to a highly dynamic "nonaloop." Increased flexibility and structural disorganization are observed for the larger terminal loop, where an entropic penalty is computed to explain the experimentally observed reduction in kissing complex formation. Overall, both SARS-CoV and SARS-CoV-2 s2m elements have a similarly pronounced L-shape due to different motif interactions. Our study establishes the atomistic three-dimensional structure and uncovers dynamic differences that arise from s2m sequence changes, which sets the stage for the interrogation of different mechanistic pathways of suspected biological function.

Insights into DNA solvation found in protein-DNA structures

The proteins that bind double-helical DNA present various microenvironments that sense and/or induce signals in the genetic material. The high-resolution structures of protein-DNA complexes reveal the nature of both the microenvironments and the conformational responses in DNA and protein. Complex networks of interactions within the structures somehow tie the protein and DNA together and induce the observed spatial forms. Here we show how the cumulative buildup of amino acid atoms around the sugars, phosphates, and bases in different protein-DNA complexes produces a binding cloud around the double helix and how different types of atoms fill that cloud. Rather than focusing on the principles of molecular binding and recognition suggested by the arrangements of amino acids and nucleotides in the macromolecular complexes, we consider the proteins in contact with DNA as organized solvents. We describe differences in the mix of atoms that come in closest contact with DNA, subtle sequence-dependent features in the microenvironment of the sugar-phosphate backbone, a direct link between the localized buildup of ionic species and the electrostatic potential surfaces of the DNA bases, and sites of atomic buildup above and below the basepair planes that transmit the unique features of the base environments along the chain backbone. The inferences about solvation that can be drawn from the survey provide new stimuli for improvement of nucleic acid force fields and fresh ideas for exploration of the properties of DNA in solution.

Structure and thermodynamics of a UGG motif interacting with Ba2+ and other metal ions: accommodating changes in the RNA structure and the presence of a G(syn)-G(syn) pair

The self-complementary triplet 5'UGG3'/5'UGG3' is a particular structural motif containing noncanonical G-G pair and two U·G wobble pairs. It constitutes a specific structural and electrostatic environment attracting metal ions, particularly Ba2+ ions. Crystallographic research has shown that two Ba2+ cations are located in the major groove of the helix and interact directly with the UGG triplet. A comparison with the unliganded structure has revealed global changes in the RNA structure in the presence of metal ions, whereas thermodynamic measurements have shown increased stability. Moreover, in the structure with Ba2+, an unusual noncanonical G(syn)-G(syn) pair is observed instead of the common G(anti)-G(syn). We further elucidate the metal binding properties of the UGG/UGG triplet by performing crystallographic and thermodynamic studies using DSC and UV melting with other metal ions. The results explain the preferences of the UGG sequence for Ba2+ cations and point to possible applications of this metal-binding propensity.

The influence of Holliday junction sequence and dynamics on DNA crystal self-assembly

The programmable synthesis of rationally engineered crystal architectures for the precise arrangement of molecular species is a foundational goal in nanotechnology, and DNA has become one of the most prominent molecules for the construction of these materials. In particular, branched DNA junctions have been used as the central building block for the assembly of 3D lattices. Here, crystallography is used to probe the effect of all 36 immobile Holliday junction sequences on self-assembling DNA crystals. Contrary to the established paradigm in the field, most junctions yield crystals, with some enhancing the resolution or resulting in unique crystal symmetries. Unexpectedly, even the sequence adjacent to the junction has a significant effect on the crystal assemblies. Six of the immobile junction sequences are completely resistant to crystallization and thus deemed “fatal,” and molecular dynamics simulations reveal that these junctions invariably lack two discrete ion binding sites that are pivotal for crystal formation. The structures and dynamics detailed here could be used to inform future designs of both crystals and DNA nanostructures more broadly, and have potential implications for the molecular engineering of applied nanoelectronics, nanophotonics, and catalysis within the crystalline context.

Engineered crystal architectures from DNA have become a foundational goal for nanotechnological precise arrangement. Here, the authors systematically investigate the structures of 36 immobile Holliday junction sequences and identify the features allowing the crystallisation of most of them, while 6 are considered fatal.

fingeRNAt—A novel tool for high-throughput analysis of nucleic acid-ligand interactions

Computational methods play a pivotal role in drug discovery and are widely applied in virtual screening, structure optimization, and compound activity profiling. Over the last decades, almost all the attention in medicinal chemistry has been directed to protein-ligand binding, and computational tools have been created with this target in mind. With novel discoveries of functional RNAs and their possible applications, RNAs have gained considerable attention as potential drug targets. However, the availability of bioinformatics tools for nucleic acids is limited. Here, we introduce fingeRNAt—a software tool for detecting non-covalent interactions formed in complexes of nucleic acids with ligands. The program detects nine types of interactions: (i) hydrogen and (ii) halogen bonds, (iii) cation-anion, (iv) pi-cation, (v) pi-anion, (vi) pi-stacking, (vii) inorganic ion-mediated, (viii) water-mediated, and (ix) lipophilic interactions. However, the scope of detected interactions can be easily expanded using a simple plugin system. In addition, detected interactions can be visualized using the associated PyMOL plugin, which facilitates the analysis of medium-throughput molecular complexes. Interactions are also encoded and stored as a bioinformatics-friendly Structural Interaction Fingerprint (SIFt)—a binary string where the respective bit in the fingerprint is set to 1 if a particular interaction is present and to 0 otherwise. This output format, in turn, enables high-throughput analysis of interaction data using data analysis techniques. We present applications of fingeRNAt-generated interaction fingerprints for visual and computational analysis of RNA-ligand complexes, including analysis of interactions formed in experimentally determined RNA-small molecule ligand complexes deposited in the Protein Data Bank. We propose interaction fingerprint-based similarity as an alternative measure to RMSD to recapitulate complexes with similar interactions but different folding. We present an application of interaction fingerprints for the clustering of molecular complexes. This approach can be used to group ligands that form similar binding networks and thus have similar biological properties. The fingeRNAt software is freely available at https://github.com/n-szulc/fingeRNAt.

Tris(2-Pyridylmethylamine)V(O)2 Complexes as Counter Ions of Diprotonated Decavanadate Anion: Potential Antineoplastic Activity

The synthesis and theoretical-experimental characterization of a novel diprotanated decavanadate is presented here due to our search for novel anticancer metallodrugs. Tris(2-pyridylmethyl)amine (TPMA), which is also known to have anticancer activity in osteosarcoma cell lines, was introduced as a possible cationic species that could act as a counterpart for the decavanadate anion. However, the isolated compound contains the previously reported vanadium (V) dioxido-tpma moieties, and the decavanadate anion appears to be diprotonated. The structural characterization of the compound was performed by infrared spectroscopy and single-crystal X-ray diffraction. In addition, DFT calculations were used to analyze the reactive sites involved in the donor-acceptor interactions from the molecular electrostatic potential maps. The level of theory mPW1PW91/6–31G(d)-LANL2DZ and ECP = LANL2DZ for the V atom was used. These insights about the compounds’ main interactions were supported by analyzing the noncovalent interactions utilizing the AIM and Hirshfeld surfaces approach. Molecular docking studies with small RNA fragments were used to assess the hypothesis that decavanadate’s anticancer activity could be attributed to its interaction with lncRNA molecules. Thus, a combination of three potentially beneficial components could be evaluated in various cancer cell lines.

The crystal structure of bromide-bound GtACR1 reveals a pre-activated state in the transmembrane anion tunnel

The crystal structure of the light-gated anion channel GtACR1 reported in our previous Research Article (Li et al., 2019) revealed a continuous tunnel traversing the protein from extracellular to intracellular pores. We proposed the tunnel as the conductance channel closed by three constrictions: C1 in the extracellular half, mid-membrane C2 containing the photoactive site, and C3 on the cytoplasmic side. Reported here, the crystal structure of bromide-bound GtACR1 reveals structural changes that relax the C1 and C3 constrictions, including a novel salt-bridge switch mechanism involving C1 and the photoactive site. These findings indicate that substrate binding induces a transition from an inactivated state to a pre-activated state in the dark that facilitates channel opening by reducing free energy in the tunnel constrictions. The results provide direct evidence that the tunnel is the closed form of the channel of GtACR1 and shed light on the light-gated channel activation mechanism.

Berberrubine Phosphate: A Selective Fluorescent Probe for Quadruplex DNA

A phosphate-substituted, zwitterionic berberine derivative was synthesized and its binding properties with duplex DNA and G4-DNA were studied using photometric, fluorimetric and polarimetric titrations and thermal DNA denaturation experiments. The ligand binds with high affinity toward both DNA forms (Kb = 2-7 × 105 M-1) and induces a slight stabilization of G4-DNA toward thermally induced unfolding, mostly pronounced for the telomeric quadruplex 22AG. The ligand likely binds by aggregation and intercalation with ct DNA and by terminal stacking with G4-DNA. Thus, this compound represents one of the rare examples of phosphate-substituted DNA binders. In an aqueous solution, the title compound has a very weak fluorescence intensity (Φfl < 0.01) that increases significantly upon binding to G4-DNA (Φfl = 0.01). In contrast, the association with duplex DNA was not accompanied by such a strong fluorescence light-up effect (Φfl < 0.01). These different fluorimetric responses upon binding to particular DNA forms are proposed to be caused by the different binding modes and may be used for the selective fluorimetric detection of G4-DNA.

Infrared Nanospectroscopy Reveals DNA Structural Modifications upon Immobilization onto Clay Nanotubes

The growing demand for innovative means in biomedical, therapeutic and diagnostic sciences has led to the development of nanomedicine. In this context, naturally occurring tubular nanostructures composed of rolled sheets of alumino-silicates, known as halloysite nanotubes, have found wide application. Halloysite nanotubes indeed have surface properties that favor the selective loading of biomolecules. Here, we present the first, to our knowledge, structural study of DNA-decorated halloysite nanotubes, carried out with nanometric spatially-resolved infrared spectroscopy. Single nanotube absorption measurements indicate a partial covering of halloysite by DNA molecules, which show significant structural modifications taking place upon loading. The present study highlights the constraints for the use of nanostructured clays as DNA carriers and demonstrates the power of super-resolved infrared spectroscopy as an effective and versatile tool for the evaluation of immobilization processes in the context of drug delivery and gene transfer.

Effects of Contamination of Freshwater Habitat With Common Heavy Metals and Anions on the Prevalence of Human Adenoviruses and Enteroviruses

The occurrence and survival of enteric viruses in open surface waters can be impacted by a host of factors including fecal emission levels, seasonal variations, virus stability and the physicochemical parameters. In this research, we aimed to document the association between contaminations of water samples with human enteric viruses (adenoviruses and enteroviruses) from a freshwater lake with variations in chemical contaminants. We collected 216 water samples from October 2010 to April 2012, from a 4 km stretch along Lake Victoria (LV) basin in Homa Bay town located in the western region of Kenya. The samples were analyzed for the existence of human adenoviruses (HAdV) and human enteroviruses (HEV), using the nested PCR (nPCR). We also assessed in the water samples the levels of twelve chemical contaminants consisting of six heavy metal elements and six anions. About 8.3 % of the samples were found to be contaminated with the enteric viruses. The concentrations of the 12 chemical contaminants were found to be largely within the WHO suggested limits. Most of the chemical contaminants were not related to the detection rates of the viruses from the statistical analysis. However, some positive and negative associations between the viral genome's detection and the chemical concentrations were established for only three metals (Fe, Pb, Cd) and the PO4³⁻ Radical. Cd had a weak positive significant relationship with HAdV (rho = 0.146, p = 0.032) while Pb and Fe had a weak positive significant relationship with HEV genome detection (rho = 0.156, p = 0.022) and (rho = 0.148 and p = 0.029) respectively. There was a modest negative relationship between phosphate ions and HEV (rho = −0.174, p = 0.010). The results of our study do not provide support for the hypothesis of an association between the presence of human enteric viruses and the levels of twelve chemical contaminants.

Deflating the RNA Mg2+ bubble. Stereochemistry to the rescue!

Proper evaluation of the ionic structure of biomolecular systems through X ray and cryo-EM techniques remains challenging but is essential for advancing our understanding of the underlying structure/activity/solvent relationships. However, numerous studies overestimate the number of Mg2+ in deposited structures due to assignment errors finding their origin in improper consideration of stereochemical rules. Herein, to tackle such issues, we re-evaluate the PDBid 6QNR and 6SJ6 models of the ribosome ionic structure. We establish that stereochemical principles need to be carefully pondered when evaluating ion binding features, even when K+ anomalous signals are available as it is the case for the 6QNR PDB entry. For ribosomes, assignment errors can result in misleading conceptions of their solvent structure. For instance, present stereochemical analysis result in a significant decrease of the number of assigned Mg2+ in 6QNR, suggesting that K+ and not Mg2+ is the prevalent ion in the ribosome 1st solvation shell. We stress that the use of proper stereochemical guidelines in combination or not with other identification techniques, such as those pertaining to the detection of transition metals, of some anions and of K+ anomalous signals, is critical for deflating the current Mg2+ bubble witnessed in many ribosome and other RNA structures. We also stress that for the identification of lighter ions such as Mg2+, Na+, …, for which no anomalous signals can be detected, stereochemistry coupled with high resolution structures (<2.4 Å) remain the best currently available option.

Extraction of Small RNAs by Titanium Dioxide Nanofibers

MicroRNAs (miRNAs) are small RNAs, that bind to mRNA targets and regulate their translation. Functional study of miRNAs and exploration of their utility as disease markers require miRNA extraction from biological samples, which contain large amounts of interfering compounds for downstream RNA identification and quantification. The most common extraction methods employ either silica columns or TRIzol reagent, but these approaches afford low recovery for small RNAs, possibly due to their short strand lengths. Here, we describe the fabrication of titanium dioxide nanofibers and the optimal extraction conditions to improve miRNA recovery from biological buffers, cell lysate, and serum.

Duplex Electrochemical DNA Sensor to Detect Bacillus anthracis CAP and PAG DNA Targets Based on the Incorporation of Tailed Primers and Ferrocene-Labeled dATP

We report the duplex amplification of two plasmid DNA markers involved in the virulence of Bacillus anthracis, CAP and PAG, and the direct electrochemical detection of these amplicons. The method consists of the simultaneous amplification of the two targets in a single-pot reaction via polymerase chain reaction (PCR) using tailed primers and ferrocene-labeled dATP. Following amplification, the PCR products hybridize to probes immobilized on electrodes in a microfabricated electrode array chip. The incorporated ferrocene labeled dATP is then detected using square wave voltammetry. We evaluated the effect of electrolyte cations, anions, and concentration to condense, bend, and shrink double-stranded DNA and their effect on the intensity of the ferrocene signal. We obtained detection limits of 0.8 and 3.4 fM for CAP and PAG targets, respectively. We successfully developed a method to detect the presence of both targets in genomic DNA extracted from real samples.

Nucleobase carbonyl groups are poor Mg2+ inner-sphere binders but excellent monovalent ion binders-a critical PDB survey

Precise knowledge of Mg²⁺ inner-sphere binding site properties is vital for understanding the structure and function of nucleic acid systems. Unfortunately, the PDB, which represents the main source of Mg²⁺ binding sites, contains a substantial number of assignment issues that blur our understanding of the functions of these ions. Here, following a previous study devoted to Mg²⁺ binding to nucleobase nitrogens, we surveyed nucleic acid X-ray structures from the PDB with resolutions ≤2.9 Å to classify the Mg²⁺ inner-sphere binding patterns to nucleotide carbonyl, ribose hydroxyl, cyclic ether, and phosphodiester oxygen atoms. From this classification, we derived a set of "prior-knowledge" nucleobase Mg²⁺ binding sites. We report that crystallographic examples of trustworthy nucleobase Mg²⁺ binding sites are fewer than expected since many of those are associated with misidentified Na⁺ or K⁺ We also emphasize that binding of Na⁺ and K⁺ to nucleic acids is much more frequent than anticipated. Overall, we provide evidence derived from X-ray structures that nucleobases are poor inner-sphere binders for Mg²⁺ but good binders for monovalent ions. Based on strict stereochemical criteria, we propose an extended set of guidelines designed to help in the assignment and validation of ions directly contacting nucleobase and ribose atoms. These guidelines should help in the interpretation of X-ray and cryo-EM solvent density maps. When borderline Mg²⁺ stereochemistry is observed, alternative placement of Na⁺, K⁺, or Ca²⁺ must be considered. We also critically examine the use of lanthanides (Yb³⁺, Tb³⁺) as Mg²⁺ substitutes in crystallography experiments.

Structure-Activity Relationship of Flavin Analogues That Target the Flavin Mononucleotide Riboswitch

The flavin mononucleotide (FMN) riboswitch is an emerging target for the development of novel RNA-targeting antibiotics. We previously discovered an FMN derivative, 5FDQD, that protects mice against diarrhea-causing Clostridium difficile bacteria. Here, we present the structure-based drug design strategy that led to the discovery of this fluoro-phenyl derivative with antibacterial properties. This approach involved the following stages: (1) structural analysis of all available free and bound FMN riboswitch structures; (2) design, synthesis, and purification of derivatives; (3) in vitro testing for productive binding using two chemical probing methods; (4) in vitro transcription termination assays; and (5) resolution of the crystal structures of the FMN riboswitch in complex with the most mature candidates. In the process, we delineated principles for productive binding to this riboswitch, thereby demonstrating the effectiveness of a coordinated structure-guided approach to designing drugs against RNA.

RNA Structural Dynamics As Captured by Molecular Simulations: A Comprehensive Overview

With both catalytic and genetic functions, ribonucleic acid (RNA) is perhaps the most pluripotent chemical species in molecular biology, and its functions are intimately linked to its structure and dynamics. Computer simulations, and in particular atomistic molecular dynamics (MD), allow structural dynamics of biomolecular systems to be investigated with unprecedented temporal and spatial resolution. We here provide a comprehensive overview of the fast-developing field of MD simulations of RNA molecules. We begin with an in-depth, evaluatory coverage of the most fundamental methodological challenges that set the basis for the future development of the field, in particular, the current developments and inherent physical limitations of the atomistic force fields and the recent advances in a broad spectrum of enhanced sampling methods. We also survey the closely related field of coarse-grained modeling of RNA systems. After dealing with the methodological aspects, we provide an exhaustive overview of the available RNA simulation literature, ranging from studies of the smallest RNA oligonucleotides to investigations of the entire ribosome. Our review encompasses tetranucleotides, tetraloops, a number of small RNA motifs, A-helix RNA, kissing-loop complexes, the TAR RNA element, the decoding center and other important regions of the ribosome, as well as assorted others systems. Extended sections are devoted to RNA-ion interactions, ribozymes, riboswitches, and protein/RNA complexes. Our overview is written for as broad of an audience as possible, aiming to provide a much-needed interdisciplinary bridge between computation and experiment, together with a perspective on the future of the field.

Theoretical Prediction of Robust Second-Row Oxyanion Clusters in the Metastable Domain of Antielectrostatic Hydrogen Bonding

We provide ab initio and density functional theory evidence for a family of surprisingly robust like-charged clusters of common HSO₄^- and H₂PO₄^- oxyanions, ranging up to tetramers of net charge 4-. Our results support other recent theoretical and experimental evidence for "antielectrostatic" hydrogen-bonded (AEHB) species that challenge conventional electrostatic conceptions and force-field modeling of closed-shell ion interactions. We provide structural and energetic descriptors of the predicted kinetic well-depths (in the range 3-10 kcal/mol) and barrier widths (in the range 2-4 Å) for simple AEHB dimers, including evidence of extremely strong hydrogen bonding in the fluoride-bisulfate dianion. For more complex polyanionic species, we employ natural-bond-orbital-based descriptors to characterize the electronic features of the cooperative hydrogen-bonding network that are able to successfully defy Coulomb explosion. The computational results suggest a variety of kinetically stable AEHB species that may be suitable for experimental detection as long-lived gas-phase species or structural units of condensed phases, despite the imposing electrostatic barriers that oppose their formation under ambient conditions.

Extraction of microRNAs from biological matrices with titanium dioxide nanofibers

MicroRNAs (miRNAs) are small RNAs that bind to mRNA targets and regulate their translation. A functional study of miRNAs and exploration of their utility as disease markers require miRNA extraction from biological samples, which contain large amounts of interfering compounds for downstream RNA identification and quantification. The most common extraction methods employ silica columns or the TRIzol reagent but give out low recovery for small RNAs probably due to their short strand lengths. Herein, we fabricated the titanium dioxide nanofibers using electrospinning to facilitate miRNA extraction and developed the optimal buffer conditions to improve miRNA recovery from biological matrices of cell lysate and serum. We found that our TiO₂ fibers could obtain a recovery of 18.0 ± 3.6% for miRNA fibers while carrying out the extraction in the more complex medium of cell lysate, much higher than the 0.02 ± 0.0001% recovery from the commercial kit. The much improved extraction of miRNAs from our fibers could be originated from the strong coordination between TiO₂ and RNA's phosphate backbone. In addition, the binding, washing, and elution buffers judiciously developed in the present study can achieve selective extraction of small RNA shorter than 500 nucleotides in length. Our results demonstrate that TiO₂ nanofibers can work as a valuable tool for extraction of miRNAs from biological samples with high recovery. Graphical abstract Schematic for extraction of small RNAs using TiO₂ nanofibers.

Mg2+ ions: do they bind to nucleobase nitrogens?

Given the many roles proposed for Mg2+ in nucleic acids, it is essential to accurately determine their binding modes. Here, we surveyed the PDB to classify Mg2+ inner-sphere binding patterns to nucleobase imine N1/N3/N7 atoms. Among those, purine N7 atoms are considered to be the best nucleobase binding sites for divalent metals. Further, Mg2+ coordination to N7 has been implied in several ribozyme catalytic mechanisms. We report that Mg2+ assigned near imine nitrogens derive mostly from poor interpretations of electron density patterns and are most often misidentified Na+, K+, NH4+ ions, water molecules or spurious density peaks. Consequently, apart from few documented exceptions, Mg2+ ions do not bind to N7 atoms. Without much of a surprise, Mn2+, Zn2+ and Cd2+, which have a higher affinity for nitrogens, may contact N7 atoms when present in crystallization buffers. In this respect, we describe for the first time a potential Zn2+ ribosomal binding site involving two purine N7 atoms. Further, we provide a set of guidelines to help in the assignment of Mg2+ in crystallographic, cryo-EM, NMR and model building practices and discuss implications of our findings related to ion substitution experiments.

Structure of a yeast activated spliceosome at 3.5 Å resolution

Pre-messenger RNA (pre-mRNA) splicing is carried out by the spliceosome, which undergoes an intricate assembly and activation process. Here, we report an atomic structure of an activated spliceosome (known as the B(act) complex) from Saccharomyces cerevisiae, determined by cryo-electron microscopy at an average resolution of 3.52 angstroms. The final refined model contains U2 and U5 small nuclear ribonucleoprotein particles (snRNPs), U6 small nuclear RNA (snRNA), nineteen complex (NTC), NTC-related (NTR) protein, and a 71-nucleotide pre-mRNA molecule, which amount to 13,505 amino acids from 38 proteins and a combined molecular mass of about 1.6 megadaltons. The 5' exon is anchored by loop I of U5 snRNA, whereas the 5' splice site (5'SS) and the branch-point sequence (BPS) of the intron are specifically recognized by U6 and U2 snRNA, respectively. Except for coordination of the catalytic metal ions, the RNA elements at the catalytic cavity of Prp8 are mostly primed for catalysis. The catalytic latency is maintained by the SF3b complex, which encircles the BPS, and the splicing factors Cwc24 and Prp11, which shield the 5' exon-5'SS junction. This structure, together with those determined earlier, outlines a molecular framework for the pre-mRNA splicing reaction.

Structure of a yeast catalytic step I spliceosome at 3.4 Å resolution

Each cycle of pre-messenger RNA splicing, carried out by the spliceosome, comprises two sequential transesterification reactions, which result in the removal of an intron and the joining of two exons. Here we report an atomic structure of a catalytic step I spliceosome (known as the C complex) from Saccharomyces cerevisiae, as determined by cryo-electron microscopy at an average resolution of 3.4 angstroms. In the structure, the 2'-OH of the invariant adenine nucleotide in the branch point sequence (BPS) is covalently joined to the phosphate at the 5' end of the 5' splice site (5'SS), forming an intron lariat. The freed 5' exon remains anchored to loop I of U5 small nuclear RNA (snRNA), and the 5'SS and BPS of the intron form duplexes with conserved U6 and U2 snRNA sequences, respectively. Specific placement of these RNA elements at the catalytic cavity of Prp8 is stabilized by 15 protein components, including Snu114 and the splicing factors Cwc21, Cwc22, Cwc25, and Yju2. These features, representing the conformation of the spliceosome after the first-step reaction, predict structural changes that are needed for the execution of the second-step transesterification reaction.

Anions in Nucleic Acid Crystallography

Nucleic acid crystallization buffers contain a large variety of chemicals fitting specific needs. Among them, anions are often solely considered for pH-regulating purposes and as cationic co-salts while their ability to directly bind to nucleic acid structures is rarely taken into account. Here we review current knowledge related to the use of anions in crystallization buffers along with data on their biological prevalence. Chloride ions are frequently identified in crystal structures but display low cytosolic concentrations. Hence, they are thought to be distant from nucleic acid structures in the cell. Sulfate ions are also frequently identified in crystal structures but their localization in the cell remains elusive. Nevertheless, the characterization of the binding properties of these ions is essential for better interpreting the solvent structure in crystals and consequently, avoiding mislabeling of electron densities. Furthermore, understanding the binding properties of these anions should help to get clues related to their potential effects in crowded cellular environments.

Using Molecular Replacement Phasing to Study the Structure and Function of RNA

In recent years a wide variety of RNA molecules regulating fundamental cellular processes has been discovered. Therefore, RNA structure determination is experiencing a boost and many more RNA structures are likely to be determined in the years to come. The broader availability of experimentally determined RNA structures implies that molecular replacement (MR) will be used more and more frequently as a method for phasing future crystallographic structures. In this report we describe various aspects relative to RNA structure determination by MR. First, we describe how to select and create MR search models for nucleic acids. Second, we describe how to perform MR searches on RNA using available crystallographic software. Finally, we describe how to refine and interpret the successful MR solutions. These protocols are applicable to determine novel RNA structures as well as to establish structural-functional relationships on existing RNA structures.

Sodium and Potassium Interactions with Nucleic Acids

Metal ions are essential cofactors for the structure and functions of nucleic acids. Yet, the early discovery in the 70s of the crucial role of Mg(2+) in stabilizing tRNA structures has occulted for a long time the importance of monovalent cations. Renewed interest in these ions was brought in the late 90s by the discovery of specific potassium metal ions in the core of a group I intron. Their importance in nucleic acid folding and catalytic activity is now well established. However, detection of K(+) and Na(+) ions is notoriously problematic and the question about their specificity is recurrent. Here we review the different methods that can be used to detect K(+) and Na(+) ions in nucleic acid structures such as X-ray crystallography, nuclear magnetic resonance or molecular dynamics simulations. We also discuss specific versus non-specific binding to different structures through various examples.

Spontaneous Formation of KCl Aggregates in Biomolecular Simulations: A Force Field Issue?

Realistic all-atom simulation of biological systems requires accurate modeling of both the biomolecules and their ionic environment. Recently, ion nucleation phenomena leading to the rapid growth of KCl or NaCl clusters in the vicinity of biomolecular systems have been reported. To better understand this phenomenon, molecular dynamics simulations of KCl aqueous solutions at three (1.0, 0.25, and 0.10 M) concentrations were performed. Two popular water models (TIP3P and SPC/E) and two Lennard-Jones parameter sets (AMBER and Dang) were combined to produce a total of 80 ns of molecular dynamics trajectories. Results suggest that the use of the Dang cation Lennard-Jones parameters instead of those adopted by the AMBER force-field produces a more accurate description of the ionic solution. In the later case, formation of salt aggregates is probably indicative of an artifact resulting from misbalanced force-field parameters. Because similar results were obtained with two different water parameter sets, the simulations exclude a water model dependency in the formation of anomalous ionic clusters. Overall, the results strongly suggest that for accurate modeling of ions in biomolecular systems, great care should be taken in choosing balanced ionic parameters even when using the most popular force-fields. These results invite a reexamination of older data obtained using available force-fields and a thorough check of the quality of current parameters sets by performing simulations at finite (>0.25 M) instead of minimal salt conditions.

Landscape of π-π and sugar-π contacts in DNA-protein interactions

There were 1765 contacts identified between DNA nucleobases or deoxyribose and cyclic (W, H, F, Y) or acyclic (R, E, D) amino acids in 672 X-ray structures of DNA-protein complexes. In this first study to compare π-interactions between the cyclic and acyclic amino acids, visual inspection was used to categorize amino acid interactions as nucleobase π-π (according to biological edge) or deoxyribose sugar-π (according to sugar edge). Overall, 54% of contacts are nucleobase π-π interactions, which involve all amino acids, but are more common for Y, F, and R, and involve all DNA nucleobases with similar frequencies. Among binding arrangements, cyclic amino acids prefer more planar (stacked) π-systems than the acyclic counterparts. Although sugar-π interactions were only previously identified with the cyclic amino acids and were found to be less common (38%) than nucleobase-cyclic amino acid contacts, sugar-π interactions are more common than nucleobase π-π contacts for the acyclic series (61% of contacts). Similar to DNA-protein π-π interactions, sugar-π contacts most frequently involve Y and R, although all amino acids adopt many binding orientations relative to deoxyribose. These DNA-protein π-interactions stabilize biological systems, by up to approximately -40 kJ mol(-1) for neutral nucleobase or sugar-amino acid interactions, but up to approximately -95 kJ mol(-1) for positively or negatively charged contacts. The high frequency and strength, despite variation in structure and composition, of these π-interactions point to an important function in biological systems.

Multiscale methods for computational RNA enzymology

RNA catalysis is of fundamental importance to biology and yet remains ill-understood due to its complex nature. The multidimensional "problem space" of RNA catalysis includes both local and global conformational rearrangements, changes in the ion atmosphere around nucleic acids and metal ion binding, dependence on potentially correlated protonation states of key residues, and bond breaking/forming in the chemical steps of the reaction. The goal of this chapter is to summarize and apply multiscale modeling methods in an effort to target the different parts of the RNA catalysis problem space while also addressing the limitations and pitfalls of these methods. Classical molecular dynamics simulations, reference interaction site model calculations, constant pH molecular dynamics (CpHMD) simulations, Hamiltonian replica exchange molecular dynamics, and quantum mechanical/molecular mechanical simulations will be discussed in the context of the study of RNA backbone cleavage transesterification. This reaction is catalyzed by both RNA and protein enzymes, and here we examine the different mechanistic strategies taken by the hepatitis delta virus ribozyme and RNase A.

Like-charge ion pairs of hydronium and hydroxide in aqueous solution?

Hydroxide can form stronger like-ion pairs than hydronium in aqueous solution mostly due to its versatile coordination ability with solvents.

Ion distributions around left- and right-handed DNA and RNA duplexes: a comparative study

The ion atmosphere around nucleic acids is an integral part of their solvated structure. However, detailed aspects of the ionic distribution are difficult to probe experimentally, and comparative studies for different structures of the same sequence are almost non-existent. Here, we have used large-scale molecular dynamics simulations to perform a comparative study of the ion distribution around (5'-CGCGCGCGCGCG-3')2 dodecamers in solution in B-DNA, A-RNA, Z-DNA and Z-RNA forms. The CG sequence is very sensitive to ionic strength and it allows the comparison with the rare but important left-handed forms. The ions investigated include Na(+), K(+) and Mg(2 +), with various concentrations of their chloride salts. Our results quantitatively describe the characteristics of the ionic distributions for different structures at varying ionic strengths, tracing these differences to nucleic acid structure and ion type. Several binding pockets with rather long ion residence times are described, both for the monovalent ions and for the hexahydrated Mg[(H2O)6](2+) ion. The conformations of these binding pockets include direct binding through desolvated ion bridges in the GpC steps in B-DNA and A-RNA; direct binding to backbone oxygens; binding of Mg[(H2O)6](2+) to distant phosphates, resulting in acute bending of A-RNA; tight 'ion traps' in Z-RNA between C-O2 and the C-O2' atoms in GpC steps; and others.

Structural basis for activity of highly efficient RNA mimics of green fluorescent protein

GFP and its derivatives revolutionized the study of proteins. Spinach is a recently reported in vitro-evolved RNA mimic of GFP, which as genetically encoded fusions makes possible live-cell, real-time imaging of biological RNAs without resorting to large RNA-binding protein-GFP fusions. To elucidate the molecular basis of Spinach fluorescence, we solved the cocrystal structure of Spinach bound to its cognate exogenous chromophore, showing that Spinach activates the small molecule by immobilizing it between a base triple, a G-quadruplex and an unpaired G. Mutational and NMR analyses indicate that the G-quadruplex is essential for Spinach fluorescence, is also present in other fluorogenic RNAs and may represent a general strategy for RNAs to induce fluorescence of chromophores. The structure guided the design of a miniaturized 'Baby Spinach', and it provides a foundation for structure-driven design and tuning of fluorescent RNAs.

Validating fragment-based drug discovery for biological RNAs: lead fragments bind and remodel the TPP riboswitch specifically

Thiamine pyrophosphate (TPP) riboswitches regulate essential genes in bacteria by changing conformation upon binding intracellular TPP. Previous studies using fragment-based approaches identified small molecule "fragments" that bind this gene-regulatory mRNA domain. Crystallographic studies now show that, despite having micromolar Kds, four different fragments bind the TPP riboswitch site-specifically, occupying the pocket that recognizes the aminopyrimidine of TPP. Unexpectedly, the unoccupied site that would recognize the pyrophosphate of TPP rearranges into a structure distinct from that of the cognate complex. This idiosyncratic fragment-induced conformation, also characterized by small-angle X-ray scattering and chemical probing, represents a possible mechanism for adventitious ligand discrimination by the riboswitch, and suggests that off-pathway conformations of RNAs can be targeted for drug development. Our structures, together with previous screening studies, demonstrate the feasibility of fragment-based drug discovery against RNA targets.

Principles of ion recognition in RNA: insights from the group II intron structures

Metal ions promote both RNA folding and catalysis, thus being essential in stabilizing the structure and determining the function of large RNA molecules, including group II introns. The latter are self-splicing metalloribozymes, containing a heteronuclear four-metal-ion center within the active site. In addition to these catalytic ions, group II introns bind many other structural ions, including delocalized ions that bind the RNA diffusively and well-ordered ions that bind the RNA tightly with high occupancy. The latter ions, which can be studied by biophysical methods, have not yet been analyzed systematically. Here, we compare crystal structures of the group IIC intron from Oceanobacillus iheyensis and classify numerous site-bound ions, which are primarily localized in the intron core and near long-range tertiary contacts. Certain ion-binding sites resemble motifs observed in known RNA structures, while others are idiosyncratic to the group II intron. Particularly interesting are (1) ions proximal to the active site, which may participate in splicing together with the catalytic four-metal-ion center, (2) organic ions that bind regions predicted to interact with intron-encoded proteins, and (3) unusual monovalent ions bound to GU wobble pairs, GA mismatches, the S-turn, the tetraloop-receptor, and the T-loop. Our analysis extends the general principles by which ions participate in RNA structural organization and it will aid in the determination and interpretation of future RNA structures.

Themes and variations in riboswitch structure and function

The complexity of gene expression control by non-coding RNA has been highlighted by the recent progress in the field of riboswitches. Discovered a decade ago, riboswitches represent a diverse group of non-coding mRNA regions that possess a unique ability to directly sense cellular metabolites and modulate gene expression through formation of alternative metabolite-free and metabolite-bound conformations. Such protein-free metabolite sensing domains utilize sophisticated three-dimensional folding of RNA molecules to discriminate between a cognate ligand from related compounds so that only the right ligand would trigger a genetic response. Given the variety of riboswitch ligands ranging from small cations to large coenzymes, riboswitches adopt a great diversity of structures. Although many riboswitches share structural principles to build metabolite-competent folds, form precise ligand-binding pockets, and communicate a ligand-binding event to downstream regulatory regions, virtually all riboswitch classes possess unique features for ligand recognition, even those tuned to recognize the same metabolites. Here we present an overview of the biochemical and structural research on riboswitches with a major focus on common principles and individual characteristics adopted by these regulatory RNA elements during evolution to specifically target small molecules and exert genetic responses. This article is part of a Special Issue entitled: Riboswitches.

Tautomerisation of thymine acts against the Hückel 4N + 2 rule. The effect of metal ions and H-bond complexations on the electronic structure of thymine

Not necessarily the π-electron delocalization is responsible for the stability of thymine tautomers.

Significant strength of charged DNA-protein π-π interactions: a preliminary study of cytosine

The present work characterized the preferred gas-phase structure and optimum interaction energy of both parallel stacked and perpendicular T-shaped dimers between cytosine (C), as a representative nucleobase, and aspartic/glutamic acid (DE), aspartate/glutamate (DE(-)) or arginine (R(+)), using detailed M06-2X/6-31+G(d,p) potential energy surface scans as a function of the relative monomer orientation. Through comparison to previous literature on the π-π interactions between the DNA nucleobases and the aromatic amino acid residues, this work will allow for comparisons between DNA-protein interactions involving aromatic and acyclic R-side chains, as well as comparisons of the relative geometric dependence and magnitude of π-π (C:DE), πcation-π (C:R(+)), and πanion-π (C:DE(-)) interactions. Our results show that the preferred relative monomer orientation is highly dependent on the monomer composition and charge, and is dictated by electrostatic-driven interactions. More importantly, for the first time, we report that the π-π interactions between cytosine and (neutral) aspartic/glutamic acid are up to approximately -40 kJ mol(-1), while the πcation-π or πanion-π interactions between cytosine and arginine or aspartate/glutamate are up to approximately -90 and -99 kJ mol(-1), respectively. An extensive investigation of the effects of the computational methodology implemented, including comparisons to detailed CCSD(T)/CBS potential energy surfaces and interaction energies, supports the use of M06-2X, as well as ωB97X-D, to study DNA-protein π-π interactions of varying composition and charge. Most importantly, the CCSD(T)/CBS results verify the strong nature of these DNA-protein π-π interactions, as well as the unique nature of the πcation-π and πanion-π counterparts. Therefore, our results emphasize that a wide variety of different types of noncovalent interactions between both cyclic and acyclic π-containing components can significantly contribute to the stability of DNA-protein complexes and likely play a larger role in biology than currently accepted.

DNA self-assembly: from chirality to evolution

Transient or long-term DNA self-assembly participates in essential genetic functions. The present review focuses on tight DNA-DNA interactions that have recently been found to play important roles in both controlling DNA higher-order structures and their topology. Due to their chirality, double helices are tightly packed into stable right-handed crossovers. Simple packing rules that are imposed by DNA geometry and sequence dictate the overall architecture of higher order DNA structures. Close DNA-DNA interactions also provide the missing link between local interactions and DNA topology, thus explaining how type II DNA topoisomerases may sense locally the global topology. Finally this paper proposes that through its influence on DNA self-assembled structures, DNA chirality played a critical role during the early steps of evolution.

What stabilizes close arginine pairing in proteins?

Close stacking of arginine residues are often observed in protein structures despite the highly repulsive nature of the close like-charged groups. Physical factors stabilizing the close guanidinium ions of arginine side-chains have been previously studied in water and in protein-like environments, and the hydration free energy has been emphasized to be an important factor. However, how close arginine pairs are stabilized in real proteins has not been fully understood yet. In this paper, we show that arginine pairs are more frequently found in the protein interior than expected from the frequency of unpaired arginines buried inside protein through a statistical analysis of the protein structure database. We then confirm that 4 selected arginine pairs buried in the protein are indeed positively charged rather than neutralized, by molecular dynamics simulations and pKa estimation with molecular mechanics-Poisson-Boltzmann calculations. Further energy decomposition analysis shows that the hydration free energy may not be strong enough to overcome the repulsive Coulomb interaction between the positively charged arginine residues buried inside the protein. Instead, a highly polar interaction network is identified around each buried arginine pair, and the electrostatic interactions within such network are strong enough to stabilize the repulsive interaction of the buried arginine pair for the 4 selected cases. The polar interaction network is highly conserved evolutionarily in some proteins, implicating their roles in protein stabilization or biochemical function.

Structural principles of RNA catalysis in a 2'-5' lariat-forming ribozyme

RNA-catalyzed lariat formation is present in both eukaryotes and prokaryotes. To date we lack structural insights into the catalytic mechanism of lariat-forming ribozymes. Here, we study an artificial 2'-5' AG1 lariat-forming ribozyme that shares the sequence specificity of lariat formation with the pre-mRNA splicing reaction. Using NMR, we solve the structure of the inactive state of the ribozyme in the absence of magnesium. The reaction center 5'-guanosine appears to be part of a helix with an exceptionally widened major groove, while the lariat-forming A48 is looped out at the apex of a pseudoknot. The model of the active state built by mutational analysis, molecular modeling, and small-angle X-ray scattering suggests that A48 is recognized by a conserved adenosine, juxtaposed to the 5'-guanosine in one base-pair step distance, while the G1-N7 coordinates a magnesium ion essential for the activation of the nucleophile. Our findings offer implications for lariat formation in RNA enzymes including the mechanism of the recognition of the branch-site adenosine.