A Molecular Dynamics Simulation of a Polyamine-Induced Conformational Change of DNA. A Possible Mechanism for the B to Z Transition

Abstract

A 75ps molecular dynamics simulation has been performed on a fully solvated complex of spermine with the B DNA decamer (dGdC)₅ · (dGdC)₅. The simulation indicates a possible mechanism by which polyamines might induce the formation of a left-handed helix, the B to Z transition. Spermine was initially located in the major groove, hydrogen bonded to the helix. During the simulation the ligand migrates deeper into the DNA, maintaining strong hydrogen bonding to the central guanine bases and destroying the Watson-Crick base pairing with their respective cytosines. Significant rotation of these and other cytosine bases was observed, in part due to interactions of the helix with the aminopropyl chains of spermine. An intermediate B_II conformation might be of importance in this process.     

Molecular dynamics of spermine-DNA interactions: sequence specificity and DNA bending for a simple ligand.

We used molecular dynamics to model interactions between the physiologically important polyamine spermine and two B-DNA oligomers, the homopolymer (dG)10-(dC)10 and the heteropolymer (dGdC)5-(dGdC)5. Water and counterions were included in the simulation. Starting coordinates for spermine-DNA complexes were structures obtained by molecular mechanics modeling of spermine with the two oligomers; in these models, spermine binding induced a bend in the heteropolymer but not in the homopolymer. During approximately 40 psec of molecular dynamics simulation, spermine moves away from the floor of the major groove and interacts nospecifically with d(G)10-d(C)10. In contrast, a spermine-induced bend in the helix of (dGdC)5-(dGdC)5 is maintained throughout the simulation and spermine remains closely associated with the major groove. These results provide further evidence that the binding of spermine to nucleic acids can be sequence specific and that bending of alternating purine-pyrimidine sequences may be a physiologically important result of spermine binding.     

Molecular basis for potentiation of bleomycin-mediated degradation of DNA by polyamines. Experimental and molecular mechanical studies

The bleomycin-mediated degradation of DNA is stimulated (amplified) by certain DNA binding compounds, such as polyamines, that distort the double helix. Computer modelling studies suggest that putrescine (1), spermidine (2), and spermine (3) bind preferentially on the floor of the major groove of (dGdC)5.(dGdC)5. This interaction results in a bend of the oligomer helix toward the major groove and enlargement of the minor groove, both effects being in the order 1 less than 2 less than 3. These polyamine-induced distortions, as obtained from theoretical studies, parallel the experimental values of the amplification activities of 1-3 in the bleomycin-mediated degradation of poly(dGdC).poly(dGdC). The amplification mechanism of non-competitive binding of amplifier molecules in the major groove, and bleomycin in the minor groove, is proposed. It is suggested that the amplifier-induced conformational changes of the DNA helix increase affinity of the activated bleomycin complex toward the DNA minor groove and, consequently, result in an increased efficiency of the bleomycin-mediated degradation of the helix.     

Interaction of spermine with different forms of DNA. A conformational study

The energy of interaction of a spermine molecule with the A - and B -forms of DNA has been calculated, assuming that the molecule of spermine is fixed in the narrow groove of the DNA helix with the formation of hydrogen bonds between the amino groups of spermine and the phosphate groups of DNA. The atom–atom potentials method was used. Optimal structures for the A-DNA–spermine and B-DNA–spermine complexes are suggested. It is shown that, in agreement with the experimental data, the interaction of the spermine molecule with the A -DNA is energetically more favorable than that with the B -DNA. Two main factors are responsible for this: (1) the distance between neighboring phosphates of the chain in A -DNA (which is about 1 Å less than that in B -DNA) corresponds better to the distance between the amino groups of the propyl part of spermine; and (2) the orientation of phosphate groups in A -DNA inside the groove is preferable for complex formation with spermine to the outside groove arrangement of the phosphates in B -DNA. These conclusions are further confirmed by the calculations for DNA–propane diamine complexes.     

Molecular dynamics simulations of excimer forming (+)-anti-BPDE-DNA adducts in aqueous solution

The chemical carcinogen (+)-anti BPDE preferentially binds covalently to the guanine base in the minor groove of DNA. Fluorescence spectroscopic studies have shown that the BPDE molecules bound to DNA can interact in their photo-excited state giving strong excimer fluorescence when bound to poly(dGdC) · poly(dGdC). It was suggested that the formation of such excited state complexes is most probable when the two (+)-anti-BPDE bind to guanines of adjacent base pairs on the two different strands of the DNA. In the present work a model for such an excimer forming DNA-BPDE double adduct system has been constructed and shown to be stable over a 300 ps molecular dynamics simulation in a water box. The model is a d(CG)₃ · d(CG)₃ molecule with two BPDE molecules bound to the guanines at the 4th position on each strand, located in the minor groove and each oriented towards the 5 end of the modified strand, respectively. The results of 300 ps MD simulation show that the two BPDE chromophores exhibited on the average a relative geometry favourable for excimer formation. The local structure at the adduct position was considerably distorted and the helix axis was bent. The modified bases were found to be paired through a stable single non-Watson Crick type of hydrogen bond. Correspondence to: A. Gräslund     

dC-dG alternating oligonucleotides: thermodynamic and kinetic aspects of the B-Z transformation

The alternating cytosine-guanine oligodeoxyribonucleotides (dCdG)n, (dGdC)n, (dCdG)ndC (n = 3,4), (dGdC)7 and dG(dCdG)3 have been studied by UV and CD spectroscopy at different temperatures and NaCl concentrations. The analysis of the melting data, assuming an all-or-none model, reveals that in the B-conformation the 5'G/C3' stacking interactions are enthalpically favoured with respect to the 5'C/G3' one. The CD investigation of the B-Z equilibrium shows that the Z-conformation is enthalpically stabilized, while the B-conformation is entropically favoured, in the range of NaCl concentration considered (1 to 5 M). The kinetic data for the B-Z transformation, obtained with a salt-jump technique for the hexamer (dCdG)3, support a mechanism by which the Watson-Crick hydrogen bonds are broken before the bases flip over separately and eventually stack, reforming the H-bonds, in the new helix.     

Fluorescence property of oxazole yellow-linked oligonucleotide. Triple helix formation and photocleavage of double-stranded DNA in the presence of spermine.

Oxazole yellow is an intercalator that shows enhanced fluorescence upon binding to DNA. We prepared an oxazole yellow-linked oligonucleotide that can form triple helix with interleukin 2 receptor alpha chain promoter. The oxazole yellow-linked oligonucleotide showed linear increase of fluorescence by the triple helix formation with double-stranded DNA and also induced photocleavage of the targeted DNA in the presence of spermine upon visible illumination. Cleavage site of one strand was 7 or 8 bases away from the site of intercalation whereas the other strand was cleaved at the intercalated site.     

Inhibition of the B to Z transition in poly(dGdC).poly(dGdC) by covalent attachment of ethidium: equilibrium studies

The effects of covalent modification of poly(dGdC).poly(dGdC) and poly(dGm5dC).poly(dGm5dC) by ethidium monoazide (a photoreactive analogue of ethidium) on the salt-induced B to Z transition are examined. Earlier studies have shown ethidium monoazide to bind DNA (in the absence of light) in a manner identical to that of the parent ethidium bromide. Photolysis of the ethidium monoazide-DNA complex with visible light results in the covalent attachment of the photoreactive analogue to the DNA. This ability to form a covalent adduct was utilized to probe the effects of an intercalating irreversibly bound adduct on the salt-induced B to Z transition of the poly(dGdC).poly(dGdC) and poly(dGm5dC).poly(dGm5dC) polynucleotides. In the absence of drug, the salt-induced transition from the B to Z structure occurs in a highly cooperative manner. In contrast, this cooperativity is diminished as the concentration of covalently attached drug is increased. The degree of inhibition of the B to Z transition is quantitated as a function of the concentration of covalently attached drug. At a concentration of one drug bound per four base pairs for poly(dGdC).poly(dGdC) and seven base pairs for poly(dGm5dC).poly(dGm5dC), total inhibition of this transition is achieved. Lower concentrations of bound drug were effective in the partial inhibition of this transition. The effects of the covalently bound intercalator on the energetics of the B to Z transition were determined and demonstrated that the adduct is effective in locking the alternating copolymer in a right-handed conformation under high salt conditions.     

Ability of spermine to differentiate between DNA sequences--preferential stabilization of A-tracts

The regulatory roles fulfilled by polyamines by governance of chromatin structure are made possible by their strong association with cellular DNA, and hence by their ability to modulate DNA structure and function. Towards this end, it is crucial to understand the manifestation of sequence-dependent polyamine binding at the secondary and tertiary structural levels of DNA. This study utilizes circular dichroism (CD) and isothermal titration calorimetry (ITC) to address this relationship by using 20bp oligonucleotides with sequences-poly(dA):poly(dT), poly(dAdT):poly(dAdT), poly(dG):poly(dC), poly(dGdC):poly(dGdC)-that yield physiologically relevant structures, and poly(dIdC):poly(dIdC). CD studies show that at physiological ionic strength (150mM NaCl), spermine preferentially stabilizes A-tracts, and increases flexibility of the G-tract oligomer; the latter is also suggested by the larger change in entropy (DeltaS) of spermine binding to G-tracts. Given the chromatin destabilizing property of these sequences, these findings suggest a role for spermine in stabilization of non-nucleosomal A-tracts, and a compensating mechanism for incorporation of G-tracts in the chromatin structure. Other implications of these findings in sequence dependent DNA packaging are discussed.     

dC-dG Alternating Oligonucleotides: Thermodynamic and Kinetic Aspects of the B-Z Transformation

Abstract

The alternating cytosine-guanine oligodeoxyribonucleotides (dCdG)_n, (dGdC)_n, (dCdG)_n, (dCdG)_ndC (n=3,4), (dGdC)₇ and dG(dCdG)₃ have been studied by UV and CD spectroscopy at different temperatures and NaCl concentrations. The analysis of the melting data, assuming an all-or-none model, reveals that in the B-conformation the 5′G/C3′ stacking interactions are enthalpically favoured with respect to the 5′C/G3′ one. The CD investigation of the B-Z equilibrium shows that the Z-conformation is enthalpically stabilized, while the B-conformation is entropically favoured, in the range of NaCl concentration considered (1 to 5M). The kinetic data for the B-Z transformation, obtained with a salt-jump technique for the hexamer (dCdG)₃, support a mechanism by which (he Watson-Crick hydrogen bonds are broken before the bases flip over separately and eventually stack, reforming the H-bonds, in the new helix.     

Triplex formation at physiological pH: comparative studies on DNA triplexes containing 5-Me-dC tethered at N4 with spermine and tetraethyleneoxyamine.

Oligodeoxynucleotides with spermine conjugation at C4 of 5-Me-dC ( sp -ODN) exhibit triple helix formation with complementary Watson-Crick duplexes, and were optimally stable at physiological pH 7.3 and low salt concentration. This was attributed to a favored reassociation of the polycationic third strand with the anionic DNA duplex. To gain further insights into the factors that contribute to the enhancement of triplex stability and for engineering improved triplex systems, the spermine appendage at C4 of 5-Me-dC was replaced with 1,11-diamino-3,6,9-trioxaundecane to create teg -ODNs. From the triple helix forming abilities of these modified ODNs studied by hysteresis behaviour and the effect of salts on triplex stability, it is demonstrated here that teg- ODNs stabilise triplexes through hydrophobic desolvation while sp -ODNs stabilise triplexes by charge effects. The results imply that factors in addition to base stacking effects and interstrand hydrogen bonds are significantly involved in modulation of triplex stability by base modified oligonucleotides.     

Inhibition of the B to Z transition in poly(dGdC).poly(dGdC) by covalent attachment of ethidium: kinetic studies

The photoaffinity analogue ethidium monoazide was used to prepare samples of poly(dGdC).poly(dGdC) containing covalently attached ethidium. The effects of both noncovalently and covalently bound ethidium on the kinetics of the NaCl-induced B to Z transition in poly(dGdC).poly(dGdC) was examined using absorbance and fluorescence spectroscopy to monitor the reaction. Covalently and noncovalently attached ethidium were equal in the extent to which they reduce the rate of the B to Z transition. By using fluorescence to selectively monitor the fate of noncovalently bound ethidium over the course of the transition, we found that ethidium completely dissociates as the reaction proceeds, but at a rate that lags behind the conversion of the polymer to the Z form. These experiments provide evidence for the redistribution of noncovalently bound ethidium over the course of the B to Z transition, leading to the development of biphasic reaction kinetics. The observed kinetics suggest that the primary effect of both covalently and noncovalently bound ethidium is on the nucleation step of the B to Z transition. The reduction in the rate of the B to Z transition by noncovalently or covalently bound ethidium may be quantitatively explained as resulting from the reduced probability of finding a drug-free length of helix long enough for nucleation to occur. As necessary ancillary experiments, the defined length deoxyoligonucleotides (dGdC)4, (dGdC)5, and (dGdC)6 were synthesized and used in kinetic experiments designed to determine the nucleation length of the B to Z transition, which was found to be 6 bp. The activation energy of the B to Z transition was demonstrated to be independent of the amount of covalently bound ethidium and was found to be 21.2 +/- 1.1 kcal mol-1. Covalent attachment of ethidium was observed to increase the rate of the reverse Z to B transition, presumably by locking regions of the polymer into a right-handed conformation and thereby providing nucleation sites from which the Z to B conversion may propagate.     

Base only binding of spermine in the deep groove of the A-DNA octamer d(GTGTACAC)

The crystal structure of a complex of spermine with the DNA octamer d(GTGTACAC) has been determined at 2.0-A resolution. The alternating sequence adopts an A-DNA conformation with a novel purine-purine extra-Watson-Crick hydrogen bond involving the central guanine G3 (G11) and adenine A13 (A5) in the deep groove. The oligocation spermine binds in the floor of the deep groove by interacting with the bases and assumes an S-shape. Its dyad is coincident with that of the DNA, reminiscent of repressor binding to B-DNA. The terminal and central ammonium groups of the top half of spermine form hydrogen-bonding interactions to the 5'-bases, GTG, of one strand; then the spermine winds across the groove to interact with the corresponding set of bases on the other strand. The methylene groups of spermine form a hydrophobic cluster with the methyl groups of the thymines and the O6 atoms of the guanines of the TGT sequences on either side of the dyad. The observed mode of binding of spermine to A-DNA can serve as a model for deep groove binding in RNA and DNA-RNA hybrids that show a propensity also for the A-conformation. It will be of interest to see if base binding of spermine to DNA is involved in the regulation of gene expression, since spermine and other oligocations are ubiquitous in cells and their concentration is coupled to stages in cell cycle.     

Structural studies of DNA fragments: the G.T wobble base pair in A, B and Z DNA; the G.A base pair in B-DNA

The crystal structures of five double helical DNA fragments containing non-Watson-Crick complementary base pairs are reviewed. They comprise four fragments containing G.T base pairs: two deoxyoctamers d(GGGGCTCC) and d(GGGGTCCC) which crystallise as A type helices; a deoxydodecamer d(CGCGAATTTGCG) which crystallises in the B-DNA conformation; and the deoxyhexamer d(TGCGCG), which crystallises as a Z-DNA helix. In all four duplexes the G and T bases form wobble base pairs, with bases in the major tautomer forms and hydrogen bonds linking N1 of G with O2 of T and O6 of G with N3 of T. The X-ray analyses establish that the G.T wobble base pair can be accommodated in the A, B or Z double helix with minimal distortion of the global conformation. There are, however, changes in base stacking in the neighbourhood of the mismatched bases. The fifth structure, d(CGCGAATTAGCG), contains the purine purine mismatch G.A where G is in the anti and A in the syn conformation. The results represent the first direct structure determinations of base pair mismatches in DNA fragments and are discussed in relation to the fidelity of replication and mismatch recognition.     

Comparison of base inclination in ribo-GC and deoxybribo-GC polymers,and synthesis of poly(rGrC)-poly(rGrC)

The inclination angle between the base normal and the helix axis, and the axes around which the bases incline, are measured for ribo-GC polymers in buffer by using flow linear dichroism (LD), and compared to measurements for deoxyribo-GC polymers in buffer and under dehydrating conditions. A new method is designed to synthesize poly(rGrC) -poly(rGrC), which is not available commercially, in large quantities. The LD of this RNA reveals inclination angles that are similar to the B-form DNA in buffer, although the axes are different. The CD of poly(dGdC)-poly(dGdC) under the dehydrating conditions is similar to poly(rGrC)-poly(rGrC), indicating it is in the A form, and the LD gives larger inclination angles than either the B form or the corresponding RNA. Poly(dG)-poly(dC) is in the A form in buffer. Comparison among poly(rG)-poly(rC) in buffer, and poly (dG)-poly(dC) in buffer and under dehydrating conditions, reveals similar inclination angles and axes, although the LD shows that the DNA has the largest inclination angles. Except for poly(rGrC)-poly(rGrC), which has a unique reduced dichroism, all the axes for G are similar, as are the axes for C. © 1995 John Wiley & Sons, Inc.     

Interactions of metallic ions with DNA. V. DNA renaturation mechanism in the presence of Cu 2+

New experimental data were obtained by means of circular dichroism, melting, renaturation, and kinetic experiments, upon Cu²⁺ binding to DNA, poly dAT, and poly dGdC. They enable us to propose a model of binding giving a satisfactory explanation to all of the data found in the literature. Two types of binding sites are proposed: (a) a “sandwich” of Cu²⁺ between two adjacent G-C pairs giving a charge-transfer complex, and (b) a chelate between a phosphate group and a nitrogen atom of the bases (N₇ of guanine and N₃ of cytosine at room temperature, N₃ of adenine after thermal opening of A-T pair). Type (a) stabilizes the helix and keeps the two strands linked. Type (b) destabilizes the helix and explains why the kinetic rate of renaturation is the same as that of copper release.     

Distamycin and penta-N-methylpyrrolecarboxamide binding sites on native DNA. A comparison of methidiumpropyl-EDTA-Fe(II) footprinting and DNA affinity cleaving

Using two direct methods we have studied the binding locations and site sizes of distamycin and penta-N-methylpyrrolecarboxamide on three DNA restriction fragments from pBR322 plasmid. We find that methidiumpropyl-EDTA.Fe(II) footprinting and DNA affinity cleaving methods report common binding locations and site sizes for the tri- and pentapeptides bound to heterogeneous DNA. The tripeptide distamycin binds 5-base-pair sites with a preference for poly(dA).poly(dT) regions. The pentapeptide binds 6-7-base-pair sites with a preference for poly(dA).poly(dT) regions. These results are consistent with distamycin binding as an isogeometric helix to the minor groove of DNA with the four carboxamide N-H's hydrogen bonding five A + T base pairs. The data supports a model where each of the carboxamide N-H's can hydrogen bond to two bases, either O(2) of thymine or N(3) of adenine, located on adjacent base pairs on opposite strands of the helix. In most (but not all) cases the tri- and pentapeptide can adopt two orientations at each A + T rich binding site.     

Triplex formation at physiological pH by 5-Me-dC-N4-(spermine) [X] oligodeoxynucleotides: non protonation of N3 in X of X*G:C triad and effect of base mismatch/ionic strength on triplex stabilities.

Oligodeoxynucleotide (ODN) directed triplex formation has therapeutic importance and depends on Hoogsteen hydrogen bonds between a duplex DNA and a third DNA strand. T*A:T triplets are formed at neutral pH and C+*G:C are favoured at acidic pH. It is demonstrated that spermine conjugation at N4 of 5-Me-dC in ODNs 1-5 (sp-ODNs) imparts zwitterionic character, thus reducing the net negative charge of ODNs 1-5. sp-ODNs form triplexes with complementary 24mer duplex 8:9 show foremost stability at neutral pH 7.3 and decrease in stability towards lower pH, unlike the normal ODNs where optimal stability is found at an acidic pH 5.5. At pH 7.3, control ODNs 6 and 7 carrying dC or 5-Me-dC, respectively, do not show any triple helix formation. The stability order of triplex containing 5-Me-dC-N4-(spermine) with normal and mismatched duplex was found to be X*G:C approximately X*A:T > X*C:G > X*T:A. The hysteresis curve of sp-ODN triplex 3*8:9 indicated a better association with complementary duplex 8:9 as compared to unmodified ODN 6 in triplex 6*8:9. pH-dependent UV difference spectra suggest that N3 protonation is not a requirement for triplex formation by sp-ODN and interstrand interaction of conjugated spermine more than compensates for loss in stability due to absence of a single Hoogsteen hydrogen bond. These results may have importance in designing oligonucleotides for antigene applications.     

An NMR structural study of deaminated base pairs in DNA.

The structurally aberrant base pairs TG, UG and TI may occur in DNA as a consequence of deamination of 5-methylcytosine, cytosine and adenine respectively. Results of NMR spectroscopic studies are reported here for these deaminated base pairs in a model seven base pair long oligonucleotide duplex. We find that in all three cases, the DNA helix is a normal B form and both mispaired bases are intrahelical and hydrogen bonded with one another in a wobble geometry. Similarly, in all three cases, all sugars are found to be normal C2' endo in conformation. Symmetric structural perturbations are observed in the helix twist on the 3' side of the mispaired pyrimidine and on the 5' side of the mispaired purine. In all three cases, the amino group of the G residue on the 3' side of the mispaired pyrimidine shows hindered rotation. Although less thermodynamically stable than helices containing only Watson-Crick base pairs, these helices melt normally from the ends and not from the mispair outwards.     

Electronic structure of DNA by DV-X alpha cluster calculations: II. d(GG).d(CC), d(CG)2, d(GC)2 A and B conformations. (Part 2) Sugars and bases

The electronic structure of d(GG).d(CC), d(CG)2, d(GC)2 which are stacked base pairs in the DNA double helix, are elucidated for both A and B conformations in detail by DV-X alpha cluster calculations. These three DNA double helix fragments are contracted from the same bases, G and C, but the electronic structures of the fragments for both A and B conformations are different from each other characteristically. There are some delicate differences in the admixture of the orbital components and the overlap populations of intra- and inter- strand stacked bases among the stacking isomers. On the other hand, the electronic states of sugars differ in the 5'-3' direction, but are not almost dependent on stacked base pairs.