Calorimetric and spectroscopic investigation of drug--DNA interactions: II. Dipyrandium binding to poly d(AT)

Solid-phase synthesis of oligodeoxynucleotides using nucleobase N-unprotected oxazaphospholidine derivatives bearing a long alkyl chain

In this study, we developed a new approach for the solid-phase synthesis of oligodeoxynucleotides (ODNs) using nucleobase-unprotected oxazaphospholidine derivatives. We tackled the problem of the difficult purification of N-unprotected monomers due to their high affinity to silica gel by introducing a tetrahydrogeranyl group into the oxazaphospholidine monomers, thereby enhancing the lipophilicity and facilitating the isolation. In addition, the cyclic structure of oxazaphospholidine enabled a hydroxy-group-selective condensation with sufficient efficiency. Unmodified and boranophosphate/phosphate chimeric ODNs were successfully synthesized using this strategy. This synthetic method can be expected to afford ODNs containing base-labile functional groups.

Forces Driving a Magic Bullet to Its Target: Revisiting the Role of Thermodynamics in Drug Design, Development, and Optimization

Drug discovery strategies have advanced significantly towards prioritizing target selectivity to achieve the longstanding goal of identifying "magic bullets" amongst thousands of chemical molecules screened for therapeutic efficacy. A myriad of emerging and existing health threats, including the SARS-CoV-2 pandemic, alarming increase in bacterial resistance, and potentially fatal chronic ailments, such as cancer, cardiovascular disease, and neurodegeneration, have incentivized the discovery of novel therapeutics in treatment regimens. The design, development, and optimization of lead compounds represent an arduous and time-consuming process that necessitates the assessment of specific criteria and metrics derived via multidisciplinary approaches incorporating functional, structural, and energetic properties. The present review focuses on specific methodologies and technologies aimed at advancing drug development with particular emphasis on the role of thermodynamics in elucidating the underlying forces governing ligand-target interaction selectivity and specificity. In the pursuit of novel therapeutics, isothermal titration calorimetry (ITC) has been utilized extensively over the past two decades to bolster drug discovery efforts, yielding information-rich thermodynamic binding signatures. A wealth of studies recognizes the need for mining thermodynamic databases to critically examine and evaluate prospective drug candidates on the basis of available metrics. The ultimate power and utility of thermodynamics within drug discovery strategies reside in the characterization and comparison of intrinsic binding signatures that facilitate the elucidation of structural-energetic correlations which assist in lead compound identification and optimization to improve overall therapeutic efficacy.

Thermodynamic Factors That Drive Sequence-Specific DNA Binding of Designed, Synthetic Minor Groove Binding Agents

Ken Breslauer began studies on the thermodynamics of small cationic molecules binding in the DNA minor groove over 30 years ago, and the studies reported here are an extension of those ground-breaking reports. The goals of this report are to develop a detailed understanding of the binding thermodynamics of pyridine-based sequence-specific minor groove binders that have different terminal cationic groups. We apply biosensor-surface plasmon resonance and ITC methods to extend the understanding of minor groove binders in two directions: (i) by using designed, heterocyclic dicationic minor groove binders that can incorporate a G•C base pair (bp), with flanking AT base pairs, into their DNA recognition site, and bind to DNA sequences specifically; and (ii) by using a range of flanking AT sequences to better define molecular recognition of the minor groove. A G•C bp in the DNA recognition site causes a generally more negative binding enthalpy than with most previously used pure AT binding sites. The binding is enthalpy-driven at 25 °C and above. The flanking AT sequences also have a large effect on the binding energetics with the -AAAGTTT- site having the strongest affinity. As a result of these studies, we now have a much better understanding of the effects of the DNA sequence and compound structure on the molecular recognition and thermodynamics of minor groove complexes.

Conformational changes in DNA upon ligand binding monitored by circular dichroism

Circular dichroism (CD) spectroscopy is an optical technique that measures the difference in the absorption of left and right circularly polarized light. This technique has been widely employed in the studies of nucleic acids structures and the use of it to monitor conformational polymorphism of DNA has grown tremendously in the past few decades. DNA may undergo conformational changes to B-form, A-form, Z-form, quadruplexes, triplexes and other structures as a result of the binding process to different compounds. Here we review the recent CD spectroscopic studies of the induction of DNA conformational changes by different ligands, which includes metal derivative complex of aureolic family drugs, actinomycin D, neomycin, cisplatin, and polyamine. It is clear that CD spectroscopy is extremely sensitive and relatively inexpensive, as compared with other techniques. These studies show that CD spectroscopy is a powerful technique to monitor DNA conformational changes resulting from drug binding and also shows its potential to be a drug-screening platform in the future.

Energetic diversity of DNA minor-groove recognition by small molecules displayed through some model ligand-DNA systems

Energetics of interactions occurring in the model ligand-DNA systems constituted from distamycin A (DST), netropsin (NET) and the oligomeric duplexes d(GCAAGTTGCGATATACG)d(CGTATATCGCAACTTGC)=D#1 and d(GCAAGTTGCGAAAAACG)d(CGTTTTTCGCAACTTGC)=D#2 was studied by spectropolarimetry, UV-absorption spectroscopy and isothermal titration calorimetry. Model analysis of the measured signals was applied to describe individual and competitive binding in terms of populations of various species in the solution. Our results reveal several unprecedented ligand-DNA binding features. DST binds to the neighboring 5'-AAGTT-3' and 5'-ATATA-3' sites of D#1 statistically in a 2:1 binding mode. By contrast, its association to D#2 appears to be a 2:1 binding event only at the DST/D#2 molar ratios between 0 and 2 while its further binding to D#2 may be considered as a step-by-step binding to the unoccupied 5'-AAAAA-3' sites resulting first in DST3D#2 and finally in DST4D#2 complex formation. Competition between DST and NET binding shows that for the most part DST displaces NET from its complexes with D#1 and D#2. In contrast to the obligatory 1:1 binding of DST to the ligand-free 5'-AAAAA-3' sites observed at DST/5'-AAAAA-3' <1 the displacement of NET bound to the 5'-AAAAA-3' sites by added DST occurs even at the smallest additions of DST in a 2:1 manner. NET can also displace DST molecules but only those bound monomerically to the 5'-AAAAA-3' sites of DST3D#2. Actually, only half of these molecules can be displaced due to the simultaneous rebinding of the displaced DST to the unreacted 5'-AAAAA-3' sites in DST3D#2. Binding of DST and NET to D#1 and D#2 is an enthalpy driven process accompanied by large unfavorable (DST), small (NET) or large favorable (NET binding to 5'-AAAAA-3') entropy contributions and negative deltaCP degrees that are reasonably close to deltaCP degrees predicted from the calculated changes in solvent-accessible surface areas that accompany complex formation. Although various modes of DST and NET binding within D#1 and D#2 are characterized by significant energetic differences they seem to be governed by the same driving forces; the hydrophobic transfer of ligand from the solution into the duplex binding site and the accompanying specific non-covalent ligand-DNA and ligand-ligand interactions occurring within the DNA minor groove.

Effects of polyamines on the thermal stability and formation kinetics of DNA duplexes with abnormal structure

The effects of ions (i.e. Na+, Mg2+ and polyamines including spermidine and spermine) on the stability of various DNA oligonucleotides in solution were studied. These synthetic DNA molecules contained sequences that mimic various cellular DNA structures, such as duplexes, bulged loops, hairpins and/or mismatched base pairs. Melting temperature curves obtained from the ultraviolet spectroscopic experiments indicated that the effectiveness of the stabilization of cations on the duplex formation follows the order of spermine > spermidine > Mg2+ > Na+ > Tris-HCl buffer alone at pH 7.3. Circular dichroism spectra showed that salts and polyamines did not change the secondary structures of those DNA molecules under study. Surface plasmon resonance (SPR) observations suggested that the rates of duplex formation are independent of the kind of cations used or the structure of the duplexes. However, the rate constants of DNA duplex dissociation decrease in the same order when those cations are involved. The enhancement of the duplex stability by polyamines, especially spermine, can compensate for the instability caused by abnormal structures (e.g. bulged loops, hairpins or mismatches). The effects can be so great as to make the abnormal DNAs as stable as the perfect duplex, both kinetically and thermodynamically. Our results may suggest that the interconversion of various DNA structures can be accomplished readily in the presence of polyamine. This may be relevant in understanding the role of DNA polymorphism in cells.

Modulation of nucleic acid structure by ligand binding: induction of a DNA.RNA.DNA hybrid triplex by DAPI intercalation

The aromatic diamidine, DAPI (4',6-diamidino-2-phenylindole), is used as an important biological and cytological tool since it forms highly fluorescent complexes with nucleic acid duplexes via minor groove-directed/intercalative modes of interaction. In this study, we find that DAPI binding can induce the formation of an RNA-DNA hybrid triplex that would not otherwise form. More specifically, through application of a broad range of spectroscopic, viscometric, and molecular modeling techniques, we demonstrate that DAPI intercalation induces the formation of the poly(dT).poly(rA).poly(dT) hybrid triple helix, a structure which does not form in the absence of the ligand. Using UV mixing studies, we demonstrate that, in the presence of DAPI, the poly(rA).poly(dT) duplex and the poly(dT) single strand form a 1:1 complex (a triplex) that does not form in the absence of DAPI. Through temperature-dependent absorbance measurements, we show that the poly(dT).poly(rA).poly(dT) triplex melts via two distinct transitions: initial conversion of the triplex to the duplex state, with the DAPI remaining bound, followed by denaturation of the duplex-DAPI complex to its component single strands and free DAPI. Using optical melting profiles, we show that DAPI binding enhances the thermal stability of the poly(dT).poly(rA).poly(dT) triplex, an observation consistent with the preferential binding of the ligand to the triplex versus the duplex and single-stranded states. Our differential scanning calorimetric measurements reveal melting of the DAPI-saturated poly(dT).poly(rA).poly(dT) triplex to be associated with a lower enthalpy but greater cooperativity than melting of the corresponding DAPI-saturated poly(rA).poly(dT) duplex. Our flow linear dichroism and viscometric data are consistent with an intercalative mode of binding when DAPI interacts with both the poly(dT).poly(rA).poly(dT) triplex and the poly(rA).poly(dT) duplex. Finally, computer modeling studies suggest that a combination of both stacking and electrostatic interactions between the intercalated ligand and the host nucleic acid play important roles in the DAPI-induced stabilization of the poly(dT).poly(rA).poly(dT) triplex. In the aggregate, our results demonstrate that ligand binding can be used to induce the formation of triplex structures that do not form in the absence of the ligand. This triplex-inducing capacity has potentially important implications in the design of novel antisense, antigene, antiviral, and diagnostic strategies.

Berenil [1,3-bis(4'-amidinophenyl)triazene] binding to DNA duplexes and to a RNA duplex: evidence for both intercalative and minor groove binding properties

Berenil is an antitrypanosomal agent that binds to nucleic acid duplexes. The generally accepted mode of berenil binding is via complexation into the minor groove of AT-rich domains of DNA double helices. We find that berenil can bind to RNA as well as DNA duplexes, while exhibiting properties characteristic of both intercalation as well as minor groove binding. More specifically, we use spectroscopic, calorimetric, and hydrodynamic techniques to characterize berenil binding to four DNA duplexes and to one RNA duplex. Our results reveal the following features: (i) Berenil binding to the poly[d(A-T)]2, poly(dA).poly(dT), poly[d(I-C)]2, poly[d(G-C)]2, and poly(rA).poly(rU) duplexes exhibits intercalative as well as minor groove binding characteristics. (ii) The apparent "site sizes" associated with berenil binding to these five duplexes range from 1 to 13 base pairs per bound berenil and depend, in part, on the host duplex. One of the site sizes common to all five duplexes is consistent with berenil binding to the minor groove. (iii) The apparent berenil binding affinity follows the hierarchy: poly(dA).poly(dT) > poly-[d(A-T)]2 approximately poly[d(I-C)]2 >> poly(rA).poly(rU) > poly[d(G-C)]2. (iv) Viscometric data reveal properties characteristic of a significant contribution from an intercalative mode of binding when berenil interacts with the poly[d(A-T)]2, poly[d(I-C)]2, poly[d(G-C)]2, and poly(rA).poly(rU) duplexes, while revealing an apparent nonintercalative mode when the drug binds to the poly(dA).poly(dT) duplex. (v) Berenil binding unwinds negative supercoils in the pBR322 plasmid, an observation consistent with an intercalative mode of binding to duplex DNA. (vi) Salt-dependent melting data suggest that both positively charged amidino groups of berenil participate in the complexation of the drug to the poly[d(I-C)]2, poly[d(A-T)]2, poly(dA).poly(dT), and poly(rA).poly(rU) duplexes, while also suggesting that the binding event is site-specific. In the aggregate, our results suggest that, in contrast to the conventional wisdom, berenil can exhibit intercalative as well as minor groove binding properties when it binds to both DNA and RNA duplexes, while also exhibiting a preference for DNA duplexes with unobstructed minor grooves. We comment on the potential correlation between drugs, such as berenil, that exhibit "mixed" binding motifs and those that express anticancer activity via inhibition of topoisomerase I activity.

A thermodynamic and spectroscopic study on the binding of berenil to poly d(AT) and to poly (dA) x poly (dT)

The complete thermodynamic profile for the non-intercalative binding of berenil to the alternating copolymer poly d(AT) and to the homopolymer poly (dA) x poly (dT) was investigated. Differential Scanning Calorimetry (DSC) and UV absorbance spectroscopy have been used to characterize and to compare the binding of berenil to the different synthetic polymers. Both double stranded DNA's show two types of binding; one stronger binding mode at low berenil concentrations and a weaker, in the case of poly d(AT)-berenil complexes slightly cooperative binding mode at higher drug to base pair ratios. For the interaction of berenil with poly d(AT) the thermodynamic data delta G(bind)0 = -33 kJ/mol drug, delta H(bind)0 = -29 kJ/mol of drug and delta S(bind)0 = +13 J/Kmol of drug were calculated. For the minor groove binding of berenil to poly (dA) x poly (dT) the following values were obtained: delta G(bind)0 = -34 kJ/mol of drug, delta H(bind)0 = -25 kJ/mol of drug and delta S(bind)0 = +30 J/Kmol of drug. Temperature-dependent UV absorbance spectroscopy revealed for both duplexes a biphasic "melting" behavior. However, the saturated nucleic acids (drug to base pair ratio 0.33) "melted" monophasically and with a decreased length of the cooperative unit. The obtained apparent equilibrium constants K(app) for the complexation with the discharged drug molecule showed to be a sensitive function of the ionic environment. But in contradiction to the expected release of two counterions into the solvent only a value of 1.0 was observed for the alternating copolymer poly d(AT). The complexation of berenil with poly (dA) x poly (dT) is followed by a release of 1.4 ions indicating stronger electrostatic interactions. For both polynucleotides the predicted release of two ions is not achieved. This is due to the presence of a binding mode, which involves less electrostatic interactions. From the complete data set it is proposed that the mode of binding is closely related to that found for the analogue minor groove binders DAPI and netropsin.

Differential scanning calorimetric study of the effect of intercalators and other kinds of DNA-binding drugs on the stepwise melting of plasmid DNA

The effect of intercalating drugs (the anthracycline group of antibiotics, ethidium bromide, actinomycin D) on stepwise melting of DNA was studied by differential scanning calorimetry (DSC). The DSC DNA melting profile of plasmid pJL3-TB5 DNA (5277 base-pairs in length) consists of seven peaks, and all the intercalators caused shifting of these peaks, particularly those formed at the high temperature ranges, to the higher temperature ranges in a characteristic manner depending upon the binding strength of the drug. The analysis of the anthracycline group of antibiotics, such as aclacinomycin A, daunomycin, adriamycin and pyrarubicin, indicates that the difference in binding is due to the sugar moiety at position O-7 of the chromophore in these antibiotics. Analysis on the basis of the helix-coil transition theory suggests that the anthracycline group of antibiotics interact preferentially with the 5'-CG-3' sequences. The effect of various DNA-binding drugs other than intercalators on stepwise melting of DNA was then studied by DSC. The representative drugs examined were distamycin A, peplomycin, cis-dichlorodiamine-platinum(II) (cis-DDP or cis-Platin) and mitomycin C, which differ in their mode of interaction with DNA; namely, minor groove binding, strand cleavage and intrastrand or interstrand cross-linking. Distamycin A caused shifting of the DSC peaks at the low temperature ranges to a higher temperature range, whereas peplomycin and cis-DDP caused shifting of all the DSC peaks to form a broad peak at a lower temperature range, suggesting that the DSC DNA melting profiles are affected in a characteristic manner depending upon the interaction mode of the drug.

Enthalpy-entropy compensations in drug-DNA binding studies

We present a comparative study of calorimetrically derived thermodynamic profiles for the binding of a series of drugs with selected DNA host duplexes. We use these data to demonstrate that comparisons between complete thermodynamic profiles (delta G zero, delta H zero, delta S zero, delta Cp) are required before drug binding can be used as a probe of DNA conformation, since enthalpy-entropy compensations can cause two drug-DNA binding events to exhibit similar binding free energies (delta G zero) despite being driven by entirely different thermodynamic forces (delta H zero, delta S zero). In this work, we employ a combination of spectroscopic and calorimetric techniques to characterize thermodynamically the DNA binding of netropsin and distamycin (two minor groove-directed ligands), ethidium (an intercalator), and daunomycin (a combined intercalator/groove binder). Our free energy data (delta G zero) show that each drug exhibits similar binding affinities at 25 degrees C for the alternating copolymer duplex poly[d(A-T)].poly[d(A-T)] and for the homopolymer duplex poly(dA).poly(dT). However, our calorimetric measurements reveal that the nature of the thermodynamic forces (delta H zero, delta S zero) that drive drug binding to these two host duplexes at 25 degrees C are entirely different, despite similar binding free energies (delta G zero) and similar salt dependencies (lnK/ln[Na+]). Specifically, the 25 degrees C binding of all four drugs to the alternating copolymer poly[d(A-T)].poly[d(A-T)] is overwhelmingly enthalpy driven, whereas the corresponding binding of each drug to the homopolymer duplex poly(dA).poly(dT) is overwhelmingly entropy driven. Thus, the similar binding free energies (delta G zero) we measure for complexation of each drug with poly[d(A-T)].poly[d(A-T)] and poly(dA).poly(dT) result from compensating changes in the enthalpy and entropy terms. Comparison with the thermodynamic profiles for the complexation of these drug molecules to other DNA host duplexes at 25 degrees C reveals that the binding of each is strongly enthalpy driven, except when the poly(dA).poly(dT) homopolymer serves as the host duplex. This comparison allows us to conclude that poly[d(A-T)].poly[d(A-T)] behaves thermodynamically as the more "normal" host duplex toward drug binding, whereas the entropy-driven binding to the poly(dA).poly(dT) duplex represents "aberrant" behavior. Furthermore, since each of the four drugs exhibits different modes of DNA binding, we conclude that the observed entropy-driven behavior for binding to poly(dA).poly(dT) reflects an intrinsic property of the homopolymer duplex that is perturbed in a common manner upon ligation rather than a common property of all four binding ligands. To rationalize the large positive entropy changes that drive drug complexation with poly(dA).poly(dT) duplex, we propose a model that emphasizes binding-induced perturbations of the more highly hydrated, altered B conformation of the homopolymer. Our results suggest that an aberrant thermodynamic binding profile may reflect an unusual DNA conformation in the host duplex. However, before such a conclusion can be reached, complete thermodynamic binding profiles must be examined, since enthalpy-entropy compensations can cause two binding events to exhibit similar binding constants even when they are driven by very different thermodynamic forces.