Stepwise binding and bending of DNA by Escherichia coli integration host factor

The effect of DNA-binding proteins on insertion sequence element transposition upstream of the bgl operon in Escherichia coli

The bglGFB operon in Escherichia coli K-12 strain BW25113, encoding the proteins necessary for the uptake and metabolism of β-glucosides, is normally not expressed. Insertion of either IS1 or IS5 upstream of the bgl promoter activates expression of the operon only when the cell is starving in the presence of a β-glucoside, drastically increasing transcription and allowing the cell to survive and grow using this carbon source. Details surrounding the exact mechanism and regulation of the IS insertional event remain unclear. In this work, the role of several DNA-binding proteins in how they affect the rate of insertion upstream of bgl are examined via mutation assays and protocols measuring transcription. Both Crp and IHF exert a positive effect on insertional Bgl+ mutations when present, active, and functional in the cell. Our results characterize IHF's effect in conjunction with other mutations, show that IHF's effect on IS insertion into bgl also affects other operons, and indicate that it may exert its effect by binding to and altering the DNA conformation of IS1 and IS5 in their native locations, rather than by directly influencing transposase gene expression. In contrast, the cAMP-CRP complex acts directly upon the bgl operon by binding upstream of the promoter, presumably altering local DNA into a conformation that enhances IS insertion.

IHF and Fis as Escherichia coli Cell Cycle Regulators: Activation of the Replication Origin oriC and the Regulatory Cycle of the DnaA Initiator

This review summarizes current knowledge about the mechanisms of timely binding and dissociation of two nucleoid proteins, IHF and Fis, which play fundamental roles in the initiation of chromosomal DNA replication in Escherichia coli. Replication is initiated from a unique replication origin called oriC and is tightly regulated so that it occurs only once per cell cycle. The timing of replication initiation at oriC is rigidly controlled by the timely binding of the initiator protein DnaA and IHF to oriC. The first part of this review presents up-to-date knowledge about the timely stabilization of oriC-IHF binding at oriC during replication initiation. Recent advances in our understanding of the genome-wide profile of cell cycle-coordinated IHF binding have revealed the oriC-specific stabilization of IHF binding by ATP-DnaA oligomers at oriC and by an initiation-specific IHF binding consensus sequence at oriC. The second part of this review summarizes the mechanism of the timely regulation of DnaA activity via the chromosomal loci DARS2 (DnaA-reactivating sequence 2) and datA. The timing of replication initiation at oriC is controlled predominantly by the phosphorylated form of the adenosine nucleotide bound to DnaA, i.e., ATP-DnaA, but not ADP-ADP, is competent for initiation. Before initiation, DARS2 increases the level of ATP-DnaA by stimulating the exchange of ADP for ATP on DnaA. This DARS2 function is activated by the site-specific and timely binding of both IHF and Fis within DARS2. After initiation, another chromosomal locus, datA, which inactivates ATP-DnaA by stimulating ATP hydrolysis, is activated by the timely binding of IHF. A recent study has shown that ATP-DnaA oligomers formed at DARS2-Fis binding sites competitively dissociate Fis via negative feedback, whereas IHF regulation at DARS2 and datA still remains to be investigated. This review summarizes the current knowledge about the specific role of IHF and Fis in the regulation of replication initiation and proposes a mechanism for the regulation of timely IHF binding and dissociation at DARS2 and datA.

Integration Host Factor Binds DNA Holliday Junctions

Integration host factor (IHF) is a nucleoid-associated protein involved in DNA packaging, integration of viral DNA and recombination. IHF binds with nanomolar affinity to duplex DNA containing a 13 bp consensus sequence, inducing a bend of ~160° upon binding. We determined that IHF binds to DNA Four-way or Holliday junctions (HJ) with high affinity regardless of the presence of the consensus sequence, signifying a structure-based mechanism of recognition. Junctions, important intermediates in DNA repair and homologous recombination, are dynamic and can adopt either an open or stacked conformation, where the open conformation facilitates branch migration and strand exchange. Using ensemble and single molecule Förster resonance energy transfer (FRET) methods, we investigated IHF-induced changes in the population distribution of junction conformations and determined that IHF binding shifts the population to the open conformation. Further analysis of smFRET dynamics revealed that even in the presence of protein, the junctions remain dynamic as fast transitions are observed for the protein-bound open state. Protein binding alters junction conformational dynamics, as cross correlation analyses reveal the protein slows the transition rate at 1 mM Mg²⁺ but accelerates the transition rate at 10 mM Mg²⁺. Stopped flow kinetic experiments provide evidence for two binding steps, a rapid, initial binding step followed by a slower step potentially associated with a conformational change. These measurements also confirm that the protein remains bound to the junction during the conformer transitions and further suggest that the protein forms a partially dissociated state that allows junction arms to be dynamic. These findings, which demonstrate that IHF binds HJs with high affinity and stabilizes junctions in the open conformation, suggest that IHF may play multiple roles in the processes of integration and recombination in addition to stabilizing bacterial biofilms.

Reciprocally rewiring and repositioning the Integration Host Factor (IHF) subunit genes in Salmonella enterica serovar Typhimurium: impacts on physiology and virulence

The Integration Host Factor (IHF) is a heterodimeric nucleoid-associated protein that plays roles in bacterial nucleoid architecture and genome-wide gene regulation. The ihfA and ihfB genes encode the subunits and are located 350 kbp apart, in the Right replichore of the Salmonella chromosome. IHF is composed of one IhfA and one IhfB subunit. Despite this 1 : 1 stoichiometry, MS revealed that IhfB is produced in 2-fold excess over IhfA. We re-engineered Salmonella to exchange reciprocally the protein-coding regions of ihfA and ihfB, such that each relocated protein-encoding region was driven by the expression signals of the other's gene. MS showed that in this 'rewired' strain, IhfA is produced in excess over IhfB, correlating with enhanced stability of the hybrid ihfB-ihfA mRNA that was expressed from the ihfB promoter. Nevertheless, the rewired strain grew at a similar rate to the wild-type and was similar in competitive fitness. However, compared to the wild-type, it was less motile, had growth-phase-specific reductions in SPI-1 and SPI-2 gene expression, and was engulfed at a higher rate by RAW macrophage. Our data show that while exchanging the physical locations of its ihf genes and the rewiring of their regulatory circuitry are well tolerated in Salmonella, genes involved in the production of type 3 secretion systems exhibit dysregulation accompanied by altered phenotypes.

Whole-Genome Analysis Reveals That the Nucleoid Protein IHF Predominantly Binds to the Replication Origin oriC Specifically at the Time of Initiation

The structure and function of bacterial chromosomes are dynamically regulated by a wide variety of nucleoid-associated proteins (NAPs) and DNA superstructures, such as DNA supercoiling. In Escherichia coli, integration host factor (IHF), a NAP, binds to specific transcription promoters and regulatory DNA elements of DNA replication such as the replication origin oriC: binding to these elements depends on the cell cycle but underlying mechanisms are unknown. In this study, we combined GeF-seq (genome footprinting with high-throughput sequencing) with synchronization of the E. coli cell cycle to determine the genome-wide, cell cycle-dependent binding of IHF with base-pair resolution. The GeF-seq results in this study were qualified enough to analyze genomic IHF binding sites (e.g., oriC and the transcriptional promoters of ilvG and osmY) except some of the known sites. Unexpectedly, we found that before replication initiation, oriC was a predominant site for stable IHF binding, whereas all other loci exhibited reduced IHF binding. To reveal the specific mechanism of stable oriC–IHF binding, we inserted a truncated oriC sequence in the terC (replication terminus) locus of the genome. Before replication initiation, stable IHF binding was detected even at this additional oriC site, dependent on the specific DnaA-binding sequence DnaA box R1 within the site. DnaA oligomers formed on oriC might protect the oriC–IHF complex from IHF dissociation. After replication initiation, IHF rapidly dissociated from oriC, and IHF binding to other sites was sustained or stimulated. In addition, we identified a novel locus associated with cell cycle-dependent IHF binding. These findings provide mechanistic insight into IHF binding and dissociation in the genome.

Analysis of nucleoid-associated protein-binding regions reveals DNA structural features influencing genome organization in Mycobacterium tuberculosis

Nucleoid-Associated Proteins (NAPs) maintain bacterial nucleoid configuration through their architectural properties of DNA bending, wrapping, and bridging. However, the contribution of DNA structural alterations to DNA-NAP recognition at the genomic scale remains unresolved. Present work dissects the DNA sequence, shape and altered structural preferences at a genomic scale for six NAPs in Mycobacterium tuberculosis. Results suggest narrower minor groove width and higher DNA rigidity are marked for the binding sites of EspR and Lsr2, while mIHF, MtHU and NapM have heterogeneous DNA structural predilections. In contrast, WhiB4-DNA binding sites were characterized by wider minor groove width, highly deformable and less curved DNA. This work provides systematic insight into NAP-mediated genome organization as a function of DNA structural features.

The Site-Specific Recombination System of the Escherichia coli Bacteriophage Φ24B

Stx bacteriophages are members of the lambdoid group of phages and are responsible for Shiga toxin (Stx) production and the dissemination of Shiga toxin genes (stx) across shigatoxigenic Escherichia coli (STEC). These toxigenic bacteriophage hosts can cause life-threatening illnesses, and Stx is the virulence determinant responsible for the severe nature of infection with enterohemorrhagic E. coli, a subset of pathogenic STEC. Stx phages are temperate, and in the present study, the identification of what is actually required for Stx phage Φ24_B and bacterial DNA recombination was tested using both in vitro and in situ recombination assays. It is well established that phage λ, which underpins most of what we understand about lambdoid phage biology, requires its own encoded phage attachment site (attP) of 250 bp, a host-encoded attachment site (attB) of 21 bp, and a host-encoded DNA binding protein known as integration host factor (IHF). The assays applied in this study enabled the manipulation of the phage attachment site (attP) and the bacterial attachment site (attB) sequences and the inclusion or exclusion of a host-encoded accessory element known as integration host factor. We were able to demonstrate that the minimal attP sequence required by Φ24_B phage is between 350 and 427 bp. Unlike phage λ, the minimal necessary flanking sequences for the attB site do not appear to be equal in size, with a total length between 62 and 93 bp. Furthermore, we identified that the Φ24_B integrase does not require IHF to drive the integration and the recombination process. Understanding how this unusual Stx phage integrase works may enable exploitation of its promiscuous nature in the context of genetic engineering.

Evidence for a bind-then-bend mechanism for architectural DNA binding protein yNhp6A

The yeast Nhp6A protein (yNhp6A) is a member of the eukaryotic HMGB family of chromatin factors that enhance apparent DNA flexibility. yNhp6A binds DNA nonspecifically with nM affinity, sharply bending DNA by >60°. It is not known whether the protein binds to unbent DNA and then deforms it, or if bent DNA conformations are ‘captured’ by protein binding. The former mechanism would be supported by discovery of conditions where unbent DNA is bound by yNhp6A. Here, we employed an array of conformational probes (FRET, fluorescence anisotropy, and circular dichroism) to reveal solution conditions in which an 18-base-pair DNA oligomer indeed remains bound to yNhp6A while unbent. In 100 mM NaCl, yNhp6A-bound DNA unbends as the temperature is raised, with no significant dissociation of the complex detected up to ∼45°C. In 200 mM NaCl, DNA unbending in the intact yNhp6A complex is again detected up to ∼35°C. Microseconds-resolved laser temperature-jump perturbation of the yNhp6a–DNA complex revealed relaxation kinetics that yielded unimolecular DNA bending/unbending rates on timescales of 500 μs−1 ms. These data provide the first direct observation of bending/unbending dynamics of DNA in complex with yNhp6A, suggesting a bind-then-bend mechanism for this protein.

Two-step interrogation then recognition of DNA binding site by Integration Host Factor: an architectural DNA-bending protein

The dynamics and mechanism of how site-specific DNA-bending proteins initially interrogate potential binding sites prior to recognition have remained elusive for most systems. Here we present these dynamics for Integration Host factor (IHF), a nucleoid-associated architectural protein, using a μs-resolved T-jump approach. Our studies show two distinct DNA-bending steps during site recognition by IHF. While the faster (∼100 μs) step is unaffected by changes in DNA or protein sequence that alter affinity by >100-fold, the slower (1–10 ms) step is accelerated ∼5-fold when mismatches are introduced at DNA sites that are sharply kinked in the specific complex. The amplitudes of the fast phase increase when the specific complex is destabilized and decrease with increasing [salt], which increases specificity. Taken together, these results indicate that the fast phase is non-specific DNA bending while the slow phase, which responds only to changes in DNA flexibility at the kink sites, is specific DNA kinking during site recognition. Notably, the timescales for the fast phase overlap with one-dimensional diffusion times measured for several proteins on DNA, suggesting that these dynamics reflect partial DNA bending during interrogation of potential binding sites by IHF as it scans DNA.

Enhanced spontaneous DNA twisting/bending fluctuations unveiled by fluorescence lifetime distributions promote mismatch recognition by the Rad4 nucleotide excision repair complex

Rad4/XPC recognizes diverse DNA lesions including ultraviolet-photolesions and carcinogen-DNA adducts, initiating nucleotide excision repair. Studies have suggested that Rad4/XPC senses lesion-induced helix-destabilization to flip out nucleotides from damaged DNA sites. However, characterizing how DNA deformability and/or distortions impact recognition has been challenging. Here, using fluorescence lifetime measurements empowered by a maximum entropy algorithm, we mapped the conformational heterogeneities of artificially destabilized mismatched DNA substrates of varying Rad4-binding specificities. The conformational distributions, as probed by FRET between a cytosine-analog pair exquisitely sensitive to DNA twisting/bending, reveal a direct connection between intrinsic DNA deformability and Rad4 recognition. High-specificity CCC/CCC mismatch, free in solution, sampled a strikingly broad range of conformations from B-DNA-like to highly distorted conformations that resembled those observed with Rad4 bound; the extent of these distortions increased with bound Rad4 and with temperature. Conversely, the non-specific TAT/TAT mismatch had a homogeneous, B-DNA-like conformation. Molecular dynamics simulations also revealed a wide distribution of conformations for CCC/CCC, complementing experimental findings. We propose that intrinsic deformability promotes Rad4 damage recognition, perhaps by stalling a diffusing protein and/or facilitating ‘conformational capture’ of pre-distorted damaged sites. Surprisingly, even mismatched DNA specifically bound to Rad4 remains highly dynamic, a feature that may reflect the versatility of Rad4/XPC to recognize many structurally dissimilar lesions.

Direct observation of the external force mediated conformational dynamics of an IHF bound Holliday junction

We have investigated the isomerization dynamics and plausible energy landscape of 4-way Holliday junctions (4WHJs) bound to integration host factor (IHF, a DNA binding protein), considering the effect of applied external force, by single-molecule FRET methods. A slowing down of the forward as well as the backward rates of the isomerization process of the protein bound 4WHJ has been observed under the influence of an external force, which indicates an imposed restriction on the conformational switching. This has also been reflected by an increase in rigidity, as observed from the increase in the single-molecule FRET (smFRET)-anisotropy values (0.270 ± 0.012 to 0.360 ± 0.008). The application of an external force has assisted the conformational transitions to share the unstacked open structure intermediate, with different rate-limiting steps and a huge induced variation in the energy landscape. Furthermore, the associated landscape of the 4WHJ is visualized in terms of rarely interconverting states embedded into the two isoforms by using nonlinear dynamics analysis, which shows that the chaoticity of the system increases at intermediate force (0.4 to 1.6 pN). The identification of chaos in our investigation provides useful information for a comprehensive explanation of the origin of the complex behavior of the system, which effectively helps us to perceive the dynamics of IHF bound 4WHJs under the influence of external force, and also demonstrates the applicability of nonlinear dynamics analysis in the field of biology.

Receptor-Ligand Rebinding Kinetics in Confinement

Rebinding kinetics of molecular ligands plays a key role in the operation of biomachinery, from regulatory networks to protein transcription, and is also a key factor in design of drugs and high-precision biosensors. In this study, we investigate initial release and rebinding of ligands to their binding sites grafted on a planar surface, a situation commonly observed in single-molecule experiments and that occurs in vivo, e.g., during exocytosis. Via scaling arguments and molecular dynamic simulations, we analyze the dependence of nonequilibrium rebinding kinetics on two intrinsic length scales: the average separation distance between the binding sites and the total diffusible volume (i.e., height of the experimental reservoir in which diffusion takes place or average distance between receptor-bearing surfaces). We obtain time-dependent scaling laws for on rates and for the cumulative number of rebinding events. For diffusion-limited binding, the (rebinding) on rate decreases with time via multiple power-law regimes before the terminal steady-state (constant on-rate) regime. At intermediate times, when particle density has not yet become uniform throughout the diffusible volume, the cumulative number of rebindings exhibits a novel, to our knowledge, plateau behavior because of the three-dimensional escape process of ligands from binding sites. The duration of the plateau regime depends on the average separation distance between binding sites. After the three-dimensional diffusive escape process, a one-dimensional diffusive regime describes on rates. In the reaction-limited scenario, ligands with higher affinity to their binding sites (e.g., longer residence times) delay entry to the power-law regimes. Our results will be useful for extracting hidden timescales in experiments such as kinetic rate measurements for ligand-receptor interactions in microchannels, as well as for cell signaling via diffusing molecules.

Characterizing Watson-Crick versus Hoogsteen Base Pairing in a DNA-Protein Complex Using Nuclear Magnetic Resonance and Site-Specifically 13C- and 15N-Labeled DNA

A( syn)-T and G( syn)-C⁺ Hoogsteen base pairs in protein-bound DNA duplexes can be difficult to resolve by X-ray crystallography due to ambiguous electron density and by nuclear magnetic resonance (NMR) spectroscopy due to poor chemical shift dispersion and size limitations with solution-state NMR spectroscopy. Here we describe an NMR strategy for characterizing Hoogsteen base pairs in protein-DNA complexes, which relies on site-specifically incorporating ¹³C- and ¹⁵N-labeled nucleotides into DNA duplexes for unambiguous resonance assignment and to improve spectral resolution. The approach was used to resolve the conformation of an A-T base pair in a crystal structure of an ∼43 kDa complex between a 34 bp duplex DNA and the integration host factor (IHF) protein. In the crystal structure (Protein Data Bank entry 1IHF ), this base pair adopts an unusual Hoogsteen conformation with a distorted sugar backbone that is accommodated by a nearby nick used to aid in crystallization. The NMR chemical shifts and interproton nuclear Overhauser effects indicate that this base pair predominantly adopts a Watson-Crick conformation in the intact DNA-IHF complex under solution conditions. Consistent with these NMR findings, substitution of 7-deazaadenine at this base pair resulted in only a small (∼2-fold) decrease in the IHF-DNA binding affinity. The NMR strategy provides a new approach for resolving crystallographic ambiguity and more generally for studying the structure and dynamics of protein-DNA complexes in solution.

Effects of electrostatic interactions on ligand dissociation kinetics

We study unbinding of multivalent cationic ligands from oppositely charged polymeric binding sites sparsely grafted on a flat neutral substrate. Our molecular dynamics simulations are suggested by single-molecule studies of protein-DNA interactions. We consider univalent salt concentrations spanning roughly a 1000-fold range, together with various concentrations of excess ligands in solution. To reveal the ionic effects on unbinding kinetics of spontaneous and facilitated dissociation mechanisms, we treat electrostatic interactions both at a Debye-Hückel (DH) (or implicit ions, i.e., use of an electrostatic potential with a prescribed decay length) level and by the more precise approach of considering all ionic species explicitly in the simulations. We find that the DH approach systematically overestimates unbinding rates, relative to the calculations where all ion pairs are present explicitly in solution, although many aspects of the two types of calculation are qualitatively similar. For facilitated dissociation (FD) (acceleration of unbinding by free ligands in solution) explicit-ion simulations lead to unbinding at lower free-ligand concentrations. Our simulations predict a variety of FD regimes as a function of free-ligand and ion concentrations; a particularly interesting regime is at intermediate concentrations of ligands where nonelectrostatic binding strength controls FD. We conclude that explicit-ion electrostatic modeling is an essential component to quantitatively tackle problems in molecular ligand dissociation, including nucleic-acid-binding proteins.

RsmW, Pseudomonas aeruginosa small non-coding RsmA-binding RNA upregulated in biofilm versus planktonic growth conditions

BACKGROUND

Biofilm development, specifically the fundamentally adaptive switch from acute to chronic infection phenotypes, requires global regulators and small non-coding regulatory RNAs (sRNAs). This work utilized RNA-sequencing (RNA-seq) to detect sRNAs differentially expressed in Pseudomonas aeruginosa biofilm versus planktonic state.

RESULTS

A computational algorithm was devised to detect and categorize sRNAs into 5 types: intergenic, intragenic, 5'-UTR, 3'-UTR, and antisense. Here we report a novel RsmY/RsmZ-type sRNA, termed RsmW, in P. aeruginosa up-transcribed in biofilm versus planktonic growth. RNA-Seq, 5'-RACE and Mfold predictions suggest RsmW has a secondary structure with 3 of 7 GGA motifs located on outer stem loops. Northern blot revealed two RsmW binding bands of 400 and 120 bases, suggesting RsmW is derived from the 3'-UTR of the upstream hypothetical gene, PA4570. RsmW expression is elevated in late stationary versus logarithmic growth phase in PB minimal media, at higher temperatures (37 °C versus 28 °C), and in both gacA and rhlR transposon mutants versus wild-type. RsmW specifically binds to RsmA protein in vitro and restores biofilm production and reduces swarming in an rsmY/rsmZ double mutant. PA4570 weakly resembles an RsmA/RsmN homolog having 49 % and 51 % similarity, and 16 % and 17 % identity to RsmA and RsmN amino acid sequences, respectively. PA4570 was unable to restore biofilm and swarming phenotypes in ΔrsmA deficient strains.

CONCLUSION

Collectively, our study reveals an interesting theme regarding another sRNA regulator of the Rsm system and further unravels the complexities regulating adaptive responses for Pseudomonas species.

How to determine local stretching and tension in a flow-stretched DNA molecule

We determine the nonuniform stretching of and tension in a mega base pairs-long fragment of deoxyribonucleic acid (DNA) that is flow stretched in a nanofluidic chip. We use no markers, do not know the contour length of the DNA, and do not have the full DNA molecule inside our field of view. Instead, we analyze the transverse thermal motion of the DNA. Tension at the center of the DNA adds up to 16 pN, giving almost fully stretched DNA. This method was devised for optical mapping of DNA, specifically, DNA denaturation patterns. It may be useful also for other studies, e.g., DNA-protein interactions, specifically, their tension dependence. Generally, wherever long strands of DNA-e.g., native DNA extracted from human cells or bacteria-must be stretched with ease for inspection, this method applies.

Twist-open mechanism of DNA damage recognition by the Rad4/XPC nucleotide excision repair complex

DNA damage repair starts with the recognition of damaged sites from predominantly normal DNA. In eukaryotes, diverse DNA lesions from environmental sources are recognized by the xeroderma pigmentosum C (XPC) nucleotide excision repair complex. Studies of Rad4 (radiation-sensitive 4; yeast XPC ortholog) showed that Rad4 "opens" up damaged DNA by inserting a β-hairpin into the duplex and flipping out two damage-containing nucleotide pairs. However, this DNA lesion "opening" is slow (˜5-10 ms) compared with typical submillisecond residence times per base pair site reported for various DNA-binding proteins during 1D diffusion on DNA. To address the mystery as to how Rad4 pauses to recognize lesions during diffusional search, we examine conformational dynamics along the lesion recognition trajectory using temperature-jump spectroscopy. Besides identifying the ˜10-ms step as the rate-limiting bottleneck towards opening specific DNA site, we uncover an earlier ˜100- to 500-μs step that we assign to nonspecific deformation (unwinding/"twisting") of DNA by Rad4. The β-hairpin is not required to unwind or to overcome the bottleneck but is essential for full nucleotide-flipping. We propose that Rad4 recognizes lesions in a step-wise "twist-open" mechanism, in which preliminary twisting represents Rad4 interconverting between search and interrogation modes. Through such conformational switches compatible with rapid diffusion on DNA, Rad4 may stall preferentially at a lesion site, offering time to open DNA. This study represents the first direct observation, to our knowledge, of dynamical DNA distortions during search/interrogation beyond base pair breathing. Submillisecond interrogation with preferential stalling at cognate sites may be common to various DNA-binding proteins.

The λ Integrase Site-specific Recombination Pathway

The site-specific recombinase encoded by bacteriophage λ (Int) is responsible for integrating and excising the viral chromosome into and out of the chromosome of its Escherichia coli host. Int carries out a reaction that is highly directional, tightly regulated, and depends upon an ensemble of accessory DNA bending proteins acting on 240 bp of DNA encoding 16 protein binding sites. This additional complexity enables two pathways, integrative and excisive recombination, whose opposite, and effectively irreversible, directions are dictated by different physiological and environmental signals. Int recombinase is a heterobivalent DNA binding protein and each of the four Int protomers, within a multiprotein 400 kDa recombinogenic complex, is thought to bind and, with the aid of DNA bending proteins, bridge one arm- and one core-type DNA site. In the 12 years since the publication of the last review focused solely on the λ site-specific recombination pathway in Mobile DNA II, there has been a great deal of progress in elucidating the molecular details of this pathway. The most dramatic advances in our understanding of the reaction have been in the area of X-ray crystallography where protein-DNA structures have now been determined for of all of the DNA-protein interfaces driving the Int pathway. Building on this foundation of structures, it has been possible to derive models for the assembly of components that determine the regulatory apparatus in the P-arm, and for the overall architectures that define excisive and integrative recombinogenic complexes. The most fundamental additional mechanistic insights derive from the application of hexapeptide inhibitors and single molecule kinetics.

Regulation of Bacterial DNA Packaging in Early Stationary Phase by Competitive DNA Binding of Dps and IHF

The bacterial nucleoid, a bacterial genome packed by nucleoid binding proteins, forms the physical basis for cellular processes such as gene transcription and DNA replication. Bacteria need to dynamically modulate their nucleoid structures at different growth phases and in response to environmental changes. At the nutrients deficient stationary phase, DNA-binding proteins from starved cells (Dps) and Integration host factors (IHF) are the two most abundant nucleoid associated proteins in E. coli. Yet, it remains unclear how the nucleoid architecture is controlled by the interplay between these two proteins, as well as the nucleoid's response to environmental changes. This question is addressed here using single DNA manipulation approach. Our results reveal that the two proteins are differentially selected for DNA binding, which can be tuned by changing environmental factors over physiological ranges including KCl (50-300 mM), MgCl2 (0-10 mM), pH (6.5-8.5) and temperature (23-37 °C). Increasing pH and MgCl2 concentrations switch from Dps-binding to IHF-binding. Stable Dps-DNA and IHF-DNA complexes are insensitive to temperature changes for the range tested. The environment dependent selection between IHF and Dps results in different physical organizations of DNA. Overall, our findings provide important insights into E. coli nucleoid architecture.

Coupled binding-bending-folding: The complex conformational dynamics of protein-DNA binding studied by atomistic molecular dynamics simulations

BACKGROUND

Protein-DNA binding often involves dramatic conformational changes such as protein folding and DNA bending. While thermodynamic aspects of this behavior are understood, and its biological function is often known, the mechanism by which the conformational changes occur is generally unclear. By providing detailed structural and energetic data, molecular dynamics simulations have been helpful in elucidating and rationalizing protein-DNA binding.

SCOPE OF REVIEW

This review will summarize recent atomistic molecular dynamics simulations of the conformational dynamics of DNA and protein-DNA binding. A brief overview of recent developments in DNA force fields is given as well.

MAJOR CONCLUSIONS

Simulations have been crucial in rationalizing the intrinsic flexibility of DNA, and have been instrumental in identifying the sequence of binding events, the triggers for the conformational motion, and the mechanism of binding for a number of important DNA-binding proteins.

GENERAL SIGNIFICANCE

Molecular dynamics simulations are an important tool for understanding the complex binding behavior of DNA-binding proteins. With recent advances in force fields and rapid increases in simulation time scales, simulations will become even more important for future studies. This article is part of a Special Issue entitled Recent developments of molecular dynamics.

Nucleoprotein architectures regulating the directionality of viral integration and excision

The virally encoded site-specific recombinase Int collaborates with its accessory DNA bending proteins IHF, Xis, and Fis to assemble two distinct, very large, nucleoprotein complexes that carry out either integrative or excisive recombination along regulated and essentially unidirectional pathways. The core of each complex consists of a tetramer of Integrase protein (Int), which is a heterobivalent DNA binding protein that binds and bridges a core-type DNA site (where strand cleavage and ligation are executed), and a distal arm-type site, that is brought within range by one or more DNA bending proteins. The recent determination of the patterns of these Int bridges has made it possible to think realistically about the global architecture of the recombinogenic complexes. Here, we combined the previously determined Int bridging patterns with in-gel FRET experiments and in silico modeling to characterize and differentiate the two 400-kDa multiprotein Holiday junction recombination intermediates formed during λ integration and excision. The results lead to architectural models that explain how integration and excision are regulated in λ site-specific recombination. Our confidence in the basic features of these architectures is based on the redundancy and self-consistency of the underlying data from two very different experimental approaches to establish bridging interactions, a set of strategic intracomplex distances from FRET experiments, and the model's ability to explain key aspects of the integrative and excisive recombination pathways, such as topological changes, the mechanism of capturing attB, and the features of asymmetry and flexibility within the complexes.

Global analysis of ion dependence unveils hidden steps in DNA binding and bending by integration host factor

Proteins that recognize and bind to specific sites on DNA often distort the DNA at these sites. The rates at which these DNA distortions occur are considered to be important in the ability of these proteins to discriminate between specific and nonspecific sites. These rates have proven difficult to measure for most protein-DNA complexes in part because of the difficulty in separating the kinetics of unimolecular conformational rearrangements (DNA bending and kinking) from the kinetics of bimolecular complex association and dissociation. A notable exception is the Integration Host Factor (IHF), a eubacterial architectural protein involved in chromosomal compaction and DNA recombination, which binds with subnanomolar affinity to specific DNA sites and bends them into sharp U-turns. The unimolecular DNA bending kinetics has been resolved using both stopped-flow and laser temperature-jump perturbation. Here we expand our investigation by presenting a global analysis of the ionic strength dependence of specific binding affinity and relaxation kinetics of an IHF-DNA complex. This analysis enables us to obtain each of the underlying elementary rates (DNA bending/unbending and protein-DNA association/dissociation), and their ionic strength dependence, even under conditions where the two processes are coupled. Our analysis indicates interesting differences in the ionic strength dependence of the bi- versus unimolecular steps. At moderate [KCl] (100-500 mM), nearly all the ionic strength dependence to the overall equilibrium binding affinity appears in the bimolecular association/dissociation of an initial, presumably weakly bent, encounter complex, with a slope SK(bi) ≈ 8 describing the loglog-dependence of the equilibrium constant to form this complex on [KCl]. In contrast, the unimolecular equilibrium constant to form the fully wrapped specific complex from the initial complex is nearly independent of [KCl], with SK(uni) < 0.5. This result is counterintuitive because there are at least twice as many ionic protein-DNA contacts in the fully wrapped complex than in the weakly bent intermediate. The following picture emerges from this analysis: in the bimolecular step, the observed [KCl]-dependence is consistent with the number of DNA counterions expected to be released when IHF binds nonspecifically to DNA whereas in the unimolecular reorganization step, the weak [KCl]-dependence suggests that two effects cancel one another. On one hand, formation of additional protein-DNA contacts in the fully wrapped complex releases bound counterions into bulk solution, which is entropically favored by decreasing [salt]. On the other hand, formation of the fully wrapped complex also releases tightly bound water molecules, which is osmotically favored by increasing [salt]. More generally, our global analysis strategy is applicable to other protein-DNA complexes, and opens up the possibility of measuring DNA bending rates in complexes where the unimolecular and bimolecular steps are not easily separable.

Multiple-binding-site mechanism explains concentration-dependent unbinding rates of DNA-binding proteins

Recent work has demonstrated concentration-dependent unbinding rates of proteins from DNA, using fluorescence visualization of the bacterial nucleoid protein Fis [Graham et al. (2011) (Concentration-dependent exchange accelerates turnover of proteins bound to double-stranded DNA. Nucleic Acids Res., 39:2249)]. The physical origin of this concentration-dependence is unexplained. We use a combination of coarse-grained simulation and theory to demonstrate that this behavior can be explained by taking into account the dimeric nature of the protein, which permits partial dissociation and exchange with other proteins in solution. Concentration-dependent unbinding is generated by this simple model, quantitatively explaining experimental data. This effect is likely to play a major role in determining binding lifetimes of proteins in vivo where there are very high concentrations of solvated molecules.

Mechanosensing of DNA bending in a single specific protein-DNA complex

Many crucial biological processes are regulated by mechanical stimuli. Here, we report new findings that pico-Newton forces can drastically affect the stability of the site-specific DNA binding of a single transcription factor, the E. coli integration host factor (IHF), by stretching a short ~150 nm DNA containing a single IHF binding site. Dynamic binding and unbinding of single IHF were recorded and analyzed for the force-dependent stability of the IHF-DNA complex. Our results demonstrate that the IHF-DNA interaction is fine tuned by force in different salt concentration and temperature over physiological ranges, indicating that, besides other physiological factors, force may play equally important role in transcription regulation. These findings have broad implications with regard to general mechanosensitivity of site-specific DNA bending proteins.

Physical organization of DNA by multiple non-specific DNA-binding modes of integration host factor (IHF)

The integration host factor (IHF) is an abundant nucleoid-associated protein and an essential co-factor for phage λ site-specific recombination and gene regulation in E. coli. Introduction of a sharp DNA kink at specific cognate sites is critical for these functions. Interestingly, the intracellular concentration of IHF is much higher than the concentration needed for site-specific interactions, suggesting that non-specific binding of IHF to DNA plays a role in the physical organization of bacterial chromatin. However, it is unclear how non-specific DNA association contributes to DNA organization. By using a combination of single DNA manipulation and atomic force microscopy imaging methods, we show here that distinct modes of non-specific DNA binding of IHF result in complex global DNA conformations. Changes in KCl and IHF concentrations, as well as tension applied to DNA, dramatically influence the degree of DNA-bending. In addition, IHF can crosslink DNA into a highly compact DNA meshwork that is observed in the presence of magnesium at low concentration of monovalent ions and high IHF-DNA stoichiometries. Our findings provide important insights into how IHF contributes to bacterial chromatin organization, gene regulation, and biofilm formation.

Single-molecule spectroscopic study of dynamic nanoscale DNA bending behavior of HIV-1 nucleocapsid protein

We have studied the conformational dynamics associated with the nanoscale DNA bending induced by human immunodeficiency virus type 1 (HIV-1) nucleocapsid (NC) protein using single-molecule Förster resonance energy transfer (SM-FRET). To gain molecular-level insights into how the HIV-1 NC locally distorts the structures of duplexed DNA segments, the dynamics, reversibility, and sequence specificity of the DNA bending behavior of NC have been systematically studied. We have performed SM-FRET measurements on a series of duplexed DNA segments with varying sequences, lengths, and local structures in the presence of the wide-type HIV-1 NC and NC mutants lacking either the basic N-terminal domain or the zinc fingers. On the basis of the SM-FRET results, we have proposed a possible mechanism for the NC-induced DNA bending in which both NC's zinc fingers and N-terminal domain are found to play crucial roles. The SM-FRET results reported here add new mechanistic insights into the biological behaviors and functions of HIV-1 NC as a retroviral DNA-architectural protein which may play critical roles in the compaction, nuclear import, and integration of the proviral DNA during the retroviral life cycle.

Structures during binding of cAMP receptor to promoter DNA: promoter search slowed by non-specific sites

The kinetics of cAMP receptor (CAP) binding to promoter DNA has been studied by stopped-flow electric-dichroism at a reduced salt concentration, where the coupling of non-specific and specific binding can be observed directly. Amplitudes, rise and decay times of dichroism transients provide detailed information about the reaction and the structure of intermediates over more than six orders of magnitude on the time scale. CAP binding during the first milliseconds after mixing is indicated by an increase of both rise- and decay-time constants. A particularly large increase of rise times reflects initial formation of non-symmetric complexes by protein binding to non-specific sites at DNA ends. The increase of the hydrodynamic dimensions continues up to ~1 s, before a decrease of time constants reflects transition to compact states with bent DNA up to the time range of ~10(3) s. The slow approach to CAP-induced DNA bending is due to non-specific complexes, which are formed initially and are converted slowly to the specific complex. At the salt concentration of 13.5 mM, conversion to specific complexes with bent DNA is completed after ~40 s at pH 8 compared to >10(3) s at pH 7, resulting from a higher affinity of CAP to non-specific sites at pH 7 than 8 by a factor of ~100. Thus, under the given conditions non-specific sites delay rather than facilitate formation of the specific complex with bent DNA. Experimental data obtained for a non-specific DNA clearly indicate the impact of pseudo-sites. The different electro-optical parameters have been combined in global fits.

Prokaryotic genome regulation: a revolutionary paradigm

After determination of the whole genome sequence, the research frontier of bacterial molecular genetics has shifted to reveal the genome regulation under stressful conditions in nature. The gene selectivity of RNA polymerase is modulated after interaction with two groups of regulatory proteins, 7 sigma factors and 300 transcription factors. For identification of regulation targets of transcription factors in Escherichia coli, we have developed Genomic SELEX system and subjected to screening the binding sites of these factors on the genome. The number of regulation targets by a single transcription factor was more than those hitherto recognized, ranging up to hundreds of promoters. The number of transcription factors involved in regulation of a single promoter also increased to as many as 30 regulators. The multi-target transcription factors and the multi-factor promoters were assembled into complex networks of transcription regulation. The most complex network was identified in the regulation cascades of transcription of two master regulators for planktonic growth and biofilm formation.

Human mitochondrial transcription factor A induces a U-turn structure in the light strand promoter

Human mitochondrial transcription factor A, TFAM, is essential for mitochondrial DNA packaging and maintenance and also has a crucial role in transcription. Crystallographic analysis of TFAM in complex with an oligonucleotide containing the mitochondrial light strand promoter (LSP) revealed two high-mobility group (HMG) protein domains that, through different DNA recognition properties, intercalate residues at two inverted DNA motifs. This induced an overall DNA bend of ~180°, stabilized by the interdomain linker. This U-turn allows the TFAM C-terminal tail, which recruits the transcription machinery, to approach the initiation site, despite contacting a distant DNA sequence. We also ascertained that structured protein regions contacting DNA in the crystal were highly flexible in solution in the absence of DNA. Our data suggest that TFAM bends LSP to create an optimal DNA arrangement for transcriptional initiation while facilitating DNA compaction elsewhere in the genome.

Basic statistical recipes for the emergence of biochemical discernment

An essential step towards understanding life would be to identify the very basic mechanisms responsible for the discerning behaviour of living biochemical systems, absent from randomly reacting chemical soups. One intuitively feels that this question goes beyond the particular nature of the biological molecules and should relate to general physical principles. The pre-eminent physicist Ludwig Boltzmann early envisioned life as a struggle for entropy, in concordance with the subsequent principle of self-organization out of equilibrium. Re-examination of elementary steady state biochemical systems from a statistical perspective supports this view and shows that sigmoidal responses arising from microstates elimination, are sufficient to explain innermost characteristics of life, including its capacity to convert random molecular interactions into accurate biological reactions. A primary operating strategy to achieve this goal is the introduction of time-irreversible transitions in molecular state conversion cycles by injection of free energy, which confers decisional capacity to single macromolecules. Selected examples from various fields of molecular biology such as enzymology and gene expression, are provided to show that these non-equilibrium steady state mechanisms remain important in contemporary biochemical systems. But in addition, information archiving allowed the emergence of the time-reversible counterparts of these mechanisms, mediated by evolutionary pre-organized macromolecular complexes capable of generating discernment in a non-dissipative manner.

Force-driven unbinding of proteins HU and Fis from DNA quantified using a thermodynamic Maxwell relation

Determining numbers of proteins bound to large DNAs is important for understanding their chromosomal functions. Protein numbers may be affected by physical factors such as mechanical forces generated in DNA, e.g. by transcription or replication. We performed single-DNA stretching experiments with bacterial nucleoid proteins HU and Fis, verifying that the force-extension measurements were in thermodynamic equilibrium. We, therefore, could use a thermodynamic Maxwell relation to deduce the change of protein number on a single DNA due to varied force. For the binding of both HU and Fis under conditions studied, numbers of bound proteins decreased as force was increased. Our experiments showed that most of the bound HU proteins were driven off the DNA at 6.3 pN for HU concentrations lower than 150 nM; our HU data were fit well by a statistical-mechanical model of protein-induced bending of DNA. In contrast, a significant amount of Fis proteins could not be forced off the DNA at forces up to 12 pN and Fis concentrations up to 20 nM. This thermodynamic approach may be applied to measure changes in numbers of a wide variety of molecules bound to DNA or other polymers. Force-dependent DNA binding by proteins suggests mechano-chemical mechanisms for gene regulation.

Concentration-dependent exchange accelerates turnover of proteins bound to double-stranded DNA

The multistep kinetics through which DNA-binding proteins bind their targets are heavily studied, but relatively little attention has been paid to proteins leaving the double helix. Using single-DNA stretching and fluorescence detection, we find that sequence-neutral DNA-binding proteins Fis, HU and NHP6A readily exchange with themselves and with each other. In experiments focused on the Escherichia coli nucleoid-associated protein Fis, only a small fraction of protein bound to DNA spontaneously dissociates into protein-free solution. However, if Fis is present in solution, we find that a concentration-dependent exchange reaction occurs which turns over the bound protein, with a rate of k_exch = 6 × 10⁴ M⁻¹s⁻¹. The bacterial DNA-binding protein HU and the yeast HMGB protein NHP6A display the same phenomenon of protein in solution accelerating dissociation of previously bound labeled proteins as exchange occurs. Thus, solvated proteins can play a key role in facilitating removal and renewal of proteins bound to the double helix, an effect that likely plays a major role in promoting the turnover of proteins bound to DNA in vivo and, therefore, in controlling the dynamics of gene regulation.

Experimental snapshots of a protein-DNA binding landscape

Protein recognition of DNA sites is a primary event for gene function. Its ultimate mechanistic understanding requires an integrated structural, dynamic, kinetic, and thermodynamic dissection that is currently limited considering the hundreds of structures of protein-DNA complexes available. We describe a protein-DNA-binding pathway in which an initial, diffuse, transition state ensemble with some nonnative contacts is followed by formation of extensive nonnative interactions that drive the system into a kinetic trap. Finally, nonnative contacts are slowly rearranged into native-like interactions with the DNA backbone. Dissimilar protein-DNA interfaces that populate along the DNA-binding route are explained by a temporary degeneracy of protein-DNA interactions, centered on "dual-role" residues. The nonnative species slow down the reaction allowing for extended functionality.

Platform for in situ real-time measurement of protein-induced conformational changes of DNA

A platform for in situ and real-time measurement of protein-induced conformational changes in dsDNA is presented. We combine electrical orientation of surface-bound dsDNA probes with an optical technique to measure the kinetics of DNA conformational changes. The sequence-specific Escherichia coli integration host factor is utilized to demonstrate protein-induced bending upon binding of integration host factor to dsDNA probes. The effects of probe surface density on binding/bending kinetics are investigated. The platform can accommodate individual spots of microarrayed dsDNA on individually controlled, lithographically designed electrodes, making it amenable for use as a high throughput assay.

Quantitative analyses of RAG-RSS interactions and conformations revealed by atomic force microscopy

During V(D)J recombination, site specific DNA excision is dictated by the binding of RAG1/2 proteins to the conserved recombination signal sequence (RSS) within the genome. The interaction between RAG1/2 and RSS is thought to involve a large DNA distortion that is permissive for DNA cleavage. In this study, using atomic force microscopy imaging (AFM), we analyzed individual RAG-RSS complexes, in which the bending angle of RAG-associated RSS substrates could be visualized and quantified. We provided the quantitative measurement on the conformations of specific RAG-12RSS complexes. Previous data indicating the necessity of RAG2 for recombination implies a structural role in the RAG-RSS complex. Surprisingly, however, no significant difference was observed in conformational bending with AFM between RAG1-12RSS and RAG1/2-12RSS. RAG1 was found sufficient to induce DNA bending, and the addition of RAG2 did not change the bending profile. In addition, a prenicked 12RSS bound by RAG1/2 proteins displayed a conformation similar to the one observed with the intact 12RSS, implying that no greater DNA bending occurs after the nicking step in the signal complex. Taken together, the quantitative AFM results on the components of the recombinase emphasize a tightly held complex with a bend angle value near 60 degrees , which may be a prerequisite step for the site-specific nicking by the V(D)J recombinase.

AAA+ ATPases in the initiation of DNA replication

All cellular organisms and many viruses rely on large, multi-subunit molecular machines, termed replisomes, to ensure that genetic material is accurately duplicated for transmission from one generation to the next. Replisome assembly is facilitated by dedicated initiator proteins, which serve to both recognize replication origins and recruit requisite replisomal components to the DNA in a cell-cycle coordinated manner. Exactly how imitators accomplish this task, and the extent to which initiator mechanisms are conserved among different organisms have remained outstanding issues. Recent structural and biochemical findings have revealed that all cellular initiators, as well as the initiators of certain classes of double-stranded DNA viruses, possess a common adenine nucleotide-binding fold belonging to the ATPases Associated with various cellular Activities (AAA+) family. This review focuses on how the AAA+ domain has been recruited and adapted to control the initiation of DNA replication, and how the use of this ATPase module underlies a common set of initiator assembly states and functions. How biochemical and structural properties correlate with initiator activity, and how species-specific modifications give rise to unique initiator functions, are also discussed.

Transition state for protein-DNA recognition

We describe the formation of protein-DNA contacts in the two-state route for DNA sequence recognition by a transcriptional regulator. Surprisingly, direct sequence readout establishes in the transition state and constitutes the bottleneck of complex formation. Although a few nonspecific ionic interactions are formed at this early stage, they mainly play a stabilizing role in the final consolidated complex. The interface is fairly plastic in the transition state, likely because of a high level of hydration. The overall picture of this two-state route largely agrees with a smooth energy landscape for binding that speeds up DNA recognition. This "direct" two-state route differs from the parallel multistep pathway described for this system, which involves nonspecific contacts and at least two intermediate species that must involve substantial conformational rearrangement in either or both macromolecules.

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(alphabeta), 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(alphabeta) 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 (> or =34 bp), at low [HU]/[DNA], HU populates a noncooperative 34 bp binding mode with a binding constant of 2.1+/-0.4x10(6) M(-1), and a binding enthalpy of +7.7+/-0.6 kcal/mol at 15 degrees C and 0.15 M Na(+). With increasing [HU]/[DNA], HU bound in the noncooperative 34 bp mode progressively converts to two cooperative (omega approximately 20) modes with site sizes of 10 bp and 6 bp. These latter modes exhibit smaller binding constants (1.1+/-0.2x10(5) M(-1) for the 10 bp mode, 3.5+/-1.4x10(4) 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 degrees 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.

New insights into the transition pathway from nonspecific to specific complex of DNA with Escherichia coli integration host factor

To elucidate the nature of the transition-state ensemble along the reaction pathway from a nonspecific protein-DNA complex to the specific complex, we have carried out measurements of DNA bending/unbending dynamics on a cognate DNA substrate in complex with integration host factor (IHF), an architectural protein from E. coli that bends its cognate site by approximately 180 degrees . We use a laser temperature jump to perturb the IHF-DNA complex and monitor the relaxation kinetics with time-resolved FRET measurements on DNA substrates end-labeled with a FRET pair. Previously, we showed that spontaneous bending/kinking of DNA, from thermal disruption of base-pairing/-stacking interactions, may be the rate-limiting step in the formation of the specific complex (Kuznetsov, S. V.; Sugimura, S.; Vivas, P.; Crothers, D. M.; Ansari, A. Proc. Natl. Acad. Sci. USA 2006, 103, 18515). Here, we probe the effect of varying [KCl], which affects the stability of the complex, on this rate-limiting step. We find that below approximately 250 mM KCl, the observed relaxation kinetics are from the unimolecular bending/unbending of DNA, and the relaxation rate kr is independent of [KCl]. Above approximately 300 mM KCl, dissociation of the IHF-DNA complex becomes significant, and the observed relaxation process includes contributions from the association/dissociation step, with kr decreasing with increasing [KCl]. The DNA bending step occurs with a positive activation enthalpy, despite the large negative enthalpy change reported for the specific IHF-DNA complex (Holbrook, J. A.; Tsodikov, O. V.; Saecker, R. M.; Record, M. T., Jr. J. Mol. Biol. 2001, 310, 379). Our conclusion from these studies is that in the uphill climb to the transition state, the DNA is kinked, but with no release of ions, as indicated by the salt-independent behavior of k(r) at low [KCl]. Any release of ions in the unimolecular process, together with conformational changes in the protein-DNA complex that facilitate favorable interactions and that contribute to the negative enthalpy change, must occur as the system leaves the transition state, downhill to the final complex.

DNA A-tracts bending: polarization effects on electrostatic interactions across their minor groove

Bending by the DNA A-tracts constitutes a contentious issue, suggesting deficiencies in the physics employed so far. Here, we inquire as to the importance in this bending of many-body polarization effects on the electrostatic interactions across their narrow minor groove. We have done this on the basis of the findings of Jarque and Buckingham who developed a procedure based on a Monte Carlo simulation for two charges of the same sign embedded in a polarizable medium. Remarkably, the present analysis reveals that for compact DNA conformations, which result from dynamic effects, an overall attractive interaction operates between the phosphate charges; this interaction is especially strong for the narrow minor groove of the A-tracts, suggesting a tendency for DNA to bend toward this groove. This tendency is in agreement with the conclusions of electrophoretic and NMR solution studies. The present analysis is also consistent with the experimental observations that the minor groove is much more easily compressible than the major groove and the bending propensity of the A-tracts is greatly reduced at "premelting" temperatures. By contrast, the dielectric screening model predicts a repulsion between the phosphate charges and is not consistent with the aforementioned bending tendency or experimental observations.