The effect of structural changes in a polyamine backbone on its DNA-binding properties

A novel analog of spermine, compound 1, 2, 6 bis(N-3-aminopropylmethanamine)-1-methoxy-4-methylbenzene, has been prepared which shows DNA binding which is altered from spermine in its base pair selectivity. A fluorescence spectroscopic assay is used to compare the complexation properties of compound 1, spermine, spermidine, putrescine, and berenil binding to calf thymus DNA, poly d(AT), and poly d(GC). The results are interpreted in terms of a major groove binding motif and compared with literature values for DNA dissociation constants.     

Does anthramycin bind to Z-DNA as well as to B-DNA? A molecular mechanics study

Molecular mechanics simulations have been presented on the covalent complexes between anthramycin and the decanucleotides d(GCGCGCGCGC)2 and d(GCGCGTGCGC).d(GCGCACGCGC) in the novel Z form, with the drug docked in the minor groove. The simulations predict that anthramycin binds to Z-DNA almost as well as it does to B-DNA, but with a few key differences in the structural aspects of binding. In particular, the 5' orientation of the drug is preferred in the Z form, contrasting with the preference for 3' orientation in the B form and the drug has a left-handed twist. While no experimental studies have been published on the binding of pyrrolo(1,4)benzodiazepines to left-handed forms of DNA, the energy refined stereochemically satisfactory models provide valuable information which will hopefully simulate high resolution 2-D NMR/NOE studies in solution.     

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.     

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.     

X-ray crystal structures of half the human papilloma virus E2 binding site: d(GACCGCGGTC).

The X-ray crystal structure of the DNA decamer d(GACCGCGGTC), containing half the human papilloma virus E2 binding site, has been solved from two crystals grown at different ionic conditions (50 mM MgCl2and 50 mM spermine or 1.56 mM MgCl2and 1.56 mM spermine). Despite the variation in salt concentration, the two DNA structures are in a very similar, A-type DNA conformation, with helical axes curving towards the major groove. Although the salt concentrations do not effect the helical parameters or hydration to a large degree, there is a change in the overall helical curvature; 18 degrees and 31 degrees for the low and high salt structures, respectively. This curvature appears to be sequence specific and biologically relevant when compared with similar DNA structures, including the E2 binding site of a protein-DNA complex.     

Influence of spermine on DNA conformation in a molecular dynamics trajectory of d(CGCGAATTCGCG)(2): major groove binding by one spermine molecule delays the A-->B transition

The effect of spermine on the A-DNA to B-DNA transition in d(CGCGAATTCGCG)(2) has been investigated by five A-start molecular dynamics simulations, using the Cornell et al. potential. In the absence of spermine an A-->B transition is initiated immediately and the DNA becomes equidistant from the A- and B-forms at 200ps. In three DNA-spermine simulations, when a spermine is located across the major groove of A-DNA in one of three different initial locations, the time taken to reach equidistance from the A- and B-forms is delayed until 800, 950 or 1000ps. In each case the A-form appears to be temporarily stabilized by spermine's electrostatic interactions with phosphates on both sides of the major groove. The onset of the A-->B transition can be correlated with the spermine losing contact with phosphates on one side of the groove and with A-like --> B-like sugar pucker transitions in the vicinity of the spermine bridge. However in the fifth trajectory, in which the spermine initially threads from the major groove via the backbone into the minor groove, the B-->A transition occurs rapidly once again and the DNA is equidistant between the A- and B-forms within 300ps. This indicates that the mere presence of spermine is insufficient to delay the transition and that major groove binding stabilizes A-DNA.     

Simulations of the B-DNA molecular dynamics of d(CGCGAATTCGCG)2 and d(GCGCGCGCGC)2: an analysis of the role of initial geometry and a comparison of united and all-atom models

Molecular dynamics simulations on the sequence d(CGCGAATTCGCG)2 have been carried out using both united atom and all-atom representations, and starting the simulations both from a regular repeating B-DNA structure and from the x-ray single crystal B-DNA structure. An all-atom B-DNA simulation on the sequence d(GCGCGCGCGC)2 has also been carried out, in order to compare it with a previous united atom simulation. The helix repeats, H-bonding, sugar pucker profiles, and average torsional angles are all in the range observed in crystallographic and nmr studies for B-DNA helices. In some of the sequences, there is a significant bend in the DNA helices. The individual helix repeats, with focus on 3'CpG5' and 3'GpC5' units, show the opposite helix repeat to that suggested by Calladine's rules.     

Flexibility of DNA in 2:1 drug-DNA complexes--simultaneous binding of two DAPI molecules to DNA.

Simultaneous binding of two DAPI molecules in the minor groove of (dA)15.(dT)15 B-DNA helix has been simulated by molecular mechanics calculations. The energy minimised structure shows some novel features in relation to binding of DAPI molecules as well as the flexibility of the grooves of DNA helices. The minor groove of the helix expands locally considerably (to 15 angstroms) to accommodate the two DAPI molecules and is achieved by positive propeller twisting of base pairs at the binding site concomitant with small variations in the local nucleotide stereochemistry. The expansion also brings forth simultaneously a contraction in the width of the major groove spread over to a few phosphates. These findings demonstrate another facet of the flexible stereochemistry of DNA helices in which the local features are significantly altered without being propagated beyond a few base pairs, and with the rest of the regions retaining the normal structure. Both the DAPI molecules are engaged in specific hydrogen bonds with the bases and non specific interactions with phosphates. Stacking interactions of DAPI molecules between themselves as well as with sugar-phosphate backbone contribute to the stability of the complex. The studies provide a stereochemical support to the experimental findings that under high drug-DNA ratio DAPI could bind in the 2:1 ratio.     

Dynamics of DNA condensation

The condensation of DNA induced by spermine and spermidine is investigated by equilibrium titrations and stopped-flow and field-jump experiments using scattered light detection. The spermine concentration required for the cooperative condensation process is measured at different DNA concentrations; these data are used to evaluate both the condensation threshold degree of spermine binding and the binding constant of spermine according to an excluded-site model. Stopped-flow measurements of the spermine-induced condensation demonstrate the existence of two processes: (1) A "fast" reaction is observed in the millisecond time range, when the reactant concentrations are around 1 microM; it is associated with a characteristic induction period and is assigned to the intramolecular condensation reaction. (2) A slow reaction with time constants of, e.g., 100 s strongly dependent upon both spermine and DNA concentrations is assigned to an intermolecular DNA association. The unusual time course of the intramolecular condensation reaction with the induction period provides evidence for a "threshold kinetics". During the induction period, spermine molecules are bound to DNA, but the degree of binding remains below the threshold value. As soon as the degree of ligand binding arrives at the threshold, the DNA is condensed in a relatively fast reaction. Model calculations of the spermine binding kinetics according to an excluded-site model demonstrate that the spermine molecules bound to DNA are mobile along the double helix. A comparison of the experimental data with the results of Monte Carlo simulations suggests a rate constant of approximately 200 s-1 for spermine movement by one nucleotide residue.(ABSTRACT TRUNCATED AT 250 WORDS)     

Influence of Spermine on DNA Conformation in a Molecular Dynamics Trajectory of d(CGCGAATTCGCG)2: Major Groove Binding by One Spermine Molecule Delays the A→B Transition

Abstract

The effect of spermine on the A-DNA to B-DNA transition in d(CGCGAATTCGCG)₂ has been investigated by five A-start molecular dynamics simulations, using the Cornell et al. potential. In the absence of spermine an A→B transition is initiated immediately and the DNA becomes equidistant from the A- and B-forms at 200ps. In three DNA-spermine simulations, when a spermine is located across the major groove of A-DNA in one of three different initial locations, the time taken to reach equidistance from the A- and B-forms is delayed until 800, 950 or 1000ps. In each case the A-form appears to be temporarily stabilized by spermine's electrostatic interactions with phosphates on both sides of the major groove. The onset of the A→B transition can be correlated with the spermine losing contact with phosphates on one side of the groove and with A-like → B-like sugar pucker transitions in the vicinity of the spermine bridge. However in the fifth trajectory, in which the spermine initially threads from the major groove via the backbone into the minor groove, the B→A transition occurs rapidly once again and the DNA is equidistant between the A- and B-forms within 300ps. This indicates that the mere presence of spermine is insufficient to delay the transition and that major groove binding stabilizes A-DNA.     

Analysis of mixed DNA‐bisnaphthalimide interactions involving groove association and intercalation with surface‐based and solution methodologies

The bisnaphthalimide cytotoxic agent elinafide exhibits a mixed DNA binding mode including groove‐association and intercalation. We have compared the interaction of elinafide and two bisnaphthalimide analogues with various natural and modified DNA sequences using solution NMR and UV‐melting methods and surface plasmon resonance (SPR) experiments at different pH conditions. The combined data obtained with these techniques established a high‐affinity binding mode comprising intercalation and strong electrostatic contacts with guanine bases in the major groove, and a weaker interaction with A·T pairs likely involving groove association. However, the SPR binding constants and the NMR and UV‐melting binding parameters responded differently to variations in DNA bases and ligand intercalating moieties. The rates and equilibrium constants determined by SPR clearly responded to changes in pH and DNA groove composition, but were rather insensitive to alterations in drug rings and DNA bases affecting the intercalation process. Conversely, the intermolecular stacking interactions detected by NMR and the ligand‐induced thermal stabilizations measured by UV depended on both sets of factors and were controlled by the sequence‐dependent properties of the DNA helices, indicating that these data were modulated by naphthalimide stacking in addition to groove association. A two‐step binding process where a groove‐bound state is required prior to intercalation is proposed as an explanation for these observations. These findings may be useful for studying other classes of DNA‐ and RNA‐binding drugs, which frequently combine groove‐binding and stacking moieties. © 2012 Wiley Periodicals, Inc. Biopolymers 97:974–987, 2012.     

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.     

Ternary interactions of spermine with DNA: 4''-epiadriamycin and other DNA: anthracycline complexes.

The recently developed anthracycline 4'-epiadriamycin, an anti-cancer drug with improved activity, differs from adriamycin by inversion of the stereochemistry at the 4'-position. We have cocrystallized 4'-epiadriamycin with the DNA hexamer d(CGATCG) and solved the structure to 1.5 A resolution using x-ray crystallography. One drug molecule binds at each d(CG) step of the hexamer duplex. The anthracycline sugar binds in the minor groove. A feature of this complex which distinguishes it from the earlier DNA:adriamycin complex is a direct hydrogen bond from the 4'-hydroxyl group of the anthracycline sugar to the adenine N3 on the floor of the DNA minor groove. This hydrogen bond results directly from inversion of the stereochemistry at the 4'-position. Spermine molecules bind in the major groove of this complex. In anthracycline complexes with d(CGATCG) a spermine molecule binds to a continuous hydrophobic zone formed by the 5-methyl and C6 of a thymidine, C5 and C6 of a cytidine and the chromophore of the anthracycline. This report discusses three anthracycline complexes with d(CGATCG) in which the spermine molecules have different conformations yet form extensive van der Waals contacts with the same hydrophobic zone. Our results suggest that these hydrophobic interactions of spermine are DNA sequence specific and provide insight into the question of whether DNA:spermine complexes are delocalized and dynamic or site-specific and static.     

Flexibility Of DNA In 2:1 Drug-DNA Complexes - Simultaneous Binding Of Two DAPI Molecules To DNA

Abstract

Simultaneous binding of two DAPI molecules in the minor groove of (dA)₁₅.(dT)₁₅ B-DNA helix has been simulated by molecular mechanics calculations. The energy minimised structure shows some novel features in relation to binding of DAPI molecules as well as the flexibility of the grooves of DNA helices. The minor groove of the helix expands locally considerably (to 15 Å) to accommodate the two DAPI molecules and is achieved by positive propeller twisting of base pairs at the binding site concomitant with small variations in the local nucleotide stereochemistry. The expansion also brings forth simultaneously a contraction in the width of the major groove spread over to a few phosphates. These findings demonstrate another facet of the flexible stereochemistry of DNA helices in which the local features are significantly altered without being propagated beyond a few base pairs, and with the rest of the regions retaining the normal structure. Both the DAPI molecules are engaged in specific hydrogen bonds with the bases and non specific interactions with phosphates. Stacking interactions of DAPI molecules between themselves as well as with sugar-phosphate backbone contribute to the stability of the complex. The studies provide a stereochemical support to the experimental findings that under high drug-DNA ratio DAPI could bind in the 2:1 ratio.     

Structure of the pure-spermine form of Z-DNA (magnesium free) at 1-A resolution.

We describe the three-dimensional X-ray structure of a complex of spermine bound to a Z-DNA duplex, [d(CGCGCG)]2, in the absence of any inorganic polyvalent cations. We have crystallized the DNA hexamer d(CGCGCG) in the exclusion of magnesium and other polyvalent ions and solved its structure at 1.0-A resolution. In the crystal of this pure-spermine form of Z-DNA, the relative orientation, position, and interactions of the DNA differ from the arrangement uniformly observed in over a dozen previously reported Z-DNA hexamers. Moreover, the conformation of the Z-DNA hexamer in this structure varies somewhat from those found in earlier structures. The DNA is compressed along the helical axis, the base pairs are shifted into the major groove, and the minor groove is more narrow. The packing of spermine-DNA complexes in crystals suggests that the molecular basis for the tendency of spermine to stabilize compact DNA structures derives from the capacity of spermine to interact simultaneously with several duplexes. This capacity is maximized by both the polymorphic nature and the length of the spermine cation. The length and flexibility of spermine and the dispersion of charge-charge, hydrogen-bonding, and hydrophobic bonding potential throughout the molecule maximize the ability of spermine to interact simultaneously with different DNA molecules.     

Structure and dynamics of ligand-template interactions of topoisomerase inhibitory analogs of Hoechst 33258: high field 1H-NMR and restrained molecular mechanics studies.

The binding characteristics of Hoechst 33258 (1), a synthetic bis-benzimidazole, and its structural analog 2, with one of the benzimidazoles replaced by a pyridoimidazole, to the self-complementary decadeoxyribonucleotide sequences d(CGCAATTGCG)2 (A) and d-(CATGGCCATG)2 (B) respectively, were examined using high field 1H-NMR techniques. Selective complexation induced chemical shift changes, the presence of exchange signals and intermolecular NOE contacts between the ligands and the minor groove protons of the oligonucleotides suggest the preferred binding sites as the centrally located AATT segment for complex A1, and the CCAT segment for complex B2. The B-type conformations of the two DNA duplexes are preserved upon complexation, as confirmed by the 2D-NOESY based sequential connectivities involving DNA base and sugar protons. Close intermolecular NOE based contacts between the ligands and their respective DNA sequences were further refined to model the ligand-DNA complexes starting from the computer generated B-type structures for the oligonucleotides. Force field calculations of ligand-DNA interaction energies indicate a more favorable contribution from the van der Waals energy component in the case of complex A1 consistent with its stronger net binding compared with the complex B2. Overall, the incorporation of a pyridinic nitrogen in Hoechst 33258 structure alters its selectivity for base pair recognition from A.T to G.C, resulting largely from the formation of a hydrogen bond between the new basic center and the 2-NH2 group of a guanosine moiety. The rates for the exchange of ligands between the two equivalent binding sites (AATT for 1, and CCAT for 2) of the self-complementary DNA sequences, are estimated from analyses of coalescence of NMR signals to be 189s-1 at 301 K for A1 and 79s-1 at 297 K for B2; which correspond to delta G++ of 13.8 and 18.6 kcal.mol-1 respectively.     

1 A crystal structures of B-DNA reveal sequence-specific binding and groove-specific bending of DNA by magnesium and calcium

The 1 A resolution X-ray crystal structures of Mg(2+) and Ca(2+) salts of the B-DNA decamers CCAACGTTGG and CCAGCGCTGG reveal sequence-specific binding of Mg(2+) and Ca(2+) to the major and minor grooves of DNA, as well as non-specific binding to backbone phosphate oxygen atoms. Minor groove binding involves H-bond interactions between cross-strand DNA base atoms of adjacent base-pairs and the cations' water ligands. In the major groove the cations' water ligands can interact through H-bonds with O and N atoms from either one base or adjacent bases, and in addition the softer Ca(2+) can form polar covalent bonds bridging adjacent N7 and O6 atoms at GG bases. For reasons outlined earlier, localized monovalent cations are neither expected nor found.Ultra-high atomic resolution gives an unprecedented view of hydration in both grooves of DNA, permits an analysis of individual anisotropic displacement parameters, and reveals up to 22 divalent cations per DNA duplex. Each DNA helix is quite anisotropic, and alternate conformations, with motion in the direction of opening and closing the minor groove, are observed for the sugar-phosphate backbone. Taking into consideration the variability of experimental parameters and crystal packing environments among these four helices, and 24 other Mg(2+) and Ca(2+) bound B-DNA structures, we conclude that sequence-specific and strand-specific binding of Mg(2+) and Ca(2+) to the major groove causes DNA bending by base-roll compression towards the major groove, while sequence-specific binding of Mg(2+) and Ca(2+) in the minor groove has a negligible effect on helix curvature. The minor groove opens and closes to accommodate Mg(2+) and Ca(2+) without the necessity for significant bending of the overall helix.The program Shelxdna was written to facilitate refinement and analysis of X-ray crystal structures by Shelxl-97 and to plot and analyze one or more Curves and Freehelix output files.     

A 5-nanosecond molecular dynamics trajectory for B-DNA: analysis of structure, motions, and solvation.

We report the results of four new molecular dynamics (MD) simulations on the DNA duplex of sequence d(CGCGAATTCGCG)2, including explicit consideration of solvent water, and a sufficient number of Na+ counterions to provide electroneutrality to the system. Our simulations are configured particularly to characterize the latest MD models of DNA, and to provide a basis for examining the sensitivity of MD results to the treatment of boundary conditions, electrostatics, initial placement of solvent, and run lengths. The trajectories employ the AMBER 4.1 force field. The simulations use particle mesh Ewald summation for boundary conditions, and range in length from 500 ps to 5.0 ns. Analysis of the results is carried out by means of time series for conformationalm, helicoidal parameters, newly developed indices of DNA axis bending, and groove widths. The results support a dynamically stable model of B-DNA for d(CGCGAATTCGCG)2 over the entire length of the trajectory. The MD results are compared with corresponding crystallographic and NMR studies on the d(CGCGAATTCGCG)2 duplex, and placed in the context of observed behavior of B-DNA by comparisons with the complete crystallographic data base of B-form structures. The calculated distributions of mobile solvent molecules, both water and counterions, are displayed. The calculated solvent structure of the primary solvation shell is compared with the location of ordered solvent positions in the corresponding crystal structure. The results indicate that ordered solvent positions in crystals are roughly twice as structured as bulk water. Detailed analysis of the solvent dynamics reveals evidence of the incorporation of ions in the primary solvation of the minor groove B-form DNA. The idea of localized complexation of otherwise mobile counterions in electronegative pockets in the grooves of DNA helices introduces an additional source of sequence-dependent effects on local conformational, helicoidal, and morphological structure, and may have important implications for understanding the functional energetics and specificity of the interactions of DNA and RNA with regulatory proteins, pharmaceutical agents, and other ligands.     

On the competition between water,sodium ions,and spermine in binding to DNA: a molecular dynamics computer simulation study

The interaction of DNA with the polyamine spermine(4+) (Spm(4+)), sodium ions, and water molecules has been studied using molecular dynamics computer simulations in a system modeling a DNA crystal. The simulation model consisted of three B-DNA decamers in a periodic hexagonal cell, containing 1200 water molecules, 8 Spm(4+), 32 Na(+), and 4 Cl(-) ions. The present paper gives a more detailed account of a recently published report of this system and compares results on this mixed Spm(4+)/Na(+)-cation system with an molecular dynamics simulation carried out for the same DNA decamer under similar conditions with only sodium counterions (Korolev et al., J. Mol. Biol. 308:907). The presence of Spm(4+) makes significant influence on the DNA hydration and on the interaction of the sodium ions with DNA. Spermine pushes water molecules out of the minor groove, whereas Na(+) attracts and organizes water around DNA. The major binding site of the Spm(4+) amino groups and the Na(+) ions is the phosphate group of DNA. The flexible polyamine spermine displays a high presence in the minor groove but does not form long-lived and structurally defined complexes. Sodium ions compete with Spm(4+) for binding to the DNA bases in the minor groove. Sodium ions also have several strong binding sites in the major groove. The ability of water molecules, Spm(4+), and Na(+) to modulate the local structure of the DNA double helix is discussed.