Structure of the 1,4-bis(2'-deoxyadenosin-N6-yl)-2R,3R-butanediol cross-link arising from alkylation of the human N-ras codon 61 by butadiene diepoxide

The solution structure of the 1,4-bis(2'-deoxyadenosin-N(6)-yl)-2R,3R-butanediol cross-link arising from N(6)-dA alkylation of nearest-neighbor adenines by butadiene diepoxide (BDO(2)) was determined in the oligodeoxynucleotide 5'-d(CGGACXYGAAG)-3'.5'-d(CTTCTTGTCCG)-3'. This oligodeoxynucleotide contained codon 61 (underlined) of the human N-ras protooncogene. The cross-link was accommodated in the major groove of duplex DNA. At the 5'-side of the cross-link there was a break in Watson-Crick base pairing at base pair X(6).T(17), whereas at the 3'-side of the cross-link at base pair Y(7).T(16), base pairing was intact. Molecular dynamics calculations carried out using a simulated annealing protocol, and restrained by a combination of 338 interproton distance restraints obtained from (1)H NOESY data and 151 torsion angle restraints obtained from (1)H and (31)P COSY data, yielded ensembles of structures with good convergence. Helicoidal analysis indicated an increase in base pair opening at base pair X(6).T(17), accompanied by a shift in the phosphodiester backbone torsion angle beta P5'-O5'-C5'-C4' at nucleotide X(6). The rMD calculations predicted that the DNA helix was not significantly bent by the presence of the four-carbon cross-link. This was corroborated by gel mobility assays of multimers containing nonhydroxylated four-carbon N(6),N(6)-dA cross-links, which did not predict DNA bending. The rMD calculations suggested the presence of hydrogen bonding between the hydroxyl group located on the beta-carbon of the four-carbon cross-link and T(17) O(4), which perhaps stabilized the base pair opening at X(6).T(17) and protected the T(17) imino proton from solvent exchange. The opening of base pair X(6).T(17) altered base stacking patterns at the cross-link site and induced slight unwinding of the DNA duplex. The structural data are interpreted in terms of biochemical data suggesting that this cross-link is bypassed by a variety of DNA polymerases, yet is significantly mutagenic [Kanuri, M., Nechev, L. V., Tamura, P. J., Harris, C. M., Harris, T. M., and Lloyd, R. S. (2002) Chem. Res. Toxicol. 15, 1572-1580].     

Conformational characteristics and correlations in crystal structures of nucleic acid oligonucleotides: evidence for sub-states

Sugar phosphate backbone conformations are a structural element inextricably involved in a complete understanding of specific recognition nucleic acid ligand interactions, from early stage discrimination of the correct target to complexation per se, including any structural adaptation on binding. The collective results of high resolution DNA, RNA and protein/DNA crystal structures provide an opportunity for an improved and enhanced statistical analysis of standard and unusual sugar-phosphate backbone conformations together with corresponding dinucleotide sequence effects as a basis for further exploration of conformational effects on binding. In this study, we have analyzed the conformations of all relevant crystal structures in the nucleic acids data base, determined the frequency distribution of all possible epsilon, zeta, alpha, beta and gamma backbone angle arrangements within four nucleic acid categories (A-RNA and A-DNA, free and bound B-DNA) and explored the relationships between backbone angles, sugar puckers and selected helical parameters. The trends in the correlations are found to be similar regardless of the nucleic acid category. It is interesting that specific structural effects exhibited by the different unusual backbone sub-states are in some cases contravariant. Certain alpha/gamma changes are accompanied by C3' endo (north) sugars, small twist angles and positive values of base pair roll, and favor a displacement of nucleotide bases towards the minor groove compared to that of canonical B form structures. Unusual epsilon/zeta combinations occur with C2' (south) sugars, high twist angles, negative values of base pair roll, and base displacements towards the major groove. Furthermore, any unusual backbone correlates with a reduced dispersion of equilibrium structural parameters of the whole double helix, as evidenced by the reduced standard deviations of almost all conformational parameters. Finally, a strong sequence effect is displayed in the free oligomers, but reduced somewhat in the ligand bound forms. The most variable steps are GpA and CpA, and, to a lesser extent, their partners TpC and TpG. The results provide a basis for considering if the variable and non-variable steps within a biological active sequence precisely determine morphological structural features as the curvature direction, the groove depth, and the accessibility of base pair for non covalent associations.     

Determination of nucleic acid backbone conformation by 1H NMR.

In DNA or RNA duplexes, the six-bond C3'-O3'-P-O5'-C5'-C4'-C3' backbone linkage connecting adjacent residues contains six torsion angles (epsilon, zeta, alpha, beta, gamma, delta) but only four protons. This seriously limits the ability to define the backbone conformation by NMR using purely 1H-1H distance geometry (DG) methods. The problem is further compounded by the inability to assign two of the four backbone protons, namely the poorly resolved H5' and H5' protons, and invariably leads to DG structures with poorly defined backbone conformations. We have developed and tested a reliable method to constrain the beta, gamma, and epsilon (and indirectly alpha and zeta) backbone torsion angles by lower-bound NOE distances to unassigned H5'/H5' resonances combined with either 1H line widths or the conservative use of sigma J measurements; the method relies only on 1H 2-D NMR data, does not involve any structural assumptions, and leads to much improved backbone convergence among DG structures. The C4'-C5' torsion angle gamma is constrained by lower-bound NOE distances from H2' and from H6/H8 to any H5'/H5', as well as by sigma JH4, coupling measurements in the 3.9-4.4 ppm region; delta is constrained by H1'-H4' NOE distances and by H3'-H4' and H3'-H2' J couplings in COSY data; epsilon is partially constrained by H3' line width and/or further constrained by subtracting the minimum possible sigma JH3'-H from the observed sigma JH3' (COSY) to arrive at the maximum possible JH3'-P, which is then converted to H3'-P distance bounds. The angle beta is partially constrained via H5'-P and H5'-P distance bounds consistent with the maximum H5'-P and H5'-P J couplings derived from the observed H5' and H5' line widths, while alpha and zeta are indirectly constrained by lower distance bounds on the observed (n)H1' to (n + 1)H5'/H5' NOEs combined with the prior partial constraints on beta, gamma, delta, and epsilon. The combined effects of these additional constraints in determining distance geometry structures have been demonstrated using a 12-base duplex, [d(GCCGTTAACGGC)]2. Coordinate RMSDs per atom between structures refined with these constraints from random-embedded DG structures, from ideal A-DNA, and from B-DNA starting structures were less than 0.4 A for the central 8 base pairs indicating good convergence. All backbone angles for the central 8 base pairs are very well constrained with less than 10 degrees variation in any of the 48 torsion angles.     

Conformational Characteristics and Correlations in Crystal Structures of Nucleic Acid Oligonucleotides: Evidence for Sub-states

Abstract

Sugar phosphate backbone conformations are a structural element inextricably involved in a complete understanding of specific recognition nucleic acid ligand interactions, from early stage discrimination of the correct target to complexation per se, including any structural adaptation on binding. The collective results of high resolution DNA, RNA and protein/DNA crystal structures provide an opportunity for an improved and enhanced statistical analysis of standard and unusual sugar-phosphate backbone conformations together with corresponding dinucleotide sequence effects as a basis for further exploration of conformational effects on binding. In this study, we have analyzed the conformations of all relevant crystal structures in the nucleic acids data base, determined the frequency distribution of all possible ?, ζ, α, β and γ backbone angle arrangements within four nucleic acid categories (A-RNA and A-DNA, free and bound B-DNA) and explored the relationships between backbone angles, sugar puckers and selected helical parameters. The trends in the correlations are found to be similar regardless of the nucleic acid category. It is interesting that specific structural effects exhibited by the different unusual backbone sub-states are in some cases contravariant. Certain α/γ changes are accompanied by C3′ endo (north) sugars, small twist angles and positive values of base pair roll, and favor a displacement of nucleotide bases towards the minor groove compared to that of canonical B form structures. Unusual ?/ζ combinations occur with C2′ (south) sugars, high twist angles, negative values of base pair roll, and base displacements towards the major groove. Furthermore, any unusual backbone correlates with a reduced dispersion of equilibrium structural parameters of the whole double helix, as evidenced by the reduced standard deviations of almost all con- formational parameters. Finally, a strong sequence effect is displayed in the free oligomers, but reduced somewhat in the ligand bound forms. The most variable steps are GpA and CpA, and, to a lesser extent, their partners TpC and TpG. The results provide a basis for considering if the variable and non-variable steps within a biological active sequence precisely determine morphological structural features as the curvature direction, the groove depth, and the accessibility of base pair for non covalent associations.     

Wobble dG X dT pairing in right-handed DNA: solution conformation of the d(C-G-T-G-A-A-T-T-C-G-C-G) duplex deduced from distance geometry analysis of nuclear Overhauser effect spectra

We report below on features of the three-dimensional structure of the d(C-G-T-G-A-A-T-T-C-G-C-G) self-complementary duplex (designated 12-mer GT) containing symmetrical G X T mismatches in the interior of the helix. The majority of the base and sugar protons in the 12-mer GT duplex were assigned by two-dimensional nuclear Overhauser effect (NOESY) spectra in H2O and D2O solution. A set of 92 short (less than 4.5-A) proton-proton distances defined by lower and upper bounds for one symmetrical half of the 12-mer GT duplex were estimated from NOESY data sets recorded as a function of mixing time. These experimental distances combined with nucleotide bond length parameters were embedded into Cartesian space; several trial structures were refined to minimize bond geometry and van der Waals and chirality error. Confidence in this approach is based on the similarity of the refined structures for the solution conformation of the 12-mer GT duplex. The G and T bases pair through two imino-carbonyl hydrogen bonds, and stacking is maintained between the G X T wobble pair and adjacent Watson-Crick G X C pairs. The experimental distance information is restricted to base and sugar protons, and hence structural features such as base pair overlap, glycosidic torsion angles, and sugar pucker are well-defined by this combination of NMR and distance geometry methods. By contrast, we are unable to define the torsion angles about the bonds C3'-O3'-P-O5'-C5'-C4' in the backbone of the nucleic acid.     

Human DNA polymerase iota incorporates dCTP opposite template G via a G.C + Hoogsteen base pair

Human DNA polymerase iota (hPoliota), a member of the Y family of DNA polymerases, differs in remarkable ways from other DNA polymerases, incorporating correct nucleotides opposite template purines with a much higher efficiency and fidelity than opposite template pyrimidines. We present here the crystal structure of hPoliota bound to template G and incoming dCTP, which reveals a G.C + Hoogsteen base pair in a DNA polymerase active site. We show that the hPoliota active site has evolved to favor Hoogsteen base pairing, wherein the template sugar is fixed in a cavity that reduces the C1'-C1' distance across the nascent base pair from approximately 10.5 A in other DNA polymerases to 8.6 A in hPoliota. The rotation of G from anti to syn is then largely in response to this curtailed C1'-C1' distance. A G.C+ Hoogsteen base pair suggests a specific mechanism for hPoliota's ability to bypass N(2)-adducted guanines that obstruct replication.     

Nucleic acid model building: the multiple backbone solutions associated with a given base morphology

A constrained model building procedure is used to generate nucleic acid structures of the familiar A-, B-, and Z-DNA duplexes. Attention is focused upon the multiple structural solutions associated with the arrangements of nucleic acid base pairs rather than the optimum sugar-phosphate structure. The glycosyl (chi) and sugar torsions (both the ring puckering and the exocyclic C5'-C4' (psi) torsion) are treated as independent variables and the resulting O3'...O5' distances are used as closure determinants. When such distances conform to the known geometry of phosphate chemical bonding, an intervening phosphorus atom with correct C-O-P valence angles can be located. Four sequential torsion angles--phi', omega', omega and phi--about the C3'-O3'-P-O5'-C5' bonds are then obtained as dependent variables. The resulting structures are categorized in terms of conformation, ranked in potential energy, and analyzed for torsional correlations. The numerical results are quite interesting with implications regarding nucleic acid models constructed to fit less than ideal experimental data. The multiple solutions to the problem are useful for comprehending the conformational complexities of the local sugar-phosphate backbone and for understanding the transitions between different helical forms. According to these studies, unique characterization of a nucleic acid duplex involves more than the determination of its base pair morphology, its sugar puckering preferences, or its groove binding features.     

Molecular mechanics simulation of the interaction of estrogen receptor with estrogen regulatory element.

A model for the interaction of 31 amino acid fragment (protein) from DNA binding domain of human estrogen receptor (hER) with a five base pair DNA sequence 5'GGTCA 3' from estrogen regulatory element (ERE) has been obtained using a step-wise procedure based on structural data on model peptides, DNA binding domain of hER, steric constrains imposed by tetrahedral coordination of the Cys sulphurs with zinc ion and classical secondary structural elements. Structure of the protein as well as its complex with DNA is obtained by energy minimization followed by refinement by molecular mechanics. The complex is stabilized by H-bonds between Lys22, Lys26 and Arg27 with DNA bases G2, T3 and T6. Lys22 also made H-bond with the backbone of G2. The backbone of Cys18 H-bonded with N7 of G1. DNA was in distorted B form and showed evidence of protein-induced conformational changes.     

Operator recognition by the phage 434 cI repressor: MD simulations of free and bound 50-bp DNA reveal important differences between the OR1 and OR2 sites

Using molecular dynamics simulations in explicit solvent, we investigated the behavior of a 50-bp DNA sequence containing the 434 bacteriophage operators OR1 and OR2 separated by an 8-bp spacer. Two simulations of 1 ns each were carried out, with DNA alone and with DNA complexed to dimers of the R1-69 DNA binding domain of the phage 434 cI repressor protein at the OR1 and OR2 sites. Strong correlations among average structural parameters are observed between our simulations and available experimental data for the bound OR1/OR2 subsites. In the free state, some differences appear between the three relevant fragments (OR1, the spacer, and OR2). Unbound OR1 exhibits a large, shallow major groove into which the base atoms protrude and is also bent toward the major groove. This structure is maintained because structural fluctuations are weak. Unbound OR2 resembles canonical B-DNA although the structural parameters show greater fluctuations, essentially due to a malleable step (the innermost CpA/TpG), absent in OR1. Complexation with the proteins slightly alters the base positions but strongly modifies the sugar and backbone motions. The most crucial repressor effects are changes in the flexibility of the OR1/OR2 sites. Structural fluctuations are enhanced for OR1, conferring a favorable energetic contribution to the OR1 binding, whereas they are reduced for OR2. Therefore, both structural and dynamic properties of DNA suggest OR1 is the most attractive site for the repressor, which may explain the different binding association constants observed for the OR1 and OR2 sites. Finally, we also investigated the impact of the protein on the DNA backbone dynamics and find that direct or indirect interactions facilitate the DNA structural variations required for achieving complementarity with the protein.     

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.     

Zipper-like Watson-Crick base-pairs

A series of DNA heptadecamers containing the DNA analogues of RNA E-like 5'-d(GXA)/(AYG)-5' motifs (X/Y is complementary T/A, A/T, C/G, or G/C pair) were studied using nuclear magnetic resonance (NMR) methodology and distance geometry (DG)/molecular dynamics (MD) approaches. Such oligomers reveal excellent resolution in NMR spectra and exhibit many unusual nuclear Overhauser effects (NOEs) that allow for good characterization of an unusual zipper-like conformation with zipper-like Watson-Crick base-pairs; the potential canonical X.Y H-bonding is not present, and the central X/Y pairs are transformed instead into inter-strand stacks that are bracketed by sheared G.A base-pairs. Such phenomenal structural change is brought about mainly through two backbone torsional angle adjustments, i.e. delta from C2'-endo to C3'-endo for the sugar puckers of unpaired residues and gamma from gauche(+) to trans for the following 3'-adenosine residues. Such motifs are analogous to the previously studied (GGA)(2) motif presumably present in the human centromeric (TGGAA)(n) tandem repeat sequence. The novel zipper-like motifs are only 4-7 deg. C less stable than the (GGA)(2) motif, suggesting that inter-strand base stacking plays an important role in stabilizing unusual nucleic acid structures. The discovery that canonical Watson-Crick G.C or A.T hydrogen-bonded pairs can be transformed into stacking pairs greatly increases the repertoire for unusual nucleic acid structural motifs.     

Solution structures of a DNA dodecamer duplex with and without a cisplatin 1,2-d(GG) intrastrand cross-link: comparison with the same DNA duplex containing an oxaliplatin 1,2-d(GG) intrastrand cross-link

Proteins that discriminate between cisplatin-DNA adducts and oxaliplatin-DNA adducts are thought to be responsible for the differences in tumor range, toxicity, and mutagenicity of these two important chemotherapeutic agents. However, the structural basis for differential protein recognition of these adducts has not been determined and could be important for the design of more effective platinum anticancer agents. We have determined high-resolution NMR structures for cisplatin-GG and undamaged DNA dodecamers in the AGGC sequence context and have compared these structures with the oxaliplatin-GG structure in the same sequence context determined previously in our laboratory. This structural study allows the first direct comparison of cisplatin-GG DNA and oxaliplatin-GG DNA solution structures referenced to undamaged DNA in the same sequence context. Non-hydrogen atom rmsds of 0.81 and 1.21 were determined for the 15 lowest-energy structures for cisplatin-GG DNA and undamaged DNA, respectively, indicating good structural convergence. The theoretical NOESY spectra obtained by back-calculation from the final average structures showed excellent agreement with the experimental data, indicating that the final structures are consistent with the NMR data. Several significant conformational differences were observed between the cisplatin-GG adduct and the oxaliplatin-GG adduct, including buckle at the 5' G6.C19 base pair, opening at the 3' G7.C18 base pair, twist at the A5G6.T20C19 base pair step, slide, twist, and roll at the G6G7.C19C18 base pair step, slide at the G7C8.C18G17 base pair step, G6G7 dihedral angle, and overall bend angle. We hypothesize that these conformational differences may be related to the ability of various DNA repair proteins, DNA binding proteins, and DNA polymerases to discriminate between cisplatin-GG and oxaliplatin-GG adducts.     

Effect of distortions in the deoxyribose phosphate backbone conformation of duplex oligodeoxyribonucleotide dodecamers containing GT, GG, GA, AC, and GU base-pair mismatches on 31P NMR spectra

We have previously suggested that variations in the 31P chemical shifts of individual phosphates in duplex oligonucleotides are attributable to torsional angle changes in the deoxyribose phosphate backbone. This hypothesis is not directly supported by analysis of the 1H/31P two-dimensional J-resolved spectra of a number of mismatch dodecamer oligonucleotide duplexes including the following sequences: d-(CGTGAATTCGCG), d(CGUGAATTCGCG), d(CGGGAATTCGCG), d(CGAGAATTCGCG), and d(CGCGAATTCACG). The 31P NMR signals of the dodecamer mismatch duplexes were assigned by 2D 1H/31P pure absorption phase constant time (PAC) heteronuclear correlation spectra. From the assigned H3' and H4' signals, the 31P signals of the base-pair mismatch dodecamers were identified. JH3'-P coupling constants for each of the phosphates of the dodecamers were obtained from 1H/31P J-resolved selective proton flip 2D spectra. By use of a modified Karplus relationship, the C4'-C3'-O3'-P torsional angles (epsilon) were obtained. JH3'-P coupling constants were measured for many of the oligonucleotides as a function of temperature. There exists a good linear correlation between 31P chemical shifts and the epsilon torsional angle. This correlation can be further extended to the C3'-O3'-P-O5' torsional angle (zeta) by using a linear relationship between epsilon and zeta obtained from crystal structure studies. The 31P chemical shifts follow the general observation that the more internally the phosphate is located within the oligonucleotide sequence, the more upfield the 31P resonance occurs. In addition, 31P chemical shifts show sequence- and site-specific variations. Analysis of the backbone torsional angle variations from the coupling constant analysis has provided additional information regarding the origin of these variations in 31P chemical shifts.     

Refinement of the solution structure of the DNA decamer 5'd(CTGGATCCAG)2: combined use of nuclear magnetic resonance and restrained molecular dynamics

The solution structure of the self-complementary DNA decamer 5'd(CTGGATCCAG)2 comprising the specific target site for the restriction endonuclease BamH1 is investigated by using nuclear magnetic resonance sectroscopy and restrained molecular dynamics. With the exception of the H5'/H5" sugar proton resonances, all the nonexchangeable proton resonances are assigned sequentially by using pure-phase absorption two-dimensional nuclear Overhauser enhancement spectroscopy. From the time dependence of the nuclear Overhauser effects a set of 160 approximate interproton distances is determined and used as the basis of a structure refinement employing restrained molecular dynamics in which the interproton distances are incorporated into the total energy function of the system in the form of an effective potential term. Two restrained dynamics simulations are carried out, starting from classical B- and A-DNA [atomic root mean square (rms) difference 5.7 A]. In both cases convergence is achieved to very similar B-type structures with an atomic rms difference of 0.9 A which is comparable to the rms fluctuations of the atoms about their average positions. In addition, the rms difference between the experimental and calculated values of the interproton distances for both average restrained dynamics structures is approximately 0.3 A. These results suggest that the converged restrained molecular dynamics structures represent reasonable approximations of the solution structure. The average restrained dynamics structures exhibit clear sequence-dependent variations of torsion angles and helical parameters. In addition, the structures exhibit a small bend of around 10-20 degrees at the second (TpG) and eighth (CpA) base pair steps. This can be attributed to the positive base roll angles and large base pair slide values at the two Pyr-Pur steps. The central core of the decamer comprising the six-base recognition site for BamH1 (GGATCC), however, is straight.     

Structural models for non-helical DNA.

Structural modelling techniques are employed to explore the energetic requirements for the transformation of classical B DNA into unwound yet double-stranded DNA structures. Structural idealization using CORELS computer program of Sussman et al. followed by energy minimization using the EREF program of Levitt, leads to two regular non-helical models. In both models, the bases are conventionally paired and stacked, yet there is no net rotation between successive base pairs. One model, N1, has a 1-bp repeating unit; the second, N2, has a 2-bp repeating unit. The dihedral angles of the backbone all have values found either in the B or the Z form of DNA, except for the P-O5'-C5'-C4' angle, which is in the unprecedented g+ or g- domains. The energy difference found between the two N form models and B form DNA are 6.6 and 3.4 kcal/mol/nucleotide for N1 and N2 respectively. These relatively low energy differences encourage the idea that non-helical forms of DNA may contribute to the alternate DNA structures found in S1 nuclease sensitive and other regulatory regions of active genes.     

Possible implications of CpG avoidance in the flatworm Schistosoma mansoni

We report the analysis of the biases of CpG, TpG, and CpA of all the DNA sequences data from the Trematode Schistosoma mansoni. Our results show CpG avoidance whereas TpG and CpA frequencies are over the expected values. These characteristics are similar to the biases displayed by methylated genomes, but in platyhelminths 5mC has not been detected by biochemical methods. The possible implications of this CpG shortage are discussed.