Structural basis of cleavage by RNase H of hybrids of arabinonucleic acids and RNA

The origins of the substrate specificity of Escherichia coli RNase H1 (termed RNase H here), an enzyme that hydrolyzes the RNA strand of DNA-RNA hybrids, are not understood at present. Although the enzyme binds double-stranded RNA, no cleavage occurs with such duplexes [Lima, W. F., and Crooke, S. T. (1997) Biochemistry 36, 390]. Therefore, the hybrid substrates may not adopt a canonical A-form geometry. Furthermore, RNase H is exquisitely sensitive to chemical modification of the DNA strands in hybrid duplexes. This is particularly relevant to the RNase H-dependent pathway of antisense action. Thus, only very few of the modifications currently being evaluated as antisense therapeutics are tolerated by the enzyme, among them phosphorothioate DNA (PS-DNA). Recently, hybrids of RNA and arabinonucleic acid (ANA) as well as the 2'F-ANA analogue were shown to be substrates of RNase H [Damha, M. J., et al. (1998) J. Am. Chem. Soc. 120, 12976]. Using X-ray crystallography, we demonstrate here that ANA analogues, such as 2'F-ANA [Berger, I., et al. (1998) Nucleic Acids Res. 26, 2473] and [3.3.0]bicyclo-ANA (bc-ANA), may not be able to adopt sugar puckers that are compatible with pure A- or a B-form duplex geometries, but rather prefer the intermediate O4'-endo conformation. On the basis of the observed conformations of these ANA analogues in a DNA dodecamer duplex, we have modeled a duplex of an all-C3'-endo RNA strand and an all-O4'-endo 2'F-ANA strand. This duplex exhibits a minor groove width that is intermediate between that of A-form RNA and B-form DNA, a feature that may be exploited by the enzyme in differentiating between RNA duplexes and DNA-RNA hybrids. Therefore, the combination of the established structural and functional properties of ANA analogues helps settle existing controversies concerning the discrimination of substrates by RNase H. Knowlegde of the structure of an analogue that exhibits enhanced RNA affinity while not interfering with RNase H activity may prove helpful in the design of future antisense modifications.     

The structure of poly(dA).poly(dT) as revealed by an X-ray fibre diffraction

X-ray diffraction in fibres revealed that the calcium salt of poly(dA).poly(dT) is a 10-fold double helix with a pitch of 3.23 nm. The opposite sugar-phosphate chains in the refined model are characterized by a complete conformational equivalence and contain sugars in a conformation close to C2'-endo. As a result a new model of the sodium salt of poly(dA).poly(dT) has been constructed, which is different from the Heteronomous DNA proposed earlier (S. Arnott et al., Nucl. Acids Res. 11, 4141 (1983)). The new model of Na-poly(dA).poly(dT) has conformationally similar opposite chains; it is a structure of the B-type, rather like that of Ca-poly(dA).poly(dT).     

1H two-dimensional nuclear Overhauser effect and relaxation studies of poly(dA).poly(dT)

The structure of poly(dA).poly(dT) in aqueous solution has been studied by using 1H two-dimensional nuclear Overhauser effect (2D NOE) spectroscopy and relaxation rate measurements on the imino and nonexchangeable protons. The assignments of the 1H resonances are determined from the observed cross-relaxation patterns in the 2D NOE experiments. The cross-peak intensities together with the measured relaxation rates show that the purine and pyrimidine strands in poly(dA).poly(dT) are equivalent in aqueous solution. The results are consistent with a right-handed B-form helix where the sugars on both strands are in the C2'-endo/anti configuration. These observations are inconsistent with a proposed heteronomous structure for poly(dA).poly(dT) [Arnott, S., Chandrasekaran, R., Hall, I. H., & Puigjaner, L. C. (1983) Nucleic Acids Res. 11, 4141-4155]. The measured relaxation rates also show that poly(dA).poly(dT) has fast, large-amplitude local internal motions (+/- 20-25 degrees) in solution and that the amplitudes of the base and sugar motions are similar. The motion of the bases in poly(dA).poly(dT) is also similar to that previously reported for poly(dA-dT).poly(dA-dT) and poly(dG-dC).poly(dG-dC) [Assa-Munt, N., Granot, J., Behling, R. W., & Kearns, D. R. (1984) Biochemistry 23, 944-955; Mirau, P. A., Behling, R. W., & Kearns, D. R. (1985) Biochemistry 24, 6200-6211].     

Complete assignment of the deoxyribose 5'/5" proton resonances of the EcoRI DNA sequence using isotropic mixing

Using two-dimensional isotropic mixing spectroscopy all 5'/5" proton resonances of the EcoRI restriction site DNA dodecamer [d(CGCGAATTCGCG)]2 have been assigned. This completes the previous assignments of 1'H to 4'H resonances of the deoxyribose spin systems (Hare et al., 1983). With mixing times of up to 500 ms, many of these resonances showed connectivities of 5'/5" protons in the two-dimensional isotropic mixing spectrum. Relying only on through-bond connectivities makes these assignments independent of assumptions about the conformation of the DNA oligonucleotide. The assignment of the 5'H/5"H resonances will allow the interpretation of intra- and interresidue NOEs to these protons, providing information about the DNA backbone conformation.     

500-MHz 1H NMR study of poly(dG).poly(dC) in solution using one-dimensional nuclear Overhauser effect

Secondary structures of poly(dG).poly(dC) and poly(dG).poly(dm5C) in solution are determined by nuclear Overhauser effect (NOE) measurements on GH8-deuterated and -nondeuterated DNAs with low presaturation pulse lengths (10-25 ms) and low-power and prolonged accumulations in the range of 50,000-72,000 scans. Under these conditions, the NOE difference spectra were free from diffusion. Primary NOEs between base protons GH8/CH6 and sugar protons H1', H2'/H2', and H3' suggest that in poly(dG).poly(dC) both guanine and cytosine nucleotides adopt a C3'-endo, low anti X = 200-220 degrees conformation. Computer modeling of the NOE data enable identification for the first time, in terms of the geometry of the nucleotide repeat, handedness, and helix geometry, of the structure of poly(dG).poly(dC) to be the A form, and the derived structure for the polymer duplex is very close to the single crystal structure of the double-helical d-GGGGCCCC [McCall, M., Brown, T., & Kennard, O. (1985) J. Mol. Biol. 183, 385-396]. Similar nuclear Overhauser effect data on poly(dG).poly(dm5C) revealed that G and m5C adopt a C2'endo, anti X = 240-260 degrees conformation, which indicates that this DNA exhibits the B form in solution. In summary, the results presented in this paper demonstrate that methylation of cytosines in poly(dG).poly(dC) causes A----B transition in the molecule.     

Duplex recognition by oligonucleotides containing 2'-deoxy-2'-fluoro-D-arabinose and 2'-deoxy-2'-fluoro-D-ribose. Intermolecular 2'-OH-phosphate contacts versus sugar puckering in the stabilization of triple-helical complexes

To gain insight into the origins of the large binding affinity of RNA toward target duplexes, 2'-deoxy-2'-fluororibonucleic acid (2'F-RNA) and 2'-deoxy-2'-fluoroarabinonucleic acid (2'F-ANA) were tested for their ability to recognize duplex DNA, duplex RNA, and RNA-DNA hybrids. 2'F-RNA, 2'F-ANA, and the corresponding control single-stranded (ss) DNA strands were shown to form triple-helical complexes only with duplex DNA and hybrid DNA (Pu)-RNA (Py), but not with duplex RNA and hybrid RNA (Pu)-DNA (Py). In contrast, an RNA third strand recognized all four possible duplexes (DD, DR, RD, and RR) as previously demonstrated by Roberts and Crothers [(1992) Science 258, 1463-1466]. The 2'F-RNA (C3'-endo) strand exhibited significantly reduced affinity for duplexes compared to an unmodified RNA (C3'-endo) strand. These findings are consistent with the intermolecular 2'-OH-phosphate contact mechanism proposed by Escudé et al. [(1993) Nucleic Acids Res. 24, 5547-5553], as a ribo 2'-F atom should not interact with a negatively charged phosphate. In addition, they emphasize the role of the 2'-OH ribose as a general recognition and binding determinant of RNA. The 2'-F arabino modification (2'F-ANA, C2'-endo) led to a considerable increase in the binding affinity for duplex DNA, as compared to those of DNA and 2'F-RNA third strands. This is likely to be the result of a greater population of C2'-endo pucker of the 2'F-ANA compared to DNA. The enhancement observed for 2'F-ANA strands toward duplex DNA is comparable to that observed with 2'-OMe RNA. Since 2'F-ANA has been shown to be more resistant to nuclease degradation than DNA, these results are likely to stimulate experimental work on arabinose derivatives in laboratories concerned with targeting DNA sequences in vivo ("antigene" strategy).     

Hydration of [d(CGC)r(aaa)d(TTTGCG)](2)

We have studied the hydration and dynamics of RNA C2'-OH in a DNA. RNA hybrid chimeric duplex [d(CGC)r(aaa)d(TTTGCG)](2). Long-lived water molecules with correlation time tau(c) larger than 0.3 ns were found close to the RNA adenine H2 and H1' protons in the hybrid segment. A possible long-lived water molecule was also detected close to the methyl group of 7T in the RNA-DNA junction but not to the other two thymine bases (8T and 9T). This result correlates with the structural studies that only DNA residue 7T in the RNA-DNA junction adopts an O4'-endo sugar conformation (intermediate between B-form and A-form), while the other DNA residues including 3C in the DNA-RNA junction, adopt C1'-exo or C2'-endo conformations (in the B-form domain). Based on the NOE cross-peak patterns, we have found that RNA C2'-OH tends to orient toward the O3' direction, forming a possible hydrogen bond with the 3'-phosphate group. The exchange rates for RNA C2'-OH were found to be around 5-20 s(-1), compared to 26.7(+/-13.8) s(-1) reported previously for the other DNA.RNA hybrid duplex. This slow exchange rate may be due to the narrow minor groove width of [d(CGC)r(aaa)d(TTTGCG)](2), which may trap the water molecules and restrict the dynamic motion of hydroxyl protons. The distinct hydration patterns of the RNA adenine H2 and H1' protons and the DNA 7T methyl group in the hybrid segment, as well as the orientation and dynamics of the RNA C2'-OH protons, may provide a molecular basis for further understanding the structure and recognition of DNA.RNA hybrid and chimeric duplexes.     

Solution structure of a trans-opened (10S)-dA adduct of (+)-(7S,8R,9S,10R)-7,8-dihydroxy-9,10-epoxy-7,8,9,10-tetrahydrobenzo[a]pyrene in a fully complementary DNA duplex: evidence for a major syn conformation

Two-dimensional NMR was used to determine the solution structure of an undecanucleotide duplex, d(CGGTCACGAGG).d(CCTCGTGACCG), in which (+)-(7S,8R,9S,10R)-7,8-dihydroxy-9,10-epoxy-7,8,9,10-tetrahydrobenzo[a]pyrene is covalently bonded to the exocyclic N(6)() amino group of the central deoxyadenosine, dA(6), through trans addition at C10 of the epoxide (to give a 10S adduct). The present study represents the first NMR structure of a benzo[a]pyrene (10S)-dA adduct in DNA with a complementary T opposite the modified dA. Exchangeable and nonexchangeable protons of the modified duplex were assigned by the use of TOCSY (in D(2)O) and NOESY spectra (in H(2)O and D(2)O). Sequential NOEs expected for a B-type DNA conformation with typical Watson-Crick base pairing are observed along the duplex, except at the lesion site. We observed a strong intraresidue NOE cross-peak between H1' and H8 of the modified dA(6). The sugar H2' and H2' ' of dC(5) lacked NOE cross-peaks with H8 of dA(6) but showed weak interactions with H2 of dA(6) instead. In addition, the chemical shift of the H8 proton (7.51 ppm) of dA(6) appears at a higher field than that of H2 (8.48 ppm). These NOE and chemical shift data for the dA(6) base protons are typical of a syn glycosidic bond at the modified base. Restrained molecular dynamics/energy minimization calculations show that the hydrocarbon is intercalated from the major groove on the 3'-side of the modified base between base pairs A(6)-T(17) and C(7)-G(16) and confirm the syn glycosidic angle (58 degrees ) of the modified dA(6). In the syn structure, a weak A-T hydrogen bond is possible between the N3-H proton of T(17) and N7 of dA(6) (at a distance of 3.11 A), whereas N1, the usual hydrogen bonding partner for N3-H of T when dA is in the anti conformation, is 6.31 A away from this proton. The 10(S)-dA modified DNA duplex remains in a right-handed helix, which bends in the direction of the aliphatic ring of BaP at about 42 degrees from the helical axis. ROESY experiments provided evidence for interconversion between the major, syn conformer and a minor, possibly anti, conformer.     

The propeller DNA conformation of poly(dA).poly(dT).

Physical properties of the DNA duplex, poly(dA).poly(dT) differ considerably from the alternating copolymer poly(dAT). A number of molecular models have been used to describe these structures obtained from fiber X-ray diffraction data. The recent solutions of single crystal DNA dodecamer structures with segments of oligo-A.oligo-T have revealed the presence of a high propeller twist in the AT regions which is stabilized by the formation of bifurcated (three-center) hydrogen bonds on the floor of the major groove, involving the N6 amino group of adenine hydrogen bonding to two O4 atoms of adjacent thymine residues on the opposite strand. Here we show that it is possible to incorporate the features of the single crystal analysis, specifically high propeller twist, bifurcated hydrogen bonds, and a narrow minor groove, as well as the close interstrand NMR signal between adenine HC2 and ribose HC1' of the opposite strand, into a model that is fully compatible with the diffraction data obtained from poly(dA).poly(dT).     

A two-dimensional NMR study of the solution structure of a DNA dodecamer comprising the concensus sequence for the specific DNA-binding sites of the glucocorticoid receptor protein

A two-dimensional 500-MHz 1H-NMR study on the non-self-complementary double-stranded DNA dodecamer 5'd(C-C-A-G-A-A-C-A-G-T-G-G)5'd(C-C-A-C-T-G-T-T-C-T-G-G), is presented. This oligonucleotide contains the consensus octanucleotide sequence 5'd(A-G-A-A-C-A-G-T) for the specific DNA-binding sites of the glucocorticoid receptor protein [Payvar, F. et al. (1984) Cell 35, 381-392]. Using a combination of two-dimensional pure phase absorption nuclear Overhauser enhancement (NOESY) and homonuclear J-correlated (COSY) spectroscopy all non-exchangeable base (with the exception fo the adenine H2 protons), methyl and deoxyribose H1', H2', H2", H3' and H4' resonances are assigned unambiguously using a sequential resonance assignment strategy. From the relative intensities of the cross peaks in the pure phase absorption NOESY spectra at two mixing times it is shown that the dodecamer adopts a B-type conformation in solution.     

Deoxyguanosine-deoxyadenosine pairing in the d(C-G-A-G-A-A-T-T-C-G-C-G) duplex: conformation and dynamics at and adjacent to the dG X dA mismatch site

Nuclear magnetic resonance (NMR) has been used to monitor the conformation and dynamics of the d-(C1-G2-A3-G4-A5-A6-T6-T5-C4-G3-C2-G1) self-complementary dodecanucleotide (henceforth called 12-mer GA) that contains a dG X dA purine-purine mismatch at position 3 in the sequence. These results are compared with the corresponding d(C-G-C-G-A-A-T-T-C-G-C-G) dodecamer duplex (henceforth called 12-mer) containing standard Watson-Crick base pairs at position 3 [Patel, D.J., Kozlowski, S.A., Marky, L.A., Broka, C., Rice, J.A., Itakura, K., & Breslauer, K.J. (1982) Biochemistry 21, 428-436]. The dG X dA interaction at position 3 was monitored at the guanosine exchangeable H-1 and nonexchangeable H-8 protons and the nonexchangeable adenosine H-2 proton. We demonstrate base-pair formation between anti orientations of the guanosine and adenosine rings on the basis of nuclear Overhauser effects (NOE) observed between the H-2 proton of adenosine 3 and the imino protons of guanosine 3 (intra base pair) and guanosines 2 and 4 (inter base pair). The dG(anti) X dA(anti) pairing should result in hydrogen-bond formation between the guanosine imino H-1 and carbonyl O-6 groups and the adenosine N-1 and NH2-6 groups, respectively. The base pairing on either side of the dG X dA pair remains intact at low temperature, but these dG X dC pairs at positions 2 and 4 are kinetically destabilized in the 12-mer GA compared to the 12-mer duplex. We have estimated the hydrogen exchange kinetics at positions 4-6 from saturation-recovery measurements on the imino protons of the 12-mer GA duplex between 5 and 40 degrees C. The measured activation energies for imino proton exchange in the 12-mer GA are larger by a factor of approximately 2 compared to the corresponding values in the 12-mer duplex. This implies that hydrogen exchange in the 12-mer GA duplex results from a cooperative transition involving exchange of several base pairs as was previously reported for the 12-mer containing a G X T wobble pair at position 3 [Pardi, A., Morden, K.M., Patel, D.J., & Tinoco, I., Jr. (1982) Biochemistry 21, 6567-6574]. We have assigned the nonexchangeable base protons by intra and inter base pair NOE experiments and monitored these assigned markers through the 12-mer GA duplex to strand transition.(ABSTRACT TRUNCATED AT 400 WORDS)     

Structure of poly(dA).poly(dT) is not identical to the AT rich regions of the single crystal structure of CGCGAATTBrCGCG. The consequence of this to netropsin binding to poly(dA).poly(dT)

The basic assumption of Dickerson and Kopka (J. Biomole. Str. Dyns. 2, 423, 1985) that the conformation of poly(dA).poly(dT) in solution is identical to the AT rich region of the single crystal structure of the Dickerson dodecamer is not supported by any experimental data. In poly(dA).poly(dT), NOE and Raman studies indicate that the dA and dT units are conformationally equivalent and display the (anti-S-type sugar)-conformation; incorporation of this nucleotide geometry into a double helix leads to a conventional regular B-helix in which the width of the minor groove is 8A. The derived structure is consistent with all available experimental data on poly(dA).poly(dT) obtained under solution conditions. In the crystal structure of the dodecamer, the dA and dT units have distinctly different conformations-dA residues adopt (anti, S-type sugar pucker), while dT residues belong to (low anti, N-type sugar pucker). These different conformations of the dA and dT units along with the large propeller twist can be accommodated in a double helix in which the minor groove is shrunk from 8A to less than 4A. In the conventional right handed B-form of poly(dA).poly(dT) with the 8A wide minor groove, netropsin has to bind asymmetrically along the dA strand to account for the NOE and chemical shift data and to generate a stereochemically sound structure (Sarma et al, J. Biomole. Str. Dyns. 2, 1085, 1985).     

Untenability of the Heteronomous DNA Model for Poly(dA) · Poly(dT) in Solution. This DNA Adopts a Right-Handed B-DNA Duplex in Which the Two Strands are Conformationally Equivalent. A 500 MHz NMR Study Using One Dimensional NOE

Abstract

ID NOE ¹H NMR spectroscopy at 500 MHz was employed to examine the structure of poly(dA)·poly(dT) in solution. NOE experiments were conducted as a function of presaturation pulse length (50, 30, 20 and 10 msec) and.power (19 and 20 db) to distinguish the primary NOEs from spin diffusion. The 10 msec NOE experiments took 49 hrs and over 55,000 scans for each case and the difference spectra were almost free from diffusion.

The spin diffused NOE difference spectra as well as difference NOE spectra in 90% H₂O + 10% D₂O in which TNH3 was presaturated enabled to make a complete assignment of the base and sugar protons. It is shown that poly(dA) ·poly(dT) melts in a fashion in which single stranded bubbles are formed with increasing temperature.

Extremely strong primary NOEs were observed at H2′/H2″ when AH8 and TH6 were presaturated. The observed NOEs at AH2′ and that AH2″ were very similar as were the NOEs at TH2′ and TH2″. The observed NOEs at AH2′ and AH2″when AH8 was presaturated were very similar to those observed at TH2′ and TH2″ when TH6 was presaturated. In addition, presaturation of H1′ of A and T residues resulted in similar NOEs at AH2′/H2″ and TH2′/H2″ region and these NOEs at H2′ and H2″ were distinctly asymmetric as expected in a C2′-endo sugar pucker. There was not a trace of NOE at AH8 and TH6 when AH3′ and TH3′ were presaturated indicating that C3′-endo, × = 30–40° conformation is not valid for this DNA. From these NOE data, chemical shift shielding calculations and stereochemistry based computer modellings, we conclude that poly(dA)·poly(dT) in solution adopts a right- handed B-DNA duplex in which both dA and dT strands are conformationally equivalent with C2′-endo sugar pucker and a glycosyl torsion, ×, of ?73°, the remaining backbone torsion angles being φ′ = 221°, ω′ = 212°, ω = 310°, φ = 149°, ψ = 42°, ψ′ = 139°. The experimental data are in total disagreement with the heteronomous DNA model of Arnott et. al. proposed for the fibrous state. (Arnott, S., Chandrasekaran, R., Hall, I.H., and Puigjaner, L.C., Nucl. Acid Res. 11, 4141, 1983).     

Crystal structure of an adenine bulge in the RNA chain of a DNA.RNA hybrid, d(CTCCTCTTC).r(gaagagagag)

Crystal structure of a DNA.RNA hybrid, d(CTCCTCTTC).r(gaagagagag), with an adenine bulge in the polypurine RNA strand was determined at 2.3 A resolution. The structure was solved by the molecular replacement method and refined to a final R-factor of 19.9% (Rfree 22.2%). The hybrid duplex crystallized in the space group I222 with unit cell dimensions, a = 46.66 A, b = 47.61 A and c = 54.05 A, and adopts the A-form conformation. All RNA and DNA sugars are in the C3'-endo conformation, the glycosyl angles in anti conformation and the majority of the C4'-C5' torsion angles in g+ except two trans angles, in conformity with the C3'-endo rigid nucleotide hypothesis. The adenine bulge is looped out and it is also in the anti C3'-endo conformation. The bulge is involved in a base-triple (C.g)*a interaction with the end base-pair (C9.g10) in the minor groove of a symmetry-related molecule. The 2' hydroxyl group of g15 is hydrogen bonded to O2P and O5' of g17, skipping the bulged adenine a16 and stabilizing the sugar-phosphate backbone of the hybrid. The hydrogen bonding and the backbone conformation at the bulged adenine site is very similar to that found in the crystal structure of a protein-RNA complex.     

Two-dimensional NMR studies on the anthramycin-d(ATGCAT)2 adduct

Two-dimensional NMR experiments were performed on the adduct of anthramycin with d(ATGCAT)2 to obtain the assignments of the nucleotide base and sugar protons as well as the anthramycin protons. Anthramycin is covalently attached to a guanine 2-amino group, forming the d(ATamGCAT).d(ATGCAT) modified duplex. The anthramycin protons in the minor groove exhibit NOEs to several nucleotide protons. The network of anthramycin-nucleotide NOEs and the measurement of the 10-Hz coupling constant between the anthramycin H11 and H11a protons shows that anthramycin is covalently attached as the S stereoisomer at the anthramycin C11 position with the side chain of anthramycin oriented toward the 5' end of the modified strand. The NOE data show that the anthramycin-modified duplex is in a right-handed conformation with all bases in an anti conformation. Analysis of the J1'-2' coupling constants for the resolved H1' resonances shows that the S-type conformation of the sugars is highly preferred.     

Sequence specific conformation of a DNA decamer containing an adenine tract studied in solution by H-NMR spectroscopy

The decanucleotide duplex d(AAAACGTTTT)2 and a variety of phase-sensitive two-dimensional (2D) NMR experiments have been used to investigate the solution conformation of an adenine-tract and its junction with another DNA sequence. 2D nuclear Overhauser effect data confirm that the oligonucleotide has a general B-type DNA morphology but an array of unusual correlations implies that the adenine tract and the 5'-ApC junction have conformations more compatible with the modified X-ray structures recently reported for DNAs of similar sequence (Nelson, H.C.M., Finch, J.T., Luisi, B.F. and Klug, A. (1987) Nature 330, 221-226). The pattern and magnitude of interstrand NOEs from the adenine H2s to the sugar H1's of the complementary base to the 5'-neighbouring residue indicate that the A-T basepairs are highly propeller twisted and that the minor groove is narrowed, showing its greatest compression at the 3'-end of the tract at the 5'-ApC step. Quantifying spin-coupling interactions within the deoxyribose rings by analysing both 1D and high-resolution 2D DQF-COSY data reveals that the conformation of the purines is predominantly C2'-endo, with the pseudorotation phase angle P lying in the range 140-180 degrees. For the pyrimidines, however, there are distortions away from this standard B-type geometry with the data being best described by P values lying in the range 90-130 degrees (i.e., O4'-endo, C1'-exo). The sugar puckers of A1, T9 and T10 are dynamically distorted no doubt as a consequence of their positions at, or close to, the ends of the duplex. Thus the conformation of the adenine and thymine sugars within the oligo(dA) and oligo(dT) strands are different with an abrupt change in sugar puckering occurring at the 5'-ApC (5'-GpT) step. Peculiar chemical shifts values for A4H2, T7CH3 and sugar C5 H1', H2' and H2", together with a number of interresidue NOEs with unusual intensities, imply that there are also substantial modifications to basepair stacking interactions at this step. Taken as a whole, our data are consistent with the view that the conformational dislocation at the 5'-ApC dinucleotide results from a combination of slide and roll manoeuvres and that the junction between the AAAA and CG sequences is a potential nucleation site for DNA bending.     

Structure of the hydration shells of oligo(dA-dT).oligo(dA-dT) and oligo(dA).oligo(dT) tracts in B-type conformation on the basis of Monte Carlo calculations

Monte Carlo simulations [(N, V, T)-ensemble] were performed for the hydration shell of poly(dA-dT).poly(dA-dT) in canonical B form and for the hydration shell of poly(dA).poly(dT) in canonical B conformation and in a conformation with narrow minor groove, highly inclined bases, but with a nearly zero-inclined base pair plane (B' conformation). We introduced helical periodic boundary conditions with a rather small unit cell and a limited number of water molecules to reduce the dimensionality of the configuration space. The coordinates of local maxima of water density and the properties of one- and two-membered water bridges between polar groups of the DNA were obtained. The AT-alternating duplex hydration mirrors the dyad symmetry of polar group distribution. At the dApdT step, a water bridge between the two carbonyl oxygens O2 of thymines is formed as in the central base-pair step of Dickerson's dodecamer. In the major groove, 5-membered water chains along the tetranucleotide pattern d(TATA).d(TATA) are observed. The hydration geometry of poly(dA).poly(dT) in canonical B conformation is distinguished by autonomous primary hydration of the base-pair edges in both grooves. When this polymer adopts a conformation with highly inclined bases and narrow minor groove, the water density distribution in the minor groove is in excellent agreement with Dickerson's spine model. One local maximum per base pair of the first layer is located near the dyad axis between adjacent base pairs, and one local maximum per base pair in the second shell lies near the dyad axis of the base pair itself. The water bridge between the two strands formed within the first layer was observed with high probability. But the water molecules of the second layer do not have a statistically favored orientation necessary for bridging first layer waters. In the major groove, the hydration geometry of the (A.T) base-pair edge resembles the main features of the AT-pair hydration derived from other sequences for the canonical B form. The preference of the B' conformation for oligo(dA).oligo(dT) tracts may express the tendency to common hydration of base-pair edges of successive base pairs in the grooves of B-type DNA. The mean potential energy of hydration of canonical B-DNA was estimated to be -60 to -80 kJ/mole nucleotides in dependence on the (G.C) contents. Because of the small system size, this estimation is preliminary.     

Structural properties of hybrid triplex of polycation deoxyribonucleic S-methylthiourea (DNmt) strands with a complementary DNA strand, probed by nanosecond molecular dynamics

The replacement of phosphodiester linkages of the polyanion DNA with S-methylthiourea linkers provides the polycation deoxyribonucleic S-methylthiourea (DNmt). Molecular dynamics studies to 1,220 ps of the hybrid triplex formed from octameric DNmt strands d(Tmt)8 with a complementary DNA oligomer strand d(Ap)8 have been carried out with explicit water solvent and Na+Cl- counterions under periodic boundary conditions using the CHARMM force field and the Ewald summation method. The Watson-Crick and Hoogsteen hydrogen-bonding patterns of the A/T tracts remained intact without any structural restraints for triplex structures throughout the simulation. The duplex portion of the triplex structure equilibrated at a B-DNA conformation in terms of the helical rise and other helical parameters. The dynamic structures of the DNmt x DNA x DNmt triplex were determined by examining histograms from the last 800 ps of the dynamics run. These included the hydrogen-bonding pattern (sequence recognition), three-centered bifurcating occurrences, minor groove width variations, and bending of tracts for the hybrid triplex structures. Together with the Watson-Crick hydrogen-bondings, the strong Hoogsteen hydrogen-bondings, the partially maintained three-centered bifurcatings in the Watson-Crick pair, and the medium-strength three-centered bifurcatings in the Hoogsteen pair suggest that the hybrid triplex is energetically favorable as compared to a duplex with similar base stacking, van der Waals interactions, and helical parameters. This is in agreement with our previously reported thermodynamic study, in which only triplex structures were observed in solution. The bending angle measured between the local axis vectors of the first and last helical axis segments is about 20 degrees for the Watson-Crick portion of the averaged structure. Propeller twist (associated with three-centered hydrogen-bonding) up to -30 degrees, native to DNA AT base pairing, was also observed for the triplex structure. The sugar pseudorotation phase angles and the ring rotation angles for the DNA strand are within the C3'-endo domain and C2'-endo domain for the DNmt strand. Water spines are observed in both minor and major grooves throughout the dynamics run. The molecular dynamics simulations of the structural properties of DNmt x DNA x DNmt hybrid triplex is compared to the DNG x DNA x DNG hybrid triplex (In DNG the -O-(PO2-)-O- linkers in DNA is replaced by -NH-C(=N+H2)-NH-).