DNA bending by charged peptides: electrophoretic and spectroscopic analyses

We are testing the idea that placement of fixed charges near one face of the DNA double helix can induce DNA bending by a purely electrostatic mechanism. If stretching forces between DNA phosphates are significant, fixed charges should induce DNA bending by asymmetrically modulating these forces. We have previously tested this hypothesis by adding charged residues to small bZIP DNA binding peptides and monitoring DNA bending using electrophoretic phasing assays. Our results were consistent with an electrostatic model of DNA bending in predicted directions. We now confirm these observations with fluorescence resonance energy transfer (FRET). Using a "U"-shaped DNA probe, we report that DNA bending by charged bZIP peptides is readily detected by FRET. We further show that charged bZIP peptides cause DNA bending rather than DNA twisting.     

DNA bending by bHLH charge variants

We wish to understand the role of electrostatics in DNA stiffness and bending. The DNA charge collapse model suggests that mutual electrostatic repulsions between neighboring phosphates significantly contribute to DNA stiffness. According to this model, placement of fixed charges near the negatively charged DNA surface should induce bending through asymmetric reduction or enhancement of these inter-phosphate repulsive forces. We have reported previously that charged variants of the elongated basic-leucine zipper (bZIP) domain of Gcn4p bend DNA in a manner consistent with this charge collapse model. To extend this result to a more globular protein, we present an investigation of the dimeric basic-helix–loop–helix (bHLH) domain of Pho4p. The 62 amino acid bHLH domain has been modified to position charged amino acid residues near one face of the DNA double helix. As observed for bZIP charge variants, DNA bending toward appended cations (away from the protein:DNA interface) is observed. However, unlike bZIP proteins, DNA is not bent away from bHLH anionic charges. This finding can be explained by the structure of the more globular bHLH domain which, in contrast to bZIP proteins, makes extensive DNA contacts along the binding face.     

An alpha-helical peptide model for electrostatic interactions of proteins with DNA. The N terminus of RecA

A series of synthetic peptides have been studied as models for non-specific protein-DNA interactions. In an alpha-helical conformation, the charged amino acid residues of the N-terminal 24 residues of RecA protein are asymmetrically distributed; at neutral pH there is a +4 charge on one face of the helix and a -3 charge on the other face. Modeling suggests that the positive face of the helix can bind five DNA phosphate groups by electrostatic interactions. Circular dichroism (c.d.) spectra indicate that the analogous peptide, Rec24 (AIDENKQKALAAALGQIEKQFGKG-amide), is largely unstructured in water but becomes highly helical in the presence of DNA. Peptide titrations of fluorescent etheno-DNA confirm that the changes in the c.d. spectrum of the peptide are associated with binding, although a dependence of the c.d. signal on the degree of DNA saturation is observed, indicating that peptide can be bound in more than one conformation. At saturation the peptide binds to 5.0(+/- 0.5) DNA phosphate groups as predicted and the electrostatic nature of the binding is confirmed by a strong dependence on salt concentration. A "mutant" peptide where an acidic glutamate residue replaces an alanine on the basic face of the Rec24 helix exhibits weaker binding to single-stranded DNA, also consistent with the electrostatic nature of the proposed peptide-DNA interaction. Extending Rec24 by ten amino acid residues, where the additional residues do not participate in the helical motif, does not noticeably affect binding. Thus, we show experimentally that an asymmetric charge distribution on an alpha-helix can represent an important element for binding nucleic acids.     

Selective DNA bending by a variety of bZIP proteins.

We have investigated DNA bending by bZIP family proteins that can bind to the AP-1 site. DNA bending is widespread, although not universal, among members of this family. Different bZIP protein dimers induced distinct DNA bends. The DNA bend angles ranged from virtually 0 to greater than 40 degrees as measured by phasing analysis and were oriented toward both the major and the minor grooves at the center of the AP-1 site. The DNA bends induced by the various heterodimeric complexes suggested that each component of the complex induced an independent DNA bend as previously shown for Fos and Jun. The Fos-related proteins Fra1 and Fra2 bent DNA in the same orientation as Fos but induced smaller DNA bend angles. ATF2 also bent DNA toward the minor groove in heterodimers formed with Fos, Fra2, and Jun. CREB and ATF1, which favor binding to the CRE site, did not induce significant DNA bending. Zta, which is a divergent member of the bZIP family, bent DNA toward the major groove. A variety of DNA structures can therefore be induced at the AP-1 site through combinatorial interactions between different bZIP family proteins. This diversity of DNA structures may contribute to regulatory specificity among the plethora of proteins that can bind to the AP-1 site.     

Peptide nucleic acid (PNA) amphiphiles: synthesis, self-assembly, and duplex stability

Peptide amphiphiles comprising a class of conjugates of peptide nucleic acid (PNA), natural amino acids, and n-alkanes were synthesized and studied. These PNA amphiphiles (PNAA) self-assemble at concentrations between 10 and 50 muM and exhibit water solubilities above 500 muM. The highly specific, stable DNA binding properties of PNAs are preserved by these modifications, with no significant differences between the thermodynamics of DNA binding of the PNA peptide and the PNA amphiphile. Proper solubilization of the PNAA required the attachment of (Lys)(2) and (Glu)(4) peptides to PNAs, which affected the PNAA-DNA duplex stability by electrostatic interactions between these charged amino acids and the negatively charged DNA backbone. These electrostatic effects did not affect the specificity of DNA binding, however. Electrostatic effects are screened with added salt, in a manner consistent with previous studies of PNA-DNA duplex stability and predictions from a charged-cylinder model for the duplex.     

Electrostatic free energy landscapes for DNA helix bending

Nucleic acids are highly charged polyanionic molecules; thus, the ionic conditions are crucial for nucleic acid structural changes such as bending. We use the tightly bound ion theory, which explicitly accounts for the correlation and ensemble effects for counterions, to calculate the electrostatic free energy landscapes for DNA helix bending. The electrostatic free energy landscapes show that DNA bending energy is strongly dependent on ion concentration, valency, and size. In a Na⁺ solution, DNA bending is electrostatically unfavorable because of the strong charge repulsion on backbone. With the increase of the Na⁺ concentration, the electrostatic bending repulsion is reduced and thus the bending becomes less unfavorable. In contrast, in an Mg²⁺ solution, ion correlation induces a possible attractive force between the different parts of the helical strands, resulting in bending. The electrostatically most favorable and unfavorable bending directions are toward the major and minor grooves, respectively. Decreasing the size of the divalent ions enhances the electrostatic bending attraction, causing an increased bending angle, and shifts the most favorable bending to the direction toward the minor groove. The microscopic analysis on ion-binding distribution reveals that the divalent ion-induced helix bending attraction may come from the correlated distribution of the ions across the grooves in the bending direction.     

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.     

Fluorescence resonance energy transfer between donor-acceptor pair on two oligonucleotides hybridized adjacently to DNA template

We use fluorescein as the energy donor and rhodamine as the acceptor to measure the efficiency of fluorescence resonance energy transfer (FRET) in a set of hybridized DNA constructs. The two fluorophores are covalently attached via linkers to two separate oligonucleotides with fluorescein at the 3' end of one oligonucleotide and rhodamine at the 5' end or in the middle of another nucleotide. For the FRET analysis both fluorophore-labeled oligonucleotides are hybridized to adjacent sections of the same DNA template to form a three-component duplex with a one base gap between the two labeled oligonucleotides. A similar configuration is implemented for a quantitative real-time polymerase chain reaction (PCR) with LightCycler technology, where a 1-5 base separation between donor and acceptor is recommended to optimize energy transfer efficiencies. Our constructs cover donor-acceptor separations from 2 to 17 base pairs (approximately 10-70 A). The results show that, when the two fluorophores are located at close distances (less than 8 base separation), FRET efficiencies are above 80%, although there may be ground-state interactions between fluorophores when the separation is under about 6 bases. Modeling calculations are used to predict the structure of these three-component constructs. The duplex mostly retains a normal double helical structure, although slight bending may occur near the unpaired base in the DNA template. Stable and reproducible energy transfer is also observed over the distance range investigated here in real-time thermal cycling. The study identifies important parameters that determine FRET response in applications such as real-time PCR.     

Charge neutralization and DNA bending by the Escherichia coli catabolite activator protein

We are interested in the role of asymmetric phosphate neutralization in DNA bending induced by proteins. We describe an experimental estimate of the actual electrostatic contribution of asymmetric phosphate neutralization to the bending of DNA by the Escherichia coli catabolite activator protein (CAP), a prototypical DNA-bending protein. Following assignment of putative electrostatic interactions between CAP and DNA phosphates based on X-ray crystal structures, appropriate phosphates in the CAP half-site DNA were chemically neutralized by methylphosphonate substitution. DNA shape was then evaluated using a semi-synthetic DNA electrophoretic phasing assay. Our results confirm that the unmodified CAP DNA half-site sequence is intrinsically curved by 26° in the direction enhanced in the complex with protein. In the absence of protein, neutralization of five appropriate phosphates increases DNA curvature to 32° (~23% increase), in the predicted direction. Shifting the placement of the neutralized phosphates changes the DNA shape, suggesting that sequence-directed DNA curvature can be modified by the asymmetry of phosphate neutralization. We suggest that asymmetric phosphate neutralization contributes favorably to DNA bending by CAP, but cannot account for the full DNA deformation.     

Mismatch-induced DNA unbending upon duplex opening

A DNA duplex can be torn open at a specific position by introducing a branch or bulge to create an asymmetric three-way junction (TWJ). The opened duplex manifests a bent conformation (bending angle approximately 60 degrees , relative to the unopened form), which leads to a dramatic decrease in gel electrophoretic mobility. In the presence of a basepair mismatch at the opening position, the DNA backbone becomes less bent and assumes a distorted T-shaped structure, resulting in an increase in polyacrylamide gel electrophoretic mobility. Both conformational changes are confirmed using fluorescence resonance energy transfer experiments and found to be similar to the signature conformational changes of DNA duplex upon MutS protein binding. Our results imply that some structural rearrangements essential for mismatch recognition are achievable without protein interference. The gel electrophoretic mobility data for DNA TWJs with and without base mismatches correlates well with rotational diffusivity, computed by taking into account the conformational change of TWJ induced by base mismatch.     

Design and calibration of a semi-synthetic DNA phasing assay

Electrophoretic assays of intrinsic DNA shape and shape changes induced by ligand binding are extremely useful because of their convenience and simplicity. The development of calibrations and empirical quantitative relationships permits highly accurate measurement of DNA shape using electrophoresis. Many conventional analyses employ the unidirectional ligation of short DNA duplexes. However, many oligonucleotides (typically more than 20) must often be synthesized for a single experiment. Additionally, the length of the DNA duplex can become limiting, preventing the analysis of certain DNA sequences. We now describe a semi-synthetic electrophoretic phasing method that offers several advantages, including a reduced number of required synthetic oligonucleotides, the ability to analyze longer DNA duplexes and a simplified approach for data analysis. We characterize semi-synthetic DNA probes in electrophoretic phasing assays by ligation of synthetic duplexes containing A₅ tracts between two longer restriction fragments. Upon electrophoresis, the gel mobility is strongly correlated with the predicted DNA curvature provided by the reference A₅ tracts. Having obtained this calibration, we show that the semi-synthetic phasing assay can be readily and economically applied to analyze DNA curvature induced by DNA charge modifications and DNA bending due to peptide binding.     

The contribution of the methyl groups on thymine bases to binding specificity and affinity by alanine-rich mutants of the bZIP motif.

We have used fluorescence anisotropy to measure in situ the thermodynamics of binding of alanine-rich mutants of the GCN4 basic region/leucine zipper (bZIP) to short DNA duplexes, in which thymines were replaced with uracils, in order to quantify the contributions of the C5 methyl group on thymines with alanine methyl side chains. We simplified the alpha-helical GCN4 bZIP by alanine substitution: 4A, 11A, and 18A contain four, 11, and 18 alanine mutations in their DNA-binding basic regions, respectively. Titration of fluorescein-labeled duplexes with increasing amounts of protein yielded dissociation constants in the low-to-mid nanomolar range for all bZIP mutants in complex with the AP-1 target site (5'-TGACTCA-3'); binding to the nonspecific control duplex was >1000-fold weaker. Small changes of <1 kcal/mol in binding free energies were observed for wild-type bZIP and 4A mutant to uracil-containing AP-1, whereas 11A and 18A bound almost equally well to native AP-1 and uracil-containing AP-1. These modest changes in binding affinities may reflect the multivalent nature of protein-DNA interactions, as our highly mutated proteins still exhibit native-like behavior. These protein mutations may compensate for changes in enthalpic and entropic contributions toward DNA-binding in order to maintain binding free energies similar to that of the native protein-DNA complex.     

DNA and Protein Requirements for Substrate Conformational Changes Necessary for Human Flap Endonuclease-1-catalyzed Reaction

Human flap endonuclease-1 (hFEN1) catalyzes the essential removal of single-stranded flaps arising at DNA junctions during replication and repair processes. hFEN1 biological function must be precisely controlled, and consequently, the protein relies on a combination of protein and substrate conformational changes as a prerequisite for reaction. These include substrate bending at the duplex-duplex junction and transfer of unpaired reacting duplex end into the active site. When present, 5′-flaps are thought to thread under the helical cap, limiting reaction to flaps with free 5′-termini in vivo. Here we monitored DNA bending by FRET and DNA unpairing using 2-aminopurine exciton pair CD to determine the DNA and protein requirements for these substrate conformational changes. Binding of DNA to hFEN1 in a bent conformation occurred independently of 5′-flap accommodation and did not require active site metal ions or the presence of conserved active site residues. More stringent requirements exist for transfer of the substrate to the active site. Placement of the scissile phosphate diester in the active site required the presence of divalent metal ions, a free 5′-flap (if present), a Watson-Crick base pair at the terminus of the reacting duplex, and the intact secondary structure of the enzyme helical cap. Optimal positioning of the scissile phosphate additionally required active site conserved residues Tyr⁴⁰, Asp¹⁸¹, and Arg¹⁰⁰ and a reacting duplex 5′-phosphate. These studies suggest a FEN1 reaction mechanism where junctions are bound and 5′-flaps are threaded (when present), and finally the substrate is transferred onto active site metals initiating cleavage.     

The GCN4 bZIP targets noncognate gene regulatory sequences: quantitative investigation of binding at full and half sites

We previously reported that a basic region/leucine zipper (bZIP) protein, a hybrid of the GCN4 basic region and C/EBP leucine zipper, not only recognizes cognate target sites AP-1 (5'-TGACTCA-3') and cAMP-response element (CRE) (5'-TGACGTCA-3') but also binds selectively to noncognate DNA sites: C/EBP (CCAAT/enhancer binding protein, 5'-TTGCGCAA), XRE1 (xenobiotic response element, 5'-TTGCGTGA), HRE (HIF response element, 5'-GCACGTAG), and E-box (5'-CACGTG). In this work, we used electrophoretic mobility shift assay (EMSA) and circular dichroism (CD) for more extensive characterization of the binding of wt bZIP dimer to noncognate sites as well as full- and half-site derivatives, and we examined changes in flanking sequences. Quantitative EMSA titrations were used to measure dissociation constants of this hybrid, wt bZIP, to DNA duplexes: Full-site binding affinities gradually decrease from cognate sites AP-1 and CRE with Kd values of 13 and 12 nM, respectively, to noncognate sites with Kd values of 120 nM to low microM. DNA-binding selectivity at half sites is maintained; however, half-site binding affinities sharply decrease from the cognate half site (Kd = 84 nM) to noncognate half sites (all Kd values > 2 microM). CD shows that comparable levels of alpha-helical structure are induced in wt bZIP upon binding to cognate AP-1 or noncognate sites. Thus, noncognate sites may contribute to preorganization of stable protein structure before binding target DNA sites. This work demonstrates that the bZIP scaffold may be a powerful tool in the design of small, alpha-helical proteins with desired DNA recognition properties.     

Location of cyanine-3 on double-stranded DNA: importance for fluorescence resonance energy transfer studies

Fluorescence resonance energy transfer provides valuable long-range distance information about macromolecules in solution. Fluorescein and Cy3 are an important donor-acceptor pair of fluorophores; the characteristic F?rster length for this pair on DNA is 56 A, so the pair can be used to study relatively long distances. Measurement of FRET efficiency for a series of DNA duplexes terminally labeled with fluorescein and Cy3 suggests that the Cy3 is close to the helical axis of the DNA. An NMR analysis of a self-complementary DNA duplex 5'-labeled with Cy3 shows that the fluorophore is stacked onto the end of the helix, in a manner similar to that of an additional base pair. This provides a known point from which distances calculated from FRET measurements are measured. Using the FRET efficiencies for the series of DNA duplexes as restraints, we have determined an effective position for the fluorescein, which is maximally extended laterally from the helix. The knowledge of the fluorophore positions can now be used for more precise interpretation of FRET data from nucleic acids.