A theoretical investigation on the sequence selective binding of daunomycin to double-stranded polynucleotides

Mechanism of Daunomycin Intercalation into DNA from Enhanced Sampling Simulations

Daunomycin is a widely used anticancer drug, yet the mechanism underlying how it binds to DNA remains contested. 469 all-atom trajectories of daunomycin binding to the DNA oligonucleotide d(GCG CAC GTG CGC) were collected using weighted ensemble (WE)-enhanced sampling. Mechanistic insights were revealed through analysis of the ensemble of trajectories. Initially, the binding process involves a ubiquitous hydrogen bond between the DNA backbone and the NH3+ group on daunomycin. During the binding process, most trajectories exhibited similar structural changes to DNA, including DNA base pair rise, bending, and minor groove width changes. Variability within the ensemble of binding trajectories illuminates differences in the orientation of daunomycin as it initially intercalates; around 10% of trajectories needed minimal rearrangement from intercalation to reaching the fully bound configuration, whereas most needed an additional 1-5 ns to rearrange. The results here emphasize the utility of generating an ensemble of trajectories to discern biomolecular binding mechanisms.

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.

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.

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.

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.

The anthracycline antibiotics: antitumor drugs that alter chromatin structure

Anthracycline antibiotics are an important group of antitumor drugs widely used in cancer chemotherapy. However, despite the increasing interest in these chemotherapeutic agents, their mechanism of action is not yet completely understood. Here, we review what is currently known about the molecular mechanisms involved with special emphasis on the interaction of these drugs with chromatin and its constitutive components: DNA and histones. The evidence suggests that one very important component of the activity of these drugs is the result of these manifold interactions that lead to a chromatin unfolding and aggregation. This chromatin structural disruption is likely to interfere with the metabolic processes of DNA (replication and transcription) and it may play an important role in the apoptosis undergone by the cells upon treatment with these drugs.

Electrostatic factors in DNA intercalation

The factors that determine the binding of a chromophore between the base pairs in DNA intercalation complexes are dissected. The electrostatic potential in the intercalation plane is calculated using an accurate ab initio based distributed multipole electrostatic model for a range of intercalation sites, involving different sequences of base pairs and relative twist angles. There will be a significant electrostatic contribution to the binding energy for chromophores with a predominantly positive electrostatic potential, but this varies significantly with sequence, and somewhat with twist angle. The usefulness of these potential maps for understanding the binding of intercalators is explored by calculating the electrostatic binding energy for 9-aminoacridine, ethidium, and daunomycin in a variety of model binding sites. The electrostatic forces play a major role in the positioning of an intercalating 9-aminoacridine and a significant stabilizing role in the binding of ethidium in its sterically constrained position, but the intercalation of daunomycin is determined by the side-chain binding. Sequence preferences are likely to be determined by a complex and subtle mixture of effects, with electrostatics being just one component. The electrostatic binding energy is also unlikely to be a major determinant of the twist angle, as its variation with angle is modest for most intercalation sites. Overall, the electrostatic potential maps give guidance on how positively charged chromophores can be chemically adapted by heteroatomic substitution to optimise their binding.

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.

Release of the cyano moiety in the crystal structure of N-cyanomethyl-N-(2-methoxyethyl)-daunomycin complexed with d(CGATCG)

Doxorubicin is among the most widely used anthracycline in cancer chemotherapy. In an attempt to avoid the cardiotoxicity and drug resistance of doxorubicin therapy, several analogues were synthesized. The cyanomorpholinyl derivative is the most cytotoxic. They differ greatly from their parent compound in their biological and pharmacological properties, inducing cross-links in drug DNA complexes. The present study concerns N-cyanomethyl-N-(2-methoxyethyl)-daunomycin (CMDa), a synthetic analogue of cyanomorpholino-daunomycin. Compared to doxorubicin, CMDa displays a cytotoxic activity on L1210 leukemia cells at higher concentration but is effective on doxorubicin resistant cells. The results of fluorescence quenching experiments as well as the melting temperature (DeltaTm = 7.5 degrees C) studies are consistent with a drug molecule which intercalates between the DNA base pairs and stabilizes the DNA double helix. The crystal structure of CMDa complexed to the hexanucleotide d(CGATCG) has been determined at 1.5 A resolution. The complex crystallizes in the space group P41212 and is similar to other anthracycline-hexanucleotide complexes. In the crystal state, the observed densities indicate the formation of N-hydroxymethyl-N-(2-methoxyethyl)-daunomycin (HMDa) with the release of the cyano moiety without DNA alkylation. The formation of this degradation compound is discussed in relation with other drug modifications when binding to DNA. Comparison with two other drug-DNA crystal structures suggests a correlation between a slight change in DNA conformation and the nature of the amino sugar substituents at the N3' position located in the minor groove.

Identification and hydropathic characterization of structural features affecting sequence specificity for doxorubicin intercalation into DNA double-stranded polynucleotides

The computer molecular modeling program HINT (Hydropathic INTeractions), an empirical hydropathic force field function that includes hydrogen bonding, coulombic and hydrophobic terms, was used to study sequence-selective doxorubicin binding/intercalation in the 64 unique CAxy, CGxy, TAxy, TGxy base pair quartet combinations. The CAAT quartet sequence is shown to have the highest binding score of the 64 combinations. Of the two regularly alternating polynucleotides, d(CGCGCG)2and d(TATATA)2, the HINT calculated binding scores reveal doxorubicin binds preferentially to d(TATATA)2. Although interactions of the chromophore with the DNA base pairs defining the intercalation site [I-1] [I+1] and the neighboring [I+2] base pair are predominant, the results obtained with HINT indicate that the base pair [I+3] contributes significantly to the sequence selectivity of doxorubicin by providing an additional hydrogen bonding opportunity for the N3' ammonium of the daunosamine sugar moiety in approximately 25% of the sequences. This observation, that interactions involving a base pair [I+3] distal to the intercalation site play a significant role in stabilizing/destabilizing the intercalation of doxorubicin into the various DNA sequences, has not been previously reported. In general terms, this work shows that molecular modeling and careful analysis of molecular interactions can have a significant role in designing and evaluating nucleotides and antineoplastic agents.

Comparison of DNA sequence selectivity of anthracycline antibiotics and their 3'-hydroxylated analogs

The sequence selectivity of three anthracyclines and their 3' hydroxylated analogs (in which an OH replaces NH3+ in the daunosamine at neutral pH) was examined in DNase I footprinting experiments on a 158-bp DNA fragment. We found that chemical modification of the daunosamine at C3' has more drastic consequences for sequence selectivity than chemical modification at C4 and C14 of the aglycone moiety. All anthracyclines and hydroxylated derivatives selectively recognize the triplet PyAPy. The importance of NH3+ in stabilizing the interaction was evidenced. First of all, comparable protection patterns require 5 times more hydroxyanthracycline than regular anthracycline. Furthermore, it is only after the replacement of NH3+ by OH that an additional protection site - CGC--appears. GGC is the site of best selectivity of the hydroxyanthracyclines. Anthracyclines can be considered both intercalators (aglycone moiety) and minor groove binders (sugar moiety). Since intercalating drugs show a slight preference for GC base pairs, we suggest hydroxylated anthracyclines to have a sequence specificity closer that of pure intercalators. Chemical modifications at C4 and C14 only modify the hydrogen bonding stabilization of the DNA-aglycone moiety complex: the more the anthracycline or its analog is lipophilic, the less it will interact with the sugar-phosphate chain.

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.

Infrared linear dichroism studies of DNA-drug complexes: quantitative determination of the drug-induced restriction of the B-A transition

The B-A transition of films or fibers of NaDNA occurs at a relative humidity of 75-85%. The fraction of DNA that changed the conformation from B to A form can be determined quantitatively by infrared linear dichroism. DNA-binding drugs can 'freeze' a fraction of DNA in the B form. This fraction of DNA is in the B form and cannot be converted to A-DNA even at a reduced relative humidity of 54%. The 'freezing' potentiality of various drugs can be described by the 'freezing' index, FI, expressed in base pairs per added drug. Drugs with a high value of FI (more than eight base pairs per drug) were observed among both intercalating and groove-binding drugs. High values of FI imply restriction of the conformational flexibility of DNA significantly going beyond the binding site of the drug. This long-range effect of drugs on the conformational flexibility of DNA may be connected with the molecular mechanism of drug action. The freezing index FI is a new quantitative parameter of drug-DNA interaction that should be considered as a valuable tool for drug design.

Transverse dipoles added to DNA chains by drug binding can induce inversion of the long-range chirality of DNA condensates

The addition of poly(ethylene glycol) (PEG) to a DNA solution induces phase separation of droplets of condensed DNA. These droplets possess liquid crystalline properties and their ordering is cholesteric. It was recently proved that daunomycin, by binding to DNA chains, inverts the long-range chirality of their tertiary packing into aggregates. The present paper suggests one possible mechanism by which this inversion can take place. Daunomycin bears a cationic group in its sugar residue. Its intercalation adds a helicoidal distribution of transverse dipoles to DNA chains. By this mechanism, in favourable cases, ionic or strongly polar groups in drugs which bind DNA can induce handedness inversion of the cholesteric ordering of its condensates. This inversion mechanism was tested experimentally using several, charged and uncharged, homologues of daunomycin. All those bearing the cationic ammonium group inverted the long-range chirality of the PEG-induced DNA mesomorphic state. The effects of the uncharged desamino homologues could not be evaluated because of their lower solubility and binding affinity for DNA.

Relative binding free energy calculations of DNA to daunomycin and its 13-dihydro analogue

We present the results of thermodynamic integration/molecular dynamics on the 13-dihydrodaunomycin.d(CGTACG) complex and 13-dihydrodaunomycin in vacuo based on the GROMOS force field. The objective is to calculate the binding free energy differences between two complexes--13-dihydrodaunomycin.d(CGTACG) and daunomycin.d(CGTACG). The results are in agreement with the experimental data, i.e. the theoretical binding free energy difference of 13-dihydrodaunomycin.d(CGTACG) and daunomycin.d(CGTACG) is 2.1 +/- 1.6 kcal/mol which accords well with the experimental value of 0.6 kcal/mol. The free energy contributions of different interactions and different groups have also been analysed.

Drug-DNA sequence-dependent interactions analysed by electric linear dichroism

The interactions between 20 drugs and a variety of synthetic DNA polymers and natural DNAs were studied by electric linear dichroism (ELD). All compounds tested, including several clinically used antitumour agents, are thought to exert their biological activities mainly by virtue of their abilities to bind to DNA. The selected drugs include intercalating agents with fused and unfused aromatic structures and several groove binders. To examine the role of base composition and base sequence in the binding of these drugs to DNA, ELD experiments were carried out with natural DNAs of widely differing base composition as well as with polynucleotides containing defined alternating and non-alternating repeating sequences, poly(dA).poly(dT), poly(dA-dT).poly(dA-dT),poly(dG).poly(dC) and poly(dG-dC).poly(dG-dC). Among intercalating agents, actinomycin D was found to be by far the most GC-selective. GC selectivity was also observed with an amsacrine-4-carboxamide derivative and to a lesser extent with methylene blue. In contrast, the binding of amsacrine and 9-aminoacridine was practically unaffected by varying the GC content of the DNAs. Ethidium bromide, proflavine, mitoxantrone, daunomycin and an ellipticine derivative were found to bind best to alternating purine-pyrimidine sequences regardless of their nature. ELD measurements provided evidence for non-specific intercalation of amiloride. A significant AT selectivity was observed with hycanthone and lucanthone. The triphenyl methane dye methyl green was found to exhibit positive and negative dichroism signals at AT and GC sites, respectively, showing that the mode of binding of a drug can change markedly with the DNA base composition. Among minor groove binders, the N-methylpyrrole carboxamide-containing antibiotics netropsin and distamycin bound to DNA with very pronounced AT specificity, as expected. More interestingly the dye Hoechst 33258, berenil and a thiazole-containing lexitropsin elicited negative reduced dichroism in the presence of GC-rich DNA which is totally inconsistent with a groove binding process. We postulate that these three drugs share with the trypanocide 4',6-diamidino-2-phenylindole (DAPI) the property of intercalating at GC-rich sites and binding to the minor groove of DNA at other sites. Replacement of guanines by inosines (i.e., removal of the protruding exocyclic C-2 amino group of guanine) restored minor groove binding of DAPI, Hoechst 33258 and berenil. Thus there are several cases where the mode of binding to DNA is directly dependent on the base composition of the polymer. Consequently the ELD technique appears uniquely valuable as a means of investigating the possibility of sequence-dependent recognition of DNA by drugs.

Anthracycline binding to DNA. High-resolution structure of d(TGTACA) complexed with 4'-epiadriamycin

Crystallographic methods have been applied to determine the high-resolution structure of the complex formed between the self-complementary oligonucleotide d(TGTACA) and the anthracycline antibiotic 4'-epiadriamycin. The complex crystallises in the tetragonal system, space group P4(1)2(1)2 with a = 2.802 nm and c = 5.293 nm, and an asymmetric unit consisting of a single DNA strand, one drug molecule and 34 solvent molecules. The refinement converged with an R factor of 0.17 for the 2381 reflections with F greater than or equal to 3 sigma F in the resolution range 0.70-0.14 nm. Two asymmetric units associate such that a distorted B-DNA-type hexanucleotide duplex is formed incorporating two drug molecules that are intercalated at the TpG steps. The amino sugar of 4'-epiadriamycin binds in the minor groove of the duplex and displays different interactions from those observed in previously determined structures. Interactions between the hydrophilic groups of the amino sugar and the oligonucleotide are all mediated by solvent molecules. Ultraviolet melting measurements and comparison with other anthracycline-DNA complexes suggest that these indirect interactions have a powerful stabilising effect on the complex.

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

The anticancer drug daunomycin has been co-crystallized with the hexanucleotide duplex sequences d(TGTACA) and d(TGATCA) and single crystal X-ray diffraction studies of these two complexes have been carried out. Structure solution of the d(TGTACA) and d(TGATCA) complexes to 1.6 and 1.7 Angstrom resolution, respectively, shows two daunomycin molecules bound to the DNA hexamer. Binding occurs via intercalation of the drug chromophore at the d(TpG) step, and hydrogen bonding interactions involving the drug, DNA and solvent molecules. The daunomycin sugar is located in the minor groove of the DNA hexamer and is stabilized by hydrogen bonds between the amino group of the sugar and functional groups on the floor of the groove. The amino sugar of the d(TGATCA) duplex interacts directly with the DNA sequence, while in the d(TGTACA) duplex, the interaction is via solvent molecules. Two other complexes d(CGTACG)-daunomycin and d(CGATCG)-daunomycin have previously been structurally characterized. Comparison of the four structures with daunomycin bound to the triplet sequences 5'TGT, 5'TGA, 5'CGT and 5'CGA reveals changes in the conformation of both the DNA hexamer and the daunomycin upon complexation, as well as the hydrogen bonding and van der Waals' interactions.

Structure of 11-deoxydaunomycin bound to DNA containing a phosphorothioate

The anthracyclines form an important family of cancer chemotherapeutic agents with a strong dependence of clinical properties on minor differences in chemical structure. We describe the X-ray crystallographic solution of the three-dimensional structure of the anthracycline 11-deoxydaunomycin plus d(CGTsACG). In this complex, two drug molecules bind to each hexamer duplex. Both the drug and the DNA are covalently modified in this complex in contrast with the three previously reported DNA-anthracycline complexes. In the 11-deoxydaunomycin complex the 11 hydroxyl group is absent and a phosphate oxygen at the TpA step has been replaced by a sulfur atom leading to a phosphorothioate with absolute stereochemistry R. Surprisingly, removal of a hydroxyl group from the 11 position does not alter the relative orientation of the intercalated chromophore. However, it appears that the phosphorothioate modification influenced the crystallization and caused the 11-deoxydaunomycin-d(CGTsACG) complex to crystallize into a different lattice (space group P2) with different lattice contacts and packing forces than the non-phosphorothioated DNA-anthracycline complexes (space group P4(1)2(1)2). In the minor groove of the DNA, the unexpected position of the amino-sugar of 11-deoxydaunomycin supports the hypothesis that in solution the position of the amino sugar is dynamic.

Biophysical chemistry of the daunomycin-DNA interaction

Daunomycin (daunorubicin) is a potent anticancer antibiotic that binds to DNA by the process of intercalation. Fundamental aspects of the physical chemistry of the daunomycin-DNA interaction are reviewed here, including the thermodynamics and kinetics of the binding reaction, and recent work that indicates that daunomycin binds preferentially to certain sites along the DNA lattice. The solution studies reviewed here combine with recent theoretical and crystallographic investigations to make the daunomycin-DNA interaction one of the best-characterized intercalation reactions. The molecular interactions that stabilize the daunomycin-DNA complex, and which contribute to its sequence preference, are discussed.

Modulation of the B-A transition of DNA by potential antitumor antibiotics. Influence of the base composition of DNA

The B-A transition of DNA in oriented films of DNA-drug complexes is more or less restricted as a consequence of drug binding as revealed by infrared linear dichroism. A fraction of DNA is irreversibly locked into the B form. This behavior is described by the number of DNA base pairs "frozen" in the B form by one drug molecule. This quantity is dependent on the DNA sequence the drug is attached to. In this paper, drug complexes of oriented films of NaDNA with a GC content of 42% from calf thymus and a GC-rich DNA from Micrococcus lysodeikticus were compared. The restriction of the B-A transition of DNA complexes with two intercalating antibiotics, aclacinomycin A and violamycin BI, is not severely influenced by the base composition of DNA. By contrast, the strong groove binding oligopeptide antibiotics netropsin and distamycin A are much less effective to restrict the B-A transition of GC-rich DNA than of AT-rich DNA. This finding is in agreement with previous results by other methods which support a model based upon a strong preference of AT clusters by these two non-intercalating drugs.

Theoretical design of novel, 4 base pair selective derivatives of mitoxantrone

Mitoxantrone (MTX) is a recently synthesized antitumor intercalative molecule, currently in use in chemotherapy. Previous theoretical computations showed that the base pair selectivity of MTX is limited to the sole two base-pair sequence making up the intercalation site. In order to further extend the recognition site, we undertook, by means of theoretical computations, the design of novel MTX derivatives, in which the terminal hydroxyl group of each side chain is esterified with oligopeptides. We compare in the present study the binding affinities of two derivatives, depsiGly-Lys(D) and depsiGly-Gly-Orn(L), for the palindromic sequences d(CCCGGG)2, d(GCCGGC)2, d(GGCGCC)2, and d(CGCGCG)2. Major groove binding of the oligopeptide arms was shown to be significantly more favourable than either minor groove binding, or binding to the sole phosphate groups. With the two arms adopting two antiparallel directions, two distinct arrangements were investigated in the major groove: (a) the two oligopeptides are brought closer together by means of two hydrogen bonds involving the backbone of their second residue in a beta-sheet like arrangement; (b) the two arms are remote from each other so as to reduce their mutual electrostatic repulsion. Whatever the disposition, the optimal binding configurations were invariably found to be those in which the cationic side chains of the terminal residues chelate N7/O6 of two successive guanines, whenever present on a given strand. A distinct energetical preference for arrangement (a) was obtained with the depsiGly-Gly-Orn(L) derivative. Replacement of the central Gly residue by a Cys one, as in the sequence depsiGly-Cys-Orn(L), was proposed subsequently, so as to further stabilize such a beta-sheet arrangement by means of a disulfide bridge between the two Cys residues. The two investigated compounds were shown to preferentially bind sequences d(CCCGGG)2 and d(GCCGGC)2, with a tetrameric core CCGG rather than sequences d(GGCGCC)2 and d(CGCGCG)2, with a tetrameric core GCGC.

Free energy calculation on base specificity of drug--DNA interactions: application to daunomycin and acridine intercalation into DNA

We present the results of free energy perturbation/molecular dynamics studies on B-DNA.daunomycin and B-DNA.9-aminoacridine complexes as well as on B-DNA itself in order to calculate the free energy differences between complexes having different base pair sequences. The results generally reproduce the trends observed experimentally, i.e., preferences of acridine and daunomycin to bind to a specific base sequence in the DNA. This is encouraging, given the simplicity of the molecular mechanical/dynamical model in which solvent is not explicitly included.

A theoretical investigation of the base sequence preferences of monointercalating polymethylene carboxamide derivatives 9-aminoacridine

Theoretical computations are performed of the comparative binding affinities of five polymethylene carboxamide derivatives of 9-aminoacridine to a series of double-stranded hexanucleotides. The purpose of this investigation is to ascertain whether minor groove recognition of a guanine base adjacent to the intercalation site can occur, and be preferentially stabilized, for a given length of the polymethylene side chain, encompassing from n = 2 up to n = 6 methylene groups. For that purpose, several representative sequences were investigated, in which intercalation of the 9-aminoacridine chromophore occurred at a central d(CpG) or d(TpA) step. Investigated were the self-complementary sequences d(CGCGCG)2, d(GCCGGC)2, d(TATATA)2 and d(ATTAAT)2, as well as the 'mixed' sequences d(ACTAAT) .d(ATTAGT) and d(TGTATA). d(TATACA). For n = 3 up to n = 6, such a recognition was enabled only when the guanine base was located downstream of the intercalation site, i.e. with steps d(CGG) and d(TAG). It occurred by means of a bidentate interaction involving, on the one hand, H(N2) and N3 of the base, and, on the other hand, the carbonyl oxygen and the cis amino hydrogen of the terminal formamide moiety of the ligand. Because of the flexibility of the side chain, however, alternative binding modes were also found to occur competitively, involving backbone-only interactions of the side chain. On the basis of the present computations, upon binding to the sequence d(GCCGGC)2, an optimal value of n = 5 could be derived, with the corresponding acridine derivative eliciting both a significant prevalence of the bidentate over backbone only binding mode, and the most favourable energy balance within the investigated series. This privileged value of n = 5 is fully consistent with the experimental results of Markovits et al. and Gaugain et al. The very flexibility of the side chain, however, hampered any preferential recognition of a triplet sequence with a downstream guanine, such as d(CGG) or d(TAG), to be elicited over sequences such as d(TAA), d(TAT) or d(TAC).

Association of anthracyclines and synthetic hexanucleotides. Structural factors influencing sequence specificity

The equilibrium and kinetic aspects of the interaction between four anthracyclines and two synthetic self-complementary hexanucleotides was investigated by fluorescence detection. Two of the studied anthracyclines are widely used antitumor drugs: doxorubicin (1, formerly adriamycin) and daunorubicin (2, formerly daunomycin). The other two, 9-deoxydoxorubicin (3) and 3'-deamino-3'-hydroxy-4'-epidoxorubicin (4), are doxorubicin analogues with modifications of the chemical groups that have been proposed as responsible for sequence specificity (Chen, K.-X., Gresh, N. and Pullman, B. (1985). J. Biomol. Struct. Dyn. 3, 445-466). One of the oligonucleotides, d(CGTACG), is identical to that used in the high resolution x-ray structure determination of the daunorubicin intercalative complex (Wang, A. H.-J., Ughetto, G., Quigley, G. J. & Rich, A. (1987). Biochemistry 26, 1152-1163). Binding to this hexanucleotide is compared with intercalation into the d(CGCGCG) duplex, revealing sequence preferences of the four anthracyclines. Taking into account the anthracycline aggregation and the dissociation of the hexanucleotide double standard form, results can be interpreted with a model that assumes complete fluorescence quenching at intercalative sites containing the CG base pair, and a large residual fluorescence after intercalation within the TpA fragment. All four anthracyclines show preferential intercalation at sites near the ends of both hexanucleotide duplexes, partly as a result of positive cooperativity in the formation of di-intercalated species at these sites. Within the limits of experimental error, complete site specificity for the CpG fragment is found in the intercalation of 1 and 2 into d(CGTACG) duplex, whereas analogues 3 and 4 give increasing evidence of intercalation at other sites including the fluorescence-preserving TpA fragment. Site specificity is less pronounced in the association with d(CGCGCG), when cooperativity is taken into account. Kinetic data corroborate the results of equilibrium studies and are interpreted with a mechanism that includes formation of an intermediate bound species followed by drug redistribution to preferential sites. Finally, from a comparison of pertinent site binding constants, approximate free energy contributions to sequence specific DNA interaction, due to C9-OH on the aglycone and -NH3+ on daunosamine, are estimated not to exceed 2 kcal/mol.

Sequence selective binding of ditrisarubicin B to DNA: comparison with daunomycin

DNase I footprinting has been used to examine the sequence selective binding of ditrisarubicin B, a novel anthracycline antibiotic, to DNA. At 37 degrees C no footprinting pattern is observed, the drug protects all sites from enzymic cleavage with equal efficiency. At 4 degrees C a footprinting pattern is induced with low drug concentrations which is different from that produced by daunomycin. The best binding sites contain the dinucleotide step GpT (ApC) and are located in regions of alternating purines and pyrimidines.

In vitro transcription analysis of the role of flanking sequence on the DNA sequence specificity of adriamycin

An in vitro transcription assay has been used to define twelve high occupancy transcriptional blockage sites in 260 bp of heterogenous DNA on three promoter-containing DNA fragments. Transcription proceeded to immediately upstream of CpA sequences in nine of these sites, and this defines the most preferred intercalation site as CpA. In almost all cases, this sequence was flanked by T on the 5' end. The consensus sequence for the highest affinity Adriamycin site is therefore 5'-TCA, with some evidence for preference of AT base pairs flanking both ends of this trinucleotide [i.e., (t)TCA(a.t)(a.t)]. In contrast, lower occupancy CpA sites were flanked on the 5' end by a GC base pair, with a preference for up to three GC base pairs on the 3' end. The low affinity consensus sequence is therefore (G.C)CA(g.c)(g.c)(g.c).

A theoretical investigation of the intercalative binding of 7-H pyrido[4.3C]carbazole chromophore into a d(CpG)2 minihelix

Theoretical computations are performed of the intercalative binding to a model d(CpG)2 minihelix of 7-H pyrido[4.3C]carbazole, the precursor of the antitumor bisintercalating drug ditercalinium. The conformations of the intercalation site are generated by the AGNAS procedure (algorithm to generate nucleic acid structures) of Miller and co-workers. The ligand-nucleotide interactions and the nucleotide conformational energies are computed with the SIBFA procedures (sum of interactions between fragments ab initio computed), which use formulas of empirical origin that reproduce ab initio SCF (self-consistent field) computations. Among the candidate intercalation sites most favored energetically, one has a pattern of conformational angles related to the one determined crystallographically by Sobell et al. in a series of x-ray structural studies of small intercalator-dinucleotide monophosphate complexes. Optimal values of the unwinding angle, found in the range of -12 degrees to -14 degrees, are consistent with available experimental data on DNA.

A 19F-NMR study of 2-fluoro-4-demethoxydaunomycin intercalation complexes with the decanucleotides d(G-C)5 and d(A-T)5

Binding configurations and equilibria of intercalation complexes formed by the novel anthracycline drug, 2-fluoro-4-demethoxydaunomycin (2FD), with the decanucleotides d(G-C)5 and d(A-T)5 have been studied by 19F-NMR spectroscopy. The 19F chemical shift of 2FD bound to d(A-T)5 was approximately 1.5 ppm downfield of that observed for 2FD bound to d(G-C)5. By mixing equimolar amounts of aqueous d(G-C)5, d(A-T)5 and 2FD, the distribution of drug between the nucleotides was followed by observing relative peak intensities and showed no G-C or A-T binding preference at room temperature. It was shown that each decanucleotide duplex bound three 2FD molecules, giving a neighbour exclusion parameter, n, of n = 3 for this drug. The stoichiometric complexes, which we denote by [d(A-T)5][2FD]3 and [d(G-C)5][2FD]3, were also purified and isolated in this study.

31P NMR study of daunorubicin-d(CGTACG) complex in solution. Evidence of the intercalation sites

The interaction of daunorubicin with the self-complementary DNA fragment d(CGTACG) was studied by 31P NMR spectroscopy. The individual phosphates have been assigned for the nucleotide and the complex and signals from bound and free species in slow exchange at 19 degrees C were detected. In solution, the hexanucleotide binds two molecules of daunorubicin, which intercalate in the d(CG) sequence at both ends of the helix. Evidence for local deformations of the backbone at the sites of C5pG6, C1pG2 and G2pT3 phosphates is given. The binding constants for the stepwise equilibrium and the rate of dissociation of the intercalated duplex were also determined.