Crystal structure of a DNA.RNA hybrid duplex with a polypurine RNA r(gaagaagag) and a complementary polypyrimidine DNA d(CTCTTCTTC)

Coarse-grained modeling of DNA-RNA hybrids

We introduce oxNA, a new model for the simulation of DNA-RNA hybrids that is based on two previously developed coarse-grained models-oxDNA and oxRNA. The model naturally reproduces the physical properties of hybrid duplexes, including their structure, persistence length, and force-extension characteristics. By parameterizing the DNA-RNA hydrogen bonding interaction, we fit the model's thermodynamic properties to experimental data using both average-sequence and sequence-dependent parameters. To demonstrate the model's applicability, we provide three examples of its use-calculating the free energy profiles of hybrid strand displacement reactions, studying the resolution of a short R-loop, and simulating RNA-scaffolded wireframe origami.

An in silico approach to evaluate the bindings of natural flavonoids and RNA-DNA hybrids

Flavonoids, low molecular weight polyphenolic compounds, are important natural products that belong to plant secondary metabolites. They have diverse biomedical applications such as antioxidative, anti-inflammatory, enzyme inhibitory, antimutagenic, anticarcinogenic, aromatase inhibitory effects, etc. Some of the flavonoids have been exported for bindings with certain DNA and tRNA structures both experimentally and computationally. RNA-DNA hybrid (RDH) falls into an important category of noncanonical nucleic acid structures that have many important biological functions. We have investigated the interaction of RDH structures with some of the dietary flavonoids with the aid of computational methods such as docking and molecular dynamics simulation. The presence of the - OH group on the ligand and the availability of a proper binding pocket in the macromolecule are the two main factors driving the binding preference. Thus, this computationally guided report explains the binding of the flavonoids with RDH structures to assist the researchers in designing noncanonical nucleic acid-targeted drug molecules.Communicated by Ramaswamy H. Sarma.

In-Cell Stability Prediction of RNA/DNA Hybrid Duplexes for Designing Oligonucleotides Aimed at Therapeutics

In cells, the formation of RNA/DNA hybrid duplexes regulates gene expression and modification. The environment inside cellular organelles is heterogeneously crowded with high concentrations of biomolecules that affect the structure and stability of RNA/DNA hybrid duplexes. However, the detailed environmental effects remain unclear. Therefore, the mechanistic details of the effect of such molecular crowding were investigated at the molecular level by using thermodynamic and nuclear magnetic resonance analyses, revealing structure-dependent destabilization of the duplexes under crowded conditions. The transition from B- to A-like hybrid duplexes due to a change in conformation of the DNA strand guided by purine-pyrimidine asymmetry significantly increased the hydration number, which resulted in greater destabilization by the addition of cosolutes. By quantifying the individual contributions of environmental factors and the bulk structure of the duplex, we developed a set of parameters that predict the stability of hybrid duplexes with conformational dissimilarities under diverse crowding conditions. A comparison of the effects of environmental conditions in living cells and in vitro crowded solutions on hybrid duplex formation using the Förster resonance energy transfer technique established the applicability of our parameters to living cells. Moreover, our derived parameters can be used to estimate the efficiency of transcriptional inhibition, genome editing, and silencing techniques in cells. This supports the usefulness of our parameters for the visualization of cellular mechanisms of gene expression and the development of nucleic acid-based therapeutics targeting different cells.

First High-Resolution Crystal Structures of DNA:2′-O-Methyl-RNA Heteroduplexes

Heteroduplexes composed of all-DNA and all-2′-OMe RNA strands do not occur in nature, but they have found application in the development of molecular beacons and could also be used as aptamers or elements of nucleic acid-based nanostructures that will contain such structural motifs. The crystallization experiments performed have shown that the introduction of overhangs at the ends of the duplex has a great influence on the success of crystallization, as well as on the DNA:2′-OMe-RNA heteroduplex crystal packing. The molecular and crystal structure of the DNA:2′-O-methyl-RNA heteroduplex in its overhanging and blunt-ended versions was determined at 100 K using synchrotron radiation with a resolution of 1.91 and 1.55 Å, respectively. The Zn-SAD method was used to resolve the original duplex structure when molecular replacement by many existing models of duplex structures failed. Both molecules analyzed adopted a conformation close to the A-RNA double helix. The presented structures provide the first insight into this type of heteroduplexes and allowed a comparative analysis with existing nucleic acid homo- and heteroduplex structures. The results of our research expand the knowledge of the structural properties of new heteroduplexes and may be useful for future applications, such as therapies using this class of compounds.

Alternative DNA Structures In Vivo: Molecular Evidence and Remaining Questions

Duplex DNA naturally folds into a right-handed double helix in physiological conditions. Some sequences of unusual base composition may nevertheless form alternative structures, as was shown for many repeated sequences in vitro However, evidence for the formation of noncanonical structures in living cells is difficult to gather. It mainly relies on genetic assays demonstrating their function in vivo or through genetic instability reflecting particular properties of such structures. Efforts were made to reveal their existence directly in a living cell, mainly by generating antibodies specific to secondary structures or using chemical ligands selected for their affinity to these structures. Among secondary structure-forming DNAs are G-quadruplexes, human fragile sites containing minisatellites, AT-rich regions, inverted repeats able to form cruciform structures, hairpin-forming CAG/CTG triplet repeats, and triple helices formed by homopurine-homopyrimidine GAA/TTC trinucleotide repeats. Many of these alternative structures are involved in human pathologies, such as neurological or developmental disorders, as in the case of trinucleotide repeats, or cancers triggered by translocations linked to fragile sites. This review will discuss and highlight evidence supporting the formation of alternative DNA structures in vivo and will emphasize the role of the mismatch repair machinery in binding mispaired DNA duplexes, triggering genetic instability.

Atypical structures of GAA/TTC trinucleotide repeats underlying Friedreich's ataxia: DNA triplexes and RNA/DNA hybrids

Expansion of the GAA/TTC repeats in the first intron of the FXN gene causes Friedreich's ataxia. Non-canonical structures are linked to this expansion. DNA triplexes and R-loops are believed to arrest transcription, which results in frataxin deficiency and eventual neurodegeneration. We present a systematic in silico characterization of the possible DNA triplexes that could be assembled with GAA and TTC strands; the two hybrid duplexes [r(GAA):d(TTC) and d(GAA):r(UUC)] in an R-loop; and three hybrid triplexes that could form during bidirectional transcription when the non-template DNA strand bonds with the hybrid duplex (collapsed R-loops, where the two DNA strands remain antiparallel). For both Y·R:Y and R·R:Y DNA triplexes, the parallel third strand orientation is more stable; both parallel and antiparallel protonated d(GA+A)·d(GAA):d(TTC) triplexes are stable. Apparent contradictions in the literature about the R·R:Y triplex stability is probably due to lack of molecular resolution, since shifting the third strand by a single nucleotide alters the stability ranking. In the collapsed R-loops, antiparallel d(TTC+)·d(GAA):r(UUC) is unstable, while parallel d(GAA)·r(GAA):d(TTC) and d(GA+A)·r(GAA):d(TTC) are stable. In addition to providing new structural perspectives for specific therapeutic aims, our results contribute to a systematic structural basis for the emerging field of quantitative R-loop biology.

GC-Content Dependence of Elastic and Overstretching Properties of DNA:RNA Hybrid Duplexes

DNA:RNA hybrid duplex plays important roles in various biological processes. Both its structural stability and interactions with proteins are highly sequence dependent. In this study, we utilize homebuilt optical tweezers to investigate how GC contents in the sequence influence the structural and mechanical properties of DNA:RNA hybrid by measuring its contour length, elasticities, and overstretching dynamics. Our results support that the DNA:RNA hybrid adopts a conformation between the A- and B-form helix, and the GC content does not affect its structural and elastic parameters obviously when varying from 40 to 60% before the overstretching transition of DNA:RNA hybrid occurs. In the overstretching transition, however, our study unravels significant heterogeneity and strong sequence dependence, suggesting the presence of a highly dynamic competition between the two processes, namely the S-form duplex formation (nonhysteretic) and the unpeeling (hysteretic). Analyzing the components left in DNA:RNA hybrid after the overstretching transition suggests that the RNA strand is more easily unpeeled than the DNA strand, whereas an increase in the GC content from 40 to 60% can significantly reduce unpeeling. Large hysteresis is observed between the stretching and relaxation processes, which is also quantitatively correlated with the percentage of unpeeling in the DNA:RNA duplex. Increasing in both the salt concentration and GC content can effectively reduce the hysteresis with the latter being more significant. Together, our study reveals that the mechanical properties of DNA:RNA hybrid duplexes are significantly different from double-stranded DNA and double-stranded RNA, and its overstretching behavior is highly sequence dependent. These results should be taken into account in the future studies on DNA:RNA-hybrid-related functional structures and motor proteins.

Stability prediction of canonical and non-canonical structures of nucleic acids in various molecular environments and cells

This review provides the biophysicochemical background and recent advances in stability prediction of canonical and non-canonical structures of nucleic acids in various molecular environments and cells.

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.

Nearest-neighbour parameters optimized for melting temperature prediction of DNA/RNA hybrids at high and low salt concentrations

Gene editing technologies sparked a renewed interest in the hybridization of DNA/RNA duplexes, yet little improvement on nearest-neighbour parameters was made over the past two decades. For low sodium concentration no parameter set was yet calculated. Here, we revised the existing experimental datasets and used an expanded set of sequences from which we recalculated the nearest-neighbour parameters, reducing the average temperature prediction uncertainty to 1.6 °C. Two experimental sets using temperatures extracted via different methods were used with similar results, with the curve-fitting method achieving a slight advantage in prediction quality over other methods. Additionally, we obtained new parameters for low salt with an average uncertainty of 0.98 °C. We also tested several types of salt correction factors and concluded that it is advisable to use those originally developed for RNA/RNA rather than for DNA/DNA.

Utilizing guanine-coordinated Zn2+ ions to determine DNA crystal structures by single-wavelength anomalous diffraction

The experimental phase determination of crystal structures of nucleic acids and nucleic acid-ligand complexes would benefit from a facile method. Even for double-stranded DNA, software-generated models are generally insufficiently accurate to serve as molecular replacement search models, necessitating experimental phasing. Here, it is demonstrated that Zn²⁺ ions coordinated to the N7 atom of guanine bases generate sufficient anomalous signal for single-wavelength anomalous diffraction (SAD) phasing of DNA crystal structures. Using zinc SAD, three crystal structures of double-stranded DNA oligomers, 5'-AGGGATCCCT-3', 5'-GGGATCCC-3' and 5'-GAGGCCTC-3', were determined. By determining the crystal structure of one of these oligomers, GAGGCCTC, in the presence of Mg²⁺ instead of Zn²⁺, it was demonstrated that Zn²⁺ is not structurally perturbing. These structures allowed the analysis of structural changes in the DNA on the binding of analogues of the natural product mithramycin to two of these oligomers, AGGGATCCCT and GAGGCCTC. Zinc SAD may become a routine approach for determining the crystal structures of nucleic acids and their complexes with small molecules.

The Mechanical Properties of RNA-DNA Hybrid Duplex Stretched by Magnetic Tweezers

RNA can anneal to its DNA template to generate an RNA-DNA hybrid (RDH) duplex and a displaced DNA strand, termed R-loop. RDH duplex occupies up to 5% of the mammalian genome and plays important roles in many biological processes. The functions of RDH duplex are affected by its mechanical properties, including the elasticity and the conformation transitions. The mechanical properties of RDH duplex, however, are still unclear. In this work, we studied the mechanical properties of RDH duplex using magnetic tweezers in comparison with those of DNA and RNA duplexes with the same sequences. We report that the contour length of RDH duplex is ∼0.30 nm/bp, and the stretching modulus of RDH duplex is ∼660 pN, neither of which is sensitive to NaCl concentration. The persistence length of RDH duplex depends on NaCl concentration, decreasing from ∼63 nm at 1 mM NaCl to ∼49 nm at 500 mM NaCl. Under high tension of ∼60 pN, the end-opened RDH duplex undergoes two distinct overstretching transitions; at high salt in which the basepairs are stable, it undergoes the nonhysteretic transition, leading to a basepaired elongated structure, whereas at low salt, it undergoes a hysteretic peeling transition, leading to the single-stranded DNA strand under force and the single-stranded RNA strand coils. The peeled RDH is difficult to reanneal back to the duplex conformation, which may be due to the secondary structures formed in the coiled single-stranded RNA strand. These results help us understand the full picture of the structures and mechanical properties of nucleic acid duplexes in solution and provide a baseline for studying the interaction of RDH with proteins at the single-molecule level.

New high copy tandem repeat in the content of the chicken W chromosome

The content of repetitive DNA in avian genomes is considerably less than in other investigated vertebrates. The first descriptions of tandem repeats were based on the results of routine biochemical and molecular biological experiments. Both satellite DNA and interspersed repetitive elements were annotated using library-based approach and de novo repeat identification in assembled genome. The development of deep-sequencing methods provides datasets of high quality without preassembly allowing one to annotate repetitive elements from unassembled part of genomes. In this work, we search the chicken assembly and annotate high copy number tandem repeats from unassembled short raw reads. Tandem repeat (GGAAA)_n has been identified and found to be the second after telomeric repeat (TTAGGG)n most abundant in the chicken genome. Furthermore, (GGAAA)_n repeat forms expanded arrays on the both arms of the chicken W chromosome. Our results highlight the complexity of repetitive sequences and update data about organization of sex W chromosome in chicken.

DNA editing in DNA/RNA hybrids by adenosine deaminases that act on RNA

Adenosine deaminases that act on RNA (ADARs) carry out adenosine (A) to inosine (I) editing reactions with a known requirement for duplex RNA. Here, we show that ADARs also react with DNA/RNA hybrid duplexes. Hybrid substrates are deaminated efficiently by ADAR deaminase domains at dA-C mismatches and with E to Q mutations in the base flipping loop of the enzyme. For a long, perfectly matched hybrid, deamination is more efficient with full length ADAR2 than its isolated deaminase domain. Guide RNA strands for directed DNA editing by ADAR were used to target six different 2΄-deoxyadenosines in the M13 bacteriophage ssDNA genome. DNA editing efficiencies varied depending on the sequence context of the editing site consistent with known sequence preferences for ADARs. These observations suggest the reaction within DNA/RNA hybrids may be a natural function of human ADARs. In addition, this work sets the stage for development of a new class of genome editing tools based on directed deamination of 2΄-deoxyadenosines in DNA/RNA hybrids.

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.

Inclusion of methoxy groups inverts the thermodynamic stabilities of DNA-RNA hybrid duplexes: A molecular dynamics simulation study

Modified nucleic acids have found profound applications in nucleic acid based technologies such as antisense and antiviral therapies. Previous studies on chemically modified nucleic acids have suggested that modifications incorporated in furanose sugar especially at 2'-position attribute special properties to nucleic acids when compared to other modifications. 2'-O-methyl modification to deoxyribose sugars of DNA-RNA hybrids is one such modification that increases nucleic acid stability and has become an attractive class of compounds for potential antisense applications. It has been reported that modification of DNA strands with 2'-O-methyl group reverses the thermodynamic stability of DNA-RNA hybrid duplexes. Molecular dynamics simulations have been performed on two hybrid duplexes (DR and RD) which differ from each other and 2'-O-methyl modified counterparts to investigate the effect of 2'-O-methyl modification on their duplex stability. The results obtained suggest that the modification drives the conformations of both the hybrid duplexes towards A-RNA like conformation. The modified hybrid duplexes exhibit significantly contrasting dynamics and hydration patterns compared to respective parent duplexes. In line with the experimental results, the relative binding free energies suggest that the introduced modifications stabilize the less stable DR hybrid, but destabilize the more stable RD duplex. Binding free energy calculations suggest that the increased hydrophobicity is primarily responsible for the reversal of thermodynamic stability of hybrid duplexes. Free energy component analysis further provides insights into the stability of modified duplexes.

The role of Mg(II) in DNA cleavage site recognition in group II intron ribozymes: solution structure and metal ion binding sites of the RNA-DNA complex

Group II intron ribozymes catalyze the cleavage of (and their reinsertion into) DNA and RNA targets using a Mg2(+)-dependent reaction. The target is cleaved 3' to the last nucleotide of intron binding site 1 (IBS1), one of three regions that form base pairs with the intron's exon binding sites (EBS1 to -3).We solved the NMR solution structure of the d3' hairpin of the Sc.ai5γ intron containing EBS1 in its 11-nucleotide loop in complex with the dIBS1 DNA 7-mer and compare it with the analogous RNA-RNA contact. The EBS1-dIBS1 helix is slightly flexible and non-symmetric. NMR data reveal two major groove binding sites for divalent metal ions at the EBS1-dIBS1 helix, and surface plasmon resonance experiments show that low concentrations of Mg2(+) considerably enhance the affinity of dIBS1 for EBS1. Our results indicate that identification of both RNA and DNA IBS1 targets, presentation of the scissile bond, and stabilization of the structure by metal ions are governed by the overall structure of EBS1-dIBS1 and the surrounding loop nucleotides but are irrespective of different EBS1-(d)IBS1 geometries and interstrand affinities.

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.

Structure-based DNA-targeting strategies with small molecule ligands for drug discovery

Nucleic acids are the molecular targets of many clinical anticancer drugs. However, compared with proteins, nucleic acids have traditionally attracted much less attention as drug targets in structure-based drug design, partially because limited structural information of nucleic acids complexed with potential drugs is available. Over the past several years, enormous progresses in nucleic acid crystallization, heavy-atom derivatization, phasing, and structural biology have been made. Many complicated nucleic acid structures have been determined, providing new insights into the molecular functions and interactions of nucleic acids, especially DNAs complexed with small molecule ligands. Thus, opportunities have been created to further discover nucleic acid-targeting drugs for disease treatments. This review focuses on the structure studies of DNAs complexed with small molecule ligands for discovering lead compounds, drug candidates, and/or therapeutics.

Synthesis of novel di-Se-containing thymidine and Se-DNAs for structure and function studies

The selenium derivatization of nucleic acids and nucleic acid-protein complexes has provided a powerful tool to solve phase problem in X-ray crystallography. Selenium atoms in the nucleotides can serve as fine scattering centers in crystal diffraction. Towards the synthesis of multiple selenium atom-containing nucleotides, which offers strong phasing power to facilitate crystal structure determination, we report here the synthesis of the thymidine analogue containing two Se atoms in one nucleobase. The novel Se-containing nucleoside and oligonucleotide DNAs were synthesized and found with the red-shifted UV spectrum and yellow color. Their unique properties are useful in phase determination, nucleic acid-based detection as well as spectroscopic studies of nucleic acids and nucleic acid-protein complexes.

Conformational flexibility influences degree of hydration of nucleic acid hybrids

Four nucleic acid duplexes-DNA/RNA hybrid, RNA/DNA hybrid, RNA duplex, and DNA duplex-were studied under molecular crowding conditions of osmolytes. Destabilization of duplexes (ΔΔG°(25)) indicated that the ΔΔG°(25) values of hybrids were intermediate between those of DNA and RNA duplexes. In the presence of polyethylene glycol 200, the ΔΔG°(25) values were estimated to be +3.0, +3.5, +3.5, and +4.1 kcal mol(-1) for the DNA duplex, DNA/RNA hybrid, RNA/DNA hybrid, and RNA duplex, respectively. Differences in the number of water molecules taken up (-Δn(w)) upon duplex formations between 0 and 37 °C (Δ(-Δn(w))) were estimated to be 44.8 and 59.7 per duplex structure for the DNA/RNA and RNA/DNA hybrids, respectively. While the Δ(-Δn(w)) value for the DNA/RNA hybrid was intermediate between those of the DNA (26.1) and RNA (59.2) duplexes, the value for RNA/DNA hybrid was close to that of RNA duplex. These differences in the thermodynamic parameters and hydration are probably a consequence of the enhanced global flexibility of the RNA/DNA hybrid structure relative to the DNA/RNA hybrid structure observed in molecular dynamics simulations. This molecular crowding study provides information not only on hydration but also on the flexibility of the conformation of nucleic acid duplexes.

Selenium derivatization of nucleic acids for X-ray crystal-structure and function studies

It is estimated that over two thirds of all new crystal structures of proteins are determined via the protein selenium derivatization (selenomethionine (Se-Met) strategy). This selenium derivatization strategy via MAD (multi-wavelength anomalous dispersion) phasing has revolutionized protein X-ray crystallography. Through our pioneer research, similarly, Se has also been successfully incorporated into nucleic acids to facilitate the X-ray crystal-structure and function studies of nucleic acids. Currently, Se has been stably introduced into nucleic acids by replacing nucleotide O-atom at the positions 2', 4', 5', and in nucleobases and non-bridging phosphates. The Se derivatization of nucleic acids can be achieved through solid-phase chemical synthesis and enzymatic methods, and the Se-derivatized nucleic acids (SeNA) can be easily purified by HPLC, FPLC, and gel electrophoresis to obtain high purity. It has also been demonstrated that the Se derivatization of nucleic acids facilitates the phase determination via MAD phasing without significant perturbation. A growing number of structures of DNAs, RNAs, and protein-nucleic acid complexes have been determined by the Se derivatization and MAD phasing. Furthermore, it was observed that the Se derivatization can facilitate crystallization, especially when it is introduced to the 2'-position. In addition, this novel derivatization strategy has many advantages over the conventional halogen derivatization, such as more choices of the modification sites via the atom-specific substitution of the nucleotide O-atom, better stability under X-ray radiation, and structure isomorphism. Therefore, our Se-derivatization strategy has great potentials to provide rational solutions for both phase determination and high-quality crystal growth in nucleic-acid crystallography. Moreover, the Se derivatization generates the nucleic acids with many new properties and creates a new paradigm of nucleic acids. This review summarizes the recent developments of the atomic site-specific Se derivatization of nucleic acids for structure determination and function study. Several applications of this Se-derivatization strategy in nucleic acid and protein research are also described in this review.

Dynamics and stability of individual base pairs in two homologous RNA-DNA hybrids

Nuclear magnetic resonance spectroscopy and proton exchange have been used to characterize two RNA-DNA hybrids from the tR2 intrinsic transcription terminator site of phage lambda. The hybrids have the same base sequence [5'-GGCGCAGGCC(T/U)(T/U)CC-3'/5'-GGAAGGCC(T/U)GCGCC-3'] but differ from each other by an interchange of DNA and RNA strands. The opening of single base pairs in the two hybrids is characterized by measuring the rates of exchange of imino protons with solvent protons as a function of the concentration of a proton acceptor (ammonia base) at 10 degrees C. The free energy change in the opening reaction provides a measure of the stability of the base pair, while the rates of opening and closing define the base pair dynamics. The results demonstrate that, within the same base sequence context, dA-rU base pairs are less stable than dT-rA base pairs. The differences in stability are enhanced when two dA-rU base pairs are located next to each other in the hybrid structure. For the G-C base pairs, the rates of opening and closing and the stability are affected by the base sequence context and by the nature of the sugar moiety attached to the guanine. The dominant feature of the base sequence is the proximity of the dA-rU base pair, which destabilizes the G-C base pair when the guanine is located on the DNA strand. Two G-C base pairs (namely, those in the fourth and 10th positions) exhibit large differences in their opening and closing rates between the two hybrids, while maintaining the same stability. These results provide the first demonstration that, for RNA-DNA hybrid structures with the same base sequence, the opening dynamics and the stability of individual base pairs are strongly influenced by the chemical nature of each strand.

A preliminary neutron crystallographic study of an A-DNA crystal

The LADI-III diffractometer at the Institut Laue-Langevin has been used to carry out a preliminary neutron crystallographic study of the self-complementary DNA oligonucleotide d(AGGGGCCCCT)(2) in the A conformation. The results demonstrate the viability of a full neutron crystallographic analysis with the aim of providing enhanced information on the ion-water networks that are known to be important in stabilizing A-DNA. This is the first account of a single-crystal neutron diffraction study of A-DNA. The study was carried out with the smallest crystal used to date for a neutron crystallographic study of a biological macromolecule.

Orientation dependence in fluorescent energy transfer between Cy3 and Cy5 terminally attached to double-stranded nucleic acids

We have found that the efficiency of fluorescence resonance energy transfer between Cy3 and Cy5 terminally attached to the 5' ends of a DNA duplex is significantly affected by the relative orientation of the two fluorophores. The cyanine fluorophores are predominantly stacked on the ends of the helix in the manner of an additional base pair, and thus their relative orientation depends on the length of the helix. Observed fluorescence resonance energy transfer (FRET) efficiency depends on the length of the helix, as well as its helical periodicity. By changing the helical geometry from B form double-stranded DNA to A form hybrid RNA/DNA, a marked phase shift occurs in the modulation of FRET efficiency with helix length. Both curves are well explained by the standard geometry of B and A form helices. The observed modulation for both polymers is less than that calculated for a fully rigid attachment of the fluorophores. However, a model involving lateral mobility of the fluorophores on the ends of the helix explains the observed experimental data. This has been further modified to take account of a minor fraction of unstacked fluorophore observed by fluorescent lifetime measurements. Our data unequivocally establish that Förster transfer obeys the orientation dependence as expected for a dipole-dipole interaction.

Selenium derivatization of nucleic acids for phase and structure determination in nucleic acid X-ray crystallography

Selenium derivatization (via selenomethionine) of proteins for crystal structure determination via MAD phasing has revolutionized protein X-ray crystallography. It is estimated that over two thirds of all new crystal structures of proteins have been determined via Se-Met derivatization. Similarly, selenium functionalities have also been successfully incorporated into nucleic acids to facilitate their structure studies and it has been proved that this Se-derivatization has advantages over halogen strategy, which was usually used as a traditional method in this field. This review reports the development of site-specific selenium derivatization of nucleic acids for phase determination since the year of 2001 (mainly focus on the 2'-position of the ribose). All the synthesis of 2'-SeMe modified phosphoramidite building blocks (U, C, T, A, G) and the according oligonucleotides are included. In addition, several structures of selenium contained nucleic acid are also described in this paper.

Theoretical analysis of antisense duplexes: determinants of the RNase H susceptibility

The structure and dynamic properties of different antisense related duplexes (DNA x RNA, 2'O-Me-DNA x RNA, 2'F-ANA x RNA, C5(Y)-propynyl-DNA x RNA, ANA x RNA, and control duplexes DNA x DNA and RNA x RNA) have been determined by means of long molecular dynamics simulations (covering more than 0.5 micros of fully solvated unrestrained MD simulation). The massive analysis presented here allows us to determine the subtle differences between the different duplexes, which in all cases pertain to the same structural family. This analysis provides information on the molecular determinants that allow RNase H to recognize and degrade some of these duplexes, whereas others with apparently similar conformations are not affected. Subtle structural and deformability features define the key properties used by RNase H to discriminate between duplexes.

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.

Crystal structures of DNA:DNA and DNA:RNA duplexes containing 5-(N-aminohexyl)carbamoyl-modified uracils reveal the basis for properties as antigene and antisense molecules

Oligonucleotides containing 5-(N-aminohexyl)carbamoyl-modified uracils have promising features for applications as antigene and antisense therapies. Relative to unmodified DNA, oligonucleotides containing 5-(N-aminohexyl)carbamoyl-2′-deoxyuridine (^NU) or 5-(N-aminohexyl)carbamoyl-2′-O-methyluridine (^NU_m), respectively exhibit increased binding affinity for DNA and RNA, and enhanced nuclease resistance. To understand the structural implications of ^NU and ^NU_m substitutions, we have determined the X-ray crystal structures of DNA:DNA duplexes containing either ^NU or ^NU_m and of DNA:RNA hybrid duplexes containing ^NU_m. The aminohexyl chains are fixed in the major groove through hydrogen bonds between the carbamoyl amino groups and the uracil O4 atoms. The terminal ammonium cations on these chains could interact with the phosphate oxygen anions of the residues in the target strands. These interactions partly account for the increased target binding affinity and nuclease resistance. In contrast to ^NU, ^NU_m decreases DNA binding affinity. This could be explained by the drastic changes in sugar puckering and in the minor groove widths and hydration structures seen in the ^NU_m containing DNA:DNA duplex structure. The conformation of ^NU_m, however, is compatible with the preferred conformation in DNA:RNA hybrid duplexes. Furthermore, the ability of ^NU_m to render the duplexes with altered minor grooves may increase nuclease resistance and elicit RNase H activity.

Selenium derivatization of nucleic acids for crystallography

The high-resolution structure of the DNA (5'-GTGTACA-C-3') with the selenium derivatization at the 2'-position of T2 was determined via MAD and SAD phasing. The selenium-derivatized structure (1.28 A resolution) with the 2'-Se modification in the minor groove is isomorphorous to the native structure (2.0 A). To directly compare with the conventional bromine derivatization, we incorporated bromine into the 5-postion of T4, determined the bromine-derivatized DNA structure at 1.5 A resolution, and found that the local backbone torsion angles and solvent hydration patterns were altered in the structure with the Br incorporation in the major groove. Furthermore, while the native and Br-derivatized DNAs needed over a week to form reasonable-size crystals, we observed that the Se-derivatized DNAs grew crystals overnight with high-diffraction quality, suggesting that the Se derivatization facilitated the crystal formation. In addition, the Se-derivatized DNA sequences crystallized under a broader range of buffer conditions, and generally had a faster crystal growth rate. Our experimental results indicate that the selenium derivatization of DNAs may facilitate the determination of nucleic acid X-ray crystal structures in phasing and high-quality crystal growth. In addition, our results suggest that the Se derivatization can be an alternative to the conventional Br derivatization.

Directed evolution of DNA polymerase, RNA polymerase and reverse transcriptase activity in a single polypeptide

DNA polymerases enable key technologies in modern biology but for many applications, native polymerases are limited by their stringent substrate recognition. Here we describe short-patch compartmentalized self-replication (spCSR), a novel strategy to expand the substrate spectrum of polymerases in a targeted way. spCSR is based on the previously described CSR, but unlike CSR only a short region (a "patch") of the gene under investigation is diversified and replicated. This allows the selection of polymerases under conditions where catalytic activity and processivity are compromised to the extent that full self-replication is inefficient. We targeted two specific motifs involved in substrate recognition in the active site of DNA polymerase I from Thermus aquaticus (Taq) and selected for incorporation of both ribonucleotide- (NTP) and deoxyribonucleotide-triphosphates (dNTPs) using spCSR. This allowed the isolation of multiple variants of Taq with apparent dual substrate specificity. They were able to synthesize RNA, while still retaining essentially wild-type (wt) DNA polymerase activity as judged by PCR. One such mutant (AA40: E602V, A608V, I614M, E615G) was able to incorporate both NTPs and dNTPs with the same catalytic efficiency as the wt enzyme incorporates dNTPs. AA40 allowed the generation of mixed RNA-DNA amplification products in PCR demonstrating DNA polymerase, RNA polymerase as well as reverse transcriptase activity within the same polypeptide. Furthermore, AA40 displayed an expanded substrate spectrum towards other 2'-substituted nucleotides and was able to synthesize nucleic acid polymers in which each base bore a different 2'-substituent. Our results suggest that spCSR will be a powerful strategy for the generation of polymerases with altered substrate specificity for applications in nano- and biotechnology and in the enzymatic synthesis of antisense and RNAi probes.

Two-dimensional strandness-dependent electrophoresis: A method to characterize single-stranded DNA, double-stranded DNA, and RNA–DNA hybrids in complex samples

We describe two-dimensional strandness-dependent electrophoresis (2D-SDE) for quantification and length distribution analysis of single-stranded (ss) DNA fragments, double-stranded (ds) DNA fragments, RNA-DNA hybrids, and nicked DNA fragments in complex samples. In the first dimension nucleic acid molecules are separated based on strandness and length in the presence of 7 M urea. After the first-dimension electrophoresis all nucleic acid fragments are heat denatured in the gel. During the second-dimension electrophoresis all nucleic acid fragments are single-stranded and migrate according to length. 2D-SDE takes about 90 min and requires only basic skills and equipment. We show that 2D-SDE has many applications in analyzing complex nucleic acid samples including (1) estimation of renaturation efficiency and kinetics, (2) monitoring cDNA synthesis, (3) detection of nicked DNA fragments, and (4) estimation of quality and in vitro damage of nucleic acid samples. Results from 2D-SDE should be useful to validate techniques such as complex polymerase chain reaction, subtractive hybridization, cDNA synthesis, cDNA normalization, and microarray analysis. 2D-SDE could also be used, e.g., to characterize biological nucleic acid samples. Information obtained with 2D-SDE cannot be readily obtained with other methods. 2D-SDE can be used for preparative isolation of ssDNA fragments, dsDNA fragments, and RNA-DNA hybrids.

Recognition of RNA by amide modified backbone nucleic acids: molecular dynamics simulations of DNA-RNA hybrids in aqueous solution

Thermodynamic and structural properties of a chemically modified DNA-RNA hybrid in which a phosphodiester linkage is replaced by a neutral amide-3 linkage (3'-CH(2)-CONH-5') were investigated using UV melting experiments, molecular dynamics simulations in explicit water, and continuum solvent models. van't Hoff analysis of the experimental UV melting curves suggests that the significant increase of the thermodynamic stability of a 15-mer DNA-RNA with seven alternated amide-3 modifications (+11 degrees C) is mainly due to an increased binding enthalpy. To further evaluate the origin in the observed affinities differences, the electrostatic contribution to the binding free energy was calculated by solving the Poisson-Boltzmann equation numerically. The nonelectrostatic contribution was estimated as the product of a hydrophobic surface tension coefficient and the surface area that is buried upon double strand formation. Structures were taken from 10 ns molecular dynamics simulations computed in a consistent fashion using explicit solvent, counterions, and the particle-mesh Ewald procedure. The present preliminary thermodynamic study suggests that the favorable binding free energy of the amide-3 DNA single strand to the complementary RNA is equally driven by electrostatic and nonpolar contributions to the binding compared to their natural analogues. In addition, molecular dynamics simulations in explicit water were performed on an amide-3 DNA single strand and the corresponding natural DNA. Results from the conformations cluster analysis of the simulated amide-3 DNA single strand ensembles suggest that the 25% of the population sampled within 10 ns has a pre-organized conformation where the sugar C3' endo pucker is favored at the 3'-flanking nucleotides. These structural and thermodynamic features contribute to the understanding of the observed increased affinities of the amide-3 DNA-RNA hybrids at the microscopic level.

Structure, Recognition Properties, and Flexibility of the DNA·RNA Hybrid

Molecular dynamics is used to investigate the properties of the DNA.RNA hybrid in aqueous solution at room temperature. The structure of the hybrid is intermediate between A and B forms but, in general, closer to the canonical A-type helix. All the riboses exhibit North puckerings, while 2'-deoxyriboses exist in North, East, and South puckerings, the latter being the most populated one. The molecular recognition pattern of the DNA.RNA hybrid is a unique combination of those of normal DNA and RNA duplexes. Finally, the results obtained from essential dynamics and stiffness analysis demonstrate the large and very asymmetric flexibility of the hybrid and the strong predilection that each strand (DNA or RNA) has on the nature of their intrinsic motions in the corresponding homoduplexes. The implications of the unique structural and dynamic properties of the DNA.RNA hybrid on the mechanism of cleavage by RNase H are discussed.

Sequence-specific artificial ribonucleases. I. Bis-imidazole-containing oligonucleotide conjugates prepared using precursor-based strategy

Antisense oligonucleotide conjugates, bearing constructs with two imidazole residues, were synthesized using a precursor-based technique employing post-synthetic histamine functionalization of oligonucleotides bearing methoxyoxalamido precursors at the 5'-termini. The conjugates were assessed in terms of their cleavage activities using both biochemical assays and conformational analysis by molecular modelling. The oligonucleotide part of the conjugates was complementary to the T-arm of yeast tRNA(Phe) (44-60 nt) and was expected to deliver imidazole groups near the fragile sequence C61-ACA-G65 of the tRNA. The conjugates showed ribonuclease activity at neutral pH and physiological temperature resulting in complete cleavage of the target RNA, mainly at the C63-A64 phosphodiester bond. For some constructs, cleavage was completed within 1-2 h under optimal conditions. Molecular modelling was used to determine the preferred orientation(s) of the cleaving group(s) in the complexes of the conjugates with RNA target. Cleaving constructs bearing two imidazole residues were found to be conformationally highly flexible, adopting no preferred specific conformation. No interactions other than complementary base pairing between the conjugates and the target were found to be the factors stabilizing the 'active' cleaving conformation(s).

Structural analysis of a eukaryotic sliding DNA clamp–clamp loader complex

Sliding clamps are ring-shaped proteins that encircle DNA and confer high processivity on DNA polymerases. Here we report the crystal structure of the five-protein clamp loader complex (replication factor-C, RFC) of the yeast Saccharomyces cerevisiae, bound to the sliding clamp (proliferating cell nuclear antigen, PCNA). Tight interfacial coordination of the ATP analogue ATP-gammaS by RFC results in a spiral arrangement of the ATPase domains of the clamp loader above the PCNA ring. Placement of a model for primed DNA within the central hole of PCNA reveals a striking correspondence between the RFC spiral and the grooves of the DNA double helix. This model, in which the clamp loader complex locks onto primed DNA in a screw-cap-like arrangement, provides a simple explanation for the process by which the engagement of primer-template junctions by the RFC:PCNA complex results in ATP hydrolysis and release of the sliding clamp on DNA.

Acyclic, achiral enamide nucleoside analogues. The importance of the C=C bond in the analogue for its ability to mimic natural nucleosides

The conformations of an acyclic, achiral enamide thymidine analogue 1 have been studied by model building and geometry calculations, as well as by NMR NOE and UV experiments. The results indicate that there are no significant barriers to rotation around any of the sigma bonds, in particular the N1-C1' enamide bond, and that the analogue should be able to accommodate conformations that mimic the conformations of natural nucleosides in A- and B-type helices quite well. For comparison the saturated analogue 2 has been prepared and built into oligonucleotides. It is shown that incorporation of 2 in oligonucleotides results in a much larger depression of the melting temperature (deltaTm -10 to -12.5 degrees C) than does incorporation of 1 (deltaTm -5 to -6.5 degrees C).

Enzymatic synthesis of phosphoroselenoate DNA using thymidine 5'-(alpha-P-seleno)triphosphate and DNA polymerase for X-ray crystallography via MAD

We report here the first study of enzymatic synthesis of two phosphoroselenoate (PSe) DNAs using the two alpha-Se-TTP diastereomers (Sp and Rp) and DNA polymerase. The experimental results indicate that Klenow equally recognizes the two individual diastereomers at the same level as natural TTP. The incorporations of the PSe groups at the expected sites have been confirmed by the digestion resistance to exonuclease III, and the different patterns of the digestion resistance of DNA I and II indicate the configurational differences of the PSe centers (Sp or Rp). Unlike chemical synthesis, which is limited to short DNAs and where the separation of the PSe DNA diastereomers is necessary, this enzymatic method can be used to prepare longer DNAs without diastereomer separation. This quantitative enzymatic approach is particular valuable for the synthesis of longer DNAs with multiple PSe groups in large scale for their X-ray crystal structure determination by the MAD phasing technique.

Synthesis of Selenium-Derivatized Cytidine and Oligonucleotides for X-ray Crystallography Using MAD

Synthesis of the novel 2'-Se-cytidine phosphoramidite was achieved via transformation of the uridine analogue to the cytidine derivative in high yield. This 2'-Se-cytidine phosphoramidite was used to synthesize selenium-derivatized DNA and RNA oligonucleotides for X-ray crystallography using MAD. The nucleotide coupling yield using this novel phosphoramidite was over 99% when 5-benzylmercaptotetrazole (5-BMT) was used as the coupling reagent.

A new type of acyclic, achiral nucleoside analogue. How does it simulate nucleosides?

The new monomer 1 seems to be an excellent mimic of nucleosides with different sugar conformations (north, south, and envelope), because of the relatively free rotation around gamma, delta, and chi. The rotation around chi is primarily controlled by the repulsion between H6 and the two hydrogen atoms on C4' and not pi conjugation between the double bond and the nucleobase. A viable synthesis of the guanine monomer 8 is described.

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

An intermolecular intercalation of base pairs was found at the CA step in the I222 crystal structure of the RNA.DNA hybrid, r(CAAAGAAAAG).d(CTTTTCTTTG), which contains two-thirds of the polypurine tract sequence of HIV-1 with a substitution of cytosine for the initial adenine. This sequence crystallized in both P212121 and I222 space groups, with an rms difference of only 0.63 A between residues 3 to 18 of the two forms. P212121 and I222 helices are both A-like, but intercalation occurs only in the I222 crystal form. The present structure shows bases stacked in parallel rather than perpendicular as in intercalated DNA (I-DNA). The base intercalation is also different from zipper-like meshing of bases seen in the center of the crystal structure of d(GCGAAAGCT), which does not have Watson-Crick base pairing. The base-step intercalation seen here is reminiscent of domain swapping in proteins; therefore, we call this phenomenon "base-pair swapping." It involves a highly mobile CA step and seems to be sequence-specific and electrostatically stable without disrupting Watson-Crick interactions. It also exhibits a large rise concurrent with unwinding of the helix (low twist). We present a base-pair swapping dimer in nucleic acids.

The solution structure of a DNA*RNA duplex containing 5-propynyl U and C; comparison with 5-Me modifications

The addition of the propynyl group at the 5 position of pyrimidine nucleotides is highly stabilising. We have determined the thermodynamic stability of the DNA.RNA hybrid r(GAAGAGAAGC)*d(GC(p)U(p)U(p)C(p)U(p) C(p)U(p)U(p)C) where p is the propynyl group at the 5 position and compared it with that of the unmodified duplex and the effects of methyl substitutions. The incorporation of the propyne group at the 5 position gives rise to a very large stabilisation of the hybrid duplex compared with the analogous 5-Me modification. The duplexes have been characterised by gel electrophoresis and NMR spectroscopy, which indicate that methyl substitutions have a smaller influence on local and global conformation than the propynyl groups. The increased NMR spectral dispersion of the propyne-modified duplex allowed a larger number of experimental restraints to be measured. Restrained molecular dynamics in a fully solvated system showed that the propyne modification leads to substantial conformational rearrangements stabilising a more A-like structure. The propynyl groups occupy a large part of the major groove and make favourable van der Waals interactions with their nearest neighbours and the atoms of the rings. This enhanced overlap may account at least in part for the increased thermodynamic stability. Furthermore, the simulations show a spine of hydration in the major groove as well as in the minor groove involving the RNA hydroxyl groups.