NMR structure determination of a modified DNA oligonucleotide containing a new intercalating nucleic acid

The intercalating nucleic acid (INA) presented in this paper is a novel 1-O-(1-pyrenylmethyl)glycerol DNA intercalator that induces high thermal affinity for complementary DNA. The duplex examined contained two INA intercalators, denoted X, inserted directly opposite each other: d(C(1)T(2)C(3)A(4)A(5)C(6)X(7)C(8)A(9)A(10)G(11)C(12)T(13)):d(A(14)G(15)C(16)T(17)-T(18)G(19)X(20)G(21)T(22)T(23)G(24)A(25)G(26)). Unlike most other nucleotide analogues, DNA with INA inserted has a lower affinity for hybridizing to complementary DNA with an INA inserted directly opposite than to complementary unmodified DNA. In this study we used two-dimensional (1)H NMR spectroscopy to determine a high-resolution solution structure of the weak INA-INA duplex. A modified ISPA approach was used to obtain interproton distance bounds from NOESY cross-peak intensities. These distance bounds were used as restraints in molecular dynamics (rMD) calculations. Twenty final structures were generated for the duplex from a B-type DNA starting structure. The root-mean-square deviation (RMSD) of the coordinates for the 20 structures of the complex was 1.95 A. This rather large value, together with broad lines in the area of insertion, reflect the high degree of internal motion in the complex. The determination of the structure revealed that both intercalators were situated in the center of the helix, stacking with each other and the neighboring nucleobases. The intercalation of the INAs caused an unwinding of the helix in the insertion area, creating a ladderlike structure. The structural changes observed upon intercalation were mainly of local character; however, a broadening of the minor groove was found throughout the helix.     

RNAPack: an integrated NMR approach to RNA structure determination

Over the last decade, a vast number of useful nuclear magnetic resonance (NMR) experiments have been developed and successfully employed to determine the structure and dynamics of RNA oligonucleotides. Despite this progress, high-resolution RNA structure determination by NMR spectroscopy still remains a lengthy process and requires programming and extensive calibrations to perform NMR experiments successfully. To accelerate RNA structure determination by NMR spectroscopy, we have designed and programmed a package of RNA NMR experiments, called RNAPack. The user-friendly package contains a set of semiautomated single, double, and triple resonance NMR experiments, which are fully optimized for high-resolution RNA solution structure determination on Varian NMR spectrometers. RNAPack provides an autocalibration feature that allows rapid calibration of all NMR experiments in a single step and thereby speeds up the NMR data collection and eliminates user errors. In our laboratory, we have successfully employed this technology to solve RNA solution structures of domains of the internal ribosome entry site of the genomic hepatitis C viral RNA in less than 3 months. RNAPack therefore makes NMR spectroscopy an attractive and rapid structural tool and allows integration of atomic resolution structural information into biochemical studies of large RNA systems.     

NMR structure of an antisense DNA.RNA hybrid duplex containing a 3'-CH(2)N(CH(3))-O-5' or an MMI backbone linker.

The solution structure of an antisense DNA.RNA hybrid duplex, d(CGCGTT-MMI-TTGCGC).r(GCGCAAAACGCG) (designated R4), containing an MMI backbone linker [3'-CH(2)N(CH(3))-O5'], is elucidated. The structural details of the MMI linker, its structural effects on the neighboring residues, and the molecular basis of the MMI effects are examined. The lipophilic N-methyl group of MMI is peripheral to the helix, assuming a conformation that is most stable with regard to the N-O torsion angle. The MMI linker promotes a 3'-endo conformation for the sugar moieties at both 3'- and 5'-adjacent positions and a backbone kink involving distant residues along the 3'-direction. Comparison of R4 with other analogous hybrid duplexes previously studied in this laboratory reveals a new family of low-energy helical conformations that can be accommodated in stable duplexes and a common feature of C3'-modified sugars for adopting a C3'-endo pucker. The results of these studies emphasize the interplay of several factors that govern the formation of stable hybrid duplexes and provide a basis for the understanding of the biological role of the MMI modifications, which are important building blocks for a family of promising chimeric antisense oligonucleotides.     

The solution structure of an oligonucleotide duplex containing a 2'-deoxyadenosine-3-(2-hydroxyethyl)- 2'-deoxyuridine base pair determined by NMR and molecular dynamics studies

Determination of the solution structure of the duplex d(GCAAGTC(HE)AAAACG)·d(CGTTTTAGACTTGC) containing a 3-(2-hydroxyethyl)-2′-deoxyuridine·deoxyadenine (HE·A) base pair is reported. The three-dimensional solution structure, determined starting from 512 models via restrained molecular mechanics using inter-proton distances and torsion angles, converged to two final families of structures. For both families the HE and the opposite A residues are intrahelical and in the anti conformation. The hydroxyethyl chain lies close to the helix axis and for one family the hydroxyl group is above the HE·A plane and in the other case it is below. These two models were used to start molecular dynamic calculations with explicit solvent to explore the hydrogen bonding possibilities of the HE·A base pair. The dynamics calculations converge finally to one model structure in which two hydrogen bonds are formed. The first is formed all the time and is between HEO4 and the amino group of A, and the second, an intermittent one, is between the hydroxyl group and the N1 of A. When this second hydrogen bond is not formed a weak interaction CH···N is possible between HEC7H2 and N1A21. All the best structures show an increase in the C1′–C1′ distance relative to a Watson–Crick base pair.     

A NMR strategy to unambiguously distinguish nucleic acid hairpin and duplex conformations applied to a Xist RNA A-repeat

All RNA sequences that fold into hairpins possess the intrinsic potential to form intermolecular duplexes because of their high self-complementarity. The thermodynamically more stable duplex conformation is favored under high salt conditions and at high RNA concentrations, posing a challenging problem for structural studies of small RNA hairpin conformations. We developed and applied a novel approach to unambiguously distinguish RNA hairpin and duplex conformations for the structural analysis of a Xist RNA A-repeat. Using a combination of a quantitative HNN-COSY experiment and an optimized double isotope-filtered NOESY experiment we could define the conformation of the 26-mer A-repeat RNA. In contrast to a previous secondary structure prediction of a double hairpin structure, the NMR data show that only the first predicted hairpin is formed, while the second predicted hairpin mediates dimerization of the A-repeat by duplex formation with a second A-repeat. The strategy employed here will be generally applicable to identify and quantify populations of hairpin and duplex conformations and to define RNA folding topology from inter- and intra-molecular base-pairing patterns.     

Probing structural factors stabilizing antisense oligonucleotide duplexes: NMR studies of a DNA.DNA duplex containing a formacetal linkage.

The duplex formed by annealing the formacetal backbone modified dodecamer d-(CGCGTTOCH2OTTGCGC) to its complementary strand, d(GCGCAAAACGCG) (duplex I), has been studied by NMR techniques and analyzed with reference to its unmodified counterpart (duplex II). Comparison of parameters such as 2D cross-peak intensities, coupling constants, and spectral patterns indicates that structural perturbations caused by the incorporation of the formacetal linkage are minimal and localized to the central T4.A4 block. Duplex I adopts a B-type helical conformation with regular Watson-Crick base pairing and normal minor groove width. The methylene group is accommodated along the phosphate backbone in a conformation similar to that of the PO2 group found in the B-form DNA family. The central T6-T7 base pairs of duplex I melt simultaneously with the duplex, indicating a cooperative transition to single strands. Although the formacetal linkage affects global melting, as evidenced by a 3 degree C reduction in Tm for duplex I with respect to duplex II, the present study indicates that this is not the result of localized premelting at the formacetal site of duplex I but rather reflects the subtle interplay of several structural and energy factors which need to be further explored.     

Preparation and characterization of a uniformly 2 H/ 15 N-labeled RNA oligonucleotide for NMR studies.

An RNA oligonucleotide that contains the binding site for Escherichia coli ribosomal protein S8 was prepared with uniform 15N isotopic enrichment and uniform deuterium enrichment at all non-exchangeable sites using enzymatic methods. The RNA binding site, which contains 44 nt, forms a hairpin in solution and requires Mg2+for proper folding. The longitudinal magnetization recovery rates of the exchangeable protons were compared for the [2H,15N]-enriched RNA molecule and for the corresponding fully [2H,15N]-enriched RNA hairpin. It was found that 1H-1H dipolar relaxation significantly contributes to the recovery of exchangeable proton longitudinal magnetization. The exchangeable proton resonance line widths were less affected by deuteration, indicating that chemical exchange with H2O remains the dominant mechanism of transverse magnetization relaxation. Nevertheless, deuteration of this RNA hairpin was found to enhance the sensitivity of NOE-based experiments relative to the fully protonated hairpin and to simplify 2D NMR spectra. The increased signal-to-noise ratio facilitated the assignment of the cytidine amino resonances and several of the purine nucleotide amino resonances and permitted the identification of NOE crosspeaks that could not be observed in spectra of the fully protonated RNA hairpin.     

Interactions of molecules with nucleic acids. I. An algorithm to generate nucleic acid structures with an application to the B-DNA structure and a counterclockwise helix.

An algorithm is developed that enables the routine determination of backbone conformations of nucleic acids. All atomic positions including hydrogen are specified in accord with experimental bond lengths and angles but with theoretically determined conformational angles. For two Watson-Crick base pairs at a separation of 3.38 Å, and perpendicular to a common helical axis, minimum energy configurations are found for all 10 combinations at helical angles of α ～ 36°–38°, corresponding to the B-DNA structure with C(2′)-endo sugar puckers. Backbone configurations exist only within the range 35.5° ? α ? 42°, which suggests the origin of the 10-fold helix. Calculated stacking energies for the B-DNA structure increases for each of the clustered groups of base pairs: G·C with G·C, G·C with A·T, and A·T with A·T, and they are in approximate agreement with experimental observations. The counter-clockwise helix is examined, and physically meaningful structures are found only when the helical axes of successive base pairs are disjointed.     

NMR studies of fully modified locked nucleic acid (LNA) hybrids: solution structure of an LNA:RNA hybrid and characterization of an LNA:DNA hybrid

LNA is a bicyclic nucleic acid analogue that contains one or more 2'-O,4'-C methylene linkage(s), which effectively locks the furanose ring in a C3'-endo conformation. We report here the NMR solution structure of a nonamer LNA:RNA hybrid and a structural characterization of a nonamer LNA:DNA hybrid, where the LNA strands are composed entirely of LNA nucleotides. This is the first structural characterization of fully modified LNA oligonucleotides. The high-resolution structure reveals that the LNA:RNA hybrid adopts an almost canonical A-type duplex morphology. The helix axis is almost straight and the duplex geometry is regular. This shows that fully modified LNA oligomers can hybridize with complementary RNA and form duplexes within the Watson-Crick framework. The LNA:DNA hybrid structurally resembles an RNA:DNA hybrid as shown by determination of deoxyribose sugar puckers and analysis of NOESY NMR spectra.     

Biophysical studies of a ruthenium(II) polypyridyl complex binding to DNA and RNA prove that nucleic acid structure has significant effects on binding behaviors

The interactions of a metal complex [Ru(phen)(2)PMIP](2+) {Ru=ruthenium, phen=1,10-phenanthroline, PMIP=2-(4-methylphenyl)imidazo[4,5-f]1,10-phenanthroline} with yeast tRNA and calf thymus DNA (CT DNA) have been investigated comparatively by UV-vis spectroscopy, fluorescence spectroscopy, viscosity measurements, isothermal titration calorimetry (ITC), as well as equilibrium dialysis and circular dichroism (CD). Spectroscopic studies together with ITC and viscosity measurements indicate that both binding modes of the Ru(II) polypyridyl complex to yeast tRNA and CT DNA are intercalation and yeast tRNA binding of the complex is stronger than CT DNA binding. ITC experiments show that the interaction of the complex with yeast tRNA is driven by a moderately favorable enthalpy decrease in combination with a moderately favorable entropy increase, while the binding of the complex to CT DNA is driven by a large favorable enthalpy decrease with a less favorable entropy increase. The results from equilibrium dialysis and CD suggest that both interactions are enantioselective and the Delta enantiomer of the complex may bind more favorably to both yeast tRNA and CT DNA than the Lambda enantiomer does, and that the complex is a better candidate for an enantioselective binder to yeast tRNA than to CT DNA. Taken together, these results indicate that the structures of nucleic acids have significant effects on the binding behaviors of metal complexes.     

Solution structures of a duplex containing an adenine opposite a gap (absence of one nucleotide). An NMR study and molecular dynamic simulations with explicit water molecules.

We investigated the behaviour of a 15mer DNA duplex, [5'd(CAGAGTCACTGGCTC)3']. [5'd(GAGCCAG)3' + 5'd(GACTCTG)3'] which contained an adenine opposite the gap. Analysis of the NMR data showed the existence of one major species, which was in equilibrium with two minor species. Their relative concentrations varied as a function of pH with a pKa of approximately 4.5. For the major species, the duplex was globally in B conformation with the central adenine stacked in the helix. The two G.C base pairs adjacent to the central adenine were well formed and a gap was present in front of this adenine. For the minor species, major structural perturbations occurred in the centre of the duplex. At neutral pH, the central adenine was involved in a G.A mismatch with G23 adjacent to the gap. Cytosine C7 was then extrahelical and no gap was observed. Under these conditions, the major neutral species corresponded to 70% of the total and the minor species to 30%. At acidic pH, the central adenine of the minor species was protonated and was involved in a G(syn).A+(anti) mismatch. The difference is that C9 is now extrahelical and G22 is implicated in the mispair. Three-dimensional models were built to initiate molecular dynamic simulations, which were in good agreement with the NMR data. Their structural stability in terms of hydrogen bonding and their flexibility are discussed and the biological significance for the interaction with DNA polymerase is evoked.     

Solution structure of an LNA hybridized to DNA: NMR study of the d(CT(L)GCT(L)T(L)CT(L)GC):d(GCAGAAGCAG) duplex containing four locked nucleotides

We have used two-dimensional (1)H NMR spectroscopy at 750 MHz to determine a high-resolution solution structure of an oligonucleotide containing restricted nucleotides with a 2'-O, 4'-C-methylene bridge (LNA) hybridized to the complementary DNA strand. The LNA:DNA duplex examined contained four thymidine LNA modifications (T(L), d(C1T(L)2G3C4T(L)5T(L)6C7T(L)8G9C10):d( G11C12A13G14A15A16G17C 18A19G20). A total relaxation matrix approach was used to obtain interproton distance bounds from NOESY cross-peak intensities. These distance bounds were used as restraints in molecular dynamics (rMD) calculations. Forty final structures were generated for the duplex from A-form and B-form DNA starting structures. The root-mean-square deviation (RMSD) of the coordinates for the 40 structures of the complex was 0.6 A. The sugar puckerings are averaged values of a dynamic interchange between N- and S-type conformation except in case of the locked nucleotides that were found to be fixed in the C3'-endo conformation. Among the other nucleotides in the modified strand, the furanose ring of C7 and G9 is predominantly in the N-type conformation whereas that of G3 is in a mixed conformation. The furanose rings of the nucleotides in the unmodified complementary strand are almost exclusively in the S-type conformation. Due to these different conformations of the sugars in the two strands, there is a structural strain between the A-type modified strand and the B-type unmodified complementary strand. This strain is relaxed by decreasing the value of rise and compensating with tip, buckle, and propeller twist. The values of twist vary along the strand but for a majority of the base pairs a value even lower than that of A-DNA is observed. The average twist over the sequence is 32+/-1 degrees. On the basis of the structure, we conclude that the high stability of LNA:DNA duplexes is caused by a local change of the phosphate backbone geometry that favors a higher degree of stacking.     

3D structure of bovine pancreatic ribonuclease A in aqueous solution: An approach to tertiary structure determination from a small basis of1H NMR NOE correlations

Summary A method is proposed to generate initial structures in cases where the distance geometry method may fail, such as when the set of¹H NMR NOE-based distance constraints is small in relation to the size of the protein. The method introduces an initial correlation between the and backbone angles (based on empirical observations) which is relaxed in later stages of the calculation. The obtained initial structures are refined by well-established methods of energy minimization and restrained molecular dynamics. The method is applied to determine the solution structure of Ribonuclease A (124 residues) from a NOE basis consisting of 467 NOE cross-correlations (97 intra-residue, 206 sequential, 23 medium-range and 141 long-range) obtained at 360 MHz. The global shape and backbone overall fold of the eight final refined structures are close to those shown by the crystal structure. A meaningful difference in the positioning of the catalytically important His¹¹⁹ side chain in the solution and crystal structures has been detected.     

Solution structure of an arabinonucleic acid (ANA)/RNA duplex in a chimeric hairpin: comparison with 2′-fluoro-ANA/RNA and DNA/RNA hybrids

Hybrids of RNA and arabinonucleic acid (ANA) as well as the 2′-fluoro-ANA analog (2′F-ANA) were recently shown to be substrates of the enzyme RNase H. Although RNase H binds to double-stranded RNA, no cleavage occurs with such duplexes. Therefore, knowledge of the structure of ANA/RNA hybrids may prove helpful in the design of future antisense oligonucleotide analogs. In this study, we have determined the NMR solution structures of ANA/RNA and DNA/RNA hairpin duplexes and compared them to the recently published structure of a 2′F-ANA/RNA hairpin duplex. We demonstrate here that the sugars of RNA nucleotides of the ANA/RNA hairpin stem adopt the C3′-endo (north, A-form) conformation, whereas those of the ANA strand adopt a ‘rigid’ O4′-endo (east) sugar pucker. The DNA strand of the DNA/RNA hairpin stem is flexible, but the average DNA/RNA hairpin structural parameters are close to the ANA/RNA and 2′F-ANA/RNA hairpin parameters. The minor groove width of ANA/RNA, 2′F-ANA/RNA and DNA/RNA helices is 9.0 ± 0.5 Å, a value that is intermediate between that of A- and B-form duplexes. These results rationalize the ability of ANA/RNA and 2′F-ANA/RNA hybrids to elicit RNase H activity.     

Molecular recognition between oligopeptides and nucleic acids. Specificity of binding of a monocationic bis-furan lexitropsin to DNA deduced from footprinting and 1H NMR studies

MPE-Fe(EDTA) footprinting of a novel monocationic bis-furan lexitropsin 6 on a HindIII/EcoRI restriction fragment of pBR322 DNA revealed a series of four-base binding sites (all 5'----3') of (primary) TGTA, TGAA, AAAT, ACAA, TTAT, and (secondary) CTAA, TCGT, TGTA, GTCA, and GGTT. Thus 6 can accept a GC pair at positions 1, 2 or 3 of the binding site with a strict 3' (4 position) AT requirement. Marked enhancement of cleavage, particularly at GC rich sequences, is observed at regions flanking or even up to 18 base pairs remote from a given binding site. The non-exchangeable and imino 1H NMR resonances of the 1:1 complex and d-[CATGGCCATG]2 were assigned using a combination of NOE differences, NOESY and COSY techniques. 1H NMR studies (ligand induced chemical shifts and NOE differences) of Lexitropsin 6 with d-[CATGGCCATG]2 show unambiguously the location and orientation of the N to C termini of 6 on the sequence 5'-G5C6C7A8-3', with the C terminus oriented to A8. This orientation of 6 in the minor groove of 5'-GCCA is confirmed by an NOE observed between H1 2a of 6 and AH8(8). This preference for binding of 6 to the sequence 5'-GCCA when challenged with d-[CATGGCCATG]2 is in accord with the conclusions of the footprinting experiments wherein GC base pairs can be accepted in the first three positions and with a strict 3' terminus AT reading requirement. Collectively the data support the inference of a GC recognizing capacity for a 2,5-substituted furan moiety within a lexitropsin. The 1H NMR data indicate that the decadeoxyribonucleotide duplex exists in the B conformation in both the 1:1 complex and the free form. The apparent binding constant of 6 to calf thymus DNA is 1.68 X 10(5) M-1 whereas netropsin under similar conditions gives a value of 1.85 X 10(7) M-1. This suggests that if advantage is to be taken of the GC recognizing property of a 2,5-substituted furan in longer lexitropsins it should be flanked by more strongly bound moieties.     

Shuffling of amino acid sequence: an important control in synthetic peptide studies of nucleic acid-binding domains. Binding properties of fragments of a conserved eukaryotic RNA binding motif.

We have used synthetic peptides to study a conserved RNA binding motif in yeast poly(A)-binding protein. Two peptides, 45 and 44 amino acids in length, corresponding to amino and carboxyl halves of a 90-amino acid RNA-binding domain in the protein were synthesized. While the amino-terminal peptide had no significant affinity for nucleic acids, the carboxyl-terminal peptide-bound nucleic acids with similar characteristics to that for the entire 577 residue yeast poly(A)-binding protein. In 100 mM NaCl, the latter peptide retained over 50% of the intrinsic binding free energy of the protein, as well as, similar RNA versus DNA binding specificity. However, shuffling of the sequence of this 44 residue peptide had surprisingly little effect on its nucleic acid binding properties suggesting the overriding importance of amino acid composition as opposed to primary sequence. Deletion studies on the 44 residue peptide with the "correct" sequence succeeded in identifying amino acids important for conferring RNA specificity and for increasing our understanding of the molecular basis for nucleic acid binding by synthetic peptides. The shuffled peptide study, however, clearly indicates that considerable caution must be exercised before extrapolating results of structure/function studies on synthetic peptide analogues to the parent protein.     

Incorporation of a complete set of deoxyadenosine and thymidine analogues suitable for the study of protein nucleic acid interactions into oligodeoxynucleotides. Application to the EcoRV restriction endonuclease and modification methylase

A complete set of dA and T analogues designed for the study of protein DNA interactions has been prepared. These modified bases have been designed by considering the groups on the dA and T bases that are accessible to proteins when these bases are incorporated into double-helical B-DNA [Seeman, N. C., Rosenberg, J. M., & Rich, A. (1976) Proc. Natl. Acad. Sci. U.S.A. 73, 804-808]. Each of the positions on the two bases, having the potential to interact with proteins, have been subject to nondisruptive, conservative change. Typically a particular group (e.g., the 6-NH2 of dA or the 5-CH3 of T) has been replaced with a hydrogen atom. Occasionally keto groups (the 2- and 4-keto oxygen atoms of T) have been replaced with sulfur. The base set has been incorporated into the self-complementary dodecamer d(GACGATATCGTC) at the central d(ATAT) sequence. Melting temperature determination shows that the modified bases do not destabilize the double helix. Additionally, circular dichroism spectroscopy shows that almost all the altered bases have very little effect on overall oligodeoxynucleotide conformation and that most of the modified oligomers have a B-DNA type structure. d(GATATC) is the recognition sequence for the EcoRV restriction modification system. Initial rate measurements (at a single oligodeoxynucleotide concentration of 20 microM) have been carried out with both the EcoRV restriction endonuclease and modification methylase. This has enabled a preliminary identification of the groups of the dA and T bases within the d(GATATC) sequence that make important contacts to both proteins.     

Crystal and molecular structure of a DNA fragment containing a 2-aminoadenine modification: the relationship between conformation, packing, and hydration in Z-DNA hexamers.

The crystal and molecular structure of d(CGUA'CG)2 (where A' is 2-aminoadenine) has been determined and refined to an R factor of 13.8% for data 8.0-1.3 A. The structure is very similar to the original Z-DNA structures with the sequence d(CGCGCG)2 [Gessner, R. V., Frederick, C. A., Quigley, G. J., Rich, A., & Wang, A. H.-J. (1989) J. Biol. Chem. 264, 7921] and shows that the substitution of 2-aminoadenine-uracil base pairs in the two central steps is consistent with Z-DNA formation. In addition, we show how waters mediating intermolecular interactions may help to explain the ZI-ZII conformational pattern found in many Z-DNA structures.