Calculations of pKa Values for a Series of Naturally Occurring Modified Nucleobases

Antisense RNA C9orf72 hexanucleotide repeat associated with amyotrophic lateral sclerosis and frontotemporal dementia forms a triplex-like structure and binds small synthetic ligand

The abnormal expansion of GGGGCC/GGCCCC hexanucleotide repeats (HR) in C9orf72 is associated with amyotrophic lateral sclerosis (ALS) and frontotemporal dementia (FTD). Structural polymorphisms of HR result in the multifactorial pathomechanism of ALS/FTD. Consequently, many ongoing studies are focused at developing therapies targeting pathogenic HR RNA. One of them involves small molecules blocking sequestration of important proteins, preventing formation of toxic nuclear foci. However, rational design of potential therapeutics is hindered by limited number of structural studies of RNA-ligand complexes. We determined the crystal structure of antisense HR RNA in complex with ANP77 ligand (1.1 Å resolution) and in the free form (0.92 and 1.5 Å resolution). HR RNA folds into a triplex structure composed of four RNA chains. ANP77 interacted with two neighboring single-stranded cytosines to form pseudo-canonical base pairs by adopting sandwich-like conformation and adjusting the position of its naphthyridine units to the helical twist of the RNA. In the unliganded structure, the cytosines formed a peculiar triplex i-motif, assembled by trans C•C+ pair and a third cytosine located at the Hoogsteen edge of the C•C+ pair. These results extend our knowledge of the structural polymorphisms of HR and can be used for rational design of small molecules targeting disease-related RNAs.

Nanopores map the acid-base properties of a single site in a single DNA molecule

Nanopores are increasingly powerful tools for single molecule sensing, in particular, for sequencing DNA, RNA and peptides. This success has spurred efforts to sequence non-canonical nucleic acid bases and amino acids. While canonical DNA and RNA bases have pKas far from neutral, certain non-canonical bases, natural RNA modifications, and amino acids are known to have pKas near neutral pHs at which nanopore sequencing is typically performed. Previous reports have suggested that the nanopore signal may be sensitive to the protonation state of an individual moiety. We sequenced ion currents with the MspA nanopore using a single stranded DNA containing a single non-canonical DNA base (Z) at various pH conditions. The Z-base has a near-neutral pKa ∼ 7.8. We find that the measured ion current is remarkably sensitive to the protonation state of the Z-base. We demonstrate how nanopores can be used to localize and determine the pKa of individual moieties along a polymer. More broadly, these experiments provide a path to mapping different protonation sites along polymers and give insight in how to optimize sequencing of polymers that contain moieties with near-neutral pKas.

Assessment of Nucleobase Protomeric and Tautomeric States in Nucleic Acid Structures for Interaction Analysis and Structure-Based Ligand Design

With increasing interest in RNA as a therapeutic and a potential target, the role of RNA structures has become more important. Even slight changes in nucleobases, such as modifications or protomeric and tautomeric states, can have a large impact on RNA structure and function, while local environments in turn affect protonation and tautomerization. In this work, the application of empirical tools for pKa and tautomer prediction for RNA modifications was elucidated and compared with ab initio quantum mechanics (QM) methods and expanded toward macromolecular RNA structures, where QM is no longer feasible. In this regard, the Protonate3D functionality within the molecular operating environment (MOE) was expanded for nucleobase protomer and tautomer predictions and applied to reported examples of altered protonation states depending on the local environment. Overall, observations of nonstandard protomers and tautomers were well reproduced, including structural C+G:C(A) and A+GG motifs, several mismatches, and protonation of adenosine or cytidine as the general acid in nucleolytic ribozymes. Special cases, such as cobalt hexamine-soaked complexes or the deprotonation of guanosine as the general base in nucleolytic ribozymes, proved to be challenging. The collected set of examples shall serve as a starting point for the development of further RNA protonation prediction tools, while the presented Protonate3D implementation already delivers reasonable protonation predictions for RNA and DNA macromolecules. For cases where higher accuracy is needed, like following catalytic pathways of ribozymes, incorporation of QM-based methods can build upon the Protonate3D-generated starting structures. Likewise, this protonation prediction can be used for structure-based RNA-ligand design approaches.

Atomic Mutagenesis of N6-Methyladenosine Reveals Distinct Recognition Modes by Human m6A Reader and Eraser Proteins

N6-methyladenosine (m6A) is an important modified nucleoside in cellular RNA associated with multiple cellular processes and is implicated in diseases. The enzymes associated with the dynamic installation and removal of m6A are heavily investigated targets for drug research, which requires detailed knowledge of the recognition modes of m6A by proteins. Here, we use atomic mutagenesis of m6A to systematically investigate the mechanisms of the two human m6A demethylase enzymes FTO and ALKBH5 and the binding modes of YTH reader proteins YTHDF2/DC1/DC2. Atomic mutagenesis refers to atom-specific changes that are introduced by chemical synthesis, such as the replacement of nitrogen by carbon atoms. Synthetic RNA oligonucleotides containing site-specifically incorporated 1-deaza-, 3-deaza-, and 7-deaza-m6A nucleosides were prepared by solid-phase synthesis and their RNA binding and demethylation by recombinant proteins were evaluated. We found distinct differences in substrate recognition and transformation and revealed structural preferences for the enzymatic activity. The deaza m6A analogues introduced in this work will be useful probes for other proteins in m6A research.

Advances in the Structural and Functional Understanding of m1A RNA Modification

ConspectusRNA modification is a co- or post-transcriptional process by which specific nucleotides are chemically altered by enzymes after their initial incorporation into the RNA chain, expanding the chemical and functional diversity of RNAs. Our understanding of RNA modifications has changed dramatically in recent years. In the past decade, RNA methyltransferases (MTases) have been highlighted in numerous clinical studies and disease models, modifications have been found to be dynamically regulated by demodification enzymes, and significant technological advances have been made in the fields of RNA sequencing, mass spectrometry, and structural biology. Among RNAs, transfer RNAs (tRNAs) exhibit the greatest diversity and density of post-transcriptional modifications, which allow for potential cross-talks and regulation during their incorporation. N1-methyladenosine (m1A) modification is found in tRNAs at positions 9, 14, 16, 22, 57, and 58, depending on the tRNA and organism.Our laboratory has used and developed a large panel of tools to decipher the different mechanisms used by m1A tRNA MTases to recognize and methylate tRNA. We have solved the structures of TrmI from Thermus thermophilus (m1A58), TrmK from Bacillus subtilis (m1A22), and human TRMT10C (m1A9). These MTases do not share the same structure or organization to recognize tRNAs, but they all modify an adenosine, forming a non-Watson-Crick (WC) interaction. For TrmK, nuclear magnetic resonance (NMR) chemical shift mapping of the binding interface between TrmK and tRNASer was invaluable to build a TrmK/tRNA model, where both domains of TrmK participate in the binding of a full-length L-shaped tRNA and where the non-WC purine 13-A22 base pair positions the A22 N1-atom close to the methyl of the S-adenosyl-l-methionine (SAM) TrmK cofactor. For TRMT10C, cryoEM structures showed the MTase poised to N1-methylate A9 or G9 in tRNA and revealed different steps of tRNA maturation, where TRMT10C acts as a tRNA binding platform for sequential docking of each maturation enzyme. This work confers a role for TRMT10C in tRNA quality control and provides a framework to understand the link between mitochondrial tRNA maturation dysfunction and diseases.Methods to directly detect the incorporation of modifications during tRNA biosynthesis are rare and do not provide easy access to the temporality of their introduction. To this end, we have introduced time-resolved NMR to monitor tRNA maturation in the cellular environment. Combined with genetic and biochemical approaches involving the synthesis of specifically modified tRNAs, our methodology revealed that some modifications are incorporated in a defined sequential order, controlled by cross-talks between modification events. In particular, a strong modification circuit, namely Ψ55 → m5U54 → m1A58, controls the modification process in the T-arm of yeast elongator tRNAs. Conversely, we showed that m1A58 is efficiently introduced on unmodified initiator tRNAiMet without the need of any prior modification. Two distinct pathways are therefore followed for m1A58 incorporation in elongator and initiator tRNAs.We are undoubtedly entering an exciting period for the elucidation of the functions of RNA modifications and the intricate mechanisms by which modification enzymes identify and alter their RNA substrates. These are promising directions for the field of epitranscriptomics.

General Principles and Limitations for Detection of RNA Modifications by Sequencing

Among the many analytical methods applied to RNA modifications, a particularly pronounced surge has occurred in the past decade in the field of modification mapping. The occurrence of modifications such as m6A in mRNA, albeit known since the 1980s, became amenable to transcriptome-wide analyses through the advent of next-generation sequencing techniques in a rather sudden manner. The term "mapping" here refers to detection of RNA modifications in a sequence context, which has a dramatic impact on the interpretation of biological functions. As a consequence, an impressive number of mapping techniques were published, most in the perspective of what now has become known as "epitranscriptomics". While more and more different modifications were reported to occur in mRNA, conflicting reports and controversial results pointed to a number of technical and theoretical problems rooted in analytics, statistics, and reagents. Rather than finding the proverbial needle in a haystack, the tasks were to determine how many needles of what color in what size of a haystack one was looking at.As the authors of this Account, we think it important to outline the limitations of different mapping methods since many life scientists freshly entering the field confuse the accuracy and precision of modification mapping with that of normal sequencing, which already features numerous caveats by itself. Indeed, we propose here to qualify a specific mapping method by the size of the transcriptome that can be meaningfully analyzed with it.We here focus on high throughput sequencing by Illumina technology, referred to as RNA-Seq. We noted with interest the development of methods for modification detection by other high throughput sequencing platforms that act directly on RNA, e.g., PacBio SMRT and nanopore sequencing, but those are not considered here.In contrast to approaches relying on direct RNA sequencing, current Illumina RNA-Seq protocols require prior conversion of RNA into DNA. This conversion relies on reverse transcription (RT) to create cDNA; thereafter, the cDNA undergoes a sequencing-by-synthesis type of analysis. Thus, a particular behavior of RNA modified nucleotides during the RT-step is a prerequisite for their detection (and quantification) by deep sequencing, and RT properties have great influence on the detection efficiency and reliability. Moreover, the RT-step requires annealing of a synthetic primer, a prerequisite with a crucial impact on library preparation. Thus, all RNA-Seq protocols must feature steps for the introduction of primers, primer landing sites, or adapters on both the RNA 3'- and 5'-ends.

Characterizing Benzo[a]pyrene Adducts in Transfer RNAs Using Liquid Chromatography Coupled with Tandem Mass Spectrometry (LC-MS/MS)

The activated forms of the environmental pollutant benzo[a]pyrene (B[a]P), such as benzo[a]pyrene diol epoxide (BPDE), are known to cause damage to genomic DNA and proteins. However, the impact of BPDE on ribonucleic acid (RNA) remains unclear. To understand the full spectrum of potential BPDE-RNA adducts formed, we reacted ribonucleoside standards with BPDE and characterized the reaction products using liquid chromatography coupled with tandem mass spectrometry (LC-MS/MS). To understand the potential types of adducts that could form with biological RNAs, eukaryotic transfer RNAs (tRNAs) were also reacted with BPDE. The isolation and analysis of the modified and adducted ribonucleosides using LC-MS/MS revealed several BPDE derivatives of post-transcriptional modifications. The approach outlined in this work enables the identification of RNA adducts from BPDE, which can pave the way for understanding the potential impacts of such adducts on the higher-order structure and function of modified RNAs.

Calculated pKa Values for a Series of Aza- and Deaza-Modified Nucleobases

A variety of synthetic modified nucleobases have been used to investigate the structure and function of RNA and DNA or act as enzyme inhibitors. A set of these modifications involves the addition or removal of a nitrogen atom in the ring. These aza and deaza modifications have garnered interest as useful biochemical tools, but information on some of their physical characteristics is lacking. In this study, the B3LYP density functional with the 6-31+G(d,p) basis set and an implicit-explicit solvent model was used to perform ab initio quantum mechanical studies to estimate pK_a values of aza- and deaza-modified nucleobases. A comparison between theoretical and known experimental pK_a values was carried out, and adjustment factors were applied to 57 pK_a values in the purine and pyrimidine data sets.

Base Modifications of Nucleosides via the Use of Peptide-Coupling Agents, and Beyond

Several naturally occurring purine and pyrimidine nucleosides contain an amide linkage as part of the heterocyclic aglycone. Enolization of the amide and conversion to leaving groups at the amide carbon atom permits base modification by addition-elimination types of processes. Although a number of methods have been developed over the years for accomplishing such conversions, the present Personal Account describes efforts from the Lakshman laboratories. Facile activation of the amido groups in nucleobases can be achieved with peptide-coupling agents. Subsequent reaction with nucleophiles then accomplishes the base modifications. In many cases, the activation and displacement steps can be done as two-step, one-pot processes, whereas in other cases, discrete storable activated nucleosides can be isolated for subsequent displacement reactions. Using such an approach a wide range of nucleoside base modifications is readily achievable. In many instances, mechanistic investigations have been conducted so as to understand the activation process.

Scaling Catalytic Contributions of Small Self-Cleaving Ribozymes

Nucleolytic ribozymes utilize general acid-base catalysis to perform phosphodiester cleavage. In most ribozyme classes, a conserved active site guanosine is positioned to act as general base, thereby activating the 2'-OH group to attack the scissile phosphate (γ-catalysis). Here, we present an atomic mutagenesis study for the pistol ribozyme class. Strikingly, "general base knockout" by replacement of the guanine N1 atom by carbon results in only 2.7-fold decreased rate. Therefore, the common view that γ-catalysis critically depends on the N1 moiety becomes challenged. For pistol ribozymes we found that γ-catalysis is subordinate in overall catalysis, made up by two other catalytic factors (α and δ). Our approach allows scaling of the different catalytic contributions (α, β, γ, δ) with unprecedented precision and paves the way for a thorough mechanistic understanding of nucleolytic ribozymes with active site guanines.

Scaling Catalytic Contributions of Small Self-Cleaving Ribozymes

Nucleolytic ribozymes utilize general acid-base catalysis to perform phosphodiester cleavage. In most ribozyme classes, a conserved active site guanosine is positioned to act as general base, thereby activating the 2'-OH group to attack the scissile phosphate (γ-catalysis). Here, we present an atomic mutagenesis study for the pistol ribozyme class. Strikingly, "general base knockout" by replacement of the guanine N1 atom by carbon results in only 2.7-fold decreased rate. Therefore, the common view that γ-catalysis critically depends on the N1 moiety becomes challenged. For pistol ribozymes we found that γ-catalysis is subordinate in overall catalysis, made up by two other catalytic factors (α and δ). Our approach allows scaling of the different catalytic contributions (α, β, γ, δ) with unprecedented precision and paves the way for a thorough mechanistic understanding of nucleolytic ribozymes with active site guanines.

Three's a crowd - stabilisation, structure, and applications of DNA triplexes

DNA is a strikingly flexible molecule and can form a variety of secondary structures, including the triple helix, which is the subject of this review. The DNA triplex may be formed naturally, during homologous recombination, or can be formed by the introduction of a synthetic triplex forming oligonucleotide (TFO) to a DNA duplex. As the TFO will bind to the duplex with sequence specificity, there is significant interest in developing TFOs with potential therapeutic applications, including using TFOs as a delivery mechanism for compounds able to modify or damage DNA. However, to combine triplexes with functionalised compounds, a full understanding of triplex structure and chemical modification strategies, which may increase triplex stability or in vivo degradation, is essential - these areas will be discussed in this review. Ruthenium polypyridyl complexes, which are able to photooxidise DNA and act as luminescent DNA probes, may serve as a suitable photophysical payload for a TFO system and the developments in this area in the context of DNA triplexes will also be reviewed.

Protein-Based Virtual Screening Tools Applied for RNA-Ligand Docking Identify New Binders of the preQ1-Riboswitch

Targeting RNA with small molecules is an emerging field. While several ligands for different RNA targets are reported, structure-based virtual screenings (VSs) against RNAs are still rare. Here, we elucidated the general capabilities of protein-based docking programs to reproduce native binding modes of small-molecule RNA ligands and to discriminate known binders from decoys by the scoring function. The programs were found to perform similar compared to the RNA-based docking tool rDOCK, and the challenges faced during docking, namely, protomer and tautomer selection, target dynamics, and explicit solvent, do not largely differ from challenges in conventional protein-ligand docking. A prospective VS with the Bacillus subtilis preQ₁-riboswitch aptamer domain performed with FRED, HYBRID, and FlexX followed by microscale thermophoresis assays identified six active compounds out of 23 tested VS hits with potencies between 29.5 nM and 11.0 μM. The hits were selected not solely based on their docking score but for resembling key interactions of the native ligand. Therefore, this study demonstrates the general feasibility to perform structure-based VSs against RNA targets, while at the same time it highlights pitfalls and their potential solutions when executing RNA-ligand docking.

DNA with Purine-Purine Base Pairs: Size and Position of Isoguanine and 8-Aza-7-deazaisoguanine Clickable Residues Control the Molecular Recognition of Guanine and 5-Aza-7-deazaguanine

Purine-purine base pairs represent an alternative recognition system to the purine-pyrimidine pairing reported by Watson and Crick. Modified purines are the source for non-canonical interactions. To mimic dG-dC interactions, 2'-deoxyisoguanosine (1a) and 8-aza-7-deaza-2'-deoxyisoguanosine (2a) are used to construct base pairs with 2'-deoxyguanosine or 5-aza-7-deaza-2'-deoxyguanosine (dZ). This work reports the chemical functionalization of 1a and its shape mimic 2a in purine-purine base pairs. Clickable rigid ethynyl and more flexible octadiynyl side chain derivatives of 1a and 2a were synthesized. They were protected and converted into phosphoramidites. Building blocks were employed in the synthesis of base-modified 12-mer oligonucleotides with clickable side chains. Pyrene azide was clicked to the linkers. After hybridization, oligonucleotides with purine-purine base pairs were constructed with linkers and pyrene adducts at position-8 of isoguanine and at position-7 of 8-aza-7-deazaisoguanine. Recognition and stability of purine-purine base pairs were explored using T_m values, thermodynamic data, and CD-spectroscopic changes. Side chains at position-7 of 8-aza-7-deazaisoguanine-guanine base pairs or with 5-aza-7-deazaguanine are well accommodated in DNA, whereas functionalization at 8-position of isoguanine makes the DNA unstable. Pyrene click adducts verified the observation. In conclusion, position-7 is the place of choice for purine-purine base pair functionalization.