Functional characterization of the Escherichia coli Fis-DNA binding sequence

The Escherichia coli protein Fis is remarkable for its ability to interact specifically with DNA sites of highly variable sequences. The mechanism of this sequence-flexible DNA recognition is not well understood. In a previous study, we examined the contributions of Fis residues to high-affinity binding at different DNA sequences using alanine-scanning mutagenesis and identified several key residues for Fis-DNA recognition. In this work, we investigated the contributions of the 15-bp core Fis binding sequence and its flanking regions to Fis-DNA interactions. Systematic base-pair replacements made in both half sites of a palindromic Fis binding sequence were examined for their effects on the relative Fis binding affinity. Missing contact assays were also used to examine the effects of base removal within the core binding site and its flanking regions on the Fis-DNA binding affinity. The results revealed that: (1) the − 7G and + 3Y bases in both DNA strands (relative to the central position of the core binding site) are major determinants for high-affinity binding; (2) the C⁵ methyl group of thymine, when present at the + 4 position, strongly hinders Fis binding; and (3) AT-rich sequences in the central and flanking DNA regions facilitate Fis-DNA interactions by altering the DNA structure and by increasing the local DNA flexibility. We infer that the degeneracy of specific Fis binding sites results from the numerous base-pair combinations that are possible at noncritical DNA positions (from − 6 to − 4, from − 2 to + 2, and from + 4 to + 6), with only moderate penalties on the binding affinity, the roughly similar contributions of − 3A or G and + 3T or C to the binding affinity, and the minimal requirement of three of the four critical base pairs to achieve considerably high binding affinities.     

The energetic contribution of induced electrostatic asymmetry to DNA bending by a site-specific protein

DNA bending can be promoted by reducing the net negative electrostatic potential around phosphates on one face of the DNA, such that electrostatic repulsion among phosphates on the opposite face drives bending toward the less negative surface. To provide the first assessment of energetic contribution to DNA bending when electrostatic asymmetry is induced by a site-specific DNA binding protein, we manipulated the electrostatics in the EcoRV endonuclease-DNA complex by mutation of cationic side chains that contact DNA phosphates and/or by replacement of a selected phosphate in each strand with uncharged methylphosphonate. Reducing the net negative charge at two symmetrically located phosphates on the concave DNA face contributes − 2.3 kcal mol^− 1 to − 0.9 kcal mol^− 1 (depending on position) to complex formation. In contrast, reducing negative charge on the opposing convex face produces a penalty of + 1.3 kcal mol^− 1. Förster resonance energy transfer experiments show that the extent of axial DNA bending (about 50°) is little affected in modified complexes, implying that modification affects the energetic cost but not the extent of DNA bending. Kinetic studies show that the favorable effects of induced electrostatic asymmetry on equilibrium binding derive primarily from a reduced rate of complex dissociation, suggesting stabilization of the specific complex between protein and markedly bent DNA. A smaller increase in the association rate may suggest that the DNA in the initial encounter complex is mildly bent. The data imply that protein-induced electrostatic asymmetry makes a significant contribution to DNA bending but is not itself sufficient to drive full bending in the specific EcoRV-DNA complex.     

Identification of Critical Residues for the Tight Binding of Both Correct and Incorrect Nucleotides to Human DNA Polymerase λ

DNA polymerase λ (Pol λ) is a novel X-family DNA polymerase that shares 34% sequence identity with DNA polymerase β. Pre-steady-state kinetic studies have shown that the Pol λ-DNA complex binds both correct and incorrect nucleotides 130-fold tighter, on average, than the DNA polymerase β-DNA complex, although the base substitution fidelity of both polymerases is 10^− 4 to 10^− 5. To better understand Pol λ's tight nucleotide binding affinity, we created single-substitution and double-substitution mutants of Pol λ to disrupt the interactions between active-site residues and an incoming nucleotide or a template base. Single-turnover kinetic assays showed that Pol λ binds to an incoming nucleotide via cooperative interactions with active-site residues (R386, R420, K422, Y505, F506, A510, and R514). Disrupting protein interactions with an incoming correct or incorrect nucleotide impacted binding to each of the common structural moieties in the following order: triphosphate ? base > ribose. In addition, the loss of Watson-Crick hydrogen bonding between the nucleotide and the template base led to a moderate increase in K_d. The fidelity of Pol λ was maintained predominantly by a single residue, R517, which has minor groove interactions with the DNA template.     

The stability of DNA-porphyrin complexes in the presence of Mn(II) ions

The complex formation of porphyrins with DNA leads to changes of stability of DNA. In the present study we investigated binding properties and the thermodynamic parameters of a water-soluble, cationic planar Cu(II)-containing meso-tetrakis(4-N-butyl-pyridiniumyl)porphyrin [CuTButPyP4] and nonplanar Co(II)-containing meso-tetrakis(4-N-butyl-pyridiniumyl)porphyrin [CoButPyP4] with calf thymus DNA in the presence of divalent manganese ions. For displaying the changes of thermodynamic parameters (T_m and ΔT) the melting curves of DNA-porphyrin complexes in the presence of Mn²⁺ ions have been obtained. The enthalpy (ΔH) of helix-coil transition has been also evaluated. It was shown that the binding of ions to DNA proceeds in two stages depending on the manganese/DNA phosphates molar ratio [Mn]/[P]. At the first stage (0.001 < [Mn]/[P] < 1), the interaction of manganese ions with DNA phosphates occurs, causing an additional screening of their negative charge and the stabilization of the double helix. As a result, the best conditions for intercalation of CuTButPyP4 or of peripheral rings of CoButPyP4 occur. The significant increase of T_m, but less changes of ΔT were observed. At the second stage (1 < [Mn]/[P] < 4), the ions interact with both the phosphates and the nitrogen bases of DNA. At this stage, it is possible for the manganese ion to coordinate simultaneously to the oxygen atom of the phosphate and the neighboring base of DNA. At a higher [Mn]/[P] ratio, the destabilization of the double helix begins, and partial breakage of the hydrogen bonds between the nitrogen bases occurs. Respectively the destabilization of DNA in the presence of both porphyrins takes place.     

Linkage of cation binding and folding in human telomeric quadruplex DNA

Formation of DNA quadruplexes requires monovalent cation binding. To characterize the cation binding stoichiometry and linkage between binding and folding, we carried out KCl titrations of Tel22 (d[A(GGGTTA)₃]), a model of the human telomere sequence, using a fluorescent indicator to determine [K⁺]_free and circular dichroism to assess the extent of folding. At [K⁺]_free = 5 mM (sufficient for > 95% folding), the apparent binding stoichiometry is 3K⁺/Tel22; at [K⁺]_free = 20 mM, it increased to 8-10K⁺/Tel22. Thermodynamic analysis shows that at [K⁺]_free = 5 mM, K⁺ binding contributes approximately − 4.9 kcal/mol for folding Tel22. The overall folding free energy is − 2.4 kcal/mol, indicating that there are energetically unfavorable contributions to folding. Thus, quadruplex folding is driven almost entirely by the energy of cation binding with little or no contribution from other weak molecular interactions.     

Diffusion of Human Replication Protein A along Single-Stranded DNA

Replication protein A (RPA) is a eukaryotic single-stranded DNA (ssDNA) binding protein that plays critical roles in most aspects of genome maintenance, including replication, recombination and repair. RPA binds ssDNA with high affinity, destabilizes DNA secondary structure and facilitates binding of other proteins to ssDNA. However, RPA must be removed from or redistributed along ssDNA during these processes. To probe the dynamics of RPA–DNA interactions, we combined ensemble and single-molecule fluorescence approaches to examine human RPA (hRPA) diffusion along ssDNA and find that an hRPA heterotrimer can diffuse rapidly along ssDNA. Diffusion of hRPA is functional in that it provides the mechanism by which hRPA can transiently disrupt DNA hairpins by diffusing in from ssDNA regions adjacent to the DNA hairpin. hRPA diffusion was also monitored by the fluctuations in fluorescence intensity of a Cy3 fluorophore attached to the end of ssDNA. Using a novel method to calibrate the Cy3 fluorescence intensity as a function of hRPA position on the ssDNA, we estimate a one-dimensional diffusion coefficient of hRPA on ssDNA of D₁ ~ 5000 nt² s^− 1 at 37 °C. Diffusion of hRPA while bound to ssDNA enables it to be readily repositioned to allow other proteins access to ssDNA.     

Structure and dynamics of a DNA-based model system for the study of electron spin-spin interactions

We report on the structure and dynamics of a model system for measuring long-range distances in biological macromolecules by saturation-recovery EPR. Four DNA duplexes that incorporate a paramagnetic dysprosium ion (Dy(III)) and a nitroxide spin-label were examined by electron paramagnetic resonance (EPR), circular dichroism (CD), and ultra-violet absorbance (UV) spectroscopy. Dy(III) is chelated by the modified base deoxythymidine-EDTA, (dT-EDTA). Electron spin-spin interactions between the Dy(III) ion and the nitroxide radical are observed at distances as great as ∼5.3 nm. A slight change in the conformation of those nucleotides lying between the EDTA(Dy(III)) complex and the nitroxide spin-label results in a “stiffening” of the DNA helix on the EPR time scale. Changes in conformation and helix dynamics are due to the binding of the EDTA(Dy(III)) complex to the phosphodiester backbone of the complementary strand. Molecular mechanics calculations indicate that binding occurs in the 5′ direction on the complementary strand, at a position 3 or 4 phosphates distant from the dT-EDTA(Dy(III)) * dA base pair.     

The Macroscopic Rate of Nucleic Acid Translocation by Hepatitis C Virus Helicase NS3h Is Dependent on Both Sugar and Base Moieties

The nonstructural protein 3 helicase (NS3h) of hepatitis C virus is a 3′-to-5′ superfamily 2 RNA and DNA helicase that is essential for the replication of hepatitis C virus. We have examined the kinetic mechanism of the translocation of NS3h along single-stranded nucleic acid with bases uridylate (rU), deoxyuridylate (dU), and deoxythymidylate (dT), and have found that the macroscopic rate of translocation is dependent on both the base moiety and the sugar moiety of the nucleic acid, with approximate macroscopic translocation rates of 3 nt s^− 1 (oligo(dT)), 35 nt s^− 1 (oligo(dU)), and 42 nt s^− 1 (oligo(rU)), respectively. We found a strong correlation between the macroscopic translocation rates and the binding affinity of the translocating NS3h protein for the respective substrates such that weaker affinity corresponded to faster translocation. The values of K_0.5 for NS3h translocation at a saturating ATP concentration are as follows: 3.3 ± 0.4 μM nucleotide (poly(dT)), 27 ± 2 μM nucleotide (poly(dU)), and 36 ± 2 μM nucleotide (poly(rU)). Furthermore, results of the isothermal titration of NS3h with these oligonucleotides suggest that differences in TΔS⁰ are the principal source of differences in the affinity of NS3h binding to these substrates. Interestingly, despite the differences in macroscopic translocation rates and binding affinities, the ATP coupling stoichiometries for NS3h translocation were identical for all three substrates (∼ 0.5 ATP molecule consumed per nucleotide translocated). This similar periodicity of ATP consumption implies a similar mechanism for NS3h translocation along RNA and DNA substrates.     

DNA binding mode transitions of Escherichia coli HU(alphabeta): evidence for formation of a bent DNA--protein complex on intact, linear duplex DNA

Escherichia coli HU_αβ, a major nucleoid-associated protein, organizes chromosomal DNA and facilitates numerous DNA transactions. Using isothermal titration calorimetry, fluorescence resonance energy transfer and a series of DNA lengths (8 bp, 15 bp, 34 bp, 38 bp and 160 bp) we established that HU_αβ interacts with duplex DNA using three different nonspecific binding modes. Both the HU to DNA molar ratio ([HU]/[DNA]) and DNA length dictate the dominant HU binding mode. On sufficiently long DNA (≥ 34 bp), at low [HU]/[DNA], HU populates a noncooperative 34 bp binding mode with a binding constant of 2.1 ± 0.4 × 10⁶ M^− 1, and a binding enthalpy of + 7.7 ± 0.6 kcal/mol at 15 °C and 0.15 M Na⁺. With increasing [HU]/[DNA], HU bound in the noncooperative 34 bp mode progressively converts to two cooperative (ω∼20) modes with site sizes of 10 bp and 6 bp. These latter modes exhibit smaller binding constants (1.1 ± 0.2 × 10⁵ M^− 1 for the 10 bp mode, 3.5 ± 1.4 × 10⁴ M^− 1 for the 6 bp mode) and binding enthalpies (4.2 ± 0.3 kcal/mol for the 10 bp mode, − 1.6 ± 0.3 kcal/mol for the 6 bp mode). As DNA length increases to 34 bp or more at low [HU]/[DNA], the small modes are replaced by the 34 bp binding mode. Fluorescence resonance energy transfer data demonstrate that the 34 bp mode bends DNA by 143 ± 6° whereas the 6 bp and 10 bp modes do not. The model proposed in this study provides a novel quantitative and comprehensive framework for reconciling previous structural and solution studies of HU, including single molecule (force extension measurement), fluorescence, and electrophoretic gel mobility-shift assays. In particular, it explains how HU condenses or extends DNA depending on the relative concentrations of HU and DNA.     

Berberine–DNA complexation: New insights into the cooperative binding and energetic aspects

The equilibrium binding of the cytotoxic plant alkaloid berberine to various DNAs and energetics of the interaction have been studied. At low ratios of bound alkaloid to base pair, the binding exhibited cooperativity to natural DNAs having almost equal proportions of AT and GC sequences. In contrast, the binding was non-cooperative to DNAs with predominantly high AT or GC sequences. Among the synthetic DNAs, cooperative binding was observed with poly(dA).poly(dT) and poly(dG).poly(dC) while non-cooperative binding was seen with poly(dA–dT).poly(dA–dT) and poly(dG–dC).poly(dG–dC). Both cooperative and non-cooperative bindings were remarkably dependent on the salt concentration of the media. Linear plots of ln K_a versus [Na⁺] for poly(dA).poly(dT) and poly(dA–dT).poly(dA–dT) showed the release of 0.56 and 0.75 sodium ions respectively per bound alkaloid. Isothermal titration calorimetry results revealed the binding to be exothermic and favoured by both enthalpy and entropy changes in all DNAs except the two AT polymers and AT rich DNA, where the same was predominantly entropy driven. Heat capacity values (ΔCp^o) of berberine binding to poly(dA).poly(dT), poly(dA–dT).poly(dA–dT), Clostridium perfringens and calf thymus DNA were − 98, − 140, − 120 and − 110 cal/mol K respectively. This study presents new insights into the binding dependent base pair heterogeneity in DNA conformation and the first complete thermodynamic profile of berberine binding to DNAs.     

Relative Domain Folding and Stability of a Membrane Transport Protein

There is a limited understanding of the folding of multidomain membrane proteins. Lactose permease (LacY) of Escherichia coli is an archetypal member of the major facilitator superfamily of membrane transport proteins, which contain two domains of six transmembrane helices each. We exploit chemical denaturation to determine the unfolding free energy of LacY and employ Trp residues as site-specific thermodynamic probes. Single Trp LacY mutants are created with the individual Trps situated at mirror image positions on the two LacY domains. The changes in Trp fluorescence induced by urea denaturation are used to construct denaturation curves from which unfolding free energies can be determined. The majority of the single Trp tracers report the same stability and an unfolding free energy of approximately + 2 kcal mol^− 1. There is one exception; the fluorescence of W33 at the cytoplasmic end of helix I on the N domain is unaffected by urea. In contrast, the equivalent position on the first helix, VII, of the C-terminal domain exhibits wild-type stability, with the single Trp tracer at position 243 on helix VII reporting an unfolding free energy of + 2 kcal mol^− 1. This indicates that the region of the N domain of LacY at position 33 on helix I has enhanced stability to urea, when compared the corresponding location at the start of the C domain. We also find evidence for a potential network of stabilising interactions across the domain interface, which reduces accessibility to the hydrophilic substrate binding pocket between the two domains.     

Kinetic control of Mg2+-dependent melting of duplex DNA ends by Escherichia coli RecBC

Escherichia coli RecBCD is a highly processive DNA helicase involved in double-strand break repair and recombination that possesses two helicase/translocase subunits with opposite translocation directionality (RecB (3′ to 5′) and RecD (5′ to 3′)). RecBCD has been shown to melt out ∼ 5-6 bp upon binding to a blunt-ended duplex DNA in a Mg²⁺-dependent, but ATP-independent reaction. Here, we examine the binding of E. coli RecBC helicase (minus RecD), also a processive helicase, to duplex DNA ends in the presence and in the absence of Mg²⁺ in order to determine if RecBC can also melt a duplex DNA end in the absence of ATP. Equilibrium binding of RecBC to DNA substrates with ends possessing pre-formed 3′ and/or 5′ single-stranded (ss)-(dT)_n flanking regions (tails) (n ranging from zero to 20 nt) was examined by competition with a fluorescently labeled reference DNA and by isothermal titration calorimetry. The presence of Mg²⁺ enhances the affinity of RecBC for DNA ends possessing 3′ or 5′-(dT)_n ssDNA tails with n < 6 nt, with the relative enhancement decreasing as n increases from zero to six nt. No effect of Mg²⁺ was observed for either the binding constant or the enthalpy of binding (ΔH_obs) for RecBC binding to DNA with ssDNA tail lengths, n ≥ 6 nucleotides. Upon RecBC binding to a blunt duplex DNA end in the presence of Mg²⁺, at least 4 bp at the duplex end become accessible to KMnO₄ attack, consistent with melting of the duplex end. Since Mg²⁺ has no effect on the affinity or binding enthalpy of RecBC for a DNA end that is fully pre-melted, this suggests that the role of Mg²⁺ is to overcome a kinetic barrier to melting of the DNA by RecBC and presumably also by RecBCD. These data also provide an accurate estimate (ΔH_obs = 8 ± 1 kcal/mol) for the average enthalpy change associated with the melting of a DNA base-pair by RecBC.     

The Osmolyte TMAO Stabilizes Native RNA Tertiary Structures in the Absence of Mg: Evidence for a Large Barrier to Folding from Phosphate Dehydration

The stabilization of RNA tertiary structures by ions is well known, but the neutral osmolyte trimethylamine oxide (TMAO) can also effectively stabilize RNA tertiary structure. To begin to understand the physical basis for the effects of TMAO on RNA, we have quantitated the TMAO-induced stabilization of five RNAs with known structures. So-called m values, the increment in unfolding free energy per molal of osmolyte at constant KCl activity, are ∼ 0 for a hairpin secondary structure and between 0.70 and 1.85 kcal mol^− 1m^− 1 for four RNA tertiary structures (30-86 nt). Further analysis of two RNAs by small-angle X-ray scattering and hydroxyl radical probing shows that TMAO reduces the radius of gyration of the unfolded ensemble to the same endpoint as seen in titration with Mg²⁺ and that the structures stabilized by TMAO and Mg²⁺ are indistinguishable. Remarkably, TMAO induces the native conformation of a Mg²⁺ ion chelation site formed in part by a buried phosphate, even though Mg²⁺ is absent. TMAO interacts weakly, if at all, with KCl, ruling out the possibility that TMAO stabilizes RNA indirectly by increasing salt activity. TMAO is, however, strongly excluded from the vicinity of dimethylphosphate (unfavorable interaction free energy, + 211 cal mol^− 1m^− 1 for the potassium salt), an ion that mimics the RNA backbone phosphate. We suggest that formation of RNA tertiary structure is accompanied by substantial phosphate dehydration (loss of 66-173 water molecules in the RNA structures studied) and that TMAO works principally by reducing the energetic penalty associated with this dehydration. The strong parallels we find between the effects of TMAO and Mg²⁺ suggest that RNA sequence is more important than specific ion interactions in specifying the native structure.     

HMGB binding to DNA: single and double box motifs

High mobility group (HMG) proteins are nuclear proteins believed to significantly affect DNA interactions by altering nucleic acid flexibility. Group B (HMGB) proteins contain HMG box domains known to bind to the DNA minor groove without sequence specificity, slightly intercalating base pairs and inducing a strong bend in the DNA helical axis. A dual-beam optical tweezers system is used to extend double-stranded DNA (dsDNA) in the absence as well as presence of a single box derivative of human HMGB2 [HMGB2(box A)] and a double box derivative of rat HMGB1 [HMGB1(box A+box B)]. The single box domain is observed to reduce the persistence length of the double helix, generating sharp DNA bends with an average bending angle of 99 ± 9° and, at very high concentrations, stabilizing dsDNA against denaturation. The double box protein contains two consecutive HMG box domains joined by a flexible tether. This protein also reduces the DNA persistence length, induces an average bending angle of 77 ± 7°, and stabilizes dsDNA at significantly lower concentrations. These results suggest that single and double box proteins increase DNA flexibility and stability, albeit both effects are achieved at much lower protein concentrations for the double box. In addition, at low concentrations, the single box protein can alter DNA flexibility without stabilizing dsDNA, whereas stabilization at higher concentrations is likely achieved through a cooperative binding mode.