A molecular dynamics study of the effect of G.T mispairs on the conformation of DNA in solution

The effect of G.T mispair incorporation into a double-helical environment was examined by molecular dynamics simulation. The 60-ps simulations performed on the two hexanucleotide duplexes d (G3C3)2 and d(G3TC2)2 included 10 Na+ counterions and first hydration shell waters. The resulting backbone torsional angle trajectories were analyzed to select time spans representative of conformational domains. The average backbone angles and helical parameters of the last time span for both duplexes are reported. During the simulation the hexamers retained B-type DNA structures that differed from typical A- or B-DNA forms. The overall helical structures for the two duplexes are vary similar. The presence of G.T mispairs did not alter the overall helical structure of the oligonucleotide duplex. Large propeller twist and buckle angles were obtained for both duplexes. The purine/pyrimidine crossover step showed a large decrease in propeller twist in the normal duplex but not in the mismatch duplex. Upon the formation of wobble mispairs in the mismatched duplex, the guanines moved into the minor groove and the thymines moved into the major groove. This helped prevent purine/purine clash and created a deformation in the relative orientation of the glycosidic bonds. It also exposed the free O4 of the thymines in the major groove and N2 of the guanines in the minor groove to interactions with solvent and counterions. These factors seemed to contribute to the apparently higher rigidity of the mismatched duplex during the simulation.     

Solution structure of an 11-mer duplex containing the 3, N(4)-ethenocytosine adduct opposite 2'-deoxycytidine: implications for the recognition of exocyclic lesions by DNA glycosylases

Lipid peroxidation products, as well as the metabolic products of vinyl chloride, react with cellular DNA producing the mutagenic adduct 3,N(4)-etheno-2'-deoxycytidine (epsilondC), along with several other exocyclic derivatives. High-resolution NMR spectroscopy and restrained molecular dynamics simulations were used to establish the solution structure of an 11-mer duplex containing an epsilondC.dC base-pair at its center. The NMR data suggested a regular right-handed helical structure having all residues in the anti orientation around the glycosydic torsion angle and Watson-Crick alignments for all canonical base-pairs of the duplex. Restrained molecular dynamics generated a three-dimensional model in excellent agreement with the spectroscopic data. The (epsilondC. dC)-duplex structure is a regular right-handed helix with a slight bend at the lesion site and no severe distortions of the sugar-phosphate backbone. The epsilondC adduct and its partner dC were displaced towards opposite grooves of the helix, resulting in a lesion-containing base-pair that was highly sheared but stabilized to some degree by the formation of a single hydrogen bond. Such a sheared base-pair alignment at the lesion site was previously observed for epsilondC.dG and epsilondC.T duplexes, and was also present in the crystal structures of duplexes containing dG.T and dG. U mismatches. These observations suggest the existence of a substrate structural motif that may be recognized by specific DNA glycosylases during the process of base excision repair.     

Atomic structures of excited state A–T Hoogsteen base pairs in duplex DNA by combining NMR relaxation dispersion,mutagenesis, and chemical shift calculations

NMR relaxation dispersion studies indicate that in canonical duplex DNA, Watson–Crick base pairs (bps) exist in dynamic equilibrium with short-lived low abundance excited state Hoogsteen bps. N1-methylated adenine (m¹A) and guanine (m¹G) are naturally occurring forms of damage that stabilize Hoogsteen bps in duplex DNA. NMR dynamic ensembles of DNA duplexes with m¹A–T Hoogsteen bps reveal significant changes in sugar pucker and backbone angles in and around the Hoogsteen bp, as well as kinking of the duplex towards the major groove. Whether these structural changes also occur upon forming excited state Hoogsteen bps in unmodified duplexes remains to be established because prior relaxation dispersion probes provided limited information regarding the sugar-backbone conformation. Here, we demonstrate measurements of C3′ and C4′ spin relaxation in the rotating frame (R1ρ) in uniformly ¹³C/¹⁵N labeled DNA as sensitive probes of the sugar-backbone conformation in DNA excited states. The chemical shifts, combined with structure-based predictions using an automated fragmentation quantum mechanics/molecular mechanics method, show that the dynamic ensemble of DNA duplexes containing m¹A–T Hoogsteen bps accurately model the excited state Hoogsteen conformation in two different sequence contexts. Formation of excited state A–T Hoogsteen bps is accompanied by changes in sugar-backbone conformation that allow the flipped syn adenine to form hydrogen-bonds with its partner thymine and this in turn results in overall kinking of the DNA toward the major groove. Results support the assignment of Hoogsteen bps as the excited state observed in canonical duplex DNA, provide an atomic view of DNA dynamics linked to formation of Hoogsteen bps, and lay the groundwork for a potentially general strategy for solving structures of nucleic acid excited states.     

The influence of helix morphology on co-operative polyamide backbone conformational flexibility in peptide nucleic acid complexes.

A systematic analysis of peptide nucleic acid (PNA) complexes deposited in the Protein Data Bank has been carried out using a set of contiguous atom torsion angle definitions. The analysis is complemented by molecular mechanics adiabatic potential energy calculations on hybrid PNA-nucleic acid model systems. Hitherto unobserved correlations in the values of the (alpha and epsilon) dihedral angles flanking the backbone secondary amide bond are found. This dihedral coupling forms the basis of a PNA backbone conformation classification scheme. Six conformations are thus characterised in experimental structures. Helix morphology is found to exert a significant influence on backbone conformation and flexibility: Watson-Crick PNA strands in complexes with DNA and RNA, that possess A-like base-pair stacking, adopt backbone conformations distinct from those in PNA.DNA-PNA triplex and PNA-PNA duplex P-helix forms. Solvation effects on Watson-Crick PNA backbone conformation in heterotriplexes are discussed and the possible involvement of inter-conformational transitions and dihedral angle uncoupling in asymmetric heteroduplex base-pair breathing is suggested.     

Fluorescent alpha-anomeric 1,N(6)etheno-deoxyadenosine in DNA duplexes. The alpha-epsilondA / dG pair

Structural properties of the fluorescent alpha-anomeric 1,N(6)ethenodeoxyadenosine residue placed in opposition to all four canonical deoxynucleotide units within 11-mer DNA duplexes have been studied. The duplex with alpha-epsilondA / dG pairing is most thermodynamically stable while the alpha-epsilondA / dC one is the least stable. Fluorescence measurements confirm the thermodynamic data and indicate base-pair dependent stacking properties of alpha-epsilondA within duplex structures. Results of molecular dynamics (MD) simulations in aqueous solution for the most stable duplex point to the presence of different conformational states of the alpha-1,N(6)etheno-deoxyadenosine residue, including formation of a hydrogen bonded pair with the dG and possible occurrence of severe kinking in the duplex.     

Interactions of DNA binding ligands with PNA-DNA hybrids.

The interactions of two representative mixed-sequence (one with an AT-stretch) PNA-DNA duplexes (10 or 15 base-pairs) and a PNA2/DNA triplex with the DNA binding reagents distamycin A, 4',6-diamidino-2-phenylindole (DAPI), ethidium bromide, 8-methoxy-psoralen and the delta and lambda enantiomers of Ru(phen)2-dppz2+ have been investigated using optical spectroscopic methods. The behaviour of these reagents versus two PNA-PNA duplexes has also been investigated. With triple helical poly(dA)/(H-T10-Lys-NH2)2 no significant intercalative binding was detected for any of the DNA intercalators, whereas DAPI, a DNA minor groove binder, was found to exhibit a circular dichroism with a positive sign and amplitude consistent with minor groove binding. Similarly, a PNA-DNA duplex containing a central AATA motif, a typical minor groove binding site for the DNA minor groove binders distamycin A and DAPI, showed binding for both of these drugs, though with strongly reduced affinity. No important interactions were found for any of the ligands with a PNA-DNA duplex consisting of a ten base-pair mixed purine-pyrimidine sequence with only two AT base-pairs in the centre. Nor did any of the ligands show any detectable binding to the PNA-PNA duplexes (one containing an AATT motif). Various PNA derivatives with extentions of the backbone, believed to increase the flexibility of the duplex to opening of an intercalation slot, were tested for intercalation of ethidium bromide or 8-methoxypsoralen into the mixed sequence PNA-DNA duplex, however, without any observation of improved binding. The importance of the ionic contribution of the deoxyribose phosphate backbone, versus interactions with the nucleobases, for drug binding to DNA is discussed in the light of these findings.     

Structure of the Wilms tumor suppressor protein zinc finger domain bound to DNA

The zinc finger domain of the Wilms tumor suppressor protein (WT1) contains four canonical Cys(2)His(2) zinc fingers. WT1 binds preferentially to DNA sequences that are closely related to the EGR-1 consensus site. We report the structure determination by both X-ray crystallography and NMR spectroscopy of the WT1 zinc finger domain in complex with DNA. The X-ray structure was determined for the complex with a cognate 14 base-pair oligonucleotide, and composite X-ray/NMR structures were determined for complexes with both the 14 base-pair and an extended 17 base-pair DNA. This combined approach allowed unambiguous determination of the position of the first zinc finger, which is influenced by lattice contacts in the crystal structure. The crystal structure shows the second, third and fourth zinc finger domains inserted deep into the major groove of the DNA where they make base-specific interactions. The DNA duplex is distorted in the vicinity of the first zinc finger, with a cytidine twisted and tilted out of the base stack to pack against finger 1 and the tip of finger 2. By contrast, the composite X-ray/NMR structures show that finger 1 continues to follow the major groove in the solution complexes. However, the orientation of the helix is non-canonical, and the fingertip and the N terminus of the helix project out of the major groove; as a consequence, the zinc finger side-chains that are commonly involved in base recognition make no contact with the DNA. We conclude that finger 1 helps to anchor WT1 to the DNA by amplifying the binding affinity although it does not contribute significantly to binding specificity. The structures provide molecular level insights into the potential consequences of mutations in zinc fingers 2 and 3 that are associated with Denys-Drash syndrome and nephritic syndrome. The mutations are of two types, and either destabilize the zinc finger structure or replace key base contact residues.     

The crystal structure of d(CCCCGGGG): a new A-form variant with an extended backbone conformation

The crystal structure of d(CCCCGGGG) has been determined at a resolution of 2.25 A. The oligomers crystallize as A-DNA duplexes occupying crystallographic two-fold axes. The backbone conformation is, in general, similar to that observed in previously reported crystal structures of A-DNA fragments, except for the central linkage, where it adopts an extended structure resulting from all trans conformation at the P-O5'-C5'-C4' bonds. This type of conformation facilitates interstrand stacking between the guanines at the C-G site. The local helix twist at this step is very small (25 degrees) compared to an overall average of 33.5 degrees. The unique structure of the C-G base-pair step, namely the extended backbone and the distinct stacking geometry, may be an important feature in the recognition mechanism between double-stranded DNA molecules and restriction endonucleases such as Msp I, which cuts the sequence CCGG very specifically with a rate unaffected by neighboring base pairs.     

Fluorescent α-Anomeric 1,N(6)Etheno- Deoxyadenosine in DNA Duplexes. the α-εdA / dG Pair

Abstract

Structural properties of the fluorescent α-anomeric 1,N(6)ethenodeoxyadenosine residue placed in opposition to all four canonical deoxynucleotide units within 11-mer DNA duplexes have been studied. The duplex with α-εedA / dG pairing is most thermodynamically stable while the α-edA / dC one is the least stable. Fluorescence measurements confirm the thermodynamic data and indicate base-pair dependent stacking properties of α-edA within duplex structures. Results of molecular dynamics (MD) simulations in aqueous solution for the most stable duplex point to the presence of different conformational states of the α-1,N(6)etheno-deoxyadenosine residue, including formation of a hydrogen bonded pair with the dG and possible occurrence of severe kinking in the duplex.     

Effect of environment, conformation, sequence and base substituents on the imino proton exchange rates in guanine and inosine-containing DNA, RNA, and DNA-RNA duplexes

High-resolution 1H nuclear magnetic resonance in H2O has been used to study the effect of sequence, conformation, environmental factors and base substituents on the exchange behavior of the hydrogen-bonded imino protons of guainine X cytosine and inosine X cytosine base-pairs in DNA, RNA, and DNA-RNA duplexes. The exchange rates were determined by measurement of the spin-lattice relaxation rates of the imino protons as a function of temperature. The exchange was not altered by the presence of high concentrations of salt, and the inability of phosphate to catalyze the exchange indicates that the exchange is limited by formation of a solvent-accessible "open" state. The exchange behavior depends on the duplex conformation and sequence. Exchange from the Z form polymers was orders of magnitude slower than the corresponding duplexes in the B conformation, and the A form RNA duplexes exchanged more slowly than the B form DNA polymers with the same sequence. The exchange behavior of the DNA-RNA hybrids was dependent on whether the purine or the pyrimidine strand contained the deoxyribose sugar. For both the guanine and inosine-containing duplexes, the homopolymer duplexes exchange more slowly than the more stable alternating copolymers. For the alternating duplexes, substitution of cytosine with 5-bromo- or 5-methylcytosine slowed the exchange and increased the activation energy for exchange. The inosine-containing duplexes exchanged more rapidly than the guanosine-containing duplexes, but both showed similar changes in exchange behavior in response to changes in sequence and base substituents. The activation energies for base-pair opening in B form DNA are correlated with the van der Waals contribution to the base-base interaction energy, suggesting that the purine base is partially unstacked in the open state. Using the relaxation measurements to set an upper limit on the exchange rate in poly(dG-dC) and the tritium exchange behavior at low temperature, we find that even though Z-DNA exchanges very slowly, the activation energy is similar to that observed in the A and B form duplexes, suggesting that exchange occurs from a similar open state.     

Assignments of 31P NMR resonances in oligodeoxyribonucleotides: origin of sequence-specific variations in the deoxyribose phosphate backbone conformation and the 31P chemical shifts of double-helical nucleic acids

It is now possible to unambiguously assign all 31P resonances in the 31P NMR spectra of oligonucleotides by either two-dimensional NMR techniques or site-specific 17O labeling of the phosphoryl groups. Assignment of 31P signals in tetradecamer duplexes, (dTGTGAGCGCTCACA)2, (dTAT-GAGCGCTCATA)2, (dTCTGAGCGCTCAGA)2, and (dTGTGTGCGCACACA)2, and the dodecamer duplex d(CGTGAATTCGCG)2 containing one base-pair mismatch, combined with additional assignments in the literature, has allowed an analysis of the origin of the sequence-specific variation in 31P chemical shifts of DNA. The 31P chemical shifts of duplex B-DNA phosphates correlate reasonably well with some aspects of the Dickerson/Calladine sum function for variation in the helical twist of the oligonucleotides. Correlations between experimentally measured P-O and C-O torsional angles and results from molecular mechanics energy minimization calculations show that these results are consistent with the hypothesis that sequence-specific variations in 31P chemical shifts are attributable to sequence-specific changes in the deoxyribose phosphate backbone. The major structural variation responsible for these 31P shift perturbations appears to be P-O and C-O backbone torsional angles which respond to changes in the local helical structure. Furthermore, 31P chemical shifts and JH3'-P coupling constants both indicate that these backbone torsional angle variations are more permissive at the ends of the double helix than in the middle. Thus 31P NMR spectroscopy and molecular mechanics energy minimization calculations appear to be able to support sequence-specific structural variations along the backbone of the DNA in solution.     

Base sequence effects in double-helical DNA. II. Configurational statistics of rodlike chains

Matrix generator techniques have been adapted to account for precise structural features of the nucleotide repeating unit and to translate the primary sequence of DNA base pairs into three-dimensional structures. Chains have been constructed to reflect the local sequence-dependent differences of bending and twisting of adjacent residues and various overall chain properties, including the average unperturbed moments of the end-to-end vector r and the mean angular orientation (〈γ〉 between base pair normals, 〈?₁〉 between long axes, and 〈?₂〉 between short axes) of terminal chain residues, have been computed. The chain backbone is treated implicitly in terms of the spatial fluctuations of successive base pairs. Motions are limited to low-energy perturbations of the standard B-DNA helix. Approximate potential energy schemes are used to represent the rules governing the patterns of local base–base morphology and flexibility. Theoretical predictions are compared with experimental observations at both the local and the macro-molecular level. Initial applications are limited to the rodlike poly(dA) · poly(dT) and poly(dG) · poly(dC) helices. The former duplex is found to be more compressed and the latter more extended than random-sequence DNA of the same chain length. The flexibility of the duplexes as a whole is described in terms of the average higher moments of the displacement vector ρ = r - 〈r〉 and the likelihood of chain cyclization is estimated from the three-dimensional Hermite series expansions of the displacement tensors. Emphasis is placed on theoretical methodology and the practical relevance of the calculated chain moments to observed physical properties.     

The inherent properties of DNA four-way junctions: comparing the crystal structures of holliday junctions

Holliday junctions are four-stranded DNA complexes that are formed during recombination and related DNA repair events. Much work has focused on the overall structure and properties of four-way junctions in solution, but we are just now beginning to understand these complexes at the atomic level. The crystal structures of two all-DNA Holliday junctions have been determined recently from the sequences d(CCGGGACCGG) and d(CCGGTACCGG). A detailed comparison of the two structures helps to distinguish distortions of the DNA conformation that are inherent to the cross-overs of the junctions in this crystal system from those that are consequences of the mismatched dG.dA base-pair in the d(CCGGGACCGG) structure. This analysis shows that the junction itself perturbs the sequence-dependent conformational features of the B-DNA duplexes and the associated patterns of hydration in the major and minor grooves only minimally. This supports the idea that a DNA four-way junction can be assembled at relatively low energetic cost. Both structures show a concerted rotation of the adjacent duplex arms relative to B-DNA, and this is discussed in terms of the conserved interactions between the duplexes at the junctions and further down the helical arms. The interactions distant from the strand cross-overs of the junction appear to be significant in defining its macroscopic properties, including the angle relating the stacked duplexes across the junction.     

Triple helices containing arabinonucleotides in the third (Hoogsteen) strand: effects of inverted stereochemistry at the 2'-position of the sugar moiety.

Arabinonucleic acid, the 2'-stereoisomer of RNA, was tested for its ability to recognize double-helical DNA, double-helical RNA and RNA-DNA hybrids. A pyrimidine oligoarabinonucleotide (ANA) was shown to form triple-helical complexes only with duplex DNA and hybrid DNA (Pu):RNA (Py) with an affinity that was slightly lower relative to the corresponding pyrimidine oligodeoxynucleotide (DNA) third strand. Neither the ANA nor DNA third strands were able to bind to duplex RNA or hybrid RNA (Pu):DNA (Py). In contrast, an RNA third strand recognized all four possible duplexes (DD, DR, RD and RR), as previously demonstrated. Such an understanding can be applied to the design of sequence-selective oligonucleotides which interact with double-stranded nucleic acids and emphasizes the role of the 2'-OH group as a general recognition and binding determinant of RNA.     

DNA-RNA hybrid secondary structures

DNA-RNA and DNA-DNA duplexes are even more polymorphic than observed previously. DNA-RNA hybrids can have secondary structures like A-DNA or A-RNA, but double helices of the synthetic DNA-RNA hybrids poly(dA) X poly(rU) and poly(dI) X poly(rC), respectively, form 11-fold and 10-fold double-helical structures in which the two chains have quite different conformations. Extensive X-ray fiber diffraction analyses show that in both structures the DNA chains have C-2'-endo-puckered furanose rings, while the anti-parallel RNA chains have C-3'-endo-puckered rings. The bidirectional properties of such duplexes may be important in the transfer of biological information from nucleic acids.     

The Crystal Structure of d(CCCCGGGG): A New A-Form Variant with an Extended Backbone Conformation

Abstract

The crystal structure of d(CCCCGGGG) has been determined at a resolution of 2.25Å. The oligomers crystallize as A-DNA duplexes occupying crystallographic two-fold axes. The backbone conformation is, in general, similar to that observed in previously reported crystal structures of A-DNA fragments, except for the central linkage, where it adopts an extended structure resulting from all trans conformation at the P-05′-C5′-C4′ bonds. This type of conformation facilitates interstrand stacking between the guanines at the C-G site. The local helix twist at this step is very small (25°) compared to an overall average of 33.5°. The unique structure of the C-G base-pair step, namely the extended backbone and the distinct stacking geometry, may be an important feature in the recognition mechanism between double- stranded DNA molecules and restriction endonucleases such as Msp I, which cuts the sequence CCGG very specifically with a rate unaffected by neighboring base pairs.     

Diversity of base-pair conformations and their occurrence in rRNA structure and RNA structural motifs

In addition to the canonical base-pairs comprising the standard Watson-Crick (C:G and U:A) and wobble U:G conformations, an analysis of the base-pair types and conformations in the rRNAs in the high-resolution crystal structures of the Thermus thermophilus 30S and Haloarcula marismortui 50S ribosomal subunits has identified a wide variety of non-canonical base-pair types and conformations. However, the existing nomenclatures do not describe all of the observed non-canonical conformations or describe them with some ambiguity. Thus, a standardized system is required to classify all of these non-canonical conformations appropriately. Here, we propose a new, simple and systematic nomenclature that unambiguously classifies base-pair conformations occurring in base-pairs, base-triples and base-quadruples that are associated with secondary and tertiary interactions. This system is based on the topological arrangement of the two bases and glycosidic bonds in a given base-pair. Base-pairs in the internal positions of regular secondary structure helices usually form with canonical base-pair groups (C:G, U:A, and U:G) and canonical conformations (C:G WC, U:A WC, and U:G Wb). In contrast, non-helical base-pairs outside of regular structure helices usually have non-canonical base-pair groups and conformations. In addition, many non-helical base-pairs are involved in RNA motifs that form a defined set of non-canonical conformations. Thus, each rare non-canonical conformation may be functionally and structurally important. Finally, the topology-based isostericity of base-pair conformations can rationalize base-pair exchanges in the evolution of RNA molecules.     

Interactions of molecules with nucleic acids. XI. Generalization of techniques to generate nucleic acid structures with applications to intercalation sites and kinked structures

A generalized procedure to generate nucleic acid structures is presented. In this procedure, the bases of a base pair are oriented first for characterization of particular DNA receptor sites. The resultant sites are then used in the study of specific molecule–DNA interactions. For example, intercalation sites, kinked DNA, and twisted and tilted bases are envisioned. Alterations of structures via anti → syn orientations of bases, as well as crankshaft motion about collinear bonds, provide additional conformations without disrupting the overall backbone structure. These approaches to the generation of nucleic acid structures are envisioned as required in studies of the intercalation phenomenon, minor adjustments of DNA to accommodate denaturation, binding of carcinogens to DNA, complex formation of transition metals with DNA, and antitumor agents as ligands. For these base-pair and base orientations, backbone orientations are calculated by the AGNAS technique to yield physically meaningful conformations, namely, those conformations for which nonbonded contacts are favourable. A procedure is presented to generate dimer duplex units that are physically meaningful and to assemble these units into a polynucleotide duplex. Double helices that begin with B-DNA, undergo a transition to one of the above-mentioned receptor sites, and return to B-DNA can be assembled from a catalog of dimer duplexes. Stereographic projections of the various receptor sites already being used to model binding to DNA are presented.