Cooperative effects on the formation of intercalation sites

Computational modeling to understand the interaction of TMPyP4 with a G-quadruplex

The potential of small molecules to bind to G-quadruplex-forming sequences in oncogene promoter regions, thereby regulating their structural equilibrium, has been explored as a promising strategy for cancer chemotherapy. The model drug 5,10,15,20-tetrakis-(N-methyl-4-pyridyl)porphine (TMPyP4) has been shown to have an affinity toward G-quadruplex DNA. However, the precise sites and modes of TMPyP4 binding to G-quadruplex DNA remain a subject of debate. In this study, we focus on identifying potential binding sites on a mutant c-MYC sequence known to fold into a single 1:2:1 loop isomer quadruplex. Our findings provide insights into the 4:1 stoichiometry reported for TMPyP4 binding to this G-quadruplex. Binding enthalpy and free energy calculations show that intercalation of a TMPyP4 molecule between the quadruplexes is thermodynamically favorable. Our calculations suggest that two of the binding sites are located at the top and bottom of the quadruplex, respectively, while the remaining two are likely intercalations.

Interactions of a Novel Anthracycline with Oligonucleotide DNA and Cyclodextrins in an Aqueous Environment

Berubicin, a chemotherapy medication belonging to the class of anthracyclines, is simulated in double-stranded DNA sequences and cyclodextrins in an aqueous environment via full-atom molecular dynamics simulations on the time scale of microseconds. The drug is studied in both the neutral and protonated states so as to better comprehend the role of its charge in the formed complexes. The noncovalent berubicin-DNA and berubicin-cyclodextrin complexes are investigated in detail, paying special attention to their thermodynamic description by employing the double decoupling method, the solvent balance method, the weighted solvent accessible surface model, and the linear interaction energy method. A novel approach for extracting the desolvation thermodynamics of the binding process is also presented. Both the binding and desolvation Gibbs energies are decomposed into entropic and enthalpic contributions so as to elucidate the nature of complexation and its driving forces. Selected structural and geometrical properties of all the complexes, which are all stable, are analyzed. Both cyclodextrins under consideration are widely utilized for drug delivery purposes, and a comparative investigation between their bound states with berubicin is carried out.

Reactions of deoxyribonucleotide bases with sulfooxymethyl or halomethyl polycyclic aromatic hydrocarbons induce unwinding of DNA supercoils

Torsional stress in double-stranded DNA enables and regulates facets of chromosomal metabolism, replication, and transcription and requires regulatory enzymatic systems including topoisomerases and histone methyltransferases. As such, this machinery may be subject to deleterious effects from reactive mutagens, including ones from carcinogenic polycyclic aromatic hydrocarbon (PAH) adduct formation with DNA. Supercoiled plasmid DNA was investigated for its torsional responses to adducts formed in vitro from PAH benzylic carbocation reactive intermediates created spontaneously by release of leaving groups. PAH sulfate esters were found to (1) unwind DNA in a concentration dependent manner, and (2) provide maximum unwinding in a pattern consistent with known carcinogenicities of the parent PAHs, that is, 6-methylbenzo[a]pyrene > 7,12-methylbenz[a]anthracene > 3-methylcholanthrene > 9-methylanthracene > 7-methylbenz[a]anthracene > 1-methylpyrene. Supercoil unwinding was demonstrated to be dependent on the presence of sulfate or chloride leaving groups such that reactive carbocations were generated in situ by hydrolysis. In silico modeling of intercalative complex topology showed PAH benzylic carbocation reactive functional groups in alignment with target nucleophiles on guanine bases in a 5'-dCdG-3' pocket in agreement with known formation of nucleotide adducts. Inhibitory or modulatory effects on PAH-induced supercoil unwinding were seen with ascorbic acid and an experimental antineoplastic agent Antineoplaston A10 in agreement with their known anticarcinogenic properties. In summary, the reactive PAH intermediates studied here undoubtedly participate in well-known mutational mechanisms such as frameshifts and apurinic site generation. However, they are also capable of random disruption of chromosomal supercoiling in a manner consistent with the known carcinogenicities of the parent compounds, and this mechanism may represent an additional detrimental motif worthy of further study for a more complete understanding of chemical carcinogenicity.

Influence of conventional hydrogen bonds in the intercalation of phenanthroline derivatives with DNA: The important role of the sugar and phosphate backbone

The influence of hydrogen bonds in model intercalated systems between guanine-cytosine and adenine-thymine DNA base pairs (bps) was analyzed with the popular intercalator 1,10-phenanthroline (phen) and derivatives obtained by substitution with OH and NH2 groups in positions 4 and 7. Semiempirical and Density Functional Theory (DFT) methods were used both including dispersion effects: PM6-DH2, M06-2X and B3LYP-D3 along with the recently developed near linear-scaling coupled cluster method DLPNO-CCSD(T) for benchmark calculations. Our results given by QTAIM and non-covalent interaction analysis confirmed the existence of hydrogen bonds created by OH and NH2 . The trends in the energy decomposition analysis for the interaction energy, ΔEint , showed that the ΔEelstat contributions are equal or even a little bit higher than the values for ΔEdisp . Such important ΔEelstat attractive contribution comes mainly from the conventional hydrogen bonds formed by OH and NH2 functional groups with DNA not only with bps but specially with the sugar and phosphate backbone. This behavior is very different from that of phen and other classical intercalators that cannot form conventional hydrogen bonds, where the ΔEdisp is the most important attractive contribution to the ΔEint . The inclusion of explicit water molecules in molecular dynamics simulations showed, as a general trend, that the hydrogen bonds with the bps disappear during the simulations but those with the sugar and phosphate backbone remain in time, which highlights the important role of the sugar and phosphate backbone in the stabilization of these systems.

From groove binding to intercalation: unravelling the weak interactions and other factors modulating the modes of interaction between methylated phenanthroline-based drugs and duplex DNA

Several antitumor drugs base their cytotoxicity on their capacity to intercalate between base pairs of DNA. Nevertheless, it has been established that the mechanism of intercalation of drugs in DNA starts with the prior groove binding mode of interaction of the drug with DNA. Sometimes, for some kind of flat small molecules, groove binding does not produce any cytotoxic effect and the fast transition of such flat small molecules to the cytotoxic intercalation mode is desirable. This is the case of methylated phenanthroline (phen) derivatives, where, changes in the substitution in the position and number of methyl groups determine their capability as cytotoxic compounds and, therefore, it is a way for the modulation of cytotoxic effects. In this work, we studied this modulation by means of the interaction of the [Pt(en)(phen)]²⁺ complex and several derivatives by methylation of phen in different number and position and the d(GTCGAC)₂ DNA hexamer via groove binding using PM6-DH2 and DFT-D methods. The analysis of the geometries, electronic structure and energetics of the studied systems was compared to experimental works to gain insight into the relation structure-interaction for the studied systems with cytotoxicity. The trends are explained by means of the Non-Covalent Interaction (NCI) index, the Energy Decomposition Analysis (EDA) and solvation contributions. Our results are in agreement with the experiments, in which the methylation of position 4 of phen seems to favour the interaction via groove binding thus making the transition to the intercalation cytotoxic mode difficult. Looking at the NCI results, these interactions come not only from the CH/π and CH/n interactions of the methyl group in position 4 but also from the ethylenediamine (en) ligand, whose orientation in the Pt complex was found in such a way that it produces a high number of weak interactions with DNA, especially with the sugar and phosphate backbone.

Thermodynamic Dissection of the Intercalation Binding Process of Doxorubicin to dsDNA with Implications of Ionic and Solvent Effects

Doxorubicin (DOX) is a cancer drug that binds to dsDNA through intercalation. A comprehensive microsecond timescale molecular dynamics study is performed for DOX with 16 tetradecamer dsDNA sequences in explicit aqueous solvent, in order to investigate the intercalation process at both binding stages (conformational change and insertion binding stages). The molecular mechanics generalized Born surface area (MM-GBSA) method is adapted to quantify and break down the binding free energy (BFE) into its thermodynamic components, for a variety of different solution conditions as well as different DNA sequences. Our results show that the van der Waals interaction provides the largest contribution to the BFE at each stage of binding. The sequence selectivity depends mainly on the base pairs located downstream from the DOX intercalation site, with a preference for (AT)₂ or (TA)₂ driven by the favorable electrostatic and/or van der Waals interactions. Invoking the quartet sequence model proved to be most successful to predict the sequence selectivity. Our findings also indicate that the aqueous bathing solution (i.e., water and ions) opposes the formation of the DOX-DNA complex at every binding stage, thus implying that the complexation process preferably occurs at low ionic strength and is crucially dependent on solvent effects.

Investigating the impacts of DNA binding mode and sequence on thermodynamic quantities and water exchange values for two small molecule drugs

Doxorubicin and nogalamycin are antitumor antibiotics that interact with DNA via intercalation and threading mechanisms, respectively. Because the importance of water, particularly its impact on entropy changes, has been established in other biological processes, we investigated the role of water in these two drug-DNA binding events. We used the osmotic stress method to calculate the number of water molecules exchanged (Δnwater), and isothermal titration calorimetry to measure Kbinding, ΔH, and ΔS for two synthetic DNAs, poly(dA·dT) and poly(dG·dC), and calf thymus DNA (CT DNA). For nogalamycin, Δnwater<0 for CT DNA and poly(dG·dC). For doxorubicin, Δnwater>0 for CT DNA and Δnwater<0 for poly(dG·dC). For poly(dA·dT), Δnwater~0 with both drugs. Net enthalpy changes were always negative, but net entropy changes depended on the drug. The effect of water exchange on the overall sign of entropy change appears to be smaller than other contributions.

A theoretical study of methylation and CH/π interactions in DNA intercalation: methylated 1,10-phenanthroline in adenine–thymine base pairs

This work shows that quality is better that quantity to estabilize the intercalation of methylated phen.

The fluorescence properties and binding mechanism of SYTOX green, a bright, low photo-damage DNA intercalating agent

DNA intercalators are widely used in cancer therapeutics, to probe protein-DNA interactions and to investigate the statistical-mechanical properties of DNA. Here, we employ single-molecule fluorescence microscopy, magnetic tweezers, and ensemble-binding assays to investigate the fluorescence properties and binding mechanism of SYTOX green, a DNA labeling dye previously used for staining dead cells and becoming of common use for single-molecule methodologies. Specifically, we show that SYTOX green presents several advantages with respect to other dyes: (1) binds DNA rapidly and with high affinity; (2) has a good signal-to-noise ratio even at low concentrations; (3) exhibits a low photobleaching rate; and (4) induces lower light-induced DNA degradation. Finally, we show that SYTOX green is a DNA intercalator that binds DNA cooperatively with a binding site of 3.5 bp, increasing the DNA length upon binding by 43%, while not affecting its mechanical properties.

Enantiomeric discrimination of chiral organic salts by chiral aza-15-crown-5 ether with C 1 symmetry: experimental and theoretical approaches

The work involves an experimental ((1)H NMR) and theoretical (MD, MM-PBSA and DFT) investigation of the molecular recognition and discrimination properties of a chiral aza-15-crown-5 against methyl esters of alanine, phenylalanine and valine hydrochloride salts. The results indicate that the receptor binds enantiomers with moderate binding constants (88-1,389 M(-1)), with phenylalanine being more discriminated. The difference in experimental binding free energies (ΔG(R) - ΔG(S)) for alanine, phenylalanine and valine enantiomers were calculated as -0.36, -1.58 and 0.80 kcal mol(-1), respectively. The differences in theoretical binding energies were calculated by MM-PBSA (ΔE(R)PB - ΔE(S)PB=) as -0.30, -1.45 and 0.88, by B3LYP/6-31+G(d) as -1.17, -0.84 and 0.74 and by M06-2X/6-31+G(d) as -1.40, -3.26 and 1.66 kcal mol(-1). The data obtained give valuable information regarding the molecular recognition mode of the organoammonium complexes of chiral aza-crown ether with C 1 symmetry, which may be relevant to biological systems.

⁹⁹mTc pyrene derivative complex causes double-strand breaks in dsDNA mainly through cluster-mediated indirect effect in aqueous solution

Radiation therapy for cancer patients works by ionizing damage to nuclear DNA, primarily by creating double-strand breaks (DSB). A major shortcoming of traditional radiation therapy is the set of side effect associated with its long-range interaction with nearby tissues. Low-energy Auger electrons have the advantage of an extremely short effective range, minimizing damage to healthy tissue. Consequently, the isotope ^99mTc, an Auger electron source, is currently being studied for its beneficial potential in cancer treatment. We examined the dose effect of a pyrene derivative ^99mTc complex on plasmid DNA by using gel electrophoresis in both aqueous and methanol solutions. In aqueous solutions, the average yield per decay for double-strand breaks is 0.011±0.005 at low dose range, decreasing to 0.0005±0.0003 in the presence of 1 M dimethyl sulfoxide (DMSO). The apparent yield per decay for single-strand breaks (SSB) is 0.04±0.02, decreasing to approximately a fifth with 1 M DMSO. In methanol, the average yield per decay of DSB is 0.54±0.06 and drops to undetectable levels in 2 M DMSO. The SSB yield per decay is 7.2±0.2, changing to 0.4±0.2 in the presence of 2 M DMSO. The 95% decrease in the yield of DSB in DMSO indicates that the main mechanism for DSB formation is through indirect effect, possibly by cooperative binding or clustering of intercalators. In the presence of non-radioactive ligands at a near saturation concentration, where radioactive Tc compounds do not form large clusters, the yield of SSB stays the same while the yield of DSB decreases to the value in DMSO. DSBs generated by ^99mTc conjugated to intercalators are primarily caused by indirect effects through clustering.

Atomistic account of structural and dynamical changes induced by small binders in the double helix of a short DNA

Nucleic acids are flexible molecules and their dynamical properties play a key role in molecular recognition events. Small binders interacting with DNA fragments induce both structural and dynamical changes in the double helix. We study the dynamics of a DNA dodecamer and of its complexes with Hoechst 33258, which is a minor groove binder, and with the ethidium cation, which is an intercalator, by molecular dynamics simulation. The thermodynamics of DNA-drug interaction is evaluated in connection with the structure and the dynamics of the resulting complexes. We identify and characterize the relevant changes in the configurational distribution of the DNA helix and relate them to the corresponding entropic contributions to the binding free energy. The binder Hoechst locks the breathing motion of the minor groove inducing a reduction of the configurational entropy of the helix, which amounts to 20 kcal mol(-1). In contrast, intercalations with the ethidium cation enhance the flexibility of the double helix. We show that the balance between the energy required to deform the helix for the intercalation and the gain in configurational entropy is the origin of cooperativity in the binding of a second ethidium and of anti-cooperativity in the binding of a third one. The results of our study provide an understanding of the relation between structure, dynamics and energetics in the interaction between DNA fragments and small binders, highlighting the role of dynamical changes and consequent variation of the configurational entropy of the DNA double helix for both types of binders.

Drug-DNA intercalation: from discovery to the molecular mechanism

The ability of small molecules to perturb the natural structure and dynamics of nucleic acids is intriguing and has potential applications in cancer therapeutics. Intercalation is a special binding mode where the planar aromatic moiety of a small molecule is inserted between a pair of base pairs, causing structural changes in the DNA and leading to its functional arrest. Enormous progress has been made to understand the nature of the intercalation process since its idealistic conception five decades ago. However, the biological functions were detected even earlier. In this review, we focus mainly on the acridine and anthracycline types of drugs and provide a brief overview of the development in the field through various experimental methods that led to our present understanding of the subject. Subsequently, we discuss the molecular mechanism of the intercalation process, free-energy landscapes, and kinetics that was revealed recently through detailed and rigorous computational studies.

Solution structure and thermodynamics of 2',5' RNA intercalation

As a means to explore the influence of the nucleic acid backbone on the intercalative binding of ligands to DNA and RNA, we have determined the solution structure of a proflavine-bound 2',5'-linked octamer duplex with the sequence GCCGCGGC. This structure represents the first NMR structure of an intercalated RNA duplex, of either backbone structural isomer. By comparison with X-ray crystal structures, we have identified similarities and differences between intercalated 3',5' and 2',5'-linked RNA duplexes. First, the two forms of RNA have different sugar pucker geometries at the intercalated nucleotide steps, yet have the same interphosphate distances. Second, as in intercalated 3',5' RNA, the phosphate backbone angle zeta at the 2',5' RNA intercalation site prefers to be in the trans conformation, whereas unintercalated 2',5' and 3',5' RNA prefer the -gauche conformation. These observations provide new insights regarding the transitions required for intercalation of a phosphodiester-ribose backbone and suggest a possible contribution of the backbone to the origin of the nearest-neighbor exclusion principle. Thermodynamic studies presented for intercalation of both structural RNA isomers also reveal a surprising sensitivity of intercalator binding enthalpy and entropy to the details of RNA backbone structure.

Electrostatic contributions drive the interaction between Staphylococcus aureus protein Efb-C and its complement target C3d

The C3-inhibitory domain of Staphylococcus aureus extracellular fibrinogen-binding protein (Efb-C) defines a novel three-helix bundle motif that regulates complement activation. Previous crystallographic studies of Efb-C bound to its cognate subdomain of human C3 (C3d) identified Arg-131 and Asn-138 of Efb-C as key residues for its activity. In order to characterize more completely the physical and chemical driving forces behind this important interaction, we employed in this study a combination of structural, biophysical, and computational methods to analyze the interaction of C3d with Efb-C and the single-point mutants R131A and N138A. Our results show that while these mutations do not drastically affect the structure of the Efb-C/C3d recognition complex, they have significant adverse effects on both the thermodynamic and kinetic profiles of the resulting complexes. We also characterized other key interactions along the Efb-C/C3d binding interface and found an intricate network of salt bridges and hydrogen bonds that anchor Efb-C to C3d, resulting in its potent complement inhibitory properties.

Analysis of actinomycin D-DNA model complexes using a quantum-chemical criterion: Mulliken overlap populations

The binding of the antitumoral drug actinomycin D to single- and double-stranded DNA was investigated using molecular modeling in the frame of MM+ molecular mechanics and AM1 semi-empirical method. Two other programs, especially conceived to analyze hydrogen-bonding patterns in biological macromolecules, HBexplore, based on geometrical criteria and SHB_interactions, based on quantum-chemical criteria (Mulliken overlap populations), were also used. The results account for the non-cooperative intercalative binding process previously investigated, and outline the contribution of specific hydrogen bonding as well as CH...O(N) and other atom-atom intermolecular interactions to the stabilization of the actinomycin D-DNA complexes. They also support the hemi-intercalation model proposed in literature for the actinomycin D-ssDNA complex.

Structure of Rhodamine 6G−DNA Complexes from Molecular Dynamics Simulations

Chromophore-DNA complexes are useful for understanding charge transport along pi-stacks once their structural properties have been clarified. We studied two rhodamine 6G semicapping complexes with 15-mer B-DNA duplexes to determine the preferred orientation of the dye with respect to the neighboring base pair. For each of these systems, two distinct chromophore alignments were identified and quantified in terms of base-step parameters. The obtained geometries agree well with those derived from an NMR structure refinement of similar complexes. Cross-correlation analysis of the base-step parameters shows that slide and twist are highly interdependent during the structural transition from one conformation to the other.

Unprecedented dual binding behaviour of acridine group of dye: A combined experimental and theoretical investigation for the development of anticancer chemotherapeutic agents

Acridine group of dyes are well known in the field of development of probes for nucleic acid structure and conformational determination because of their relevance in the development of novel chemotherapeutic agents, footprinting agents and for gene manipulation in biotechnology and medicine. Here, we report the interaction of 9-N,N-dimethylaniline decahydroacridinedione (DMAADD), a new class of dye molecule with calf thymus DNA (CT-DNA) which has been studied extensively by means of traditional experimental and theoretical techniques. The changes in the base stacking of CT-DNA upon the binding of DMAADD are reflected in the circular dichroic (CD) spectral studies. Competitive binding study shows that the enhanced emission intensity of ethidium bromide (EB) in presence of DNA was quenched by the addition of DMAADD indicating that it displaces EB from its binding site in DNA and the apparent binding constant has been estimated to be (3.3+/-0.2)x10(5) M(-1). This competitive binding study and further fluorescence experiments reveal that DMAADD is a moderate binder of CT-DNA, while viscosity measurements show that the mode of binding is partial intercalation. Generally, one would expect increase in the melting temperature (T(m)) of DNA in presence of intercalators. Interestingly, an unusual decrease in melting temperature (DeltaT(m) of -4+/-0.2 degrees C) of DNA by the addition of DMAADD was observed. From our knowledge such a decreasing trend in melting point was not reported before for all the possible modes of binding. Molecular modeling gave the pictorial view of the binding model which clearly shows that of the various mode of binding, the dye prefers the major groove binding to the sites rich in GC residues and to the sites rich in AT residues it prefers intercalation mode of binding either through major or minor groove with the inclusion of the N,N-dimethylaniline (DMA) group inside the double helix which has been stacked in between the bases, under physiological relevant pH of 7.5.

Docking, triggering, and biological activity of dynemicin A in DNA: a computational study

The triggering and biological activity of the naturally occurring enediyne dynemicin A (1) was investigated, both inside and outside the minor groove of the duplex 10-mer B-DNA sequence d(CTACTACTGG).d(CCAGTAGTAG), using density functional theory (B3LYP with the 3-21G and 6-31G(d) basis set), BD(T)/cc-pVDZ (Brueckner doubles with a perturbative treatment of triple excitations), and the ONIOM approach. Enediyne 1 is triggered by NADPH in a strongly exothermic reaction (-88 kcal/mol), which involves a number of intermediate steps. Untriggered 1 has a high barrier for the Bergman cyclization (52 kcal/mol) that is lowered after triggering to 16.7 kcal/mol due to an epoxide opening and the accompanying strain relief. The Bergman reaction of triggered 1 is slightly exothermic by 2.8 kcal/mol. The singlet biradical formed in this reaction is kinetically stable (activation enthalpies of 19.5 and 21.8 kcal/mol for retro-Bergman reactions) and is as reactive as para-benzyne. The activity-relevant docking mode is an edge-on insertion into the minor groove, whereas the intercalation between base pairs, although leading to larger binding energies, excludes a triggering of 1 and the development of its biological activity. Therefore, an insertion-intercalation model is developed, which can explain all known experimental observations made for 1. On the basis of the insertion-intercalation model it is explained why large intercalation energies suppress the biological activity of dynemicin and why double-strand scission can be achieved only in a two-step mechanism that involves two enediyne molecules, explaining thus the high ratio of single-strand to double-strand scission observed for 1.