Molecular mechanics studies of dinucleoside methylphosphonates: influence of methylphosphonate and its chirality on the phosphodiester conformation.

Molecular mechanics studies are performed on single stranded as well as base paired forms of dinucleoside methylphosphonates comprising different base sequences for both the S- and R-isomers of methylphosphonate (MP). S-MP produces noticeable distortions in the geometry, locally at the phosphate center, and this enables the stereochemical feasibility of compact g- g- phosphodiester. Besides, it tends to perturb the conformations around the P-O3' and glycosyl bonds, causing minor variations in stacking interactions. In single stranded dinucleosides, the gain in adjacent base stacking interaction energies seems to be sufficient to overcome the barrier to P-O3' bond rotation arising due to S-MP...sugar interaction, and this results in transition to a compact phosphodiester (BI-type) from an initial extended phosphodiester (BII-type) conformation. Such a thing seems rather difficult in base pair constrained duplexes. Dinucleosides with R-MP behave analogous to normal phosphate duplexes as the methyl group is away from the sugar. It is found that dinucleoside methylphosphonates are energetically less favoured than the corresponding dinucleoside phosphates mainly due to the depletion of contributions from electrostatic attractive interactions involving the base and sugar with the methylphosphonate consequent to the nonionic nature of the latter. Neither S-MP nor R-MP seem to significantly alter the stereochemistry of duplex structure.     

Molecular Mechanics Studies of Dinucleoside Methylphosphonates: Influence of Methylphosphonate and its Chirality on the Phosphodiester Conformation

Abstract

Molecular mechanics studies are performed on single stranded as well as base paired forms of dinucleoside methylphosphonates comprising different base sequences for both the Sand R-isomers of methylphosphonate (MP). S-MP produces noticeable distortions in the geometry, locally at the phosphate center, and this enables the stereochemical feasibility of compact g^? g^? phosphodiester. Besides, it tends to perturb the conformations around the P- 03′ and glycosyl bonds, causing minor variations in stacking interactions. In single stranded dinucleosides, the gain in adjacent base stacking interaction energies seems to be sufficient to overcome the barrier to P-03′ bond rotation arising due to S-MP…sugar interaction, and this results in transition to a compact phosphodiester (B_I-type) from an initial extended phosphodiester (B_II-type) conformation. Such a thing seems rather difficult in base pair constrained duplexes. Dinucleosides with R-MP behave analogous to normal phosphate duplexes as the methyl group is away from the sugar. It is found that dinucleoside methylphosphonates are energetically less favoured than the corresponding dinucleoside phosphates mainly due to the depletion of contributions from electrostatic attractive interactions involving the base and sugar with the methylphosphonate consequent to the nonionic nature of the latter. Neither S-MP nor R-MP seem to significantly alter the stereochemistry of duplex structure.     

Ab-initio quantum mechanical calculations of NMR chemical shifts in nucleic acids constituents. III. Chemical shift variations due to base stacking

Ab inito computations of the different contributions to chemical shift variations due to intra and interstrand stacking are reported for the GC, CG, AT and TA sequences of a B DNA helix. The results obtained for the non hydrogen atoms of the GC stacks show that the chemical shift variations are mainly due to the polarization contribution, the term which decreases slowly with the intermolecular distance. Because of the weaker polarity of adenine and thymine the geometric and polarization contributions are of closer absolute magnitude for the non hydrogen atoms of the intrastrand stacks but the polarization term is the determining contribution in the corresponding interstrand stacks. For the protons which undergo smaller shifts due to the polarization (or electric field effects) the role of the geometric contribution is more important and is even the leading one for the hydrogens of cytosine and thymine in the case of intrastrand stacking. The charge transfer plus exchange term has a non negligeable value for a limited number of cases corresponding to the shortest intermolecular interatomic distances. These results are discussed in relation with the qualitative differences observed between the proton and carbon spectra of dinucleotides and B-DNA duplexes.     

Different intrastrand and interstrand contributions to stacking account for roll variations at the alternating purine-pyrimidine sequences in A-DNA and A-RNA.

An explanation is suggested for the roll alternation between low and high values in A-type nucleic acid duplexes containing alternating sequences of purine and pyrimidine residues. The explanation combines two points. (1) Roll inevitably occurs in A-type duplexes due to geometrical reasons. (2) Intrastrand base stacking is much more impaired by roll than interstrand base stacking in A-type duplexes. Therefore purine-pyrimidine steps, whose bases mainly exhibit an intrastrand stacking, resist roll and decrease it. By contrast, bases at pyrimidine-purine steps exhibit a significant interstrand stacking that is tolerant to roll in A-type nucleic acid duplexes. In consequence, it is favourable if the purine-pyrimidine and pyrimidine-purine steps adopt low and high rolls, respectively in A-conformations of DNA and RNA molecules containing alternating purine-pyrimidine sequences. This is actually observed in the relevant molecular crystal structures.     

Effects of base substituents on the hydration of B- and Z-DNA: correlations to the B- to Z-DNA transition.

We present a study of how substituent groups of naturally occurring and modified nucleotide bases affect the degree of hydration of right-handed B-DNA and left-handed Z-DNA. A comparison of poly(dG-dC) and poly(dG-dm5C) titrations with the lipotropic salts of the Hofmeister series infers that the methyl stabilization of cytosines as Z-DNA is primarily a hydrophobic effect. The hydration free energies of various alternating pyrimidine-purine sequences in the two DNA conformations were calculated as solvent free energies from solvent accessible surfaces. Our analysis focused on the N2 amino group of purine bases that sits in the minor groove of the double helix. Removing this amino group from guanine to form inosine (I) destabilizes Z-DNA, while adding this group to adenines to form 2-aminoadenine (A') stabilizes Z-DNA. These predictions were tested by comparing the salt concentrations required to crystallize hexanucleotide sequences that incorporate d(CG), d(CI), d(TA) and d(TA') base pairs as Z-DNA. Combining the current results with our previous analysis of major groove substituents, we derived a thermodynamic cycle that relates the systematic addition, deletion, or substitution of each base substituent to the B- to Z-DNA transition free energy.     

Enantioselective binding of chiral 1,14-dimethyl[5]helicene–spermine ligands with B- and Z-DNA

Duplex DNA adopts a right-handed B-DNA conformation under physiological conditions. Z-DNA, meanwhile, has a left-handed helical structure and is in equilibrium with right-handed B-DNA. We recently reported that the bisnaphthyl maleimide–spermine conjugate (1) induced a B- to Z-DNA transition with high efficiency at low salt concentrations. It was also found that the bisnaphthyl ligand (1) spontaneously transformed into the corresponding [5]helicene derivative (2). Because [5]helicene 2 can potentially be chiral and because the chiral discrimination of B- and Z-DNA is also of interest, we became interested in whether enatiomerically pure [5]helicene–spermine conjugates might discriminate the chirality of B- or Z-DNA. In this study, we have demonstrated an efficient synthesis of chiral DNA-binding ligands by the conjugation of a [5]helicene unit with a spermine unit. These chiral helicene ligands exhibited recognition of B- and Z-DNA, with (P)-3 displaying preference for B-DNA and (M)-3 for Z-DNA. The characteristic features of the helicene–spermine ligands developed in this study include two points: the cationic spermine portion produces electrostatic interactions along the phosphate backbone of the minor groove, and the helicene forms complexes in an end-stacking mode. Such binding modes, together with the thermodynamic parameters, account for the mode of chiral recognition of (P)- and (M)-3 for B- and Z-DNA.     

Theoretical study of the geometrical arrangement of GT and GA wobble pairs in two short duplexes, Proton NMR chemical shifts and interaction energy calculations

NMR shielding constants are calculated for the base protons of duplexes formed by the dodecamer d(CGTGAATTCGCG) and the decamer d(CCAAGATTGG). A good agreement with experimental data is obtained for B-DNA helices in which the wobble GT and GA pairs are in the plane of the corresponding GC pairs of the parent duplexes formed by d(CGCGAATTCGCG) and d(CCAAGCTTGG), if the glycosyl bonds of T and G or A and G are symmetrical with respect to the dyad axis of the Watson-Crick GC pair. Interaction energy calculations show that this type of geometrical arrangement, which implies a distortion of the ribonphosphate backbone of both strands of the duplexes are more stable than those in which only one strand has its conformation modified by the presence of the wobble pair. For the duplex containing the GA pair, NMR chemical shifts as well as interaction energy computations favour the Watson-Crick hydrogen bonding scheme. The variation of the different contributions (intrastrand, interstrand, pair-pair) to the interaction energy between the bases of the duplexes, with the geometrical arrangement of the wobble pairs, is reported.     

Non-contiguous regions of Z-DNA in a DNA dodecamer.

The conformation of the self-complimentary DNA dodecamer d(br5CGbr5CGAATTbr5CGbr5CG) has been investigated in a variety of salt and solvent conditions by one and two-dimensional 1H NMR. In low salt aqueous solutions, the molecule forms a regular B-DNA structure similar to the unmodified dodecamer. However, in aqueous solution containing high salt concentration and methanol, the dodecamer adopts a structure in which the br5CGbr5CG ends of the molecule are in a Z-DNA like conformation and the AATT region is neither standard B-DNA nor Z-DNA. The implications of these results for the structure of junctions between B and Z-DNA and the sequence specificity of Z-DNA are discussed.     

The zab domain of the human RNA editing enzyme ADAR1 recognizes Z-DNA when surrounded by B-DNA

The Zab domain of the editing enzyme ADAR1 binds tightly and specifically to Z-DNA stabilized by bromination or supercoiling. A stoichiometric amount of protein has been shown to convert a substrate of suitable sequence to the Z form, as demonstrated by a characteristic change in the CD spectrum of the DNA. Now we show that Zab can bind not only to isolated Z-forming d(CG)(n) sequences but also to d(CG)(n) embedded in B-DNA. The binding of Zab to such sequences results in a complex including Z-DNA, B-DNA, and two B-Z junctions. In this complex, the d(CG)(n) sequence, but not the flanking region, is in the Z conformation. The presence of Z-DNA was detected by cleavage with a Z-DNA specific nuclease, by undermethylation using Z-DNA sensitive SssI methylase, and by circular dichroism. It is possible that Zab binds to B-DNA with low affinity and flips any favorable sequence into Z-DNA, resulting in a high affinity complex. Alternatively, Zab may capture Z-DNA that exists transiently in solution. The binding of Zab to potential as well as established Z-DNA segments suggests that the range of biological substrates might be wider than previously thought.     

Raman spectroscopy of DNA modified by intrastrand cross-links of antitumor cisplatin

Raman spectroscopy was employed to characterize the perturbations to DNA conformation induced in DNA by two different intrastrand adducts of antitumor cis-diamminedichloroplatinum(II) (cisplatin), namely by its 1,2-GG or 1,3-GTG intrastrand cross-links. We examined short deoxyribooligonucleotide duplexes containing single, site-specific cross-link by Raman spectroscopy and assigned the spectral alterations to conformational changes induced in DNA by 1,2-GG or 1,3-GTG intrastrand CLs determined earlier by other biochemical and biophysical methods. The results confirmed significant perturbations to the B-form DNA backbone due to the intrastrand lesions and that several nucleotides changed their conformation from C2'-endo to C3'-endo. Evidence for a partial transition from B- to A-form was found in several regions of the Raman spectra as well. The spectra also confirmed the different and more extensive distortion induced in B-DNA by 1,3-GTG in comparison with 1,2-GG intrastrand CLs, consistent with their already known high resolution structures. The results of the present work demonstrate that Raman spectroscopy represents a suitable tool to provide insights into structural factors involved in the mechanisms underlying antitumor effects of platinum drugs.     

The role of methylation in the intrinsic dynamics of B- and Z-DNA

Methylation of cytosine at the 5-carbon position (5 mC) is observed in both prokaryotes and eukaryotes. In humans, DNA methylation at CpG sites plays an important role in gene regulation and has been implicated in development, gene silencing, and cancer. In addition, the CpG dinucleotide is a known hot spot for pathologic mutations genome-wide. CpG tracts may adopt left-handed Z-DNA conformations, which have also been implicated in gene regulation and genomic instability. Methylation facilitates this B-Z transition but the underlying mechanism remains unclear. Herein, four structural models of the dinucleotide d(GC)(5) repeat sequence in B-, methylated B-, Z-, and methylated Z-DNA forms were constructed and an aggregate 100 nanoseconds of molecular dynamics simulations in explicit solvent under physiological conditions was performed for each model. Both unmethylated and methylated B-DNA were found to be more flexible than Z-DNA. However, methylation significantly destabilized the BII, relative to the BI, state through the Gp5mC steps. In addition, methylation decreased the free energy difference between B- and Z-DNA. Comparisons of α/γ backbone torsional angles showed that torsional states changed marginally upon methylation for B-DNA, and Z-DNA. Methylation-induced conformational changes and lower energy differences may contribute to the transition to Z-DNA by methylated, over unmethylated, B-DNA and may be a contributing factor to biological function.     

Luminescence of ruthenium(II) polypyridyls: evidence for intercalative binding to Z-DNA.

Photophysical studies have been undertaken to characterize the binding interactions of enantiomers of Ru(phen)3(2+), Ru(DIP)3(2+), and racemic Ru(bpy)2dppz2+ (where phen = 1,10-phenanthroline, DIP = 4,7-diphenylphenanthroline, and dppz = dipyridophenazine) with Z-form poly d(GC). Parallel enhancements in steady state luminescent intensity and a lengthening of luminescent lifetimes are seen for ruthenium enantiomers with Z-DNA as for B-DNA but with enantioselectivities reversed. Greater enhancements are seen for delta-isomers with the right-handed helix but for lambda-isomers with the left-handed helix. Ru(bpy)2dppz2+, an avid intercalator in B-DNA, displays no luminescence free in aqueous solution, but luminesces brightly bound to either B- or Z-poly d(GC). Stern-Volmer quenching studies also support the enantioselective preference in binding to B-DNA by delta-isomers and a reversal with binding to Z-DNA preferentially by the lambda-isomers. Steady state polarization studies indicate a rigid association of the complexes with both B- and Z-DNA on the time-scale of their emission and again with symmetrical enantioselectivities for the left and right-handed helices. Given the well characterized intercalative association of the complexes with B-DNA, the parallel results seen here with Z-DNA point strongly to a comparable intercalative association with the Z-form helix. That molecules may interact with Z-DNA through intercalation has not been demonstrated previously and now requires consideration in describing the range of interactions of small molecules and proteins with Z-DNA.     

Left-handed Z-DNA binding by the recA protein of Escherichia coli

recA binding to left-handed Z-DNA was measured using nitrocellulose filter binding assays with four DNA polymers with defined nucleotide sequences and four recombinant plasmids. Two to 7-fold preferential binding of recA to Z-DNA polymers was observed. Left-handed Z-DNA polymer binding by recA required ATP or its nonhydrolyzable analog, ATP(gamma S), while ADP inhibited binding. Complex formation with both B- and Z-forms was influenced by polymer length; recA bound longer DNAs better. recA binding to recombinant plasmids containing supercoil-stabilized Z-DNA was essentially similar to that found for the control vector; thus, no preferential binding of recA to the Z-form was observed. Comparative experiments with the rec1 protein of Ustilago maydis and the Escherichia coli recA protein were performed. In our hands, recA and rec1 have a similar capacity for binding left-handed Z-DNA polymers and for binding recombinant plasmids containing B- and/or Z-regions. recA contains a left-handed Z-DNA-stimulated ATPase activity. This activity differs from the right-handed B-DNA-stimulated activity since it is less sensitive to increasing pH. The kinetics of ATP hydrolysis in B-DNA/Z-DNA mixing experiments showed that the turnover of the Z-DNA recA complex was slower than for B-DNA suggesting that left-handed Z-DNA is more stably bound by recA. Our results are consistent with the postulate that left-handed Z-DNA is involved in genetic recombination.     

NMR study of hexanucleotide d(CCGCGG)2 containing two triplet repeats of fragile X syndrome

Long repeated stretches of d(CCG) and tri-nucleotide are crucial mutations that cause hereditary forms of mental retardation (fragile X-syndrome). Moreover, the alternating (CG) di-nucleotide is one of the candidates for Z-DNA conformation. Solution NMR structure of d(CCGCGG)(2) has been solved and is discussed. The determined NMR solution structure is a distorted highly bent B-DNA conformation with increased flexibility in both terminal residues. This conformation differs significantly from the Z-DNA tetramer structure reported for the same hexamer in the crystal state at similar ionic strength by Malinina and co-workers. Crystal structure of d(CCGCGG)(2) at high salt concentration includes a central alternating tetramer in Z-DNA conformation, while the initial cytosine swings out and forms a Watson-Crick base-pair with the terminal guanine of a symmetry-related molecule. In solution, NMR data for sugar ring puckering combined with restrained molecular dynamics simulations starting from a Z-DNA form show that terminal furanose residues could adopt the conformation required for aromatic bases swinging out. Therefore, tetramer formation could be considered possible once the hexanucleotide had previously adopted the Z-DNA form. This work gives some insight into correlations between anomalous crystal structures and their accessibility in the solution state.     

Mutual interactions of the phosphate groups in locally deformed backbones of various DNA double helices at high salt concentrations

Changes in the free energy of mutual phosphate group interactions are calculated that accompany bending of the A-, B- and Z-DNA backbones in 0.7, 2.1 and 4.2 mol/l NaCl aqueous solutions. The bending is often found to be favoured in the direction of the double helix grooves; B-DNA prefers bending into the major groove while minor groove is the preferred bending direction of A-DNA in the presence of 0.7 mol/l NaCl. Interestingly, the preferences are reversed in 4.2 mol/l NaCl. Further stabilization of A-DNA and B-DNA backbones is achieved in some cases if bending is combined with suitable local double helix twist alterations. Bending tendencies of Z-DNA backbone are generally weaker and they decrease, in contrast to B-DNA and A-DNA, with the increasing ionic strength.     

Structural variability and new intermolecular interactions of Z-DNA in crystals of d(pCpGpCpGpCpG).

We have determined the single crystal x-ray structure of the synthetic DNA hexamer d(pCpGpCpGpCpG) in two different crystal forms. The hexamer pCGCGCG has the Z-DNA conformation and in both cases the asymmetric unit contains more than one Z-DNA duplex. Crystals belong to the space group C222(1) with a = 69.73, b = 52.63, and c = 26.21 A, and to the space group P2(1) with a = 49.87, b = 41.26, c = 21.91 A, and gamma = 97.12 degrees. Both crystals show new crystal packing modes. The molecules also show striking new features when compared with previously determined Z-DNA structures: 1) the bases in one duplex have a large inclination with respect to the helical axis, which alters the overall shape of the molecule. 2) Some cytosine nitrogens interact by hydrogen bonding with phosphates in neighbor molecules. Similar base-phosphate interactions had been previously detected in some B-DNA crystals. 3) Basepair stacking between the ends of neighbor molecules is variable and no helical continuity is maintained between contiguous hexamer duplexes.     

Unusual frequencies of certain alternating purine-pyrimidine runs in natural DNA sequences: relation to Z-DNA

Prokaryotic, eukaryotic and mitochondrial DNA sequences of total Length 300 000 nucleotides have been analyzed to find out whether stretches of alternating purines and pyrimidines are unusual in terms of occurrence, composition and base sequence. Alternating runs longer than 5 nucleotides are significantly under-represented in the natural sequences as compared to random ones. Octanucleotides are the most deficient, occurring at only 60% of the frequency expected in random sequences. An unexpectedly high proportion of these octamers consists of alternating tetramers with the repeat structure (PuPyPuPy)2 or (PyPuPyPu)2. DNA stretches containing such sequences can potentially form a S1 nuclease sensitive slippage (staggered loop) structure, which might serve as a locally unstacked intermediate in the B- to Z-DNA conformational transition.     

Base dynamics of local Z-DNA conformations as detected by electron paramagnetic resonance with spin-labeled deoxycytidine analogues

Conformation detection and base dynamics of spin-labeled Z-DNA have been investigated by electron paramagnetic resonance (EPR) spectroscopy. The two synthesized and characterized probes used in this study were C(5)-nitroxide-labeled 2'-deoxycytidine 5'-triphosphates, pppDCAT and pppDCAVAT, which serve as suitable substrates for Micrococcus luteus DNA polymerase. Enzymatic incorporation of these probes into (dG-dC)n yields the EPR-active alternating copolymers (dG-dC,DCAT)n and (dG-dC,-DCAVAT)n. These polymers assume typical B- and Z-DNA conformations under respective low (0.1 M NaCl) and high (4.5 M NaCl) salt conditions, as evidenced by their UV-circular dichroism spectra. The EPR line shape of (dG-dC,DCAT)n in Z-form is unique and significantly different from the B-form EPR spectrum. A similar observation is made for (dG-dC,DCAVAT)n. Thus, the EPR line shapes of these spin-labeled DNAs are indicative of their local conformations. The EPR spectra, analyzed with a previously published motional model [Kao, S.-C., Polnaszek, C.F., Toppin, C.R., & Bobst, A.M. (1983) Biochemistry 22, 5563-5568], indicate tau perpendicular values of 4 and 7 ns for the B- and Z-forms, respectively. Therefore, the base dynamics of Z-DNA are about two times slower than in B-DNA.