Binding of cationic and neutral phenanthridine intercalators to a DNA oligomer is controlled by dispersion energy: quantum chemical calculations and molecular mechanics simulations

Correlated ab initio as well as semiempirical quantum chemical calculations and molecular dynamics simulations were used to study the intercalation of cationic ethidium, cationic 5-ethyl-6-phenylphenanthridinium and uncharged 3,8-diamino-6-phenylphenanthridine to DNA. The stabilization energy of the cationic intercalators is considerably larger than that of the uncharged one. The dominant energy contribution with all intercalators is represented by dispersion energy. In the case of the cationic intercalators, the electrostatic and charge-transfer terms are also important. The DeltaG of ethidium intercalation to DNA was estimated at -4.5 kcal mol(-1) and this value agrees well with the experimental result. Of six contributions to the final free energy, the interaction energy value is crucial. The intercalation process is governed by the non-covalent stacking (including charge-transfer) interaction while the hydrogen bonding between the ethidium amino groups and the DNA backbone is less important. This is confirmed by the evaluation of the interaction energy as well as by the calculation of the free energy change. The intercalation affects the macroscopic properties of DNA in terms of its flexibility. This explains the easier entry of another intercalator molecule in the vicinity of an existing intercalation site.     

Intercalators. 1. Nature of stacking interactions between intercalators (ethidium, daunomycin, ellipticine, and 4',6-diaminide-2-phenylindole) and DNA base pairs. Ab initio quantum chemical, density functional theory, and empirical potential study

Properties of isolated intercalators (ethidium (E), daunomycin (D), ellipticine (EL), and 4,6'-diaminide-2-phenylindole (DAPI)) and their stacking interactions with adenine...thymine (AT) and guanine...cytosine (GC) nucleic acid base pairs were investigated by means of a nonempirical correlated ab initio method. All intercalators exhibit large charge delocalization, and none of them (including the DAPI dication) exhibits a site with dominant charge. All intercalators have large polarizability and are good electron acceptors, while base pairs are good electron donors. MP2/6-31G*(0.25) stabilization energies of intercalator...base pair complexes are large (E...AT, 22.4 kcal/mol; D...GC, 17.8 kcal/mol; EL...GC, 18.2 kcal/mol; DAPI...GC, 21.1 kcal/mol) and are well reproduced by modified AMBER potential (van der Waals radii of intercalator atoms are enlarged and their energy depths are increased). Standard AMBER potential underestimates binding, especially for DAPI-containing complexes. Because the DAPI dication is the best electron acceptor (among all intercalators studied), this difference is explained by the importance of the charge-transfer term, which is not included in the AMBER potential. For the neutral EL molecule, the standard AMBER force field provides correct results. The Hartree-Fock and DFT/B3LYP methods, not covering the dispersion energy, fail completely to reveal any energy minimum at the potential energy curve of the E...AT complex, and these methods thus cannot be recommended for a study of intercalation process. On the other hand, an approximate version of the DFT method, which was extended to cover London dispersion energy, yields for all complexes very good stabilization energies that are well comparable with referenced ab initio data. Besides the vertical dependence of the interaction, an energy twist dependence of the interaction energy was also investigated by a reference correlated ab initio method and empirical potentials. It is concluded that, despite the cationic (E +1, D +1, DAPI +2) or polar (EL) character of the intercalators investigated, it is the dispersion energy which predominantly contributes to the stability of intercalator...base pair complexes. Any procedure which does not cover dispersion energy is thus not suitable for studying the process of intercalation.     

Physical nature of ethidium and proflavine interactions with nucleic acid bases in the intercalation plane

On the basis of the crystallographic structures of three nucleic acid intercalation complexes involving ethidium and proflavine, we have analyzed the interaction energies between intercalator chromophores and their four nearest bases, using a hybrid variation-perturbation method at the second-order M?ller-Plesset theory level (MP2) with a 6-31G(d,p) basis set. A total MP2 interaction energy minimum precisely reproduces the crystallographic position of the ethidium chromophore in the intercalation plane between UA/AU bases. The electrostatic component constitutes the same fraction of the total energy for all three studied structures. The multipole electrostatic interaction energy, calculated from cumulative atomic multipole moments (CAMMs), was found to converge only after including components above the fifth order. CAMM interaction surfaces, calculated on grids in the intercalation planes of these structures, reasonably reproduce the alignment of intercalators in crystal structures; they exhibit additional minima in the direction of the DNA grooves, however, which also need to be examined at higher theory levels if no crystallographic data are given.     

Hydration changes for DNA intercalation reactions

The hydration changes that accompany the DNA binding of five intercalators (ethidium, propidium, proflavine, daunomycin, and 7-aminoactinomycin D) were measured by the osmotic stress method with use of the osmolytes betaine, sucrose, and triethylene glycol. Water uptake was found to accompany complex formation for all intercalators except ethidium. The difference in the number of bound water molecules between the complex and the free reactants (Deltan(w)) was different for each intercalator. The values found for Deltan(w) were the following: propidium, +6; daunomycin, +18; proflavine, +30; and 7-aminoactinomycin D, +32. For ethidium binding to DNA a value of Deltan(w) = +0.25(+/-0.3) was found, indicating that within experimental error no water was released or taken up upon complex formation. Intercalation association constants measured in D2O were found to increase relative to values measured in H2O for all compounds except ethidium. A positive correlation between the ratio of binding constants (K(D2O)/K(H2O)) and Deltan(w) was found. These combined studies identify water as an important thermodynamic participant in the formation of certain intercalation complexes.     

Micelle-sequestered dissociation of cationic DNA-intercalated drugs: unexpected surfactant-induced rate enhancement

Detergent sequestration using micelles as a hydrophobic sink for dissociated drug molecules is an established technique for determination of dissociation rates. The anionic surfactant molecules are generally assumed not to interact with the anionic DNA and thereby not to affect the rate of dissociation. By contrast, we here demonstrate that the surfactant molecules sodium dodecyl sulfate (SDS), sodium decyl sulfate, and sodium octyl sulfate all induce substantial rate enhancements of the dissociation of intercalators from DNA. Four different cationic DNA intercalators are studied with respect to surfactant-induced dissociation. Except for the smallest intercalator, ethidium, the dissociation rate constants increase monotonically with surfactant concentration both below cmc and (more strongly) above cmc, much more than expected from electrostatic effects of increased counterion concentration. The rate enhancement, most pronounced for the bulky, multicationic, hydrophobic DNA ligands in this study, indicates a reduction of the activation energy for the ligand to pass out from a deeply penetrating intercalation site of DNA. The discovery that surfactants enhance the rate of dissociation of cationic DNA-intercalators implies that rate constants previously determined by micelle-sequestered dissociation may have been overestimated. As an alternative, more reliable method, we suggest instead the addition of excess of dummy DNA as an absorbent for dissociated ligand.     

Ethidium and proflavine binding to a 2',5'-linked RNA duplex

Despite over 40 years of physical investigations, fundamental questions persist regarding the energetics of RNA and DNA intercalation. The dramatic unwinding of a nucleic acid duplex upon intercalation immediately suggests that the nucleic acid backbone should play a significant role in dictating the free energy of intercalation. However, the contribution of the backbone to intercalation free energy is difficult to appreciate given the intertwined energetics associated with intercalation (e.g., pi-pi stacking and solvent effects). Fluorescence titrations were used to determine the association constants of two known intercalators, proflavine and ethidium, for duplex 2',5'-linked RNA. Proflavine was found to bind 2',5' RNA with an association constant 25-fold greater than that measured for standard, 3',5'-linked RNA. In contrast, ethidium binds 2',5' RNA less favorably than standard RNA.     

Unmodified "GNP-oligonucleotide" nanobiohybrids: a simple route for emission enhancement of DNA intercalators

We present herein a simple method for enhancing the emission of DNA intercalators in homogeneous nanobiohybrids of unlabeled oligonucleotides and unmodified gold nanoparticles (GNPs). Pristine single‐stranded DNA (ss‐DNA) has been wrapped around unmodified GNPs to induce metal‐enhanced fluorescence (MEF) of DNA intercalators, such as ethidium bromide and propidium iodide. The thickness of the ss‐DNA layer on the gold nanosurface determines the extent of MEF, since this depends on the position of the intercalator in relation to the metal surface. Presumably, at a suitable thickness of this DNA layer, more of the intercalator is localized at the optimum distance from the nanoparticle to give rise to MEF. Importantly, no external spacer or coating agent was needed to induce the MEF effect of the GNPs. The concentration ratios of Au to DNA in the nanohybrids, as well as the capping agents applied to the GNPs, play key roles in enhancing the emission of the intercalators. The dimensions of both components of the nanobiohybrids, that is, the size of the GNPs and the length of the oligonucleotide, have considerable influences on the emission enhancement of the intercalators. Emission intensity increased with increasing size of the GNPs and length of the oligonucleotide only when the DNA efficiently wrapped the nanoparticles. An almost 100 % increment in the quantum yield of ethidium bromide was achieved with the GNP–DNA nanobiohybrid compared with that with DNA alone (in the absence of GNP), and the fluorescence emission was enhanced by 50 % even at an oligonucleotide concentration of 2 nM . The plasmonic effect of the GNPs in the emission enhancement was also established by the use of similar nanobioconjugates of ss‐DNA with nonmetallic carbon nanoparticles and TiO₂ nanoparticles, with which no increase in the fluorescence emission of ethidium bromide was observed.     

FLUORESCENCE ENERGY TRANSFER BETWEEN DIMETHYLDIAZAPEROPYRENIUM DICATION AND ETHIDIUM INTERCALATED IN POLY d(A-T)

Dimethyldiazaperopyrenium is one of the largest known DNA intercalators. Fluorescence energy transfer occurred between dimethyldiazaperopyrenium (donor) and ethidium (acceptor) when these dyes were bound to a double-stranded polynucleotide such as poly d(A-T). The addition of increasing amounts of ethidium bromide led to a marked shortening of the fluorescence lifetime of the donor, whereas the excited state of the acceptor was progressively populated via energy transfer from the donor. Critical Förster distance between these two chromophores was calculated to be 3.8 nm. The observed transfer efficiency was lower than that calculated on the basis of this critical distance and a statistical distribution of bound drugs. These results are discussed taking into account the conformational change induced by intercalation of dimethyldiazaperopyrenium in the double-stranded polynucleotide.     

Importance of model size in quantum mechanical studies of DNA intercalation

The convergence of DFT-computed interaction energies with increasing binding site model size was assessed. The data show that while accurate intercalator interaction energies can be derived from binding site models featuring only the flanking nucleotides for uncharged intercalators that bind parallel to the DNA base pairs, errors remain significant even when including distant nucleotides for intercalators that are charged, exhibit groove-binding tails that engage in noncovalent interactions with distant nucleotides, or that bind perpendicular to the DNA base pairs. Consequently, binding site models that include at least three adjacent nucleotides are required to consistently predict converged binding energies. The computationally inexpensive HF-3c method is shown to provide reliable interaction energies and can be routinely applied to such large models.     

Bridging Ligand Length Controls AT Selectivity and Enantioselectivity of Binuclear Ruthenium Threading Intercalators

The slow dissociation of DNA threading intercalators makes them interesting as model compounds in the search for new DNA targeting drugs, as there appears to be a correlation between slow dissociation and biological activity. Thus, it would be of great value to understand the mechanisms controlling threading intercalation, and for this purpose we have investigated how the length of the bridging ligand of binuclear ruthenium threading intercalators affects their DNA binding properties. We have synthesised a new binuclear ruthenium threading intercalator with slower dissociation kinetics from ct‐DNA than has ever been observed for any ruthenium complex with any type of DNA, a property that we attribute to the increased distance between the ruthenium centres of the new complex. By comparison with previously studied ruthenium complexes, we further conclude that elongation of the bridging ligand reduces the sensitivity of the threading interaction to DNA flexibility, resulting in a decreased AT selectivity for the new complex. We also find that the length of the bridging ligand affects the enantioselectivity with increasing preference for the ΔΔ enantiomer as the bridging ligand becomes longer.     

Intercalator-DNA interactions revealed by NIR surface-enhanced Raman spectroscopy

Using 1064 nm excited surface-enhanced Raman spectroscopy (SERS) a well known intercalator, ethidium bromide (EB), and a structurally related compound, 4-methyl-2,7-diamino-5,10-diphenyl-4,9-diazapyrenium hydrogensulfate (ADAP), have been studied. Concentration dependent SERS spectra of both aromatic species (1 × 10⁻⁷-5 × 10⁻⁵ M) indicated existence of dimeric associates at high concentration and an equilibrium shift towards monomers with a concentration decrease. Interactions of the intercalating molecules with DNA have been studied for various intercalator/DNA (base pair) molar ratios ranging from 10/1 to 1/10. In colloidal samples containing an intercalator in excess relative to DNA binding sites (from 10/1 to 2/1) enhancement of the Raman scattering gradually weakened, indicating a decrease in a number of free molecules adsorbed on the metal surface due to binding with DNA. At the drug/DNA ratios of 1/2 and 1/5 weaker but observable SERS bands indicated insertion of the drug molecules between the base pairs (intercalation strongly diminished interaction of the drug molecules with metal surface) as well as non-intercalative binding of the drug molecules able to stay in closer contact with a metal surface. A total intercalation of EB and ADAP molecules (intercalator/DNA of 1/7 and 1/10) resulted in almost complete loss of the SERS signal. Intensity of the SERS spectra of the intercalator/DNA complexes relative to the SERS intensity of the free intercalating molecules diminished to a lesser degree for ADAP/DNA than for EB/DNA. The obtained difference was attributed to a larger aromatic surface of the ADAP molecules which, although intercalated, could be positioned near the enhancing nanoparticles, unlike the smaller EB molecules which were deeply inserted within the DNA helix.     

Stability of Hoogsteen‐Type Triplexes – Electrostatic Attraction between Duplex Backbone and Triplex‐Forming Oligonucleotide (TFO) Using an Intercalating Conjugate

Syntheses are described for two novel twisted intercalating nucleic acid (TINA) monomers where the intercalator comprises a benzene ring linked to a naphthalimide moiety via an ethynediyl bridge. The intercalators Y and Z have a 2‐(dimethylamino)ethyl and a methyl residue on the naphthalimide moiety, respectively. When used as triplex‐forming oligonucleotides (TFOs), the novel naphthalimide TINAs show extraordinary high thermal stability in Hoogsteen‐type triplexes and duplexes with high discrimination of mismatch strands. DNA Strands containing the intercalator Y show higher thermal triplex stability than DNA strands containing the intercalator Z . This observation can be explained by the ionic interaction of the protonated dimethylamino group under physiological conditions, targeting the negatively charged phosphate backbone of the duplex. This interaction leads to an extra binding mode between the TFO and the duplex, in agreement with molecular‐modeling studies. We believe that this is the first example of an intercalator linking the TFO to the phosphate backbone of the duplex by an ionic interaction, which is a promising tool to achieve a higher triplex stability.     

二氢杨梅素-锌配合物的合成及其与DNA相互作用研究

合成了二氢杨梅素-锌配合物(DMY-Zn),采用紫外可见光谱、红外光谱、元素分析及热重分析(TG-DTA)法对其结构进行了表征,结果表明二氢杨梅素与Zn2+离子形成了配合物,其组成为[C15H10O8Zn].2H2O。采用EB为荧光探针利用荧光滴定法和粘度法进一步研究了二氢杨梅素-锌配合物与DNA的相互作用,发现二氢杨梅素-锌配合物与DNA有较强的相互作用,能和EB竞争与DNA结合,插入到ctDNA相邻的碱基对中,其作用方式为插入作用,Stern-Volmer线性猝灭常数Ksq为1.01。     

Designed Intercalators for Modification of DNA Origami Surface Properties

The modification of the backbone properties of DNA origami nanostructures through noncovalent interactions with designed intercalators, based on acridine derivatized with side chains containing esterified fatty acids or oligo(ethylene glycol) residues is reported. Spectroscopic analyses indicate that these intercalators bind to DNA origami structures. Atomic force microscopy studies reveal that intercalator binding does not affect the structural intactness but leads to altered surface properties of the highly negatively charged nanostructures, as demonstrated by their interaction with solid mica or graphite supports. Moreover, the noncovalent interaction between the intercalators and the origami structures leads to alteration in cellular uptake, as shown by confocal microscopy studies using two different eukaryotic cell lines. Hence, the intercalator approach offers a potential means for tailoring the surface properties of DNA nanostructures.     

DNA Binding Studies and Cytotoxicity of a Dinuclear PtII Diazapyrenium‐Based Metallo‐supramolecular Rectangular Box

The interaction with native DNA of a 2,7‐diazapyrenium‐based ligand 1 and its Pt^II rectangular metallacycle 2 is explored through circular and linear dichroism and fluorescence spectroscopies. The metal‐free ligand 1 binds through intercalation, with a binding constant of approximately 5×10⁵ M ^?1, whereas the metallacycle 2 binds and bends the DNA with a binding constant of 7×10⁶ M ^?1. PCR assays show that metallo‐supramolecular box 2 interferes with DNA transactions in vitro whereas the intercalator 1 does not. The metallacycle is active against four human cancer cell lines, with IC₅₀ values ranging between 3.1 and 19.2 μM and shows similar levels of efficacy, but a different spectrum of activity, to cisplatin.     

Effect of headgroup on DNA-cationic surfactant interactions

The interaction behavior of DNA with different types of hydroxylated cationic surfactants has been studied. Attention was directed to how the introduction of hydroxyl substituents at the headgroup of the cationic surfactants affects the compaction of DNA. The DNA-cationic surfactant interaction was investigated at different charge ratios by several methods like UV melting, ethidium bromide exclusion, and gel electrophoresis. Studies show that there is a discrete transition in the DNA chain from extended coils (free chain) to a compact form and that this transition does not depend substantially on the architecture of the headgroup. However, the accessibility of DNA to ethidium bromide is preserved to a significantly larger extent for the more hydrophilic surfactants. This was discussed in terms of surfactant packing. Observations are interpreted to reflect that the surfactants with more substituents have a larger headgroup and therefore form smaller micellar aggregates; these higher curvature aggregates lead to a less efficient, "patch-like" coverage of DNA. The more hydrophilic surfactants also presented a significantly lower cytotoxicity, which is important for biotechnological applications.     

Quadruplex-to-duplex transition of G-rich oligonucleotides probed by cationic water-soluble conjugated polyelectrolytes

G-quartet DNA converts to duplex form in the presence of its complementary strand. This conformational change can be detected in real time by a homogeneous assay method based on the signal amplification of conjugated polyelectrolytes and the specific interaction of intercalating dyes with double-stranded DNA (dsDNA). The probe solution contains a cationic, conjugated polymer (CCP), G-quadruplex labeled with a fluorescein at the 5'-terminus (G-quadruplex-Fl), and ethidium bromide (EB). The addition of a complementary target results in the transition from G-quadruplex to duplex (dsDNA-Fl) and EB intercalation within the duplex structure. Excitation of the CCP leads to energy transfer from CCP to dsDNA-Fl (FRET-1) and then energy transfer from dsDNA-Fl to EB (FRET-2). Increasing the number of mismatched bases discourages dsDNA formation, which is detected in the assay.     

Initial DNA Interactions of the Binuclear Threading Intercalator Λ,Λ‐[μ‐bidppz(bipy)4Ru2]4+: An NMR Study with [d(CGCGAATTCGCG)]2

Binuclear polypyridine ruthenium compounds have been shown to slowly intercalate into DNA, following a fast initial binding on the DNA surface. For these compounds, intercalation requires threading of a bulky substituent, containing one Ru^II, through the DNA base‐pair stack, and the accompanying DNA duplex distortions are much more severe than with intercalation of mononuclear compounds. Structural understanding of the process of intercalation may greatly gain from a characterisation of the initial interactions between binuclear Ru^II compounds and DNA. We report a structural NMR study on the binuclear Ru^II intercalator Λ,Λ‐B (Λ,Λ‐[μ‐bidppz(bipy)₄Ru₂]⁴⁺; bidppz=11,11′‐bis(dipyrido[3,2‐a:2′,3′‐c]phenazinyl, bipy = 2,2′‐bipyridine) mixed with the palindromic DNA [d(CGCGAATTCGCG)]₂. Threading of Λ,Λ‐B depends on the presence and length of AT stretches in the DNA. Therefore, the latter was selected to promote initial binding, but due to the short stretch of AT base pairs, final intercalation is prevented. Structural calculations provide a model for the interaction: Λ,Λ‐B is trapped in a well‐defined surface‐bound state consisting of an eccentric minor‐groove binding. Most of the interaction enthalpy originates from electrostatic and van der Waals contacts, whereas intermolecular hydrogen bonds may help to define a unique position of Λ,Λ‐B. Molecular dynamics simulations show that this minor‐groove binding mode is stable on a nanosecond scale. To the best of our knowledge, this is the first structural study by NMR spectroscopy on a binuclear Ru compound bound to DNA. In the calculated structure, one of the positively charged Ru²⁺ moieties is near the central AATT region; this is favourable in view of potential intercalation as observed by optical methods for DNA with longer AT stretches. Circular dichroism (CD) spectroscopy suggests that a similar binding geometry is formed in mixtures of Λ,Λ‐B with natural calf thymus DNA. The present minor‐groove binding mode is proposed to represent the initial surface interactions of binuclear Ru^II compounds prior to intercalation into AT‐rich DNA.     

Oxidative charge transport through DNA in nucleosome core particles

Eukaryotic DNA is packaged into nucleosomes, made up of 146 bp of DNA wrapped around a core of histone proteins. We used photoexcited rhodium intercalators to explore DNA charge transport within these assemblies. Although histone proteins inhibit intercalation of the rhodium complex within the core particle, they do not prevent 5'-GG-3' oxidation, the signature of oxidative charge transport through DNA. Moreover, using rhodium intercalators tethered to the 5' terminus of the DNA, we found that guanine bases within the nucleosome can be oxidized from a distance of 24 bp. Histone binding did not affect the pattern and extent of this oxidation. Therefore, although the structure of the nucleosome core particle generally protects DNA from damage by solution-borne molecules, packaging within the nucleosome does not protect DNA from charge transfer damage through the base pair stack.     

Voltammetric and spectroscopic studies on binding of antitumor Morin, Morin-Cu complex and Morin-beta-cyclodextrin with DNA

A systematic comparative study of the binding of antitumor Morin and its complexes with DNA has been investigated in the Britton-Robison (BR) buffer solutions using voltammetric and spectroscopic methods. The results show that Morin molecule, acting as an intercalator, is inserted into the cavity of the beta-cyclodextrin (beta-CD) as well as into the base stacking domain of the DNA double helix. The interaction of Morin-Cu complex or the inclusion complex of Morin-beta-CD with ds-DNA causes hypochromism in the absorption spectra, along with pronounced changes in the electrochemical behavior of the Morin complexes. An isobestic point and a new spectrum band appeared indicating the formation of the new system of Morin-Cu-DNA at lambda(m)=391 nm and Morin-beta-CD-DNA at lambda(m)=375 nm. The intercalation of Morin-Cu and Morin-beta-CD complexes with DNA produces an electrochemically inactive supramolecular complex. The binding constants were calculated from the increase of the solubility, the strong hypochromism, and the decrease in peak current of Morin and its complexes upon the addition of the host molecules. Calculation of the thermodynamic parameters of the interaction of the inclusion complex of Morin-beta-CD with DNA, including Gibbs free energy change, Helmholz free energy and entropy change shows that the complexation is a spontaneous process of association.