DNA-drug interactions. The crystal structures of d(TGTACA) and d(TGATCA) complexed with daunomycin

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.

Anti-Inflammatory and Antioxidant Pyrrolo[3,4-d]pyridazinone Derivatives Interact with DNA and Bind to Plasma Proteins-Spectroscopic and In Silico Studies

From the point of view of the search for new pharmaceuticals, pyridazinone derivatives are a very promising group of compounds. In our previous works, we have proved that newly synthesized ligands from this group have desirable biological and pharmacokinetic properties. Therefore, we decided to continue the research evaluating the activity of pyrrolo[3,4-dpyridazinone derivatives. In this work, we focused on the interactions of five pyridazinone derivatives with the following biomolecules: DNA and two plasma proteins: orosomucoid and gamma globulin. Using several of spectroscopic methods, such as UV-Vis, CD, and fluorescence spectroscopy, we proved that the tested compounds form stable complexes with all biomacromolecules selected for analysis. These findings were also confirmed by the results obtained by molecular modeling. All tested pyridazinone derivatives bind to the ctDNA molecule via groove binding mechanisms. All these molecules can also be bound and transported by the tested plasma proteins; however, the stability of the complexes formed is lower than those formed with serum albumin.

Doxorubicin and other anthracyclines in cancers: Activity, chemoresistance and its overcoming

Anthracyclines have been important and effective treatments against a number of cancers since their discovery. However, their use in therapy has been complicated by severe side effects and toxicity that occur during or after treatment, including cardiotoxicity. The mode of action of anthracyclines is complex, with several mechanisms proposed. It is possible that their high toxicity is due to the large set of processes involved in anthracycline action. The development of resistance is a major barrier to successful treatment when using anthracyclines. This resistance is based on a series of mechanisms that have been studied and addressed in recent years. This work provides an overview of the anthracyclines used in cancer therapy. It discusses their mechanisms of activity, toxicity, and chemoresistance, as well as the approaches used to improve their activity, decrease their toxicity, and overcome resistance.

Preparation and Evaluation of Thermosensitive Liposomes Encapsulating I-125-Labeled Doxorubicin Derivatives for Auger Electron Therapy

Auger electrons (AEs) are very low-energy electrons emitted by radionuclides such as I-125 (¹²⁵I). This energy is deposited across a small distance (<0.5 μm), resulting in high linear energy transfer that is potent for causing lethal damage to cancer cells. Thus, AE-emitting radiotherapeutic agents have great potential for cancer treatment. In this study, thermosensitive liposomes (TSLs) encapsulating ¹²⁵I-labeled doxorubicin (DOX) derivatives were developed for Auger electron therapy, targeting the DNA of cancer cells. A radioiodinated DOX derivative [¹²⁵I]5 highly accumulated in the nuclei of cancer cells and showed potent cytotoxicity against Colon 26 cancer cells by AEs. Subsequently, [¹²⁵I]5 was loaded into the TSLs with high encapsulation efficiency. Potent release of [¹²⁵I]5 from TSLs was achieved with heating, whereas a decreased release was observed without heating. Furthermore, TSLs encapsulating [¹²⁵I]5 showed a high uptake in the nuclei at 42 °C for 1 h. We supposed that [¹²⁵I]5 was released by heating at 42 °C and accumulated in the nuclei in the cells. These results suggest that the combination of TSLs encapsulating [¹²⁵I]5 and hyperthermia is an effective cancer therapy.

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.

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.

Four steps for revealing and adjusting the 3D structure of aptamers in solution by small-angle X-ray scattering and computer simulation

Nucleic acid (NA) aptamers bind to their targets with high affinity and selectivity. The three-dimensional (3D) structures of aptamers play a major role in these non-covalent interactions. Here, we use a four-step approach to determine a true 3D structure of aptamers in solution using small-angle X-ray scattering (SAXS) and molecular structure restoration (MSR). The approach consists of (i) acquiring SAXS experimental data of an aptamer in solution, (ii) building a spatial distribution of the molecule's electron density using SAXS results, (iii) constructing a 3D model of the aptamer from its nucleotide primary sequence and secondary structure, and (iv) comparing and refining the modeled 3D structures with the experimental SAXS model. In the proof-of-principle we analyzed the 3D structure of RE31 aptamer to thrombin in a native free state at different temperatures and validated it by circular dichroism (CD). The resulting 3D structure of RE31 has the most energetically favorable conformation and the same elements such as a B-form duplex, non-complementary region, and two G-quartets which were previously reported by X-ray diffraction (XRD) from a single crystal. More broadly, this study demonstrates the complementary approach for constructing and adjusting the 3D structures of aptamers, DNAzymes, and ribozymes in solution, and could supply new opportunities for developing functional nucleic acids. Graphical abstract.

DNA-binding mechanism of spiropyran photoswitches: the role of electrostatics

The binding mechanism of the protonated open form of three spiropyran derivatives into a 12-mer (poly-dAT)₂ has been unveiled by means of computational methods.

The binding mechanism of the protonated open form of three spiropyran derivatives into a 12-mer (poly-dAT)₂ has been unveiled by means of computational methods. It is found that the electrostatic term in the probe:DNA binding energy, modulates the binding mode, providing new guidelines for the design of spiropyran photoswitches with specific binding modes to DNA.

Molecular mechanism and binding free energy of doxorubicin intercalation in DNA

The intercalation process of binding doxorubicin (DOX) in DNA is studied by extensive MD simulations. Many molecular factors that control the binding affinity of DOX to DNA to form a stable complex are inspected and quantified by employing continuum solvation models for estimating the binding free energy. The modified MM-PB(GB)SA methodology provides a complete energetic profile of ΔGele, ΔGvDW, ΔGpolar, ΔGnon-polar, TΔStotal, ΔGdeform, ΔGcon, and ΔGion. To identify the sequence specificity of DOX, two different DNA sequences, d(CGATCG) or DNA1 and d(CGTACG) or DNA2, with one molecule (1 : 1 complex) or two molecule (2 : 1 complex) configurations of DOX were selected in this study. Our results show that the DNA deformation energy (ΔGdeform), the energy cost from translational and rotational entropic contributions (TΔStran+rot), the total electrostatic interactions (ΔGpolar-PB/GB + ΔGele) of incorporation, the intramolecular electrostatic interactions (ΔGele) and electrostatic polar solvation interactions (ΔGpolar-PB/GB) are all unfavorable to the binding of DOX to DNA. However, they are overcome by at least five favorable interactions: the van der Waals interactions (ΔGvDW), the non-polar solvation interaction (ΔGnon-polar), the vibrational entropic contribution (TΔSvib), and the standard concentration dependent free energies of DOX (ΔGcon) and the ionic solution (ΔGion). Specifically, the van der Waals interaction appears to be the major driving force to form a stable DOX-DNA complex. We also predict that DOX has stronger binding to DNA1 than DNA2. The DNA deformation penalty and entropy cost in the 2 : 1 complex are less than those in the 1 : 1 complex, thus they indicate that the 2 : 1 complex is more stable than the 1 : 1 complex. We have calculated the total binding free energy (BFE) (ΔGt-sim) using both MM-PBSA and MM-GBSA methods, which suggests a more stable DOX-DNA complex at lower ionic concentration. The calculated BFE from the modified MM-GBSA method for DOX-DNA1 and DOX-DNA2 in the 1 : 1 complex is -9.1 and -5.1 kcal mol-1 respectively. The same quantities from the modified MM-PBSA method are -12.74 and -8.35 kcal mol-1 respectively. The value of the total BFE ΔGt-sim in the 1 : 1 complex is in reasonable agreement with the experimental value of -7.7 ± 0.3 kcal mol-1.

Structure-based DNA-targeting strategies with small molecule ligands for drug discovery

Nucleic acids are the molecular targets of many clinical anticancer drugs. However, compared with proteins, nucleic acids have traditionally attracted much less attention as drug targets in structure-based drug design, partially because limited structural information of nucleic acids complexed with potential drugs is available. Over the past several years, enormous progresses in nucleic acid crystallization, heavy-atom derivatization, phasing, and structural biology have been made. Many complicated nucleic acid structures have been determined, providing new insights into the molecular functions and interactions of nucleic acids, especially DNAs complexed with small molecule ligands. Thus, opportunities have been created to further discover nucleic acid-targeting drugs for disease treatments. This review focuses on the structure studies of DNAs complexed with small molecule ligands for discovering lead compounds, drug candidates, and/or therapeutics.

Early stage intercalation of doxorubicin to DNA fragments observed in molecular dynamics binding simulations

The intercalation mode between doxorubicin (an anticancer drug) and two 6-base-pair DNA model fragments (d(CGATCG)₂ and d(CGTACG)₂) has been well studied by X-ray crystallography and NMR experimental methods. Yet, the detailed intercalation pathway at molecular level remains elusive. In this study, we conducted molecular dynamics binding simulations of these two systems using AMBER DNA (parmbsc0) and drug (GAFF) force fields starting from the unbound state. We observed outside binding (minor groove binding or end-binding) in all six independent binding simulations (three for each DNA fragment), followed by the complete intercalation of a drug molecule in two simulations (one for each DNA fragment). First, our data directly supported that the minor groove binding is the dominant pre-intercalation step. Second, we observed that the opening and flipping of a local base pair (A3-T10 for d(CGATCG)₂ and C1-G12 for d(CGTACG)₂) in the two intercalation trajectories. This locally cooperative flipping-intercalation mechanism was different from the previously proposed rise-insertion mechanism by which the distance between two neighboring intact base pairs increases to create a space for the drug insertion. Third, our simulations provided the first set of data to support the applicability of the AMBER DNA and drug force fields in drug-DNA atomistic binding simulations. Implications on the kinetics pathway and drug action are also discussed.

A novel sensitive DNA-biosensor for detection of a carcinogen, Sudan II, using electrochemically treated pencil graphite electrode by voltammetric methods

A simple and inexpensive methodology was used to develop a novel electrochemical sensor for the determination of Sudan II. The interaction of Sudan II with salmon sperm ds-DNA on the surface of salmon sperm ds-DNA-modified pencil graphite electrode (PGE) and in solution phase was studied, using differential pulse voltammetry. The difference between adenine and guanine signals of the ds-DNA after and before interaction with Sudan II was directly proportional to Sudan II concentration, which used for quantitative inspections. Using PGE, a linear calibration curve (R(2)=0.9958) was observed with 0.5-6.0 μg mL(-1) Sudan II. Furthermore, the LOD of 0.4 μg mL(-1) and linear range between 0.5 and 4.0 μg mL(-1) were achieved in solution phase. In the second part, Sudan II was determined on a pretreated pencil graphite electrode by means of adsorptive stripping differential pulse voltammetry. The peak current was linearly dependent on Sudan II concentration over the range of 0.0015-0.30 μg mL(-1), with a detection limit of 0.00007 μg mL(-1) Sudan II. Both ds-DNA-modified PGE and PPGE were applied to analyze Sudan II in real samples.

Quinoline alkaloids as intercalative topoisomerase inhibitors

Quinoline alkaloids are abundant in the Rutaceae, and many have exhibited cytotoxic activity. Because structurally related antitumor alkaloids such as camptothecin and fagaronine are known to function as intercalative topoisomerase poisons, it is hypothesized that cytotoxic Stauranthus alkaloids may also serve as intercalative topoisomerase inhibitors. To test this hypothesis theoretically, ten Stauranthus quinoline alkaloids were examined for potential intercalation into DNA using a molecular docking approach. Four of the alkaloids (stauranthine, skimmianine, 3',6'-dihydroxy-3',6'-dihydrostauranthine, and trans-3',4'-dihydroxy-3',4'-dihydrostauranthine) were able to intercalatively dock consistently into DNA. In order to probe the intermolecular interactions that may be responsible for intercalation of these quinoline alkaloids, density functional calculations have been carried out using both the B3LYP and M06 functionals. M06 calculations indicated favorable pi-pi interactions between either skimmianine or stauranthine and the guanine-cytosine base pair. Furthermore, the lowest-energy face-to-face orientation of stauranthine with guanine is consistent with favorable dipole-dipole orientations, favorable electrostatic interactions, and favorable frontier molecular orbital interactions. Likewise, the lowest-energy face-to-face orientation of stauranthine with the guanine-cytosine base pair reveals favorable electrostatic interactions as well as frontier molecular orbital interactions. Thus, not only can quinoline alkaloids dock intercalatively into DNA, but the docked orientations are also electronically favorable.

Studies on drug-DNA complexes, adriamycin-d-(TGATCA)(2) and 4'-epiadriamycin-d-(CGATCG)(2), by phosphorus-31 nuclear magnetic resonance spectroscopy

The complexes of adriamycin-d-(TGATCA)(2) and 4'-epiadriamycin-d-(CGATCG)(2) are studied by one- and two-dimensional (31)P nuclear magnetic resonance spectroscopy (NMR) at 500 MHz in the temperature range 275-328 K and as a function of drug to DNA ratio (0.0-2.0). The binding of drug to DNA is clearly evident in (31)P-(31)P exchange NOESY spectra that shows two sets of resonances in slow chemical exchange. The phosphate resonances at the intercalating steps, T1pG2/C1pG2 and C5pA6/C5pG6, shift downfield up to 1.7 ppm and that at the adjacent step shift downfield up to 0.7 ppm, whereas the central phosphate A3pT4 is relatively unaffected. The variations of chemical shift with drug to DNA ratio and temperature as well as linewidths are different in each of the two complexes. These observations reflect change in population of B(I)/B(II) conformation, stretching of backbone torsional angle zeta, and distortions in O-P-O bond angles that occur on binding of drug to DNA. To the best of our knowledge, there are no solution studies on 4'-epiadriamycin, a better tolerated drug, and binding of daunomycin or its analogue to d-(TGATCA)(2) hexamer sequence. The studies report the use of (31)P NMR as a tool to differentiate various complexes. The specific differences may well be the reasons that are responsible for different antitumor action of these drugs due to different binding ability and distortions in DNA.

Solution studies on the complex of 4'-epiadriamycin-d-(CGATCG)2 followed by time-resolved fluorescence measurement, diffusion ordered spectroscopy and restrained molecular dynamics simulations

4'-Epiadriamycin is a better-tolerated anthracycline drug, due to lesser cardiotoxicity. We report here a study of the 2:1 complex of 4'-epiadriamycin-d-(CGATCG)(2) by proton Nuclear Magnetic Resonance Spectroscopy which show the absence of sequential connectivities between C1pG2 and C5pG6 base pair steps and presence of intermolecular cross peaks of the drug and DNA. Our studies establish the role of 9OH, NH3+, 7O, 4OCH(3) groups in binding to DNA. Time-resolved fluorescence measurement and diffusion ordered spectroscopic studies reveal the formation of complex. The nonspecific interactions as well as those essential for biological activity are discussed along with its medicinal importance.

Volume and hydration changes of DNA–ligand interactions

We report the volumetric and other thermodynamic properties of ethidium bromide (EB), propidium iodide (PI) and daunomycin (DAU) intercalating with poly(dA).poly(dT), poly[d(A-T)].poly[d(A-T)], and poly[d(G-C)].poly[d(G-C)], respectively, as well as minor groove binder Hoechst 33258 binding with poly[d(A-T)].poly[d(A-T)]. The data were obtained using fluorescence titration and hydrostatic pressure measurements. Our thermodynamic data are combined with enthalpies from literature reports to analyze the thermodynamic characteristics of the different interactions. The differences are interpreted based on three processes related to hydration: I. burial of non-polar hydrophobic solvent accessible surface, II. burial of polar surface and formation of solute-solute H-bonds, and III. disruption of "structural" hydration. Sequence dependent conformational changes may also be important when comparing ligand binding to different DNA sequences. We conclude that a combination of different thermodynamic parameters, especially volume change, is essential in order to understand the role of hydration in the energetics of DNA-ligand interactions.

Intercalation of daunomycin into d(CG)4 oligomer duplex containing G x T mismatches by vibrational circular dichroism and infrared absorption spectroscopy

The vibrational circular dichroism (VCD) and infrared absorption (IR) spectra of the mismatched octamer oligonucleotides d(CGTGCGCG)(2) (CGT) and d(CGCGTGCG)(2) (CGC) and their complexes with the antitumor drug daunomycin were measured in D(2)O, interpreted, and compared to the octamer d(CGCGCGCG)(2) (CG). The IR spectra of the mismatched octamers in the carbonyl-stretching region are similar to those of the parent CG, whereas the VCD spectra differ in several respects between each other. The main VCD feature due to carbonyl stretching is informative for the mismatches and CG. Vibrational modes in the sugar-phosphate region remain essentially unchanged especially for PO(2) (-) symmetric stretching. Differences between the free and complexed mismatch octamers occurred mainly in the carbonyl-stretching region (1,700-1,600 cm(-1)). The absorption intensity of the C==O peak of G is more prominent for CGC than CGT and resembles CG in this respect. The detailed composition of this doublet is clearly visible, indicating the geometric rearrangement of the base pairs in the presence of the mismatch and upon forming the daunomycin complex.

Restrained molecular dynamics studies on complex of adriamycin with DNA hexamer sequence d-CGATCG

Adriamycin is an anthracycline anticancer drug used widely for solid tumors in spite of its adverse side effects. The solution structure of 2:1 adriamycin-d-(CGATCG)(2) complex has been studied by restrained molecular dynamics simulations. The restraint data set consists of several intramolecular and intermolecular nuclear Overhauser enhancement cross-peaks obtained from two-dimensional nuclear magnetic resonance spectroscopy data. The drug is found to intercalate between CG and GC base pairs at two d-CpG sites. The drug-DNA complex is stabilized via specific hydrogen bonding and van der Waal's interactions involving 4OCH(3), O5, 6OH, and NH(3)(+) moiety of daunosamine sugar, and rings A protons. The O-glycosidic bond C7-O7-C1'-C2' lies in the range 138 degrees -160 degrees during the course of simulations. The O6-H6...O5 hydrogen bond is stable while O11-H11...O12 hydrogen bond is not favored. The intercalating base pairs are buckled and minor groove is wider in the complex. The phosphate on one strand at intercalation site C1pG2 is in B(I) conformation and the phosphates directly lying on opposite strand is in B(II) conformation. The phosphorus on adjacent site G2pA3 is in B(II) conformation and hence a distinct pattern of B(I) and B(II) conformations is induced and stabilized. The role of various functional groups by which the molecular action is mediated has been discussed and correlated to the available biochemical evidence.

Structure of DNA sequence d-TGATCA by two-dimensional nuclear magnetic resonance spectroscopy and restrained molecular dynamics

The 5' d-TpG 3' element is a part of DNA sequences involved in regulation of gene expression and is also a site for intercalation of several anticancer drugs. Solution conformation of DNA duplex d-TGATCA containing this element has been investigated by two-dimensional NMR spectroscopy. Using a total of 12 torsional angles and 121 distance constraints, structural refinement has been carried out by restrained molecular dynamics (rMDs) in vacuum up to 100 ps. The structure is characterized by a large positive roll at TpG/CpA base pair step and large negative propeller twist for AT and TA base pairs. The backbone torsional angle, gamma(O5'-C5'-C4'-C3'), of T1 residue adopts a trans-conformation which is corroborated by short intra nucleotide T1H6-T1H5' (3.7A) distance in nuclear overhauser effect spectroscopy (NOESY) spectra while the backbone torsional angle, beta(P-O5'-C5'-C4'), exists in trans as well as gauche state for T1 and C5 residues. There is evidence of significant flexibility of the sugar-phosphate backbone with rapid inter-conversion between two different conformers at TpG/CpA base pair step. The base sequence dependent variations and local structural heterogeneity have important implications in specific recognition of DNA by ligands.

Sequence selective binding of bis-daunorubicin WP631 to DNA

We have used footprinting techniques on a wide range of natural and synthetic footprinting substrates to examine the sequence-selective interaction of the bis-daunorubicin antibiotic WP631 with DNA. The ligand produces clear DNase I footprints that are very different from those seen with other anthracycline antibiotics such as daunorubicin and nogalamycin. Footprints are found in a diverse range of sequences, many of which are rich in GT (AC) or GA (TC) residues. As expected, the ligand binds well to the sequences CGTACG and CGATCG, but clear footprints are also found at hexanucleotide sequences such GCATGC and GCTAGC. The various footprints do not contain any particular unique di-, tri- or tetranucleotide sequences, but are frequently contain the sequence (G/C)(A/T)(A/T)(G/C). All sequences with this composition are protected by the ligand, though it can also bind to some sites that differ from this consensus by one base pair.

Electrochemical genosensor for Mitomycin C–DNA interaction based on guanine signal

The interaction of mitomycin C (MC) with fish sperm or calf thymus DNA immobilized onto carbon screen-printed electrodes (CSPE) and carbon paste electrode (CPE) have been studied by using electrochemical techniques as square wave voltammetry (SWV) and differential pulse voltammetry (DPV). After the interaction was occurred between DNA and MC on electrode surface, it was observed that the guanine signal was higher with bare electrode than DNA-modified one. The changes in the experimental parameters such as the concentration of MC, and the accumulation time of MC were studied by using SWV and DPV. In addition, reproducibility, and detection limit parameters were determined using both electrodes. The partition coefficient of MC was also calculated before and after interaction of MC with dsDNA at CPE surface. These results showed that these two different DNA biosensors could be used for the sensitive, rapid and cost effective detection of MC-DNA interaction.

Daunomycin Intercalation Stabilizes Distinct Backbone Conformations of DNA

Daunomycin is a widely used antibiotic of the anthracycline family. In the present study we reveal the structural properties and important intercalator-DNA interactions by means of molecular dynamics. As most of the X-ray structures of DNA-daunomycin intercalated complexes are short hexamers or octamers of DNA with two drug molecules per doublehelix we calculated a self complementary 14-mer oligodeoxyribonucleotide duplex d(CGCGCGATCGCGCG)2 in the B-form with two putative intercalation sites at the 5'-CGA-3' step on both strands. Consequently we are able to look at the structure of a 1:1 complex and exclude crystal packing effects normally encountered in most of the X-ray crystallographic studies conducted so far. We performed different 10 to 20 ns long molecular dynamics simulations of the uncomplexed DNA structure, the DNA-daunomycin complex and a 1:2 complex of DNA-daunomycin where the two intercalator molecules are stacked into the two opposing 5'-CGA-3' steps. Thereby--in contrast to X-ray structures--a comparison of a complex of only one with a complex of two intercalators per doublehelix is possible. The chromophore of daunomycin is intercalated between the 5'-CG-3' bases while the daunosamine sugar moiety is placed in the minor groove. We observe a flexibility of the dihedral angle at the glycosidic bond, leading to three different positions of the ammonium group responsible for important contacts in the minor groove. Furthermore a distinct pattern of BI and BII around the intercalation site is induced and stabilized. This indicates a transfer of changes in the DNA geometry caused by intercalation to the DNA backbone.

Disposable electrochemical biosensor for the detection of the interaction between DNA and lycorine based on guanine and adenine signals

The interaction of lycorine (LYC) with calf thymus double stranded DNA (dsDNA) and calf thymus single stranded DNA (ssDNA) was studied electrochemically based on the oxidation signals of guanine and adenine by using differential pulse voltammetry (DPV) at carbon paste electrode (CPE) and pencil graphite electrode (PGE). As a result of the interaction of LYC with DNA, the voltammetric signal of guanine and adenine greatly decreased. The changes in the experimental parameters such as the concentration of LYC, and the accumulation time of LYC were studied by using DPV with PGE. The interaction of LYC with synthetic polynucleotides, such as poly[G] was also observed. The interaction of LYC with dsDNA was also observed at PGE in solution phase. In addition, the detection limit and the reproducibility was determined by using both electrochemical transducers. The application of electrochemical methods on the interactions between DNA and DNA targeted agent were explored.

Interaction of daunomycin antibiotic with histone H(1): ultraviolet spectroscopy and equilibrium dialysis studies

Using ultraviolet spectroscopy and equilibrium dialysis techniques, we have investigated the interaction of anticancer drug, daunomycin with calf thymus histone H(1) chromosomal protein in 20 mM phosphate buffer, pH 7.0, 1 mM EDTA at room temperature. The UV spectroscopy results show that daunomycin (5.0-100 microM) decreases the absorbance of histone H(1) at 210-230 nm and induces hypochromicity in the absorption spectrum of the protein. The equilibrium dialysis data show that daunomycin binds to histone H(1) and the binding process is positive cooperative with two binding sites as Scatchard plot and Hill coefficient confirm it. The results suggest that daunomycin binds to histone H(1) and changes its conformation.

The effects of daunomycin antibiotic on histone H(1): thermal denaturation and fluorescence spectroscopy studies

Using thermal denaturation and fluorescence spectroscopy, we have investigated the interaction of antitumor antibiotic, daunomycin, with calf thymus histone H(1) under several ionic strengths. The results show that daunomycin binds to histone H(1) and increases its melting temperature. Increasing ionic strength elevates this effect. Fluorescence emission data show that the interaction of daunomycin with histone H(1) decreases the emission intensity at 325 nm and induces hyperchromicity in the emission spectrum of the drug. The results suggest that histone H(1) can be considered as a new target for drug action at the chromatin level.

NMR investigation of the complexation of daunomycin with deoxytetranucleotides of different base sequence in aqueous solution

500 MHz NMR spectroscopy has been used to investigate the complexation of the anthracycline antibiotic daunomycin (DAU) with self-complementary deoxytetranucleotides, 5'-d(CGCG), 5'-d(GCGC), 5'-d(TGCA), 5'-d(ACGT) and 5'-d(AGCT), of different base sequence in aqueous salt solution. 2D homonuclear 1H NMR spectroscopy (TOCSY and NOESY) and heteronuclear 1H - 31P NMR spectroscopy (HMBC) have been used for complete assignment of the non-exchangeable protons and the phosphorus resonance signals, respectively, and for a qualitative determination of the preferred binding sites of the drug. Analysis shows that DAU intercalates preferentially into the terminal sites of each of the tetranucleotides and that the aminosugar of the antibiotic is situated in the minor groove of the tetramer duplex, partly eclipsing the third base pair. A quantitative determination of the complexation of DAU with the deoxytetranucleotides has been made using the experimental concentration and temperature dependences of the drug proton chemical shifts; these have been analysed in terms of the equilibrium reaction constants, limiting proton chemical shifts and thermodynamical parameters (enthalpies deltaH, entropies deltaS) of different drug-DNA complexes (1:1, 1:2, 2:1, 2:2) in aqueous solution. It is found that DAU interacts with sites containing three adjacent base pairs but does not show any significant sequence specificity of binding with either single or double-stranded tetranucleotides, in contrast with other intercalating drugs such as proflavine, ethidium bromide and actinomycin D. The most favourable structures of the 1:2 complexes have been derived from the induced limiting proton chemical shifts of the drug in the intercalated complexes with the tetranucleotide duplex, in conjunction with 2D NOE data. It has been found that the conformational parameters of the double helix and the orientation of the DAU chromophore in the intercalated complexes depend on base sequence at the binding site of the tetramer duplexes in aqueous solution.

A continuous transition from A-DNA to B-DNA in the 1:1 complex between nogalamycin and the hexamer dCCCGGG

The antibiotic nogalamycin, a drug with high specificity for TG and CG steps in double-stranded DNA, has been crystallized as a 1:1 complex with the hexamer d(CCCGGG). The antibiotic is inserted at the central CG step of the duplex, with the two sugars oriented in the same direction and with strong interactions with the DNA within the grooves. The amino-glucose residue makes an integral part of a well defined major groove hydration network with van der Waals contacts and several strong hydrogen bonds to the duplex. The nogalose residue resides in the minor groove, making primarily van der Waals contacts. The single site allows an accurate molecular description of the intercalation, without perturbations from end effects observed previously. The local unwinding induced by nogalamycin is completely relaxed 2 base pairs away from the intercalation site. The two strands of the DNA show a continuous deformation from the A to the B form: 1) the cytosines toward the 5' end of the nogalomycin site in each strand have c3'-endo conformations while 5 guanosines toward the 3' ends have c2'-endo conformations; 2) within each strand, the phosphate-phosphate distances increase in a continuous manner from 5.7 A (A-form) to 7.1 A (B-form).

A proton nuclear magnetic resonance investigation of the conformation of daunomycin

The 500 MHz proton NMR spectra of 4.95 mM daunomycin in D2O have been investigated in the temperature range 277-350 K. Down field shifts of approximately 0.15 to 0.20 ppm in 1H, 2H and 3H protons with increasing temperature indicate that daunomycin exists in aggregated form at 277 K which is stabilized by stacking of aromatic rings. The 7H, 10axH and 10eqH protons show a change in chemical shift of 0.12 to 0.16 ppm, while other ring A/sugar protons shift by 0.0 to 0.08 ppm due to self association. The drug exists in monomer state at about 350 K. The conformational features have been ascertained from NOESY spectra at 297 K. The 7H proton is found to be strongly coupled to 8axH (J = 5 Hz) as compared to 8eqH (J = 2 Hz), while the 4'H-5'H connectivity is not observed in the COSY spectra. Besides the NOESY cross peaks between the spin-spin coupled protons, several intramolecular NOEs are seen. The 8axH and 8eqH are equally distant from 5'H proton. The distances of 3'H and 4'H daunosamine sugar protons from the ring A protons--7H, 8eqH, 9COCH3--are in the range 2.9-3.1 A, giving moderate cross peaks in NOESY spectra. The observed results imply the existence of predominantly 9H8 half chair conformation of ring A in aqueous solution, which is marginally different from that obtained by X-ray crystal analysis.

Binding of biologically important molecules to DNA, probed using electro-fluorescence polarization spectroscopy

When small fluorescent tag molecules are bound with a single geometry to larger carrier molecules the fluorescence excited using polarized light may itself be partially polarized. If this is the case then alignment of the carrier molecules changes the fluorescence polarization. Electro-fluorescence polarization spectrometry (EFPS) makes use of the fluorescence polarization changes which occur when the carrier molecules are aligned by an electric field. EFPS allows the determination of the binding geometry of the small tag molecules to the larger carrier molecules. In this study, EFPS has been used to probe the interactions between DNA and a series of chromosomal stains and anti-neoplastic (anti-cancer) agents. Results are presented for calf thymus DNA, treated with the dyes acridine orange, ethidium bromide, proflavine, Hoechst 33258, Hoechst 33342, 4',6-diamidino-2-phenylindole and the anti-tumour agents, doxorubicin, daunomycin, actinomycin C and actinomycin V.

Molecular structure of the halogenated anti-cancer drug iododoxorubicin complexed with d(TGTACA) and d(CGATCG)

4'-Deoxy-4'-iododoxorubicin, a halogenated anthracycline derivative, is an anticancer agent currently under Phase II clinical trials. In preclinical studies, it has demonstrated significantly reduced levels of cardiotoxicity compared to currently employed anthracyclines. It also has modified pharmacological properties resulting in an altered spectrum of experimental antitumor activity. The iodine atom at the 4' position of the sugar ring reduces the basicity and enhances the lipophilicity of this compound as compared to related anthracycline drugs. We report here single crystal X-ray diffraction studies of the complexes of 4'-deoxy-4'-iododoxorubicin with the hexanucleotide duplex sequences d(TGTACA) and d(CGATCG) at 1.6 and 1.5 A, respectively. The iodine substituent does not alter the geometry of intercalation as compared to previously solved anthracycline complexes, but appears to markedly affect the solvent environment of the structures. This could have consequences for the interaction of this drug with DNA and DNA binding proteins in cells.

Hedamycin intercalates the DNA helix and, through carbohydrate-mediated recognition in the minor groove, directs N7-alkylation of guanine in the major groove in a sequence-specific manner

BACKGROUND

The pluramycins are a class of antitumor antibiotics that exert their biological activity through interaction with DNA. Recent studies with the analog altromycin B have determined that these agents intercalate into the DNA molecule, position carbohydrate substituents into both major and minor grooves, and alkylate the DNA molecule by epoxide-mediated electrophilic attack on N7 of guanine located to the 3' side of the drug molecule. Alkylation is sequence dependent and appears to be modulated by glycoside substituents attached at the corners of a planar chromophore. The altromycin B-like analogs preferentially alkylate 5'AG sequences; hedamycin-like analogs prefer 5'TG and 5'CG sequences. Although the mechanism of guanine modification by altromycin B has been extensively studied, the mechanism of action of hedamycin has not been previously determined.

RESULTS

Using high-field NMR, we have shown that hedamycin stacks to the 5' side of the guanine nucleotide at the site of intercalation in a DNA decamer, positioning both aminosaccharides into the minor groove to direct alkylation by the epoxide moiety on N7 of guanine. The C10 linked N,N-dimethylvancosamine sugar moiety interacts to the 5' side of the intercalation site, while the C8 linked anglosamine moiety interacts to the 3' side. The binding interactions of the two aminosugars steer the C2 double epoxide located in the major groove into the proximity of N7 of guanine. Unexpectedly, it is not the first epoxide that undergoes electrophilic addition to N7 of guanine, which would correspond to altromycin B, but the second, terminal epoxide.

CONCLUSIONS

We have used two-dimensional NMR to elucidate the sequence-selective recognition of DNA by hedamycin and the mechanism of covalent modification of guanine by this antibiotic. Characterization of the intermolecular interactions between both hedamycin and altromycin B and their targeted DNA sequences has yielded a better understanding of the reasons for variations in sequence selectivity and alkylation reactivity among the pluramycin compounds.