DNA sequence-specific recognition of peptides incorporating the HPRK and polyamide motifs

Three peptide amides, HPRK(Py)(4)HPRK-NH(2) (PyH-12), HPRK(Py)(3)HPRK-NH(2) (PyH-11) and HPRK(Py)(2)HPRK-NH(2) (PyH-10), incorporating two HPRK motifs and various 4-amino-1-methylpyrrole-2-carboxylic acid residues (Py) were synthesized by solid-phase peptide methodology. The binding of these three peptides to a 5'-32P-labeled 158-mer DNA duplex (Watson fragment) and to a 5'-32P-labeled 135-mer DNA duplex (complementary Crick fragment) was investigated by quantitative DNase I footprinting. On the 158-mer Watson strand, the most distinctive DNase I blockages seen with all three peptides occur around positions 105-112 and 76-79, corresponding to the sequences 5'-GAGAAAAT-3' and 5'-CGGT-3', respectively. However, on the complementary Crick strand, only PyH-12 strongly discriminates the 5'-TTT-3' site around positions 108-110 whereas both PyH-11 and PyH-10 have moderate binding around positions 102-112 comprising the sequence 5'-ATTTTCTCCTT-3'. Possible bidentate and single interactions of the side-chain functions and alpha-amino protons of the peptides with DNA bases are discussed.     

Search for the optimal linker in tandem hairpin polyamides

In order to target specific DNA sequences >or=10 base pairs in size by minor groove binding ligands, a search for the optimal linker in dimers of hairpin polyamides was initiated. Two series of tandem polyamides ImPyIm-(R)[ImPyIm-(R)(H2N)gamma-PyPyPy-L](HN)gamma-PyPyPy-beta-Dp (1a-e), where L represents a series of 4-8 carbon long aliphatic amino acid linkers, and ImPyIm-(R)[ImPyIm-(R)(H2N)gamma-PyPyPyIm-L](HN)gamma-PyPyPy-beta-Dp (2a-e), where L represents a series of 2-6 carbon long aliphatic amino acid linkers, were synthesized and characterized by quantitative DNase I footprinting. beta, gamma and Dp represents beta-alanine, gamma-aminobutyric acid, and 3-(dimethylamino)propylamine, respectively. It was found that the five-carbon 5-aminovaleric acid (delta), is suitable to span one base-pair (bp) of DNA when incorporated into a tandem polyamide. ImPyIm-(R)[ImPyIm-(R)(H2N)gamma-PyPyPy-delta](HN)gamma-PyPyPy-beta-Dp (1b) binds the 10 bp binding-site 5'-AGTGAAGTGA-3' with equilibrium association constant K(a)=3.2 x 10(10) M(-1) and ImPyIm-(R)[ImPyIm-(R)(H2N)gamma-PyPyPyIm-delta](HN)gamma-PyPyPy-beta-Dp (2d) binds the 11 bp binding-site 5'-AGTGATAGTGA-3' with K(a)=9.7 x 10(9) M(-1). Tandem 1b also bind the 11 bp site but with lower affinity affording a 15-fold specificity for the shorter binding site. Replacing a methylene group in the amino acid linker with an oxygen atom to form tandem polyamide ImPyIm-(R)[ImPyIm-(R)(H2N)gamma-PyPyPy-E](HN)gamma-PyPyPy-beta-Dp (4) where E represents the ether linker, resulted in that an 80-fold specificity for the 10 bp binding site over the 11 bp site.     

Quantitative footprinting analysis of the chromomycin A3--DNA interaction.

Chromomycin A3 (CHR) binding to the duplex d(CAAGTCTGGCCATCAGTC).d(GACTGATGGCCAGACTTG) has been studied using quantitative footprinting methods. Previous NMR studies indicated CHR binds as a dimer in the minor groove. Analysis of autoradiographic spot intensities derived from DNase I cleavage of the 18-mer in the presence of various amounts of CHR revealed that the drug binds as a dimer to the sequence 5'-TGGCCA-3',3'-ACCGGT-5' in the 18-mer with a binding constant of (2.7 +/- 1.4) x 10(7) M-1. Footprinting and fluorescence data indicate that the dimerization constant for the drug in solution is approximately 10(5) M-1. Since it has been suggested that CHR binding alters DNA to the A configuration, quantitative footprinting studies using dimethyl sulfate, which alkylates at N-7 of guanine in the major groove, were also carried out. Apparently, any drug-induced alteration in DNA structure does not affect cleavage by DMS enough to be observed by these experiments.     

Effects of a hairpin polyamide on DNA melting: comparison with distamycin and Hoechst 33258

We have used DNase I footprinting and fluorescence melting studies to study the interaction of the hairpin polyamide Im-Py-Py-Py-(R)H2Ngamma-Im-Py-Py-Py-beta-Dp with its preferred binding sites (5'-WGWWCW; W=A or T) and other sequences. DNase I footprinting confirmed that the ligand binds to the sequence AGAACA at nanomolar concentrations and that changing the terminal A to G causes a dramatic decrease in affinity, while there was no interaction with the reverse sequence WCWWGW. Fluorescence melting studies with 11-mer duplexes showed that the polyamide had very different effects on the forward (TGWWCT) and reverse (TCTAGT) sequences. At low concentrations, the polyamide produced biphasic melting curves with TGATCT, TGTACT and TGAACT, suggesting a strong interaction. In contrast, the melting profiles with TCTAGT were always monophasic and showed much smaller concentration dependent changes in Tm. The polyamide also showed weak binding to the sequence TGATCT when one of the central AT pairs was replaced with an AC mismatch. These melting profiles were compared with those produced by the AT-selective minor groove binding agents distamycin and Hoechst 33258 at the same sites and at similar sequences containing A5 and (AT)3, which are expected to bind distamycin in the 1:1 and 2:1 modes, respectively. These ligands produced simple monophasic melting curves in which the Tm steadily increased as the ligand concentration was raised.     

Site-specific binding constants for actinomycin D on DNA determined from footprinting studies.

We report site-specific binding constants for the intercalating anticancer drug actinomycin D (Act-D), binding to a 139-base-pair restriction fragment from pBR 322 DNA. The binding constants are derived from analysis of footprinting experiments, in which the radiolabeled 139-mer is cleaved using DNase I, the cleavage products undergo gel electrophoresis, and, from the gel autoradiogram, spot intensities, proportional to amounts of cleaved fragments, are measured. A bound drug prevents DNase I from cleaving at approximately 7 bonds, leading to decreased amounts of corresponding fragments. With the radiolabel on the 3' end of the noncoding strand (A-label), we measured relative amounts of 54 cleavage products at 25 Act-D concentrations. For cleavage of the 139-mer with the label on the 3' end of the coding strand (G-label), relative amounts of 43 cleavage products at 11 Act-D concentrations were measured. These measurements give information about approximately 120 base pairs of the restriction fragment (approximately 12 turns of the DNA helix); in this region, 14 strong and weak Act-D binding sites were identified. The model used to interpret the footprinting plots is derived in detail. Binding constants for 14 sites on the fragment are obtained simultaneously. It is important to take into account the effect of drug binding at its various sites on the local concentration of probe elsewhere. It is also necessary to include in the model weak as well as strong Act-D sites on the carrier DNA which is present, since the carrier DNA controls the free-drug concentration. As expected, the strongest sites are those with the sequence (all sequences are 5'----3') GC, with TGCT having the highest binding constant, 6.4 x 10(6) M-1. Sites having the sequence GC preceded by G are weak binding sites, having binding constants approximately 1 order of magnitude lower than those of the strong sites. Also, the non-GC-containing sequences CCG and CCC bind Act-D with a binding constant comparable to those of the weak GGC sites. The analysis may reveal drug-induced structural changes on the DNA, which are discussed in terms of the mechanism of Act-D binding.     

Hydration pattern of A4T4 and T4A4 DNA: a molecular dynamics study

Hydration pattern and energetics of 'A-tract' containing duplexes have been studied using molecular dynamics on 12-mer self-complementary sequences 5'-d(GCA4T4GC)-3' and 5'-d(CGT4A4CG)-3'. The structural features for the simulated duplexes showed correlation with the corresponding experimental structures. Analysis of the hydration pattern confirmed that water network around the simulated duplexes is more conformation specific rather than sequence specific. The calculated heat capacity change upon duplex formation showed that the process is entropically driven for both the sequences. Furthermore, the theoretical free energy estimates calculated using MMPBSA approach showed a higher net electrostatic contribution for A4T4 duplex formation than for T4A4, however, energetically both the duplexes are observed to be equally stable.     

Detection of multiple network-based allosteric interactions between peptides and arrays of DNA binding sites

Quantitative DNase I footprinting shows that three designed peptides containing N-methylpyrrole (Py) moieties display different types of network-based allosteric communication in binding to DNA: circuit type, incomplete-circuit type, and non-circuit type characterized by interstrand bidentate interactions. Positive cooperative binding of all three peptides to individual DNA binding sites is commonly observed. CD spectral characterization of the interaction between peptides and model undecanucleotide duplexes is consistent with the footprinting results and supports the allosteric model. This study provides insights relating to the interaction network nature of allostery in complex DNA–small molecule interactions.     

Gibberellin treatment stimulates nuclear factor binding to the gibberellin response complex in a barley alpha-amylase promoter.

The promoters of a majority of cereal alpha-amylase genes contain three highly conserved sequences (gibberellin response element, box I, and pyrimidine box). Recent studies have demonstrated the functional importance of four regions that either coincide with or are immediately proximal to these three conserved elements as well as an upstream Opaque-2 binding sequence. In this study, we describe the characterization of nuclear protein factors from barley aleurone layers whose binding activity toward gibberellin response complex sequences from the barley low-pl alpha-amylase gene (Amy32b) promoter is stimulated by gibberellin A3 (GA3) treatment. Barley proteins isolated from crude nuclear extracts prepared from aleurone layers incubated with or without GA3 were fractionated by anion exchange fast protein liquid chromatography and studied using band shift assays, sequence-specific competitions, and DNase I footprinting. A GA3-dependent binding activity eluting at 210 mM KCl was shown to bind specifically to the gibberellin response element and the closely associated box I. DNase I footprinting with the proteins in this fraction indicated interactions with sequences in the gibberellin response element and box I. A second DNA binding activity eluting at 310 mM KCl was present constitutively in extracts prepared from tissues incubated both in the absence and in the presence of hormone. Proteins in this fraction were able to bind to many DNA sequences and, in general, were largely nonspecific. DNase I footprinting with the proteins in this fraction indicated a large area of protection with a single unoccupied region located at the 3' end of box I. The possible function of such an activity in hormone regulation of the alpha-amylase genes is discussed.     

Efficient conjugation and characterization of distamycin-based peptides with selected oligonucleotide stretches

Selected sequences of oligodeoxyribonucleotides (ODNs) have been conjugated efficiently with distamycin-based peptides containing reactive cysteine and oxyamine functionalities at the C-terminus. The conjugation was performed easily within 30-60 min, using individual modified oligonucleotide stretches having sequences of 5'-d(GCTTTTTTCG)-3', 5'-d(GCTATATACG)-3', and 5'-AGCGCGCGCA-3'. Two types of linkages were used for making the covalent connection: (i) a five-membered thiazolidine ring and (ii) an oxime. These distamycin-like polyamide-ODN conjugates were then converted to the corresponding DNA duplexes using complementary oligonucleotide sequences. To elucidate the binding specificity of the distamycin-oligonucleotide conjugates, UV-melting temperature measurements were performed. These studies indicated that the distamycin-ODN conjugate favored binding with the duplex with sequence 5'-d(GCTTTTTTCG)-3' rather than 5'-d(GCTATATACG)-3'. On the other hand, no stabilization of the duplex with sequence 5'-d(AGCGCGCGCA)-3' was observed. UV results also suggest that the thiazolidine and oxime linkages do not significantly influence the process of distamycin binding to the minor groove surface of the DNA duplex. The results obtained from duplex UV-melting studies were further corroborated by a temperature-dependent study of the circular dichroism spectra of the conjugates and a fluorescence displacement titration assay using Hoechst 33258 fluorophore as a competitive binder for the minor groove. All these studies reinforce the fact that the specific stabilization of A/T rich DNA-DNA duplexes by distamycin was preserved upon conjugation with oligonucleotide stretches.     

Binding of distamycin A and netropsin to the 12mer DNA duplexes containing mixed AT.GC sequences with at most five or three successive AT base pairs

Circular dichroism (CD), isothermal calorimetric titrations (ITC), and temperature-dependent UV spectroscopy were used to investigate binding of the minor groove-directed ligands distamycin A (Dst) and netropsin (Net) to the following duplexes: d(GTTAGTATTTGG). d(CCAAATACTAAC), d(GTTAGTATATGG).d(CCATATACTAAC), d(GTTAGTACTTGG). d(CCAAGTACTAAC), and d(GTTAGTAGTTGG).d(CCAACTACTAAC). Our results reveal that Dst binds within the minor grooves of these dodecamers that contain five-AT and/or four-AT.GC binding sites exclusively in a dimeric high-affinity 2:1 binding mode (K approximately 10(16) M(-)(2)). By contrast, Net exhibits high-affinity binding only when it binds in a 1:1 mode (K(1) approximately 10(9) M(-)(1)) to the two duplexes that contain five-AT sites (5'-TATTT-3' and 5'-TATAT-3'). Its further binding to these two duplexes occurs in a low-affinity mode (K(2) approximately 10(6) M(-)(1)) and results in the formation of 2:1 Net-DNA complexes. To the other two duplexes that contain sequences with at most three AT consecutive base pairs Net binds in two distinctive low-affinity 1:1 binding modes (K(1) approximately 10(7) M(-)(1), K(2) approximately 10(6) M(-)(1)). Competition experiments (CD and ITC titrations) reveal that Dst entirely displaces Net from its 1:1 and 2:1 complexes with any of the four duplexes. We discuss and interpret our optical and calorimetric results in the context of the available structural information about the complexes between DNA and the sequence-specific minor groove binders Dst and Net.     

Deoxyribonuclease I footprinting reveals different DNA binding modes of bifunctional platinum complexes

Deoxyribonuclease I (DNase I) footprinting methodology was used to analyze oligodeoxyribonucleotide duplexes containing unique and single, site-specific adducts of trinuclear bifunctional platinum compound, [{trans-PtCl(NH3)2}2 mu-trans-Pt(NH3)2{H2N(CH2)6NH2}2]4+ (BBR3464) and the results were compared with DNase I footprints of some adducts of conventional mononuclear cis-diamminedichloroplatinum(II) (cisplatin). These examinations took into account the fact that the local conformation of the DNA at the sites of the contacts of DNase I with DNA phosphates, such as the minor groove width and depth, sequence-dependent flexibility and bendability of the double helix, are important determinants of sequence-dependent binding to and cutting of DNA by DNase I. It was shown that various conformational perturbations induced by platinum binding in the major groove translated into the minor groove, allowing their detection by DNase I probing. The results also demonstrate the very high sensitivity of DNase I to DNA conformational alterations induced by platinum complexes so that the platinum adducts which induce specific local conformational alterations in DNA are differently recognized by DNase I.     

Cooperative binding of oligonucleotides to adjacent sites of single-stranded DNA: sequence composition dependence at the junction

The quantitative parameters of cooperative binding of deoxyribooligonucleotides to adjacent sites by double helix formation have been determined as a function of sequence composition at the junction. The base stacks 5'-Py/p-Py-3', 5'-Pu/p-Py-3' and 5'-Pu/p-Pu-3' (p is phosphate group, Py and Pu are pyrimidine and purine nucleoside, respectively) including mismatches on the 3'-side of the junction were studied using complementary addressed modification titration (CAMT) at 25 degrees C and pH 7.5, 0.16 M NaCl, 0.02 M Na2HPO4, 0.1 mM EDTA. The equilibrium binding constants of alkylating derivatives of 8-mer oligonucleotides (reagents) with 22-mer oligonucleotides (targets) were determined using the dependence of the target limit modification extents on the concentrations of the reagents. The parameters of cooperativity were calculated as the ratio of binding constants of reagents in the presence and the absence of a second 8-mer oligonucleotides (effectors) occupying the adjacent site on the 22-mer targets. For the stacks 5'-Py/p-Py-3' the parameters of cooperativity were around unity both for matched and mismatched nucleotides at the junction indicating the absence of cooperativity. The parameters of cooperativity for the stacks 5'-Pu/p-Pu-3' were higher than for the stacks 5'-Pu/p-Py-3' in perfect and non-perfect duplexes. Discrimination of mismatches was higher in nicked than in normal duplexes.     

Structural constraints for the binding of short peptides to claudin-4 revealed by surface plasmon resonance

Claudin family transmembrane proteins play an important role in tight junction structure and function in epithelial cells. Among the 24 isoforms identified in mice and humans, claudin-4 and -3 serve as the receptor for Clostridium perfringens enterotoxin (Cpe). The second extracellular loop (Ecl2) of claudin-4 is responsible for the binding to the C-terminal 30 amino acids of Cpe (Cpe30). To define the structural constraints for the claudin-4/Cpe30 interaction, a surface plasmon resonance (SPR) method was developed. GST fusions with claudin-4 revealed that Ecl2 with the downstream transmembrane domain of claudin-4 reconstituted the basic structural requirement for optimal binding activity to Cpe30, with affinity in the nanomolar range. Two 12-mer peptides selected by phage display against claudin-4-transfected CHO cells and a 12-mer Cpe mutant peptide also showed significant affinity for claudin-4 with this SPR assay, suggesting that a short peptide can establish stable contact with Ecl2 with nanomolar affinity. Alignment of these short peptides unveiled a common Ecl2 binding motif: . Whereas the short peptides bound native claudin-4 on transfected CHO cells in pull-down assays, only the larger Cpe30 peptide affected trans-epithelial electrical resistance (TER) in peptide-treated Caco-2BBe monolayers. Importantly, Cpe30 retained its binding to claudin-4 when fused to the C terminus of influenza hemagglutinin, demonstrating that its binding activity can be maintained in a different biochemical context. These studies may help in the design of assays for membrane receptor interactions with soluble ligands, and in applying new targeting ligands to delivering attached "cargo" proteins.     

Binding of f-PIP, a pyrrole- and imidazole-containing triamide, to the inverted CCAAT box-2 of the topoisomerase IIalpha promoter and modulation of gene expression in cells

An N-formamido pyrrole- and imidazole-containing triamide (f-PIP) has been shown by DNase I footprinting, SPR, and CD studies to bind as a stacked dimer to its cognate sequences: 5'-TACGAT-3' (5'-flank of the inverted CCAAT box-2 of the human topoisomerase IIalpha promoter) and 5'-ATCGAT-3'. A gel shift experiment provided evidence for f-PIP to inhibit protein-DNA interaction at the ICB2 site. Western blot studies showed that expression of the topoisomerase IIalpha gene in confluent NIH 3T3 cells was induced by treatment with f-PIP. The results suggested that the triamide was able to enter the nucleus, interacted with the target site within ICB2, inhibited NF-Y binding, and activated gene expression.     

Characterization of major recognition sequences for a herpes simplex virus type 1 origin-binding protein.

To investigate early initiation events in the replication of herpes simplex virus type 1, we analyzed interactions of proteins from infected cell extracts with the small origin of herpes simplex virus type 1 (oris1). Using the mobility shift assay, we detected two origin-specific binding interactions. We characterized the more prominent interaction on both strands of the DNA duplex with DNase I protection and methylation interference assays. Protein binding protects 17 bases of DNA on each strand from DNase I. These sequences are located at the left end of the central palindrome and are shifted four bases relative to one another. On the basis of the DNase protection pattern, we believe this protein to be related to the origin-binding protein defined by Elias et al. (P. Elias, M.E. O'Donnell, E.S. Mocarski, and I.R. Lehman, Proc. Natl. Acad. Sci. 83:6322-6326, 1986). Our DNase I footprint shows both strong and weak areas of protection. The regions strongly protected from DNase I align with the essential contact residues identified by interference footprinting. Methylation interference defines a small binding domain of 8 base pairs: 5'-GTTCGCAC-3'/3'-CAAGCGTG-5'. This recognition sequence contains two inverted 5'-GT(T/G)CG-3' repeats which share a 2-base overlap; thus, the origin-binding protein probably binds to the inverted repeats as a dimer.     

Modulation of the chemical and biological properties of trans platinum complexes: monofunctional platinum complexes containing one nucleobase as potential antiviral chemotypes

Replacement of one of the chloride leaving groups in trans-[PtCl2(NH3)(L)] by the nucleobase 9-ethylguanine gives the nucleobase cations [SP-4-2]-[PtCl(9-ethylguanine)(NH3)(L)]+ (L = NH3, 1; L = quinoline, 3), which are models for the monofunctional adduct on DNA. Displacement of Cl- in 1 and 3 by either 5'-guanosine monophosphate (5'-GMP) or N-acetyl-L-methionine (N-AcMet) showed clear kinetic preference for the sulfur (estimated half-lives of 1.5 and 4 h with N-AcMet against 7 and 17 h for 5'-GMP for 1 and 3, respectively). To further examine the kinetic preference, 1-methylcytosine (1-MeCyt) analogs were prepared, [SP-4-2]-[PtCl(1-Me-Cyt)(NH3)(L)]+ (L=NH3, 2; L=quinoline, 4). The -MeCyt compounds, 2 and 4, resulted in slower rates of substitution by both 5'-GMP and N-AcMet in comparison to 1 and 3 (estimated half-lives for N-AcMet of 5 and 13.5 h and for 5'-GMP of 6 and 14 h for 2 and 4, respectively). Interestingly in this case, however, no selectivity for the sulfur site was observed, a possible explanation being that molecular recognition across the square plane enhances the rate of reaction with 5'-GMP. The affinity of 3 towards S-donor ligands was exploited to remove zinc from the zinc-finger site of the C-terminal finger of the HIV-nucleocapsid protein, NCp7. The ability to eject zinc further suggested the biological antiviral application of [SP-4-2]-[PtCl(nucleobase)(NH3)(L)]+. A preliminary survey against HIV and herpes viruses indeed showed encouraging results with some antiviral specificity, dependent on the exact nature of the compound. The initial results suggest consideration of [SP-4-2]-[PtCl(nucleobase)(NH3)(L)]+ as a novel antiviral chemotype.     

Covalent carcinogenic lesions in DNA: NMR studies of O6-methylguanosine containing oligonucleotide duplexes

We report on proton and phosphorus high resolution NMR investigations of the self-complementary dodecanucleotide d(C1-G2-N3-G4-A5-A6-T7-T8-C9-O6meG10-C11-G12) duplexes (henceforth called O6 meG.N 12-mers), N = C, T, A and G, which contain N3.O6meG10 interactions in the interior of the helix. These sequences containing a single modified O6meG per strand were prepared by phosphoamidite synthesis and provide an excellent model for probing the structural basis for covalent carcinogenic lesions in DNA. Distance dependent nuclear Overhauser effect (NOE) measurements and line widths of imino protons demonstrate that the N3 and O6meG.10 bases stack into the duplex and are flanked by stable Watson-Crick base pairs at low temperature for all four O6meG.N 12-mer duplexes. The imino proton of T3 in the O6meG.T 12-mer and G3 in the O6meG.N 12-mer helix, which are associated with the modification site, resonate at unusually high field (8.5 to 9.0 ppm) compared to imino protons in Watson-Crick base pairs (12.5 to 14.5 ppm). The nonexchangeable base and sugar protons have been assigned from two dimensional correlated (COSY) and nuclear Overhauser effect (NOESY) measurements on the O6meG.N 12-mer helices. The directionality of the distance dependent NOEs establish all O6meG.N duplexes to be right-handed helices in solution. The glycosidic torsion angles are in the anti range at the N3.O6meG10 modification site except for O6meG10 in the O6meG.G 12-mer duplex which adopts a syn configuration. This results in altered NOEs between the G3 (anti).O6meG10 (syn) pair and flanking G2.C11 and G4.C9 base pairs in the O6meG.G 12-mer duplex. We observe pattern reversal for cross peaks in the COSY spectrum linking the sugar H1' protons with the H2',2" protons at the G2 and O6meG10 residues in the O6meG.N 12-mer duplexes with the effect least pronounced for the O6meG.T 12-mer helix. The proton chemical shift and NOE data have been analyzed to identify regions of conformational perturbations associated with N3.O6meG10 modification sites in the O6meG.N 12-mer duplexes. The proton decoupled phosphorus spectrum of O6meG.T 12-mer duplex exhibits an unperturbed phosphodiester backbone in contrast to the phosphorus spectra of the O6meG.C 12-mer, O6meG.G 12-mer and O6meG.A 12-mer duplexes which exhibit phosphorus resonances dispersed over 2 ppm characteristic of altered phosphodiester backbones at the modification site. Tentative proposals are put forward for N3.O6meG10 pairing models based on the available NMR data and serve as a guide for the design of future experiments.