Drug-DNA recognition: energetics and implications for design

In this article we review thermodynamic studies designed to examine the interaction of low molecular weight ligands or drugs with DNA. Over the past 10 years there has been an increase in the number of rigorous biophysical studies of DNA-drug interactions and considerable insight has been gained into the energetics of these binding reactions. The advent of high-sensitivity calorimetric techniques has meant that the energetics of DNA-drug association reactions can be probed directly and enthalpic and entropic contributions to the binding free energy established. There are two principal consequences arising from this type of work, firstly three-dimensional structures of DNA-drug complexes from X-ray and NMR studies can be put into a thermodynamic context and the energetics responsible for stabilizing the observed structures can be more fully understood. Secondly, any rational approach to structure-based drug design requires a fundamental base of knowledge where structural detail and thermodynamic data on complex formation are intimately linked. Therefore these types of studies allow a set of general guidelines to be established, which can then be used to develop drug design algorithms. In this review we describe recent breakthroughs in duplex DNA-directed drug design and also discuss how similar principles are now being used to target higher-order DNA molecules, for example, triplex (three-stranded) and tetraplex (four-stranded) structures.     

Comparative Thermodynamics of Benzodiazepine Receptor Ligand Interactions in Rat Neuronal Membranes

The effects of temperature on the interaction of various ligands with the benzodiazepine receptor were studied in rat brain membrane preparations. The affinities of all ligands studied were reduced on raising the temperature from 4 to 37 degrees C. The variation of affinity constant with temperature deviated from the classical relationship for both the anticonvulsant ligand [3H]flunitrazepam and the proconvulsant ligand [3H]ethyl beta-carboline-3-carboxylate. This implies a variation of observed enthalpy change of binding with temperature. Possible reasons for this are discussed. Gamma-Aminobutyric acid and sodium chloride both enhance the binding of [3H]flunitrazepam--the former by an increase in the entropic component of the binding energy, and the latter by an increase in the enthalpic component. In a series of ligands of different biological activities, no simple correlation was observed between biological activity and temperature dependence of binding.     

Molecular Mechanisms,Thermodynamics, and Dissociation Kinetics of Knob-Hole Interactions in Fibrin

Polymerization of fibrin, the primary structural protein of blood clots and thrombi, occurs through binding of knobs ‘A’ and ‘B’ in the central nodule of fibrin monomer to complementary holes ‘a’ and ‘b’ in the γ- and β-nodules, respectively, of another monomer. We characterized the A:a and B:b knob-hole interactions under varying solution conditions using molecular dynamics simulations of the structural models of fibrin(ogen) fragment D complexed with synthetic peptides GPRP (knob ‘A’ mimetic) and GHRP (knob ‘B’ mimetic). The strength of A:a and B:b knob-hole complexes was roughly equal, decreasing with pulling force; however, the dissociation kinetics were sensitive to variations in acidity (pH 5–7) and temperature (T = 25–37 °C). There were similar structural changes in holes ‘a’ and ‘b’ during forced dissociation of the knob-hole complexes: elongation of loop I, stretching of the interior region, and translocation of the moveable flap. The disruption of the knob-hole interactions was not an “all-or-none” transition as it occurred through distinct two-step or single step pathways with or without intermediate states. The knob-hole bonds were stronger, tighter, and more brittle at pH 7 than at pH 5. The B:b knob-hole bonds were weaker, looser, and more compliant than the A:a knob-hole bonds at pH 7 but stronger, tighter, and less compliant at pH 5. Surprisingly, the knob-hole bonds were stronger, not weaker, at elevated temperature (T = 37 °C) compared with T = 25 °C due to the helix-to-coil transition in loop I that helps stabilize the bonds. These results provide detailed qualitative and quantitative characteristics underlying the most significant non-covalent interactions involved in fibrin polymerization.     

Ca2+ and Mg2+ Influence the Thermodynamics of Peptide-Membrane Interactions

Accurate quantitative estimates of protein-membrane interactions are critical to studies of membrane proteins. Here, we demonstrate that thermodynamic analyses based on current hydropathy scales do not account for the significant and experimentally determined effects that Ca²⁺ or Mg²⁺ have on protein-membrane interactions. We examined distinct modes of interaction (interfacial partitioning and folding and transmembrane insertion) by studying three highly divergent peptides: Bid-BH3 (derived from apoptotic regulator Bid), peripherin-2-derived prph2-CTER, and the cancer-targeting pH-Low-Insertion-Peptide (pHLIP). Fluorescence experiments demonstrate that adding 1–2 mM of divalent cations led to a substantially more favorable bilayer partitioning and insertion, with free energy differences of 5–15 kcal/mol.     

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.     

Kinetics and Thermodynamics of 1-anilino-8-naphthalene Sulfonate Interactions with Urinary Trypsin Inhibitor

The interactions between Urinary Trypsin Inhibitor (UTI) and 1-anilino-8-naphthalene sulfonate (ANS) were investigated by fluorescence spectra, isothermal titration calorimetry and molecular modeling. The results revealed the presence of four specific binding sites for ANS on UTI, with interactions driven mainly by electrostatic forces. The four specific binding sites indicated the involvement of four hydrophobic patches on UTI. Experimental data also confirmed the presence of a further five nonspecific binding sites that interacted mainly by the formation of salt bridges between the sulfonates of ANS and positive residues on the surface of UTI.     

BglⅠ限制性核酸内切酶与环状pBR322-DNA专一性和非专一性相互作用的动力学及热力学

 本文选择限制性核酸内切酶BglⅠ和pBR322-DNA为试验系统,用酶促反应的动力学和热力学方法来研究内切酶对环状DNA分子的专一性和非专一性结合及切割过程,求得了各限制位点的切割速度k及活化能E,各限制位点催化速度常数k_c,酶同限制位点专一性结合的平衡常数k_S非专一性结合的平衡常数k_N及其热力学参数△H,△S。研究表明:不论底物的构型如何(线状还是环状),内切酶都以相似的动力学和热力过程对其进行结合与切割。     

Drug-DNA interactions: spectroscopic and footprinting studies of site and sequence specificity of elliptinium.

The binding of the antitumoral ellipticine derivative 2-methyl-9-hydroxyellipticinium acetate (elliptinium; NMHE) to DNA was analyzed by the combined use of DNase I footprinting and spectroscopic methods. Using two fragments of pBR322 DNA, five discrete NMHE binding sites of 5-7 protected base pairs (bp) were detected by footprinting at 4 degrees C on the analyzed regions. These corresponded to alternating pyrimidines and purines. The inactive derivative 2-methyl ellipticinium acetate L(NME) lacking a hydroxy group failed to demonstrate DNA protection even at low temperature. Ultraviolet-absorption and 1H-nmr analysis was performed using two autocomplementary octanucleotides d(TGACGTCA) (I) and d(ACTGCAGT) (II). The uv-absorption titrations resulted in an intercalative binding mode for NMHE in the oligomers. Analysis of the derived biphasic Scatchard plots yielded two binding sites corresponding to approximately 6-bp and 2-bp sizes and characterized by apparent association constants K1 approximately 10(8) M-1 and K2 approximately 10(6) M-1, respectively. The 1H-nmr analysis of exchangeable (imino) protons and nonexchangeable protons performed in the one- and two-dimensional modes confirmed the intercalation of NMHE, and further revealed the existence of multiple sites on DNA. Assuming that imino resonance line width concerned the sole kinetic effects, 10-ms order lifetimes were estimated for the drug-oligonucleotide complexes at 7 degrees C, pH 7, and 0.1 ionic strength. Finally, examination of every drug-DNA spectra in the light of the footprinting results indicated that there was a preference for binding of NMHE to the CpG (octamer I) and TpG (octamers I and II) steps.     

Thermodynamics of mechanosensitivity

Mechanosensitivity of ion channels is conventionally interpreted as being driven by a change of their in-plane cross-sectional area A(msc). This, however, does not include any factors relating to membrane stiffness, thickness, spontaneous curvature or changes in channel shape, length or stiffness. Because the open probability of a channel is sensitive to all these factors, we constructed a general thermodynamic formalism. These equations provide the basis for the analysis of the behaviour of mechanosensitive channels in lipids of different geometric and chemical properties such as the hydrophobic mismatch at the boundary between the protein and lipid or the effects of changes in the bilayer intrinsic curvature caused by the adsorption of amphipaths. This model predicts that the midpoint gamma(1/2) and the slope(1/2) of the gating curve are generally not independent. Using this relationship, we have predicted the line tension at the channel/lipid border of MscL as approximately 10 pN, and found it to be much less than the line tension of aqueous pores in pure lipid membranes. The MscL channel appears quite well matched to its lipid environment. Using gramicidin as a model system, we have explained its observed conversion from stretch-activated to stretch-inactivated gating as a function of bilayer thickness and composition.     

Mutation Induced Conformational Changes in Genomic DNA from Cancerous K562 Cells Influence Drug-DNA Binding Modes

Normal human genomic DNA (N-DNA) and mutated DNA (M-DNA) from K562 leukemic cells show different thermodynamic properties and binding affinities on interaction with anticancer drugs; adriamycin (ADR) and daunomycin (DNM). Isothermal calorimetric thermograms representing titration of ADR/DNM with N-DNA and M-DNA on analysis best fitted with sequential model of four and three events respectively. From Raman spectroscopy it has been identified that M-DNA is partially transformed to A form owing to mutations and N-DNA on binding of drugs too undergoes transition to A form of DNA. A correlation of thermodynamic contribution and structural data reveal the presence of different binding events in drug and DNA interactions. These events are assumed to be representative of minor groove complexation, reorientation of the drug in the complex, DNA deformation to accommodate the drugs and finally intercalation. Dynamic light scattering and zeta potential data also support differences in structure and mode of binding of N and M DNA. This study highlights that mutations can manifest structural changes in DNA, which may influence the binding efficacy of the drugs. New generation of drugs can be designed which recognize the difference in DNA structure in the cancerous cells instead of their biochemical manifestation.     

Thermodynamics of protein cross-links

The thermal transitions of native lysozyme and a well-characterized cross-linked derivative of lysozyme [Imoto, T., and Rupley, J. A. (1973), J. Mol. Biol. 80, 657] have been studied in 1.94 M guanidine hydrochloride at pH 2. The observed increase in the melting temperature from 32.4 degrees C for native lysozyme to 61.8 degrees C for the cross-linked derivative corresponds to a calculated 5.2 kcal/mol increase in the free energy of denaturation. This free-energy change is attributed to the decreased entropy of the unfolded polypeptide chain following introduction of a cross-link and is shown to compare well with theoretical predictions. The possibility that an introduction of a cross-link could also affect the enthalpy of an unfolded protein was investigated. The heats of reduction of bovine serum albumin and lysozyme by dithioerythritol in 6 M guanidine hydrochloride were determined and compared to that for the model peptide, oxidized glutathione. The near identity of the observed heats was taken as evidence that the introduction of cross-links into a random-coil protein does not, in general, introduce strain.     

Thermodynamics of phospholipid self-assembly

Negatively charged phospholipids are an important component of biological membranes. The thermodynamic parameters governing self-assembly of anionic phospholipids are deduced here from isothermal titration calorimetry. Heats of demicellization were determined for dioctanoyl phosphatidylglycerol (PG) and phosphatidylserine (PS) at different ionic strengths, and for dioctanoyl phosphatidic acid at different pH values. The large heat capacity (ΔC^o_P ∼ −400 J.mol⁻¹ K⁻¹ for PG and PS), and zero enthalpy at a characteristic temperature near the physiological range (T^∗ ∼ 300 K for PG and PS), demonstrate that the driving force for self-assembly is the hydrophobic effect. The pH and ionic-strength dependences indicate that the principal electrostatic contribution to self-assembly comes from the entropy associated with the electrostatic double layer, in agreement with theoretical predictions. These measurements help define the thermodynamic effects of anionic lipids on biomembrane stability.     

Thermodynamics of primary photosynthesis

The thermodynamics of photosynthesis has been much discussed, but recent articles have pointed to some confusion on the subject. The aim of this review is to clarify a limited part of this state of affairs.     

Flexibility Of DNA In 2:1 Drug-DNA Complexes - Simultaneous Binding Of Two DAPI Molecules To DNA

Abstract

Simultaneous binding of two DAPI molecules in the minor groove of (dA)₁₅.(dT)₁₅ B-DNA helix has been simulated by molecular mechanics calculations. The energy minimised structure shows some novel features in relation to binding of DAPI molecules as well as the flexibility of the grooves of DNA helices. The minor groove of the helix expands locally considerably (to 15 Å) to accommodate the two DAPI molecules and is achieved by positive propeller twisting of base pairs at the binding site concomitant with small variations in the local nucleotide stereochemistry. The expansion also brings forth simultaneously a contraction in the width of the major groove spread over to a few phosphates. These findings demonstrate another facet of the flexible stereochemistry of DNA helices in which the local features are significantly altered without being propagated beyond a few base pairs, and with the rest of the regions retaining the normal structure. Both the DAPI molecules are engaged in specific hydrogen bonds with the bases and non specific interactions with phosphates. Stacking interactions of DAPI molecules between themselves as well as with sugar-phosphate backbone contribute to the stability of the complex. The studies provide a stereochemical support to the experimental findings that under high drug-DNA ratio DAPI could bind in the 2:1 ratio.     

Thermodynamics of aminoglycoside-rRNA recognition

2-Deoxystreptamine (2-DOS) aminoglycosides are a family of structurally related broad-spectrum antibiotics that are used widely in the treatment of infections caused by aerobic Gram-negative bacilli. Their antibiotic activities are ascribed to their abilities to bind a highly conserved A site in the 16 S rRNA of the 30 S ribosomal subunit and interfere with protein synthesis. The abilities of the 2-DOS aminoglycosides to recognize a specific subdomain of a large RNA molecule make these compounds archetypical models for RNA-targeting drugs. This article presents a series of calorimetric, spectroscopic, osmotic stress, and computational studies designed to evaluate the thermodynamics (DeltaG, DeltaH, DeltaS, DeltaCp) of aminoglycoside-rRNA interactions, as well as the hydration changes that accompany these interactions. In conjunction with the current structural database, the results of these studies provide important insights into the molecular forces that dictate and control the rRNA binding affinities and specificities of the aminoglycosides. Significantly, identification of these molecular driving forces [which include binding-linked drug protonation reactions, polyelectrolyte contributions from counterion release, conformational changes, hydration effects, and molecular interactions (e.g., hydrogen bonds and van der Waals interactions)], as well as the relative magnitudes of their contributions to the binding free energy, could not be achieved by consideration of structural data alone, highlighting the importance of acquiring both thermodynamic and structural information for developing a complete understanding of the drug-RNA binding process. The results presented here begin to establish a database that can be used to predict, over a range of conditions, the relative affinity of a given aminoglycoside or aminoglycoside mimetic for a targeted RNA site vs binding to potential competing secondary sites. This type of predictive capability is essential for establishment of a rational design approach to the development of new RNA-targeted drugs.     

Thermodynamics of drug-DNA interactions

Many anticancer, antibiotic, and antiviral drugs exert their primary biological effects by reversibly interacting with nucleic acids. Therefore, these biomolecules represent a major target in drug development strategies designed to produce next generation therapeutics for diseases such as cancer. In order to improve the clinical efficacy of existing drugs and also to design new ones it is necessary to understand the molecular basis of drug-DNA interactions in structural, thermodynamic, and kinetic detail. The past decade has witnessed an increase in the number of rigorous biophysical studies of drug-DNA systems and considerable knowledge has been gained in the energetics of these binding reactions. This is, in part, due to the increased availability of high-sensitivity calorimetric techniques, which have allowed the thermodynamics of drug-DNA interactions to be probed directly and accurately. The focus of this article is to review thermodynamic approaches to examining drug-DNA recognition. Specifically, an overview of a recently developed method of analysis that dissects the binding free energy of these reactions into five component terms is presented. The results of applying this analysis to the DNA binding interactions of both minor groove drugs and intercalators are discussed. The solvent water plays a key role in nucleic acid structure and consequently in the binding of ligands to these biomolecules. Any rational approach to DNA-targeted drug design requires an understanding of how water participates in recognition and binding events. Recent studies examining hydration changes that accompany DNA binding by intercalators will be reviewed. Finally some aspects of cooperativity in drug-DNA interactions are described and the importance of considering cooperative effects when examining these reactions is highlighted.     

Thermodynamics of lipid-peptide interactions

This review is focused on peptide molecules which exhibit a limited solubility in the aqueous phase and bind to the lipid membrane from the aqueous medium. Surface adsorption, membrane insertion, and specific binding are usually accompanied by changes in the heat content of the system and can be measured conveniently with isothermal titration calorimetry, avoiding the necessity of peptide labeling. The driving forces for peptide adsorption and binding are hydrophobicity, electrostatics, and hydrogen bonding. An exclusively hydrophobic interaction is exemplified by the immunosuppressant drug cyclosporine A. Its insertion into the membrane can be described by a simple partition equilibrium X(b)=K(0)C(eq). If peptide and membrane are both charged, electrostatic interactions are dominant leading to nonlinear binding curves. The concentration of the peptide near the membrane interface can then be much larger than its bulk concentration. Electrostatic effects must be accounted for by means of the Gouy-Chapman theory before conventional binding models can be applied. A small number of peptides and proteins bind with very high affinity to a specific lipid species only. This is illustrated for the lantibiotic cinnamycin (Ro 09-0198) which forms a 1:1 complex with phosphatidyethanolamine with a binding constant of 10(8) M(-1). Membrane adsorption and insertion can be accompanied by conformational transitions facilitated, in part, by hydrogen bonding mechanisms. The two membrane-induced conformational changes to be discussed are the random coil-to-alpha-helix transition of amphipathic peptides and the random coil-to-beta-structure transition of Alzheimer peptides.