The effect of intercalator structure on binding strength and base-pair specificity in DNA interactions

The interaction of naphthothiophene, phenanthrene and anthracene ring systems, which have amide and ester side chains with cationic groups (synthesized from the aromatic acid chlorides and appropriate amines and alcohols), with calf thymus DNA has been investigated by using viscometric titrations, spectrophotometric binding experiments and 1H-, 31P- and 17O-NMR methods. The viscosity and NMR experiments suggest that all of these compounds bind to DNA by intercalation. These experiments and spectrophotometric binding studies, however, indicate that there is considerable variation in the interaction of these compounds with DNA. These variations can all be explained by the geometry of the ring systems, the position of protons adjacent to the side chains, and the relative sizes of the amide and ester side chains. With the naphthothiophene ester and amide, for example, the planar amide cannot rotate into the plane of the naphthothiophene ring whereas the smaller planar ester can. With this ring system the ester has a significantly higher binding constant than the amide derivative. Additional binding studies with poly[d(A-T)2] and poly[d(G-C)2] have shown that all of these compounds bind more strongly to the A-T- than the G-C-containing polymer. Since the ester compounds do not have hydrogen bond donating groups proximate to the aromatic ring, these results suggest a model for the A-T specificity of these compounds that involves a solvent-mediated hydrogen bond between the C-2 carbonyl of thymine and the carbonyl group of the intercalators.     

Kinetic and equilibrium analysis of a threading intercalation mode: DNA sequence and ion effects

The interaction of a symmetric naphthalene diimide with alkylamino substituents at each imide position was investigated with the alternating sequence polymers, poly[d(A-T)]2 and poly[d(G-C)]2. Spectrophotometric binding studies indicate strong binding of the diimide to both sequences although the GC binding constant is 20-25 times larger than the AT binding constant. Analysis of the effects of salt concentration on the binding equilibria shows that the diimide forms two ion pairs in its complex with both polymers as expected for a simple dication. Stopped-flow kinetics experiments demonstrate that the diimide both associates and dissociates from DNA more slowly than classical intercalators with similar binding constants. Analysis of salt concentration effects on dissociation kinetics rate constants (kd) reveals that slopes in log kd versus log [Na+] plots are only approximately half the value obtained for classical dicationic intercalators that have both charged groups in the same groove. These kinetics results support a threading intercalation model, with one charged diimide substituent in each of the DNA grooves rather than with both side chains in the same groove, for the diimide complex with DNA. In the rate-determining step of the mechanism for dissociation of a threading complex only one ion pair is broken; the free side chain can then slide between base pairs to put both diimide side chains in the same groove, and this is followed by rapid full dissociation of the diimide. This sequential release of ion pairs makes the dissociation slope for dicationic threading intercalators more similar to the slope for classical monocationic intercalating ligands.(ABSTRACT TRUNCATED AT 250 WORDS)     

Evaluating the Effect of Ionic Strength on Duplex Stability for PNA Having Negatively or Positively Charged Side Chains

The enhanced thermodynamic stability of PNA:DNA and PNA:RNA duplexes compared with DNA:DNA and DNA:RNA duplexes has been attributed in part to the lack of electrostatic repulsion between the uncharged PNA backbone and negatively charged DNA or RNA backbone. However, there are no previously reported studies that systematically evaluate the effect of ionic strength on duplex stability for PNA having a charged backbone. Here we investigate the role of charge repulsion in PNA binding by synthesizing PNA strands having negatively or positively charged side chains, then measuring their duplex stability with DNA or RNA at varying salt concentrations. At low salt concentrations, positively charged PNA binds more strongly to DNA and RNA than does negatively charged PNA. However, at medium to high salt concentrations, this trend is reversed, and negatively charged PNA shows higher affinity for DNA and RNA than does positively charged PNA. These results show that charge screening by counterions in solution enables negatively charged side chains to be incorporated into the PNA backbone without reducing duplex stability with DNA and RNA. This research provides new insight into the role of electrostatics in PNA binding, and demonstrates that introduction of negatively charged side chains is not significantly detrimental to PNA binding affinity at physiological ionic strength. The ability to incorporate negative charge without sacrificing binding affinity is anticipated to enable the development of PNA therapeutics that take advantage of both the inherent benefits of PNA and the multitude of charge-based delivery technologies currently being developed for DNA and RNA.     

NMR characterization of the electrostatic interaction of the basic residues in HDGF and FGF2 during heparin binding

Electrostatic interaction is a major driving force in the binding of proteins to highly acidic glycosaminoglycan, such as heparin. Although NMR backbone chemical shifts have generally been used to identify the heparin-binding site on a protein, however, there is no correlation between the binding free energies and the perturbed backbone chemical shifts for individual residues. The binding event occurs at the end of a side chain of basic residue, and does not require causing significant alterations in the backbone environment at a distance of multiple bonds. We used the H₂CN NMR pulse sequence to detect heparin binding through the side-chain resonances Hε–Cε–Nζ of Lys and Hδ–Cδ–Nε of Arg in the two proteins of hepatoma-derived growth factor (HDGF) and basic fibroblast growth factor (FGF2). H₂CN titration experiments revealed chemical shift perturbations in the side chains, which were correlated with the free energy changes in various mutants. The residues K19 in HDGF and K125 in FGF2 demonstrated the most significant perturbations, consistent with our previous observation that the two residues are crucial for binding. The result suggests that H₂CN NMR provides a precise evaluation for the electrostatic interactions. The discrepancy observed between backbone and side chain chemical shifts is correlated to the solvent accessibility of residues that the K19 and K125 backbones are highly buried with the restricted backbone conformation and are not strongly affected by the events at the end of the side chains.     

Differential DNA and RNA sequence discrimination by PNA having charged side chains

PNA sequences modified with charged side chains were evaluated for base-pairing sequence selectivity under physiological conditions. PNA having negatively charged aspartic acid side chains shows higher selectivity with RNA, while PNA having positively charged lysine side chains shows higher selectivity with DNA. These observations provide insight into the binding selectivity of modified PNA in antisense and antigene applications.     

Novel antitumor agent family of 1H-benzo[c,d]indol-2-one with flexible basic side chains: synthesis and biological evaluation

A series of mono-1H-benzo[c,d]indol-2-one with different amine side chains and bis-1H-benzo[c,d]indol-2-one as novel family of DNA intercalators were designed and synthesized, the contributions of aromatic chromophores and amine side chains for DNA binding properties, for example, intercalation and electrostatic binding, respectively, were evaluated. Among them, A3 tailed with N,N-dimethylamino-ethyl-ethane-1,2-diamine showed selective anti-tumor activities against cell lines A549 and P388 with IC(50) 0.428microm and 1.69microm.     

Binding of DNA to large fragment of DNA polymerase I: identification of strong and weak electrostatic forces and their biological implications.

Examination of the electrostatic potential of a modeled complex, consisting of the Klenow fragment of E. coli DNA polymerase I and DNA template-primer, suggested the presence of two distinct interacting regions. The one displaying a strong electropositive potential field is generated by side chains of basic amino acid pairs and is directed towards the major groove site in DNA. The second electrostatic potential field around DNA is somewhat weaker and appears to be exerted by a pair of vicinal side chains of acidic and basic amino acids. The distribution of charges in this manner appears well suited for the binding of enzyme to the template-primer required in the enzymatic synthesis of DNA.     

Structural basis for the binding affinity of a homologous series of synthetic phenoxazone drugs with DNA: NMR and molecular mechanics analysis

The molecular basis of the marked structure-activity relationship for a homologous series of DNA-binding phenoxazone drugs (ActII-ActIV) has been investigated by NMR spectroscopy and molecular mechanics. The spatial structures of the complexes between the drugs and a model deoxytetranucleotide, 5'-d(TpGpCpA), have been determined by molecular mechanics methods using homonuclear (1)H-(1)H 2D-NOESY and heteronuclear (1)H-(31)P (HMBC) NMR spectroscopic data. Observed intermolecular NOE contacts and equilibrium binding studies confirm that the binding affinity of the synthetic phenoxazone derivatives with d(TGCA) decreases with an increase in the number of CH(2) groups in the dimethylaminoalkyl side chains, i.e., ActII > ActIII > ActIV, in agreement with the observed biological activity of these compounds. Molecular mechanics calculations of the spatial structures of the intercalated complexes of ActII-ActIV with d(TGCA) indicate that the different binding constants of the phenoxazone derivatives with the DNA oligomer are due to the different degrees of intercalation of the chromophore and the different steric arrangements of aminoalkyl side chains in the minor groove of the tetramer duplex; this results in different distances between the negatively-charged phosphates of the DNA duplex and the terminal positively-charged N(CH(3))(2) groups of the side chains.     

Side-chain flexibility in proteins upon ligand binding

Ligand binding may involve a wide range of structural changes in the receptor protein, from hinge movement of entire domains to small side-chain rearrangements in the binding pocket residues. The analysis of side chain flexibility gives insights valuable to improve docking algorithms and can provide an index of amino-acid side-chain flexibility potentially useful in molecular biology and protein engineering studies. In this study we analyzed side-chain rearrangements upon ligand binding. We constructed two non-redundant databases (980 and 353 entries) of "paired" protein structures in complexed (holo-protein) and uncomplexed (apo-protein) forms from the PDB macromolecular structural database. The number and identity of binding pocket residues that undergo side-chain conformational changes were determined. We show that, in general, only a small number of residues in the pocket undergo such changes (e.g., approximately 85% of cases show changes in three residues or less). The flexibility scale has the following order: Lys > Arg, Gln, Met > Glu, Ile, Leu > Asn, Thr, Val, Tyr, Ser, His, Asp > Cys, Trp, Phe; thus, Lys side chains in binding pockets flex 25 times more often then do the Phe side chains. Normalizing for the number of flexible dihedral bonds in each amino acid attenuates the scale somewhat, however, the clear trend of large, polar amino acids being more flexible in the pocket than aromatic ones remains. We found no correlation between backbone movement of a residue upon ligand binding and the flexibility of its side chain. These results are relevant to 1. Reduction of search space in docking algorithms by inclusion of side-chain flexibility for a limited number of binding pocket residues; and 2. Utilization of the amino acid flexibility scale in protein engineering studies to alter the flexibility of binding pockets.     

Naphthalimide intercalators with chiral amino side chains: effects of chirality on DNA binding, photodamage and antitumor cytotoxicity

Several novel heterocyclic-fused naphthalimides intercalators with chiral amino side chains were investigated. Their side chains' chiral configuration determines DNA binding activities in the order: S-enantiomers > R-enantiomers. And their DNA photodamaging activities were in good agreement with their DNA binding constants, the S-enantiomers could photocleave circular supercoiled pBR322 DNA more efficiently than their R-enantiomers. S-enantiomer B(3) could photodamage DNA at 0.2 microM and cleave supercoiled plasmid DNA from form I to form II completely at 50 microM. Almost all of these intercalators showed effective cytoxicities against human lung cancer cells and murine leukemia cells. S-enantiomers showed different antitumor cytotoxicity by comparison with R-enantiomers. This work may provide additional information for the role of amino side chains on intercalators as antitumor agents.     

Xeroderma pigmentosum complementation group A protein is driven to nucleotide excision repair sites by the electrostatic potential of distorted DNA

The presumed DNA-binding cleft of xeroderma pigmentosum group A (XPA) protein, a key regulatory subunit of the eukaryotic nucleotide excision repair complex, displays a distinctive array of 6 positively charged amino acid side chains. Here, the molecular function of these closely spaced electropositive residues has been tested by systematic site-directed mutagenesis. After the introduction of single amino acid substitutions, the mutants were probed for protein-DNA interactions in electrophoretic mobility shift and photochemical crosslinking assays. This analysis led to the identification of a critical hot-spot for DNA substrate recognition composed of two neighboring lysines at codons 141 and 179 of the human XPA sequence. The replacement of other basic side chains in the DNA interaction domain conferred more moderate defects of substrate binding. When the function of XPA was tested as a fusion product with either mCherry or green-fluorescent protein, a glutamate substitution of one of the positively charged residues at positions 141 and 179 was sufficient to decrease DNA repair activity in human fibroblasts. Thus, the removal of a single cationic side chain abolished DNA-binding activity and significant excision repair defects could be induced by single charge inversions on the XPA surface, indicating that this molecular sensor participates in substrate recognition by monitoring the electrostatic potential of distorted DNA repair sites.     

DNA binding properties of the yeast Msh2-Msh6 and Mlh1-Pms1 heterodimers

We describe here our recent studies of the DNA binding properties of Msh2-Msh6 and Mlh1-Pms1, two protein complexes required to repair mismatches generated during DNA replication. Mismatched DNA binding by Msh2-Msh6 was probed by mutagenesis based on the crystal structure of the homologous bacterial MutS homodimer bound to DNA. The results suggest that several amino acid side chains inferred to interact with the DNA backbone near the mismatch are critical for repair activity. These contacts, which are different in Msh2 and Msh6, likely facilitate stacking and hydrogen bonding interactions between side chains in Msh6 and the mismatched base, thus stabilizing a kinked DNA conformation that permits subsequent repair steps coordinated by the Mlh1-Pms1 heterodimer. Mlh1-Pms1 also binds to DNA, but independently of a mismatch. Mlh1-Pms1 binds short DNA substrates with low affinity and with a slight preference for single-stranded DNA. It also binds longer duplex DNA molecules, but with a higher affinity indicative of cooperative binding. Indeed, imaging by atomic force microscopy reveals cooperative DNA binding and simultaneous interaction with two DNA duplexes. The novel DNA binding properties of Mlh1-Pms1 may be relevant to signal transduction during DNA mismatch repair and to recombination, meiosis and cellular responses to DNA damage.     

Identification of functionally important residues in the DNA binding region of the mnt repressor

Single amino acid substitutions have been introduced throughout the N-terminal DNA binding region of the Mnt repressor, and the operator binding properties of the resulting mutant repressors have been assayed. These studies show that the side chains of Arg2, His6, Asn8, and Arg10 are critical for high affinity binding to operator DNA. Other side chains in the N-terminal region do not appear to play major roles in DNA recognition and binding. Specific alterations in the pattern of methylation protection afforded by the Arg2----Lys mutant protein suggest that Arg2 contacts the N7 groups of guanines 10 and 12 in the operator. In conjunction with previous results, these findings suggest that part of the N-terminal region of Mnt binds as an extended polypeptide strand within the major groove of the mnt operator.     

Structural Basis for the Binding Affinity of a Homologous Series of Synthetic Phenoxazone Drugs with DNA: NMR and Molecular Mechanics Analysis

Abstract

The molecular basis of the marked structure-activity relationship for a homologous series of DNA-binding phenoxazone drugs (ActII-ActIV) has been investigated by NMR spectroscopy and molecular mechanics. The spatial structures of the complexes between the drugs and a model deoxytetranucleotide, 5′-d(TpGpCpA), have been determined by molecular mechanics methods using homonuclear ¹H-¹H 2D-NOESY and heteronuclear ¹H-³¹P (HMBC) NMR spectroscopic data. Observed intermolecular NOE contacts and equilibrium binding studies confirm that the binding affinity of the synthetic phenoxazone derivatives with d(TGCA) decreases with an increase in the number of CH₂ groups in the dimethylami- noalkyl side chains, i.e., ActII > ActIII > ActIV, in agreement with the observed biological activity of these compounds. Molecular mechanics calculations of the spatial structures of the intercalated complexes of ActII-ActIV with d(TGCA) indicate that the different binding constants of the phenoxazone derivatives with the DNA oligomer are due to the different degrees of intercalation of the chromophore and the different steric arrangements of aminoalkyl side chains in the minor groove of the tetramer duplex; this results in different distances between the negatively-charged phosphates of the DNA duplex and the terminal positively-charged N(CH₃)₂ groups of the side chains.     

Effect of mutation of the Sac7d intercalating residues on the temperature dependence of DNA distortion and binding thermodynamics

Sac7d is a small chromatin protein from the hyperthermophile Sulfolobus acidocaldarius which kinks duplex DNA by approximately 66 degrees at a single base pair step with intercalation of V26 and M29 side chains. Site-directed mutagenesis coupled with calorimetric and spectroscopic data has been used to characterize the influence of the intercalating side chains on the structure and thermodynamics of the DNA complex from 5 to 85 degrees C. Two single-alanine substitutions (V26A and M29A) and five double-glycine, -alanine, -leucine, -phenylalanine, and -tryptophan substitutions of the surface residues have been created. NMR and fluorescence titrations indicated that the substitutions had little effect on the structure of the protein or DNA binding site size. Each of the mutant proteins demonstrated a temperature-dependent binding enthalpy which was correlated with a similar temperature dependence in the structure of the complex reflected by changes in fluorescence and circular dichroism. A positive heat capacity change (DeltaC(p)) for DNA binding was observed for only those mutants which also demonstrated a thermotropic structural transition in the complex, and the temperature range for the positive DeltaC(p) coincided with that observed for the structural transition. The thermodynamic data are interpreted using a model in which binding is linked to an endothermic distortion of the DNA in the complex. The results support the proposal that the unfavorable enthalpy of binding of Sac7d at 25 degrees C is due in part to the distortion of DNA.     

Linear diffusion of the restriction endonuclease EcoRV on DNA is essential for the in vivo function of the enzyme.

Linear diffusion along DNA is a mechanism of enhancing the association rates of proteins to their specific recognition sites on DNA. It has been demonstrated for several proteins in vitro, but to date in no case in vivo. Here we show that the restriction endonuclease EcoRV slides along the DNA, scanning approximately 1000 bp in one binding event. This process is critically dependent on contacts between amino acid residues of the protein and the backbone of the DNA. The disruption of single hydrogen bonds and, in particular, the alteration of electrostatic interactions between amino acid side chains of the protein and phosphate groups of the DNA interfere with or abolish effective sliding. The efficiency of linear diffusion is dependent on salt concentration, having a maximum at 50 mM NaCl. These results suggest that a nonspecific and mobile binding mode capable of linear diffusion is dependent on a subtle balance of forces governing the interaction of the enzyme and the DNA. A strong correlation between the ability of EcoRV mutants to slide along the DNA in vitro and to protect Escherichia coli cells from phage infection demonstrates that linear diffusion occurs in vivo and is essential for effective phage restriction.     

Interaction of a macrocyclic bisacridine with DNA

The binding of the macrocycle SDM to DNA was investigated by visible spectroscopy, stopped-flow kinetics, and NMR spectroscopy. SDM is composed of two 9-aminoacridines linked via the amino groups by a spermine side chain and via the 4-positions by a N,N'-[(methylthio)ethyl]succinamide side chain [Zimmerman, S. C., Lamberson, C. R., Cory, M., & Fairley, T. A. (1989) J. Am. Chem. Soc. 111, 6805-6809]. The visible spectrum of SDM bound to poly[d(A-T)]2 or poly[d(G-C)]2 is red-shifted relative to the spectrum of SDM alone and displays considerable hypochromicity. Results from titrations of SDM with polymer indicate a binding site size of three base pairs per macrocycle. The dissociation constant for SDM bound to either poly[d(A-T)]2 or poly[d(G-C)]2 is an order of magnitude lower than that for a similar bisacridine linked only by a spermine side chain. In addition, the dependence of the dissociation constant on ionic strength is significantly reduced. NMR studies of SDM complexes with poly[d(A-T)]2 or a tetramer, d(CGCG)2, show that intercalation is the mode of binding. The magnitudes of the chemical shift differences for SDM aromatic protons in the free and bound states support intercalation with the acridine ring systems essentially parallel to the long axis of the base pairs. Cross peaks from NOESY spectra of the SDM complex with d(CGCG)2 further support this mode of binding and provide information on the structure of the complex. The results are analyzed for consistency with each of three binding models: (i) bisintercalation with the two side chains in the same groove; (ii) bisintercalation according to the neighbor-exclusion principle with the two side chains in opposite grooves; and (iii) bisintercalation with two side chains in opposite grooves but with violation of the neighbor-exclusion principle. Model i is found to be unlikely on the basis of all evidence obtained, including preliminary modeling studies. Both models ii and iii can be reconciled with the experimental evidence and from a modeling standpoint are energetically feasible.     

Interaction of DNA with connective tissue matrix proteins reveals preferential binding to type V collagen

The interaction of DNA with type I to VI collagens and laminin was studied in vitro in systems in which the connective tissue components were immobilized, as well as when in solution. In studies on immobilized components, significant binding of DNA was observed only for type V collagen, and the binding of radiolabeled DNA to this component could be effectively inhibited in a concentration-dependent manner by the addition of unlabeled DNA. Similar results were observed in solution assays in which it was observed that DNA binding to type V collagen was dependent on the native triple-helical conformation of the collagen. The preferential binding of DNA to native type V collagen may be due to the relative basicity of type V collagen chains, as well as the unique spatial arrangement of amino acid side chains in the native molecules. The data are of potential clinical relevance in that binding of DNA to type V collagen may represent at least one component of the mechanism whereby DNA and its immune complexes are deposited in connective tissues in certain pathologic conditions.     

Intercalating residues determine the mode of HMG1 domains A and B binding to cisplatin-modified DNA

Cisplatin exerts its anticancer activity by forming covalent adducts with DNA. High-mobility group (HMG)-domain proteins recognize the major 1,2-intrastrand cisplatin-DNA cross-links and can mediate cisplatin cytotoxicity. The crystal structure of HMG1 domain A bound to cisplatin-modified DNA, further analyzed here, reveals intercalation of a key Phe37 residue. Other published structures of HMG domains bound to DNA, including NHP6A and HMG-D, similarly indicate amino acid side chains intercalating into linear DNA to form a bend. To delineate the importance of such side chain intercalations and further to explore the binding modes of different HMG domains toward prebent DNA structures, site-directed mutagenesis was used to generate HMG1 domain A and domain B mutants. The affinities of these mutant proteins for cisplatin-modified DNA were determined in gel electrophoresis mobility shift assays. The results indicate that intercalating residues at positions 16 or 37 can both contribute to the binding affinity. The data further reveal that the length of the loop between helices I and II is not critical for binding affinity. Footprinting analyses indicate that the position of the intercalating residue dictates the binding mode of the domain toward platinated DNA. Both congruent and offset positioning of the HMG domain with respect to the locus of the cisplatin-induced bend in the DNA were encountered. Packing interactions in the crystal structure suggest how full-length HMG1 might bind to DNA by contacting more than one duplex simultaneously. Taken together, these results demonstrate that cisplatin modification of DNA provides an energetically favorable, prebent target for HMG domains, which bind to these targets through one or more side chain and favorable hydrophobic surface interactions.