Crystal structure of an RNA.DNA hybrid reveals intermolecular intercalation: dimer formation by base-pair swapping

RPA transforms RNase H1 to a bidirectional exoribonuclease for processive RNA-DNA hybrid cleavage

RNase H1 has been acknowledged as an endoribonuclease specializing in the internal degradation of the RNA moiety within RNA-DNA hybrids, and its ribonuclease activity is indispensable in multifaceted aspects of nucleic acid metabolism. However, the molecular mechanism underlying RNase H1-mediated hybrid cleavage remains inadequately elucidated. Herein, using single-molecule approaches, we probe the dynamics of the hybrid cleavage by Saccharomyces cerevisiae RNase H1. Remarkably, a single RNase H1 enzyme displays 3'-to-5' exoribonuclease activity. The directional RNA degradation proceeds processively and yet discretely, wherein unwinding approximately 6-bp hybrids as a prerequisite for two consecutive 3-nt RNA excisions limits the overall rate within each catalytic cycle. Moreover, Replication Protein A (RPA) reinforces RNase H1's 3'-to-5' nucleolytic rate and processivity and stimulates its 5'-to-3' exoribonuclease activity. This stimulation is primarily realized through the pre-separation of the hybrids and consequently transfers RNase H1 to a bidirectional exoribonuclease, further potentiating its cleavage efficiency. These findings unveil unprecedented characteristics of an RNase and provide a dynamic view of RPA-enhanced processive hybrid cleavage by RNase H1.

Comprehensive Assessment of Force-Field Performance in Molecular Dynamics Simulations of DNA/RNA Hybrid Duplexes

Mixed double helices formed by RNA and DNA strands, commonly referred to as hybrid duplexes or hybrids, are essential in biological processes like transcription and reverse transcription. They are also important for their applications in CRISPR gene editing and nanotechnology. Yet, despite their significance, the hybrid duplexes have been seldom modeled by atomistic molecular dynamics methodology, and there is no benchmark study systematically assessing the force-field performance. Here, we present an extensive benchmark study of polypurine tract (PPT) and Dickerson-Drew dodecamer hybrid duplexes using contemporary and commonly utilized pairwise additive and polarizable nucleic acid force fields. Our findings indicate that none of the available force-field choices accurately reproduces all the characteristic structural details of the hybrid duplexes. The AMBER force fields are unable to populate the C3'-endo (north) pucker of the DNA strand and underestimate inclination. The CHARMM force field accurately describes the C3'-endo pucker and inclination but shows base pair instability. The polarizable force fields struggle with accurately reproducing the helical parameters. Some force-field combinations even demonstrate a discernible conflict between the RNA and DNA parameters. In this work, we offer a candid assessment of the force-field performance for mixed DNA/RNA duplexes. We provide guidance on selecting utilizable force-field combinations and also highlight potential pitfalls and best practices for obtaining optimal performance.

Structure and Dynamics of DNA and RNA Double Helices Formed by d(CTG), d(GTC), r(CUG), and r(GUC) Trinucleotide Repeats and Associated DNA-RNA Hybrids

Myotonic dystrophy type 1 is the most frequent form of muscular dystrophy in adults caused by an abnormal expansion of the CTG trinucleotide. Both the expanded DNA and the expanded CUG RNA transcript can fold into hairpins. Co-transcriptional formation of stable RNA·DNA hybrids can also enhance the instability of repeat tracts. We performed molecular dynamics simulations of homoduplexes associated with the disease, d(CTG)n and r(CUG)n, and their corresponding r(CAG)n:d(CTG)n and r(CUG)n:d(CAG)n hybrids that can form under bidirectional transcription and of non-pathological d(GTC)n and d(GUC)n homoduplexes. We characterized their conformations, stability, and dynamics and found that the U·U and T·T mismatches are dynamic, favoring anti-anti conformations inside the helical core, followed by anti-syn and syn-syn conformations. For DNA, the secondary minima in the non-expanding d(GTC)n helices are deeper, wider, and longer-lived than those in d(CTG)n, which constitutes another biophysical factor further differentiating the expanding and non-expanding sequences. The hybrid helices are closer to A-RNA, with the A-T and A-U pairs forming two stable Watson-Crick hydrogen bonds. The neutralizing ion distribution around the non-canonical pairs is also described.

Targeting the ALS/FTD-associated A-DNA kink with anthracene-based metal complex causes DNA backbone straightening and groove contraction

The use of a small molecule compound to reduce toxic repeat RNA transcripts or their translated aberrant proteins to target repeat-expanded RNA/DNA with a G4C2 motif is a promising strategy to treat C9orf72-linked disorders. In this study, the crystal structures of DNA and RNA–DNA hybrid duplexes with the -GGGCCG- region as a G4C2 repeat motif were solved. Unusual groove widening and sharper bending of the G4C2 DNA duplex A-DNA conformation with B-form characteristics inside was observed. The G4C2 RNA–DNA hybrid duplex adopts a more typical rigid A form structure. Detailed structural analysis revealed that the G4C2 repeat motif of the DNA duplex exhibits a hydration shell and greater flexibility and serves as a ‘hot-spot’ for binding of the anthracene-based nickel complex, Ni^II(Chro)₂ (Chro = Chromomycin A3). In addition to the original GGCC recognition site, Ni^II(Chro)₂ has extended specificity and binds the flanked G:C base pairs of the GGCC core, resulting in minor groove contraction and straightening of the DNA backbone. We have also shown that Chro-metal complexes inhibit neuronal toxicity and suppresses locomotor deficits in a Drosophila model of C9orf72-associated ALS. The approach represents a new direction for drug discovery against ALS and FTD diseases by targeting G4C2 repeat motif DNA.

Design, Optimization, and Study of Small Molecules That Target Tau Pre-mRNA and Affect Splicing

Approximately 95% of human genes are alternatively spliced, and aberrant splicing events can cause disease. One pre-mRNA that is alternatively spliced and linked to neurodegenerative diseases is tau (microtubule-associated protein tau), which can cause frontotemporal dementia and parkinsonism linked to chromosome 17 (FTDP-17) and can contribute to Alzheimer's disease. Here, we describe the design of structure-specific lead small molecules that directly target tau pre-mRNA from sequence. This was followed by hit expansion and analogue synthesis to further improve upon these initial lead molecules. The emergent compounds were assessed for functional activity in a battery of assays, including binding assays and an assay that mimics molecular recognition of tau pre-mRNA by a U1 small nuclear ribonucleoprotein (snRNP) splicing factor. Compounds that emerged from these studies had enhanced potency and selectivity for the target RNA relative to the initial hits, while also having significantly improved drug-like properties. The compounds are shown to directly target tau pre-mRNA in cells, via chemical cross-linking and isolation by pull-down target profiling, and to rescue disease-relevant splicing of tau pre-mRNA in a variety of cellular systems, including primary neurons. More broadly, this study shows that lead, structure-specific compounds can be designed from sequence and then further optimized for their physicochemical properties while at the same time enhancing their activity.

Structural Flexibility of DNA-RNA Hybrid Duplex: Stretching and Twist-Stretch Coupling

DNA-RNA hybrid (DRH) duplexes play essential roles during the replication of DNA and the reverse transcription of RNA viruses, and their flexibility is important for their biological functions. Recent experiments indicated that A-form RNA and B-form DNA have a strikingly different flexibility in stretching and twist-stretch coupling, and the structural flexibility of DRH duplex is of great interest, especially in stretching and twist-stretch coupling. In this work, we performed microsecond all-atom molecular dynamics simulations with new AMBER force fields to characterize the structural flexibility of DRH duplex in stretching and twist-stretch coupling. We have calculated all the helical parameters, stretch modulus, and twist-stretch coupling parameters for the DRH duplex. First, our analyses on structure suggest that the DRH duplex exhibits an intermediate conformation between A- and B-forms and closer to A-form, which can be attributed to the stronger rigidity of the RNA strand than the DNA strand. Second, our calculations show that the DRH duplex has the stretch modulus of 834 ± 34 pN and a very weak twist-stretch coupling. Our quantitative analyses indicate that, compared with DNA and RNA duplexes, the different flexibility of the DRH duplex in stretching and twist-stretch coupling is mainly attributed to its apparently different basepair inclination in the helical structure.

Crystal structure of RNA-DNA duplex provides insight into conformational changes induced by RNase H binding

RNA-DNA hybrids play essential roles in a variety of biological processes, including DNA replication, transcription, and viral integration. Ribonucleotides incorporated within DNA are hydrolyzed by RNase H enzymes in a removal process that is necessary for maintaining genomic stability. In order to understand the structural determinants involved in recognition of a hybrid substrate by RNase H we have determined the crystal structure of a dodecameric non-polypurine/polypyrimidine tract RNA-DNA duplex. A comparison to the same sequence bound to RNase H, reveals structural changes to the duplex that include widening of the major groove to 12.5 Å from 4.2 Å and decreasing the degree of bending along the axis which may play a crucial role in the ribonucleotide recognition and cleavage mechanism within RNase H. This structure allows a direct comparison to be made about the conformational changes induced in RNA-DNA hybrids upon binding to RNase H and may provide insight into how dysfunction in the endonuclease causes disease.

DNA–RNA hybrid duplexes with decreasing pyrimidine content in the DNA strand provide structural snapshots for the A- to B-form conformational transition of nucleic acids

A gradual increase in the deoxypyrimidine content in DNA–RNA hybrids leads to B- to A-form nucleic acid transition. Possible factors that govern nuclease activity on hybrid duplexes are presented.

Identification and role of the homodimerization interface of the glycosylphosphatidylinositol-anchored membrane type 6 matrix metalloproteinase (MMP25)

The membrane type (MT) 6 matrix metalloproteinase (MMP) (MMP25) is a glycosylphosphatidylinositol-anchored matrix metalloproteinase (MMP) that is highly expressed in leukocytes and in some cancer tissues. We previously showed that natural MT6-MMP is expressed on the cell surface as a major reduction-sensitive form of M(r) 120, likely representing enzyme homodimers held by disulfide bridges. Among the membrane type-MMPs, the stem region of MT6-MMP contains three cysteine residues at positions 530, 532, and 534 which may contribute to dimerization. A systematic site-directed mutagenesis study of the Cys residues in the stem region shows that Cys(532) is involved in MT6-MMP dimerization by forming an intermolecular disulfide bond. The mutagenesis data also suggest that Cys(530) and Cys(534) form an intramolecular disulfide bond. The experimental observations on cysteines were also investigated by computational studies of the stem peptide, which validate these proposals. Dimerization is not essential for transport of MT6-MMP to the cell surface, partitioning into lipid rafts or cleavage of alpha-1-proteinase inhibitor. However, monomeric forms of MT6-MMP exhibited enhanced autolysis and metalloprotease-dependent degradation. Collectively, these studies establish the stem region of MT6-MMP as the dimerization interface, an event whose outcome imparts protease stability to the protein.

Molecular dynamics of a kappaB DNA element: base flipping via cross-strand intercalative stacking in a microsecond-scale simulation

The sequence-dependent structural variability and conformational dynamics of DNA play pivotal roles in many biological milieus, such as in the site-specific binding of transcription factors to target regulatory elements. To better understand DNA structure, function, and dynamics in general, and protein···DNA recognition in the ‘κB’ family of genetic regulatory elements in particular, we performed molecular dynamics simulations of a 20-bp DNA encompassing a cognate κB site recognized by the proto-oncogenic ‘c-Rel’ subfamily of NF-κB transcription factors. Simulations of the κB DNA in explicit water were extended to microsecond duration, providing a broad, atomically detailed glimpse into the structural and dynamical behavior of double helical DNA over many timescales. Of particular note, novel (and structurally plausible) conformations of DNA developed only at the long times sampled in this simulation—including a peculiar state arising at ≈0.7 μs and characterized by cross-strand intercalative stacking of nucleotides within a longitudinally sheared base pair, followed (at ≈1 μs) by spontaneous base flipping of a neighboring thymine within the A-rich duplex. Results and predictions from the microsecond-scale simulation include implications for a dynamical NF-κB recognition motif, and are amenable to testing and further exploration via specific experimental approaches that are suggested herein.

RNase H activity: structure, specificity, and function in reverse transcription

This review compares the well-studied RNase H activities of human immunodeficiency virus, type 1 (HIV-1) and Moloney murine leukemia virus (MoMLV) reverse transcriptases. The RNase H domains of HIV-1 and MoMLV are structurally very similar, with functions assigned to conserved subregions like the RNase H primer grip and the connection subdomain, as well as to distinct features like the C-helix and loop in MoMLV RNase H. Like cellular RNases H, catalysis by the retroviral enzymes appears to involve a two-metal ion mechanism. Unlike cellular RNases H, the retroviral RNases H display three different modes of cleavage: internal, DNA 3' end-directed, and RNA 5' end-directed. All three modes of cleavage appear to have roles in reverse transcription. Nucleotide sequence is an important determinant of cleavage specificity with both enzymes exhibiting a preference for specific nucleotides at discrete positions flanking an internal cleavage site as well as during tRNA primer removal and plus-strand primer generation. RNA 5' end-directed and DNA 3' end-directed cleavages show similar sequence preferences at the positions closest to a cleavage site. A model for how RNase H selects cleavage sites is presented that incorporates both sequence preferences and the concept of a defined window for allowable cleavage from a recessed end. Finally, the RNase H activity of HIV-1 is considered as a target for anti-virals as well as a participant in drug resistance.

Atomic detail investigation of the structure and dynamics of DNA.RNA hybrids: a molecular dynamics study

DNA.RNA hybrid duplexes are biologically important molecules and are shown to have potential therapeutic properties. To investigate the relationship between structures, energetics, solvation and RNase H activity of hybrid duplexes in comparison with pure DNA and RNA duplexes, a molecular dynamics study using the CHARMM27 force field was undertaken. The structural properties of all four nucleic acids considered are in very good agreement with the experimental data. The backbone dihedral angles and the puckering of the (deoxy)ribose indicate that the purine rich strands retain their A-/B-like properties but the pyrimidine rich DNA strand undergoes A-B conformational transitions. The minor groove widths of the hybrid structures are narrower than those in the RNA duplex, a requirement for RNase H binding. In addition, sampling of noncanonical phosphodiester backbone dihedrals by the DNA strands, differential solvation properties and helical properties, most notably rise, are suggested to contribute to hybrids being RNase H substrates. Differential RNase H activity toward hybrids containing purine versus pyrimidine rich RNA strands is suggested to be due to sampling of values of the phosphodiester backbone dihedrals in the DNA strands. Notably, the present results indicate that hybrids have decreased flexibility as compared to RNA, in contrast to previous reports.

Base-tetrad swapping results in dimerization of RNA quadruplexes: implications for formation of the i-motif RNA octaplex

Nucleic acids adopt different multistranded helical architectures to perform various biological functions. Here, we report a crystal structure of an RNA quadruplex containing "base-tetrad swapping" and bulged nucleotide at 2.1-Angstroms resolution. The base-tetrad swapping results in a dimer of quadruplexes with an intercalated octaplex fragment at the 5' end junction. The intercalated base tetrads provide the basic repeat unit for constructing a model of intercalated RNA octaplex. The model we obtained shows fundamentally different characteristics from duplex, triplex, and quadruplex. We also observed two different orientations of bulged uridine residues that are related to the interaction with surroundings. This structural evidence reflects the conformational flexibility of bulged nucleotides in RNA quadruplexes and implies the potential roles of bulged nucleotides as recognition and interaction sites in RNA-protein and RNA-RNA interactions.

Sequence, distance, and accessibility are determinants of 5'-end-directed cleavages by retroviral RNases H

The RNase H activity of reverse transcriptase is essential for retroviral replication. RNA 5'-end-directed cleavages represent a form of RNase H activity that is carried out on RNA/DNA hybrids that contain a recessed RNA 5'-end. Previously, the distance from the RNA 5'-end has been considered the primary determinant for the location of these cleavages. Employing model hybrid substrates and the HIV-1 and Moloney murine leukemia virus reverse transcriptases, we demonstrate that cleavage sites correlate with specific sequences and that the distance from the RNA 5'-end determines the extent of cleavage. An alignment of sequences flanking multiple RNA 5'-end-directed cleavage sites reveals that both enzymes strongly prefer A or U at the +1 position and C or G at the -2 position, and additionally for HIV-1, A is disfavored at the -4 position. For both enzymes, 5'-end-directed cleavages occurred when sites were positioned between the 13th and 20th nucleotides from the RNA 5'-end, a distance termed the cleavage window. In examining the importance of accessibility to the RNA 5'-end, it was found that the extent of 5'-end-directed cleavages observed in substrates containing a free recessed RNA 5'-end was most comparable to substrates with a gap of two or three bases between the upstream and downstream RNAs. Together these finding demonstrate that the selection of 5'-end-directed cleavage sites by retroviral RNases H results from a combination of nucleotide sequence, permissible distance, and accessibility to the RNA 5'-end.

Nucleic acid crystallography: current progress

Fifty years after the publication of the DNA double helix model by Watson and Crick, new nucleic acid structures keep emerging at an ever-increasing rate. The past three years have brought a flurry of new oligonucleotide structures, including those of a Hoogsteen-paired DNA duplex, Holliday junctions, DNA-drug complexes, quadruplexes, a host of RNA motifs and various nucleic acid analogues. Major advances were also made in terms of the structure and function of catalytic RNAs. These range from improved models of the phosphodiester cleavage reactions catalyzed by the hairpin and hepatitis delta virus ribozymes to the visualization of a complete active site of a group I self-splicing intron with bound 5'- and 3'-exons. These triumphs are complemented by a refined understanding of cation-nucleic-acid interactions and new routes to the generation of derivatives for phasing of DNA and RNA structures.

Recognition of internal cleavage sites by retroviral RNases H

The RNase H activity of reverse transcriptase is essential to complete retroviral replication. Many studies have characterized how reverse transcriptase associates with recessed and exposed DNA 3' ends or RNA 5' ends to position the RNase H domain for cleavage, but little is known about how a nick might affect RNase H cleavages, or how RNase H carries out internal cleavages, which do not require positioning by a nucleic acid end. We have addressed these issues using model hybrid substrates and the reverse transcriptases of Moloney murine leukemia virus (M-MuLV) and human immunodeficiency virus type 1 (HIV-1). Our results show that a nick separating an upstream RNA and a downstream RNA annealed to DNA is essentially ignored by RNase H, indicating that the RNA 5' end at a nick is not sufficient to position 5' end-directed cleavages. Cleavage sites that are located close to the 5' end of the downstream RNA are not recognized in the absence of the upstream RNA, and the 5' ends of the shorter upstream RNAs enhance cleavage at these sites. The recognition of an internal cleavage site depends on local sequence features found both upstream and downstream of the cleavage site, designated as the -1/+1 position. By analyzing the nucleotide frequencies in the sequence surrounding strong internal cleavage sites, preferred nucleotides have been identified in the flanking sequences spanning positions -14 to +1 for HIV-1 and -11 to +1 for M-MuLV. These data reveal that general degradation of the retroviral genome after minus-strand synthesis can occur through sequence-specific cleavages.