An innate twist between Crick's wobble and Watson-Crick base pairs

RNA thermometers are widespread upstream of ABC transporter genes in bacteria

RNA thermometers are temperature-sensing non-coding RNAs that regulate the expression of downstream genes. A well-characterized RNA thermometer motif discovered in bacteria is the ROSE-like element (repression of heat shock gene expression). ATP-binding cassette (ABC) transporters are a superfamily of transmembrane proteins that harness ATP hydrolysis to facilitate the export and import of substrates across cellular membranes. Through structure-guided bioinformatics, we discovered that ROSE-like RNA thermometers are widespread upstream of ABC transporter genes in bacteria. X-ray crystallography, biochemistry, and cellular assays indicate that these RNA thermometers are functional regulatory elements. This study expands the known biological role of RNA thermometers to these key membrane transporters.

Clinical, genetic and structural delineation of RPL13-related spondyloepimetaphyseal dysplasia suggest extra-ribosomal functions of eL13

Spondyloepimetaphyseal dysplasia with severe short stature, RPL13-related (SEMD-RPL13), MIM#618728), is a rare autosomal dominant disorder characterized by short stature and skeletal changes such as mild spondylar and epimetaphyseal dysplasia affecting primarily the lower limbs. The genetic cause was first reported in 2019 by Le Caignec et al., and six disease-causing variants in the gene coding for a ribosomal protein, RPL13 (NM_000977.3) have been identified to date. This study presents clinical and radiographic data from 12 affected individuals aged 2-64 years from seven unrelated families, showing highly variable manifestations. The affected individuals showed a range from mild to severe short stature, retaining the same radiographic pattern of spondylar- and epi-metaphyseal dysplasia, but with varying severity of the hip and knee deformities. Two new missense variants, c.548 G>A, p.(Arg183His) and c.569 G>T, p.(Arg190Leu), and a previously known splice variant c.477+1G>A were identified, confirming mutational clustering in a highly specific RNA binding motif. Structural analysis and interpretation of the variants' impact on the protein suggests that disruption of extra-ribosomal functions of the protein through binding of mRNA may play a role in the skeletal phenotype of SEMD-RPL13. In addition, we present gonadal and somatic mosaicism for the condition.

G·U base pairing motifs in long non-coding RNAs

Long non-coding RNAs (lncRNAs) are recently-discovered transcripts involved in gene expression regulation and associated with diseases. Despite the unprecedented molecular complexity of these transcripts, recent studies of the secondary and tertiary structure of lncRNAs are starting to reveal the principles of lncRNA structural organization, with important functional implications. It therefore starts to be possible to analyze lncRNA structures systematically. Here, using a set of prototypical and medically-relevant lncRNAs of known secondary structure, we specifically catalogue the distribution and structural environment of one of the first-identified and most frequently occurring non-canonical Watson-Crick interactions, the G·U base pair. We compare the properties of G·U base pairs in our set of lncRNAs to those of the G·U base pairs in other well-characterized transcripts, like rRNAs, tRNAs, ribozymes, and riboswitches. Furthermore, we discuss how G·U base pairs in these targets participate in establishing interactions with proteins or miRNAs, and how they enable lncRNA tertiary folding by forming intramolecular or metal-ion interactions. Finally, by identifying highly-G·U-enriched regions of yet unknown function in our target lncRNAs, we provide a new rationale for future experimental investigation of these motifs, which will help obtain a more comprehensive understanding of lncRNA functions and molecular mechanisms in the future.

The structure and mechanism of action of a distinct class of dicistrovirus intergenic region IRESs

Internal ribosomal entry sites (IRESs) engage with the eukaryotic translation apparatus to promote end-independent initiation. We identified a conserved class of ∼150 nt long intergenic region (IGR) IRESs in dicistrovirus genomes derived from members of the phyla Arthropoda, Bryozoa, Cnidaria, Echinodermata, Entoprocta, Mollusca and Porifera. These IRESs, exemplified by Wenling picorna-like virus 2, resemble the canonical cricket paralysis virus (CrPV) IGR IRES in comprising two nested pseudoknots (PKII/PKIII) and a 3'-terminal pseudoknot (PKI) that mimics a tRNA anticodon stem-loop base-paired to mRNA. However, they are ∼50 nt shorter than CrPV-like IRESs, and PKIII is an H-type pseudoknot that lacks the SLIV and SLV stem-loops that are primarily responsible for the affinity of CrPV-like IRESs for the 40S ribosomal subunit and that restrict initial binding of PKI to its aminoacyl (A) site. Wenling-class IRESs bound strongly to 80S ribosomes but only weakly to 40S subunits. Whereas CrPV-like IRESs must be translocated from the A site to the peptidyl (P) site by elongation factor 2 for elongation to commence, Wenling-class IRESs bound directly to the P site of 80S ribosomes, and decoding begins without a prior translocation step. A chimeric CrPV clone containing a Wenling-class IRES was infectious, confirming that the IRES functioned in cells.

Alu RNA fold links splicing with signal recognition particle proteins

Transcriptomic diversity in primates was considerably expanded by exonizations of intronic Alu elements. To better understand their cellular mechanisms we have used structure-based mutagenesis coupled with functional and proteomic assays to study the impact of successive primate mutations and their combinations on inclusion of a sense-oriented AluJ exon in the human F8 gene. We show that the splicing outcome was better predicted by consecutive RNA conformation changes than by computationally derived splicing regulatory motifs. We also demonstrate an involvement of SRP9/14 (signal recognition particle) heterodimer in splicing regulation of Alu-derived exons. Nucleotide substitutions that accumulated during primate evolution relaxed the conserved left-arm AluJ structure including helix H1 and reduced the capacity of SRP9/14 to stabilize the closed Alu conformation. RNA secondary structure-constrained mutations that promoted open Y-shaped conformations of the Alu made the Alu exon inclusion reliant on DHX9. Finally, we identified additional SRP9/14 sensitive Alu exons and predicted their functional roles in the cell. Together, these results provide unique insights into architectural elements required for sense Alu exonization, identify conserved pre-mRNA structures involved in exon selection and point to a possible chaperone activity of SRP9/14 outside the mammalian signal recognition particle.

Structural analysis of mitochondrial rRNA gene variants identified in patients with deafness

The last few years have witnessed dramatic advances in our understanding of the structure and function of the mammalian mito-ribosome. At the same time, the first attempts to elucidate the effects of mito-ribosomal fidelity (decoding accuracy) in disease have been made. Hence, the time is right to push an important frontier in our understanding of mitochondrial genetics, that is, the elucidation of the phenotypic effects of mtDNA variants affecting the functioning of the mito-ribosome. Here, we have assessed the structural and functional role of 93 mitochondrial (mt-) rRNA variants thought to be associated with deafness, including those located at non-conserved positions. Our analysis has used the structural description of the human mito-ribosome of the highest quality currently available, together with a new understanding of the phenotypic manifestation of mito-ribosomal-associated variants. Basically, any base change capable of inducing a fidelity phenotype may be considered non-silent. Under this light, out of 92 previously reported mt-rRNA variants thought to be associated with deafness, we found that 49 were potentially non-silent. We also dismissed a large number of reportedly pathogenic mtDNA variants, 41, as polymorphisms. These results drastically update our view on the implication of the primary sequence of mt-rRNA in the etiology of deafness and mitochondrial disease in general. Our data sheds much-needed light on the question of how mt-rRNA variants located at non-conserved positions may lead to mitochondrial disease and, most notably, provide evidence of the effect of haplotype context in the manifestation of some mt-rRNA variants.

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.

Uncovering translation roadblocks during the development of a synthetic tRNA

Ribosomes are remarkable in their malleability to accept diverse aminoacyl-tRNA substrates from both the same organism and other organisms or domains of life. This is a critical feature of the ribosome that allows the use of orthogonal translation systems for genetic code expansion. Optimization of these orthogonal translation systems generally involves focusing on the compatibility of the tRNA, aminoacyl-tRNA synthetase, and a non-canonical amino acid with each other. As we expand the diversity of tRNAs used to include non-canonical structures, the question arises as to the tRNA suitability on the ribosome. Specifically, we investigated the ribosomal translation of allo-tRNAUTu1, a uniquely shaped (9/3) tRNA exploited for site-specific selenocysteine insertion, using single-molecule fluorescence. With this technique we identified ribosomal disassembly occurring from translocation of allo-tRNAUTu1 from the A to the P site. Using cryo-EM to capture the tRNA on the ribosome, we pinpointed a distinct tertiary interaction preventing fluid translocation. Through a single nucleotide mutation, we disrupted this tertiary interaction and relieved the translation roadblock. With the continued diversification of genetic code expansion, our work highlights a targeted approach to optimize translation by distinct tRNAs as they move through the ribosome.

Structural Insights Into the 5′UG/3′GU Wobble Tandem in Complex With Ba2+ Cation

G•U wobble base pair frequently occurs in RNA structures. The unique chemical, thermodynamic, and structural properties of the G•U pair are widely exploited in RNA biology. In several RNA molecules, the G•U pair plays key roles in folding, ribozyme catalysis, and interactions with proteins. G•U may occur as a single pair or in tandem motifs with different geometries, electrostatics, and thermodynamics, further extending its biological functions. The metal binding affinity, which is essential for RNA folding, catalysis, and other interactions, differs with respect to the tandem motif type due to the different electrostatic potentials of the major grooves. In this work, we present the crystal structure of an RNA 8-mer duplex r[UCGUGCGA]₂, providing detailed structural insights into the tandem motif I (5′UG/3′GU) complexed with Ba²⁺ cation. We compare the electrostatic potential of the presented motif I major groove with previously published structures of tandem motifs I, II (5′GU/3′UG), and III (5′GG/3′UU). A local patch of a strongly negative electrostatic potential in the major groove of the presented structure forms the metal binding site with the contributions of three oxygen atoms from the tandem. These results give us a better understanding of the G•U tandem motif I as a divalent metal binder, a feature essential for RNA functions.

Conformational distortions induced by periodically recurring A…A in d(CAG).d(CAG) provide stereochemical rationale for the trapping of MSH2.MSH3 in polyQ disorders

CAG repeat instability causes a number of neurodegenerative disorders. The unusual hairpin stem structure formed by the CAG repeats in DNA traps the human mismatch repair MSH2.MSH3 (Mutsβ) complex. To understand the mechanism behind the abnormal binding of Mutsβ with the imperfect hairpin stem structure formed by CAG repeats, molecular dynamics simulations have been carried out for Mutsβ-d(CAG)₂(CAG)(CAG)₂.d(CTG)₂(CAG)(CTG)₂ (1 A…A mismatch) and Mutsβ-d(CAG)₅.d(CAG)₅ (5 mismatches, wherein, A…A occurs periodically) complexes. The interaction of MSH3 residue Tyr₂₄₅ at the minor groove side of A…A, an essential interaction responsible for the recognition by Mutsβ, are retained in both the cases. Nevertheless, the periodic unwinding caused by the nonisostericity of A…A with the flanking canonical base pairs in d(CAG)₅.d(CAG)₅ distorts the regular B-form geometry. Such an unwinding exposes one of the A…A mismatches (that interacts with Tyr₂₄₅) at the major groove side and also facilitates the on and off hydrogen bonding interaction with Lys₅₄₆ sidechain (MSH2-domain-IV). In contrast, kinking of the DNA towards the major groove in Mutsβ-d(CAG)₂(CAG)(CAG)₂.d(CTG)₂(CAG)(CTG)₂ doesn’t facilitate such an exposure of the bases at the major groove. Further, the unwinding of the helix in d(CAG)₅.d(CAG)₅ enhances the tighter binding between MSH2-domain-I and d(CAG)₅.d(CAG)₅ at the major groove side as well as between MSH3-domain-I and MSH3-domain-IV. Markedly, such enhanced interactions are absent in Mutsβ-d(CAG)₂(CAG)(CAG)₂.d(CTG)₂(CAG)(CTG)₂ that has a single A…A mismatch. Thus, the above-mentioned enhancement in intra- and inter- molecular interactions in Mutsβ-d(CAG)₅.d(CAG)₅ provide the stereochemical rationale for the trapping of Mutsβ in CAG repeat expansion disorders.

Characterization of the complete mitochondrial genome of the many-lined sun skink (Eutropis multifasciata) and comparison with other Scincomorpha species

Characterizating the complete mitochondrial genome (mitogenome) of an organism allows detailed genomic studies in systematics and evolution. The present study decodes the mitogenome (17,062 bp) of the many-lined sun skink, Eutropis multifasciata, using next-generation sequencing. To compare the diversity of mitogenomic structure and investigate intraspecific evolutionary relationships among the Asian Scincomorpha, the mitogenomes of 46 other species were examined concurrently. Within the group, the size of mitogenomes varied predominantly in the length at their control regions. The Ka/Ks ratios of 12 protein codon genes (PCGs) were lower than 1.00, demonstrating that they were under relaxed or moderate purifying selection. However, the ND5 had a Ka/Ks ratio >1, and was considered to be under positive selection. Currently there are two superfamilies in Scincomorpha (i.e. Scincoidea and Lacertoidea), but phylogenetic analysis using Bayesian Inference and Maximum-Likelihood Estimations produced phylogenetic trees with three clades in Scincomorpha ((Scincoidea + Lacertoidea (part)) + Gymnophthalmidae)).

Secondary structural choice of DNA and RNA associated with CGG/CCG trinucleotide repeat expansion rationalizes the RNA misprocessing in FXTAS

CGG tandem repeat expansion in the 5'-untranslated region of the fragile X mental retardation-1 (FMR1) gene leads to unusual nucleic acid conformations, hence causing genetic instabilities. We show that the number of G…G (in CGG repeat) or C…C (in CCG repeat) mismatches (other than A…T, T…A, C…G and G…C canonical base pairs) dictates the secondary structural choice of the sense and antisense strands of the FMR1 gene and their corresponding transcripts in fragile X-associated tremor/ataxia syndrome (FXTAS). The circular dichroism (CD) spectra and electrophoretic mobility shift assay (EMSA) reveal that CGG DNA (sense strand of the FMR1 gene) and its transcript favor a quadruplex structure. CD, EMSA and molecular dynamics (MD) simulations also show that more than four C…C mismatches cannot be accommodated in the RNA duplex consisting of the CCG repeat (antisense transcript); instead, it favors an i-motif conformational intermediate. Such a preference for unusual secondary structures provides a convincing justification for the RNA foci formation due to the sequestration of RNA-binding proteins to the bidirectional transcripts and the repeat-associated non-AUG translation that are observed in FXTAS. The results presented here also suggest that small molecule modulators that can destabilize FMR1 CGG DNA and RNA quadruplex structures could be promising candidates for treating FXTAS.

Sequence dependent influence of an A…A mismatch in a DNA duplex: An insight into the recognition by hZαADAR1 protein

Base pair mismatches can erroneously be incorporated in the DNA. An adenine pairing with another adenine is one of the eight possible mismatches. The atomistic insights about the structure and dynamics of an A…A mismatch in a DNA (unbound form) is not yet accessible to any experimental technique. Earlier molecular dynamics (MD) simulations have shown that A…A mismatch in the midst of 5'CAG/3'GAC, 5'GAC/3'CAG and 5'CAA/3'GAT (underline represents the mismatch) are highly dynamic in nature. By employing MD simulation, the influence of an A…A mismatch in the midst of 5'GAA/3'CAT, 5'GAG/3'CAC, 5'AAC/3'TAG, 5'AAG/3'TAC, 5'TAA/3'AAT, 5'TAT/3'AAA and 5'AAT/3'TAA sequences have been investigated here. The results indicate that irrespective of the flanking sequences, the mismatch samples a variety of transient conformations, including a B-Z junction. Further, circular dichroism studies have been carried out to explore the ability of these sequences to bind with hZα_ADAR1 which specifically recognizes B-Z junction/Z-DNA. The results indicate that hZα_ADAR1 could not lead to a complete B to Z transition in the above sequences. Notably, a complete transition to Z-form has been reported earlier for 5'GAC/3'CAG upon titrating with hZα_ADAR1. Intriguingly, 5'AAC/3'TAG, 5'AAG/3'TAC and 5'GAA/3'CAT exhibit a B-Z junction formation rather than a complete transition to Z-form, similar to the situation of 5'CAA/3'GAT. These indicate that although A…A mismatch could induce a local B-Z junction transiently, hZα_ADAR1 requires the presence of a G…C/C…G base pair adjacent to the A…A mismatch for the binding. Additionally, the extent of B-Z junction has enhanced upon binding with hZα_ADAR1 in the presence of the A…A mismatch (specifically when CG, CA, AC, GA and AG steps occur), but not in the presence of the canonical base pairs. These confirm the inclination of A…A mismatch towards the B-Z junction.

Solution NMR readily reveals distinct structural folds and interactions in doubly 13C- and 19F-labeled RNAs

RNAs form critical components of biological processes implicated in human diseases, making them attractive for small-molecule therapeutics. Expanding the sites accessible to nuclear magnetic resonance (NMR) spectroscopy will provide atomic-level insights into RNA interactions. Here, we present an efficient strategy to introduce ¹⁹F-¹³C spin pairs into RNA by using a 5-fluorouridine-5'-triphosphate and T7 RNA polymerase-based in vitro transcription. Incorporating the ¹⁹F-¹³C label in two model RNAs produces linewidths that are twice as sharp as the commonly used ¹H-¹³C spin pair. Furthermore, the high sensitivity of the ¹⁹F nucleus allows for clear delineation of helical and nonhelical regions as well as GU wobble and Watson-Crick base pairs. Last, the ¹⁹F-¹³C label enables rapid identification of a small-molecule binding pocket within human hepatitis B virus encapsidation signal epsilon (hHBV ε) RNA. We anticipate that the methods described herein will expand the size limitations of RNA NMR and aid with RNA-drug discovery efforts.

Stacking geometry between two sheared Watson-Crick basepairs: Computational chemistry and bioinformatics based prediction

BACKGROUND

Molecular modeling of RNA double helices is possible using most probable values of basepair parameters obtained from crystal structure database. The A:A w:wC non-canonical basepair, involving Watson-Crick edges of two Adenines in cis orientation, appears quite frequently in database. Bimodal distribution of its Shear, due to two different H-bonding schemes, introduces the confusion in assigning most the probable value. Its effect is pronounced when the A:A w:wC basepair stacks on Sheared wobble G:U W:WC basepairs.

METHODS

We employed molecular dynamics simulations of three possible double helices with GAG, UAG and GAU sequence motifs at their centers and quantum chemical calculation for non-canonical A:A w:wC basepair stacked on G:U W:WC basepair.

RESULTS

We noticed stable structures of GAG motif with specifically negative Shear of the A:A basepair but stabilities of the other motifs were not found with A:A w:wC basepairing. Hybrid DFT-D and MP2 stacking energy analyses on dinucleotide step sequences, A:A w:wC::G:U W:WC and A:A w:wC::U:G W:WC reveal that viable orientation of A:A::G:U prefers one of the H-bonding modes with negative Shear, supported by crystal structure database. The A:A::U:G dinucleotide, however, prefers structure with only positive Shear.

CONCLUSIONS

The quantum chemical calculations explain why MD simulations of GAG sequence motif only appear stable. In the cases of the GAU and UAG motifs "tug of war" situation between positive and negative Shears of A:A w:wC basepair induces conformational plasticity.

GENERAL SIGNIFICANCE

We have projected comprehensive reason behind the promiscuous nature of A:A w:wC basepair which brings occasional structural plasticity.

Adenine·cytosine substitutions are an alternative pathway of compensatory mutation in angiosperm ITS2

Compensatory mutations are crucial for functional RNA because they maintain RNA configuration and thus function. Compensatory mutation has traditionally been considered to be a two-step substitution through the GU-base-pair intermediate. We tested for an alternative AC-mediated compensatory mutation (ACCM). We investigated ACCMs by using a comprehensive sampling of ribosomal internal transcribed spacer 2 (ITS2) from 3934 angiosperm species in 80 genera and 55 families. We predicted ITS2 consensus secondary structures by using LocARNA for structure-based alignment and partitioning paired and unpaired regions. We examined and compared the substitution rates and frequencies among base pairs by using RNA-specific models. Base-pair states of ACCMs were mapped onto the inferred phylogenetic trees to infer their evolution. All types of compensatory mutations involving the AC intermediate were observed, but the most frequent substitutions were with AU or GC pairs, which are part of the AU-AC-GC pathway. Compared with the GU intermediate, AC had a lower frequency and higher mutability. Within the AU-AC-GC pathway, the AU-AC substitution rate was much slower than the AC-GC substitution rate. No consistently higher overall rate was identified for either pathway among all 80 sampled lineages, though compensatory mutations through the AC intermediate averaged about half that through the GU intermediate. These results demonstrate an alternative compensatory mutation between AU and GC that helps address the controversial inference of inferred simultaneous double substitutions.

The multiple flavors of GoU pairs in RNA

Wobble GU pairs (or GoU) occur frequently within double‐stranded RNA helices interspersed within the standard G═C and A─U Watson‐Crick pairs. However, other types of GoU pairs interacting on their Watson‐Crick edges have been observed. The structural and functional roles of such alternative GoU pairs are surprisingly diverse and reflect the various pairings G and U can form by exploiting all the subtleties of their electronic configurations. Here, the structural characteristics of the GoU pairs are updated following the recent crystallographic structures of functional ribosomal complexes and the development in our understanding of ribosomal translation.

Modulation of Bacterial sRNAs Activity by Epigenetic Modifications: Inputs from the Eukaryotic miRNAs

Trans-encoded bacterial regulatory RNAs (sRNAs) are functional analogues of eukaryotic microRNAs (miRNAs). These RNA classes act by base-pairing complementarity with their RNA targets to modulate gene expression (transcription, half-life and/or translation). Based on base-pairing, algorithms predict binding and the impact of small RNAs on targeted-RNAs expression and fate. However, other actors are involved such as RNA binding proteins and epigenetic modifications of the targeted and small RNAs. Post-transcriptional base modifications are widespread in all living organisms where they lower undesired RNA folds through conformation adjustments and influence RNA pairing and stability, especially if remodeling their ends. In bacteria, sRNAs possess RNA modifications either internally (methylation, pseudouridinylation) or at their ends. Nicotinamide adenine dinucleotide were detected at 5′-ends, and polyadenylation can occur at 3′-ends. Eukaryotic miRNAs possess N⁶-methyladenosine (m⁶A), A editing into I, and non-templated addition of uridines at their 3′-ends. Biological functions and enzymes involved in those sRNA and micro RNA epigenetic modifications, when known, are presented and challenged.

A B-Z junction induced by an A … A mismatch in GAC repeats in the gene for cartilage oligomeric matrix protein promotes binding with the hZαADAR1 protein

GAC repeat expansion from five to seven in the exonic region of the gene for cartilage oligomeric matrix protein (COMP) leads to pseudoachondroplasia, a skeletal abnormality. However, the molecular mechanism by which GAC expansions in the COMP gene lead to skeletal dysplasias is poorly understood. Here we used molecular dynamics simulations, which indicate that an A … A mismatch in a d(GAC)6·d(GAC)6 duplex induces negative supercoiling, leading to a local B-to-Z DNA transition. This transition facilitates the binding of d(GAC)7·d(GAC)7 with the Zα-binding domain of human adenosine deaminase acting on RNA 1 (ADAR1, hZαADAR1), as confirmed by CD, NMR, and microscale thermophoresis studies. The CD results indicated that hZαADAR1 recognizes the zigzag backbone of d(GAC)7·d(GAC)7 at the B-Z junction and subsequently converts it into Z-DNA via the so-called passive mechanism. Molecular dynamics simulations carried out for the modeled hZαADAR1-d(GAC)6d(GAC)6 complex confirmed the retention of previously reported important interactions between the two molecules. These findings suggest that hZαADAR1 binding with the GAC hairpin stem in COMP can lead to a non-genetic, RNA editing-mediated substitution in COMP that may then play a crucial role in the development of pseudoachondroplasia.

Structure of RNA Stem Loop B from the Picornavirus Replication Platform

The presumptive RNA cloverleaf at the start of the 5'-untranslated region of the picornavirus genome is an essential element in replication. Stem loop B (SLB) of the cloverleaf is a recognition site for the host polyC-binding protein, which initiates a switch from translation to replication. Here we present the solution structure of human rhinovirus isotype 14 SLB using nuclear magnetic resonance spectroscopy. SLB adopts a predominantly A-form helical structure. The stem contains five Watson-Crick base pairs and one wobble base pair and is capped by an eight-nucleotide loop. The wobble base pair introduces perturbations into the helical parameters but does not appear to introduce flexibility. However, the helix major groove appears to be accessible. Flexibility is seen throughout the loop and in the terminal nucleotides. The pyrimidine-rich region of the loop, the apparent recognition site for the polyC-binding protein, is the most disordered region of the structure.

Copper(II)-Thymine Coordination Polymer Nanoribbons as Potential Oligonucleotide Nanocarriers

The direct reaction between copper nitrate, thymine-1-acetic acid, and 4,4'-bipyridine in water leads to the formation of a blue colloid comprising uniform crystalline nanoribbons (length >1 μm; width ca. 150-185 nm; diameter ca. 15-60 nm) of a coordination polymer. The polymer displays a thymine-based structure freely available for supramolecular interactions. These nanostructures show significant selective interaction with single-stranded oligonucleotides based on adenine. Remarkably, they present low cell toxicity in three cell lines-despite the copper(II) content-and can be used as nanocarriers of oligonucleotides. These results suggest the potential of these types of nanostructures in several biological applications.

Common miR-590 Variant rs6971711 Present Only in African Americans Reduces miR-590 Biogenesis

MicroRNAs (miRNAs) are recognized as important regulators of cardiac development, hypertrophy and fibrosis. Recent studies have demonstrated that genetic variations which cause alterations in miRNA:target interactions can lead to disease. We hypothesized that genetic variations in miRNAs that regulate cardiac hypertrophy/fibrosis might be involved in generation of the cardiac phenotype in patients diagnosed with hypertrophic cardiomyopathy (HCM). To investigate this question, we Sanger sequenced 18 miRNA genes previously implicated in myocyte hypertrophy/fibrosis and apoptosis, using genomic DNA isolated from the leukocytes of 199 HCM patients. We identified a single nucleotide polymorphism (rs6971711, C57T SNP) at the 17th position of mature miR-590-3p (= 57th position of pre-miR-590) that is common in individuals of African ancestry. SNP frequency was higher in African American HCM patients (n = 55) than ethnically-matched controls (n = 100), but the difference was not statistically significant (8.2% vs. 6.5%; p = 0.5). Using a cell culture system, we discovered that presence of this SNP resulted in markedly lower levels of mature miR-590-5p (39 ± 16%, p<0.003) and miR-590-3p (20 ± 2%, p<0.003), when compared with wild-type (WT) miR-590, without affecting levels of pri-miR-590 and pre-miR-590. Consistent with this finding, the SNP resulted in reduced target suppression when compared to WT miR-590 (71% suppression by WT vs 60% suppression by SNP, p<0.03). Since miR-590 can regulate TGF-β, Activin A and Akt signaling, SNP-induced reduction in miR-590 biogenesis could influence cardiac phenotype by de-repression of these signaling pathways. Since the SNP is only present in African Americans, population studies in this patient population would be valuable to investigate effects of this SNP on myocyte function and cardiac physiology.

Selective Preference of Parallel DNA Triplexes Is Due to the Disruption of Hoogsteen Hydrogen Bonds Caused by the Severe Nonisostericity between the G*GC and T*AT Triplets

Implications of DNA, RNA and RNA.DNA hybrid triplexes in diverse biological functions, diseases and therapeutic applications call for a thorough understanding of their structure-function relationships. Despite exhaustive studies mechanistic rationale for the discriminatory preference of parallel DNA triplexes with G*GC & T*AT triplets still remains elusive. Here, we show that the highest nonisostericity between the G*GC & T*AT triplets imposes extensive stereochemical rearrangements contributing to context dependent triplex destabilisation through selective disruption of Hoogsteen scheme of hydrogen bonds. MD simulations of nineteen DNA triplexes with an assortment of sequence milieu reveal for the first time fresh insights into the nature and extent of destabilization from a single (non-overlapping), double (overlapping) and multiple pairs of nonisosteric base triplets (NIBTs). It is found that a solitary pair of NIBTs, feasible either at a G*GC/T*AT or T*AT/G*GC triplex junction, does not impinge significantly on triplex stability. But two overlapping pairs of NIBTs resulting from either a T*AT or a G*GC interruption disrupt Hoogsteen pair to a noncanonical mismatch destabilizing the triplex by ~10 to 14 kcal/mol, implying that their frequent incidence in multiples, especially, in short sequences could even hinder triplex formation. The results provide (i) an unambiguous and generalised mechanistic rationale for the discriminatory trait of parallel triplexes, including those studied experimentally (ii) clarity for the prevalence of antiparallel triplexes and (iii) comprehensive perspectives on the sequence dependent influence of nonisosteric base triplets useful in the rational design of TFO's against potential triplex target sites.

Coarse-Grained Modeling of Nucleic Acids Using Anisotropic Gay-Berne and Electric Multipole Potentials

In this work, we attempt to apply a coarse-grained (CG) model, which is based on anisotropic Gay-Berne and electric multipole (EMP) potentials, to the modeling of nucleic acids. First, a comparison has been made between the CG and atomistic models (AMBER point-charge model) in the modeling of DNA and RNA hairpin structures. The CG results have demonstrated a good quality in maintaining the nucleic acid hairpin structures, in reproducing the dynamics of backbone atoms of nucleic acids, and in describing the hydrogen-bonding interactions between nucleic acid base pairs. Second, the CG and atomistic AMBER models yield comparable results in modeling double-stranded DNA and RNA molecules. It is encouraging that our CG model is capable of reproducing many elastic features of nucleic acid base pairs in terms of the distributions of the interbase pair step parameters (such as shift, slide, tilt, and twist) and the intrabase pair parameters (such as buckle, propeller, shear, and stretch). Finally, The GBEMP model has shown a promising ability to predict the melting temperatures of DNA duplexes with different lengths.

Twisting right to left: A…A mismatch in a CAG trinucleotide repeat overexpansion provokes left-handed Z-DNA conformation

Conformational polymorphism of DNA is a major causative factor behind several incurable trinucleotide repeat expansion disorders that arise from overexpansion of trinucleotide repeats located in coding/non-coding regions of specific genes. Hairpin DNA structures that are formed due to overexpansion of CAG repeat lead to Huntington’s disorder and spinocerebellar ataxias. Nonetheless, DNA hairpin stem structure that generally embraces B-form with canonical base pairs is poorly understood in the context of periodic noncanonical A…A mismatch as found in CAG repeat overexpansion. Molecular dynamics simulations on DNA hairpin stems containing A…A mismatches in a CAG repeat overexpansion show that A…A dictates local Z-form irrespective of starting glycosyl conformation, in sharp contrast to canonical DNA duplex. Transition from B-to-Z is due to the mechanistic effect that originates from its pronounced nonisostericity with flanking canonical base pairs facilitated by base extrusion, backbone and/or base flipping. Based on these structural insights we envisage that such an unusual DNA structure of the CAG hairpin stem may have a role in disease pathogenesis. As this is the first study that delineates the influence of a single A…A mismatch in reversing DNA helicity, it would further have an impact on understanding DNA mismatch repair.

When a set of 3 nucleotides in a DNA sequence repeats beyond a certain number, it leads to incurable neurological or neuromuscular disorders. Such DNA sequences tend to form unusual DNA structures comprising of base pairing schemes different from the canonical A…T/G…C base pairs. Influence of such unusual base pairing on the overall 3-dimensional structure of DNA and its impact on the pathogenesis of disorder is not well understood. CAG repeat overexpansion that leads to Huntington’s disorder and several spinocerebellar ataxias forms noncanonical A…A base pair in between canonical C…G and G…C base pairs. However, no detailed structural information is available on the influence of an A…A mismatch on a DNA structure under any sequence context. Here, we have shown for the first time that A…A base pairing in a CAG repeat provokes the formation of left-handed Z-DNA due to the pronounced structural dissimilarity of A…A base pair with G…C base pair, leading to periodic B-Z junction. Thus, these results suggest that formation of periodic B-Z junction may be one of the molecular bases for CAG repeat instability.

Targeting vertebrate intron-encoded box C/D 2'-O-methylation guide RNAs into the Cajal body

Post-transcriptional pseudouridylation and 2'-O-methylation of splicesomal small nuclear ribonucleic acids (snRNAs) is mediated by box H/ACA and box C/D small Cajal body (CB)-specific ribonucleoproteins (scaRNPs), respectively. The WD-repeat protein 79 (WDR79) has been proposed to interact with both classes of modification scaRNPs and target them into the CB. The box H/ACA scaRNAs carry the common CAB box motif (consensus, ugAG) that is required for both WDR79 binding and CB-specific accumulation. Thus far, no cis-acting CB-localization element has been reported for vertebrate box C/D scaRNAs. In this study, systematic mutational analysis of the human U90 and another newly identified box C/D scaRNA, mgU2-47, demonstrated that the CB-specific accumulation of vertebrate intron-encoded box C/D scaRNAs relies on GU- or UG-dominated dinucleotide repeat sequences which are predicted to form the terminal stem-loop of the RNA apical hairpin. While the loop nucleotides are unimportant, the adjacent terminal helix that is composed mostly of consecutive G.U and U.G wobble base-pairs is essential for CB-specific localization of box C/D scaRNAs. Co-immunoprecipitation experiments confirmed that the newly identified CB localization element, called the G.U/U.G wobble stem, is crucial for in vivo association of box C/D scaRNPs with WDR79.

Is the DPT tautomerization of the long A·G Watson-Crick DNA base mispair a source of the adenine and guanine mutagenic tautomers? A QM and QTAIM response to the biologically important question

Herein, we first address the question posed in the title by establishing the tautomerization trajectory via the double proton transfer of the adenine·guanine (A·G) DNA base mispair formed by the canonical tautomers of the A and G bases into the A*·G* DNA base mispair, involving mutagenic tautomers, with the use of the quantum-mechanical calculations and quantum theory of atoms in molecules (QTAIM). It was detected that the A·G ↔ A*·G* tautomerization proceeds through the asynchronous concerted mechanism. It was revealed that the A·G base mispair is stabilized by the N6H···O6 (5.68) and N1H···N1 (6.51) hydrogen bonds (H-bonds) and the N2H···HC2 dihydrogen bond (DH-bond) (0.68 kcal·mol(-1) ), whereas the A*·G* base mispair-by the O6H···N6 (10.88), N1H···N1 (7.01) and C2H···N2 H-bonds (0.42 kcal·mol(-1) ). The N2H···HC2 DH-bond smoothly and without bifurcation transforms into the C2H···N2 H-bond at the IRC = -10.07 Bohr in the course of the A·G ↔ A*·G* tautomerization. Using the sweeps of the energies of the intermolecular H-bonds, it was observed that the N6H···O6 H-bond is anticooperative to the two others-N1H···N1 and N2H···HC2 in the A·G base mispair, while the latters are significantly cooperative, mutually strengthening each other. In opposite, all three O6H···N6, N1H···N1, and C2H···N2 H-bonds are cooperative in the A*·G* base mispair. All in all, we established the dynamical instability of the А*·G* base mispair with a short lifetime (4.83·10(-14) s), enabling it not to be deemed feasible source of the A* and G* mutagenic tautomers of the DNA bases. The small lifetime of the А*·G* base mispair is predetermined by the negative value of the Gibbs free energy for the A*·G* → A·G transition. Moreover, all of the six low-frequency intermolecular vibrations cannot develop during this lifetime that additionally confirms the aforementioned results. Thus, the A*·G* base mispair cannot be considered as a source of the mutagenic tautomers of the DNA bases, as the A·G base mispair dissociates during DNA replication exceptionally into the A and G monomers in the canonical tautomeric form.